| 5 | 1/1 | 返回列表 |
| 查看: 3405 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
justinzhouke新蟲 (初入文壇)
|
[交流]
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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +4 | Ymlll 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +5 | S.H1 2026-03-18 | 5/250 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動化所實習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|