| 1 | 1/1 | 返回列表 |
| 查看: 492 | 回復(fù): 0 | ||
kndtmsl新蟲 (初入文壇)
|
[求助]
用espresso計算硫化鎢的能帶結(jié)構(gòu),無法有效設(shè)置k路徑
|
|
大家好,我初學(xué)espresso,老師讓我計算硫化鎢的能帶結(jié)構(gòu)。 第一個輸入文件是ws2.cry.scf.in,用來完成scf計算。 &control calculation = 'scf' restart_mode='from_scratch', prefix='ws2cry' tstress = .true. tprnfor = .true. pseudo_dir = '.', outdir='.' / &system ibrav=4, celldm(1)= 3.17, celldm(3) =3.9093, nat=6, ntyp= 2, ecutwfc =60.0, / &electrons diagonalization='david' mixing_mode = 'plain' mixing_beta = 0.7 conv_thr = 1.0d-8 / ATOMIC_SPECIES W 183.84 W.pbe-nsp-van.UPF S 32.065 S.pbe-van_bm.UPF ATOMIC_POSITIONS S 0.33333 0.66667 0.6225 S 0.33333 0.66667 -0.1225 S -0.33333 -0.66667 -0.6225 S -0.33333 -0.66667 0.1225 W -0.33333 -0.66667 -0.25 W 0.33333 0.66667 0.25 K_POINTS (automatic) 15 15 3 1 1 1 EOF 第二個輸入文件是ws2.cry.band.in,用來計算能帶。 &control calculation='bands' pseudo_dir = '.', outdir='.', prefix='ws2cry', / &system ibrav= 4, celldm(1)= 3.17, celldm(3) =3.9093, nat= 6, ntyp= 2, ecutwfc =60.0, ecutrho = 400.0, nbnd = 28, / &electrons diagonalization='david', / ATOMIC_SPECIES W 183.84 W.pbe-nsp-van.UPF S 32.065 S.pbe-van_bm.UPF ATOMIC_POSITIONS S 0.33333 0.66667 0.6225 S 0.33333 0.66667 -0.1225 S -0.33333 -0.66667 -0.6225 S -0.33333 -0.66667 0.1225 W -0.33333 -0.66667 -0.25 W 0.33333 0.66667 0.25 K_POINTS 2 0.0 0.0 0.0 4.0 0.36 0.0 0.0 1.0 EOF 還有bands.x的輸入文件ws2.cry.bands.in: &bands prefix = 'ws2cry', outdir = '.', filband = 'ws2cry.dat', lsym=.true, / 先執(zhí)行pw.x <ws2.cry.scf.in> ws2.cry.scf.out,然后執(zhí)行pw.x <ws2.cry.band.in> ws2.cry.band.out,最后執(zhí)行bands.x <ws2.cry.bands.in> ws2.cry.bands.out,其中ws2.cry.band.in的k points被完全無視了,我試過各種不同組合的k points,但輸出文件ws2cry.dat中的k points數(shù)量非常多,經(jīng)過比對是ws2.cry.scf.in中自動生成的k points。另外輸出文件有26個能帶,和我的設(shè)定不一樣(nbnd=28)。我之前用espresso算過銅的能帶結(jié)構(gòu),輸出文件的k points與能帶的輸入文件完全一致,為什么這一次輸出文件直接無視了能帶的輸入文件,而用的是能帶的輸入文件呢?還有一個細(xì)節(jié)就是做scf計算的時間很長,而計算能帶結(jié)構(gòu)用的的時間卻很短,幾乎是秒出結(jié)果。 |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +3 | 炸呀炸呀炸薯條 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 08工學(xué)調(diào)劑 +7 | 用戶573181 2026-03-20 | 11/550 |
|
|
[考研] 350求調(diào)劑 +6 | weudhdk 2026-03-19 | 6/300 |
|
|
[考研] 306求調(diào)劑 +9 | chuanzhu川燭 2026-03-18 | 9/450 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
[考研]
|
2117205181 2026-03-21 | 8/400 |
|
|
[考研] 315分,誠求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +5 | 我愛生物生物愛?/a> 2026-03-17 | 5/250 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
[考研] 一志愿華南師大 070300(化學(xué))304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|