计算化学公社

标题: Gaussian与VASP处理自旋的差异？ [打印本页]

作者
Author: archer 时间: 2019-12-5 22:14
标题: Gaussian与VASP处理自旋的差异？
对于分子计算，Gaussian需要设置自旋多重度，而VASP只需要设置自旋极化就行了。VASP能够自行收敛到能量最低的自旋态吗？

作者
Author: 卡开发发 时间: 2019-12-6 07:27
VASP不能，一般能做变磁矩（VMM）计算的程序也不能。这主要由算法决定，大部分优化或求根的算法通常找到的都是局部解，要找其他的解常规手段就是调整初猜，因此进行VMM计算猜测磁矩是必须的。

一些例子：
（1）分子FeO+的就有Fe高、低自旋两种电子组态，如果猜高自旋的情况，低自旋的结果显然就得不到了。
（2）不设置初始磁矩的情况，某些程序会按照铁磁性（比如每个原子磁矩为1）处理，有些程序则默认每个原子初始磁矩为0，这样计算铁磁性的体系后者磁性做出来可能不对或者SCF收敛不了。

作者
Author: archer 时间: 2019-12-6 10:46

卡开发发发表于 2019-12-6 07:27
VASP不能，一般能做变磁矩（VMM）计算的程序也不能。这主要由算法决定，大部分优化或求根的算法通常找到的 ...

您好！我现在做大量的结构计算，只有一个金属原子，无法一个一个去找最稳定自旋态，不偏向于研究磁性。如果用VASP不设置磁矩，能够收敛就行，请问这样做是否有价值？

作者
Author: 卡开发发 时间: 2019-12-6 11:50

archer 发表于 2019-12-6 10:46
您好！我现在做大量的结构计算，只有一个金属原子，无法一个一个去找最稳定自旋态，不偏向于研究磁性。如 ...

看啥样的体系吧，一个原子的话，也没有上面出现那种要区分高低自旋态的情况，应该能简单处理。上面提到过如果不手动猜测会是什么样的机制了。

作者
Author: archer 时间: 2019-12-6 12:15

卡开发发发表于 2019-12-6 11:50
看啥样的体系吧，一个原子的话，也没有上面出现那种要区分高低自旋态的情况，应该能简单处理。上面提到过 ...

老师，您好！如果是过渡金属被两个分子的夹心结构，能否猜测是最低自旋态，用Gaussian算这样可以吗？

作者
Author: 卡开发发 时间: 2019-12-6 12:55

archer 发表于 2019-12-6 12:15
老师，您好！如果是过渡金属被两个分子的夹心结构，能否猜测是最低自旋态，用Gaussian算这样可以吗？

不好说。如果你过渡金属周围化学环境都差不多，可以简单测试几个磁矩，其他照做。

作者
Author: Q-Chembio-llg 时间: 2021-7-8 20:54

卡开发发发表于 2019-12-6 12:55
不好说。如果你过渡金属周围化学环境都差不多，可以简单测试几个磁矩，其他照做。

老师好，想请教两个问题：
1.我最近在用VASP优化包含g-c3n4和双Fe，O2的体系时，INCAR中并未设置MAGMOM，但是最后算出来的mag=-2，想知道这个正常吗？
2.我在用VASP计算一个氧气在双铁上解离的路径时，发现反应物、过渡态、产物的mag分别为4、-2和8，这样净自旋不同这个反应路径还有意义吗？
谢谢老师，还请老师不吝赐教。

作者
Author: 卡开发发 时间: 2021-7-8 22:16

Q-Chembio-llg 发表于 2021-7-8 20:54
老师好，想请教两个问题：
1.我最近在用VASP优化包含g-c3n4和双Fe，O2的体系时，INCAR中并未设置MAGMOM ...

1、正常，总磁矩的正负没关系。一般还是建议设置，有时候可能能算出来某个磁矩，但可能未必是基态，复杂情况还是要多考虑一些。
2、O2的解离情况是比较复杂，我认为确实很难讲清楚这里面的意义，可能这些结构不在同一势能面上了。

作者
Author: Q-Chembio-llg 时间: 2021-7-8 22:41

卡开发发发表于 2021-7-8 22:16
1、正常，总磁矩的正负没关系。一般还是建议设置，有时候可能能算出来某个磁矩，但可能未必是基态，复杂 ...

谢谢老师！针对第2点还想继续请教一下，那不在一个势能面上是不是就是说这条反应路径和中间的过渡态都没有意义呢？我有尝试过在INCAR中强行把MAGMOM设成和反应物的总磁矩相同，但是好像对输出的总磁矩没有明显的影响，还是和反应物的总磁矩不一样。

作者
Author: 卡开发发 时间: 2021-7-9 10:38

Q-Chembio-llg 发表于 2021-7-8 22:41
谢谢老师！针对第2点还想继续请教一下，那不在一个势能面上是不是就是说这条反应路径和中间的过渡态都没 ...

不同磁矩就是不同的自旋多重度，不过如果对于多参考效应比较强的情况这点可能不太容易讨论。

作者
Author: jiangning198511 时间: 2021-7-9 11:23

Q-Chembio-llg 发表于 2021-7-8 22:41
谢谢老师！针对第2点还想继续请教一下，那不在一个势能面上是不是就是说这条反应路径和中间的过渡态都没 ...

MAGMOM只是初始猜测，最后的结果不一定能保持这个设定

作者
Author: jiangning198511 时间: 2021-7-9 11:24

Q-Chembio-llg 发表于 2021-7-8 22:41
谢谢老师！针对第2点还想继续请教一下，那不在一个势能面上是不是就是说这条反应路径和中间的过渡态都没 ...

VASP里面可以强制保持自旋不变进行计算

作者
Author: sunjk 时间: 2022-6-20 21:34

jiangning198511 发表于 2021-7-9 11:24
VASP里面可以强制保持自旋不变进行计算

你好，请问VASP强制保持自旋不变的参数是哪个？

作者
Author: jiangning198511 时间: 2022-6-23 14:46

sunjk 发表于 2022-6-20 21:34
你好，请问VASP强制保持自旋不变的参数是哪个？

NUPDOWN

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