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模块计算结果基本一致。

zhaixiaoyi001 · 发表于 Post on 2025-12-16 09:17:16

陈AG 发表于 2025-12-13 16:32
lz，请问我gjf都是geom=read的，里面没有坐标信息，是不是需要修改一下代码，从fch里读元素呀

感觉也可以不改，把计算的log文件重新保存一下就有坐标信息了。

陈AG · 发表于 Post on 2025-12-13 16:32:12

lz，请问我gjf都是geom=read的，里面没有坐标信息，是不是需要修改一下代码，从fch里读元素呀

wzkchem5 · 发表于 Post on 2024-11-7 16:42:59

qaqfdmmj 发表于 2024-11-7 03:48
老师，我这里还是想确认一下是否写对了……我是计算振动分辨的荧光光谱的，已经提前对S0基态进行了opt+fr ...

我感觉没问题，但是nroots不用设那么大，一般如果优化第一激发态的话，nroots=3足够了

qaqfdmmj · 发表于 Post on 2024-11-7 10:48:27

wzkchem5 发表于 2024-6-28 15:41
1. 是基态的
2. 应该用.spectrum文件作图。out文件不包含绘制振动分辨光谱所需的信息，因此不管是avogad ...

老师，我这里还是想确认一下是否写对了……我是计算振动分辨的荧光光谱的，已经提前对S0基态进行了opt+freq的计算，然后代入手册中AHAS的写法如下：
! PBE0 def2-TZVP def2/J def2-TZVP/C RIJCOSX tightSCF noautostart miniprint nopop ESD(FLUOR)
%maxcore    4000
%pal nprocs  48 end
%TDDFT
NROOTS 10
IROOT 1
END
%ESD
GSHESSIAN "6_gs.hess"
DOHT TRUE
HESSFLAG AHAS
END
* XYZFILE 0 1  6_gs.xyz
其中6_gs对应的都是S0的结构和hessian，想问一下老师：这样得到的spectrum文件是荧光光谱数据吗？

zhaixiaoyi001 · 发表于 Post on 2024-7-25 11:26:48

heroooo 发表于 2024-7-23 16:43
请问把高斯计算的.fchk转成.hess文件的原理是什么呢？不同的字段是怎么对应的呢，有没有坐标单位的转换？

就是把fchk文件中的hessian提取出来，转换成了ORCA .hess文件输出的格式而已，单位什么的都处理好了。

heroooo · 发表于 Post on 2024-7-23 16:43:25

zhaixiaoyi001 发表于 2023-1-10 09:00
最近测试了80原子体系，双路Xeon Platinum 8175M, PBE0-D3/def2-SVP计算水平。 ORCA计算S1 hessian用时约26 ...

请问把高斯计算的.fchk转成.hess文件的原理是什么呢？不同的字段是怎么对应的呢，有没有坐标单位的转换？

wzkchem5 · 发表于 Post on 2024-6-28 15:41:42

qaqfdmmj 发表于 2024-6-28 07:36
是的老师，在手册里有VG和AHAS的示例，不过是在吸收光谱的章节中。
q1：如果是计算荧光光谱这里esd板块 ...

1. 是基态的
2. 应该用.spectrum文件作图。out文件不包含绘制振动分辨光谱所需的信息，因此不管是avogadro还是其他软件，通过读取out文件作的图应该都是非振动分辨的谱

qaqfdmmj · 发表于 Post on 2024-6-28 14:36:11

wzkchem5 发表于 2024-6-27 20:11
对，但这是AH的写法。VG或AHAS只需要基态结构和基态Hessian，手册另有其他示例

是的老师，在手册里有VG和AHAS的示例，不过是在吸收光谱的章节中。
q1：如果是计算荧光光谱这里esd板块中需要的结构和hessian是基态还是激发态的呢？
q2：绘制振动分辨荧光光谱时，请问我用官方推荐的avogrado打开out文件->extensions->spectra还是将.spectrum文件中的数据导入origin自己作图呢？因为我看到avogadro中绘制的是uv和cd，不太明白是不是所需要的荧光光谱

wzkchem5 · 发表于 Post on 2024-6-27 20:11:28

qaqfdmmj 发表于 2024-6-27 10:52
wzk老师你好，我之前是用gaussian16算振动分辨荧光但是显示 Total convergence = 0.0%.查阅了论坛发现您 ...

对，但这是AH的写法。VG或AHAS只需要基态结构和基态Hessian，手册另有其他示例

qaqfdmmj · 发表于 Post on 2024-6-27 17:52:23

wzkchem5 发表于 2024-6-23 02:36
结构优化和频率计算的软件应当统一。如果频率打算用高斯算，那么结构可以用orca预优化，但最终要过一遍高 ...

wzk老师你好，我之前是用gaussian16算振动分辨荧光但是显示 Total convergence = 0.0%.查阅了论坛发现您曾在别的贴子说过这是因为基态和激发态结构差异较大所致，因此推荐使用orca的esd的vg或ahas方法进行计算。我翻阅了orca5.0.4的manual，请问计算振动分辨荧光前是先得到S0和S1的hessian，然后再输入以下的文件吗（manual中的示例）？
ps：这里如果要使用ahas而不是默认的vg是不是就在esd板块加上hessflag ahas即可？

wzkchem5 · 发表于 Post on 2024-6-23 02:36:04

qaqfdmmj 发表于 2024-6-22 04:07
你好，我最近刚刚接触orca（算振动分辨），想问一下这个方法利用gaussian计算得到ben-s1-freq.hess，是不是 ...

结构优化和频率计算的软件应当统一。如果频率打算用高斯算，那么结构可以用orca预优化，但最终要过一遍高斯的结构优化。

qaqfdmmj · 发表于 Post on 2024-6-22 11:07:21

你好，我最近刚刚接触orca（算振动分辨），想问一下这个方法利用gaussian计算得到ben-s1-freq.hess，是不是说就不用orca进行频率计算而仅进行结构优化即可？

zhaixiaoyi001 · 发表于 Post on 2023-6-25 18:32:01

dacaoyang2008 发表于 2023-6-24 08:40
请问为什么orca_vib t1_opt.hess时总是提示Error-Message:error while reading $act_energy

可能是文件内容有问题，你把所有需要的文件发上来看看

dacaoyang2008 · 发表于 Post on 2023-6-24 08:40:59

请问为什么orca_vib t1_opt.hess时总是提示Error-Message:error while reading $act_energy

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

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

评分 Rate

评分 Rate

评分 Rate

浏览过的版块