计算化学公社

标题: 求助cp2k计算拉曼活性只给出了星号 [打印本页]

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

作者
Author: Uus/pMeC6H4-/キ 时间: 5 hour ago
本帖最后由 Uus/pMeC6H4-/キ于 2026-4-27 18:21 编辑

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

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

作者
Author: Richael 时间: 4 hour ago

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

但是，为了确 ...

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

作者
Author: Uus/pMeC6H4-/キ 时间: 3 hour ago

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上折腾不太值得。

作者
Author: Richael 时间: 2 hour ago

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

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

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