quantumespresso:phonon:フォノンの計算
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision | Next revision Both sides next revision | ||
|
quantumespresso:phonon:フォノンの計算 [2020/11/18 19:50] koudai [Γ点のフォノン] |
quantumespresso:phonon:フォノンの計算 [2020/11/18 20:05] koudai [Γ点のフォノン] |
||
|---|---|---|---|
| Line 94: | Line 94: | ||
| * 振動モードの既約表現はSi.phG.out内に出力されている。 | * 振動モードの既約表現はSi.phG.out内に出力されている。 | ||
| + | |||
| + | フォノンの振動に対応する変位ベクトルはXcrysDenを使って見ることができます。 | ||
| + | 次のファイルを用意します。 | ||
| + | |||
| + | <file - Si.dynmat.in> | ||
| + | & | ||
| + | dyn =' | ||
| + | asr =' | ||
| + | lperm = .true., | ||
| + | q(1)=1.0, | ||
| + | q(2)=0.0, | ||
| + | q(3)=0.0 | ||
| + | / | ||
| + | </ | ||
| + | |||
| + | Γ点のフォノンは q->0 の極限で定義されます。lpermは q=(q(1), q(2), q(3)) に沿って q->0 の極限をとることを指定します。 | ||
| + | |||
| + | 実行は次のようにします。 | ||
| + | < | ||
| + | $ dynmat.x < Si.dynmat.in > Si.dynmat.out | ||
| + | </ | ||
| + | |||
| + | すると dynmat.axsf というファイルができるので、以下の手順で変位ベクトルを表示します。 | ||
| + | |||
| + | - XCrysDenを起動して、[File] -> [Open Structure] -> [Open AXSF] よりdynmat.axsfを開く | ||
| + | - ウィンドウが開くので、矢印ボタンを押して変位ベクトルを表示したいフォノンを選ぶ（エネルギーの低い順です） | ||
| + | - いったん [Hide] ボタンを押してウィンドウを最小化する | ||
| + | - [Display] から [Forces] を選ぶ | ||
| + | - このままだと矢印が長すぎるので、[Modify] -> [Force Settings] でLength Factorの大きさを調節する（デフォルトは200だが30くらいがちょうどよい） | ||
| ===== X点のフォノン ===== | ===== X点のフォノン ===== | ||
quantumespresso/phonon/フォノンの計算.txt · Last modified: 2021/06/27 22:04 (external edit)
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
