| 5 | 1/1 | 返回列表 |
| 查看: 1499 | 回復(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)問:問題出現(xiàn)在哪。课叶技彼懒 |
鐵桿木蟲 (正式寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 330分求調(diào)劑 +4 | qzenlc 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +8 | 15188958825 2026-03-25 | 8/400 |
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +3 | Feisty。 2026-03-28 | 7/350 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
[考研]
材料學(xué)碩,求調(diào)劑
6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|