24小时热门版块排行榜    

CyRhmU.jpeg
查看: 514  |  回复: 0
【悬赏金币】回答本帖问题,作者护发专家将赠送您 10 个金币

护发专家

铜虫 (小有名气)

[求助] 求助!算吸附,离子步能量和电子步能量差距好大

N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.887594668591E+03   -0.13338E-03   -0.20532E-01  3648   0.614E-01    0.102E-01
RMM:   2    -0.887595006596E+03   -0.33801E-03   -0.41688E-03  3648   0.981E-02    0.868E-02
RMM:   3    -0.887595007513E+03   -0.91721E-06   -0.34736E-04  3645   0.268E-02    0.805E-02
RMM:   4    -0.887594989210E+03    0.18303E-04   -0.97734E-05  3490   0.126E-02    0.607E-02
RMM:   5    -0.887594988735E+03    0.47464E-06   -0.75309E-05  3350   0.100E-02    0.533E-02
RMM:   6    -0.887594981209E+03    0.75265E-05   -0.46260E-05  2924   0.804E-03
388 F= -.90430791E+03 E0= -.90434664E+03  d E =-.632607E-03  mag=     0.4967
我用MD的结构算表面吸附,结果用离子步收敛能量-904eV算吸附能竟然-20ev,明显不对,仔细找了之后发现,用电子步收敛能量-887ev算就刚好,请问为何电子步能量和离子步能量差距会有这么大。是不是INCAR参数设置有问题。

LCHARG  =.FALSE.
   LWAVE   =.FALSE.
   PREC    = Normal
   LREAL   = Auto
   ALGO    = Fast
   EDIFF   = 1e-5
   EDIFFG  =-0.03
   NSW     = 500
   NELM    = 100
   NELMIN  = 6
   ENCUT   = 400
   IVDW    = 11
   ISPIN   = 2
   IDIPOL  = 3
   NGX     =122
   NGY     =70
   NGZ     =146
Ionic relaxation
   ISIF    = 2
   IBRION  = 2
DOS related values
   ISMEAR  = 1
   SIGMA   = 0.2
Parallelization
   NSIM    = 4
   NPAR    = 4
   LPLANE  = .TRUE.
   LSCALU  = .FALSE.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 护发专家 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见