24小时热门版块排行榜    

查看: 3802  |  回复: 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
回复此楼

» 猜你喜欢

已阅   回复此楼   关注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的关系,哪一个到了收敛精度都会被认为收敛了,你可以把能量的收敛标准再提高一个数量级试试?
4楼2012-08-05 16:07:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料277求调剂 +5 min3 2026-03-24 5/250 2026-03-26 15:13 by zzll406
[考研] 总分293求调剂 +6 加一一九 2026-03-25 8/400 2026-03-26 13:30 by yujianx
[考研] 316求调剂 +9 梁茜雯 2026-03-19 9/450 2026-03-25 23:01 by peike
[考研] 08工学调剂 +13 用户573181 2026-03-20 18/900 2026-03-25 22:00 by zbssa
[考研] 26考研-291分-厦门大学(085601)-柔性电子学院材料工程专业求调剂 +3 min3 2026-03-24 4/200 2026-03-25 18:22 by xcjcqu
[考研] 招08考数学 +8 laoshidan 2026-03-20 17/850 2026-03-25 17:52 by 一个红太阳
[考研] 329求调剂 +3 钮恩雪 2026-03-25 3/150 2026-03-25 14:43 by 糖加冰
[考研] 085602 289分求调剂 +7 WWW西西弗斯 2026-03-24 7/350 2026-03-25 14:28 by 3Strings
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 11/550 2026-03-25 10:23 by userper
[考研] 085601求调剂总分293英一数二 +3 钢铁大炮 2026-03-24 3/150 2026-03-24 22:03 by bingxueer79
[考研] 一志愿南航材料专317分求调剂 +5 炸呀炸呀炸薯条 2026-03-23 5/250 2026-03-24 16:52 by 星空星月
[考研] 300求调剂,材料科学英一数二 +5 leaflight 2026-03-24 5/250 2026-03-24 16:25 by laoshidan
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
信息提示
请填处理意见