justinzhouke

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 3.5
紅花: 1
帖子: 24
在線: 6.5小時
蟲號: 3513863
注冊: 2014-11-02
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[交流] 石墨烯聲子譜Gamma點附近虛頻，怎么處理已有4人參與

計算石墨烯的聲子譜，計算發(fā)現(xiàn)Gamma點附近出現(xiàn)虛頻（如下)。不知道怎么處理。

  1 0.000000    0.0000 0.0000 0.0000  879.5294 1560.4160 1560.4160
  2 0.022223    -3.7499 41.7646 65.5168  879.2473 1559.8164 1561.6153
  3 0.044447    -6.5222 83.4849  130.6515  878.3968 1557.9969 1565.0765
  4 0.066670    -6.7089  125.1109  195.0346  876.9649 1554.9026 1570.4069
  5 0.088893    6.6043  166.5984  258.3457  874.9310 1550.4615 1577.0073
  6 0.111117    17.6058  207.9030  320.3191  872.2675 1544.6098 1584.1435
  7 0.133340    29.2947  248.9746  380.7432  868.9412 1537.3218 1591.0273

基本設(shè)置包括PP，是按照論文里面來的（PHYSICAL REVIEW B 80, 033406 2009）。
我做過ecut和q-mesh的收斂性測試，q-mesh在16的時候就收斂了。Ecut是55Ry收斂。

In the present analysis, we take advantage of the PWSCF
package of the QUANTUM-ESPRESSO distribution13 to model
the characteristic graphene lattice dynamics. Specifically,
an ultrasoft pseudopotential in the local-density
approximation14 is used with a plane-wave expansion up to
55 Ry and 32x32x1 Monkhorst-Pack meshes adopted for
Brillouin-zone sampling with 0.02 Ry Fermi-Dirac smearing
in the electronic occupation. The calculation also employs
the theoretical lattice parameter a=2.46 Å, which is in excellent
agreement with the experimental value at 300 K. For
bilayer graphene, the Bernal stacking is considered with the
interlayer distance set to the theoretical minimum of
3.36 Å.

輸入文件
Gra.scf.in:
1  &control
  2       calculation = 'scf',
  3       restart_mode= 'from_scratch',
  4       pseudo_dir  = '/home/xuzp/bin/lib/pwscf_pot',
  5       outdir    = './tmp',
  6       prefix='graphene',
  7       !disk_io    = 'none',
  8       tstress = .true.
  9       tprnfor = .true.
10       wf_collect=.true.
11 /
12 &system
13    ibrav          = 4
14    celldm(1)    = 4.64873
15    celldm(3)    = 8.13
16    nat          = 2
17    ntyp          = 1
18    ecutwfc       = 55
19    ecutrho       = 550
20    occupations    = 'smearing'
21    smearing       = 'fd'
22    degauss       = 0.02
23 /
24 &electrons
25    diagonalization = 'david'
26    mixing_beta    = 0.7
27    conv_thr       = 1.0d-10
28 /
29 ATOMIC_SPECIES
30 C  12.01078 C.pbe-rrkjus.UPF  !C.pz-vbc.UPF
31 ATOMIC_POSITIONS crystal
32 C 0.6666666666666667 0.3333333333333333 0.5000000000000000
33 C 0.3333333333333333 0.6666666666666667 0.5000000000000000
34 K_POINTS automatic
35 32 32 1 0 0 0

Gra.ph.in:
  1 phonon for graphene
  2 &inputph
  3       tr2_ph       =1.0d-18,
  4       prefix       ='graphene',
  5       epsil          =.false.,
  6       amass(1)       =12.01078,
  7       outdir       = './tmp',
  8       fildyn='graphene.dynG',
  9       !elph=.false.,
10       !trans=.true.,
11       ldisp=.true.
12       nq1=6, nq2=6, nq3=1
13 /

回復(fù)此樓

» 猜你喜歡

材料考研調(diào)劑已經(jīng)有7人回復(fù)
270調(diào)劑已經(jīng)有7人回復(fù)
275學(xué)碩081000服從調(diào)劑到其他專業(yè)，保不住本專業(yè)了已經(jīng)有6人回復(fù)
0710生物學(xué)，325求調(diào)劑已經(jīng)有3人回復(fù)
296材料專碩求調(diào)劑已經(jīng)有5人回復(fù)
385分生物學(xué)（071000）求調(diào)劑已經(jīng)有5人回復(fù)
377求調(diào)劑已經(jīng)有3人回復(fù)
【求調(diào)劑】新能源材料本科，一志愿211，初試321 已經(jīng)有6人回復(fù)
316求調(diào)劑已經(jīng)有11人回復(fù)
321求調(diào)劑已經(jīng)有9人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

vasp_DFTP+phonopy 出現(xiàn)負(fù)頻且有一個聲學(xué)支斜率為負(fù)值已經(jīng)有19人回復(fù)
用VASP-DFPT+phonopy計算石墨烯聲子譜求助已經(jīng)有37人回復(fù)
vasp+dfpt算聲子性質(zhì)用phonopy軟件處理數(shù)據(jù)時得到band圖真的不忍直視，它怎么了？？已經(jīng)有17人回復(fù)
phonopy聲子譜怎么列出gamma點頻率值（&qpoints) 已經(jīng)有10人回復(fù)
vasp5.2可以用linear response theory方法計算聲子譜，聲子態(tài)密度嗎？？已經(jīng)有15人回復(fù)
Vasp+phonopy 聲子譜出現(xiàn)離奇虛頻！已經(jīng)有25人回復(fù)
聲子譜虛頻求助已經(jīng)有10人回復(fù)
CASTEP 怎么計算石墨烯的聲子譜？已經(jīng)有8人回復(fù)
求助用MS計算聲子譜出錯的問題已經(jīng)有15人回復(fù)
為何用晶格常數(shù)的實驗值計算聲子？已經(jīng)有18人回復(fù)
用siesta畫聲子譜的問題已經(jīng)有12人回復(fù)
怎么消除虛頻，細(xì)節(jié) 已經(jīng)有10人回復(fù)
【求助】castep計算聲子譜，為什么不能進(jìn)行振動分析已經(jīng)有9人回復(fù)
【求助】有關(guān)聲子軟化的處理已經(jīng)有12人回復(fù)
【求助】聲子與相變已經(jīng)有28人回復(fù)
【求助】彈性不穩(wěn)定，聲子會穩(wěn)定嗎？已經(jīng)有18人回復(fù)

whoisyourdady

1樓 2015-05-05 14:11:11

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

KalaShayminS

鐵桿木蟲 (著名寫手)

應(yīng)助: 432 (碩士)
金幣: 6174.2
紅花: 56
帖子: 1088
在線: 379小時
蟲號: 2408111
注冊: 2013-04-08

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

聲子計算總是需要先結(jié)構(gòu)優(yōu)化，達(dá)到足夠高的精度，否則由于計算條件略有區(qū)別（贗勢等）會產(chǎn)生虛頻。

贊一下

回復(fù)此樓

2樓2015-05-05 14:32:33

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

lchchl

銅蟲 (小有名氣)

應(yīng)助: 7 (幼兒園)
金幣: 2462.4
散金: 128
紅花: 1
帖子: 229
在線: 160.4小時
蟲號: 1677234
注冊: 2012-03-09
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

Gamma點附近的一點虛頻，也許改一下聲學(xué)聲子求和規(guī)則就可以消除了。

贊一下

回復(fù)此樓

為什么小木蟲每天都是扣四個金幣�。。。。�！

3樓2015-05-18 15:10:06

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

magicmonk

至尊木蟲 (著名寫手)

應(yīng)助: 172 (高中生)
金幣: 11451.3
紅花: 24
帖子: 1124
在線: 3108.8小時
蟲號: 1191315
注冊: 2011-01-17
專業(yè): 理論和計算化學(xué)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

3樓: Originally posted by lchchl at 2015-05-18 15:10:06
Gamma點附近的一點虛頻，也許改一下聲學(xué)聲子求和規(guī)則就可以消除了。

此話怎講，求詳細(xì)說明求和規(guī)則

贊一下

回復(fù)此樓

我們最深的恐懼不是我們能力不夠，我們最深的恐懼是能力超越了極限。

4樓2015-05-18 15:33:50

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

lchchl

銅蟲 (小有名氣)

應(yīng)助: 7 (幼兒園)
金幣: 2462.4
散金: 128
紅花: 1
帖子: 229
在線: 160.4小時
蟲號: 1677234
注冊: 2012-03-09
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

4樓: Originally posted by magicmonk at 2015-05-18 15:33:50
此話怎講，求詳細(xì)說明求和規(guī)則...

有個參數(shù)叫zasr，這個是控制求和規(guī)則的，印象中可以設(shè)置的求和規(guī)則包括，no，simple，crystal等。

贊一下

回復(fù)此樓

為什么小木蟲每天都是扣四個金幣！�。。。�！

5樓2015-05-18 19:44:30

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

三葉蟲的腳印

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 2.5
帖子: 5
在線: 7.1小時
蟲號: 4299870
注冊: 2015-12-20
專業(yè): 實驗地球化學(xué)和計算地球化

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

樓主你好。請問你Gamma點的虛頻消掉了嗎？我最近計算聲子譜時也出現(xiàn)了相同的問題。Gamma點處總是有3個小虛頻。。

贊一下(1人)

回復(fù)此樓

6樓2016-12-22 20:17:47

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

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

[考研] 一志愿南昌大學(xué)324求調(diào)劑 +12	hanamiko 2026-03-27	12/600	2026-04-01 22:34 by 運氣yunqi
[考研] 266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑 +8	哇呼哼呼哼 2026-04-01	8/400	2026-04-01 21:41 by cal0306
[考研] 求調(diào)劑0703 +5	周嘉堯 2026-03-31	8/400	2026-04-01 20:32 by ltltkkk
[考研] 0856初試324分求調(diào)劑 +5	想上學(xué)求調(diào) 2026-04-01	5/250	2026-04-01 15:28 by 考研上岸呀一定?/a>
[考研] 311求調(diào)劑 +10	李芷新1 2026-03-31	10/500	2026-04-01 14:38 by chenqifeng666
[考研] 環(huán)境工程調(diào)劑 +9	hyzzzzzzz. 2026-04-01	9/450	2026-04-01 14:20 by salamander`
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +10	美味可樂雞翅 2026-03-26	12/600	2026-04-01 11:17 by syh9288
[考研] 化學(xué)0703 調(diào)劑 306分一志愿211 +12	26要上岸 2026-03-28	12/600	2026-04-01 11:10 by chemdavid
[考研] 環(huán)境工程 085701，267求調(diào)劑 +15	minht 2026-03-29	16/800	2026-04-01 10:13 by li_sujuan99
[考研] 339求調(diào)劑 +5	zjjkt 2026-03-31	5/250	2026-04-01 09:18 by JourneyLucky
[考研] 0856求調(diào)劑 +9	楒桉 2026-03-28	9/450	2026-03-31 19:06 by 暮澤12
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6	hanamiko 2026-03-29	6/300	2026-03-31 16:35 by hypershenger
[考研] 一志愿浙江大學(xué)工科動力工程370,數(shù)一121,專業(yè)課135，現(xiàn)在能去哪里 +3	080700調(diào)劑 2026-03-30	4/200	2026-03-31 12:00 by KLMY666
[考研] 085600，材料與化工321分求調(diào)劑 +10	大饞小子 2026-03-28	10/500	2026-03-29 23:35 by 飛行日記西
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6	零八# 2026-03-28	6/300	2026-03-29 21:20 by nanaliuyun
[考研] 332求92調(diào)劑 +8	蕉蕉123 2026-03-28	8/400	2026-03-29 10:46 by 周梓丹
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 315調(diào)劑 +4	0860求調(diào)劑 2026-03-26	5/250	2026-03-27 11:23 by wangjy2002
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 085602化學(xué)工程求調(diào)劑。 +4	平樂樂樂 2026-03-26	4/200	2026-03-26 17:57 by fmesaito