CP2K计算NH3在氧化物表面解离过程难收敛

喝酸奶的巧乐兹 · 发表于 Post on 2025-11-5 10:17:43

本帖最后由喝酸奶的巧乐兹于 2025-11-5 10:19 编辑

各位老师好！我用cp2k计算NH3在Fe2O3 001表面解离的过渡态，目前在优化初末态。
初态的NH3吸附已经优化完成，但末态N*+3H*的吸附怎么都优化不收敛，默认精度下SCF难收敛，降低精度到EPS_SCF 3.0E-05后离子步则不收敛，甚至体块中Fe出现较大位移
对角化、OT均尝试过，吸附位点是根据文献选择的，DFT+U的U值则是与计算初态NH3吸附时保持一致
四种cif文件已给出，输入文件和输出文件已给出，四种结构的输入文件仅原子坐标不同，
如何修改使末态几何优化顺利进行呢？恳请各位老师给出建议！

spinel · 发表于 Post on 2025-11-5 11:10:37

同学，可否上传你已经优化收敛的初态输入输出文件，我做的Fe2O3(001)表面一直无法收敛，可以参考一下你的设置吗？

喝酸奶的巧乐兹 · 发表于 Post on 2025-11-5 15:49:43

spinel 发表于 2025-11-5 11:10
同学，可否上传你已经优化收敛的初态输入输出文件，我做的Fe2O3(001)表面一直无法收敛，可以参考一下你的设 ...

你可以把上传的inp文件里的N和H删去试试，删去就是Fe2O3 001表面了

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

[CP2K] CP2K计算NH3在氧化物表面解离过程难收敛

浏览过的版块