| 5 | 1/1 | 返回列表 |
| 查看: 1497 | 回復(fù): 8 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
求助ZnO單層聲子計(jì)算
|
||
|
剛剛接觸PWSCF,用ZnO單層做嘗試,計(jì)算的聲子總是不對(duì),在gamma點(diǎn)處的計(jì)算結(jié)果如下: omega( 1) = -4.186512 [THz] = -139.647010 [cm-1] ( 0.666757 0.000000 0.019032 0.000000 0.000000 0.000000 ) ( 0.744729 0.000000 0.021257 0.000000 0.000000 0.000000 ) omega( 2) = -4.186512 [THz] = -139.647010 [cm-1] ( 0.019032 0.000000 -0.666757 0.000000 -0.000000 0.000000 ) ( 0.021257 0.000000 -0.744729 0.000000 0.000000 0.000000 ) omega( 3) = 1.330295 [THz] = 44.373869 [cm-1] ( -0.000000 0.000000 -0.000000 0.000000 -0.815519 0.000000 ) ( 0.000000 0.000000 0.000000 0.000000 -0.578730 0.000000 ) omega( 4) = 7.850469 [THz] = 261.863452 [cm-1] ( -0.000000 0.000000 -0.000000 0.000000 -0.171074 0.000000 ) ( 0.000000 0.000000 0.000000 0.000000 0.985258 0.000000 ) omega( 5) = 17.606727 [THz] = 587.297190 [cm-1] ( 0.000819 0.000000 0.263621 0.000000 -0.000000 0.000000 ) ( -0.002996 0.000000 -0.964621 0.000000 0.000000 0.000000 ) omega( 6) = 17.606727 [THz] = 587.297190 [cm-1] ( 0.263621 0.000000 -0.000819 0.000000 -0.000000 0.000000 ) ( -0.964621 0.000000 0.002996 0.000000 0.000000 0.000000 ) 計(jì)算過程是,首先做自洽計(jì)算(原子位置已經(jīng)經(jīng)過嚴(yán)格的優(yōu)化) &control calculation = 'scf' restart_mode='from_scratch', prefix='zno-001', tstress = .true. tprnfor = .true. pseudo_dir = '/home/LAB/byqu/Pseudpo-pwscf', outdir='./tmp' / &system ibrav= 4, a= 3.28 , c=15.0, nat= 2, ntyp= 2, ecutwfc = 325.0 , / &electrons diagonalization='david' mixing_mode = 'plain' mixing_beta = 0.7 conv_thr = 1.0d-8 / ATOMIC_SPECIES Zn 65.39 Zn.pbe-d-hgh.UPF O 15.9994 O.pbe-hgh.UPF ATOMIC_POSITIONS {crystal} Zn 0.6666666699999979 0.3333333300000021 0.5000000000000000 O 0.3333333300000021 0.6666666699999979 0.5000000000000000 K_POINTS {automatic} 9 9 1 0 0 0 然后做聲子計(jì)算 phonons for ZnO-monolayer &inputph tr2_ph = 1.0d-12 prefix = 'zno-001' outdir = './tmp' amass(1) = 65.39 amass(2) = 15.9994 fildyn = 'zno-001.dyn' trans = .true. ldisp = .true. nq1=15, nq2=15, nq3=1 / 換了幾個(gè)贗勢(shì),計(jì)算結(jié)果都問題,請(qǐng)問:?jiǎn)栴}出現(xiàn)在哪啊?我都急死了 |
木蟲 (著名寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-28 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | 爭(zhēng)取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 297求調(diào)劑 +6 | 田洪有 2026-03-26 | 6/300 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生(材料、化學(xué)、電化學(xué),環(huán)境) +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|