| 查看: 2840 | 回复: 17 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 shijing879 的 9 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会 | |||
[交流]
BSE+GW计算光学吸收谱(激子结合能)
|
|||
| 大家好，请问大家谁有BSE+GW计算光学吸收谱(激子结合能)的INCAR设置步骤，KPOINTS和注意事项，要求能计算成功，我要计算的是二维材料，官网和网上的一些教程就不要推荐了，我已经试过！ |
回复此楼» 猜你喜欢
交叉科学部支持青年基金，对三无青椒是个机会吗？
已经有5人回复
-
国家基金申请书模板内插入图片不可调整大小?
已经有6人回复
-
招博士
已经有6人回复
-
限项规定
已经有8人回复
-
国家级人才课题组招收2026年入学博士
已经有5人回复
-
Fe3O4@SiO2合成
已经有6人回复
-
青年基金C终止
已经有4人回复
-
青椒八年已不青，大家都被折磨成啥样了？
已经有7人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 抢金币啦！回帖就可以得到:
-
邀您投稿 Polymers 特刊
+1/477
-
“超分子材料交叉研究团队”联合诚聘博士后 [清华/吉大/复旦/北大]
+1/79
-
[求助] 半导体物理PPT课件 东华理工机电学院 彭新村老师
+1/71
-
关于本子打包
+1/67
-
北京理工大学 珠海校区全职院士招数名博士生--申请考核制-半导体、光学、微电子
+2/64
-
大连工业大学 超临界流体技术团队（纺材学院）招收2026级“申请-考核制”博士生
+1/58
-
东北大学-招收2026年硕士研究生2-3名（金属材料3D打印方向）
+5/40
-
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生
+1/30
-
招聘农用化学产品销售一名，须具备良好的英语口语，以便拓展海外市场。
+1/14
-
辽宁材料实验室框架复合材料课题组招收联合培养研究生（长期有效）
+2/12
-
深圳大学信息功能聚合物电介质方向“申请-考核制”博士生招生
+2/10
-
电子科技大学崔春华课题组招收物理化学背景博士生1名-申请考核制
+2/8
-
华南理工大学宋波教授招聘材料和化学方向博士后（长期有效）
+1/6
-
苏州大学医学院纳米生物医学方向招收申请考核制博士生1名
+1/4
-
澳科大诚招2026年秋季全奖博士研究生（药剂学/生物材料方向）
+1/4
-
广东以色列理工学院博士/硕士招生-通过稀疏观测用数据驱动方法预测湍流
+1/3
-
【博士全奖 3–4 名｜材料 / 化学 / 机械｜双校资源】澳大利亚西悉尼大学
+1/3
-
【博士后/科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+1/2
-
华南理工大学宋波教授招收2026年博士生（二氧化碳有机、高分子和光电生物催化还原）
+1/2
-
尊敬的基金委领导，您好！
+1/1
★
shijing879(金币+1): 谢谢参与
shijing879(金币+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. |
6楼2023-08-28 10:16:03
|
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 |
8楼2023-08-29 02:56:10
|
•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. |
9楼2023-08-29 02:58:01
|
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. |
10楼2023-08-29 03:03:57
|
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 |
11楼2023-08-29 03:07:17
|
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. |
12楼2023-08-29 03:10:12
7楼2023-08-28 15:48:51
13楼2023-08-29 15:19:18
14楼2023-08-30 09:32:08
15楼2023-08-31 14:55:35
16楼2023-09-03 02:23:57
17楼2024-11-19 18:37:07
简单回复
2023-08-27 15:44
回复
shijing879(金币+1): 谢谢参与
toppaper3楼
2023-08-27 16:11
回复
tzynew4楼
2023-08-27 16:53
回复
shijing879(金币+1): 谢谢参与
2023-08-27 19:53
回复
shijing879(金币+1): 谢谢参与
XG-WUST18楼
2024-11-24 07:11
回复
up 发自小木虫IOS客户端