求助对ZnO表面吸附O2结构优化时能量一直震荡

申野会

我使用vasp对ZnO表面吸附O2结构优化时，能量一直震荡怎么办。下面是我的incar，原子共有80多个。
# ZnO(100)H钝化表面吸附O₂ - 弛豫+能量计算SYSTEM = ZnO_100_H_O2_ads
# 基础设置（精度改为中等）
ISTART = 0        # 从头计算
ICHARG = 2        # 电荷密度初始化
PREC = Medium     # 中等精度（平衡速度与精度）
LREAL = .FALSE.   # 全平面波（保证界面精度）
# 交换关联泛函
GGA = PE          # PBE泛函（匹配Zn/O/H PAW_PBE赝势）
# 电子步收敛
ENCUT = 400       # 截断能400 eV（匹配Zn/O赝势）
EDIFF = 1E-4      # 电子步精度（保证吸附能准确）
NELM = 100        # 电子步最大迭代次数
# 离子步弛豫
NSW = 100         # 最大弛豫步数
IBRION = 2        # 共轭梯度法弛豫
ISIF = 2          # 仅弛豫原子位置（固定晶胞）
EDIFFG = -0.02    # 离子步力阈值
# 磁性（O₂顺磁必需）
ISPIN = 2         # 开启自旋极化
# 对称性
ISYM = 0          # 关闭对称性（吸附后破缺）
# 输出控制
ISMEAR = 0        # Gaussian展宽
SIGMA = 0.02      # 展宽系数（平衡精度）
LCHARG = .TRUE.   # 保存电荷密度（分析界面）
LWAVE = .FALSE.   # 不保存波函数（省空间）

阮文升-河内国大

你这是AI生成的输入文件？仔细看：
计算新手不要用AI
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=52382
(出处: 计算化学公社)

EDIFF = 1E-4 精度太差，建议 1E-6 起步，条件允许的话 1E-8
你这里面磁矩最好用MAGMOM手动设，被吸附的O2有磁矩，其它原子全都没有磁矩。否则默认初猜是所有原子都有磁矩，严重偏离正解，自然难以收敛。

申野会

阮文升-河内国大 发表于 2026-3-22 17:51
你这是AI生成的输入文件？仔细看：
计算新手不要用AI
http://bbs.keinsci.com/forum.php?mod=viewthread& ...

好的，谢谢大佬

nuonuo123456

PREC=Medium不太够用，LREAL用auto比较好吧，sigma可以用0.05或者0.1，把色散校正打开，EDIFF用-5也可以，算是够用，条件允许就-6，-7就可以，vasp官网上well converged calculations也就建议到-6的精度，不建议到-8这么高的精度。