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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 为什么力没有达到收敛标准，就提示已经达到了收敛精度，并且停止了呢？
411/5返回列表12345下一页
查看: 3616  |  回复: 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
回复此楼

» 收录本帖的淘帖专辑推荐

DFT 第一性计算

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-08-05 07:47:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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的回帖

c03185

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

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

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

谢谢区长爷爷现身，我按照您说的只驰豫上3层和下3层，结果还是那样的。这是我的POSCAR：
AO(001)
4.2487171130865562   -2.9342791501140022    0.0000000000000010
     0.0000000000000000    4.5431245416429560   -2.4538970306330690
     6.8981356599745016    9.9882204546856510   18.4921082297066661
   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
赞一下(1人) 回复此楼
16楼2012-08-05 23:00:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
★ ★
liliangfang: 金币+2, 谢谢指教 2012-08-07 09:20:59
怀疑问题不是出在算法上，而是作业输入文件有可能写错了，建议认真检查输入参数和输入文件
如果书写正确的话，不应该出现这样的问题
EDIFFG is negative it has a different meaning: In this case the relaxation will stop if all forces are smaller than |EDIFFG |
如果可以的话，可以将INCAR,POSCAR和KPOINTS的内容都贴上来，让大家一起找找原因

[ Last edited by uuv2010 on 2012-8-6 at 11:33 ]
赞一下(1人) 回复此楼
26楼2012-08-06 11:32:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

成广

铜虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fzx2008: 金币+2, 谢谢交流！ 2012-08-17 11:19:17
c03185: 金币+10, ★★★很有帮助 2012-08-18 02:25:52
vasp的官方论坛有讨论过相关问题：
http://cms.mpi.univie.ac.at/vasp-forum/forum_viewtopic.php?4.4717
这个pdf的第28也有解决办法：
http://cms.mpi.univie.ac.at/vasp-workshop/slides/accuracy.pdf
但是一般来说，如果你是受限弛豫的话，一般不让动的维度力常常达不到精度的，lz可以考虑你的模型的合理性
赞一下(1人) 回复此楼
38楼2012-08-16 23:29:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖

c03185

新虫 (小有名气)
引用回帖:
4楼: Originally posted by itpfeng at 2012-08-05 16:07:21
力和能量大概是OR而不是AND的关系，哪一个到了收敛精度都会被认为收敛了，你可以把能量的收敛标准再提高一个数量级试试？...

还是不行阿
回复此楼
5楼2012-08-05 16:28:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c03185

新虫 (小有名气)
继续求助阿。同样的参数用在其它的晶格结构里边，力可以达到收敛的标准，能量也达到了收敛的标准。
可用在这个结构里，力就是始终没有达到收敛标准，就提示收敛了。
继续求帮助。。。。
谢谢阿谢谢
赞一下 回复此楼
6楼2012-08-05 16:31:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c03185

新虫 (小有名气)
引用回帖:
6楼: Originally posted by c03185 at 2012-08-05 16:31:26
继续求助阿。同样的参数用在其它的晶格结构里边，力可以达到收敛的标准，能量也达到了收敛的标准。
可用在这个结构里，力就是始终没有达到收敛标准，就提示收敛了。
继续求帮助。。。。
谢谢阿谢谢

换了赝势还是不行
赞一下 回复此楼
7楼2012-08-05 18:22:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c03185

新虫 (小有名气)
引用回帖:
6楼: Originally posted by c03185 at 2012-08-05 16:31:26
继续求助阿。同样的参数用在其它的晶格结构里边，力可以达到收敛的标准，能量也达到了收敛的标准。
可用在这个结构里，力就是始终没有达到收敛标准，就提示收敛了。
继续求帮助。。。。
谢谢阿谢谢

ENCUT，KPOINTS都修改了，结果还是一样，力没有收敛就提示收敛了
赞一下 回复此楼
8楼2012-08-05 18:22:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangll99

木虫 (职业作家)
★ ★
franch: 金币+2, 谢谢回帖交流。 2012-08-05 22:08:20
收敛标准看的是你选的为T的自由度，F的可以大于收敛标准

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
9楼2012-08-05 18:30:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

VASP_2012

新虫 (小有名气)
★
franch: 金币+1, 谢谢回帖交流。 2012-08-05 22:08:29
引用回帖:
9楼: Originally posted by huangll99 at 2012-08-05 18:30:48
收敛标准看的是你选的为T的自由度，F的可以大于收敛标准

谢谢，可是T的也没有达到呀
赞一下 回复此楼
10楼2012-08-05 18:37:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 c03185 的主题更新
411/5返回列表12345下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭