计算化学公社

标题: 利用CP2K做了振动频率分析，之后Shermo计算热力学数据报错 [打印本页]
作者
Author: 桂の猿    时间: 2025-4-3 17:45
标题: 利用CP2K做了振动频率分析，之后Shermo计算热力学数据报错
老师好

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

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

谢谢
桂
作者
Author: Stardust0831    时间: 2025-4-3 18:24
本帖最后由 Stardust0831 于 2025-4-3 18:44 编辑

输出文件里的hessian矩阵部分有很异常的"*"，可能是超出了cp2k程序能输出的位数。 (, 下载次数 Times of downloads: 5)
请确保结构在计算振动分析前经过了结构优化，我质疑当前结构未被结构优化过的原因是初始结构的受力也很大，而不是接近零：
(, 下载次数 Times of downloads: 5)
如果是孤立体系，请使用量子化学软件（如orca）来做结构优化+振动分析。


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



作者
Author: 丁越    时间: 2025-4-3 20:27
  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读取第一个频率时候出现了问题。

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



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3