| 5 | 1/1 | 返回列表 |
| 查看: 2537 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
求助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算出能帶的。懇請各位使用過GW方法的能夠給個(gè)建議,謝謝大家啦! |
鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[基金申請] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|