计算化学公社

标题: 求助：表面吸附分子后优化问题 [打印本页]

作者
Author: 涨涨涨 时间: 2021-9-10 11:03
标题: 求助：表面吸附分子后优化问题
vasp初学者：在面心立方Au（111）表面吸附N2分子然后做优化，无论是top，bridge，还是hollow位，算着算着N2分子就会远离表面跑到真空层中，没有达到设置的NSW就计算结束，OUTCAR中也没有reached required accuracy - stopping structural energy minimisation。
OSZICAR最后为
  N    E                   dE          d eps    ncg    rms       rms(c)DAV: 1 -0.126410230295E+03 -0.86555E-07 -0.22073E-06  3020 0.524E-03 0.853E-04
RMM: 2 -0.126410230314E+03 -0.18190E-07 0.00000E+00 89 0.101E-03
148 F= -.12641023E+03 E0= -.12640915E+03  d E =-.136229E-03
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.126410230344E+03 -0.48385E-07 -0.47422E-07  2800 0.280E-03 0.159E-03
RMM: 2 -0.126410230576E+03 -0.23236E-06 0.00000E+00 185 0.590E-04
149 F= -.12641023E+03 E0= -.12640915E+03  d E =-.136491E-03
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.126410230366E+03 -0.22119E-07 -0.16157E-06  3020 0.586E-03 0.134E-03
RMM: 2 -0.126410230797E+03 -0.43080E-06 0.00000E+00 248 0.696E-04
150 F= -.12641023E+03 E0= -.12640915E+03  d E =-.136712E-03
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.126410230342E+03 0.23705E-07 -0.20677E-06  3060 0.661E-03 0.692E-04
RMM: 2 -0.126410230466E+03 -0.12408E-06 0.00000E+00 120 0.511E-04
151 F= -.12641023E+03 E0= -.12640915E+03  d E =-.136382E-03
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.126410230336E+03 0.64247E-08 -0.58386E-07  2780 0.356E-03 0.422E-04
RMM: 2 -0.126410230364E+03 -0.28755E-07 0.00000E+00 78 0.428E-04
152 F= -.12641023E+03 E0= -.12640915E+03  d E =-.136280E-03
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.126409864245E+03 0.36609E-03 -0.24175E-01  5240 0.756E-01 0.550E-02
RMM: 2 -0.126410107672E+03 -0.24343E-03 -0.26896E-03  5175 0.105E-01 0.387E-02
RMM: 3 -0.126410105548E+03 0.21240E-05 -0.78674E-05  2888 0.183E-02
153 F= -.12641011E+03 E0= -.12640907E+03  d E =-.114629E-04

我的INCAR为
SYSTEM = Aucalculations
ISMEAR = 1
SIGMA = 0.1 #改为0.2还是上面的情况
ENCUT = 450
EDIFF = 1E-4 #减小还是上面情况
EDIFFG = -0.02
LDIPOL = .TRUE.
IDIPOL = 3
IBRION = 2
POTIM = 0.5#减小还是上面情况
NSW = 500
ISIF = 2
LREAL = Auto
ALGO = F

之后查找收敛解决办法又将ibrion改为1后，很快收敛（约三十多个离子步）OUTCAR中出现reached required accuracy - stopping structural energy minimisation，但查看CONTCAR发现结构与大牛计算的结果相差较大，并且与物理事实不符。

我的问题是：在使用ibrion=2时对于这种吸附分子在计算过程中远离表面并且最后不收敛的情况该如何计算相关能量？出现上面情况是否说明吸附分子不能在表面产生化学吸附？如果要计算吸附能OUTCAR没有收敛如何计算？

作者
Author: wsz 时间: 2021-9-11 09:11
加上D3校正试试（IVDW=11or12）

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