求助，CP2K对无定型碳模型结构优化不收敛

haqipaqi · 发表于 Post on 6 hour ago

本帖最后由 haqipaqi 于 2026-4-26 10:27 编辑

本人是cp2k的计算新手。使用cp2k进行文献中无定型碳模型的结构优化。使用CG算法时，除了Maximum gradient其他的都收敛了，后换成了 BFGS方法续算，但是发现跑了500步，更差了。步长0.02.怎么样才能使得结构收敛呢。原模型是文献中使用LAMMPS的Tersoff力场熔融淬火生成的。做结构优化是想要后面跑AIMD来生成机器学习势的数据。

Uus/pMeC6H4-/キ · 发表于 Post on 6 hour ago

“原模型是文献中使用LAMMPS的Tersoff力场熔融淬火生成的”具体是哪篇文献（请给出标题和doi号），从什么晶体结构出发做的熔融？这一帧结构是高温熔融的时候取的，还是足够缓慢地退火降温到0 K再取的？优化完的结果要用来做什么，预计优化精度对目标量的影响有多大？

粗略来说，分子力场只考虑原子层面的问题，而量子化学计算得考虑原子核和电子的问题，二者对势能面的描述不完全相符，不能指望分子力场下能量极小点的结构在量子化学看来也一定有化学意义/能量较小，所以分子力场模拟产生的扭曲复杂结构拿来做量子化学计算不收敛（无论是SCF还是几何优化方面）并不奇怪。

sobereva · 发表于 Post on 1 hour ago

本来CG就更为robust，能用CG顺利算出想要的结果就没必要改成BFGS
注意监控SCF收敛情况，确保优化跑起来后SCF都充分收敛了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助，CP2K对无定型碳模型结构优化不收敛