计算化学公社

 找回密码 Forget password
 注册 Register

求助:用关键词guess=(mix,always)计算对称破缺态遇到了一些疑惑

查看数: 1145 | 评论数: 5 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2023-11-27 21:21

正文摘要:

本帖最后由 哇卡卡卡卡卡 于 2023-11-28 09:01 编辑 各位老师好!最近在进行对称破缺态的一些计算,根据文献中提到的TS-1的结构,根据帖子谈谈片段组合波函数与自旋极化单重态 - 思想家公社的门口:量子化学·分 ...

回复 Reply

zjxitcc 发表于 Post on 2023-12-2 20:42:00
哇卡卡卡卡卡 发表于 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章补充一下基础知识。
哇卡卡卡卡卡 发表于 Post on 2023-11-28 09:03:15
北大-陶豫 发表于 2023-11-28 08:59
对称破缺单重态写 ub3lyp/6-31g(d) guess=(mix,always),电荷多重度分别设为0和1,这没问题,但是你算的 ...

关键词在图片上,那这个自旋密度图还有自旋密度<S^^2>值,应该怎么理解呢?
北大-陶豫 发表于 Post on 2023-11-28 08:59:55
哇卡卡卡卡卡 发表于 2023-11-28 08:53
好的,谢谢老师的指导!所以如果我要研究对称破缺单重态,关键词应该这样写是吗(#p opt=(ts,calcfc,noei ...

对称破缺单重态写 ub3lyp/6-31g(d) guess=(mix,always),电荷多重度分别设为0和1,这没问题,但是你算的是氧气分子为什么要 opt=TS,这恐怕没有过渡态吧。
哇卡卡卡卡卡 发表于 Post on 2023-11-28 08:53:27
本帖最后由 哇卡卡卡卡卡 于 2023-11-28 09:03 编辑
zjxitcc 发表于 2023-11-27 21:52
guess=mix是专门给单重态UHF/UDFT计算用的(构造对称破缺初猜),不是给三重态计算用的。对于 非单重态,由 ...

好的,谢谢老师的指导!所以如果我要研究对称破缺单重态,关键词应该这样写是吗(#p opt freq ub3lyp/6-31g(d) guess=(mix,always),电荷多重度分别设为0和1),我前面有尝试去计算O2,上面同样的关键词,多重度分别设置为1和3,最后绘制出来的自旋密度图,这两者之间的差别我看不太懂,图我上传到附件了,请老师给予指导,谢谢!
zjxitcc 发表于 Post on 2023-11-27 21:52:14
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.

快速回复 返回顶部 返回列表 Return to list