求助cp2k计算拉曼活性只给出了星号

Richael

求助熟悉cp2k计算的大佬，我在利用cp2k计算拉曼活性时，最后输出只给出了*******，无法得到拉曼信息，请问该如何解决？这里给出了我的输入文件，希望有大佬帮忙指出问题。

Uus/pMeC6H4-/キ

我有一个修改源代码重新编译的解决方法，预计下个CP2K版本可以生效，也可以在这里说一说。

但是，为了确保这个计算首先有意义，我必须问一问：你用的CP2K是什么版本，输入文件是怎么生成的，为什么振动分析/计算拉曼的任务会&MD字段（还是Langevin系综下2 fs/步跑300万步这么奇怪的设定），./GST.qbs基组文件和类似TZVP-GTH-q4-very-confined这样的基组是又是哪里来的，为什么PERIODIC_DIPOLE_OPERATOR设为了F？

Richael

Uus/pMeC6H4-/キ 发表于 2026-4-27 18:15
我有一个修改源代码重新编译的解决方法，预计下个CP2K版本可以生效，也可以在这里说一说。

但是，为了确 ...

是2023.1版本，输入文件是在课题组之前用来跑MD的输入文件基础上参考了Multiwfn生成的拉曼计算指令。这个基组文件是课题组内部使用的，用了快十年了，是稳定精确的。PERIODIC_DIPOLE_OPERATOR我也不太清楚更多的是Multiwfn生成的拉曼计算指令默认项。

Uus/pMeC6H4-/キ

Richael 发表于 2026-4-27 19:11
是2023.1版本，输入文件是在课题组之前用来跑MD的输入文件基础上参考了Multiwfn生成的拉曼计算指令。这个 ...

PERIODIC_DIPOLE_OPERATOR设T表示周期性体系、F表示孤立体系，T是选项省略不写时的默认以及multiwfn生成输入文件时的默认，所以设F肯定是有意改过的，我不知道这是否符合你的体系。2023.1已经是三年前的老版本了，能更新还是尽量更新，毕竟有好多故障修复和特性增添了。

对用&VIBRATIONAL_ANALYSIS/&PRINT/&MOLDEN_VIB输出振动信息molden文件的情形，去molden找到这个振动模式的拉曼强度（红外同理）然后把输出文件的星号替换掉即可。要不然就得去源代码里找输出这部分信息的处理逻辑、修改相应的fortran格式，然后再编译CP2K了，但是我觉得在这么老版本的CP2K上折腾不太值得。

Richael

Uus/pMeC6H4-/キ 发表于 2026-4-27 20:01
PERIODIC_DIPOLE_OPERATOR设T表示周期性体系、F表示孤立体系，T是选项省略不写时的默认以及multiwfn生成 ...

感谢您的回答，我先进行一些尝试

Richael

Uus/pMeC6H4-/キ 发表于 2026-4-27 20:01
PERIODIC_DIPOLE_OPERATOR设T表示周期性体系、F表示孤立体系，T是选项省略不写时的默认以及multiwfn生成 ...

想请问怎么z在molden下找到拉曼强度，只找到了红外强度

Uus/pMeC6H4-/キ

Richael 发表于 2026-4-29 11:45
想请问怎么z在molden下找到拉曼强度，只找到了红外强度

哦我记错了，目前输出molden是不包含拉曼强度的（这个格式的原始定义似乎只有红外强度）。要么就等2026.2新版本发布，要么就修改2026.1的源代码来重新编译了，控制振动分析信息打印的部分应该在src/motion/vibrational_analysis.F的SUBROUTINE vib_out里面。