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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 收敛精度
51/1返回列表
查看: 2659  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

bjmaowei25

铜虫 (小有名气)

[求助] 收敛精度

设置运行的NSW=1000, 可运行60步后自动停止，并告诉精度到达要求，但明显dE 很大没收敛，求
大家的点评，谢谢！
INCAR

SYSTEM   = M2O3
ENCUT    = 500
ISTART   = 0
ICHARG   = 2
ISMEAR   = 1
PREC     = M
EDIFF    = 1E-3
NSW      = 1000
ISIF     = 2
IBRION   = 2
IALGO    = 48
ALGO     = very fast
POTIM    = 0.1
LREAL    = Auto


OSZICAR
N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1     0.923112380146E+03    0.92311E+03   -0.50333E+04    74   0.357E+03
RMM:   2     0.175369114913E+03   -0.74774E+03   -0.10738E+04    74   0.394E+02
RMM:   3     0.174267355228E+03   -0.11018E+01   -0.15431E+03    74   0.342E+02
RMM:   4    -0.542524666320E+01   -0.17969E+03   -0.86607E+02    74   0.123E+02
RMM:   5    -0.373367180864E+02   -0.31911E+02   -0.29624E+02    74   0.151E+02
RMM:   6    -0.677436920725E+02   -0.30407E+02   -0.17804E+02    74   0.658E+01
RMM:   7    -0.800862495786E+02   -0.12343E+02   -0.74106E+01    74   0.524E+01
RMM:   8    -0.882359395166E+02   -0.81497E+01   -0.51421E+01    74   0.285E+01
RMM:   9    -0.955797020611E+02   -0.73438E+01   -0.58969E+01   248   0.225E+01
RMM:  10    -0.963124520811E+02   -0.73275E+00   -0.92937E+00   283   0.642E+00
RMM:  11    -0.964586535670E+02   -0.14620E+00   -0.14783E+00   277   0.327E+00
RMM:  12    -0.965482989444E+02   -0.89645E-01   -0.90834E-01   280   0.220E+00    0.216E+01
RMM:  13    -0.990880302596E+02   -0.25397E+01   -0.72834E+01   212   0.238E+01    0.466E+01
RMM:  14    -0.988030396415E+02    0.28499E+00   -0.71966E+01   192   0.442E+01    0.294E+01
RMM:  15    -0.901814394107E+02    0.86216E+01   -0.84060E+00   205   0.176E+01    0.184E+01
RMM:  16    -0.885813824920E+02    0.16001E+01   -0.98976E-01   235   0.470E+00    0.134E+01
RMM:  17    -0.906172378817E+02   -0.20359E+01   -0.11675E+00   233   0.458E+00    0.146E+01
RMM:  18    -0.916309571854E+02   -0.10137E+01   -0.14801E+00   211   0.617E+00    0.991E+00
RMM:  19    -0.889660004866E+02    0.26650E+01   -0.57077E+00   192   0.128E+01    0.581E+00
RMM:  20    -0.887204715274E+02    0.24553E+00   -0.11635E+00   211   0.564E+00    0.551E+00
RMM:  21    -0.864728355587E+02    0.22476E+01   -0.50053E-01   194   0.416E+00    0.479E+00
RMM:  22    -0.881407375740E+02   -0.16679E+01   -0.49187E-01   199   0.354E+00    0.315E+00
RMM:  23    -0.870375191558E+02    0.11032E+01   -0.34042E-01   159   0.352E+00    0.548E+00
RMM:  24    -0.900241807071E+02   -0.29867E+01   -0.71858E-01   153   0.443E+00    0.335E+00
RMM:  25    -0.871196437618E+02    0.29045E+01   -0.48867E-01   173   0.412E+00    0.515E+00
RMM:  26    -0.904592699743E+02   -0.33396E+01   -0.74835E-01   160   0.450E+00    0.380E+00
RMM:  27    -0.872593727356E+02    0.31999E+01   -0.52872E-01   159   0.449E+00    0.491E+00
RMM:  28    -0.908610841605E+02   -0.36017E+01   -0.67250E-01   155   0.461E+00    0.423E+00
RMM:  29    -0.876984696886E+02    0.31626E+01   -0.63891E-01   152   0.492E+00    0.464E+00
RMM:  30    -0.909369033777E+02   -0.32384E+01   -0.69348E-01   147   0.468E+00    0.460E+00
RMM:  31    -0.881459421380E+02    0.27910E+01   -0.85276E-01   150   0.536E+00    0.423E+00
RMM:  32    -0.906263470891E+02   -0.24804E+01   -0.68250E-01   148   0.472E+00    0.483E+00
RMM:  33    -0.884082531914E+02    0.22181E+01   -0.88297E-01   149   0.547E+00    0.392E+00
RMM:  34    -0.903375932469E+02   -0.19293E+01   -0.66028E-01   145   0.466E+00    0.494E+00
RMM:  35    -0.885481926372E+02    0.17894E+01   -0.89977E-01   148   0.547E+00    0.373E+00
RMM:  36    -0.901607324062E+02   -0.16125E+01   -0.67139E-01   145   0.464E+00    0.500E+00
RMM:  37    -0.886275118245E+02    0.15332E+01   -0.94080E-01   148   0.547E+00    0.362E+00
RMM:  38    -0.900544748007E+02   -0.14270E+01   -0.69192E-01   145   0.464E+00    0.503E+00
RMM:  39    -0.886695168003E+02    0.13850E+01   -0.94613E-01   149   0.549E+00    0.355E+00
RMM:  40    -0.899878387757E+02   -0.13183E+01   -0.68633E-01   146   0.463E+00    0.505E+00
RMM:  41    -0.886901460457E+02    0.12977E+01   -0.94424E-01   145   0.547E+00    0.350E+00
RMM:  42    -0.899453742520E+02   -0.12552E+01   -0.68720E-01   145   0.464E+00    0.507E+00
RMM:  43    -0.887034990119E+02    0.12419E+01   -0.95275E-01   145   0.545E+00    0.346E+00
RMM:  44    -0.899179872721E+02   -0.12145E+01   -0.69601E-01   147   0.464E+00    0.510E+00
RMM:  45    -0.887090983630E+02    0.12089E+01   -0.94088E-01   142   0.547E+00    0.342E+00
RMM:  46    -0.898958996580E+02   -0.11868E+01   -0.68101E-01   147   0.462E+00    0.512E+00
RMM:  47    -0.887119752044E+02    0.11839E+01   -0.93817E-01   142   0.543E+00    0.337E+00
RMM:  48    -0.898810921160E+02   -0.11691E+01   -0.68506E-01   149   0.464E+00    0.515E+00
RMM:  49    -0.887149032964E+02    0.11662E+01   -0.94807E-01   141   0.545E+00    0.333E+00
RMM:  50    -0.898688771465E+02   -0.11540E+01   -0.69777E-01   146   0.462E+00    0.518E+00
RMM:  51    -0.887154505611E+02    0.11534E+01   -0.94228E-01   141   0.544E+00    0.329E+00
RMM:  52    -0.898563194566E+02   -0.11409E+01   -0.68267E-01   151   0.464E+00    0.521E+00
RMM:  53    -0.887138792528E+02    0.11424E+01   -0.93836E-01   140   0.544E+00    0.324E+00
RMM:  54    -0.898450068695E+02   -0.11311E+01   -0.68500E-01   146   0.462E+00    0.524E+00
RMM:  55    -0.887150853364E+02    0.11299E+01   -0.95461E-01   140   0.543E+00    0.319E+00
RMM:  56    -0.898350463963E+02   -0.11200E+01   -0.69105E-01   148   0.462E+00    0.527E+00
RMM:  57    -0.887128917073E+02    0.11222E+01   -0.94546E-01   141   0.545E+00    0.314E+00
RMM:  58    -0.898213736552E+02   -0.11085E+01   -0.66544E-01   147   0.462E+00    0.529E+00
RMM:  59    -0.887064903571E+02    0.11149E+01   -0.93655E-01   142   0.542E+00    0.312E+00
RMM:  60    -0.898133215169E+02   -0.11068E+01   -0.66742E-01   148   0.459E+00
   1 F= -.89813322E+02 E0= -.89816821E+02  d E =-.898133E+02
curvature:   0.00 expect dE= 0.000E+00 dE for cont linesearch  0.000E+00
trial: gam= 0.00000 g(F)=  0.144E-58 g(S)=  0.000E+00 ort = 0.000E+00 (trialstep = 0.100E+01)
search vector abs. value=  0.100E-09
reached required accuracy - stopping structural energy minimisation
writing wavefunctions
回复此楼

» 猜你喜欢
追求真理
1楼 2012-02-28 17:00:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bjmaowei25

铜虫 (小有名气)
引用回帖:
楼: Originally posted by thinkllh at 2012-02-29 23:57:54:
你只是贴出来了一个离子步的过程，一个离子步电子自恰计算步数默认为60步。很显然，你的参数设置太粗糙了，导致体系在两个态之间震荡，没有收敛的趋势。建议如楼上，改ISMEAR＝0，即采用Gaussian smearing方法，S ...

谢谢你的建议，由于我初步的更改SIGMA值, 0..1，0.2，0.3，0.4，0.5，0.6
发现在0.6很快收敛，得到这样的结果
  free  energy   TOTEN  =       -88.553907 eV

  energy  without entropy=      -82.267707  energy(sigma->0) =      -85.410807

误差有点是不是有点大？
等会再精试一下！
energy  without entropy 之后的能量值，是指energy(sigma->0)？
这个值要在1meV~10meV,根据侯老师的入门指南上说。
赞一下 回复此楼
追求真理
7楼2012-03-01 03:24:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

uuv2010

荣誉版主 (职业作家)

优秀版主

【答案】应助回帖

感谢参与，应助指数 +1
liliangfang(金币+1): 多谢指教 2012-02-29 11:10:34
NSW=1000设置的是粒子优化的最大步数，现在你所关心的是电子优化，所以需要自己改变一下NELM,这个是电子优化的最大步数，不写的话则是默认值60
赞一下 回复此楼
2楼2012-02-29 10:26:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bjmaowei25

铜虫 (小有名气)
引用回帖:
楼: Originally posted by uuv2010 at 2012-02-29 10:26:37:
NSW=1000设置的是粒子优化的最大步数，现在你所关心的是电子优化，所以需要自己改变一下NELM,这个是电子优化的最大步数，不写的话则是默认值60

谢谢你的回复，我可以增加NELM设置，但dE 还是很大E+01左右，运行结果怎么会提示是收敛呢eached required accuracy - stopping structural energy minimisation？
赞一下 回复此楼
追求真理
3楼2012-02-29 10:46:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
liliangfang(金币+1): 多谢指教 2012-02-29 20:50:55
可以试试下面的：
SYSTEM   = M2O3
ENCUT    = 500
ISTART   = 0
ICHARG   = 2
ISMEAR   = 0; SIGMA = 0.1
PREC     = Accurate
EDIFF    = 1E-5
NSW      = 100
ISIF     = 2
IBRION   = 2
#IALGO    = 48
#ALGO     = very fast
POTIM    = 0.1
LREAL    = Auto
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

bjmaowei25
4楼2012-02-29 19:02:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 4/200 2026-02-23 06:46 by jsjzfl
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 w89i99eaeh 2026-02-22 4/200 2026-02-23 06:36 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:34 by w4l55oybr1
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:29 by w4l55oybr1
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 8/400 2026-02-23 06:24 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +6 3dfhjxgsh7 2026-02-22 8/400 2026-02-23 06:21 by w4l55oybr1
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 10/500 2026-02-23 04:58 by 5jlh3qtdvx
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 6/300 2026-02-23 02:08 by 5jlh3qtdvx
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +5 3dfhjxgsh7 2026-02-22 6/300 2026-02-23 02:04 by 5jlh3qtdvx
[教师之家] 版面费该交吗 +7 苹果在哪里 2026-02-22 8/400 2026-02-22 22:37 by otani
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +5 successhe 2026-02-16 6/300 2026-02-22 21:38 by 山西悬空寺空悬�
[基金申请] 面上可以超过30页吧？ +4 阿拉贡aragon 2026-02-22 4/200 2026-02-22 21:22 by 山西悬空寺空悬�
[教师之家] 为什么中国大学教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？ +5 欢乐颂叶蓁 2026-02-21 5/250 2026-02-22 21:15 by 山西悬空寺空悬�
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 6/300 2026-02-22 19:48 by w89i99eaeh
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 3dfhjxgsh7 2026-02-22 4/200 2026-02-22 16:52 by khieu8v8m0
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-22 3/150 2026-02-22 16:37 by khieu8v8m0
[基金申请] “人文社科而论，许多学术研究还没有达到民国时期的水平” +4 苏东坡二世 2026-02-18 5/250 2026-02-22 16:07 by liangep1573
[基金申请] 什么是人一生最重要的？ +4 瞬息宇宙 2026-02-21 4/200 2026-02-22 11:44 by huagongfeihu
[基金申请] 今年春晚有几个节目很不错，点赞！ +11 瞬息宇宙 2026-02-16 12/600 2026-02-21 21:14 by lq493392203
[基金申请] 体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低 +9 瞬息宇宙 2026-02-20 12/600 2026-02-21 10:39 by 欢乐颂叶蓁
信息提示
关闭
请填处理意见
关闭