计算化学公社

 找回密码 Forget password
 注册 Register

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

查看数: 1436 | 评论数: 4 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2024-7-15 16:51

正文摘要:

各位老师,目前我在学习sob老师写的基于fch中的Hessian矩阵计算振动频率的简单程序Hess2freq(http://sobereva.com/328)。我想将ORCA输出的Hessian矩阵采用类似的方法求得谐振频率和正则坐标。我采用ORCA计算了一个水 ...

回复 Reply

lucky1999 发表于 Post on 2024-7-15 18:09:52
王二葛 发表于 2024-7-15 17:22
数值问题,通常对称化之后再使用

收到,谢谢老师的解答。
lucky1999 发表于 Post on 2024-7-15 18:09:18
wdcqx 发表于 2024-7-15 17:41
因为在构建交换相关项的静态二阶导数时仅包含了格点权重的一阶导数,导致了Hessian的对称性破缺。可以参考F ...

收到,谢谢老师的解答。
wdcqx 发表于 Post on 2024-7-15 17:41:00
因为在构建交换相关项的静态二阶导数时仅包含了格点权重的一阶导数,导致了Hessian的对称性破缺。可以参考Frank Neese的这篇综述(10.1080/00268976.2015.1025114)。
王二葛 发表于 Post on 2024-7-15 17:22:07
数值问题,通常对称化之后再使用
  1. H = (H + H.T) / 2.
复制代码

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 03:11 , Processed in 0.170138 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list