计算化学公社

标题: Ni水界面跑aimd的scf不收敛 [打印本页]

作者
Author: Wednesday 时间: 2025-11-6 17:15
标题: Ni水界面跑aimd的scf不收敛
本帖最后由 Wednesday 于 2025-11-6 20:38 编辑

各位老师，我目前在用cp2k跑Ni水界面的aimd，结构是用MS和packmol生成的，目前遇到一些问题想请教一下：
1.首先没有对整个体系进行结构优化就直接进行aimd了，因为之前做类似的根据博文修改后cp2k仍然不能收敛，实在是不想用vasp去做400多原子的aimd，但是cp2k做Ni表面又不好收敛，所以跳过了这个步骤直接加上水层开始跑aimd，out文件也可以看到前几轮跑的不太合理。我参考了别人做的Ni/NiOH2水界面体系aimd，发现out文件中大部分时候scf都不能收敛，只有少部分可以达到收敛标准，但跑到一定步数后scf收敛前convergence最后都在0.00001000左右也就是收敛限上下，想请教一下，这样跑完全程的话结果合理么？
2.一开始是用abacus的补偿电荷方法去模拟恒电位，想请教下cp2k中有没有类似的方法？之前看论坛中有用匀强电场去模拟电极电位的，这么做合理么？
感谢各位老师！

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