 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2540 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
求助vasp GW方法計(jì)算過(guò)渡金屬錯(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沒(méi)有提示錯(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ì)算過(guò)程中保持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 問(wèn)題出在第三步： 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é)果就沒(méi)問(wèn)題。但是這時(shí)候的結(jié)果應(yīng)該是沒(méi)有意義的，應(yīng)該是無(wú)法用wannier90算出能帶的。懇請(qǐng)各位使用過(guò)GW方法的能夠給個(gè)建議，謝謝大家啦！ |
回復(fù)此樓» 猜你喜歡
285化工學(xué)碩求調(diào)劑（081700）
已經(jīng)有10人回復(fù)
-
求材料調(diào)劑
已經(jīng)有10人回復(fù)
-
本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑
已經(jīng)有4人回復(fù)
-
國(guó)自科面上基金字體
已經(jīng)有8人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有11人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有9人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
0703化學(xué) 305求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求各位老師收留
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- pwscf收斂速度，參數(shù)調(diào)節(jié)與贗勢(shì)
已經(jīng)有19人回復(fù)
- CESTEP計(jì)算第五周期過(guò)渡金屬氧化物參數(shù)設(shè)置的一些疑問(wèn)
已經(jīng)有11人回復(fù)
- 求助quantum-espresso使用問(wèn)題！
已經(jīng)有5人回復(fù)
- 新手求助一個(gè)quantum-espresso的使用問(wèn)題！
已經(jīng)有13人回復(fù)
- vasp接口wannier90計(jì)算問(wèn)題
已經(jīng)有24人回復(fù)
- 手冊(cè)中計(jì)算能帶時(shí)k點(diǎn)設(shè)置的問(wèn)題
已經(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ù)
|
大神你好，我用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)。 |
7樓2015-11-21 10:04:12
KalaShayminS
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 432 (碩士)
- 金幣: 6174.2
- 紅花: 56
- 帖子: 1088
- 在線: 379小時(shí)
- 蟲號(hào): 2408111
- 注冊(cè): 2013-04-08
2樓2015-07-12 19:07:54
3樓2015-07-12 22:33:06
KalaShayminS
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 432 (碩士)
- 金幣: 6174.2
- 紅花: 56
- 帖子: 1088
- 在線: 379小時(shí)
- 蟲號(hào): 2408111
- 注冊(cè): 2013-04-08
4樓2015-07-13 10:03:28
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
