计算化学公社

标题: 为什么ORCA计算得到的Hessian矩阵不是对称的 [打印本页]

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

作者
Author: 王二葛 时间: 2024-7-15 17:22
数值问题，通常对称化之后再使用

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

复制代码

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

作者
Author: lucky1999 时间: 2024-7-15 18:09

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

收到，谢谢老师的解答。

作者
Author: lucky1999 时间: 2024-7-15 18:09

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

收到，谢谢老师的解答。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)