关于单线态氧原子的计算关键词的选择

Mr.zhen

想要计算某一有机分子与单线态氧原子（O1D）的过渡态，考虑到O1D自旋多重度也是1，是不是可以不用td关键词，写成：#p opt=(calcfc,ts,noeigen) freq m062x/6-31g(d) scrf ，然后自旋多重度设为1进行计算；还是要写成：#p td opt=(calcfc,ts,noeigen) freq m062x/6-31g(d) scrf ，然后自旋多重度设为3呢？之前试过把自旋多重度设为1，但没能找到满意的过渡态，使用td的话耗时极多，在等待期间希望大家能指点一二

wzkchem5

单线态氧原子用单参考方法几乎不可能定性正确地描述（即使用SF-TDDFT，也只能正确描述自旋简并性而不能正确描述空间简并性）。但鉴于你说TDDFT已经耗时极多了，用多参考态方法估计也跑不动，所以我的建议是用ORCA或Q-Chem，以三线态为参考态做SF-TDA计算

Mr.zhen

wzkchem5 发表于 2023-4-20 16:03
单线态氧原子用单参考方法几乎不可能定性正确地描述（即使用SF-TDDFT，也只能正确描述自旋简并性而不能正确 ...

谢谢老师，但是因为手头只有Gaussian可用，并且不加td的时候算得还比较快，所以可以用#p opt=(calcfc,ts,noeigen) freq m062x/6-31g(d) scrf ，然后自旋多重度设为1，“凑合”用吗，只是不知道这样写关键词行不行

wzkchem5

Mr.zhen 发表于 2023-4-20 09:25
谢谢老师，但是因为手头只有Gaussian可用，并且不加td的时候算得还比较快，所以可以用#p opt=(calcfc,ts, ...

不能，这样得到的电子结构是定性错误的。
如果只买了高斯的话，倒可以算是不用Q-Chem的理由，但是不能作为不用ORCA的理由，毕竟ORCA对学术界是免费的，想下载随时注册账号下载，而且也不难上手

wjc404

求助时可以描述得具体一点，比如要和氧原子反应的有机分子长啥样，反应位点在哪里。