为什么ORCA计算得到的Hessian矩阵不是对称的

lucky1999 · 发表于 Post on 2024-7-15 16:51:16

各位老师，目前我在学习sob老师写的基于fch中的Hessian矩阵计算振动频率的简单程序Hess2freq(http://sobereva.com/328)。我想将ORCA输出的Hessian矩阵采用类似的方法求得谐振频率和正则坐标。我采用ORCA计算了一个水分子的Hessian矩阵，请问ORCA计算得到的Hessian矩阵为什么不是完全对称的呢？附件中为ORCA的输入文件与.hess文件。谢谢各位老师的解答。

王二葛 · 发表于 Post on 2024-7-15 17:22:07

数值问题，通常对称化之后再使用

H = (H + H.T) / 2.

复制代码

wdcqx · 发表于 Post on 2024-7-15 17:41:00

因为在构建交换相关项的静态二阶导数时仅包含了格点权重的一阶导数，导致了Hessian的对称性破缺。可以参考Frank Neese的这篇综述（10.1080/00268976.2015.1025114）。

lucky1999 · 发表于 Post on 2024-7-15 18:09:18

wdcqx 发表于 2024-7-15 17:41
因为在构建交换相关项的静态二阶导数时仅包含了格点权重的一阶导数，导致了Hessian的对称性破缺。可以参考F ...

收到，谢谢老师的解答。

lucky1999 · 发表于 Post on 2024-7-15 18:09:52

王二葛发表于 2024-7-15 17:22
数值问题，通常对称化之后再使用

收到，谢谢老师的解答。

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

[ORCA] 为什么ORCA计算得到的Hessian矩阵不是对称的