Differences

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

--- quantumespresso:半導体 [2020/07/23 18:48] – [K_POINTS] koudai
+++ quantumespresso:半導体 [2021/07/13 02:13] (current) – [K_POINTS] koudai
@@ Line 52: / Line 52: @@
 /
 &electrons
+  conv_thr = 1.0d-8
 /
 ATOMIC_SPECIES
@@ Line 84: / Line 85: @@
 ^変数^初期値^説明^
 |ibrav     |なし|結晶がどのブラベー格子に属するかを指定。2は面心立方格子である。|
-|celldm(1) |なし|ボーア半径を1としたときのブラベー格子の格子定数a。立方晶系の場合は一つだけ指定します。オングストロームで指定したければ代わりに変数Aを使います|
+|celldm(1) |なし|ボーア半径を1としたときのブラベー格子の格子定数a。オングストロームで指定したければ代わりに変数Aを使います|
 |nat       |なし|単位胞に含まれる原子の総数。Siはダイヤモンド構造をとるので2です。|
 |ntype     |なし|原子の種類の数。今の場合は単体のSiなので1です。|
@@ Line 91: / Line 92: @@
 ==== &electrons ====
-実際の計算にあたってのオプションです。 今はデフォルトのままで計算することにするので何も指定しません。
+実際の計算にあたってのオプションです。
+^変数^初期値^説明^
+|conv_thr|1.0d-6|SCF計算におけるエネルギーの収束判定。デフォルトだとややゆるいので厳しめにします。|
@@ Line 107: / Line 111: @@
 計算に使うk点を指定します。
-  * automaticオプションにより自動的に逆格子を6×6×6に分割したものを使います。 分割の数が多ければ多きほど計算が正確になりますが、それだけ計算に時間がかかるようになります。
+  * automaticオプションにより自動的に逆格子を6×6×6に分割したものを使います (Monkhorst-Pack grid)。分割の数が多ければ多きほど計算が正確になりますが、それだけ計算に時間がかかるようになります。
   * 後ろの3つの数字は、分割した点のk点を使うなら0、分割した点と点の中間点を使うなら1を指定します（つまり0を指定したらガンマ点が含まれますが、1を指定すると含まれません）。
+    * 半分だけずらした方が、少ないk点数でより精度が高くなることが知られています
 ===== 計算の実行 =====
@@ Line 117: / Line 121: @@
   $ pw.x < Si.scf.in > Si.scf.out
-実行したら Si.scf.out というファイルと、計算の途中でできたファイルがつまった silicon.save というディレクトリができます。
+実行したら、計算の途中経過が出力された Si.scf.out というファイルと、電荷密度などの計算結果のファイルが入った tmp/Si.save/ というディレクトリができます。
 計算に失敗したら crash という名前のファイルもできるので、それを開いて原因を確認します。
 計算が収束する様子は Si.scf.out の estimated scf accuracy を見ることで確認できます。
@@ Line 165: / Line 169: @@
 /
 &electrons
+  conv_thr = 1.0d-8
 /
 ATOMIC_SPECIES
@@ Line 232: / Line 237: @@
 /
 &electrons
+  conv_thr = 1.0d-8
 /
 ATOMIC_SPECIES
@@ Line 294: / Line 300: @@
   * Range ... どこまでのエネルギーの範囲をプロットするか、エネルギーの最高値と最低値の表示を参考に決めます。
   * output file (xmgr) ... Graceというグラフ描画ソフト用の出力ファイルの名前を指定します。 また、output file (ps)は、Postscript用のファイルです。
-  * Efermi ... フェルミ・エネルギーを指定しますが、今の場合は半導体なので、さっき見積もった価電子帯のトップのエネルギーを入力します。
+  * Efermi ... フェルミ・エネルギーを指定しますが、今の場合は半導体なので、さっき調べた価電子帯のトップのエネルギーを入力します。
   * deltaE ... グラフのエネルギーの刻み幅
   * reference E ... 刻みの基準となるエネルギーの値