dantevinsky 发表于 2023-3-3 17:59
有可能是数值误差的原因，也有可能是没有投影掉平动转动的原因，具体得看ASE package代码中间有没有这一步 ...

我去看了下代码，应该是ASE没有去掉平动转动项。

ahxb 发表于 2023-3-3 19:11
在Gaussian输出文件中给出频率及振动模式之前的部分可以看到如下信息：

而计算出的频率却是：

谢谢！我又去算了几个体系，确实都是少了平动和转动项。

在Gaussian输出文件中给出频率及振动模式之前的部分可以看到如下信息：

1. Full mass-weighted force constant matrix:
   
2. Low frequencies --- -576.0604  -13.2694   -5.9568   -3.0655    0.0009    0.0012
   
3. Low frequencies ---    0.0013  186.9810  331.0460

复制代码

而计算出的频率却是：

1. Frequencies --   -576.0603               186.9797               331.0456
   
2. Frequencies --    348.3210               396.3831               454.3103

复制代码

因此对应于平动和转动的6个自由度被忽略了。具体是用什么算法忽略的我不清楚。

有可能是数值误差的原因，也有可能是没有投影掉平动转动的原因，具体得看ASE package代码中间有没有这一步。