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:23] koudai [電子状態計算] |
quantumespresso:phonon:フォノンの計算 [2020/11/18 20:32] koudai [電子状態計算] |
||
---|---|---|---|
Line 1: | Line 1: | ||
- | ===== 電子状態計算 ===== | + | ===== SCF計算 ===== |
- | Siを例に、フォノンの簡単な計算を実行してみる。 事前にpw.xによる電子状態計算を行う。 | + | Siを例に、フォノンの簡単な計算を実行してみる。 事前にpw.xによるSCF計算を行う。 |
<file - Si.scf.in> | <file - Si.scf.in> | ||
Line 56: | Line 56: | ||
^変数^初期値^説明^ | ^変数^初期値^説明^ | ||
|tr2_ph|1.0d-12|フォノンの計算の収束条件。経験的に1.0d-14を使ったほうがよい。| | |tr2_ph|1.0d-12|フォノンの計算の収束条件。経験的に1.0d-14を使ったほうがよい。| | ||
- | |prefix|pwscf|si.scf.inで使用したプレフィックスと同じにする。| | + | |prefix|pwscf|SCF計算で使用したプレフィックスと同じにする。| |
|epsil|.false.|誘電率を計算する。Γ点の計算に必要。ただし.true.にする場合は、半導体かつΓ点の計算でないとエラーが出る。| | |epsil|.false.|誘電率を計算する。Γ点の計算に必要。ただし.true.にする場合は、半導体かつΓ点の計算でないとエラーが出る。| | ||
|outdir|./ | |outdir|./ | ||
Line 69: | Line 69: | ||
</ | </ | ||
- | * 金属の計算でSi.phG.out | + | * 金属の計算で |ddv_scf|^2 が発散してしまう場合は、SCF計算の際にdegauss=0.001など小さいものを使ってください。計算がうまく行けば、収束がうまくいく範囲でdegaussを大きくします。 |
結果はfildynで指定したもの(今の場合はSi.dynG)に出力されます。 | 結果はfildynで指定したもの(今の場合はSi.dynG)に出力されます。 |