24小时热门版块排行榜    

查看: 1582  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gleerat

木虫 (正式写手)

[交流] 【求助成功】EDIFFG设为-1E-4,但是TOTFORCE到0.001eV/A就退出了 已有8人参与

官方指南上说道:
引用回帖:
EDIFFG defines the break condition for the ionic relaxation loop. If the change in the total (free) energy is smaller than EDIFFG between two ionic steps relaxation will be stopped. If EDIFFG is negative it has a different meaning: In this case the relaxation will stop if all forces are smaller than |EDIFFG|. This is usually a more convenient setting.

我设置的EDIFFG=-1E-4,可是TOTAL-FORCE在meV/A量级时就退出了:
引用回帖:
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
     -1.39378      0.00000      2.82373         0.000491      0.000000     -0.000582
     -4.20826      0.00000      8.45886         0.000000      0.000000     -0.000424
      1.40275      0.00000      8.48186         0.000000      0.000000      0.001014
     -1.41173      0.00000     14.11700        -0.000491      0.000000     -0.000582
     -1.37499     -5.50242      2.81574         0.001671      0.000000      0.000416
     -4.20826     -5.50242      8.47020         0.000000      0.000000     -0.002063
      1.40275     -5.50242      8.49494         0.000000      0.000000     -0.000786
     -1.43052     -5.50242     14.10902        -0.001671      0.000000      0.000416
      1.40275     -2.74676      2.82412         0.000000     -0.000593     -0.000412
     -1.40165     -2.76141      8.47504         0.000417      0.000661     -0.000208
      4.20715     -2.76141      8.47504        -0.000417      0.000661     -0.000208
      1.40275     -2.72309     14.10978         0.000000     -0.000324      0.000215
      1.40275     -8.25807      2.82412         0.000000      0.000593     -0.000412
     -1.40165     -8.24343      8.47504         0.000417     -0.000661     -0.000208
      4.20715     -8.24343      8.47504        -0.000417     -0.000661     -0.000208
      1.40275     -8.28175     14.10978         0.000000      0.000324      0.000215
      0.00288     -1.37439     -0.01437         0.003051     -0.000512      0.001832
     -2.80228     -1.37924      5.63469        -0.000181     -0.001731     -0.001290
      2.80385     -1.37653      5.66093        -0.000879     -0.000745     -0.002184
     -0.00406     -1.37259     11.30601         0.001224     -0.000351      0.000525
      0.00472     -6.88185      0.00131        -0.004483      0.001564     -0.001495
     -2.81690     -6.86016      5.58381         0.001325      0.000220     -0.001528
      2.80560     -6.87768      5.65247         0.001363     -0.000520      0.000544
     -0.00194     -6.87960     11.34442        -0.004149     -0.000214     -0.001433
      0.00472     -4.12298      0.00131        -0.004483     -0.001564     -0.001495
     -2.81690     -4.14468      5.58381         0.001325     -0.000220     -0.001528
      2.80560     -4.12716      5.65247         0.001363      0.000520      0.000544
     -0.00194     -4.12523     11.34442        -0.004149      0.000214     -0.001433
      0.00288     -9.63045     -0.01437         0.003051      0.000512      0.001832
     -2.80228     -9.62560      5.63469        -0.000181      0.001731     -0.001290
      2.80385     -9.62831      5.66093        -0.000879      0.000745     -0.002184
     -0.00406     -9.63225     11.30601         0.001224      0.000351      0.000525
     -0.00010     -4.12716      5.65247        -0.001363      0.000520      0.000544
     -2.81023     -4.12298     11.29458         0.004483     -0.001564     -0.001495
      2.80744     -4.12523     11.34442         0.004149      0.000214     -0.001433
      0.01139     -4.14468     16.87708        -0.001325     -0.000220     -0.001528
      0.00166     -9.62831      5.66093         0.000879      0.000745     -0.002184
     -2.80839     -9.63045     11.27890        -0.003051      0.000512      0.001832
      2.80956     -9.63225     11.30601        -0.001224      0.000351      0.000525
     -0.00323     -9.62560     16.92796         0.000181      0.001731     -0.001290
      0.00166     -1.37653      5.66093         0.000879     -0.000745     -0.002184
     -2.80839     -1.37439     11.27890        -0.003051     -0.000512      0.001832
      2.80956     -1.37259     11.30601        -0.001224     -0.000351      0.000525
     -0.00323     -1.37924     16.92796         0.000181     -0.001731     -0.001290
     -0.00010     -6.87768      5.65247        -0.001363     -0.000520      0.000544
     -2.81023     -6.88185     11.29458         0.004483      0.001564     -0.001495
      2.80744     -6.87960     11.34442         0.004149     -0.000214     -0.001433
      0.01139     -6.86016     16.87708        -0.001325      0.000220     -0.001528
      1.40275     -5.50242     12.67755         0.000000      0.000000      0.005588
     -1.41293      0.00000      6.95033        -0.000273      0.000000      0.002822
     -4.20826      0.00000     12.58582         0.000000      0.000000     -0.004215
      1.40275      0.00000     12.59997         0.000000      0.000000     -0.000004
     -1.39258      0.00000     18.24360         0.000273      0.000000      0.002822
     -1.42630     -5.50242      6.93186         0.001306      0.000000     -0.000019
     -4.20826     -5.50242     12.59332         0.000000      0.000000     -0.000022
     -1.37920     -5.50242     18.22513        -0.001306      0.000000     -0.000019
      1.40275      0.00000      4.35735         0.000000      0.000000      0.004561
     -1.39442      0.00000      9.98792         0.003392      0.000000     -0.002411
      4.19993      0.00000      9.98792        -0.003392      0.000000     -0.002411
      1.40275      0.00000     15.63891         0.000000      0.000000      0.000105
      1.40275     -5.50242      4.35102         0.000000      0.000000      0.001704
     -1.38562     -5.50242     10.02199        -0.000590      0.000000      0.002788
      4.19113     -5.50242     10.02199         0.000590      0.000000      0.002788
      1.40275     -5.50242     15.50333         0.000000      0.000000      0.002125
      1.40275     -2.75728      6.96646         0.000000      0.003008      0.002751
     -1.41894     -2.74070     12.60271        -0.000265     -0.001679     -0.002589
      4.22444     -2.74070     12.60271         0.000265     -0.001679     -0.002589
      1.40275     -2.77262     18.20247         0.000000      0.003103      0.001452
      1.40275     -8.24756      6.96646         0.000000     -0.003008      0.002751
     -1.41894     -8.26414     12.60271        -0.000265      0.001679     -0.002589
      4.22444     -8.26414     12.60271         0.000265      0.001679     -0.002589
      1.40275     -8.23222     18.20247         0.000000     -0.003103      0.001452
     -1.38752     -2.74600      4.33772        -0.001263     -0.000785      0.003097
     -4.20826     -2.75662      9.97505         0.000000      0.002716     -0.001933
      1.40275     -2.76734     10.03288         0.000000      0.003262      0.000582
     -1.41799     -2.74600     15.63100         0.001263     -0.000785      0.003097
     -1.38752     -8.25884      4.33772        -0.001263      0.000785      0.003097
     -4.20826     -8.24822      9.97505         0.000000     -0.002716     -0.001933
      1.40275     -8.23750     10.03288         0.000000     -0.003262      0.000582
     -1.41799     -8.25884     15.63100         0.001263      0.000785      0.003097
-----------------------------------------------------------------------------------
    total drift:                                0.000000      0.000005     -0.000805

[ Last edited by gleerat on 2011-3-15 at 23:20 ]
回复此楼
明察、慎思、笃行
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)

引用回帖:
Originally posted by 小木虫于 at 2011-03-02 15:54:52:
你看看你的OUTCAR中是否有FORCES: max atom, RMS     0.268526    0.061907
这样一行,ediffg设为负值,应该是在最大原子受力小于你的设置时计算就结束了,如上我的最大原子受力为0.268526.    还有就是我用不同 ...

好的,谢谢。
明察、慎思、笃行
11楼2011-03-03 16:53:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

ustbwgs

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1): 谢谢回复 2011-03-02 08:32:06
不太可能吧,我从来没碰到过,应该是每一个数都达到你设置的那个精度才收敛啊
doingwhatishoulddo!
2楼2011-03-01 23:11:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1): 谢谢交流 2011-03-02 08:32:19
-1E-4这个精度非常非常高了,VASP很难达到哦。一般小数点后面两位,比如0.02eV/A,某些计算可以选择小数点后三位,你都4位了,为什么设这么高的精度。
3楼2011-03-01 23:22:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)


zzy870720z(金币+1): 谢谢参与交流 2011-03-02 08:32:37
这个标准好高啊,是不是达到最大离子步数了,默认是NSW=200
估计你这个是跑了200步,你这个标准估计200步很难收敛
心之所向无惧无悔
4楼2011-03-02 00:29:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见