This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
quantumespresso:phonon:フォノンの計算 [2020/11/18 20:06] koudai [Γ点のフォノン] |
quantumespresso:phonon:フォノンの計算 [2020/11/18 20:16] koudai [X点のフォノン] |
||
---|---|---|---|
Line 126: | Line 126: | ||
===== X点のフォノン ===== | ===== X点のフォノン ===== | ||
- | 今度はX点(2π/ | + | 今度はX点(2π/ |
+ | 手順はさきほどと同じですが、誘電率の計算が不要です。 | ||
- | <file - si.phX.in> | + | <file - Si.phX.in> |
phonons of Si at X | phonons of Si at X | ||
& | & | ||
tr2_ph=1.0d-14, | tr2_ph=1.0d-14, | ||
- | prefix=' | + | prefix=' |
amass(1)=28.08, | amass(1)=28.08, | ||
- | outdir=' | + | outdir=' |
- | fildyn=' | + | fildyn=' |
/ | / | ||
1.0 0.0 0.0 | 1.0 0.0 0.0 | ||
</ | </ | ||
- | $ ph.x < si.phX.in > si.phX.out | + | $ ph.x < Si.phX.in > Si.phX.out |
- | 結果はsi.dynXに出力される。 | + | 結果は |
- | <file -si.dynX> | + | <file -Si.dynX> |
(中略) | (中略) | ||
| | ||
Line 168: | Line 169: | ||
</ | </ | ||
- | 2重に縮退したものが3つあることがわかる。 | + | 2重に縮退したものが3つあることがわかります。 |
+ | |||
+ | 変位ベクトルを計算する場合は次のファイルを用意して、Γ点と同じように計算してください。 | ||
+ | |||
+ | |||
+ | <file - Si.dynmat.in> | ||
+ | & | ||
+ | dyn =' | ||
+ | asr =' | ||
+ | / | ||
+ | </ |