| 5 | 1/1 | 返回列表 |
| 查看: 3404 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +7 | llllkkkhh 2026-03-18 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-18 | 3/150 |
|
|
[考研] 求材料調(diào)劑 +7 | 隔壁陳先生 2026-03-12 | 7/350 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|