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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 收敛精度
81/1返回列表
查看: 2197  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者bjmaowei25将赠送您 2 个金币

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的回帖
回帖置顶 ( 共有4个 )

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

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

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

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

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

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 的主题更新
81/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 请问大家，我有希望进浙大普博么 +5 高艺文 2024-09-28 8/400 2024-09-29 00:08 by 小红豆
[论文投稿] applied catalysis B 投稿模版 5+4 午后夏日 2024-09-27 5/250 2024-09-28 22:59 by 莱茵润色
[硕博家园] 希望有师兄师姐帮忙引荐申博 +6 zhanyaqian 2024-09-25 11/550 2024-09-28 22:24 by 自强2019
[教师之家] 东南大学电气学院电力电子方向招收硕士博士 +4 蜡笔小鑫1989 2024-09-26 7/350 2024-09-28 21:16 by 伍文童
[基金申请] QB怎么还没消息 +21 wuxiacl 2024-09-24 57/2850 2024-09-28 14:48 by xiaoxiao270
[考博] 25有机合成申博 +3 再等你啊 2024-09-24 4/200 2024-09-28 13:09 by majun521
[论文投稿] 想跟大家聊聊对国产SCI期刊的看法 +14 wleizl 2024-09-26 17/850 2024-09-28 11:59 by wleizl
[教师之家] 西湖大学教授：我可以自由地选择讲课的方式, 讲授的角度 +4 zju2000 2024-09-22 4/200 2024-09-28 09:40 by bio-polymer
[电化学] 刚装好的电池电压掉很快怎么解决？？ +3 花花你最棒 2024-09-26 6/300 2024-09-28 08:20 by cnofc
[考博] 材料/电信/生物-2025普博生自荐-985本双非硕一区一作 +4 enowei0127 2024-09-23 8/400 2024-09-27 22:40 by enowei0127
[育儿交流] 熊孩子睡觉前很多事，啥原因？ +7 quan2153 2024-09-25 14/700 2024-09-27 16:46 by quan2153
[基金申请] 请问大家的计划书填写列表中状态更新了吗？ +7 Laker610 2024-09-25 9/450 2024-09-27 16:36 by 田田hj
[教师之家] 压力山大 +17 暮色恋伊人 2024-09-24 25/1250 2024-09-27 16:10 by 领先于123
[考博] 电催化，电池方向学生自荐 30+3 Freelancedd 2024-09-24 10/500 2024-09-27 10:27 by 莱茵润色
[论文投稿] 期刊论文发表了但查询到sci未收录 50+5 yibuxiao 2024-09-25 11/550 2024-09-27 10:24 by 莱茵润色
[论文投稿] 为什么审稿人一直拒绝审稿啊？ +7 潇湘雪月 2024-09-24 13/650 2024-09-27 07:53 by weizhi111
[考博] 有机转材料申博难度 30+3 昕散 2024-09-22 10/500 2024-09-26 22:42 by 824282658
[论文投稿] NC投稿多久有消息呀？ +4 Jordanblood 2024-09-24 6/300 2024-09-25 16:58 by FZX_2024
[有机交流] chemdraw结构式复制到word变形 +5 笑看人生1993 2024-09-22 11/550 2024-09-25 14:44 by 冰蓝夜游神
[论文投稿] 爱斯维尔旗下Postharvest Biology and Technology，采用什么系统查查 20+4 sdsdsxh 2024-09-22 9/450 2024-09-25 07:33 by sdsdsxh
信息提示
关闭
请填处理意见
关闭