vasp dimer跑一个甲基环己烷在Pt111表面上的脱氢的过渡态散掉不收敛

赵非命

各位老师，我用vasp在跑一个甲基环己烷在Pt111表面上的脱氢第一步过渡态搜索，就是C7H14-Pt111变成C7H13+H-Pt111，两个优化完后在过渡态搜索中先用nebmake.pl IS/POSCAR FS/POSCAR 7然后采用了其中文件夹04的poscar，再用modemake.pl IS/POSCAR FS/POSCAR生成MODECAR开始跑SYSTEM = default


Start parameter for this Run:
  ISTART = 0
  ICHARG = 2
  INIWAV = 1
# NCORE=26
  NPAR= 1




Electronic Relaxation:
  ENCUT = 400
  ENAUG = 650
# PREC = Accurate
  IALGO = 48
# ALGO=Normal
  NELM = 100
  NELMIN = 8 # surface = 8, bulk = 4
  NELMDL = -12 # When IALGO = 38, 0
  EDIFF = 1E-04
# NBANDS = 1600
  GGA = PE
# VOSKOWN = 1
  LREAL = Auto
  WEIMIN = 0


Ionic Relaxation:
  EDIFFG = -0.03
  NSW = 400
  IBRION = 3
# IBRION=1：对于初始结构接近局域极小时是最好的选择
# IBRION=2：比较可靠，弛豫困难时，推荐使用
# IBRION=3：对非常糟糕的初始结构进行弛豫时可考虑使用
  ISIF = 2
# 在每个压力点下优化晶格常数，保持原子位置不动     ISIF=7
# 在每个压力点下优化原子位置，保持晶格常数不动     ISIF=2
  POTIM = 0
# POTIM代表步长，步长太大，收敛的不是极值；POTIM 太小是否会影响优化速度
IOPT = 2 #（VTST的代码）
  ISYM = 0
# ISYM：决定是否使用对称性。取值：-1~3（-1,0-不使用，1,2,3-使用）
# 使用US-PP 1；使用PAW 2；ISYM=3时仅考虑力和应力张量的对称性，而电荷密度是非对称的。


DOS related values:
  SIGMA = 0.05
# 越小精度越高，vasp默认为0.2
  ISMEAR = 0
# 计算能带：导体的话，用ISMEAR=1；半导体或绝缘体，用ISMEAR=0 。
# 计算 DOS： ISMEAR一般取-5都可以算；金属可用0、1，非金属不要大过0，体材料可用－4、－5（面的话就用－1、0吧？）




Spin polarized:
  ISPIN = 1
# 1非自旋，得到up的dos；2自旋，比如Fe Co等的氧化物，得到双向的dos。
#  MAGMOM = 1(amount)*2(initial moment) 3*2


File writing
  LWAVE = .FALSE.
  LCHARG = .FALSE.
  
Calculation of DOS
#  RWIGS = 1.286 1.000 0.370 # It is not necessary to use it for PAW potential.
#  NPAR = 1
#  LORBIT = 1


# vdW Calculation
# PARAM1 = 0.1234
# PARAM2 = 1.0000
# LUSE_VDW = .TRUE.
# AGGAC = 0.0000
#  IVDW = 1


# U Calculation
dimer Calculation
  ICHAIN = 2
  DdR = 0.005
  DRotMax = 10
  DFNMin = 0.01
  DFNMax = 1.0
但是跑完原子Pt放开的表面两层原子乱了 dimer也是越跑越大，这个情况应该只能解决，我还能继续复制contcar继续跑吗，是不是已经散掉了

后到咸阳斯蒂庞

大概率是MODECAR指定的初始方向不太合适，正常来讲dimer的斜率绝大多数应该是负值。不如先用CI-NEB的结果做个振动分析，然后把沿反应方向振动的虚频输入到MODECAR里。具体操作步骤详见白皛白老师的帖子http://bbs.keinsci.com/thread-41457-1-1.html，这个方法基本上百试不爽且效率很高。