计算化学公社

标题: Fe原胞优化不收敛，能量一直震荡，望大佬相助~ [打印本页]

作者
Author: TaylorS 时间: 2023-4-20 21:03
标题: Fe原胞优化不收敛，能量一直震荡，望大佬相助~
尝试了BFGS和CG，能量一直处于震荡状态，不知道是什么原因...以下是我的输入文件
(, 下载次数 Times of downloads: 25)

作者
Author: WilliamH 时间: 2023-4-20 21:38
cp2k处理金属体系比较难。论坛里之前有关cp2k加速收敛的讨论，可以试试 http://bbs.keinsci.com/forum.php ... 630&fromuid=246

作者
Author: sobereva 时间: 2023-4-20 23:25
你这体系哪是Fe的原胞？

要考虑反铁磁性耦合状态的可能，恰当设置初猜磁矩

说清楚怎么个震荡法，又没输出文件和轨迹文件，别人根本无从判断。如果确实是几何优化震荡而没有SCF不收敛的事，把BFGS步长上限设小再试。
如果是SCF不收敛，尝试各种帮助解决SCF不收敛的做法，诸如减小ALPHA、加大NBROYDEN、加大CUTOFF、加大MAX_SCF等等。对于几何优化，靠前阶段SCF不收敛可以暂时性无视，后期能SCF收敛就行。另外

第一个晶胞矢方向k点用5太浪费，根本没必要这么大。

作者
Author: TaylorS 时间: 2023-4-21 19:34

sobereva 发表于 2023-4-20 23:25
你这体系哪是Fe的原胞？

要考虑反铁磁性耦合状态的可能，恰当设置初猜磁矩

感谢卢老师！在尝试了几个初猜磁矩后最终将磁矩设为6，跑了55步后结构收敛了！

作者
Author: TaylorS 时间: 2023-4-21 19:35

WilliamH 发表于 2023-4-20 21:38
cp2k处理金属体系比较难。论坛里之前有关cp2k加速收敛的讨论，可以试试 http://bbs.keinsci.com/forum.php? ...

嗯嗯，谢谢~

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