关于在CP2K使用Metadynamics计算能量的不合理问题

Wguan

使用Metadynamics计算一个反应的能垒，CV设置的是其中的一个关键键长的变化，输入文件已经上传，最终得出的结果画图也已上传，但是能量的变化和文献以及使用高斯算出的能量变化（能垒7 kcal/mol）不太相符，主要有以下几个问题：1.是不是因为参数设置问题，导致计算结果出错；2.是否是因为没有联用Plumed导致整个模拟过程出错（用VMD看了整个过程，符合要求的反应）；3.使用的是graph.psmp -ndim 1 -ndw 1  -file MTD-10000.restart -cp2k 输出的CV与能量的输出文件，是否不太正确。

ljh123

不知道体系没人能给出很好的建议。
文献应该发一下，因为很多文献算能垒都是乱算的，不一定准。
建议用plumed，参数管理也更清晰些。另外你SCF部分参数写的也不合适，可以用Multiwfn转换一个替换上去

Wguan

ljh123 发表于 2025-12-8 13:45
不知道体系没人能给出很好的建议。
文献应该发一下，因为很多文献算能垒都是乱算的，不一定准。
建议用pl ...

已将我的体系上传
因为我现在也是在试验阶段，所以只是先用三分子来试试能不能出现恰当的结果，根据我自己使用高斯算出来的能垒是在7kcal/mol，而CP2K出来的却是一个无势垒过程，这就是我一直疑惑的地方。