版塊導(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í)熱門版塊排行榜

返回列表

justinzhouke

新蟲 (初入文壇)

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

[交流] gamma點(diǎn)聲子虛頻已有5人參與

我在計(jì)算閃鋅礦結(jié)構(gòu)的InAs的gamma點(diǎn)頻率的時(shí)候，結(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點(diǎn)虛頻居然如此之大。我以及做了收斂性測(cè)試和結(jié)構(gòu)優(yōu)化。
但是為何gamma點(diǎn)虛頻居然如此之大？。誰能幫助我一下。

輸入文件如下：
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ù)此樓

» 猜你喜歡

085600材料與化工調(diào)劑 324分已經(jīng)有5人回復(fù)
344求調(diào)劑已經(jīng)有6人回復(fù)
266求調(diào)劑已經(jīng)有9人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有18人回復(fù)
314求調(diào)劑已經(jīng)有8人回復(fù)
08工科 320總分求調(diào)劑已經(jīng)有5人回復(fù)
085601專碩，總分342求調(diào)劑，地區(qū)不限已經(jīng)有5人回復(fù)
311求調(diào)劑已經(jīng)有6人回復(fù)
化工學(xué)碩306求調(diào)劑已經(jīng)有10人回復(fù)
286求調(diào)劑已經(jīng)有8人回復(fù)

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

有誰用過 phonon 這個(gè)軟件，希望交流下使用經(jīng)驗(yàn) 已經(jīng)有13人回復(fù)
聲子譜計(jì)算在Gamma點(diǎn)出0.2THz虛頻，可以忽略嗎？已經(jīng)有52人回復(fù)
石墨烯聲子譜Gamma點(diǎn)附近虛頻，怎么處理已經(jīng)有5人回復(fù)
用VASP-DFPT+phonopy計(jì)算石墨烯聲子譜求助已經(jīng)有37人回復(fù)
聲子譜��！謝謝��！已經(jīng)有4人回復(fù)
Gamma點(diǎn)處出現(xiàn)比較大的虛頻說明什么已經(jīng)有3人回復(fù)
你們?cè)趺刺幚鞧amma的虛頻的？已經(jīng)有4人回復(fù)
聲子譜虛頻求助已經(jīng)有10人回復(fù)
degauss的收斂已經(jīng)有13人回復(fù)
哭了，各種虛頻！已經(jīng)有23人回復(fù)
用自己生成的贗勢(shì)算能帶效果很好，算聲子卻一團(tuán)糟已經(jīng)有11人回復(fù)
關(guān)于用超胞方法計(jì)算聲子的一些經(jīng)驗(yàn) 已經(jīng)有21人回復(fù)
怎么消除虛頻，細(xì)節(jié) 已經(jīng)有10人回復(fù)
求助如何分析castep聲子計(jì)算的結(jié)果已經(jīng)有7人回復(fù)
自旋體系的聲子計(jì)算應(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小時(shí)
蟲號(hào): 3028351
注冊(cè): 2014-03-08
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

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

樓主是用什么軟件做的？我用Phonopy做聲子譜時(shí)，曾經(jīng)因?yàn)閷?duì)稱性，超胞的奇偶等很多原因出現(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小時(shí)
蟲號(hào): 3513863
注冊(cè): 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做聲子譜時(shí)，曾經(jīng)因?yàn)閷?duì)稱性，超胞的奇偶等很多原因出現(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小時(shí)
蟲號(hào): 1191315
注冊(cè): 2011-01-17
專業(yè): 理論和計(jì)算化學(xué)

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

引用回帖:

對(duì)稱性和超胞的奇偶？求樓主詳細(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小時(shí)
蟲號(hào): 1760414
注冊(cè): 2012-04-17
專業(yè): 物理學(xué)I

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

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

發(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小時(shí)
蟲號(hào): 2978560
注冊(cè): 2014-02-19
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

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

引用回帖:

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

贊一下

回復(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小時(shí)
蟲號(hào): 4437346
注冊(cè): 2016-02-26
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

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

贊一下

回復(fù)此樓

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

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

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

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 一志愿華中師范071000，325求調(diào)劑 +6	RuitingC 2026-03-12	6/300	2026-03-16 14:50 by 可淡不可忘
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz

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

[交流] gamma點(diǎn)聲子虛頻 已有5人參與

» 猜你喜歡

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

[交流] gamma點(diǎn)聲子虛頻已有5人參與

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