| 5 | 1/1 | 返回列表 |
| 查看: 2000 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
zw_tju木蟲 (小有名氣)
|
[求助]
【求助】GW或HSE雜化泛函帶結(jié)構(gòu)計(jì)算的問題
|
|
|
請教,我在使用HSE06或GW計(jì)算能帶時出現(xiàn)問題:error in IBZKPT_HF: two k-points are equivalent 1 60 this will cause problems in the HF routine。請問這是什么問題,如何解決啊。。。。 使用的計(jì)算過程是網(wǎng)上給出的計(jì)算Si的模版: GW Band Structure Plots To obtain band structure plots from the GW method, first you have to perform a DFT calculation. The INCAR file I used for DFT on Si looks like this: ISTART = 0 ISYM = 0 ICHARG = 2 ISMEAR = 0 SIGMA = 0.1 LORBIT = 11 NSW = 200 IBRION = 2 ISIF = 3 EDIFF = 0.0001 EDIFFG = -0.01 KPOINTS file for the run is: Monkhorst Pack 0 Gamma 4 4 4 0 0 0 POSCAR file is: Si 5.38936000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 0.5000000000000000 2 Direct 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.2500000000000000 0.2500000000000000 0.2500000000000000 POTCAR file is supplied with VASP 5.2, remember to use paw pseudopotentials so that in the band structure calculation the Wigner-Seitz radius for Si need not be directly specified in INCAR. After this run completes you are ready to perform the HF-type (HSE06) calculation. Change the INCAR file to: ISTART = 1 ICHARG = 2 EDIFF = 0.0001 EDIFFG = -0.01 ENCUT = 400 ENAUG = 800 LREAL = .FALSE. LWAVE = .TRUE. LCHARG = .TRUE. NELM = 200 NSW = 0 IBRION = -1 LMAXMIX = 4 ISMEAR = 0 SIGMA = 0.1 NSIM = 4 IALGO = 48 ISYM = 0 LHFCALC = .TRUE. HFSCREEN = 0.2 ALGO = Damped TIME = 0.4 ENCUTFOCK = 0 NKRED = 2 AEXX = 0.25 Keep the other files constant for this run (e.g. KPOINTS, POSCAR, POTCAR need not be changed). Once this run has completed, add to the bottom of the INCAR file from the HSE06 calculation: LOPTICS = .TRUE. This step calculates the dielectric matrix Once this run has completed, change ALGO = Damped to ALGO = GWo and ICHARG = 2 to ICHARG = 11 and add: NOMEGA = 50 ENCUTGW = 100 LORBIT = 11 and change the KPOINTS file to: k-points along high symmetry lines 20 ! 20 intersections Line-mode rec 0 0 0 ! gamma 0.5 0.0 0 ! X 0.5 0.0 0 ! X 0.5 0.5 0.0 ! M 0.5 0.5 0.0 ! M 0.0 0.0 0.0 ! gamma 0.0 0.0 0.0 ! gamma 0.5 0.5 0.0 ! M 0.5 0.5 0.0 ! M 0.5 0.5 0.5 ! R 0.5 0.5 0.5 ! R 0.0 0.0 0.0 ! gamma These values are the points of high symmetry along the lines of high symmetry within the first Brillouin Zone of the Si unit cell structure in reciprocal space. And that’s it! Once complete you can plot the PROCAR, EIGENVAL, and DOSCAR files from the output of this run as described in the previous posts. The band structure plot produced from the plotting utility is shown below. Quantum Quibbles » GW Band Structure Plots [ Last edited by zw_tju on 2011-5-1 at 09:02 ] |
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +7 | 妮妮ninicgb 2026-03-15 | 11/550 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
[考研] 085600材料與化工 309分請求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|