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 19:50] koudai [Γ点のフォノン] |
quantumespresso:phonon:フォノンの計算 [2020/11/18 20:06] koudai [電子状態計算] |
||
---|---|---|---|
Line 7: | Line 7: | ||
& | & | ||
| | ||
- | | + | |
- | | + | |
- | | + | |
/ | / | ||
&system | &system | ||
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点のフォノン ===== |