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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

junfeiyan

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 52.2
散金: 30
紅花: 1
帖子: 149
在線: 229.1小時
蟲號: 2419400
注冊: 2013-04-15
專業(yè): 固體力學(xué)

[求助] 求助大神， VASP 使用雜化泛函結(jié)構(gòu)優(yōu)化一直報錯已有5人參與

小弟最近使用雜化泛函進(jìn)行結(jié)構(gòu)優(yōu)化，原因是我算的是含f電子稀土化合物， vasp 給的贗勢用GGA+U 收斂極其困難，用雜化泛函好些，已經(jīng)有發(fā)表在nature 子篇的論文，做過類似的體系，用的雜化泛函優(yōu)化，求助各位大神，看看是小弟的輸入文件有問題，還是天河超算中心的VASP 編譯問題，感謝大神
小弟的INCAR

Name of system
SYSTEM =    a

Start parameter
ISTART =    0
INIWAV =    1
ICHARG =    2

Write flags
NWRITE =    2
LWAVE    =    T
LCHARG =    T

Parallelization
NPAR    =    72
LPLANE =    T
LSCALU =    F
LSCALAPACK = F

Electronic Relaxation 1
ENCUT    =    500.00  eV
EDIFF    =    1.0E-05  eV
PREC    =    High
NELM    =    200
NELMIN =    2

Electronic Relaxation 2
ALGO    =    Damped
TIME    =    0.4
LREAL    =    Auto

Ionic Relaxation
NSW    =    150
IBRION =    2
ISIF    =    3
EDIFFG =    -0.05
POTIM    =    0.3
ISYM    =    3
IWAVPR =    1
LCORR    =    T

DOS related values
ISMEAR =    0
SIGMA    =    0.025
NEDOS    =    3000
# NBANDS =    3958

Magnetic parameter
LORBIT =    11
Exchange-correlation functional
GGA    =    PS
ISPIN = 2
  MAGMOM = 4*8 4*4 12*0
HSE06 parameter
LHFCALC =    T
HFSCREEN  =    0.2
PRECFOCK  =    Fast
# NKRED    =    2

還有一個問體就是我想用HSE 考慮稀土元素的f 電子，那么我選擇贗勢的時候是否可以用PBE 里面不含f 電子的贗勢，不知道雜化泛函會不會自己考慮f 電子
前提是我用GGA+U 算不含f電子的贗勢，很快收斂，用含f電子的贗勢極難收斂
報錯一
ROT:  1 0.159495398721E+04 0.17831E+04BRMIX: very serious problems
the old and the new charge density differ
old charge density: 160.00000 new***********
0.603E+06
ERROR FEXCF: supplied exchange-correlation table
is too small, maximal index :    912695

這個報錯有人說應(yīng)該減小POTIM。有成功的經(jīng)驗嗎
報錯2
ROT:  1  -0.155058592914E+03 -0.30390E+02 0.284E+01
ROT:  2  -0.141575035880E+03 0.13484E+02 0.160E+01
ROT:  3  -0.139046123851E+03 0.25289E+01 0.114E+01
ROT:  4  -0.139495088906E+03 -0.44897E+00 0.453E+00
ROT:  5  -0.139199449050E+03 0.29564E+00 0.133E+00
gam= 0.382 g(H,U,f)=  0.305E+03 0.218E+02 0.320E+00 ort(H,U,f) =-0.188E+03-0.426E+02-0.465E+00
DMP: 7 -0.112154452901E+03 -0.50536E+02 -0.95574E+02  2160 0.327E+03-0.231E+03
      N    E                   dE          rms(c)
ROT:  1  -0.159501125680E+03 -0.20006E+02 0.795E+01
ROT:  2  -0.157770040931E+03 0.17311E+01 0.495E+01
ROT:  3  -0.155549182836E+03 0.22209E+01 0.173E+01
ROT:  4  NaN                   NaN          NaN
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image             PC             Routine          Line       Source
vasp             00000000005A769A  Unknown             Unknown  Unknown
vasp             000000000070A100  Unknown             Unknown  Unknown

U31~GV(QLRG}V{D75RS6[]8.png

回復(fù)此樓

1樓 2015-08-28 09:52:29

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

回帖支持 ( 顯示支持度最高的前 50 名 )

golddoushi

木蟲 (正式寫手)

應(yīng)助: 36 (小學(xué)生)
金幣: 4007.2
散金: 25
紅花: 15
帖子: 337
在線: 447.7小時
蟲號: 2546370
注冊: 2013-07-15
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

ALGO=Normal
就行了

贊一下(1人)

回復(fù)此樓

2樓2015-08-28 11:36:41

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

lshok

鐵桿木蟲 (著名寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 5536.1
紅花: 4
沙發(fā): 1
帖子: 1537
在線: 262.5小時
蟲號: 1225253
注冊: 2011-03-08
專業(yè): 物理學(xué)II

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

★
感謝參與，應(yīng)助指數(shù) +1
liliangfang: 金幣+1, 謝謝交流 2015-09-05 18:29:29

NPAR=72?
VASP 推薦NPAR＝sqrt（cores），你設(shè)置的NPAR＝72對應(yīng)5184個核。。。
瞅瞅是不是你這個參數(shù)設(shè)置過大，導(dǎo)致內(nèi)存不夠引起的錯誤

贊一下(1人)

回復(fù)此樓

5樓2015-08-29 16:44:55

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

emilyoyang

木蟲 (正式寫手)

應(yīng)助: 232 (大學(xué)生)
金幣: 4340.4
紅花: 42
帖子: 667
在線: 189.9小時
蟲號: 1630291
注冊: 2012-02-20
性別: MM
專業(yè): 理論和計算化學(xué)

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

★
liliangfang: 金幣+1, 謝謝交流 2015-09-05 18:29:36

引用回帖:

3樓: Originally posted by junfeiyan at 2015-08-28 15:15:50
謝謝您，麻煩在問一下，你說的這個參數(shù)，是對應(yīng)哪個報錯的呢，還有幫我看看其他參數(shù)有沒有需要修改的呢，小弟沒怎么用過雜化泛函...

上面的人肯定沒算過雜化泛函。ALGO=Normal不用于雜化泛函計計算，對于雜化泛函需要ALGO=ALL or Damped
仔細(xì)看看vasp manual http://cms.mpi.univie.ac.at/vasp ... k_calculations.html

贊一下(1人)

回復(fù)此樓

7樓2015-09-01 13:53:05

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

普通回帖

junfeiyan

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 52.2
散金: 30
紅花: 1
帖子: 149
在線: 229.1小時
蟲號: 2419400
注冊: 2013-04-15
專業(yè): 固體力學(xué)

引用回帖:

2樓: Originally posted by golddoushi at 2015-08-28 11:36:41
ALGO=Normal
就行了

謝謝您，麻煩在問一下，你說的這個參數(shù)，是對應(yīng)哪個報錯的呢，還有幫我看看其他參數(shù)有沒有需要修改的呢，小弟沒怎么用過雜化泛函

贊一下

回復(fù)此樓

3樓2015-08-28 15:15:50

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

zuocuiping

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

應(yīng)助: 14 (小學(xué)生)
金幣: 4914.9
散金: 966
紅花: 60
帖子: 4023
在線: 1597.1小時
蟲號: 890422
注冊: 2009-11-01
性別: GG
專業(yè): 半導(dǎo)體光電子器件

能把你說的參考文獻(xiàn)的題目發(fā)給我嗎
我想看一下
謝謝

贊一下

回復(fù)此樓

4樓2015-08-29 12:50:43

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

junfeiyan

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 52.2
散金: 30
紅花: 1
帖子: 149
在線: 229.1小時
蟲號: 2419400
注冊: 2013-04-15
專業(yè): 固體力學(xué)

引用回帖:

5樓: Originally posted by lshok at 2015-08-29 16:44:55
NPAR=72?
VASP 推薦NPAR＝sqrt（cores），你設(shè)置的NPAR＝72對應(yīng)5184個核。。。
瞅瞅是不是你這個參數(shù)設(shè)置過大，導(dǎo)致內(nèi)存不夠引起的錯誤

謝謝，回復(fù)，我也了解前輩說的這個參數(shù)，很多人推薦設(shè)置您說的核數(shù)的平方根，也有人說，直接設(shè)置為核數(shù)也行，我再測試測試

贊一下

回復(fù)此樓

6樓2015-08-29 17:53:56

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

emilyoyang

木蟲 (正式寫手)

應(yīng)助: 232 (大學(xué)生)
金幣: 4340.4
紅花: 42
帖子: 667
在線: 189.9小時
蟲號: 1630291
注冊: 2012-02-20
性別: MM
專業(yè): 理論和計算化學(xué)

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

★
liliangfang: 金幣+1, 謝謝交流 2015-09-05 18:29:43

引用回帖:

6樓: Originally posted by junfeiyan at 2015-08-29 17:53:56
謝謝，回復(fù)，我也了解前輩說的這個參數(shù)，很多人推薦設(shè)置您說的核數(shù)的平方根，也有人說，直接設(shè)置為核數(shù)也行，我再測試測試...

這回復(fù)。。。。瞎說！雜化泛函要要求NPAR=#of cores
你的NPAR設(shè)置沒問題

贊一下

回復(fù)此樓

8樓2015-09-01 13:54:07

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

emilyoyang

木蟲 (正式寫手)

應(yīng)助: 232 (大學(xué)生)
金幣: 4340.4
紅花: 42
帖子: 667
在線: 189.9小時
蟲號: 1630291
注冊: 2012-02-20
性別: MM
專業(yè): 理論和計算化學(xué)

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

★
liliangfang: 金幣+1, 謝謝交流 2015-09-05 18:29:51

你的問題我來回答一下吧

1. ERROR FEXCF: supplied exchange-correlation table
is too small, maximal index : 912695
這個問題是由于intel編譯器的問題，（我猜）天河上用intel可能是2013 or 2015。你用intel2011編譯的vasp運行是不會出現(xiàn)這個錯誤的。
BTW，這個問題偶爾通過減小SIGMA 還是POTIM可以得到解決(以前試過，忘了，anyway, 這不是終極解決方法）

2. ROT: 4 NaN NaN NaN
forrtl: severe (174): SIGSEGV, segmentation fault occurred
這個錯誤明顯是你體系energy 不收斂。你看energy 都NaN了！所以才死掉。
4f電子計算有很多tricks。你慢慢研究研究吧。

贊一下

回復(fù)此樓

9樓2015-09-01 13:59:55

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

junfeiyan

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 52.2
散金: 30
紅花: 1
帖子: 149
在線: 229.1小時
蟲號: 2419400
注冊: 2013-04-15
專業(yè): 固體力學(xué)

引用回帖:

9樓: Originally posted by emilyoyang at 2015-09-01 13:59:55
你的問題我來回答一下吧

1. ERROR FEXCF: supplied exchange-correlation table
is too small, maximal index : 912695
這個問題是由于intel編譯器的問題，（我猜）天河上用intel可能是2013 or 2015。 ...

謝謝，謝謝，希望前輩介紹一下關(guān)于4f電子計算的經(jīng)驗，讓晚生少有些彎路！

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

10樓2015-09-01 21:09:09

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

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

返回列表

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

[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8	皮皮皮123456 2026-03-22	8/400	2026-03-23 23:47 by 熱情沙漠
[考研] 299求調(diào)劑 +5	某某某某位 2026-03-21	5/250	2026-03-23 23:36 by 熱情沙漠
[考研] 材料專業(yè)求調(diào)劑 +11	hanamiko 2026-03-18	11/550	2026-03-23 23:12 by peike
[考研] 327求調(diào)劑 +5	prayer13 2026-03-23	5/250	2026-03-23 22:11 by 星空星月
[考研] 化學(xué)308分求調(diào)劑 +3	你好明天你好 2026-03-23	3/150	2026-03-23 20:11 by macy2011
[考研] 0703化學(xué)求調(diào)劑 +4	奶油草莓. 2026-03-22	5/250	2026-03-23 19:37 by pswait
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +6	zxts12 2026-03-21	9/450	2026-03-23 18:37 by macy2011
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 08工學(xué)調(diào)劑 +7	用戶573181 2026-03-20	11/550	2026-03-23 15:47 by 我愛學(xué)習(xí)學(xué)習(xí)使?/a>
[考研] 化學(xué)調(diào)劑 +5	yzysaa 2026-03-21	5/250	2026-03-21 22:12 by peike
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO
[考研] 301求調(diào)劑 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-21 00:39 by JourneyLucky
[考研] 274求調(diào)劑 +10	S.H1 2026-03-18	10/500	2026-03-20 23:51 by JourneyLucky
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 329求調(diào)劑 +9	想上學(xué)吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨

24小時熱門版塊排行榜

[求助] 求助大神， VASP 使用雜化泛函結(jié)構(gòu)優(yōu)化一直報錯 已有5人參與

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

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

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

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

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

[求助] 求助大神， VASP 使用雜化泛函結(jié)構(gòu)優(yōu)化一直報錯已有5人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: