當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

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ù)此樓

» 猜你喜歡

能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有11人回復(fù)
求材料調(diào)劑已經(jīng)有9人回復(fù)
材料專碩英一數(shù)二306 已經(jīng)有5人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有9人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有12人回復(fù)
0703化學(xué) 305求調(diào)劑已經(jīng)有4人回復(fù)
0703化學(xué)調(diào)劑，求各位老師收留已經(jīng)有10人回復(fù)
271材料工程求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(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的回帖

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的回帖

查看全部 7 個(gè)回答

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的回帖

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

[考研] 085600材料與化工調(diào)劑 324分 +9	llllkkkhh 2026-03-18	9/450	2026-03-19 07:40 by BruceLiu320
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 0703化學(xué)求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9	柴郡貓_ 2026-03-12	9/450	2026-03-17 10:18 by Sammy2
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu