| 查看: 2173 | 回復(fù): 18 | |||
soiu2000銀蟲 (初入文壇)
|
[求助]
ZnO的聲子計(jì)算
|
|
最近新學(xué)pwscf,想計(jì)算ZnO Gamma點(diǎn)聲子。但是怎么算都和文獻(xiàn)上差很遠(yuǎn),我有的文獻(xiàn)上看到別人可以算的很接近。不知道是什么原因。還請(qǐng)各位高手指點(diǎn)一二。 附上輸入文件。 omega( 1 - 1) = -15.8 [cm-1] --> A_1 I+R omega( 2 - 3) = 41.6 [cm-1] --> E_1 I+R omega( 4 - 5) = 103.5 [cm-1] --> E_2 R omega( 6 - 6) = 234.2 [cm-1] --> B_2 omega( 7 - 8) = 315.2 [cm-1] --> E_1 I+R omega( 9 - 9) = 323.1 [cm-1] --> A_1 I+R omega( 10 - 11) = 369.1 [cm-1] --> E_2 R omega( 12 - 12) = 493.8 [cm-1] --> B_2 而文獻(xiàn)上 omega LO(A1)=574 [cm-1] omega LO(E1)=583 [cm-1] omega TO(A1)=380 [cm-1] omega TO(E1)=407 [cm-1] omega (E2)=437 [cm-1] omega (E2)=101 [cm-1] |
鐵桿木蟲 (正式寫手)
|
&electrons conv_thr = 1.0d-10 mixing_beta = 0.7 electron_maxstep=200 / ATOMIC_SPECIES O 15.999 O.pz-van_ak.UPF Zn 65.39 Zn.pz-van_ak.UPF ATOMIC_POSITIONS {crystal} Zn 0.000000000 0.000000000 -0.000969838 Zn 0.666666667 0.333333333 0.499030162 O 0.000000000 0.000000000 0.378963856 O 0.666666667 0.333333333 0.878963856 ATOMIC_SPECIES 和 ATOMIC_POSITIONS {crystal}中的原子的順序,顛倒了。是否可以顛倒呢? |
|
首先,我不明白你的pw的in文件中ecutwfc要設(shè)置這么高,K點(diǎn)網(wǎng)格又相對(duì)比較低?ph的in文件中你只算Gamma點(diǎn),為什么還要寫q網(wǎng)格?還有,你的結(jié)果中,有一點(diǎn)點(diǎn)虛頻的存在,前面的relax完全了嗎? 最后,我想問(wèn)你下,你用pwscf 4.3.1的并行效率這么樣?我自己最近也剛裝4.3.1 怎么感覺并行效率很低啊。CPU使用率都只有40%左右。你的呢? |
鐵桿木蟲 (正式寫手)
鐵桿木蟲 (正式寫手)
銀蟲 (初入文壇)
榮譽(yù)版主 (文壇精英)
super-emuch

銅蟲 (小有名氣)
|
首先確保你的優(yōu)化是對(duì)的,先做收斂測(cè)試,這是很嚴(yán)謹(jǐn)?shù)模矣X得你截?cái)嗳∧敲锤撸赡苡袃蓚(gè)原因1,沒(méi)做收斂測(cè)試,1,你可能以為截?cái)嗟膯挝皇莈V吧,應(yīng)該Ry。還有這些參數(shù):wf_collect=.true. forc_conv_thr= 1.0d-7 etot_conv_thr= 1.0d-7 nstep=100, electron_maxstep=200 一般是默認(rèn)值,就可以。 原子坐標(biāo)有問(wèn)題,閑心礦可以用MS生成坐標(biāo),在對(duì)其優(yōu)化。 晶格常數(shù)不用精確這么多為吧,三位足以。 |

銅蟲 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +5 | 十三加油 2026-03-21 | 5/250 |
|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +6 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-17 | 8/400 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變,給編輯又寫郵件了,沒(méi)回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|