计算化学公社

标题: 关于NiFe掺杂AlOOH模型做吸附计算结构优化力不收敛问题求助 [打印本页]

作者
Author: FalseAndRight 时间: 2025-2-16 11:58
标题: 关于NiFe掺杂AlOOH模型做吸附计算结构优化力不收敛问题求助
体系是：NiFe掺杂AlOOH，晶面切的是100（参考文献），由于实验测得Ni、Fe是+2、+3价态都存在，觉得都直接替换无法体现，因此做了一个氧空位，可以看poscar文件，然后结构已经优化成功了，添加了OH，做小分子吸附。
输入文件内容如下：（先做粗收敛，再做精收敛）
# ------------- DFT-D3 矫正模块 -------------
IVDW = 11 (DFT-D3 method of method with no damping)

# Global Parameters
ENCUT = 400 (Cut-off energy for plane wave basis set, in eV)
ALGO = N
# --------- 自旋极化模块 start -------------
ISPIN = 2 # 2表示考虑自旋极化
AMIX = 0.05 # 这个参数需要进行梯度矫正来判断是否合适
BMIX=0.0001
AMIX MAG=0.8
BMIX MAG=0.0001
# --------- 是否输出电磁轨道相关的信息 -------------
LORBIT = 11
# --------- 材料本身决定的参数 -------------
ISMEAR= 0 (Gaussian smearing, metals:1)
SIGMA=0.05 (Smearing value in eV, metals:0.2)

# Electronic Relaxation
NELM =  200          (Max electronic SCF steps)
EDIFF  =  1E-03       (SCF energy convergence, in eV)
# Ionic Relaxation
NSW =  200       (Max ionic steps)
IBRION =  2          (Algorithm: 0-MD, 1-Quasi-New, 2-CG)
ISIF =  2          (Stress/relaxation: 2-Ions, 3-Shape/Ions/V, 4-Shape/Ions)
EDIFFG = -5E-01       (Ionic convergence, eV/AA)

# ---------- 偶极矫正 -------------
LDIPOL = .TRUE.
IDIPOL = 3    # 针对z轴

为了计算得快一点，K点采用的Gamma点，内容如下：
K-Spacing Value to Generate K-Mesh: 0.000
0
Gamma
1 1 1
0.0  0.0  0.0

目前的情况是：无论怎么修改输入文件，包括更改IBRION=1，ISMEAR改为正数等等都没有效果，我在思考是不是因为结构本身就不太好，导致无法收敛，但是对于这方面没有经验，不知道如何判断是否结构存在问题，只能不停尝试修改。有没有高手指点我一下，万分感谢。
下面是力的收敛情况：基本上每次修改INCAR文件再次运行都是这样的结果，收敛到比较大的值，比如2.0/1.0/0.5左右就反复震荡不继续向下收敛了。但是只优化结构，是可以正常收敛的，一加上小分子做吸附计算，就无法收敛。
(, 下载次数 Times of downloads: 1)
三维结构图如下：
(, 下载次数 Times of downloads: 3)
(, 下载次数 Times of downloads: 3)
(, 下载次数 Times of downloads: 3)

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