求助ORCA算双自由基的方法

ljc050512

各位老师好，最近在研究计算双自由基体系。已经拜读了sob老师的《谈谈片段组合波函数与自旋极化单重态 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)》，里面是关于高斯软件计算双自由基体系的方法。根据帖子中的方法1，混合HOMO与LUMO，对benzyne进行了计算练习。现在想用免费且强大的ORCA进行计算。ORCA5.0.3手册795页的SF-TDA可以计算benzyne。
我的问题是，ORCA中能不能实现sob老师帖子中提到的针对高斯的方法1（混合HOMO与LUMO）或者方法(3)利用stable=opt关键词，输入文件应该怎么写呢？

zjxitcc

对称破缺单重态各种计算方式：
（1）在ORCA中使用BrokenSym或FlipSpin关键词，例子见《双金属中心化合物的磁交换耦合常数的计算》，《量子化学软件基础》习题 （3）》
（2）从别的程序产生对称破缺初猜、甚至获得收敛、稳定的UHF波函数，然后传给ORCA，例子见《利用MOKIT从Gaussian向其他量化程序传轨道》，《利用MOKIT从PySCF向其他量化程序传轨道》

检验波函数稳定性
%scf
  STABPerform true
  STABRestartUHFifUnstable true
  STABMaxIter 500
  STABDTol 1e-5
  STABRTol 1e-5
end

另外，这些只是对称破缺单重态UHF/UDFT的计算技巧，与SF-TDA关系不大，SF-TDA的参考态是三重态。

ljc050512

zjxitcc 发表于 2023-10-26 13:38
对称破缺单重态各种计算方式：
（1）在ORCA中使用BrokenSym或FlipSpin关键词，例子见《双金属中心化合物的 ...

非常感谢您的指点！参考您推荐的博文进行了学习和练习。