Differences

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

--- quantumespresso:phonon:フォノンの計算 [2020/11/18 20:12]
koudai [X点のフォノン]
+++ quantumespresso:phonon:フォノンの計算 [2020/11/18 20:24]
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 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>
+  * 金属の計算で |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>
@@ Line 142: / Line 152: @@
   $ ph.x < Si.phX.in > Si.phX.out
-結果は Si.dynX に出力される。
+結果は Si.dynX に出力されます。
 <file -Si.dynX>
@@ Line 168: / Line 178: @@
 </file>
-重に縮退したものが3つあることがわかる。
+重に縮退したものが3つあることがわかります。
 変位ベクトルを計算する場合は次のファイルを用意して、Γ点と同じように計算してください。

雑多な記録

User Tools

Site Tools

Differences

Page Tools