应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】EDIFFG设为-1E-4，但是TOTFORCE到0.001eV/A就退出了
151/2返回列表12下一页
查看: 2105  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
明察、慎思、笃行
1楼 2011-03-01 16:19:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbwgs

木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1): 谢谢回复 2011-03-02 08:32:06
不太可能吧，我从来没碰到过，应该是每一个数都达到你设置的那个精度才收敛啊
赞一下(2人) 回复此楼
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位了，为什么设这么高的精度。
赞一下(2人) 回复此楼
3楼2011-03-01 23:22:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★
zzy870720z(金币+1): 谢谢参与交流 2011-03-02 08:32:37
这个标准好高啊，是不是达到最大离子步数了，默认是NSW=200
估计你这个是跑了200步，你这个标准估计200步很难收敛
赞一下(1人) 回复此楼
心之所向无惧无悔
4楼2011-03-02 00:29:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbwgs

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
版主和专家说的对啊，如果原子个数比较多的话，根本就不可能达到这么高的精度，你能量收敛设置个-4也足够了。
赞一下(1人) 回复此楼
doingwhatishoulddo!
5楼2011-03-02 09:56:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by stractor at 2011-03-01 23:22:35:
-1E-4这个精度非常非常高了，VASP很难达到哦。一般小数点后面两位，比如0.02eV/A,某些计算可以选择小数点后三位，你都4位了，为什么设这么高的精度。

VASP默认的原子力收敛标准为-1E-3，我-1E-4的收敛标准不算太离谱吧？
但是同样的标准，有时候原子力的收敛精度就会很低，这里就只能达到0.01量级。
引用回帖:
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
     -0.00000      0.00000      0.00000        -0.000000     -0.000000      0.000000
      1.98805     -0.00001      2.59903        -0.000005      0.000005      0.000005
     -0.00000      1.98805      2.59903         0.000005     -0.000005     -0.000006
     -0.00001     -0.00001      3.76228         0.000009      0.000011      0.083150
     -0.00001     -0.00002      6.63385        -0.000008     -0.000010     -0.083150
-----------------------------------------------------------------------------------
    total drift:                                0.000000     -0.000000      0.000000

我看过我以前的结果，发现基本上以达到meV/A量级居多。
赞一下 回复此楼
明察、慎思、笃行
6楼2011-03-02 10:23:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by sunyang1988 at 2011-03-02 00:29:19:
这个标准好高啊，是不是达到最大离子步数了，默认是NSW=200
估计你这个是跑了200步，你这个标准估计200步很难收敛

我曾经使用的是NSW=80，但是看到原子力达不到我定的精度，就把NSW调整为了120。
赞一下 回复此楼
明察、慎思、笃行
7楼2011-03-02 10:25:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
120不可能达到这个标准啊，再多的步数也难说
赞一下(1人) 回复此楼
心之所向无惧无悔
8楼2011-03-02 10:29:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
bingmou(金币+2): 感谢参与 2011-03-02 15:14:24
VASP默认的是以能量作为收敛标准的，而你EDIFFG=-1e-4是表示力的收敛标准，前面说了这是相当相当高的标准。你可以改为EDIFFG= 1e-4。

[ Last edited by stractor on 2011-3-2 at 14:17 ]
赞一下(2人) 回复此楼
9楼2011-03-02 14:14:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫于

金虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
youzhizhe(金币+3): 谢谢交流。 2011-03-02 16:07:13
youzhizhe(金币+3): 谢谢交流。 2011-03-02 16:07:16
youzhizhe(金币+3): 谢谢交流。 2011-03-02 16:07:23
youzhizhe(金币+3): 谢谢交流。 2011-03-02 16:07:33
youzhizhe(金币+3): 谢谢交流。 2011-03-02 16:07:40
youzhizhe(金币+3): 谢谢交流。 2011-03-02 16:07:45
youzhizhe(金币+3): 谢谢交流。 2011-03-02 16:07:55
youzhizhe(金币+3): 谢谢交流。 2011-03-02 16:08:01
youzhizhe(金币+3): 谢谢交流。 2011-03-02 16:08:08
youzhizhe(金币+3): 谢谢交流。 2011-03-02 16:08:14
youzhizhe(金币-27): 由于网络错误导致的重复评贴，对此造成的误会，版主诚挚向您道歉。谢谢交流。 2011-03-02 16:35:56
你看看你的OUTCAR中是否有FORCES: max atom, RMS     0.268526    0.061907
这样一行，ediffg设为负值，应该是在最大原子受力小于你的设置时计算就结束了，如上我的最大原子受力为0.268526.    还有就是我用不同地方的机器算，有的有这一行，有的没有这样一行，不明白为什么。如果有这行，只要这个数字小于了你的设置ediffg（负值），计算就达到了精度，计算就结束。
赞一下(12人) 回复此楼
10楼2011-03-02 15:54:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gleerat 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭