计算化学公社

标题: CP2K计算周期性氧化铈和三氧化二铈体系结果不收敛 [打印本页]

作者
Author: zhifanfan 时间: 2025-3-24 11:13
标题: CP2K计算周期性氧化铈和三氧化二铈体系结果不收敛
各位老师好，我想使用CP2K计算周期性氧化铈和三氧化二铈体系的波函数，使用Multiwfn软件生成输入文件，并更改所用势函数，结果CeO2体系收敛成功计算出波函数，而Ce2O3却无法收敛，不知道改=该怎么解决，氢各位老师点拨。

作者
Author: kilig 时间: 2025-3-25 08:57
我看到你的inp文件里面只给Ce设置了U值而没有给初始化磁矩，可能是因为这个原因吧，我之前跑含有Co和Fe体系的MTD，合适的磁矩是非常有利与体系收敛。
我之前常用的方法都是vasp跑一遍得到原子的磁化率和磁化强度为参考，从而进行设置初始化磁矩。

作者
Author: zhifanfan 时间: 2025-3-25 16:02

kilig 发表于 2025-3-25 08:57
我看到你的inp文件里面只给Ce设置了U值而没有给初始化磁矩，可能是因为这个原因吧，我之前跑含有Co和Fe体系 ...

这个我也不是很清楚，但是二氧化铈就能很轻松收敛

作者
Author: xmuqyj 时间: 11 hour ago

zhifanfan 发表于 2025-3-25 16:02
这个我也不是很清楚，但是二氧化铈就能很轻松收敛

楼主，我开OT算CeO2死活不收敛，开对角化反而能更快收敛吗

作者
Author: sobereva 时间: 5 hour ago

xmuqyj 发表于 2026-6-2 12:01
楼主，我开OT算CeO2死活不收敛，开对角化反而能更快收敛吗

此体系没必要刻意开OT
对角化能收敛用对角化就完了

北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里讲这个例子就用的Multiwfn创建CP2K输入文件默认的对角化方式算的，SCF收敛非常顺利
(, 下载次数 Times of downloads: 2)

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