计算化学公社

标题: 使用cp2k计算拉曼光谱时，峰强度太大，显示 ************ [打印本页]

作者
Author: liuzhanning09 时间: 2025-9-6 19:07
标题: 使用cp2k计算拉曼光谱时，峰强度太大，显示 ************
各位老师，我使用cp2k对MOF体系进行了Raman光谱计算，得到的结果导入到multiwfn中绘图时，报如下错误，请问该如何调整？

(, 下载次数 Times of downloads: 30)

作者
Author: DuAn111 时间: 2025-9-7 07:45
你的输入文件异常，拉曼活性数据1和4都太大了，需要你手动输入，建议检查你的计算结果输入文件。

作者
Author: sobereva 时间: 2025-9-7 15:14
别把Multiwfn极其清楚的英文提示当成乱码看待

作者
Author: liuzhanning09 时间: 2025-9-7 16:37

sobereva 发表于 2025-9-7 15:14
别把Multiwfn极其清楚的英文提示当成乱码看待

不好意思社长，没太理解

，您是指我截图的地方Multiwfn提示让我手动输入拉曼强度吗？我用文本器打开out文件中里边显示的也是*******，所以不知道是什么原因？请问是我的inp文件命令有错误吗？

作者
Author: sobereva 时间: 2025-9-7 22:16

liuzhanning09 发表于 2025-9-7 16:37
不好意思社长，没太理解，您是指我截图的地方Multiwfn提示让我手动输入拉曼强度吗？我用文本器打开ou ...

要么自己随便猜一个数值，要么自己改CP2K源代码，修改输出格式，使得能够输出很大的值。数值超过当前格式能显示的上限就会显示为星号

作者
Author: liuzhanning09 时间: 2025-9-8 08:10

sobereva 发表于 2025-9-7 22:16
要么自己随便猜一个数值，要么自己改CP2K源代码，修改输出格式，使得能够输出很大的值。数值超过当前格式 ...

了解了，感谢社长

作者
Author: sjyy 时间: 2026-2-3 12:02
请问这个解决了吗，我也遇到了这个问题

作者
Author: liuzhanning09 时间: 2026-3-7 16:36
没有解决。放弃了。。用的MS算的。

作者
Author: Uus/pMeC6H4-/キ 时间: 2026-5-3 01:01

sobereva 发表于 2025-9-7 22:16
要么自己随便猜一个数值，要么自己改CP2K源代码，修改输出格式，使得能够输出很大的值。数值超过当前格式 ...

社长，开发版的CP2K把这部分的输出格式从3(1X,F12.6,8X)改成3(1X,ES18.10E3,2X)了，这应该不会影响multiwfn的读取吧？在(A^4/amu)单位下，超过原先F12.6输出格式的拉曼强度数值得至少十几万了，此类结果如何确认是合理的，而不是如几何优化/振动分析的收敛限不够严或原子冻结不当等因素造成的异常呢？

作者
Author: sobereva 时间: 2026-5-3 08:03

Uus/pMeC6H4-/キ发表于 2026-5-3 01:01
社长，开发版的CP2K把这部分的输出格式从3(1X,F12.6,8X)改成3(1X,ES18.10E3,2X)了，这应该不会影响multiw ...

Multiwfn目前是按自由格式读取的，只要CP2K数值输出的格式宽度不变（看起来都是21个字符），就没问题

红外活性正比于极化率对正则坐标导数的平方，并且体系gap越小极化率倾向越低，如果凭借这些物理原因看起来大得有理有据，并且收敛精度看起来已经足够充分、收敛过程顺畅，就不至于有问题

作者
Author: sobereva 时间: 2026-5-3 12:40

Uus/pMeC6H4-/キ发表于 2026-5-3 01:01
社长，开发版的CP2K把这部分的输出格式从3(1X,F12.6,8X)改成3(1X,ES18.10E3,2X)了，这应该不会影响multiw ...

强调一下，是拉曼活性Raman activity（而不是拉曼强度Raman intensity）正比于极化率对正则坐标导数的平方。后者需要由前者进一步转化，看下文的相关说明
使用Multiwfn绘制红外、拉曼、UV-Vis、ECD、VCD和ROA光谱图
http://sobereva.com/224

北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）这页幻灯片也说了
(, 下载次数 Times of downloads: 23)

欢迎光临计算化学公社 (http://bbs.keinsci.com/)