计算化学公社

标题: CP2K计算NH3在氧化物表面解离过程难收敛 [打印本页]

作者
Author: 喝酸奶的巧乐兹 时间: 2025-11-5 10:17
标题: CP2K计算NH3在氧化物表面解离过程难收敛
本帖最后由喝酸奶的巧乐兹于 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文件已给出，输入文件和输出文件已给出，四种结构的输入文件仅原子坐标不同，
如何修改使末态几何优化顺利进行呢？恳请各位老师给出建议！

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

作者
Author: 喝酸奶的巧乐兹 时间: 2025-11-5 15:49

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

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

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