求助：用CP2K计算液体水的红外光谱特征峰显示不全

THz

用CP2K跑分子动力学计算了64个水分子，时间步0.5 fs，总时长20 ps，然后对输出的偶极矩做自相关后再傅里叶变换，得到的频谱缺少bend模式（1800 cm^-1左右）的特征峰，而且整个谱线看起来噪声比较多。请教一下是什么原因。（本人感觉或许20 ps有点短，输入文件中相关参数设置的精度不够）

CP2K的输入文件和频谱后处理的MATLAB脚本放在附件里了。