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:12] koudai [X点のフォノン] |
quantumespresso:phonon:フォノンの計算 [2020/11/18 20:23] koudai [電子状態計算] |
||
---|---|---|---|
Line 4: | Line 4: | ||
Siを例に、フォノンの簡単な計算を実行してみる。 事前にpw.xによる電子状態計算を行う。 | Siを例に、フォノンの簡単な計算を実行してみる。 事前にpw.xによる電子状態計算を行う。 | ||
- | <file - si.scf.in> | + | <file - Si.scf.in> |
& | & | ||
| | ||
Line 26: | Line 26: | ||
</ | </ | ||
- | $ pw.x < si.scf.in > si.scf.out | + | |
+ | |||
+ | < | ||
+ | $ pw.x < Si.scf.in > Si.scf.out | ||
+ | </ | ||
Line 57: | Line 61: | ||
|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 のようにΓ点から少しずらしてやるとうまくいきます。 |
実行は次のようにします。 | 実行は次のようにします。 | ||
- | | + | < |
+ | $ ph.x < Si.phG.in > Si.phG.out | ||
+ | </ | ||
+ | |||
+ | * 金属の計算でSi.phG.out |ddv_scf|^2 が発散してしまう場合は、SCF計算の際にdegauss=0.001など小さいものを使ってください。計算がうまく行けば、収束がうまくいく範囲でdegaussを大きくします。 | ||
結果はfildynで指定したもの(今の場合はSi.dynG)に出力されます。 | 結果はfildynで指定したもの(今の場合はSi.dynG)に出力されます。 | ||
Line 126: | Line 135: | ||
===== X点のフォノン ===== | ===== X点のフォノン ===== | ||
- | 今度はX点(2π/ | + | 今度はX点(2π/ |
+ | 手順はさきほどと同じですが、誘電率の計算が不要です。 | ||
<file - Si.phX.in> | <file - Si.phX.in> | ||
Line 142: | Line 152: | ||
$ ph.x < Si.phX.in > Si.phX.out | $ ph.x < Si.phX.in > Si.phX.out | ||
- | 結果は Si.dynX に出力される。 | + | 結果は Si.dynX に出力されます。 |
<file -Si.dynX> | <file -Si.dynX> | ||
Line 168: | Line 178: | ||
</ | </ | ||
- | 2重に縮退したものが3つあることがわかる。 | + | 2重に縮退したものが3つあることがわかります。 |
変位ベクトルを計算する場合は次のファイルを用意して、Γ点と同じように計算してください。 | 変位ベクトルを計算する場合は次のファイルを用意して、Γ点と同じように計算してください。 |