24小时热门版块排行榜    

查看: 3754  |  回复: 40

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 第一性计算

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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的关系,哪一个到了收敛精度都会被认为收敛了,你可以把能量的收敛标准再提高一个数量级试试?
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
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 ]
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可以考虑你的模型的合理性
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 4/200 2026-02-08 08:46 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 4/200 2026-02-08 08:27 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:12 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:07 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 08:06 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:52 by vs90ilomwc
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:46 by vs90ilomwc
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:32 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:27 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:26 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:07 by vs90ilomwc
[硕博家园] 博士延得我,科研能力直往上蹿 +8 偏振片 2026-02-02 8/400 2026-02-08 06:52 by liyeqik
[教师之家] 有院领导为了换新车,用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[基金申请] 同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生,欢迎大家推荐交流(博导是本人) +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生|MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
请填处理意见