开壳层体系计算SOC

13672068720

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

zjxitcc

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

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

1. 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方法，手册里有详细示例。

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

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