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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

rainxin2008

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1718.2
散金: 120
帖子: 265
在線: 163.3小時
蟲號: 1619644
注冊: 2012-02-15
專業(yè): 金屬功能材料

引用回帖:

18樓: Originally posted by mousekingadv at 2012-08-12 19:31:31
表面的晶格大小是固定的，弛豫的是離子位置...

你說的我不太懂，我的是最下面一層固定F，其余的全是放松的T，師姐們的都是ISIF=0。不管我改ISIF=0還是2，結(jié)果都是不收斂，你能幫幫我嗎？謝謝你

贊一下

回復(fù)此樓

21樓2012-08-12 19:41:47

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

mousekingadv

金蟲 (小有名氣)

應(yīng)助: 51 (初中生)
金幣: 1141.8
紅花: 1
帖子: 140
在線: 28.7小時
蟲號: 1721002
注冊: 2012-03-27
專業(yè): 理論和計算化學(xué)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
WDD880227: 金幣+2, 多謝指教 2012-08-13 19:04:45
rainxin2008: 金幣+10, ★有幫助 2012-08-14 09:21:10

引用回帖:

21樓: Originally posted by rainxin2008 at 2012-08-12 19:41:47
你說的我不太懂，我的是最下面一層固定F，其余的全是放松的T，師姐們的都是ISIF=0。不管我改ISIF=0還是2，結(jié)果都是不收斂，你能幫幫我嗎？謝謝你...

IBRION=2在搜尋能量最低點的時候是以下過程：
計算原子受力，乘以POTIM后讓離子移動一定的距離，
{再次計算受力、能量；
if 當(dāng)前位置是這條路徑上的最低點->計算原子受力，乘以POTIM后讓離子移動一定的距離
else 根據(jù)當(dāng)前受力，估計可能的最低點（僅限在此一維路徑上繼續(xù)跑）
}
你所遇到的錯誤是因為計算誤差或最低點附近能量曲線過于平坦之類的原因無法找到此路徑上的最低點。如果確認(rèn)輸入文件沒有問題的話，沒有關(guān)系，照他說的從當(dāng)前CONTCAR繼續(xù)就行了（即停止尋找最低點，直接跳出循環(huán)外找下一個路徑）
如果反復(fù)出現(xiàn)這樣的問題，不妨換IBRION=1試試。

贊一下

回復(fù)此樓

22樓2012-08-13 01:20:32

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

rainxin2008

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1718.2
散金: 120
帖子: 265
在線: 163.3小時
蟲號: 1619644
注冊: 2012-02-15
專業(yè): 金屬功能材料

引用回帖:

22樓: Originally posted by mousekingadv at 2012-08-13 01:20:32
IBRION=2在搜尋能量最低點的時候是以下過程：
計算原子受力，乘以POTIM后讓離子移動一定的距離，
{再次計算受力、能量；
if 當(dāng)前位置是這條路徑上的最低點->計算原子受力，乘以POTIM后讓離子移動一定的距離 ...

恩，謝謝你，我把IBRION改為1了，ISIF=0.現(xiàn)在在重算，看結(jié)果能不能收斂了。。。

贊一下

回復(fù)此樓

23樓2012-08-13 09:29:16

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

rainxin2008

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1718.2
散金: 120
帖子: 265
在線: 163.3小時
蟲號: 1619644
注冊: 2012-02-15
專業(yè): 金屬功能材料

引用回帖:

你好，我把IBRION改為1算的結(jié)果如下： 195 F= -.17800925E+03 E0= -.17800569E+03  d E =-.497371E-02
BRION: g(F)=  0.255E-01 g(S)=  0.000E+00 retain N=  4 mean eig= 1.98
eig: 7.417  0.028  0.240  0.240
bond charge predicted
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.178021946329E+03 -0.12684E-01 -0.37661E+00  3273 0.301E+00 0.898E-02
RMM: 2 -0.178025353881E+03 -0.34076E-02 -0.38320E-02  3261 0.298E-01 0.686E-02
RMM: 3 -0.178025396033E+03 -0.42152E-04 -0.13120E-03  3267 0.684E-02 0.635E-02
RMM: 4 -0.178025320041E+03 0.75992E-04 -0.10529E-04  1793 0.188E-02
196 F= -.17802532E+03 E0= -.17802177E+03  d E =-.160742E-01
BRION: g(F)=  0.644E+00 g(S)=  0.000E+00 retain N=  5 mean eig= 1.59
eig: 7.444  0.235  0.235  0.022  0.022
bond charge predicted
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.178025456026E+03 -0.59993E-04 -0.41268E-03  3225 0.988E-02 0.842E-03
RMM: 2 -0.178025458590E+03 -0.25644E-05 -0.44000E-05  1737 0.115E-02 0.621E-03
RMM: 3 -0.178025457360E+03 0.12308E-05 -0.97267E-06  1648 0.719E-03
197 F= -.17802546E+03 E0= -.17802193E+03  d E =-.137319E-03
BRION: g(F)=  0.143E-01 g(S)=  0.000E+00 retain N=  1 mean eig= 0.00
eig: 0.005
bond charge predicted
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.178034835671E+03 -0.93771E-02 -0.36337E+00  3274 0.296E+00 0.937E-02
RMM: 2 -0.178038051320E+03 -0.32156E-02 -0.37021E-02  3262 0.293E-01 0.694E-02
RMM: 3 -0.178038089780E+03 -0.38460E-04 -0.12707E-03  3270 0.666E-02 0.629E-02
RMM: 4 -0.178038012742E+03 0.77038E-04 -0.10441E-04  1793 0.188E-02
198 F= -.17803801E+03 E0= -.17803453E+03  d E =-.125554E-01
BRION: g(F)=  0.175E+00 g(S)=  0.000E+00 retain N=  2 mean eig= 0.01
eig: 0.005  0.019
bond charge predicted
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.178038061183E+03 0.28597E-04 -0.51879E-02  3264 0.343E-01 0.137E-02
RMM: 2 -0.178038107135E+03 -0.45953E-04 -0.49736E-04  2332 0.336E-02 0.854E-03
RMM: 3 -0.178038108580E+03 -0.14447E-05 -0.34037E-05  1706 0.123E-02
199 F= -.17803811E+03 E0= -.17803467E+03  d E =-.958382E-04
BRION: g(F)=  0.746E-02 g(S)=  0.000E+00 retain N=  1 mean eig= 0.03
eig: 0.032
bond charge predicted
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.178043166425E+03 -0.50593E-02 -0.17555E+00  3274 0.208E+00 0.772E-02
RMM: 2 -0.178044778060E+03 -0.16116E-02 -0.17710E-02  3262 0.204E-01 0.560E-02
RMM: 3 -0.178044782047E+03 -0.39871E-05 -0.59850E-04  2607 0.463E-02
200 F= -.17804478E+03 E0= -.17804144E+03  d E =-.667347E-02
BRION: g(F)=  0.937E+00 g(S)=  0.000E+00 retain N=  2 mean eig= 0.02
eig: 0.032  0.004
start time  Mon Aug 13 14:54:12 CST 2012
end time  Tue Aug 14 01:51:52 CST 2012
。。。。。。。。。。。。。。。。
是不是說還沒有算完？算到200步就結(jié)束了，我的NSW設(shè)的是200.
還是就是我的是結(jié)構(gòu)優(yōu)化，IBRION設(shè)為1和2算的結(jié)果對我的結(jié)構(gòu)有影響沒？是不是為

贊一下

回復(fù)此樓

24樓2012-08-14 09:15:21

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

mousekingadv

金蟲 (小有名氣)

應(yīng)助: 51 (初中生)
金幣: 1141.8
紅花: 1
帖子: 140
在線: 28.7小時
蟲號: 1721002
注冊: 2012-03-27
專業(yè): 理論和計算化學(xué)

引用回帖:

24樓: Originally posted by rainxin2008 at 2012-08-14 09:15:21
你好，我把IBRION改為1算的結(jié)果如下： 195 F= -.17800925E+03 E0= -.17800569E+03  d E =-.497371E-02
BRION: g(F)=  0.255E-01 g(S)=  0.000E+00 retain N=  4 mean eig= 1.98
eig: 7.417  0.028  0.240  0 ...

確實沒算完。只是因為NSW到了就停了
IBRION和ISMEAR以及各種MIX，只要最后收斂了，收斂到的都是局域的能量最低點，不會因為方法不同而使基態(tài)變化。IBRION=2（CG方法）常用于大范圍地尋找能量低點，IBRION=1（quasi-Newton）在搜尋范圍較小但能得到精確的結(jié)果（CG基本EDIFFG=-0.01就是極限了，當(dāng)然0.01已經(jīng)很精確了）
為什么能量比你之前貼出來的下降了那么多（2eV）？五個原子吸附上去不可能有這么大的能量差吧。你的襯底原子間距是怎么得到的？
另外建議過一段時間就用VESTA看看CONTCAR，確認(rèn)在跑的東西是你想要的。

贊一下

回復(fù)此樓

25樓2012-08-14 10:30:14

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

rainxin2008

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1718.2
散金: 120
帖子: 265
在線: 163.3小時
蟲號: 1619644
注冊: 2012-02-15
專業(yè): 金屬功能材料

引用回帖:

25樓: Originally posted by mousekingadv at 2012-08-14 10:30:14
確實沒算完。只是因為NSW到了就停了
IBRION和ISMEAR以及各種MIX，只要最后收斂了，收斂到的都是局域的能量最低點，不會因為方法不同而使基態(tài)變化。IBRION=2（CG方法）常用于大范圍地尋找能量低點，IBRION=1（quas ...

我的襯底原子間距就是直接在MS里導(dǎo)入的Au,然后建的模型，間距也沒有手動修改，不知道為什么能量下降了那么多。怎么這么多問題啊，哎，太多不懂了

贊一下

回復(fù)此樓

26樓2012-08-14 10:37:15

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

mousekingadv

金蟲 (小有名氣)

應(yīng)助: 51 (初中生)
金幣: 1141.8
紅花: 1
帖子: 140
在線: 28.7小時
蟲號: 1721002
注冊: 2012-03-27
專業(yè): 理論和計算化學(xué)

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

★
liliangfang: 金幣+1, 謝謝交流 2012-08-14 18:59:27

引用回帖:

26樓: Originally posted by rainxin2008 at 2012-08-14 10:37:15
我的襯底原子間距就是直接在MS里導(dǎo)入的Au,然后建的模型，間距也沒有手動修改，不知道為什么能量下降了那么多。怎么這么多問題啊，哎，太多不懂了...

晶格常數(shù)不能使用實驗數(shù)據(jù)，必須用相同的方法和贗勢優(yōu)化一遍塊材Au。不需要真空層。ENCUT和K點不同也會有影響，不過差距不是很大。

贊一下

回復(fù)此樓

27樓2012-08-14 12:39:51

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

夢里花

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 140
紅花: 3
帖子: 26
在線: 41.7小時
蟲號: 1657927
注冊: 2012-03-02
專業(yè): 金屬材料的微觀結(jié)構(gòu)

改IBRION=1，就可以了，算到200步還不收斂，就先把POTIM改成0.5，收斂后再改成0.3左右繼續(xù)算，試試吧

贊一下

回復(fù)此樓

28樓2012-11-15 16:09:52

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

李加偉

榮譽版主 (文壇精英)

不過凡人

應(yīng)助: 258 (大學(xué)生)
貴賓: 1.257
金幣: 24680.9
散金: 8829
紅花: 146
沙發(fā): 13
帖子: 12823
在線: 3507.4小時
蟲號: 2359263
注冊: 2013-03-19
性別: GG
專業(yè): 可再生與替代能源利用中的
管轄: 論文道賀祈福

遇到相同問題，請問樓主后來怎么處理的

贊一下

回復(fù)此樓

人生海海，在于內(nèi)心

29樓2014-03-25 08:26:13

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

微6涼

鐵蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 202.8
帖子: 89
在線: 35.6小時
蟲號: 2724445
注冊: 2013-10-14
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶

引用回帖:

6樓: Originally posted by rainxin2008 at 2012-08-11 11:03:37
我的用CONTCAR替換POSCAR又算了一遍，還是不收斂。怎么辦啊�。考{了悶了...

請問樓主問題解決了嗎？我也遇到了同樣的問題

贊一下

回復(fù)此樓

30樓2015-08-11 16:47:37

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

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

返回列表

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

[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +3	beefly 2026-03-14	3/150	2026-03-17 16:45 by RRRKKK
[考研] 277調(diào)劑 +3	自由煎餅果子 2026-03-16	4/200	2026-03-17 16:36 by 無際的草原
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 289求調(diào)劑 +6	步川酷紫123 2026-03-11	6/300	2026-03-17 10:23 by Sammy2
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +3	張vvvv 2026-03-15	5/250	2026-03-16 20:25 by 張vvvv
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 279求調(diào)劑 +3	抓著星星的女孩 2026-03-10	3/150	2026-03-13 23:47 by userper
[考研] 0703化學(xué)調(diào)劑 +4	快樂的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 0856化學(xué)工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor
[考研] 296求調(diào)劑 +3	大口吃飯身體健 2026-03-13	3/150	2026-03-13 10:31 by 學(xué)員8dgXkO