| 5 | 1/1 | 返回列表 |
| 查看: 2539 | 回復: 7 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
求助vasp GW方法計算過渡金屬錯誤CHECK_OCCUPATIONS 已有1人參與
|
||
|
流程參考網(wǎng)上vasp gw計算方法帖子。 原帖地址: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沒有提示錯誤) 第三步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計算過程中保持KPOINTS等其他的文件不變。使用的計算腳本如下: #!/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文件錯誤提示如下: 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ù),然而始終提示這個錯誤。 此外,如果把KPOINTS改為如下: KPOINTS 0 G 1 1 1 其他參數(shù)保持不變,結果就沒問題。但是這時候的結果應該是沒有意義的,應該是無法用wannier90算出能帶的。懇請各位使用過GW方法的能夠給個建議,謝謝大家啦! |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 0703化學調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 東南大學364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 080500,材料學碩302分求調(diào)劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|