| 5 | 1/1 | 返回列表 |
| 查看: 5568 | 回復(fù): 9 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||||
[交流]
【求助】PBE0 計(jì)算Si、Ge能帶結(jié)構(gòu) 已有6人參與
|
|||||
|
請教高手,用vasp5.2的pbe0計(jì)算硅鍺能帶結(jié)構(gòu)時(shí),所得的能帶結(jié)構(gòu)圖上總會(huì)出現(xiàn)尖銳的峰值,顯然計(jì)算數(shù)據(jù)有問題,用4.6版計(jì)算的能帶結(jié)構(gòu)和教科書上比較溫和,為何同樣的設(shè)置做PBE0時(shí)就出錯(cuò)了。請做過PBE0計(jì)算的大蝦指點(diǎn)一下 以下是我的個(gè)文件: 以ge為例,做優(yōu)化時(shí)的設(shè)置 INCAR: SYSTEM = Azoc MD PREC = Medium EDIFF = 1e-5 EDIFFG = -2e-2 IALGO = 48 NPAR = 2 LPLANE = TRUE LCHARG = .TRUE. NELMIN = 3 LREAL = F ROPT = 1E-3 ISYM = 1 NWRITE = 1 SPIN POLARIZATION SPECIFICATION ISPIN = 1 IONIC RELAXATION NSW = 200 NBLOCK = 1 IBRION = 2 POTIM = 0.5 LWAVE = .F. LCHG = .T. ISIF = 2 LCHARG = .TRUE. ELECTRON OCCUPATION ISMEAR = 0 occupation forced by INCAR DOS Stuff SIGMA = 0.1 check T*S to make sure LHFCALC = .TRUE. AEXX = 0.2 KPOINTS: 0 Monkhorst-Pack 4 4 4 0. 0. 0. POSCAR: Ge 5.6575 0.5 0.5 0.0 0.0 0.5 0.5 0.5 0.0 0.5 2 Cartesian 0.0000000000000 0.00000000000 0.0000000000000 0.2500000000000 0.25000000000 0.2500000000000 算能帶時(shí)設(shè)置 INCAR: SYSTEM = Azoc MD PREC = Medium EDIFF = 1e-5 EDIFFG = -2e-2 IALGO = 48 NPAR = 2 LPLANE = TRUE LCHARG = .TRUE. NELMIN = 3 LREAL = F ROPT = 1E-3 ISYM = 1 NWRITE = 1 ICHARG = 11 SPIN POLARIZATION SPECIFICATION ISPIN = 1 IONIC RELAXATION NSW = 1 NBLOCK = 1 IBRION = 1 POTIM = 0.5 LWAVE = .F. LCHG = .T. ISIF = 2 LCHARG = .TRUE. ELECTRON OCCUPATION ISMEAR = 0 occupation forced by INCAR DOS Stuff SIGMA = 0.1 check T*S to make sure LHFCALC = .TRUE. AEXX = 0.2 KPOINTS: Automatic mesh 19 rec 0.50000000 0.50000000 0.50000000 0.050 0.44444444 0.44444444 0.44444444 0.050 0.38888889 0.38888889 0.38888889 0.050 0.33333333 0.33333333 0.33333333 0.050 0.27777778 0.27777778 0.27777778 0.050 0.22222222 0.22222222 0.22222222 0.050 0.16666667 0.16666667 0.16666667 0.050 0.11111111 0.11111111 0.11111111 0.050 0.05555556 0.05555556 0.05555556 0.050 0.00000000 0.00000000 0.00000000 0.050 0.05555556 0.05555556 0.00000000 0.050 0.11111111 0.11111111 0.00000000 0.050 0.16666667 0.16666667 0.00000000 0.050 0.22222222 0.22222222 0.00000000 0.050 0.27777778 0.27777778 0.00000000 0.050 0.33333333 0.33333333 0.00000000 0.050 0.38888889 0.38888889 0.00000000 0.050 0.44444444 0.44444444 0.00000000 0.050 0.50000000 0.50000000 0.00000000 0.050 懇請高手指點(diǎn),非常急,拜托了! [ Last edited by 1984mjx on 2010-6-4 at 16:29 ] |
VASP | 計(jì)算與催化 |
金蟲 (正式寫手)
|
看到了樓上的帖子了,實(shí)際上官方論壇的VASP001 和現(xiàn)在的賬號(hào)是一個(gè)人,都是我提的。對于Si這個(gè)體系而言,你用HSE06 得出的帶隙和實(shí)驗(yàn)值僅差0.1ev,而且,計(jì)算時(shí)間比PBE0 少很多, 具體原因,自己看看,這個(gè)庫侖截?cái)噙是起的作用很大。此外,不一定用AlGO=Normal ,實(shí)際上ALGO=Damped 或者ALL 更快,我測試的結(jié)果是DAMPED最快。 |

至尊木蟲 (職業(yè)作家)
|
請參考vasp官方論壇提到的如何進(jìn)行hybrid functional來計(jì)算band structure的步驟: 下面的是從上面給拷貝下來的: band structures using hybrid functionals have to be calculated the following way: -) First perform a selfconsistent Hartree-Fock/HSE calculation using a conventional KPOINTS file. -) Copy the IBZKPT file to KPOINTS, and explicitely add all desired k-points along high-symmetry lines of the BZ that are needed for the bandstructure plot. Add the points at the end of the KPOINTS file, but set the weights of these added k-points to 0. -) Do not forget to set the number of k-points in KPOINTS correctly (to the number of the k-points used in the standard mesh PLUS the number of the k-points along the lines), such that all k-points are used for the calculations. -) NKRED can not be used. -) Perform a second VASP run: It is recommended to use the Davidson algorithm, since it converges that eigen energies at the new k-points fastest. Since VASP terminates when the total energy is converged to a certain threshold, it is important to force VASP to do a minimum number of steps, so that the orbitals at the new k-points are fully converged (note: since their weight is zero, they do not contribute to the total energy). This can be done using e.g. ALGO = N ; NELMIN = 5 ! Davidson, minimum 5 scf-steps IBRION = 1 ! Use simple charge mixer, since Pulay might blow up -) The KS-eigenvalues of the states along the high-symmetry lines are written in OUTCAR, EIGENVAL, vasprun.xml; please cut the k-points required for the bandstructure from one of these files and proceed as usual (using p4vasp or any other graphics tool you usually use to produce bandstructure-plots) Note: A Hartee-Fock calculation can NOT be continued from an existing CHGCAR file, since the non-local exchange is not determined by the charge density but by the density matrix and/or the KS-orbitals. |
| 首先非常感謝樓上的同仁的幫助,我也是自己在摸索pbe0的計(jì)算,不太懂里面的具體設(shè)置,花了不少時(shí)間也沒解決好這個(gè)問題,用PAE0做能帶計(jì)算的步驟和gga的不一樣嗎,我都是先做結(jié)構(gòu)優(yōu)化,然后將產(chǎn)生的CHGCAR/CHG文件考到新文件夾下,將ICHARG = 11設(shè)定后計(jì)算,最后在產(chǎn)生的文件EIGENVAL讀取能帶信息,請樓上的幫個(gè)忙,能把你做pBE0計(jì)算能帶的INCAR 、kpoints、poscar、貼出來共享一下嗎,不勝感激額,先行謝過了! |
金蟲 (正式寫手)
|
計(jì)算能帶的方法見VASP 官方論壇,有別于常規(guī)的DFT 能帶計(jì)算。。你有問題,再提問吧。 Band like this : > # Startparameter for this run: > ISPIN = 2 > ISTART = 0 > ICHARG = 2 > NWRITE = 2 > PREC = Normal > > # Electronic Relaxation > ENCUT = 400.0 > LREAL = AUTO > EDIFF = 1E-04 > NELM = 80 > NELMIN = 5 > > # Ionic relaxation > EDIFFG = -1E-02 > NSW = 0 > POTIM = 0.50 > IBRION = 1 > ISIF = 2 > ISYM = 2 > > > # DOS related values: > ISMEAR = 0; SIGMA = 0.20 > > > # Electronic relaxation 2 (details) > ALGO = Normal > > # Write flags > LWAVE =.TRUE. > LCHARG =.FALSE. > LORBIT = 11 > NBANDS=48 > # Accelorate computation > NSIM = 4 > NPAR = 64 > LPLANE =.TRUE. > LSCALU =.FALSE. > MAXMIX = 80 > > # Exchange correlation treatment: > LHFCALC = .TRUE. > HFSCREEN = 0.2 > AEXX = 0.25 > AGGAX = 0.75 > ALDAC = 1.0 > AGGAc = 1.0 > > > ENCUTFOCK= 0 > LMAXFOCK= 4 > # NKRED = 2 > TIME = 0.4 > > KPOINTS is like this : > > Automatically generated mesh > 56 > Reciprocal lattice > 0.05000000000000 0.05000000000000 0.12500000000000 8 0.15000000000000 > 0.05000000000000 0.12500000000000 16 0.25000000000000 0.05000000000000 > 0.12500000000000 16 0.35000000000000 0.05000000000000 0.12500000000000 > 16 0.45000000000000 0.05000000000000 0.12500000000000 16 > 0.15000000000000 0.15000000000000 0.12500000000000 8 0.25000000000000 > 0.15000000000000 0.12500000000000 16 0.35000000000000 0.15000000000000 > 0.12500000000000 16 0.45000000000000 0.15000000000000 0.12500000000000 > 16 0.25000000000000 0.25000000000000 0.12500000000000 8 > 0.35000000000000 0.25000000000000 0.12500000000000 16 0.45000000000000 > 0.25000000000000 0.12500000000000 16 0.35000000000000 0.35000000000000 > 0.12500000000000 8 0.45000000000000 0.35000000000000 0.12500000000000 > 16 0.45000000000000 0.45000000000000 0.12500000000000 8 > 0.05000000000000 0.05000000000000 0.37500000000000 8 0.15000000000000 > 0.05000000000000 0.37500000000000 16 0.25000000000000 0.05000000000000 > 0.37500000000000 16 0.35000000000000 0.05000000000000 0.37500000000000 > 16 0.45000000000000 0.05000000000000 0.37500000000000 16 > 0.15000000000000 0.15000000000000 0.37500000000000 8 0.25000000000000 > 0.15000000000000 0.37500000000000 16 0.35000000000000 0.15000000000000 > 0.37500000000000 16 0.45000000000000 0.15000000000000 0.37500000000000 > 16 0.25000000000000 0.25000000000000 0.37500000000000 8 > 0.35000000000000 0.25000000000000 0.37500000000000 16 0.45000000000000 > 0.25000000000000 0.37500000000000 16 0.35000000000000 0.35000000000000 > 0.37500000000000 8 0.45000000000000 0.35000000000000 0.37500000000000 > 16 0.45000000000000 0.45000000000000 0.37500000000000 8 0.0 0.0 0.0 0 > 0.1 0.1 0.0 0 > 0.2 0.2 0.0 0 > 0.3 0.3 0.0 0 > 0.4 0.4 0.0 0 > 0.5 0.5 0.0 0 > 0.5 0.5 0.1 0 > 0.5 0.5 0.2 0 > 0.5 0.5 0.3 0 > 0.5 0.5 0.4 0 > 0.5 0.5 0.5 0 > 0.5 0.4 0.5 0 > 0.5 0.3 0.5 0 > 0.5 0.2 0.5 0 > 0.5 0.1 0.5 0 > 0.5 0.0 0.5 0 > 0.4 0.0 0.5 0 > 0.3 0.0 0.5 0 > 0.2 0.0 0.5 0 > 0.1 0.0 0.5 0 > 0.0 0.0 0.5 0 > 0.0 0.0 0.4 0 > 0.0 0.0 0.3 0 > 0.0 0.0 0.2 0 > 0.0 0.0 0.1 0 > 0.0 0.0 0.0 0 這是HSE06 的,PBE0 ,要去掉那個(gè)HFSCREEN 截?cái)唷?br /> [ Last edited by vasp001 on 2010-6-9 at 14:28 ] |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +8 | 冬十三 2026-03-15 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級(jí)475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 材料專碩288分求調(diào)劑 一志愿211 +4 | 在家想你 2026-03-11 | 4/200 |
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 快樂的香蕉 2026-03-11 | 4/200 |
|
|
[考研] 一志愿中科院,化學(xué)方向,295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 301求調(diào)劑 +6 | Liyouyumairs 2026-03-11 | 6/300 |
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123! 2026-03-11 | 4/200 |
|
|
[考研] 270求調(diào)劑 085600材料與化工專碩 +3 | YXCT 2026-03-11 | 3/150 |
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|