| 8 | 1/1 | 返回列表 |
| 查看: 2542 | 回復(fù): 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ù)保持不變,結(jié)果就沒問題。但是這時候的結(jié)果應(yīng)該是沒有意義的,應(yīng)該是無法用wannier90算出能帶的。懇請各位使用過GW方法的能夠給個建議,謝謝大家啦! |
鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
|
這取決于計算的初始態(tài),如果HSE自洽結(jié)果不是金屬,就可以算,否則和LDA/GGA并無區(qū)別。GW自洽也需要初始態(tài),在這個問題上和G0W0沒有區(qū)別。 ABINIT可以算,但結(jié)果也需要謹(jǐn)慎。 |
|
大神你好,我用GW 算到第三步的時候,一直卡在這是是什么意思? 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ǎo)。 |

捐助貴賓 (小有名氣)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工(085600)總分338 +4 | 蔡大美女 2026-03-13 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +5 | 小材化本科 2026-03-18 | 5/250 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
[考研] 268求調(diào)劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 土木第一志愿276求調(diào)劑,科研和技能十分豐富,求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|