计算化学公社

标题: 求助对ZnO表面吸附O2结构优化时能量一直震荡 [打印本页]

作者
Author: 申野会 时间: 2026-3-22 16:51
标题: 求助对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. # 不保存波函数（省空间）

作者
Author: 阮文升-河内国大 时间: 2026-3-22 17:51
本帖最后由阮文升-河内国大于 2026-3-22 20:10 编辑

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

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

作者
Author: 申野会 时间: 2026-3-28 13:51

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

好的，谢谢大佬

作者
Author: nuonuo123456 时间: 2026-3-28 15:18
本帖最后由 nuonuo123456 于 2026-3-28 15:22 编辑

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

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