CP2K和VASP进行CI-NEB计算结果差别，以及如何选择在CP2K的中的插点数量

Ying_Zhang

想请各位老师给点建议，因为接触CP2K有一段时间，所以为了后续跑大模型选择CP2K，之前一直跑VASP。所以，最近拿同一个模型在CP2K和VASP计算，检验一下我CP2K掌握的能力，即计算Mg和Bi金属化合物的空位过渡态计算（CI-NEB）。发现两个问题：

1）VASP可以根据给定的初末结构，通过VTST脚本，确定插入的过渡态的数量，但是CP2K如何找合适的NEB数量？因为我跑CP2K的CI-NEB，发现不同的数量，很大程度影响计算，插入较少和过多，计算会卡在某一步，一直循环，但是一直尝试又较浪费时间。有什么经验和方法，减少时间确定较为适合的NEB数量。

2）VASP和CP2K结果出现较大差别，VASP结果过渡态为0.773607 eV，CP2K为1.382307142 eV，采用相同的surpercell的结构且K-point相同，CP2K采用DZVP-MOLOPT-SR-GTH级别计算。VASP结果更接近文献，当然文献也是VASP结果。所以，这差别是什么造成的？还是我CP2K计算过渡态仍存在问题？

3）附上我计算的NEB.inp和NEB.out结果以及VASP的相关结果，麻烦各位老师看看，有什么问题，哪里可以修改。

nuonuo123456

1.可以用VTST脚本确定数量，CP2K和vasp用一样的NEB数量
2.vasp我个人认为计算过渡态是更稳定和精准的，当然如何CP2K输入文件合理的话应该差异也不大，DZVP计算能量基组有点小了，可以提到TZV2P再试试

Ying_Zhang

nuonuo123456 发表于 2026-3-12 19:59
1.可以用VTST脚本确定数量，CP2K和vasp用一样的NEB数量
2.vasp我个人认为计算过渡态是更稳定和精准的，当 ...

感谢建议
1.第一次跑CP2K过渡态，我就是拿vtst决定插入数量，结果每个点一直停留在Step1，一直持续计算能量和受力，不进入到Step2。发了帖子询问，后面加上自己尝试，发现提升NEB插点数量，可以使CI-NEB跑起来，只不过也存在，有的跑起来中间又停止在某一步，不继续进行，有的可以顺利收敛。
2. 是的我也觉得VASP是更准确的，但是我后续体系上千原子，拿VASP是很难算的，所以才尝试CP2K。
3. 老师，我附上我的NEB.inp和NEB.out，麻烦您看一下，除了基组选择以外，我哪里还需要改进一下？

nuonuo123456

Ying_Zhang 发表于 2026-3-12 20:54
感谢建议
1.第一次跑CP2K过渡态，我就是拿vtst决定插入数量，结果每个点一直停留在Step1，一直持续计算 ...

不用喊老师，大家都是好社友。我个人感觉inp没什么问题，就是CUTOFF是不是有点设的过高了，具体每个元素设多少我不太了解，你可以让其他坛友再看看。但是这个NEB卡住的情况我之前试过插点设多设少都会有卡住的情况，我用的不是最新的版本，不知道最新版是不是会有这样的情况，我感觉是软件自己的问题。你可以换一下版本再试试。

Ying_Zhang

nuonuo123456 发表于 2026-3-12 22:02
不用喊老师，大家都是好社友。我个人感觉inp没什么问题，就是CUTOFF是不是有点设的过高了，具体每个元素 ...

嗯，好的。Mg元素这个值存在这个问题，它的值较高，我也做了cutoff/recutoff测试。我现在用的就是最新的2026.1版本，请问当您出现卡顿如何解决？我会从头计算相同内容，但是后续仍会卡顿，但是有人我看到会继续卡的那步进行续算，但是说不准确？您是如何处理的，哪个方法更准确？

northern_mar

Ying_Zhang 发表于 2026-3-12 20:54
感谢建议
1.第一次跑CP2K过渡态，我就是拿vtst决定插入数量，结果每个点一直停留在Step1，一直持续计算 ...

关于1，会不会可能是有些副本没有同时SCF收敛呢？因为他要所有副本同时收敛完才会到下一步。当然也有可能是软件问题，个别问题个别分析。
关于3，可以取消掉EXTRAPOLATION USE_PREV_P的注释
能量不准的问题：可以使用MN15L搭配TZVPP
Mg的截断值过高可以使用全电子基组，可以对比一下提高cutoff和使用全电子基组哪个更节约耗时

Ying_Zhang

northern_mar 发表于 2026-3-13 08:57
关于1，会不会可能是有些副本没有同时SCF收敛呢？因为他要所有副本同时收敛完才会到下一步。当然也有可能 ...

关于1 我后面也觉得是存在的插入NEB的点结构在某一步较难收敛 所以导致所有副本停留在某一步 一直循环
关于3 感谢建议 我会马上进行尝试