This shows you the differences between two versions of the page.
Both sides previous revision Previous revision | Next revision Both sides next revision | ||
quantumespresso:phonon:フォノンの計算 [2020/11/18 20:16] koudai [X点のフォノン] |
quantumespresso:phonon:フォノンの計算 [2020/11/18 20:21] koudai [Γ点のフォノン] |
||
---|---|---|---|
Line 57: | Line 57: | ||
|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)に出力されます。 |