quantumespresso:phonon:フォノンの計算
Differences
This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| quantumespresso:phonon:フォノンの計算 [2020/11/20 02:39] – [X点のフォノン] koudai | quantumespresso:phonon:フォノンの計算 [2021/06/27 22:04] (current) – external edit 127.0.0.1 | ||
|---|---|---|---|
| Line 57: | Line 57: | ||
| ^変数^初期値^説明^ | ^変数^初期値^説明^ | ||
| - | |tr2_ph|1.0d-12|フォノンの計算の収束条件。経験的に1.0d-14を使ったほうがよい。| | + | |tr2_ph|1.0d-12|フォノンの計算の収束条件。経験的に1.0d-14を使ったほうがよい| |
| - | |prefix|pwscf|SCF計算で使用したプレフィックスと同じにする。| | + | |prefix|pwscf|SCF計算で使用したプレフィックスと同じにする| |
| - | |epsil|.false.|誘電率を計算する。Γ点の計算に必要。ただし.true.にする場合は、半導体かつΓ点の計算でないとエラーが出る。| | + | |epsil|.false.|誘電率を計算する。Γ点の計算に必要。ただし.true.にする場合は、半導体かつΓ点の計算でないとエラーが出る| |
| - | |outdir|./ | + | |outdir|./ |
| - | |fildyn|matdyn|結果であるdynamical matrixを出力するファイル名。| | + | |fildyn|matdyn|結果であるdynamical matrixを出力するファイル名| |
| * 金属の場合は、epsil=.false.として、計算するk点を 0.01 0.0 0.0 のようにΓ点から少しずらしてやるとうまくいきます。 | * 金属の場合は、epsil=.false.として、計算するk点を 0.01 0.0 0.0 のようにΓ点から少しずらしてやるとうまくいきます。 | ||
| Line 120: | Line 120: | ||
| * Γ点のフォノンは q->0 の極限で定義されます。lpermは q=(q(1), q(2), q(3)) に沿って q->0 の極限をとることを指定します。 | * Γ点のフォノンは q->0 の極限で定義されます。lpermは q=(q(1), q(2), q(3)) に沿って q->0 の極限をとることを指定します。 | ||
| - | | + | * Si.phG.inで波数を 0.01 0.0 0.0 とした場合はasrおよびlpermの行を削除して q(1) = 0.01 としてください。 |
| + | | ||
| + | * zero-dim ... 分子の計算で使います | ||
| + | * one-dim ... 1次元系の計算(カーボンナノチューブなど)に使います | ||
| + | * crystal ... simpleよりも計算精度があがりますが、その分コストも増えます。simpleでうまく行かない場合はこちらにします。 | ||
| 実行は次のようにします。 | 実行は次のようにします。 | ||
quantumespresso/phonon/フォノンの計算.1605807573.txt.gz · Last modified: 2021/06/27 21:59 (external edit)
