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

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

quby7503

銀蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 485.7
散金: 3713
紅花: 15
帖子: 1488
在線: 263.7小時
蟲號: 638688
注冊: 2008-10-28
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[求助] 求助ZnO單層聲子計算

剛剛接觸PWSCF，用ZnO單層做嘗試，計算的聲子總是不對，在gamma點(diǎn)處的計算結(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ī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
然后做聲子計算
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
/
換了幾個贗勢，計算結(jié)果都問題，請問：問題出現(xiàn)在哪啊？我都急死了

回復(fù)此樓

» 猜你喜歡

085601材料工程找調(diào)劑已經(jīng)有11人回復(fù)
總分293求調(diào)劑已經(jīng)有9人回復(fù)
334分一志愿武理材料求調(diào)劑已經(jīng)有5人回復(fù)
327求調(diào)劑已經(jīng)有6人回復(fù)
329求調(diào)劑已經(jīng)有10人回復(fù)
275求調(diào)劑已經(jīng)有15人回復(fù)
291求調(diào)劑已經(jīng)有6人回復(fù)
309求調(diào)劑已經(jīng)有7人回復(fù)
297求調(diào)劑已經(jīng)有11人回復(fù)
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有6人回復(fù)

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

單層和多層石墨電子譜和聲子譜研究已經(jīng)有66人回復(fù)
【文獻(xiàn)學(xué)習(xí)】Li,Na摻雜ZnO空穴的雜化密度泛函計算(PRB，2009) 已經(jīng)有21人回復(fù)
ZnO的聲子計算已經(jīng)有18人回復(fù)

1樓 2013-11-03 10:01:28

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

huai

木蟲 (職業(yè)作家)

應(yīng)助: 4 (幼兒園)
金幣: 303.1
散金: 104
紅花: 6
帖子: 4289
在線: 550.6小時
蟲號: 420388
注冊: 2007-07-09
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

怎么不對了？沒看出問題。

贊一下

回復(fù)此樓

2樓2013-11-03 11:09:43

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

quby7503

銀蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 485.7
散金: 3713
紅花: 15
帖子: 1488
在線: 263.7小時
蟲號: 638688
注冊: 2008-10-28
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

2樓: Originally posted by huai at 2013-11-03 11:09:43
怎么不對了？沒看出問題。

前三個頻率對應(yīng)著聲學(xué)支，在gamma點(diǎn)時，頻率應(yīng)該為0，即使算上數(shù)值誤差的話，也應(yīng)該是個很小很小的數(shù)值才是。但我計算的頻率都太大了

贊一下

回復(fù)此樓

3樓2013-11-03 11:14:13

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

huai

木蟲 (職業(yè)作家)

應(yīng)助: 4 (幼兒園)
金幣: 303.1
散金: 104
紅花: 6
帖子: 4289
在線: 550.6小時
蟲號: 420388
注冊: 2007-07-09
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

★ ★
franch: 金幣+2, 謝謝回帖交流 2013-11-04 18:19:18

引用回帖:

3樓: Originally posted by quby7503 at 2013-11-03 11:14:13
前三個頻率對應(yīng)著聲學(xué)支，在gamma點(diǎn)時，頻率應(yīng)該為0，即使算上數(shù)值誤差的話，也應(yīng)該是個很小很小的數(shù)值才是。但我計算的頻率都太大了...

我沒有太好的建議，我認(rèn)為這是數(shù)值誤差。前三個比起后三個要小很多。我只是做過對角化動力學(xué)矩陣得到聲子譜。對pwscf不了解，你問一下他們用過的人，看虛頻是多大的量級。

贊一下

回復(fù)此樓

4樓2013-11-03 16:41:39

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

future_wl

木蟲 (著名寫手)

應(yīng)助: 87 (初中生)
金幣: 3295.4
散金: 1881
紅花: 5
帖子: 1404
在線: 676.6小時
蟲號: 950891
注冊: 2010-01-30
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
quby7503: 金幣+3 2013-11-04 12:24:31
franch: 金幣+2, 謝謝回帖交流 2013-11-04 18:19:29

計算表面聲子比較復(fù)雜，你是否考慮了不同真空層厚度。

贊一下

回復(fù)此樓

未來就是現(xiàn)在

5樓2013-11-04 08:28:04

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

quby7503

銀蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 485.7
散金: 3713
紅花: 15
帖子: 1488
在線: 263.7小時
蟲號: 638688
注冊: 2008-10-28
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

5樓: Originally posted by future_wl at 2013-11-04 08:28:04
計算表面聲子比較復(fù)雜，你是否考慮了不同真空層厚度。

謝謝，真空層是經(jīng)過嚴(yán)格優(yōu)化的。

贊一下

回復(fù)此樓

6樓2013-11-04 12:24:02

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

倔強(qiáng)的堅果墻

鐵桿木蟲 (正式寫手)

應(yīng)助: 39 (小學(xué)生)
金幣: 8437.1
散金: 52
紅花: 12
帖子: 463
在線: 878.6小時
蟲號: 1251958
注冊: 2011-04-01
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
quby7503: 金幣+3 2013-11-05 12:01:58

沒有計算過ZnO，不過我建議你從以下幾個方面入手：1.設(shè)置encut的值（需要測試）2.檢查ecutwfc設(shè)置是否足夠（印象中g(shù)raphene的設(shè)置達(dá)到500Ry）3.K點(diǎn)是否測試達(dá)到收斂 4.結(jié)構(gòu)是否充分弛豫

贊一下

回復(fù)此樓

像小強(qiáng)一樣活著

7樓2013-11-04 21:05:45

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

quby7503

銀蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 485.7
散金: 3713
紅花: 15
帖子: 1488
在線: 263.7小時
蟲號: 638688
注冊: 2008-10-28
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

7樓: Originally posted by 倔強(qiáng)的堅果墻 at 2013-11-04 21:05:45
沒有計算過ZnO，不過我建議你從以下幾個方面入手：1.設(shè)置encut的值（需要測試）2.檢查ecutwfc設(shè)置是否足夠（印象中g(shù)raphene的設(shè)置達(dá)到500Ry）3.K點(diǎn)是否測試達(dá)到收斂 4.結(jié)構(gòu)是否充分弛豫

謝謝，都是經(jīng)過測試的

贊一下

回復(fù)此樓

8樓2013-11-05 12:01:50

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

倔強(qiáng)的堅果墻

鐵桿木蟲 (正式寫手)

應(yīng)助: 39 (小學(xué)生)
金幣: 8437.1
散金: 52
紅花: 12
帖子: 463
在線: 878.6小時
蟲號: 1251958
注冊: 2011-04-01
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

8樓: Originally posted by quby7503 at 2013-11-05 12:01:50
謝謝，都是經(jīng)過測試的...

那運(yùn)用聲學(xué)支求和規(guī)則畫聲子譜看看，假如虛頻不是很大，也是可以接受的

贊一下

回復(fù)此樓

像小強(qiáng)一樣活著

9樓2013-11-05 16:36:05

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

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

返回列表

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

[考研] 327求調(diào)劑 +6	汲亦昊 2026-03-29	6/300	2026-03-29 13:40 by peike
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 352分-085602-一志愿985 +5	海納百川Ly 2026-03-29	5/250	2026-03-29 09:57 by Sjndkwm
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 085600 材料與化工 329分求調(diào)劑 +11	Mr. Z 2026-03-25	11/550	2026-03-29 08:14 by 無際的草原
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4	IceGuo 2026-03-26	5/250	2026-03-29 01:25 by griffith2014
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4	小張zxy張 2026-03-26	8/400	2026-03-28 19:25 by lbsjt
[考研] 一志愿北化085600材料專碩275\|有文章專利｜求調(diào)劑 +7	Micky11223 2026-03-25	7/350	2026-03-28 18:34 by 無際的草原
[考研] 394求調(diào)劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒
[考研] 085602 307分求調(diào)劑 +7	不知道叫什么！ 2026-03-26	7/350	2026-03-28 09:57 by 神馬都不懂
[考研] 286求調(diào)劑 +4	丟掉懶惰 2026-03-27	7/350	2026-03-28 08:07 by baoball
[考研] 285求調(diào)劑 +4	AZMK 2026-03-27	7/350	2026-03-27 20:59 by AZMK
[考研] 307求調(diào)劑 +8	超級伊昂大王 2026-03-24	9/450	2026-03-27 15:34 by 超級伊昂大王
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +4	美味可樂雞翅 2026-03-26	4/200	2026-03-27 15:27 by 星空星月
[考研] 求調(diào)劑 +3	劉柯@ 2026-03-24	4/200	2026-03-27 11:28 by shangxh
[考研] 286求調(diào)劑 +4	lim0922 2026-03-26	4/200	2026-03-27 10:28 by guoweigw
[考研] 271求調(diào)劑 +6	生如夏花… 2026-03-22	6/300	2026-03-26 16:48 by 張凱十八號
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4	鋼鐵大炮 2026-03-24	4/200	2026-03-26 16:28 by dick_runner
[考研] 化工專碩求調(diào)劑 +3	question挽風(fēng) 2026-03-24	3/150	2026-03-24 18:48 by jhhcooi
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106

24小時熱門版塊排行榜

[求助] 求助ZnO單層聲子計算

» 猜你喜歡

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

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

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

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