| 24小時(shí)熱門(mén)版塊排行榜 |
| 10 | 1/1 | 返回列表 |
| 查看: 5566 | 回復(fù): 9 | |||||
[交流]
【求助】PBE0 計(jì)算Si、Ge能帶結(jié)構(gòu) 已有6人參與
|
|
請(qǐng)教高手,用vasp5.2的pbe0計(jì)算硅鍺能帶結(jié)構(gòu)時(shí),所得的能帶結(jié)構(gòu)圖上總會(huì)出現(xiàn)尖銳的峰值,顯然計(jì)算數(shù)據(jù)有問(wèn)題,用4.6版計(jì)算的能帶結(jié)構(gòu)和教科書(shū)上比較溫和,為何同樣的設(shè)置做PBE0時(shí)就出錯(cuò)了。請(qǐng)做過(guò)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 懇請(qǐng)高手指點(diǎn),非常急,拜托了! [ Last edited by 1984mjx on 2010-6-4 at 16:29 ] |
VASP | 計(jì)算與催化 |
至尊木蟲(chóng) (職業(yè)作家)
|
請(qǐng)參考vasp官方論壇提到的如何進(jìn)行hybrid functional來(lái)計(jì)算band structure的步驟: 下面的是從上面給拷貝下來(lái)的: 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. |
金蟲(chóng) (正式寫(xiě)手)
|
看到了樓上的帖子了,實(shí)際上官方論壇的VASP001 和現(xiàn)在的賬號(hào)是一個(gè)人,都是我提的。對(duì)于Si這個(gè)體系而言,你用HSE06 得出的帶隙和實(shí)驗(yàn)值僅差0.1ev,而且,計(jì)算時(shí)間比PBE0 少很多, 具體原因,自己看看,這個(gè)庫(kù)侖截?cái)噙是起的作用很大。此外,不一定用AlGO=Normal ,實(shí)際上ALGO=Damped 或者ALL 更快,我測(cè)試的結(jié)果是DAMPED最快。 |

| 首先非常感謝樓上的同仁的幫助,我也是自己在摸索pbe0的計(jì)算,不太懂里面的具體設(shè)置,花了不少時(shí)間也沒(méi)解決好這個(gè)問(wèn)題,用PAE0做能帶計(jì)算的步驟和gga的不一樣嗎,我都是先做結(jié)構(gòu)優(yōu)化,然后將產(chǎn)生的CHGCAR/CHG文件考到新文件夾下,將ICHARG = 11設(shè)定后計(jì)算,最后在產(chǎn)生的文件EIGENVAL讀取能帶信息,請(qǐng)樓上的幫個(gè)忙,能把你做pBE0計(jì)算能帶的INCAR 、kpoints、poscar、貼出來(lái)共享一下嗎,不勝感激額,先行謝過(guò)了! |
金蟲(chóng) (正式寫(xiě)手)
|
計(jì)算能帶的方法見(jiàn)VASP 官方論壇,有別于常規(guī)的DFT 能帶計(jì)算。。你有問(wèn)題,再提問(wèn)吧。 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 ] |

金蟲(chóng) (小有名氣)

金蟲(chóng) (小有名氣)

木蟲(chóng) (正式寫(xiě)手)

|
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +4 | yy要上岸呀 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +8 | 小熊joy 2026-03-14 | 8/400 |
|
|
[考研] 材料工程專(zhuān)碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +11 | 大米飯! 2026-03-15 | 14/700 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 277調(diào)劑 +4 | 自由煎餅果子 2026-03-16 | 5/250 |
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
[考研] 085601專(zhuān)碩,總分342求調(diào)劑,地區(qū)不限 +3 | share_joy 2026-03-16 | 3/150 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 332求調(diào)劑 +3 | zjy101327 2026-03-11 | 6/300 |
|
|
[考研] 308求調(diào)劑 +5 | 是Lupa啊 2026-03-11 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考博] 2026年博士申請(qǐng) +3 | QwQwQW10 2026-03-11 | 3/150 |
|