版塊導(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í)熱門版塊排行榜

返回列表

onion2440

銀蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 95.5
散金: 650
紅花: 4
帖子: 838
在線: 263.8小時(shí)
蟲號(hào): 1790016
注冊(cè): 2012-05-02
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[求助] 求助vasp GW方法計(jì)算過渡金屬錯(cuò)誤CHECK_OCCUPATIONS 已有1人參與

流程參考網(wǎng)上vasp gw計(jì)算方法帖子。
原帖地址：http://blog.sina.com.cn/s/blog_7c63531f0101ibx9.html
第一步的INCAR：
Start parameter for this run:
  ISTART=0
  ICHARG=2
#  ISIF =3
  NSW=0
  IBRION=-1
#  LWAVE=F
#  LCHARG=F
ENCUT = 400.0
EDIFF = 1E-4
#EDIFFG=-1E-2
ISMEAR=0
LREAL=A
SIGMA=0.1
#PREC=A
POTIM=0.35
#NPAR=4
#NELMIN=4
NBANDS=64
第二步的INCAR：
Start parameter for this run:
  ISTART=1
  ICHARG=2
#  ISIF =3
  NSW=0
  IBRION=-1
#  LWAVE=F
#  LCHARG=F
ENCUT = 400.0
EDIFF = 1E-4
#EDIFFG=-1E-2
ISMEAR=0
LREAL=A
SIGMA=0.1
#PREC=A
POTIM=0.35
#NPAR=4
#NELMIN=4
NBANDS=64
ALGO = Exact
LOPTICS=.TRUE.
（這兩步vasp沒有提示錯(cuò)誤）
第三步INCAR：
Start parameter for this run:
  ISTART=1
  ICHARG=1
#  ISIF =3
  NSW=0
  IBRION=-1
#  LWAVE=F
#  LCHARG=F
ENCUT = 400.0
EDIFF = 1E-4
#EDIFFG=-1E-2
ISMEAR=0
LREAL=A
SIGMA=0.1
#PREC=A
POTIM=0.35
#NPAR=4
#NELMIN=4
NBANDS=64
ALGO = GW0
LOPTICS=.TRUE.
LSPECTRAL = T
NOMEGA = 64
NBANDSGW=64
ENCUTGW=100
PRECFOCK=F
LWANNIER90=.TRUE.
LRPA=.FALSE.
在以上三步vasp計(jì)算過程中保持KPOINTS等其他的文件不變。使用的計(jì)算腳本如下：
#!/bin/bash
mpirun -np 16 vaspwannier9012 >out
wait
mkdir optic
cp INCAR* KPOINTS POTCAR POSCAR WAVE* CHGCAR ./optic
cd optic
cp INCAR-optic INCAR
mpirun -np 16  vaspwannier9012 >out
wait
mkdir gw
cp INCAR-gw KPOINTS POTCAR POSCAR WAVE* CHGCAR ./gw
cd gw
cp INCAR-gw INCAR
mpirun -np 16 vaspwannier9012 >out

問題出在第三步：
out文件錯(cuò)誤提示如下：
NQ= 168 -0.0909 0.4545 0.4000,
|.........|.........|.........|.........|.........|.........|.........|.........
performed    923544000 updates of chi_q(r,r)
total number of BLAS operations 102922.89 Gflops
resolving degeneracies of screened two electron integrals

calculate QP shifts <psi_nk| G(iteration)W_0 |psi_nk>: iteration 1
writing wavefunctions
entering main loop
optical routines
imaginary and real dielectric function
recalculating local potential from charge density
CHECK_OCCUPATIONS: ERROR, system is not insulating.
CHECK_OCCUPATIONS: ERROR, system is not insulating.
CHECK_OCCUPATIONS: ERROR, system is not insulating.
根據(jù)提示，我不斷增加NBANDS空帶的條數(shù)，然而始終提示這個(gè)錯(cuò)誤。
此外，如果把KPOINTS改為如下：
KPOINTS
0
G
1 1 1
其他參數(shù)保持不變，結(jié)果就沒問題。但是這時(shí)候的結(jié)果應(yīng)該是沒有意義的，應(yīng)該是無法用wannier90算出能帶的。懇請(qǐng)各位使用過GW方法的能夠給個(gè)建議，謝謝大家啦！

回復(fù)此樓

» 猜你喜歡

321求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中海洋材料工程專碩330分求調(diào)劑已經(jīng)有7人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有5人回復(fù)
288求調(diào)劑已經(jīng)有3人回復(fù)
材料080500調(diào)劑求收留已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有9人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有8人回復(fù)
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑已經(jīng)有10人回復(fù)
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑已經(jīng)有3人回復(fù)
332求調(diào)劑已經(jīng)有3人回復(fù)

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

pwscf收斂速度，參數(shù)調(diào)節(jié)與贗勢(shì) 已經(jīng)有19人回復(fù)
CESTEP計(jì)算第五周期過渡金屬氧化物參數(shù)設(shè)置的一些疑問已經(jīng)有11人回復(fù)
求助quantum-espresso使用問題！已經(jīng)有5人回復(fù)
新手求助一個(gè)quantum-espresso的使用問題！已經(jīng)有13人回復(fù)
vasp接口wannier90計(jì)算問題已經(jīng)有24人回復(fù)
手冊(cè)中計(jì)算能帶時(shí)k點(diǎn)設(shè)置的問題已經(jīng)有12人回復(fù)
新手求助pwscf的進(jìn)行磁性運(yùn)算和LDA+U運(yùn)算已經(jīng)有35人回復(fù)
Problem - nscf - c_bands: eigenvalues not converged 已經(jīng)有20人回復(fù)
vasp中利用GW算法算能帶結(jié)構(gòu) 已經(jīng)有9人回復(fù)
ABINIT和VASP相比，計(jì)算GW，哪個(gè)更快啊？已經(jīng)有13人回復(fù)
vasp跨節(jié)點(diǎn)運(yùn)行出錯(cuò)，mpiexec_node-1 (handle_stdin_input 1089) 已經(jīng)有5人回復(fù)
ms中吸附能的計(jì)算已經(jīng)有40人回復(fù)
第一原理軟件的兩種磁性計(jì)算方法已經(jīng)有14人回復(fù)
【求助】結(jié)構(gòu)優(yōu)化不收斂已經(jīng)有9人回復(fù)
【求助】VASP計(jì)算光學(xué)性質(zhì)時(shí)的GW方法（RPA, no local field effects) 已經(jīng)有5人回復(fù)
【求助成功】Vasp 5.2.11編譯出錯(cuò) 已經(jīng)有10人回復(fù)
【求助】對(duì)于難收斂的體系已經(jīng)有9人回復(fù)
【求助】pwscf結(jié)構(gòu)優(yōu)化中結(jié)果出的怪異啊已經(jīng)有16人回復(fù)
【求助】VASP計(jì)算電子態(tài)密度已經(jīng)有34人回復(fù)
【求助】vasp5.2的GW怎么算？已經(jīng)有14人回復(fù)

1樓 2015-07-11 21:27:20

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

KalaShayminS

鐵桿木蟲 (著名寫手)

應(yīng)助: 432 (碩士)
金幣: 6174.2
紅花: 56
帖子: 1088
在線: 379小時(shí)
蟲號(hào): 2408111
注冊(cè): 2013-04-08

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

★
感謝參與，應(yīng)助指數(shù) +1
liliangfang: 金幣+1, 謝謝交流 2015-11-22 11:17:23

VASP的GW不能計(jì)算金屬，有的能算結(jié)果無意義，有的直接禁止計(jì)算。

贊一下(1人)

回復(fù)此樓

2樓2015-07-12 19:07:54

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

onion2440

銀蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 95.5
散金: 650
紅花: 4
帖子: 838
在線: 263.8小時(shí)
蟲號(hào): 1790016
注冊(cè): 2012-05-02
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

2樓: Originally posted by KalaShayminS at 2015-07-12 19:07:54
VASP的GW不能計(jì)算金屬，有的能算結(jié)果無意義，有的直接禁止計(jì)算。

你好，謝謝回復(fù)！你的意思是，如果能帶是金屬型，那么GW是不能計(jì)算的，如果能帶是半金屬性或者半導(dǎo)體性，這時(shí)候就可以用GW計(jì)算，是嗎？

贊一下

回復(fù)此樓

3樓2015-07-12 22:33:06

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

KalaShayminS

鐵桿木蟲 (著名寫手)

應(yīng)助: 432 (碩士)
金幣: 6174.2
紅花: 56
帖子: 1088
在線: 379小時(shí)
蟲號(hào): 2408111
注冊(cè): 2013-04-08

引用回帖:

3樓: Originally posted by onion2440 at 2015-07-12 22:33:06
你好，謝謝回復(fù)！你的意思是，如果能帶是金屬型，那么GW是不能計(jì)算的，如果能帶是半金屬性或者半導(dǎo)體性，這時(shí)候就可以用GW計(jì)算，是嗎？...

vasp里是這樣。

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

回復(fù)此樓

4樓2015-07-13 10:03:28

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

onion2440

銀蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 95.5
散金: 650
紅花: 4
帖子: 838
在線: 263.8小時(shí)
蟲號(hào): 1790016
注冊(cè): 2012-05-02
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

4樓: Originally posted by KalaShayminS at 2015-07-13 10:03:28
vasp里是這樣。
...

如果換一種方法，例如基于HSE的波函數(shù)做GW計(jì)算，或者直接采用GW波函數(shù)自洽計(jì)算，是不是會(huì)好一些？
另外，abinit等其他的第一性原理code做GW計(jì)算，也會(huì)遇到同樣的問題嗎？

贊一下

回復(fù)此樓

5樓2015-07-13 10:20:32

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

KalaShayminS

鐵桿木蟲 (著名寫手)

應(yīng)助: 432 (碩士)
金幣: 6174.2
紅花: 56
帖子: 1088
在線: 379小時(shí)
蟲號(hào): 2408111
注冊(cè): 2013-04-08

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
onion2440: 金幣+50, ★★★★★最佳答案, 回復(fù)很詳細(xì)，很有幫助，非常感謝！ 2015-07-13 10:39:16
liliangfang: 金幣+1, 謝謝交流 2015-11-22 11:17:29

引用回帖:

5樓: Originally posted by onion2440 at 2015-07-13 10:20:32
如果換一種方法，例如基于HSE的波函數(shù)做GW計(jì)算，或者直接采用GW波函數(shù)自洽計(jì)算，是不是會(huì)好一些？
另外，abinit等其他的第一性原理code做GW計(jì)算，也會(huì)遇到同樣的問題嗎？...

這取決于計(jì)算的初始態(tài)，如果HSE自洽結(jié)果不是金屬，就可以算，否則和LDA/GGA并無區(qū)別。GW自洽也需要初始態(tài)，在這個(gè)問題上和G0W0沒有區(qū)別。
ABINIT可以算，但結(jié)果也需要謹(jǐn)慎。

贊一下

回復(fù)此樓

6樓2015-07-13 10:32:02

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

gaozhibin

銅蟲 (正式寫手)

應(yīng)助: 9 (幼兒園)
金幣: 100.9
散金: 202
紅花: 7
帖子: 391
在線: 149.5小時(shí)
蟲號(hào): 2617861
注冊(cè): 2013-08-26
性別: MM
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

6樓: Originally posted by KalaShayminS at 2015-07-13 10:32:02
這取決于計(jì)算的初始態(tài)，如果HSE自洽結(jié)果不是金屬，就可以算，否則和LDA/GGA并無區(qū)別。GW自洽也需要初始態(tài)，在這個(gè)問題上和G0W0沒有區(qū)別。
ABINIT可以算，但結(jié)果也需要謹(jǐn)慎。...

大神你好，我用GW 算到第三步的時(shí)候，一直卡在這是是什么意思？

responsefunction array rank=6432
LDA part: xc-table for Pade appr. of Perdew
allocating 2 responsefunctions rank=  6432
resonant part calculated as well
allocating 2 responsefunctions rank=  6432
allocating 1 responsefunctions rank=  6432
Doing          2  frequencies on each core in blocks of          1
NQ= 1 0.0000 0.0000 0.0000,

下面一直不動(dòng)了。。。請(qǐng)大神指導(dǎo)。

贊一下

回復(fù)此樓

加油努力，堅(jiān)持爬行的小烏龜

7樓2015-11-21 10:04:12

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

gy_ustc

捐助貴賓 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 162.6
帖子: 70
在線: 32.4小時(shí)
蟲號(hào): 3120049
注冊(cè): 2014-04-08
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

7樓: Originally posted by gaozhibin at 2015-11-21 10:04:12
大神你好，我用GW 算到第三步的時(shí)候，一直卡在這是是什么意思？

responsefunction array rank=6432
LDA part: xc-table for Pade appr. of Perdew
allocating 2 responsefunctions rank= 6432
resona ...

前輩您好！您已經(jīng)解決了嗎？我現(xiàn)在遇到一模一樣的問題，我也卡在這個(gè)地方了

贊一下

回復(fù)此樓

8樓2018-08-11 10:25:58

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

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

返回列表

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

[考研] 271材料工程求調(diào)劑 +5	.6lL 2026-03-18	5/250	2026-03-19 03:07 by 無懈可擊111
[考研] 085700資源與環(huán)境308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-18 22:35 by bingxueer79
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4	守候夕陽CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	9/450	2026-03-18 15:31 by vgtyfty
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 考研調(diào)劑 +4	芬達(dá)46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7	dtdxzxx 2026-03-12	8/400	2026-03-13 14:43 by jxchenghu
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike