计算化学公社

标题: 请问CP2K进行VIBRATIONAL_ANALYSIS计算的原理 [打印本页]

作者
Author: wqqxxy 时间: 2022-11-25 21:10
标题: 请问CP2K进行VIBRATIONAL_ANALYSIS计算的原理
本人半导体工艺新手一枚，想请问一下站内老师们，可否详细说一下CP2K进行VIBRATIONAL_ANALYSIS进行计算时的原理。目前仅了解它是通过对体系内每个原子向+x, -x, +y, -y, +z, -z 6个方向分别进行移动，用数值的方法得到能量的二阶导即力常数的一个Hessian矩阵，然后计算频率。它是怎样通过移动原子位置求得力常数矩阵。力常数矩阵又和原子的振动频率又有什么关系，是如何通过力常数求得原子的振动频率的。提前感谢各位老师的解答。

作者
Author: sobereva 时间: 2022-11-25 22:54
Hessian矩阵就是能量对核坐标的二阶导数矩阵，去Google了解有限差分（finite difference）算导数的数值方法。CP2K里是对解析梯度做一阶有限差分来得到半数值Hessian的。

Hessian转化到振动频率的原理看
基于fch中的Hessian矩阵计算振动频率的简单程序Hess2freq
http://sobereva.com/328（http://bbs.keinsci.com/thread-3368-1-1.html）

作者
Author: wqqxxy 时间: 2022-11-26 09:40
十分感谢sob老师

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