求助：CsPbI3的slab模型计算离子迁移能末态能量很高

josy

我想算下CsPbI3表面的I迁移到外界的迁移能垒，参考了别人算的进行了建模，初态是就是一个CsPbI3的slab，末态是将表面中心的I挪到外界。在优化末态的结构时，发现不论怎么调外界这个I的位置，都要比初态的能量高0.6eV左右。然后我就尝试了一下把优化好后的末态和初态用neb插image，neb结束之后发现迁移能曲线是递增的，各位老师请问这是为什么呢，是因为末态的位置找的不对吗？我也算过表面的I迁移到他表面临近I空位的迁移能，这个是正常的，但实验的人要求要算从表面迁移到外界的迁移能。（第一张图是我的计算；第二张图是别人算的，他的末态的位置我也试过，我优化出来的末态结构是跑出来的I会把临近的Pb扯出表面，也没有算出只比初态高0.2eV左右；结构优化的INCAR和NEB计算的INCAR贴在后面的）




结构优化的INCAR：
ISYM = 0
IVDW = 11
# I/O设置(读入、读出)
ISTART = 0  # 0代表一个全新的计算
ICHARG = 2
LWAVE = .FALSE.  
LCHARG = .F.
LVOT = .FALSE.
LELF = .FALSE.  #以上四个文件一般在需要计算准确电子结构（自洽运算）时才会出现

# Electronic Relaxation(电子步)
ENCUT = 450    #平面阶段能，由赝势决定（POTCAR中的ENMAX*1.25~1.50）（数值越大，精度越高）
NELM = 100    #最大电子步数量，默认60步（难收敛体系，如过渡金属，可设置为200、300~）
ALGO = Normal    #自洽循环算法（Fast,Normal,VaryFast)
PREC = Normal    #精度设置
ISMEAR = 0    # 金属0+，非金属0-（金属不可以设置为负数）
SIGMA = 0.05  # 划分布里渊区,展宽
EDIFF = 1E-4    #电子步收敛精度
#AMIX = 0.1    #
#BMIX = 0.01   #加快收敛

# Ionic Relaxation (离子步)
IBRION = 2     #离子弛豫的算法（2较稳定）
NSW = 300      #最大离子步的个数
EDIFFG = -0.05   #离子步收敛标准（默认是EDIFF*10）
ISIF = 2    # 对晶胞的弛豫方法(3:全弛豫 2:固定体积的弛豫 4:固定体积但允许形状改变)

# Polarization (含磁性体系)
#ISPIN = 2       # 自旋极化   
#LOSRBIT = .False.  #考虑非线性磁性（默认不考虑）
#MAGMOM = 18*1 6*3 6*1

# Parallization(并行计算的相关设定)
LREAL = Auto   #计算进行在实空间或是导空间
NPAR = 8    #取节点数量或核心数量的根号值
#LPLANE = .FALSE. #针对平面波的设置


NEB计算的INCAR：
#for neb
  IBRION =   3
  POTIM  =   0
  IOPT = 1
  ICHAIN = 0
  LCLIMB = .TRUE.
  SPRING =  -5
  IMAGES = 5

ISYM = 0
IVDW = 11
# I/O设置(读入、读出)
ISTART = 0  # 0代表一个全新的计算
ICHARG = 2
LWAVE = .FALSE.  
LCHARG = .F.
LVOT = .FALSE.
LELF = .FALSE.  #以上四个文件一般在需要计算准确电子结构（自洽运算）时才会出现

# Electronic Relaxation(电子步)
ENCUT = 450    #平面阶段能，由赝势决定（POTCAR中的ENMAX*1.25~1.50）（数值越大，精度越高）
NELM = 100    #最大电子步数量，默认60步（难收敛体系，如过渡金属，可设置为200、300~）
ALGO = Normal    #自洽循环算法（Fast,Normal,VaryFast)
PREC = Normal    #精度设置
ISMEAR = 0    # 金属0+，非金属0-（金属不可以设置为负数）
SIGMA = 0.05  # 划分布里渊区,展宽
EDIFF = 1E-4    #电子步收敛精度
#AMIX = 0.1    #
#BMIX = 0.01   #加快收敛

# Ionic Relaxation (离子步)
#IBRION = 2     #离子弛豫的算法（2较稳定）
NSW = 300      #最大离子步的个数
EDIFFG = -0.05   #离子步收敛标准（默认是EDIFF*10）
ISIF = 2    # 对晶胞的弛豫方法(3:全弛豫 2:固定体积的弛豫 4:固定体积但允许形状改变)

# Polarization (含磁性体系)
#ISPIN = 2       # 自旋极化   
#LOSRBIT = .False.  #考虑非线性磁性（默认不考虑）
#MAGMOM = 18*1 6*3 6*1

# Parallization(并行计算的相关设定)
LREAL = Auto   #计算进行在实空间或是导空间
NPAR = 10    #取节点数量或核心数量的根号值
#LPLANE = .FALSE. #针对平面波的设置