zjxitcc 发表于 2022-3-24 11:02 邹老师,您好!请教您一个问题,想计算基态氧气,得到Sz=0,1,-1的三个子态,从结果看salter行列式的差别,应该用什么方法?先谢谢您! |
zjxitcc 发表于 2022-3-24 11:02 老师,请问楼上贴的那张电子结构的图可以从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)的电子结构应该是什么呢? |
喵星大佬 发表于 2022-3-24 10:57 好的 十分感谢您 |
zjxitcc 发表于 2022-3-24 11:02 好的 好的 谢谢您 受益匪浅 |
|
氧气基态是三重态,称为T0,它的一些激发态S1和S2都是多参考态性质及其显著的。这张图恰如中文、英文网上很多氧气电子排布图一样都有很大的错误(往好听点说,不严谨之处):Delta_g态至少需要2个行列式,最高的2个电子可以都填充左边的pi*,也可以都填充在右边的pi*;而Sigma_g同样至少需要2个行列式,一个是alpha, beta,另一个是beta, alpha。两行列式比重是一样的,缺了一个行列式就是定性错误。单参考方法没法描述这样的情形。 你先搜一下论坛上以往关于氧气计算的帖子,不计其数,先看了再说。 |
|
显然应该是第三种。。。。。 |
sobereva 发表于 2022-3-24 09:25 老师 ,您好,我想请问在我附件这个图里,第二种单重态氧气,所有电子都是自旋配对的,不是闭壳层嘛?我可能这个问题问的有点傻 |
| 好的,谢谢各位老师,氧气在我这个体系中是作为反应物的,我想请问一般在反应里面我们都是用三重态去考虑嘛?因为他们实验是用单重态做的,所以我有这个疑问,还是非常感谢各位老师 |
| B3LYP当闭壳层算单重态氧气,电子结构完全是错的 |
本帖最后由 喵星大佬 于 2022-3-24 07:15 编辑 wzkchem5 发表于 2022-3-24 02:08 用TDDFT应该是可以描述单重态氧气的 CAS的话最大也就(12,8),没啥算不了的,实际上只需要(8,6)甚至(6,4),还有就是用MC-PDFT比NEVPT2/MRCISD啥的便宜得多还有梯度,如果需要后续算反应也可以考虑,不过精度不太好说,而且OpenMolcas用起来也有点麻烦 |
| 参与人数Participants 1 | eV +4 | 收起 理由Reason |
|---|---|---|
|
| + 4 |
|
单重态氧气是多参考态体系,本来就不应该用DFT算,而应该用NEVPT2、MRCISD等多参考态方法算。用单参考态方法强行算多参考态体系,不仅难收敛,即使收敛了结果也不可信。 至于为什么B3LYP碰巧能收敛,应该是因为B3LYP的HF交换项成分大,所以HOMO-LUMO gap大,收敛还不至于那么困难。但是还是那句话,虽然B3LYP收敛了,结果也是不能用的,因为这是多参考态体系。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-18 09:10 , Processed in 0.184933 second(s), 26 queries , Gzip On.