计算化学公社

标题: 求助单原子铁力不收敛的原因 [打印本页]
作者
Author: WVzzz    时间: 2023-2-28 14:48
标题: 求助单原子铁力不收敛的原因
本帖最后由 WVzzz 于 2023-2-28 16:22 编辑

求助各位老师，Fe单原子结构优化过程出现电子步收敛，离子步不收敛的情况，目前已经尝试更换IBRION=3却不行，直接散掉，目前的想法是1.真空层原因：目前为14A-15A
2.POTIM原因，默认0.5打算改小一点
希望各位老师给予建议，谢谢
作者
Author: 含光君    时间: 2023-2-28 16:01
1. IBRION可改2，若初始结构差应选择CG法；
2.  EDIFF增加收敛精度，可能因为电子步精度低而离子步收敛要求高而导致后者不收敛；
3. 建模上如非必要不需要叠四层graphene；
4. 提问时POTCAR没必要传，但最好给OUTCAR和OSZICAR。
作者
Author: WVzzz    时间: 2023-2-28 16:22
本帖最后由 WVzzz 于 2023-2-28 16:38 编辑
含光君 发表于 2023-2-28 16:01
1. IBRION可改2，若初始结构差应选择CG法；
2.  EDIFF增加收敛精度，可能因为电子步精度低而离子步收敛要 ...

好的，多谢，再请问一下，在表面结构优化过程需要加上偶极矫正的LDIPOL、IDIPOL等参数嘛
作者
Author: 含光君    时间: 2023-3-1 00:57
对于你目前建立的体系，z方向不对称应该加上偶极矫正。
作者
Author: WVzzz    时间: 2023-3-14 16:02
含光君 发表于 2023-3-1 00:57
对于你目前建立的体系，z方向不对称应该加上偶极矫正。

多谢老师



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3