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

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

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

返回列表

soiu2000

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 58.2
帖子: 7
在線: 39.1小時(shí)
蟲號(hào): 233903
注冊(cè): 2006-03-30
專業(yè): 半導(dǎo)體材料

[求助] 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]

回復(fù)此樓

» 猜你喜歡

326求調(diào)劑已經(jīng)有8人回復(fù)
333求調(diào)劑已經(jīng)有7人回復(fù)
一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑已經(jīng)有3人回復(fù)
考研化學(xué)學(xué)碩調(diào)劑，一志愿985 已經(jīng)有7人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有21人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有15人回復(fù)
0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷已經(jīng)有14人回復(fù)
326求調(diào)劑已經(jīng)有4人回復(fù)
311求調(diào)劑已經(jīng)有13人回復(fù)
311求調(diào)劑已經(jīng)有7人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【文獻(xiàn)學(xué)習(xí)】Li,Na摻雜ZnO空穴的雜化密度泛函計(jì)算(PRB，2009) 已經(jīng)有21人回復(fù)
新手求教pw算聲子譜提示缺少.dyn1文件已經(jīng)有9人回復(fù)
自己使用QHA時(shí)計(jì)算退出已經(jīng)有23人回復(fù)
【討論】abint計(jì)算ZnO帶隙差的太大已經(jīng)有26人回復(fù)
【求助】pwscf計(jì)算ZnO聲子時(shí)找不到對(duì)稱性已經(jīng)有7人回復(fù)
【求助】重金求計(jì)算ZnO或GaN等材料方面形成能，結(jié)合能，激活能的計(jì)算過程已經(jīng)有29人回復(fù)

1樓 2011-10-28 13:59:38

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

ridge8888

鐵桿木蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 8843.8
散金: 23
帖子: 411
在線: 137.4小時(shí)
蟲號(hào): 212298
注冊(cè): 2006-03-08
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

【答案】應(yīng)助回帖

★ ★
franch(金幣+2): 謝謝回帖交流 2011-10-31 21:27:34

&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}中的原子的順序，顛倒了。是否可以顛倒呢？

贊一下(1人)

回復(fù)此樓

2樓2011-10-28 14:43:10

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

fzuvivi

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 1168
散金: 93
紅花: 1
帖子: 224
在線: 330.1小時(shí)
蟲號(hào): 881678
注冊(cè): 2009-10-23
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

【答案】應(yīng)助回帖

★ ★
franch(金幣+2): 謝謝回帖交流 2011-10-31 21:27:21

首先，我不明白你的pw的in文件中ecutwfc要設(shè)置這么高，K點(diǎn)網(wǎng)格又相對(duì)比較低？ph的in文件中你只算Gamma點(diǎn)，為什么還要寫q網(wǎng)格？還有，你的結(jié)果中，有一點(diǎn)點(diǎn)虛頻的存在，前面的relax完全了嗎？

最后，我想問你下，你用pwscf 4.3.1的并行效率這么樣？我自己最近也剛裝4.3.1 怎么感覺并行效率很低啊。CPU使用率都只有40%左右。你的呢？

贊一下(1人)

回復(fù)此樓

3樓2011-10-28 18:54:28

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

huazhorg

鐵桿木蟲 (正式寫手)

應(yīng)助: 21 (小學(xué)生)
金幣: 4953.8
散金: 5
紅花: 2
帖子: 381
在線: 363.2小時(shí)
蟲號(hào): 629899
注冊(cè): 2008-10-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★
franch(金幣+2): 謝謝回帖交流 2011-10-31 21:27:57

引用回帖:

1樓: Originally posted by soiu2000 at 2011-10-28 13:59:38:
最近新學(xué)pwscf，想計(jì)算ZnO Gamma點(diǎn)聲子。但是怎么算都和文獻(xiàn)上差很遠(yuǎn)，我有的文獻(xiàn)上看到別人可以算的很接近。不知道是什么原因。還請(qǐng)各位高手指點(diǎn)一二。
附上輸入文件。

omega( 1 - 1) = -15.8 [cm ...

感覺你的 ecutwfc = 500.0, ecutrho = 1200.0設(shè)置太高了，一般來說要做個(gè)收斂測(cè)試，獲取適當(dāng)?shù)慕財(cái)�。就你的情況，以我的經(jīng)驗(yàn)，你分別設(shè)置成30，350試試。

贊一下(1人)

回復(fù)此樓

4樓2011-10-31 19:24:27

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

huazhorg

鐵桿木蟲 (正式寫手)

應(yīng)助: 21 (小學(xué)生)
金幣: 4953.8
散金: 5
紅花: 2
帖子: 381
在線: 363.2小時(shí)
蟲號(hào): 629899
注冊(cè): 2008-10-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★
franch(金幣+2): 謝謝回帖交流 2011-10-31 21:28:08

引用回帖:

3樓: Originally posted by fzuvivi at 2011-10-28 18:54:28:
首先，我不明白你的pw的in文件中ecutwfc要設(shè)置這么高，K點(diǎn)網(wǎng)格又相對(duì)比較低？ph的in文件中你只算Gamma點(diǎn)，為什么還要寫q網(wǎng)格？還有，你的結(jié)果中，有一點(diǎn)點(diǎn)虛頻的存在，前面的relax完全了嗎？

最后，我想問你下 ...

顯然是你安裝的問題。就我的情況看，CPU利用率一般在97～99%。

贊一下(1人)

回復(fù)此樓

5樓2011-10-31 19:25:56

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

soiu2000

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 58.2
帖子: 7
在線: 39.1小時(shí)
蟲號(hào): 233903
注冊(cè): 2006-03-30
專業(yè): 半導(dǎo)體材料

感謝大家的幫助。ecutwfc給的太大了知識(shí)因?yàn)橄朐囋嚫呓財(cái)嗄苁欠衲芴岣呗曌拥挠?jì)算精度。
問題在于不論我怎樣提高自洽計(jì)算的精度，聲子計(jì)算的結(jié)果都沒有顯著改善。是不是因?yàn)槌涄I勢(shì)的原因？但看文檔中給的Si的例子也是用超軟贗勢(shì)，結(jié)果很好。所以很迷惑。已經(jīng)困在這個(gè)問題上太久了。始終找不到問題所在。

贊一下

回復(fù)此樓

6樓2011-10-31 23:15:34

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

w76990480

榮譽(yù)版主 (文壇精英)

super-emuch

應(yīng)助: 47 (小學(xué)生)
貴賓: 13.729
金幣: 109757
散金: 29386
紅花: 337
沙發(fā): 50
帖子: 23155
在線: 3166.8小時(shí)
蟲號(hào): 942733
注冊(cè): 2010-01-15
專業(yè): 固體力學(xué)
管轄: 文獻(xiàn)求助

引用回帖:

6樓: Originally posted by soiu2000 at 2011-10-31 23:15:34:
感謝大家的幫助。ecutwfc給的太大了知識(shí)因?yàn)橄朐囋嚫呓財(cái)嗄苁欠衲芴岣呗曌拥挠?jì)算精度。
問題在于不論我怎樣提高自洽計(jì)算的精度，聲子計(jì)算的結(jié)果都沒有顯著改善。是不是因?yàn)槌涄I勢(shì)的原因？但看文檔中給的Si的例 ...

請(qǐng)檢查一下你的vc-realx的過程，首先確保你scf的結(jié)構(gòu)是最穩(wěn)定的，不穩(wěn)定很容易出現(xiàn)虛頻！

贊一下

回復(fù)此樓

專業(yè)文獻(xiàn)求助

7樓2011-11-01 08:13:58

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

hanyanli0475

銅蟲 (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 173.1
散金: 32
紅花: 3
帖子: 296
在線: 130.5小時(shí)
蟲號(hào): 1160645
注冊(cè): 2010-12-02
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★
賀儀(金幣+2): 多謝指教！ 2011-11-03 21:57:34

首先確保你的優(yōu)化是對(duì)的，先做收斂測(cè)試，這是很嚴(yán)謹(jǐn)?shù)�，我覺得你截?cái)嗳∧敲锤撸赡苡袃蓚€(gè)原因1，沒做收斂測(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)有問題，閑心礦可以用MS生成坐標(biāo)，在對(duì)其優(yōu)化。
晶格常數(shù)不用精確這么多為吧，三位足以。

贊一下(1人)

回復(fù)此樓

什么時(shí)候能算完�。。�！

8樓2011-11-01 10:45:39

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

hanyanli0475

銅蟲 (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 173.1
散金: 32
紅花: 3
帖子: 296
在線: 130.5小時(shí)
蟲號(hào): 1160645
注冊(cè): 2010-12-02
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

贗勢(shì)可能會(huì)影響，我現(xiàn)在也是這個(gè)問題，用模守恒聲子圖還可以，但是用超軟就是負(fù)的。

贊一下

回復(fù)此樓

什么時(shí)候能算完�。。�！

9樓2011-11-01 10:48:34

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

fzuvivi

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 1168
散金: 93
紅花: 1
帖子: 224
在線: 330.1小時(shí)
蟲號(hào): 881678
注冊(cè): 2009-10-23
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

4樓: Originally posted by huazhorg at 2011-10-31 19:24:27:
感覺你的 ecutwfc = 500.0, ecutrho = 1200.0設(shè)置太高了，一般來說要做個(gè)收斂測(cè)試，獲取適當(dāng)?shù)慕財(cái)�。就你的情況，以我的經(jīng)驗(yàn)，你分別設(shè)置成30，350試試。

我原來用的是4.1.2的，直接configure然后make all，其CPU利用率能達(dá)到95以上。同樣，我4.3.1也是直接configure然后make all，并行環(huán)境是有的，就是利用率很低。望多指教。。我的并行是mpich2.

贊一下

回復(fù)此樓

10樓2011-11-03 09:57:21

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

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

返回列表

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

[考研] 311求調(diào)劑 +12	冬十三 2026-03-15	13/650	2026-03-21 19:07 by ColorlessPI
[考研] 材料學(xué)碩301分求調(diào)劑 +6	Liyouyumairs 2026-03-21	6/300	2026-03-21 17:42 by JourneyLucky
[考研] 336求調(diào)劑 +5	rmc8866 2026-03-21	5/250	2026-03-21 17:24 by 學(xué)員8dgXkO
[考研] 0805材料320求調(diào)劑 +3	深海物語 2026-03-20	3/150	2026-03-21 15:46 by 無際的草原
[考研] 307求調(diào)劑 +10	冷笙123 2026-03-17	10/500	2026-03-21 01:54 by JourneyLucky
[考研] 321求調(diào)劑 +9	何潤(rùn)采123 2026-03-18	11/550	2026-03-20 23:19 by JourneyLucky
[考研] 287求調(diào)劑 +7	晨昏線與星海 2026-03-19	8/400	2026-03-20 22:19 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡(jiǎn)木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	3/150	2026-03-20 20:35 by 學(xué)員8dgXkO
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82

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

[求助] ZnO的聲子計(jì)算

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: