版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

bjmaowei25

銅蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 362.1
紅花: 1
帖子: 121
在線: 293.5小時(shí)
蟲號(hào): 1531554
注冊(cè): 2011-12-10
性別: GG
專業(yè): 結(jié)構(gòu)陶瓷

[求助] 收斂精度

設(shè)置運(yùn)行的NSW=1000, 可運(yùn)行60步后自動(dòng)停止，并告訴精度到達(dá)要求，但明顯dE 很大沒(méi)收斂，求
大家的點(diǎn)評(píng)，謝謝！
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

回復(fù)此樓

» 猜你喜歡

268求調(diào)劑已經(jīng)有3人回復(fù)
302求調(diào)劑已經(jīng)有5人回復(fù)
311求調(diào)劑已經(jīng)有10人回復(fù)
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復(fù)
303求調(diào)劑已經(jīng)有5人回復(fù)
考研求調(diào)劑已經(jīng)有4人回復(fù)
085601專碩，總分342求調(diào)劑，地區(qū)不限已經(jīng)有4人回復(fù)
085601材料工程專碩求調(diào)劑已經(jīng)有5人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有6人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有6人回復(fù)

追求真理

1樓 2012-02-28 17:00:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖置頂 ( 共有4個(gè) )

bjmaowei25

銅蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 362.1
紅花: 1
帖子: 121
在線: 293.5小時(shí)
蟲號(hào): 1531554
注冊(cè): 2011-12-10
性別: GG
專業(yè): 結(jié)構(gòu)陶瓷

引用回帖:

樓: Originally posted by uuv2010 at 2012-02-29 10:26:37:
NSW=1000設(shè)置的是粒子優(yōu)化的最大步數(shù)，現(xiàn)在你所關(guān)心的是電子優(yōu)化，所以需要自己改變一下NELM,這個(gè)是電子優(yōu)化的最大步數(shù)，不寫的話則是默認(rèn)值60

謝謝你的回復(fù)，我可以增加NELM設(shè)置，但dE 還是很大E+01左右，運(yùn)行結(jié)果怎么會(huì)提示是收斂呢eached required accuracy - stopping structural energy minimisation？

贊一下

回復(fù)此樓

追求真理

3樓2012-02-29 10:46:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bjmaowei25

銅蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 362.1
紅花: 1
帖子: 121
在線: 293.5小時(shí)
蟲號(hào): 1531554
注冊(cè): 2011-12-10
性別: GG
專業(yè): 結(jié)構(gòu)陶瓷

送鮮花一朵

引用回帖:

樓: 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 ...

謝謝侯老師！

贊一下

回復(fù)此樓

追求真理

5樓2012-02-29 19:55:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bjmaowei25

銅蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 362.1
紅花: 1
帖子: 121
在線: 293.5小時(shí)
蟲號(hào): 1531554
注冊(cè): 2011-12-10
性別: GG
專業(yè): 結(jié)構(gòu)陶瓷

引用回帖:

樓: Originally posted by thinkllh at 2012-02-29 23:57:54:
你只是貼出來(lái)了一個(gè)離子步的過(guò)程，一個(gè)離子步電子自恰計(jì)算步數(shù)默認(rèn)為60步。很顯然，你的參數(shù)設(shè)置太粗糙了，導(dǎo)致體系在兩個(gè)態(tài)之間震蕩，沒(méi)有收斂的趨勢(shì)。建議如樓上，改ISMEAR＝0，即采用Gaussian smearing方法，S ...

謝謝你的建議，由于我初步的更改SIGMA值, 0..1，0.2，0.3，0.4，0.5，0.6
發(fā)現(xiàn)在0.6很快收斂，得到這樣的結(jié)果
  free  energy TOTEN  =    -88.553907 eV

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

誤差有點(diǎn)是不是有點(diǎn)大？
等會(huì)再精試一下！
energy  without entropy 之后的能量值，是指energy(sigma->0)？
這個(gè)值要在1meV~10meV,根據(jù)侯老師的入門指南上說(shuō)。

贊一下

回復(fù)此樓

追求真理

7樓2012-03-01 03:24:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bjmaowei25

銅蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 362.1
紅花: 1
帖子: 121
在線: 293.5小時(shí)
蟲號(hào): 1531554
注冊(cè): 2011-12-10
性別: GG
專業(yè): 結(jié)構(gòu)陶瓷

引用回帖:

侯老師，你好！

由于收斂的比較慢，我初步的更改SIGMA值, 0..1，0.2，0.3，0.4，0.5，0.6
發(fā)現(xiàn)在0.6很快收斂，得到這樣的結(jié)果
  free  energy TOTEN  =    -88.553907 eV
  energy  without entropy=    -82.267707  energy(sigma->0) =    -85.410807

根據(jù)你的入門指南談到， OUTCAR 中 energy  without entropy 之后的能量值，是指energy(sigma->0)？
這個(gè)值要在1meV~10meV。
還是指兩者之差在1meV~10meV，也就是 energy(sigma->0)  減去 energy  without entropy 呢？

謝謝了!

還是根據(jù)侯老師的入門指南上說(shuō)。

贊一下

回復(fù)此樓

追求真理

8樓2012-03-01 03:29:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

uuv2010

榮譽(yù)版主 (職業(yè)作家)

1ST強(qiáng)帖: 6
應(yīng)助: 160 (高中生)
貴賓: 3.14
金幣: 3681.7
散金: 5802
紅花: 79
沙發(fā): 35
帖子: 4560
在線: 916小時(shí)
蟲號(hào): 966216
注冊(cè): 2010-03-09
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1
liliangfang(金幣+1): 多謝指教 2012-02-29 11:10:34

NSW=1000設(shè)置的是粒子優(yōu)化的最大步數(shù)，現(xiàn)在你所關(guān)心的是電子優(yōu)化，所以需要自己改變一下NELM,這個(gè)是電子優(yōu)化的最大步數(shù)，不寫的話則是默認(rèn)值60

贊一下

回復(fù)此樓

2樓2012-02-29 10:26:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲 (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線: 873.8小時(shí)
蟲號(hào): 1007127
注冊(cè): 2010-04-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +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

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

bjmaowei25

4樓2012-02-29 19:02:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

thinkllh

鐵桿木蟲 (正式寫手)

應(yīng)助: 34 (小學(xué)生)
金幣: 6569.1
散金: 50
紅花: 1
帖子: 538
在線: 157小時(shí)
蟲號(hào): 715800
注冊(cè): 2009-03-05
性別: GG
專業(yè): 極端條件下使用的金屬材料

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1
WDD880227(金幣+1): 感謝交流，很有幫助 2012-03-01 08:28:54

你只是貼出來(lái)了一個(gè)離子步的過(guò)程，一個(gè)離子步電子自恰計(jì)算步數(shù)默認(rèn)為60步。很顯然，你的參數(shù)設(shè)置太粗糙了，導(dǎo)致體系在兩個(gè)態(tài)之間震蕩，沒(méi)有收斂的趨勢(shì)。建議如樓上，改ISMEAR＝0，即采用Gaussian smearing方法，SIGMA值的設(shè)置對(duì)收斂與否及收斂的速度影響很大，可從0.1測(cè)試，若仍不收斂，增大SIGMA，繼續(xù)測(cè)試，祝好運(yùn)。

贊一下

回復(fù)此樓

Thinkwiththeworld.

6樓2012-02-29 23:57:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 bjmaowei25 的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 268求調(diào)劑 +3	簡(jiǎn)單點(diǎn)0 2026-03-17	3/150	2026-03-17 23:15 by wangjihu
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	7/350	2026-03-17 20:20 by 花125533
[考研] 301求調(diào)劑 +5	yy要上岸呀 2026-03-17	5/250	2026-03-17 20:20 by peike
[考研] 08工科 320總分求調(diào)劑 +4	梨花珞晚風(fēng) 2026-03-17	4/200	2026-03-17 13:38 by houyaoxu
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6	SUICHILD 2026-03-12	6/300	2026-03-17 09:24 by 霧散后相遇lc
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 085600材料與化工求調(diào)劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級(jí)475 +5	Becho359 2026-03-11	5/250	2026-03-14 11:35 by 哦哦123
[考研] 0703化學(xué)調(diào)劑 +4	快樂(lè)的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 求材料調(diào)劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655