ORCA5.0.3利用G16解析Hessian进行ESD计算

zhaixiaoyi001 · 发表于 Post on 2022-12-27 18:07:53

本帖最后由 zhaixiaoyi001 于 2025-12-16 09:40 编辑

1. 前言
      ORCA 5.0.3能够高效优化激发态结构，但使用ESD模块需要用到激发态Hessian。由于ORCA 5.0.3仍然不支持激发态解析Hessian，本文尝试使用能够计算激发态解析Hessian的Gaussian 16生成相应计算水平下的Hessian，让ORCA5.0.3读取该Hessian完成ESD模块计算。

2. 测试步骤
      1、首先使用ORCA 5.0.3优化苯分子S0和S1态结构并进行频率计算，得到优化后的结构：ben-s0-opt.out 和 ben-s1-opt.out，以及优化后的基态和激发态Hessian文件：ben-s0-opt.hess和 ben-s1-opt.hess。
      2、将 ben-s1-opt.out (S1结构)存为gjf文件 ben-s1-freq.gjf，使用Gaussian 16进行S1态频率计算。
      3、通过本文的python脚本将Gaussian 16 fchk文件中的Hessian转换为ORCA 5.0.3的.hess文件，得到 ben-s1-freq.hess。
      4、使用Gaussian 16 fchk文件转换得到的ben-s1-freq.hess与ORCA直接得到的ben-s1-opt.hess同时进行ESD荧光发射速率的计算，对比计算结果，观察结果可靠性。
注意：在结构优化过程中使用了RI加速，关键词见附件。ESD的计算过程完全参考ORCA 5.0.3手册。Gaussian 16直接使用freq TD=(nstates=5,root=1)关键词计算ORCA 5.0.3优化后S1结构的Hessian，使用了nosymm关键词。为了保证计算级别的一致性，ORCA 5.0.3计算使用 B3LYP/G 泛函，Gaussian 16计算使用B3LYP泛函，都使用了def2-SVP基组，未加色散校正。

3. 测试结果
      如图所示，使用Gaussian 16 fchk文件转化得到的S1态.hess进行ESD计算，与使用ORCA 5.0.3生成的S1态.hess文件计算得到的结果几乎一致，证明该方法较为可靠，通过Gaussian16进行激发态解析Hessian计算能够减少ORCA ESD模块的计算成本。

4. 脚本使用方法
      Python 环境 Python 3.7.6, numpy
      将fchk文件和Gaussian 16频率计算的gjf文件放入脚本目录，运行脚本，输入fchk文件名和gjf文件名，得到初步完成格式转化的.hess文件。将.hess文件复制到安装有ORCA环境的机器上，使用orca_vib xxx.hess 补全 Hessian文件内容，随后可以将该 Hessian 文件用于 ESD计算。

5. 注意事项
      本脚本只进行了少数测试，可能有bug，并且Hessian文件中的 $atom 项需要内置标准原子质量，本脚本中没有内置太多元素，使用脚本时需要进一步修改。本文涉及的输入输出文件打包放在了附件中。

zhaixiaoyi001 · 发表于 Post on 2023-1-10 09:00:19

最近测试了80原子体系，双路Xeon Platinum 8175M, PBE0-D3/def2-SVP计算水平。 ORCA计算S1 hessian用时约26 h，G16A 用时约1 h 30min，二者ESD模块计算结果基本一致。

chrinide · 发表于 Post on 2022-12-27 19:53:11

都在薅高斯的羊毛啊

话说高斯的羊毛薅起来真爽

Novice · 发表于 Post on 2022-12-28 09:26:33

这里我有个小小的疑问，G16直接使用freq关键词计算的S1结构的Hessian，和td + freq的Hessian是一样的吗？

冰释之川 · 发表于 Post on 2022-12-28 10:25:09

Novice 发表于 2022-12-28 09:26
这里我有个小小的疑问，G16直接使用freq关键词计算的S1结构的Hessian，和td + freq的Hessian是一样的吗？

前者是基于S1结构算该结构下基态的freq，后者是算激发态结构下激发态对应的freq

Novice · 发表于 Post on 2022-12-28 11:05:28

冰释之川发表于 2022-12-28 10:25
前者是基于S1结构算该结构下基态的freq，后者是算激发态结构下激发态对应的freq

嗯嗯，谢谢。所以我觉得楼主此处的步骤中应该采用的是td + freq。

zhaixiaoyi001 · 发表于 Post on 2022-12-28 11:14:49

Novice 发表于 2022-12-28 09:26
这里我有个小小的疑问，G16直接使用freq关键词计算的S1结构的Hessian，和td + freq的Hessian是一样的吗？

这里用的freq关键词也用了TD关键词，写的有点歧义，我做下修改

Novice · 发表于 Post on 2022-12-28 11:38:13

看了下，对于苯这么小的体系，ORCA的S1频率计算耗时17.343min，Gaussian不到5min，的确太有用了。不过话说等ORCA支持激发态解析频率了，那时候应该就不用这么麻烦了

flyingchow · 发表于 Post on 2023-1-19 22:54:10

请问群里有人测试过Turbomole吗？

Kalius · 发表于 Post on 2023-1-28 21:36:38

Orca做结构优化最好使用TightOPT+verytightSCF罢，要不然精度可能不够

RAL · 发表于 Post on 2023-1-28 21:56:36

Kalius 发表于 2023-1-28 21:36
Orca做结构优化最好使用TightOPT+verytightSCF罢，要不然精度可能不够

至少ORCA5的opt默认优化收敛限比Gaussian16严格，只是ORCA喜欢接近收敛限就停，用DefGrid3基本不会有小虚频的问题，没必要进一步提高收敛限

wzkchem5 · 发表于 Post on 2023-1-28 23:30:47

Kalius 发表于 2023-1-28 14:36
Orca做结构优化最好使用TightOPT+verytightSCF罢，要不然精度可能不够

默认收敛限和TightOPT+verytightSCF的差别基本上远小于格点误差了，所以影响精度不至于，如果遇到小虚频再把收敛限设严消虚频就行了

Kalius · 发表于 Post on 2023-1-29 17:30:12

wzkchem5 发表于 2023-1-28 23:30
默认收敛限和TightOPT+verytightSCF的差别基本上远小于格点误差了，所以影响精度不至于，如果遇到小虚频 ...

多谢，谨受教

tianmafei · 发表于 Post on 2023-5-7 16:00:49

老师好！这个脚本是 ./H_gTo.py这样运行吗？

zhaixiaoyi001 · 发表于 Post on 2023-5-18 09:39:50

tianmafei 发表于 2023-5-7 16:00
老师好！这个脚本是 ./H_gTo.py这样运行吗？

我一般是使用pycharm直接打开工程配置好运行库直接run的。
也可以 python3 H_gTo.py
现在看来这个脚本已经没什么用处了，ORCA光物理参数这块算的不太好，要算这些光物理参数可以用开源的fcclasses3

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] ORCA5.0.3利用G16解析Hessian进行ESD计算

评分 Rate

评分 Rate

浏览过的版块