计算化学公社

标题: 单重态氧气计算求助 [打印本页]

作者
Author: cmd27 时间: 2022-3-23 21:40
标题: 单重态氧气计算求助
大家好，我用BLYP/6-31g(d)优化单重态的氧气时，SCF不收敛，我之前用的B3LYP同一基组是可以收敛的，请教大家这是为什么？我附了输入和输出，感谢

作者
Author: wzkchem5 时间: 2022-3-24 02:08
单重态氧气是多参考态体系，本来就不应该用DFT算，而应该用NEVPT2、MRCISD等多参考态方法算。用单参考态方法强行算多参考态体系，不仅难收敛，即使收敛了结果也不可信。
至于为什么B3LYP碰巧能收敛，应该是因为B3LYP的HF交换项成分大，所以HOMO-LUMO gap大，收敛还不至于那么困难。但是还是那句话，虽然B3LYP收敛了，结果也是不能用的，因为这是多参考态体系。

作者
Author: 喵星大佬 时间: 2022-3-24 07:02
本帖最后由喵星大佬于 2022-3-24 07:15 编辑

wzkchem5 发表于 2022-3-24 02:08
单重态氧气是多参考态体系，本来就不应该用DFT算，而应该用NEVPT2、MRCISD等多参考态方法算。用单参考态方 ...

用TDDFT应该是可以描述单重态氧气的
CAS的话最大也就(12,8)，没啥算不了的，实际上只需要(8,6)甚至(6,4)，还有就是用MC-PDFT比NEVPT2/MRCISD啥的便宜得多还有梯度，如果需要后续算反应也可以考虑，不过精度不太好说，而且OpenMolcas用起来也有点麻烦

作者
Author: sobereva 时间: 2022-3-24 09:25
B3LYP当闭壳层算单重态氧气，电子结构完全是错的

作者
Author: cmd27 时间: 2022-3-24 10:47
好的，谢谢各位老师，氧气在我这个体系中是作为反应物的，我想请问一般在反应里面我们都是用三重态去考虑嘛？因为他们实验是用单重态做的，所以我有这个疑问，还是非常感谢各位老师

作者
Author: cmd27 时间: 2022-3-24 10:52

sobereva 发表于 2022-3-24 09:25
B3LYP当闭壳层算单重态氧气，电子结构完全是错的

老师，您好，我想请问在我附件这个图里，第二种单重态氧气，所有电子都是自旋配对的，不是闭壳层嘛？我可能这个问题问的有点傻

作者
Author: cmd27 时间: 2022-3-24 10:52

sobereva 发表于 2022-3-24 09:25
B3LYP当闭壳层算单重态氧气，电子结构完全是错的

C:\Users\lixiaojiao27\Desktop

作者
Author: 喵星大佬 时间: 2022-3-24 10:57

cmd27 发表于 2022-3-24 10:52

显然应该是第三种。。。。。

作者
Author: zjxitcc 时间: 2022-3-24 11:02

cmd27 发表于 2022-3-24 10:52

氧气基态是三重态，称为T0，它的一些激发态S1和S2都是多参考态性质及其显著的。这张图恰如中文、英文网上很多氧气电子排布图一样都有很大的错误（往好听点说，不严谨之处）：Delta_g态至少需要2个行列式，最高的2个电子可以都填充左边的pi*，也可以都填充在右边的pi*；而Sigma_g同样至少需要2个行列式，一个是alpha, beta，另一个是beta, alpha。两行列式比重是一样的，缺了一个行列式就是定性错误。单参考方法没法描述这样的情形。

你先搜一下论坛上以往关于氧气计算的帖子，不计其数，先看了再说。

作者
Author: cmd27 时间: 2022-3-24 11:05

zjxitcc 发表于 2022-3-24 11:02
氧气基态是三重态，称为T0，它的一些激发态S1和S2都是多参考态性质及其显著的。这张图恰如中文、英文网上 ...

好的好的谢谢您受益匪浅

作者
Author: cmd27 时间: 2022-3-24 11:06

喵星大佬发表于 2022-3-24 10:57
显然应该是第三种。。。。。

好的十分感谢您

作者
Author: 小西瓜123 时间: 2022-8-15 10:00

zjxitcc 发表于 2022-3-24 11:02
氧气基态是三重态，称为T0，它的一些激发态S1和S2都是多参考态性质及其显著的。这张图恰如中文、英文网上 ...

老师，请问楼上贴的那张电子结构的图可以从Gaussian或者Molpro的相关输出给出吗？那第三个或者第四个激发态的电子结构应该是怎么样呢？最近看文献，OH+HO2(A')生成H2O+O2(b1)，这个O2(b1)就是图中的第三个结构，然后H+HO2(A')=H2+O2(a1)和H+HO2=H2+O2(A3Δ)，这个O2(a1)就是上图中第二个结构，为什么不同反应物生成的氧的状态不同呢？此外，O2(A3Δ)是有一个电子跃迁到了更高能级了吗？还有一个电子态，比O2(A3Δ)低一点点，为O2(C1)的电子结构应该是什么呢？

作者
Author: tianmafei 时间: 2024-12-18 19:41

zjxitcc 发表于 2022-3-24 11:02
氧气基态是三重态，称为T0，它的一些激发态S1和S2都是多参考态性质及其显著的。这张图恰如中文、英文网上 ...

邹老师，您好！请教您一个问题，想计算基态氧气，得到Sz=0,1,-1的三个子态，从结果看salter行列式的差别，应该用什么方法？先谢谢您！

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