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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 为什么力没有达到收敛标准，就提示已经达到了收敛精度，并且停止了呢？
53/1返回列表
查看: 3665  |  回复: 40
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

c03185

新虫 (小有名气)

[求助] 为什么力没有达到收敛标准，就提示已经达到了收敛精度，并且停止了呢？

我的EDIFFG = -2E-2，为什么力没有达到收敛标准，就提示已经达到了收敛精度(reached required accuracy - stopping structural energy minimisation)，并且停止了呢？
POSCAR是对每个原子的进行（F F T）驰豫。

这是我的INCAR：
SYSTEM = CuO
ENCUT = 520
ISTART = 0; ICHARG = 2
ISMEAR = 1; SIGMA = 0.20
NSW = 600; IBRION = 2
POTIM = 0.1
EDIFF = 1E-5; EDIFFG = -2E-2
PREC= Accurate
IALGO=48

这是收敛后OUTCAR中的力的情况：
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.50768     -0.85473      0.11467        -0.669668     -0.711244      0.550448
      4.61607      3.84534      3.38384         0.428006      0.878143     -0.550448
      3.18038      1.48121     -0.65192         1.032349     -0.325156     -0.292996
      3.28877      6.18128      2.61724        -1.032349      0.066748      0.432572
      1.85308     -0.44066      0.88127        -0.428006      0.941812     -0.432572
      1.96147      4.25940      4.15044         0.669668     -0.850303      0.292996
      4.03164     -0.41006      2.97621         0.320228     -0.221158      0.000000
      1.11022      4.73140      1.28891         0.000000      0.342418     -0.184952
      2.43752      2.11015      2.20961        -0.320228     -0.121260      0.184952
      3.63869      0.72555      1.34328        -0.012267     -0.017762     -0.032885
      3.59276      5.20217      4.19875         0.002868      0.004152      0.007687
      2.87640     -0.37822     -0.70023        -0.002868     -0.004152     -0.007687
      2.83046      4.09840      2.15524         0.012267      0.017762      0.032885
      0.66516      0.96312      1.78311        -0.000486     -0.000704     -0.001304
      1.55528      2.25198      4.16930         0.000486      0.000704      0.001304
回复此楼

» 猜你喜欢
1楼 2012-08-05 07:47:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c03185

新虫 (小有名气)
送鲜花一朵
引用回帖:
14楼: Originally posted by cenwanglai at 2012-08-05 21:02:08
分析一下你的几个原子最后最后一列的受力数据，你看看。受力是成对的，一个正A，一个负A，两相抵消恰好是0.

为什么不试着把某层原子全部固定呢？

能把POSCAR贴出来不？也可以直接发给我。
cwlzju@zju.edu.cn....

结果还是那样，找不到原因。同样的参数用在其他的晶格结构里边，驰豫原子对应的力是收敛的。在这个晶格结构就是不对。版主爷爷给指条活路吧。。。

INCAR：
SYSTEM = AO
ENCUT = 520
ISTART = 0; ICHARG = 2
ISMEAR = 1; SIGMA = 0.20
NSW = 600; IBRION = 2
POTIM = 0.1
EDIFF = 1E-7
EDIFFG = -5E-2
PREC= Accurate
IALGO=48

POSCAR：
AO(001)_9
   1.00000000000000
     4.2487171130865562   -2.9342791501140022    0.0000000000000010
     0.0000000000000000    4.5431245416429560   -2.4538970306330690
     6.8981356599745016    9.9882204546856510   18.4921082297066661
   O    Ti
   9   6
Selective dynamics
Direct
  0.9790666575233118  0.3333333333333357  0.0505438132934409   F   F   T
  0.6457333241899832  0.6666666666666643  0.2713454927088639   F   F   T
  0.6666666666666643  0.6457333241899832  0.0505438132934409   F   F   T
  0.3333333333333357  0.9790666575233118  0.2713454927088639   F   F   T
  0.3542666758100168  0.0209333424766882  0.0505438132934409   F   F   T
  0.0209333424766882  0.3542666758100168  0.2713454927088639   F   F   T
  0.6876000091433525  0.0000000000000000  0.1609446530011525   F   F   F
  0.0000000000000000  0.6876000091433525  0.1609446530011525   F   F   F
  0.3123999908566475  0.3123999908566475  0.1609446530011525   F   F   F
  0.6666666666666643  0.3333333333333357  0.1168738138540851   F   F   F
  0.3333333333333357  0.6666666666666643  0.3154820665970479   F   F   T
  0.6666666666666643  0.3333333333333357  0.0064072394052569   F   F   T
  0.3333333333333357  0.6666666666666643  0.2050154921482271   F   F   F
  0.0000000000000000  0.0000000000000000  0.0963564137588521   F   F   T
  0.0000000000000000  0.0000000000000000  0.2255328922434601   F   F   T


最终结果:
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.50830     -0.85384      0.11633        -0.671847     -0.681135      0.568972
      4.61545      3.84444      3.38219         0.399076      0.869518     -0.568972
      3.18100      1.48211     -0.65027         1.032169     -0.327554     -0.257660
      3.28815      6.18039      2.61559        -1.032169      0.035882      0.415202
      1.85370     -0.43977      0.88293        -0.399076      0.952574     -0.415202
      1.96085      4.25851      4.14878         0.671847     -0.849285      0.257660
      4.03164     -0.41006      2.97621         0.308830     -0.213287      0.000000
      1.11022      4.73140      1.28891         0.000000      0.330230     -0.178369
      2.43752      2.11015      2.20961        -0.308830     -0.116944      0.178369
      3.63869      0.72555      1.34328        -0.005537     -0.008018     -0.014844
      3.59252      5.20183      4.19812         0.005220      0.007559      0.013995
      2.87663     -0.37788     -0.69961        -0.005220     -0.007559     -0.013995
      2.83046      4.09840      2.15524         0.005537      0.008018      0.014844
      0.66483      0.96265      1.78224         0.011722      0.016973      0.031424
      1.55561      2.25245      4.17017        -0.011722     -0.016973     -0.031424
-----------------------------------------------------------------------------------
    total drift:                                0.000075      0.000108      0.000200
d Force = 0.1619927E-03[ 0.151E-03, 0.173E-03]  d Energy =-0.4911406E-03 0.653E-03
d Force = 0.2019964E+01[ 0.202E+01, 0.202E+01]  d Ewald  = 0.2019964E+01-0.245E-07
--------------------------------------------------------------------------------------------------------
    POTLOK:  cpu time    1.34: real time    0.69
--------------------------------------------------------------------------------------------------------
Steepest descent step on ions:
trial-energy change:    0.000491  1 .order   -0.000162   -0.000173   -0.000151
  (g-gl).g = 0.173E-03      g.g   = 0.173E-03  gl.gl    = 0.000E+00
g(Force)  = 0.173E-03   g(Stress)= 0.000E+00 ortho     = 0.000E+00
gamma     =   0.00000
trial     =   1.00000
opt step  =   4.00000  (harmonic =   7.63264) maximal distance =0.00162217
next E    =  -127.127400   (d E  =  -0.00066)
--------------------------------------------------------------------------------------------------------
reached required accuracy - stopping structural energy minimisation
writing wavefunctions
赞一下 回复此楼
22楼2012-08-06 01:04:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 41 个回答

itpfeng

木虫 (著名写手)
★
xueht987: 金币+1, 鼓励交流！ 2012-08-05 11:44:42
大概是能量的收敛已经到给定的精度了吧
赞一下 回复此楼
2楼2012-08-05 09:08:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c03185

新虫 (小有名气)
引用回帖:
2楼: Originally posted by itpfeng at 2012-08-05 09:08:26
大概是能量的收敛已经到给定的精度了吧

可是力没有达到收敛的精度，为什么就停止了呢？
赞一下 回复此楼
3楼2012-08-05 15:59:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

itpfeng

木虫 (著名写手)
★ ★
franch: 金币+2, 谢谢回帖交流。 2012-08-05 22:08:09
引用回帖:
3楼: Originally posted by c03185 at 2012-08-05 15:59:47
可是力没有达到收敛的精度，为什么就停止了呢？...

力和能量大概是OR而不是AND的关系，哪一个到了收敛精度都会被认为收敛了，你可以把能量的收敛标准再提高一个数量级试试？
赞一下(1人) 回复此楼
4楼2012-08-05 16:07:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 41 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭