Differences

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

--- quantumespresso:phonon:フォノンの計算 [2020/11/18 20:06]
koudai [電子状態計算]
+++ quantumespresso:phonon:フォノンの計算 [2020/11/18 20:23]
koudai [電子状態計算]
@@ Line 4: / Line 4: @@
 Siを例に、フォノンの簡単な計算を実行してみる。 事前にpw.xによる電子状態計算を行う。
-<file - si.scf.in>
+<file - Si.scf.in>
 &control
    calculation='scf',
@@ Line 26: / Line 26: @@
 </file>
-  $ pw.x < si.scf.in > si.scf.out
+  * 金属の場合はsmearingを指定する必要がありますが、テトラヘドロン法はフォノン計算に使えません。
+<code>
+$ pw.x < Si.scf.in > Si.scf.out
+</code>
@@ Line 39: / Line 43: @@
   prefix='Si',
   epsil=.true.,
-  outdir='./',
+  outdir='./tmp/',
   fildyn='Si.dynG',
  /
@@ Line 57: / Line 61: @@
 |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
+<code>
+$ ph.x < Si.phG.in > Si.phG.out
+</code>
+  * 金属の計算でSi.phG.out |ddv_scf|^2 が発散してしまう場合は、SCF計算の際にdegauss=0.001など小さいものを使ってください。計算がうまく行けば、収束がうまくいく範囲でdegaussを大きくします。
 結果はfildynで指定したもの（今の場合はSi.dynG）に出力されます。
@@ Line 126: / Line 135: @@
 ===== X点のフォノン =====
-今度はX点(2π/a,0,0)のフォノンを調べる。 手順はさきほどと同じであるが、誘電率が計算できないことに注意する。
+今度はX点(2π/a,0,0)のフォノンを調べます。
+手順はさきほどと同じですが、誘電率の計算が不要です。
-<file - si.phX.in>
+<file - Si.phX.in>
 phonons of Si at X
  &inputph
   tr2_ph=1.0d-14,
-  prefix='si',
+  prefix='Si',
   amass(1)=28.08,
-  outdir='./',
+  outdir='./tmp/',
-  fildyn='si.dynX',
+  fildyn='Si.dynX',
  /
 .0 0.0 0.0
 </file>
-  $ 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: @@
 </file>
-重に縮退したものが3つあることがわかる。
+重に縮退したものが3つあることがわかります。
+変位ベクトルを計算する場合は次のファイルを用意して、Γ点と同じように計算してください。
+<file - Si.dynmat.in>
+&inputfil
+  dyn ='Si.dynX',
+  asr ='no'
+/
+</file>

雑多な記録

User Tools

Site Tools

Differences

Page Tools