计算化学公社

标题: 请问在计算吸附是否开自旋的问题 [打印本页]

作者
Author: xmuqyj 时间: 2024-3-22 21:42
标题: 请问在计算吸附是否开自旋的问题
本帖最后由 xmuqyj 于 2024-3-22 21:46 编辑

我的体系是氧化亚铜上吸附小分子，体系肯定是开了磁距的，但是请问计算零点能时是否开自旋呢？
我在几个体系计算零点能时不开自旋，发现部分结构电子步能收敛，但是部分结构电子步不能收敛，是不是证明一定要开自旋呢？
我的另一个体系，在二氧化硅负载铜单原子，在上面吸附分子，做了磁性测试，发现吸附H原子的时候必须开磁性，不然不收敛。而吸附大的有机分子，开磁性，在OUTCAR里看原子的磁性都是0.000x，不开磁性也能收敛。
在二氧化硅表面负载铜簇的时候，开磁性在OUTCAR查看磁距，发现所有原子也是0。
所以不太理解，到底这种体系要不要开磁性呢？（当然希望能够节约时间）

作者
Author: wzkchem5 时间: 2024-3-23 16:45
应当根据化学常识，确认你的体系有没有单电子。很显然，即使多一个氢、少一个氢，都会改变体系电子数的奇偶性，而电子数为偶数的情况下，很小的结构变化就可能导致两个本来该配对的电子不再配对，变成单电子。
一定一定不要以是否收敛为标准判断是否应该开磁性，设置错误反倒导致收敛更容易的例子多如牛毛。从化学常识判断是否应当开磁性以后，就应该咬准这个结论不松口，即使不收敛也去通过加其他关键词的方法帮助收敛，而不要再动磁性设置。

作者
Author: xmuqyj 时间: 2024-3-23 21:18

wzkchem5 发表于 2024-3-23 16:45
应当根据化学常识，确认你的体系有没有单电子。很显然，即使多一个氢、少一个氢，都会改变体系电子数的奇偶 ...

感谢老师的回复，我还请教老师，在铜slab上做吸附催化反应，大家一般都是不加磁性的，我做测试也是没有磁性。而在铜单原子上是不是向您说这这种有单电子需要加磁性呢？那么铜团簇呢？有什么需要注意的吗

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