利用CP2K做了振动频率分析，之后Shermo计算热力学数据报错

桂の猿

老师好

我利用CP2K做了振动频率分析，之后，希望利用Shermo输出相应的热力学数据。

但是运行过程中发生报错，
报错信息如下：
severe(59): list-directed I/O syntax error, unit 10, file /...../1.out
因为文件类型限制，我将.out文件后缀改为了.dat，才可以上传。
如果老师要用这个文件运行Shermo，烦请改回.out

谢谢
桂

Stardust0831

输出文件里的hessian矩阵部分有很异常的"*"，可能是超出了cp2k程序能输出的位数。
请确保结构在计算振动分析前经过了结构优化，我质疑当前结构未被结构优化过的原因是初始结构的受力也很大，而不是接近零：

如果是孤立体系，请使用量子化学软件（如orca）来做结构优化+振动分析。


限制out文件格式的原因是因为文件上传前需要被压缩，见：《！！！！！计算化学公社新社员必读！！！！！》

关于帖子附件：发帖时点击“附件”按钮就可以上传附件。本论坛的服务器空间是很有限的，所以请大家珍惜使用。分享小文件为了方便可以用附件，大文件的话请尽量用网盘。如果是提问的时候需要发计算程序输出文件（或者其它各种文本型文件）给别人看，且文件大于500 KB，一定要压缩后再发，以节约论坛空间、减少下载耗时！如果压缩后还超过5MB，请用网盘。如果你对文件压缩毫无概念的话，认真看《计算化学研究者应当了解的文件压缩的常识》（http://sobereva.com/672）。

丁越

1. VIB|                        NORMAL MODES - CARTESIAN DISPLACEMENTS
   
2. VIB|
   
3. VIB|                         1                    2                    3
   
4. VIB|Frequency (cm^-1) ************          -648.356471          -522.899214
   
5. VIB|IR int (KM/Mole)   6684.348464             0.300352             9.253806
   
6. VIB|Red.Masses (a.u.)    12.087560             1.131195             1.064454
   
7. VIB|Frc consts (a.u.) ************            -0.000512            -0.000204

复制代码

很明显Shermo读取第一个频率时候出现了问题。

一定要对结构做了优化之后再计算频率。