justinzhouke

新蟲 (初入文壇)

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

[交流] gamma點聲子虛頻已有5人參與

我在計算閃鋅礦結(jié)構(gòu)的InAs的gamma點頻率的時候，結(jié)算結(jié)果出奇的奇怪。具體如下：
   Diagonalizing the dynamical matrix

   q = ( 0.000000000 0.000000000 0.000000000 )

**************************************************************************
   omega( 1) =    -6.801498 [THz] = -226.873553 [cm-1]
( -0.502968  0.000000 -0.207252  0.000000 -0.106323  0.000000 )
(  0.755263  0.000000  0.311212  0.000000  0.159655  0.000000 )
   omega( 2) =    -6.801498 [THz] = -226.873553 [cm-1]
(  0.147751  0.000000 -0.088266  0.000000 -0.526891  0.000000 )
( -0.221864  0.000000  0.132541  0.000000  0.791185  0.000000 )
   omega( 3) =    -6.801498 [THz] = -226.873553 [cm-1]
(  0.180077  0.000000 -0.506449  0.000000  0.135339 -0.000000 )
( -0.270406  0.000000  0.760490 -0.000000 -0.203226  0.000000 )
   omega( 4) =    -0.136566 [THz] =    -4.555354 [cm-1]
(  0.006948  0.000000 -0.160492  0.000000  0.690824  0.000000 )
(  0.006906  0.000000 -0.159518  0.000000  0.686631  0.000000 )
   omega( 5) =    -0.136566 [THz] =    -4.555354 [cm-1]
(  0.617358  0.000000  0.341420  0.000000  0.073109  0.000000 )
(  0.613612  0.000000  0.339348  0.000000  0.072665  0.000000 )
   omega( 6) =    -0.136566 [THz] =    -4.555354 [cm-1]
( -0.349090  0.000000  0.600599 -0.000000  0.143042  0.000000 )
( -0.346972  0.000000  0.596954 -0.000000  0.142174  0.000000 )
**************************************************************************

gamma點虛頻居然如此之大。我以及做了收斂性測試和結(jié)構(gòu)優(yōu)化。
但是為何gamma點虛頻居然如此之大？。誰能幫助我一下。

輸入文件如下：
scf：
&control
      calculation = 'scf',
    restart_mode= 'from_scratch',
    pseudo_dir  = './',
    outdir    = './tmp',
      prefix='inas',
      tstress = .true.
      tprnfor = .true.
/
&system
ibrav=  0,
celldm(1)=11.450120305668,
nat=  2,
ntyp= 2,
ecutwfc =48,
/
&electrons
diagonalization='cg'
mixing_mode = 'plain'
mixing_beta = 0.7
conv_thr =  5.0d-12
/
ATOMIC_SPECIES
In  111.818  In.pbe-hgh.UPF
As 74.920  As.pbe-hgh.UPF

ATOMIC_POSITIONS crystal
In 0.00 0.00 0.00
As 0.25 0.25 0.25
K_POINTS automatic
9 9 9 0 0 0
CELL_PARAMETERS alat
0.000000000 0.508881713 0.508881713
0.508881713 0.000000000 0.508881713
0.508881713 0.508881713  -0.000000000

ph：
phonons of InAs
&inputph
   tr2_ph=1.0d-14,
   prefix='inas',
   epsil=.true.
      amass(1) =111.818,
       amass(2) =74.920,
   outdir='./tmp',
   fildyn='inas.dyn',
/
0.0 0.0 0.0

自恰輸出結(jié)果為：
。。。。
! total energy             =    -16.65528050 Ry
   Harris-Foulkes estimate =    -16.65528050 Ry
   estimated scf accuracy <       4.8E-12 Ry

   The total energy is the sum of the following terms:

   one-electron contribution =    1.42930748 Ry
   hartree contribution    =    1.95307784 Ry
   xc contribution          =    -4.59653835 Ry
   ewald contribution       =    -15.44112747 Ry

   convergence has been achieved in  34 iterations

   Forces acting on atoms (Ry/au):

   atom 1 type  1 force =    0.00000000 0.00000000 0.00000000
   atom 2 type  2 force =    0.00000000 0.00000000 0.00000000

   Total force =    0.000000    Total SCF correction =    0.000000

   entering subroutine stress ...

      total stress  (Ry/bohr**3)                (kbar)    P= 0.00
0.00000001 0.00000000 0.00000000       0.00    0.00    0.00
0.00000000 0.00000001  -0.00000000       0.00    0.00    -0.00
0.00000000  -0.00000000 0.00000001       0.00    -0.00    0.00

。。。。

回復(fù)此樓

» 猜你喜歡

材料考研調(diào)劑已經(jīng)有3人回復(fù)
材料調(diào)劑已經(jīng)有12人回復(fù)
英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑已經(jīng)有14人回復(fù)
085410 一志愿211 22408分?jǐn)?shù)359求調(diào)劑已經(jīng)有4人回復(fù)
271求調(diào)劑已經(jīng)有19人回復(fù)
385分生物學(xué)（071000）求調(diào)劑已經(jīng)有3人回復(fù)
一志愿安徽大學(xué)計算機(jī)科學(xué)與技術(shù)學(xué)碩，331分求調(diào)劑已經(jīng)有3人回復(fù)
318求調(diào)劑，計算材料方向已經(jīng)有8人回復(fù)
291求調(diào)劑已經(jīng)有25人回復(fù)
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑已經(jīng)有6人回復(fù)

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

有誰用過 phonon 這個軟件，希望交流下使用經(jīng)驗已經(jīng)有13人回復(fù)
聲子譜計算在Gamma點出0.2THz虛頻，可以忽略嗎？已經(jīng)有52人回復(fù)
石墨烯聲子譜Gamma點附近虛頻，怎么處理已經(jīng)有5人回復(fù)
用VASP-DFPT+phonopy計算石墨烯聲子譜求助已經(jīng)有37人回復(fù)
聲子譜�。≈x謝��！已經(jīng)有4人回復(fù)
Gamma點處出現(xiàn)比較大的虛頻說明什么已經(jīng)有3人回復(fù)
你們怎么處理Gamma的虛頻的？已經(jīng)有4人回復(fù)
聲子譜虛頻求助已經(jīng)有10人回復(fù)
degauss的收斂已經(jīng)有13人回復(fù)
哭了，各種虛頻！已經(jīng)有23人回復(fù)
用自己生成的贗勢算能帶效果很好，算聲子卻一團(tuán)糟已經(jīng)有11人回復(fù)
關(guān)于用超胞方法計算聲子的一些經(jīng)驗已經(jīng)有21人回復(fù)
怎么消除虛頻，細(xì)節(jié) 已經(jīng)有10人回復(fù)
求助如何分析castep聲子計算的結(jié)果已經(jīng)有7人回復(fù)
自旋體系的聲子計算應(yīng)該用什么軟件呢？已經(jīng)有5人回復(fù)

whoisyourdady

1樓 2015-09-02 13:59:03

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

afan0012

木蟲 (小有名氣)

應(yīng)助: 7 (幼兒園)
金幣: 4328
散金: 7
紅花: 6
帖子: 228
在線: 65.2小時
蟲號: 3028351
注冊: 2014-03-08
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

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

樓主是用什么軟件做的？我用Phonopy做聲子譜時，曾經(jīng)因為對稱性，超胞的奇偶等很多原因出現(xiàn)過虛頻

贊一下

回復(fù)此樓

2樓2015-09-04 10:55:04

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

justinzhouke

新蟲 (初入文壇)

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

引用回帖:

2樓: Originally posted by afan0012 at 2015-09-04 10:55:04
樓主是用什么軟件做的？我用Phonopy做聲子譜時，曾經(jīng)因為對稱性，超胞的奇偶等很多原因出現(xiàn)過虛頻

我用pwscf

回復(fù)此樓

whoisyourdady

3樓2015-09-04 14:31:11

已閱回復(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, 給個紅包，謝謝回帖

引用回帖:

對稱性和超胞的奇偶？求樓主詳細(xì)解釋一下

贊一下

回復(fù)此樓

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

4樓2015-09-05 11:28:12

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

yyyu200

金蟲 (小有名氣)

應(yīng)助: 27 (小學(xué)生)
金幣: 713
散金: 790
紅花: 4
帖子: 225
在線: 129.5小時
蟲號: 1760414
注冊: 2012-04-17
專業(yè): 物理學(xué)I

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

需要用超胞，原胞計算的結(jié)果由于平移對稱性破壞會有很大誤差

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

5樓2015-09-05 18:00:07

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

呂小強(qiáng)強(qiáng)

新蟲 (正式寫手)

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

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

引用回帖:

最近也遇到了同樣的問題，是不是用tolerance定義對稱性的時候遇到的？

贊一下

回復(fù)此樓

6樓2018-01-08 16:40:21

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

liuxinbao

鐵蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 38.5
紅花: 2
帖子: 46
在線: 10.5小時
蟲號: 4437346
注冊: 2016-02-26
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

我也碰到類似問題,有時候碰到不了,有時候碰到,同求

贊一下

回復(fù)此樓

7樓2018-02-01 21:14:49

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

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

[考研] 085801 總分275 本科新能源求調(diào)劑 +10	bradoner 2026-04-01	11/550	2026-04-01 22:57 by 來看流星雨10
[考研] 江蘇科技大學(xué)招材料研究生 +4	Su032713. 2026-04-01	5/250	2026-04-01 22:03 by cccchenso
[考研] 311求調(diào)劑 +11	藍(lán)月亮亮 2026-03-30	11/550	2026-04-01 16:33 by 七度不信任
[考研] 0703求調(diào)劑 +4	zizimo 2026-03-31	4/200	2026-04-01 16:04 by yanflower7133
[考研] 08工科275求調(diào)劑，可跨考。 +5	AaAa7420 2026-03-31	5/250	2026-04-01 15:21 by 159357hjz
[考研] 材料0856 英一數(shù)二 323 求調(diào)劑 +9	袁sy 2026-04-01	9/450	2026-04-01 14:30 by wangjy2002
[考研] 336求調(diào)劑 +3	kiyy 2026-04-01	3/150	2026-04-01 12:46 by wxiongid
[碩博家園] 博一被送出聯(lián)培感覺不適應(yīng)怎么辦 +3	全村的狗 2026-03-31	3/150	2026-04-01 10:44 by 328838485
[考研] 318求調(diào)劑 +8	七憶77 2026-04-01	8/400	2026-04-01 10:37 by Jaylen.
[考研] 【調(diào)劑】一志愿廈大生物與醫(yī)藥調(diào)劑 +3	Echo蝦米 2026-03-31	3/150	2026-04-01 08:40 by JourneyLucky
[考研] 289求調(diào)劑 +7	BrightLL 2026-03-29	7/350	2026-03-31 22:05 by 544594351
[考研] 375求調(diào)劑 +7	雨夏整夜 2026-03-29	7/350	2026-03-31 18:52 by xhai2011
[考研] 070300化學(xué)354求調(diào)劑 +15	101次希望 2026-03-28	15/750	2026-03-31 17:58 by jp9609
[考研] 254材料與化工求調(diào)劑 +3	翰冬林楠 2026-03-30	4/200	2026-03-31 17:53 by yishunmin
[考研] 一志愿西電085401數(shù)一英一299求調(diào)劑六級521 +4	愛吃大鴨梨 2026-03-31	4/200	2026-03-31 11:51 by 搏擊518
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-29	3/150	2026-03-31 10:06 by cal0306
[考研] 282求調(diào)劑 +3	呼吸都是減肥 2026-03-31	3/150	2026-03-31 10:00 by wangjy2002
[考研] 吉大生物學(xué)326分求調(diào)劑 +3	sunnyupup 2026-03-31	3/150	2026-03-31 09:28 by longlotian
[考研] 材料與化工304求B區(qū)調(diào)劑 +4	邱gl 2026-03-26	7/350	2026-03-30 08:39 by 探123
[考研] 265求調(diào)劑11408 +3	劉小鹿lu 2026-03-27	3/150	2026-03-27 20:53 by nihaoar