This shows you the differences between two versions of the page.
Next revision | Previous revision Next revision Both sides next revision | ||
quantumespresso:phonon:フォノンの計算 [2020/11/18 19:42] koudai 作成 |
quantumespresso:phonon:フォノンの計算 [2020/11/18 19:50] koudai [Γ点のフォノン] |
||
---|---|---|---|
Line 33: | Line 33: | ||
電子状態の計算が終われば、フォノンの振動数を計算できるようになる。 まず、k=(0, | 電子状態の計算が終われば、フォノンの振動数を計算できるようになる。 まず、k=(0, | ||
- | <file - si.phG.in> | + | <file - Si.phG.in> |
phonons of Si at Gamma | phonons of Si at Gamma | ||
& | & | ||
tr2_ph=1.0d-14, | tr2_ph=1.0d-14, | ||
- | prefix=' | + | prefix=' |
epsil=.true., | epsil=.true., | ||
outdir=' | outdir=' | ||
- | fildyn=' | + | fildyn=' |
/ | / | ||
0.0 0.0 0.0 | 0.0 0.0 0.0 | ||
</ | </ | ||
- | 第1行はコメント行とみなされ無視される。 また、最後の行は計算する波数(デカルト座標で、単位は格子定数aを使って2π/ | + | 第1行はコメント行とみなされ無視されます。 |
- | 各引数の意味は次のとおりである。 | + | また、最後の行は計算する波数(デカルト座標で、単位は格子定数aを使って2π/ |
- | 詳細はマニュアルを参照のこと。 | + | |
- | |tr2_ph|フォノンの計算の収束条件。デフォルトは1.0d-12であるが、経験的に1.0d-14を使ったほうがよい。| | + | 各引数の意味は次のとおりです。 |
- | |prefix|si.scf.inで使用したプレフィックスと同じにする。| | + | |
- | |epsil|デフォルトは.false.だが、誘電率を計算するならば.true.を指定する。ただし.true.にする場合は、半導体かつΓ点の計算でないとエラーが出る。| | + | |
- | |outdir|出力ファイルの場所。| | + | |
- | |fildyn|結果であるdynamical matrixを出力するファイル名。| | + | |
- | 実行は次のようにする。 | + | ^変数^初期値^説明^ |
+ | |tr2_ph|1.0d-12|フォノンの計算の収束条件。経験的に1.0d-14を使ったほうがよい。| | ||
+ | |prefix|pwscf|si.scf.inで使用したプレフィックスと同じにする。| | ||
+ | |epsil|.false.|誘電率を計算する。Γ点の計算に必要。ただし.true.にする場合は、半導体かつΓ点の計算でないとエラーが出る。| | ||
+ | |outdir|./ | ||
+ | |fildyn|matdyn|結果であるdynamical matrixを出力するファイル名。| | ||
- | $ ph.x < si.phG.in > si.phG.out | + | 金属の場合は、epsil=.false.として、計算するk点を 0.01 0.0 0.0 のようにΓ点から少しずらしてやるとうまくいきます。 |
+ | 実行は次のようにします。 | ||
- | 結果はsi.dynGで見ることができる。 | + | $ ph.x < Si.phG.in > Si.phG.out |
- | <file - si.dynG> | + | 結果はfildynで指定したもの(今の場合はSi.dynG)に出力されます。 |
+ | |||
+ | <file - Si.dynG> | ||
(中略) | (中略) | ||
| |