计算化学公社

标题: 求助对开壳层体系如何计算不同电子态之间的自旋-轨道耦合常数 [打印本页]

作者
Author: 13672068720 时间: 2026-7-2 16:41
标题: 求助对开壳层体系如何计算不同电子态之间的自旋-轨道耦合常数
计算自旋-轨道耦合矩阵元时，对开壳层体系，例如自由基-生色团（二重态-四重态间的系间窜越）或双自由基（单重态-三重态间的系间窜越），存在自旋污染，电子态无严格意义的自旋多重度，如何计算不同电子态之间的自旋-轨道耦合常数

作者
Author: zjxitcc 时间: 2026-7-2 17:11
本帖最后由 zjxitcc 于 2026-7-3 08:25 编辑

以下建议假设你已经优化好了几何结构，只考虑一个静态结构下的SOC计算：
（1）通用，说服力强，但计算量较大。在NEVPT2下计算旋轨耦合，用ORCA/Molpro都可以，自己跟着ORCA手册学。注意，虽然在CASSCF下也可以计算旋轨耦合，但没考虑动态相关，精度可能一般；万一碰到多参考领域审稿人，可能会遭受质疑。

（2）有局限性，但经济实惠。目前只适用于Sn-Tn，不适用于二重态-四重态。先用Gaussian/PySCF进行三重态ROKS计算，然后使用MOKIT的fch2inp小程序，例如运行

fch2inp O2_triplet.fch -mrsf

复制代码

产生MRSF-TDA输入文件（包含坐标，轨道系数，基组数据，MRSF关键词，计算精度控制），绝大部分内容已写好。如果直接使用免费量化软件GAMESS对该输入文件做计算，可以在MRSF-TDA下获得Sn和Tn能量，基态可以是开壳层单重态或三重态，该方法都没问题。但MRSF-SOC尚未合并至GAMESS主分支，需要找韩国庆北国立大学的Cheol Ho Choi 索要他们组里的GAMESS。然后修改inp文件中$TDDFT那一行为如下两行
$TDDFT NSTATE=5 NRAD=99 NLEB=590 MRSOC=.T. $END
$TRANST $END
注意$符号之前必须有且仅有一个空格，不可以顶格写。

（3）尝试ORCA的ROCIS/DFT方法，手册里有详细示例。

（4）这里夹带私货：当然也可以找我合作:(

如果你还没找到合适的方法优化几何结构，则应当先解决结构问题。

作者
Author: sobereva 时间: 2026-7-3 00:53
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “开壳层体系计算SOC” 改了，以后务必注意！

作者
Author: sobereva 时间: 2026-7-3 00:54
ORCA的ROCIS/DFT对开壳层体系可以计算出没有自旋污染的激发态，并计算SOC

北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）里专有一部分讲了具体用法
(, 下载次数 Times of downloads: 3)

对于参考态是自旋极化单重态的情况也有办法
(, 下载次数 Times of downloads: 5)

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