版塊導(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ù)此樓

» 猜你喜歡

0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287 已經(jīng)有5人回復(fù)
070300化學(xué)319求調(diào)劑已經(jīng)有6人回復(fù)
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑已經(jīng)有6人回復(fù)
一志愿西南交大，求調(diào)劑已經(jīng)有3人回復(fù)
286求調(diào)劑已經(jīng)有6人回復(fù)
焦慮已經(jīng)有11人回復(fù)
331求調(diào)劑（0703有機(jī)化學(xué) 已經(jīng)有7人回復(fù)
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有5人回復(fù)
297求調(diào)劑已經(jīng)有5人回復(fù)
302求調(diào)劑已經(jīng)有10人回復(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ā)表

[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +6	Liwangman 2026-03-15	6/300	2026-03-18 13:21 by 盡舜堯1
[考研] 344求調(diào)劑 +5	knight344 2026-03-16	5/250	2026-03-18 12:39 by 學(xué)員HHvs9Y
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3	QGZDSYS 2026-03-13	5/250	2026-03-18 11:13 by QGZDSYS
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8	秋有木北 2026-03-14	8/400	2026-03-16 15:21 by 哦哦123
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-12	7/350	2026-03-13 23:24 by 邱gl
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5	瑪卡巴卡啊哈 2026-03-11	5/250	2026-03-13 21:40 by JourneyLucky
[考研] 工科，求調(diào)劑 +3	我887 2026-03-11	3/150	2026-03-13 21:39 by JourneyLucky
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 【考研調(diào)劑求收留】 +3	Ceciilia 2026-03-11	3/150	2026-03-13 20:18 by JourneyLucky
[考研] 求調(diào)劑 +7	18880831720 2026-03-11	7/350	2026-03-13 16:10 by JourneyLucky
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky

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

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

» 猜你喜歡

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

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

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