计算化学公社

标题: 求助ORCA算双自由基的方法 [打印本页]

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



作者
Author:
zjxitcc    时间: 2023-10-26 13:38
本帖最后由 zjxitcc 于 2023-10-26 13:47 编辑

对称破缺单重态各种计算方式:
(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的参考态是三重态。

作者
Author:
ljc050512    时间: 2023-10-26 17:23
zjxitcc 发表于 2023-10-26 13:38
对称破缺单重态各种计算方式:
(1)在ORCA中使用BrokenSym或FlipSpin关键词,例子见《双金属中心化合物的 ...

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





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3