版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請盡快對帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

【有獎(jiǎng)交流】積極回復(fù)本帖子，參與交流，就有機(jī)會(huì)分得作者 shijing879 的 9 個(gè)金幣，回帖就立即獲得 1 個(gè)金幣，每人有 1 次機(jī)會(huì)

shijing879

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 477
帖子: 14
在線: 5.2小時(shí)
蟲號(hào): 13793021

[交流] BSE+GW計(jì)算光學(xué)吸收譜(激子結(jié)合能)

大家好，請問大家誰有BSE+GW計(jì)算光學(xué)吸收譜(激子結(jié)合能)的INCAR設(shè)置步驟，KPOINTS和注意事項(xiàng)，要求能計(jì)算成功，我要計(jì)算的是二維材料，官網(wǎng)和網(wǎng)上的一些教程就不要推薦了，我已經(jīng)試過！

回復(fù)此樓

» 猜你喜歡

307求調(diào)劑已經(jīng)有4人回復(fù)
收調(diào)劑已經(jīng)有3人回復(fù)
315食品工程求調(diào)劑已經(jīng)有6人回復(fù)
研究生招生已經(jīng)有6人回復(fù)
材料與化工，291，求調(diào)劑已經(jīng)有12人回復(fù)
0817一志愿蘇州大學(xué)280復(fù)試調(diào)劑已經(jīng)有8人回復(fù)
求調(diào)劑一志愿蘇州大學(xué)，0856化工323分 | 本科應(yīng)化 | 有專利/競賽/科研助手經(jīng)歷 | 已經(jīng)有9人回復(fù)
085602化工求調(diào)劑已經(jīng)有4人回復(fù)
083000環(huán)境科學(xué)與工程調(diào)劑已經(jīng)有6人回復(fù)
求材料讀博院校已經(jīng)有9人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2023-08-27 13:28:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

wangyikeco

木蟲 (正式寫手)

應(yīng)助: 84 (初中生)
金幣: 7026
帖子: 542
在線: 111.3小時(shí)
蟲號(hào): 25598800

★
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.

贊一下(1人)

回復(fù)此樓

6樓2023-08-28 10:16:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangyikeco

木蟲 (正式寫手)

應(yīng)助: 84 (初中生)
金幣: 7026
帖子: 542
在線: 111.3小時(shí)
蟲號(hào): 25598800

引用回帖:

7樓: originally posted by shijing879 at 2023-08-28 15:48:51
這些我都知道，請問是否有每一步具體的incar，kpoints， thanks in advance!...

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

贊一下(1人)

回復(fù)此樓

8樓2023-08-29 02:56:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangyikeco

木蟲 (正式寫手)

應(yīng)助: 84 (初中生)
金幣: 7026
帖子: 542
在線: 111.3小時(shí)
蟲號(hào): 25598800

•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.

贊一下(1人)

回復(fù)此樓

9樓2023-08-29 02:58:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangyikeco

木蟲 (正式寫手)

應(yīng)助: 84 (初中生)
金幣: 7026
帖子: 542
在線: 111.3小時(shí)
蟲號(hào): 25598800

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.

贊一下(1人)

回復(fù)此樓

10樓2023-08-29 03:03:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangyikeco

木蟲 (正式寫手)

應(yīng)助: 84 (初中生)
金幣: 7026
帖子: 542
在線: 111.3小時(shí)
蟲號(hào): 25598800

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

贊一下(1人)

回復(fù)此樓

11樓2023-08-29 03:07:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangyikeco

木蟲 (正式寫手)

應(yīng)助: 84 (初中生)
金幣: 7026
帖子: 542
在線: 111.3小時(shí)
蟲號(hào): 25598800

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.

贊一下(1人)

回復(fù)此樓

12樓2023-08-29 03:10:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shijing879

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 477
帖子: 14
在線: 5.2小時(shí)
蟲號(hào): 13793021

引用回帖:

6樓: Originally posted by wangyikeco at 2023-08-28 10:16:03
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
...

這些我都知道，請問是否有每一步具體的INCAR，KPOINTS， Thanks in advance!

贊一下

回復(fù)此樓

7樓2023-08-28 15:48:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shijing879

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 477
帖子: 14
在線: 5.2小時(shí)
蟲號(hào): 13793021

引用回帖:

12樓: Originally posted by wangyikeco at 2023-08-29 03:10:12
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.0 ...

這些我也知道，都是官網(wǎng)上面的，我已經(jīng)按照這個(gè)算過了，沒成功。請問你有沒自己親身計(jì)算成功的一些算例。

贊一下

回復(fù)此樓

13樓2023-08-29 15:19:18

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangyikeco

木蟲 (正式寫手)

應(yīng)助: 84 (初中生)
金幣: 7026
帖子: 542
在線: 111.3小時(shí)
蟲號(hào): 25598800

引用回帖:

13樓: Originally posted by shijing879 at 2023-08-29 15:19:18
這些我也知道，都是官網(wǎng)上面的，我已經(jīng)按照這個(gè)算過了，沒成功。請問你有沒自己親身計(jì)算成功的一些算例。...

What version of VASP are you doing?

贊一下

回復(fù)此樓

14樓2023-08-30 09:32:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shijing879

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 477
帖子: 14
在線: 5.2小時(shí)
蟲號(hào): 13793021

引用回帖:

14樓: Originally posted by wangyikeco at 2023-08-30 09:32:08
What version of VASP are you doing?...

我用的是vasp6.3.2

贊一下

回復(fù)此樓

15樓2023-08-31 14:55:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangyikeco

木蟲 (正式寫手)

應(yīng)助: 84 (初中生)
金幣: 7026
帖子: 542
在線: 111.3小時(shí)
蟲號(hào): 25598800

引用回帖:

15樓: Originally posted by shijing879 at 2023-08-31 14:55:35
我用的是vasp6.3.2...

我的老版本上似乎沒問題
你升級(jí)到最新的6.4.2
或者把你的script發(fā)給我, zrli3 AT hotmail dot com

贊一下

回復(fù)此樓

16樓2023-09-03 02:23:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

理想的美好

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 89.3
帖子: 14
在線: 16.8小時(shí)
蟲號(hào): 29250578

★
shijing879(金幣+1): 謝謝參與

您好，樓主，我想請教您如何用 GW+BSE 來計(jì)算激子結(jié)合能，同時(shí)給它畫成能帶形式

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

17樓2024-11-19 18:37:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

簡單回復(fù)

shenrenren2樓

2023-08-27 15:44 回復(fù)

shijing879(金幣+1): 謝謝參與

toppaper3樓

2023-08-27 16:11 回復(fù)

shijing879(金幣+1): 謝謝參與

發(fā)自小木蟲Android客戶端

tzynew4樓

2023-08-27 16:53 回復(fù)

shijing879(金幣+1): 謝謝參與

miaojiabing5樓

2023-08-27 19:53 回復(fù)

shijing879(金幣+1): 謝謝參與

XG-WUST18樓

2024-11-24 07:11 回復(fù)

up 發(fā)自小木蟲IOS客戶端

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 shijing879 的主題更新

返回列表

提示: 如果您在30分鐘內(nèi)回復(fù)過其他散金貼，則可能無法領(lǐng)取此貼金幣

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 求調(diào)劑一志愿蘇州大學(xué)，0856化工323分 \| 本科應(yīng)化 \| 有專利/競賽/科研助手經(jīng)歷 \| +7	橙子cyx 2026-03-06	9/450	2026-03-10 12:03 by kiokin
[考研] 085602化工求調(diào)劑 +4	董boxing 2026-03-10	4/200	2026-03-10 12:02 by kiokin
[考研] 278求調(diào)劑 +7	Gale1314 2026-03-06	7/350	2026-03-10 10:20 by allen-yin
[考研] 一志愿天津大學(xué)材料與化工275求調(diào)劑 10+4	穿只靴子 2026-03-07	15/750	2026-03-10 10:07 by 30660438
[考研] 085600材料與化工 326 求調(diào)劑 +4	熱愛生活ing 2026-03-09	4/200	2026-03-10 08:25 by Linda Hu
[考研] 085600材料與化工，一志愿廣州985，求調(diào)劑 +15	qqyyaill 2026-03-05	15/750	2026-03-10 08:23 by Linda Hu
[考研] 264求調(diào)劑 +5	thext 2026-03-03	5/250	2026-03-10 08:12 by 一切OK
[考研] 調(diào)劑 +4	調(diào)劑的考研學(xué)生 2026-03-09	4/200	2026-03-09 20:22 by wsxw
[考研] 醫(yī)學(xué)檢驗(yàn)專業(yè)26考研求調(diào)劑 +3	琪ya 2026-03-03	5/250	2026-03-08 22:56 by 清風(fēng)月
[考研] 267化工調(diào)劑求助 +7	聰少OZ 2026-03-04	7/350	2026-03-07 17:33 by 無際的草原
[考研] 085600材料與化工 292分求調(diào)劑 +6	程晴之 2026-03-06	6/300	2026-03-07 09:22 by 斬魂滴兔子！
[考研] 材料調(diào)劑 +10	ounce. 2026-03-04	12/600	2026-03-07 09:14 by Mornach1988
[考研] 085600材料與化工（高分子）290分求調(diào)劑 +9	wengyujian 2026-03-04	9/450	2026-03-06 15:59 by 聰明的大松鼠
[考研] 289求調(diào)劑 +3	馳光絕景 2026-03-04	4/200	2026-03-06 12:34 by 馳光絕景
[考研] 085600，一志愿鄭州大學(xué)，280分求調(diào)劑 +7	Wuqi725 2026-03-05	7/350	2026-03-06 08:12 by 楊楊楊紫
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-04	3/150	2026-03-05 10:39 by Iveryant
[考研] 材料085600 303求調(diào)劑 +7	1bygone 2026-03-04	7/350	2026-03-05 09:39 by oxidpl
[考研] 一志愿西交化工專碩288專業(yè)課93求調(diào)劑 +6	好運(yùn)好運(yùn)接接 2026-03-04	6/300	2026-03-05 09:36 by houyaoxu
[考研] 307求調(diào)劑 +4	xuyuan111 2026-03-04	5/250	2026-03-04 20:39 by kakakapanpan
[考研] 0857求調(diào)劑 +4	CDDD1013 2026-03-04	4/200	2026-03-04 16:08 by zhukairuo

24小時(shí)熱門版塊排行榜

[交流] BSE+GW計(jì)算光學(xué)吸收譜(激子結(jié)合能)

» 猜你喜歡

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

» 搶金幣啦！回帖就可以得到:

查看全部散金貼