计算化学公社

标题: 计算频率相关的G矩阵和F矩阵怎么获得？ [打印本页]

作者
Author: 找镜 时间: 2021-10-23 23:31
标题: 计算频率相关的G矩阵和F矩阵怎么获得？
各位老师好！
最近一直在研究关于力常数相关的东西，最初是借助Compliance计算了柔性力常数。
后来有老师给我提到了G矩阵和F矩阵之类的东西，但是我理论学习得不好，一直没有搞明白是什么东西。

虽然也看过了关于卢老师对柔性力常数和Hess2Freq的相关介绍，但是对原理还是一知半解，这些地方似乎也没有提到G矩阵和F矩阵。

请问各位老师，G矩阵和F矩阵的原理是什么？
怎么才能基于Gaussian中的Hessian矩阵计算得到？
是否有相关的软件可以直接输出呢？

如果老师们有空，是否可以顺便指点一下:
Compliance中的逆矩阵能否基于Gaussian中的Hessian矩阵获得？
比如三原子的9×9的Hessian矩阵，通过质权得到力常数矩阵，然后有办法直接通过对角元之类的看到键、键角的力常数吗？
如果不可以，笛卡尔坐标形式的Hessian矩阵，可以转换成内坐标下的吗？

作者
Author: dantevinsky 时间: 2022-7-13 11:39
推荐看Wilson的Molecular Vibrations: The Theory of Infrared and Raman Vibrational Spectra；
其实就是将笛卡尔坐标系下的Hessian转换至内坐标下求解。首先构造将笛卡尔坐标转换至内坐标的B矩阵，再之后根据GF method就能得到相关的简谐频率，这跟高斯里面基于笛卡尔坐标系下算出来的简谐频率结果是一样的。(Sob老师的Hess2Freq省略了笛卡尔坐标系中投影平动与转动的过程，但为了保证最终简谐频率中低频的准确，这种投影是必要的)
原理上面已经说了。要想得到你得自己先定义一套内坐标，然后构造B矩阵，然后计算G与F矩阵。
相关软件的话，我不清楚，我是自己手写的轮子。。。。

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