哇卡卡卡卡卡 发表于 2023-11-28 09:03 三重态<S^2>理想值为2.0,你算出来2.0066,说明几乎没有自旋污染,算是自旋纯态,图表示两个O原子各有1个alpha单电子。 单重态<S^2>理想值为0,你算出来1.0033,说明自旋污染很大,有一个alpha未成对电子,有一个beta未成对电子,它们分布在两个不同的空间轨道上(可能你会好奇,那什么是“alpha和beta电子分布在同一个空间轨道上”?因为分布在一个空间轨道上=闭壳层单重态,与计算结果矛盾)。 如果你完全不理解我说的,建议阅读Szabo的《Modern Quantum Chemistry》一书前3章补充一下基础知识。 |
北大-陶豫 发表于 2023-11-28 08:59 关键词在图片上,那这个自旋密度图还有自旋密度<S^^2>值,应该怎么理解呢? |
哇卡卡卡卡卡 发表于 2023-11-28 08:53 对称破缺单重态写 ub3lyp/6-31g(d) guess=(mix,always),电荷多重度分别设为0和1,这没问题,但是你算的是氧气分子为什么要 opt=TS,这恐怕没有过渡态吧。 |
本帖最后由 哇卡卡卡卡卡 于 2023-11-28 09:03 编辑 zjxitcc 发表于 2023-11-27 21:52 好的,谢谢老师的指导!所以如果我要研究对称破缺单重态,关键词应该这样写是吗(#p opt freq ub3lyp/6-31g(d) guess=(mix,always),电荷多重度分别设为0和1),我前面有尝试去计算O2,上面同样的关键词,多重度分别设置为1和3,最后绘制出来的自旋密度图,这两者之间的差别我看不太懂,图我上传到附件了,请老师给予指导,谢谢! |
guess=mix是专门给单重态UHF/UDFT计算用的(构造对称破缺初猜),不是给三重态计算用的。对于 非单重态,由于alpha/beta电子数不相等,体系可以自然收敛到UHF/UDFT解(而非ROHF/RODFT解),因此无需guess=mix,当然always就更不需要了。如果你偏要写,那就是扰乱了初猜,但往往也能收敛到正确的解,所以你写不写最后结果都没差别,只是写了会暴露出自己理解不够。 对于②中的计算,本来关键词就固定了是闭壳层计算,显然不需要看“自旋密度图”,不用算也知道必然没有单电子,这是理论方法决定的。 如果你想研究的是对称破缺单重态,目前的文字描述中没有任何一个是。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-26 00:46 , Processed in 0.179243 second(s), 25 queries , Gzip On.