请问CP2K进行VIBRATIONAL_ANALYSIS计算的原理

wqqxxy · 发表于 Post on 2022-11-25 21:10:51

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

sobereva · 发表于 Post on 2022-11-25 22:54:56

Hessian矩阵就是能量对核坐标的二阶导数矩阵，去Google了解有限差分（finite difference）算导数的数值方法。CP2K里是对解析梯度做一阶有限差分来得到半数值Hessian的。

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

wqqxxy · 发表于 Post on 2022-11-26 09:40:34

十分感谢sob老师

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

[CP2K] 请问CP2K进行VIBRATIONAL_ANALYSIS计算的原理