User Tools

Site Tools


quantumespresso:phonon:フォノンの計算

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

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 4: Line 4:
 Siを例に、フォノンの簡単な計算を実行してみる。 事前にpw.xによる電子状態計算を行う。 Siを例に、フォノンの簡単な計算を実行してみる。 事前にpw.xによる電子状態計算を行う。
  
-<file - si.scf.in>+<file - Si.scf.in>
 &control &control
    calculation='scf',    calculation='scf',
Line 29: Line 29:
  
 <code> <code>
-$ pw.x < si.scf.in > si.scf.out+$ pw.x < Si.scf.in > Si.scf.out
 </code> </code>
  
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:
 </code> </code>
  
-  * 金属の計算でSi.phG.out |ddv_scf|^2 が発散してしまう場合は、SCF計算の際にdegauss=0.001など小さいものを使ってください。計算がうまく行けば、収束がうまくいく範囲でdegaussを大きくします。+  * 金属の計算で |ddv_scf|^2 が発散してしまう場合は、SCF計算の際にdegauss=0.001など小さいものを使ってください。計算がうまく行けば、収束がうまくいく範囲でdegaussを大きくします。
  
 結果はfildynで指定したもの(今の場合はSi.dynG)に出力されます。 結果はfildynで指定したもの(今の場合はSi.dynG)に出力されます。
quantumespresso/phonon/フォノンの計算.txt · Last modified: 2021/06/27 22:04 (external edit)