计算化学公社

标题: 求助：金属团簇的磁矩计算 [打印本页]

作者
Author: yiranfengbai 时间: 2019-3-14 10:37
标题: 求助：金属团簇的磁矩计算
最近测试了Aln(n=2-7)小团簇的磁矩，发现了比较奇怪的结果：n=偶数的团簇总磁矩为0 ub，而n=奇数的团簇总磁矩为1ub。而参考文献中说:n=偶数的团簇总磁矩为2ub，而n=奇数的团簇总磁矩为1ub。
这是为什么呢？计算在20A*20A*20A的盒子中，参数：CASTEP,GGA-PBE， cutoff=300eV, k-point=1。

作者
Author: sobereva 时间: 2019-3-14 11:16
这种孤立体系，而且还是小团簇，用量化程序诸如Gaussian算比用CASTEP强得多得多，而且还能充分利用对称性，结合Multiwfn考察电子结构特征
那两种磁矩值相当于计算时设的自旋多重度不同，只有其中能量最低的才是基态，才对应真实情况

作者
Author: 卡开发发 时间: 2019-3-14 14:38
在Modify-Electronic Configuration里面把每个原子的磁矩都猜+1看看，Al默认可能是0就会算成反铁磁的，还有动能截断最好测试下。

作者
Author: yiranfengbai 时间: 2019-3-14 19:48

sobereva 发表于 2019-3-14 11:16
这种孤立体系，而且还是小团簇，用量化程序诸如Gaussian算比用CASTEP强得多得多，而且还能充分利用对称性， ...

社长，Gaussian计算磁矩是用命令NBO=npa，然后通过alpha与beta电子做差吗?

作者
Author: yiranfengbai 时间: 2019-3-14 19:49

卡开发发发表于 2019-3-14 14:38
在Modify-Electronic Configuration里面把每个原子的磁矩都猜+1看看，Al默认可能是0就会算成反铁磁的，还有 ...

谢谢老师，我给Al都赋予磁矩+1之后，计算出来的磁矩终于正常了。通过这次测试，我感觉以后再计算类似磁性团簇的时候都要先赋予磁矩才行。

作者
Author: 卡开发发 时间: 2019-3-14 21:00
本帖最后由卡开发发于 2019-3-15 04:41 编辑

yiranfengbai 发表于 2019-3-14 19:49
谢谢老师，我给Al都赋予磁矩+1之后，计算出来的磁矩终于正常了。通过这次测试，我感觉以后再计算类似磁性 ...

磁性系统吃不准默认给出的磁矩是啥样的情况建议都手工做初猜。这里CASTEP一部分原子默认是0，很容易算成反铁磁或无磁性，而VASP开自旋极化默认是1，很容易算成铁磁。很多人觉得开了自旋极化就能算出来想要的结果，这个还不一定。

作者
Author: sobereva 时间: 2019-3-15 06:44

yiranfengbai 发表于 2019-3-14 19:48
社长，Gaussian计算磁矩是用命令NBO=npa，然后通过alpha与beta电子做差吗?

跟pop=NPA毫无直接关系

仔细看了此文就懂了，专门说了磁矩
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353（http://bbs.keinsci.com/thread-4422-1-1.html）

作者
Author: yiranfengbai 时间: 2019-3-15 09:32

卡开发发发表于 2019-3-14 21:00
磁性系统吃不准默认给出的磁矩是啥样的情况建议都手工做初猜。这里CASTEP一部分原子默认是0，很容易算成 ...

明白了，谢谢老师

作者
Author: yiranfengbai 时间: 2019-3-15 09:32

sobereva 发表于 2019-3-15 06:44
跟pop=NPA毫无直接关系

仔细看了此文就懂了，专门说了磁矩

谢谢社长了

欢迎光临计算化学公社 (http://ccc.keinsci.com/)