各位老师们好，最近在学习用python去做normal mode采样，目前的思路是调用ASE第三方库中的Vibrations函数去计算分子的振动模式，然后沿着每个振动模式的方向进行位移。在验证计算结果是否准确的时候，选取了一个分子结构，在HF/Def2-SVP的计算水平下，用高斯和代码都算了一遍，发现由于ase进行振动分析时采用的是数值差分的方法计算hessian，即对体系中的每个原子分别进行各个方向0.01埃的移动，计算原子力的变化，所以在数值上与高斯的计算结果有所差异，同时ase输出的振动模式中，还存在一些很小的实频和虚频（一般在100i-100cm^-1），我不太清楚这些多出来的频率是由于数值计算所带来的误差，还是高斯通过某种规则忽略了这些频率，希望有老师可以为我解惑。 (, 下载次数 Times of downloads: 6)