| 18 | 1/1 | 返回列表 |
| 查看: 2956 | 回復(fù): 17 | |||
| 【有獎交流】積極回復(fù)本帖子,參與交流,就有機會分得作者 shijing879 的 9 個金幣 ,回帖就立即獲得 1 個金幣,每人有 1 次機會 | |||
[交流]
BSE+GW計算光學(xué)吸收譜(激子結(jié)合能)
|
|||
| 大家好,請問大家誰有BSE+GW計算光學(xué)吸收譜(激子結(jié)合能)的INCAR設(shè)置步驟,KPOINTS和注意事項,要求能計算成功,我要計算的是二維材料,官網(wǎng)和網(wǎng)上的一些教程就不要推薦了,我已經(jīng)試過! |
» 搶金幣啦!回帖就可以得到:
+1/171
+1/88
+1/83
+1/82
+1/49
+1/37
+1/36
+1/34
+1/34
+1/33
+1/29
+1/11
+1/6
+1/6
+1/6
+1/4
+1/4
+1/3
+1/1
+1/1
|
VASP sucks! Goal: account for excitonic effects in the frequency dependent dielectric function • GW0+BSE: Solve the Bethe-Salpeter-Equation on top of GW0 WORKFLOW of a GW0+BSE calculation The workflow of the GW0+BSE calculation in job.sh consists of five consecutive steps: Step 1: a “standard" DFT groundstate calculation Step 2: compute additional DFT “virtual” orbitals (empty states): • Needs the WAVECAR from Step 1. Step 3: the GW0 calculation: • Needs the WAVECAR and WAVEDER files from Step 2. Step 4: (optional) use LOPTICS=.TRUE. to plot the IPA dielectric function using the GW0-QP energies instead of DFT eigenenergies • Needs the WAVECAR file from Step 3. Step 5: the BSE calculation: • Needs WAVECAR file from Step 3 and WAVEDER file from Step 2. |
|
Step 1: a “standard” DFT groundstate calculation KPOINTS (e.g. KPOINTS.6): 6x6x6 0 G 6 6 6 0 0 0 INCAR (copy INCAR.DFT to INCAR): PREC = Normal ENCUT = 250 Default Set cutoff energy to 250 eV ISMEAR = 0 SIGMA = 0.01 Gaussian smearing Set small smearing width KPAR = 2 Additional parallelization EDIFF = 1.E-8 Set tight convergence criterium |
|
•Step 2: compute additional DFT “virtual” orbitals (empty states) INCAR (copy INCAR.DIAG to INCAR): PREC = Normal ; ENCUT = 250 ALGO = Exact NELM = 1 Exact diagonalisation of the Hamiltonian Only 1 electronic step ISMEAR = 0 ; SIGMA = 0.01 NBANDS = 128 Keep 128 bands after diagonalisation LOPTICS = .TRUE. LPEAD = .TRUE. Compute derivative of the orbitals w.r.t. the Bloch wave vector (stored in the WAVEDER file) using the PEAD formalism OMEGAMAX = 40 N.B.: This calculation needs the orbitals (WAVECAR file) written in Step 1. |
|
Step 3: RPA quasiparticles: single-shot GW (G0W0) INCAR (copy INCAR.GW0 INCAR): PREC = Normal ; ENCUT = 250 ALGO = GW0 ; NELM = 1 Request G0W0 calculation ISMEAR = 0 ; SIGMA = 0.01 ENCUTGW = 150 NOMEGA = 50 OMEAGTL = 280 PRECFOCK = Normal Eenergy cutoff for response functions Number of points for frequency integration NBANDS = 128 Use the same #-of-bands as in Step 2., otherwise the WAVEDER file can not be read correctly NBANDSGW = 12 Compute QP energies only for the first 12 bands KPAR = 2 LWAVE = .TRUE. Write WAVECAR (contains QP-energies) N.B.: This calculation needs the orbitals (WAVECAR file) , and the derivative of the orbitals w.r.t. the Bloch wave vectors (WAVEDER file) written in Step 2. |
|
Step 4: (optional) plot IPA dielectric function using G0W0QP-energies INCAR (copy INCAR.NONE INCAR): ALGO = Nothing ; NELM = 1 Do nothing except read WAVECAR and proceed to post-processing NBANDS = 128 LWAVE = .FALSE. Do not write WAVECAR LOPTICS = .TRUE. ; LPEAD = .TRUE. Compute IPA dielectric function |
|
Step 5: the BSE calculation INCAR (copy INCAR.BSE INCAR): PREC = Normal ; ENCUT = 250 ALGO = BSE ANTIRES = 0 Request BSE calculation Use the Tamm-Dancoff approximation ISMEAR = 0 ; SIGMA = 0.01 ENCUTGW = 150 Energy cutoff for response functions NBANDS = 128 Use the same #-of-bands as in Step 2., otherwise the WAVEDER file can not be read correctly NBANDSO = 4 NBANDSV = 8 Setup BSE matrix for NBANDSO HOMOs and NBANDSV LUMOs OMEGAMAX = 20 PRECFOCK = Normal N.B.: This calculation needs the orbitals (WAVECAR file) of Step 3, and the derivative of the orbitals w.r.t. the Bloch wave vectors (WAVEDER file) written in Step 2. |
| 18 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854總分272 +5 | 打小就是老實人 2026-03-02 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 接收調(diào)劑 +15 | 津萌津萌 2026-03-02 | 23/1150 |
|
|
[考研] 304分材料專碩求調(diào)劑 +3 | qiuzhigril 2026-03-03 | 5/250 |
|
|
[基金申請]
|
xhuama 2026-03-02 | 10/500 |
|
|
[考研] 中國林科院林化所(南京)2026年招收化學(xué)/材料/環(huán)境工程等背景碩士研究生3名 +3 | realstar2006 2026-02-27 | 3/150 |
|
|
[考研] 085600材料與化工調(diào)劑 280分 +10 | yyqqhh 2026-03-03 | 10/500 |
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +4 | 青科11 2026-03-02 | 4/200 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +11 | February_Feb 2026-03-01 | 11/550 |
|
|
[考研] 284求調(diào)劑 +6 | 天下熯 2026-03-02 | 6/300 |
|
|
[考研] 沒上岸的看過來 +3 | tangxiaotian 2026-03-01 | 3/150 |
|
|
[考研]
材料工程專碩283求調(diào)劑
5+8
|
,? 2026-03-02 | 10/500 |
|
|
[考研] 080500材料科學(xué)與工程 +4 | 202114020319 2026-03-03 | 4/200 |
|
|
[考研]
材料325求調(diào)劑
30+3
|
mariusuki 2026-03-02 | 7/350 |
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
[論文投稿]
通訊作者寫誰,問題是你意想不到的問題
15+3
|
阿爾法啊 2026-03-01 | 3/150 |
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學(xué)生 2026-02-28 | 9/450 |
|
|
[考研] 材料工程274求調(diào)劑 +5 | Lilithan 2026-03-01 | 5/250 |
|
|
[考研] 261求調(diào)劑 +3 | 陸lh 2026-03-01 | 3/150 |
|
|
[考研] 292求調(diào)劑 +7 | yhk_819 2026-02-28 | 7/350 |
|
|
[考研] 調(diào)劑 +3 | 簡木ChuFront 2026-02-28 | 3/150 |
|