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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 收敛精度
北京石油化工学院2026年研究生招生接收调剂公告
81/1返回列表
查看: 2715  |  回复: 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的回帖

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

bjmaowei25

铜虫 (小有名气)
送鲜花一朵
引用回帖:
楼: Originally posted by valenhou001 at 2012-02-29 19:02:18:
可以试试下面的：
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    ...

谢谢侯老师！
回复此楼
追求真理
5楼2012-02-29 19:55:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

thinkllh

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
WDD880227(金币+1): 感谢交流，很有帮助 2012-03-01 08:28:54
你只是贴出来了一个离子步的过程，一个离子步电子自恰计算步数默认为60步。很显然，你的参数设置太粗糙了，导致体系在两个态之间震荡，没有收敛的趋势。建议如楼上，改ISMEAR＝0，即采用Gaussian smearing方法，SIGMA值的设置对收敛与否及收敛的速度影响很大，可从0.1测试，若仍不收敛，增大SIGMA，继续测试，祝好运。
赞一下 回复此楼
Thinkwiththeworld.
6楼2012-02-29 23:57:54
已阅   回复此楼   关注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的回帖

bjmaowei25

铜虫 (小有名气)
引用回帖:
楼: Originally posted by valenhou001 at 2012-02-29 19:02:18:
可以试试下面的：
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    ...

侯老师，你好！

由于收敛的比较慢，我初步的更改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

根据你的入门指南谈到， OUTCAR 中 energy  without entropy 之后的能量值，是指energy(sigma->0)？
这个值要在1meV~10meV。
还是指两者之差在1meV~10meV， 也就是 energy(sigma->0)  减去 energy  without entropy 呢？

谢谢了!

还是根据侯老师的入门指南上说。
赞一下 回复此楼
追求真理
8楼2012-03-01 03:29:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bjmaowei25 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 334分 一志愿武理 材料求调剂 +7 李李不服输 2026-03-26 7/350 2026-03-29 16:25 by 学员8dgXkO
[考研] 085600 286分 材料求调剂 +8 麻辣鱿鱼 2026-03-27 9/450 2026-03-29 15:36 by Serene1974
[考研] 311求调剂 +5 冬十三 2026-03-24 5/250 2026-03-29 08:55 by qingfeng258
[考研] 315求调剂 +4 akie... 2026-03-28 5/250 2026-03-28 21:05 by zhq0425
[考研] 299求调剂 +7 嗯嗯嗯嗯2 2026-03-27 7/350 2026-03-28 13:09 by 唐沐儿
[考研] 283求调剂 +7 A child 2026-03-28 7/350 2026-03-28 12:05 by zllcz
[考研] 一志愿南京航空航天大学材料学硕求调剂 +3 @taotao 2026-03-28 3/150 2026-03-28 10:26 by JourneyLucky
[考研] 291求调剂 +6 HanBeiNingZC 2026-03-24 6/300 2026-03-28 07:55 by baoball
[考研] 一志愿上海理工能源动力（085800）310分求调剂 +3 zhangmingc 2026-03-27 4/200 2026-03-27 19:01 by 给你你注意休息
[考研] 322求调剂 +4 我真的很想学习 2026-03-23 4/200 2026-03-27 13:51 by 杨杨杨紫
[考研] 一志愿211，335分，0856，求调剂院校和导师 +4 倾____萧 2026-03-27 5/250 2026-03-27 11:52 by zhshch
[考研] 304材料求调剂 +4 钟llll 2026-03-26 4/200 2026-03-27 03:42 by wxiongid
[考研] 341求调剂 +7 青柠檬1 2026-03-26 7/350 2026-03-27 00:19 by wxiongid
[考研] 327求调剂 +7 prayer13 2026-03-23 7/350 2026-03-26 20:48 by 不吃魚的貓
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-26 3/150 2026-03-26 19:57 by nihaoar
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 309求调剂 +4 gajsj 2026-03-25 5/250 2026-03-26 00:27 by Dyhoer
[考研] 318求调剂 +3 plum李子 2026-03-23 3/150 2026-03-25 09:42 by 雾散后相遇lc
[考研] 一志愿南航材料专317分求调剂 +5 炸呀炸呀炸薯条 2026-03-23 5/250 2026-03-24 16:52 by 星空星月
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
信息提示
关闭
请填处理意见
关闭