Differences

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

--- wien2k:scf計算 [2020/07/09 13:09]
koudai [計算の設定の初期化]
+++ wien2k:scf計算 [2021/07/14 21:48]
koudai
@@ Line 202: / Line 202: @@
 実行後、いろいろなファイルができますが、計算の経過や重要な情報は TiC.scf で確認できます。
+<code>
-  $ grep ":ENE" TiC.scf     # 全エネルギー
+$ grep ":ENE" TiC.scf     # 全エネルギー
-  $ grep ":FER" TiC.scf     # フェルミ準位。単位は Ry です。
+（略）
-  $ grep ":GAP" TiC.scf     # エネルギーギャップ。金属の場合は0になります。
+:ENE  : ********** TOTAL ENERGY IN Ry =        -1783.96056300
+$ grep ":FER" TiC.scf     # フェルミ準位。単位は Ry です。
+（略）
+:FER  : F E R M I - ENERGY(TETRAH.M.)=   0.7492912035
+$ grep ":GAP" TiC.scf     # エネルギーギャップ。金属の場合は0になります。
+（略）
+:GAP (global)   :    0.0    Ry =     0.0   eV (metal)
+</code>
 これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。
@@ Line 212: / Line 219: @@
 ===== SCF計算の再実行 =====
-SCF計算を再度実行したい場合は、フラグ-NIをつけて run_lapw を実行します。
+SCF計算を再度実行したい場合は、ブロイデン法の履歴を削除してから再実行します。
+（電子密度の分布は前回の計算のものが引き継がれます）
 <code>
-$ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -NI
+$ rm *.broyd*
+$ run_lapw -cc 0.0001 -ec 0.00001 -i 100
 </code>
-もしrun_lapwのフラグに変更がある場合は、ブロイデン法を使ってSCF計算を収束させているので、その履歴を削除してから再実行します。
+フラグ-NIをつけて再実行することも可能ですが、計算条件を変更する場合は使えません（エラーが出ます）。
-（電子密度の分布は前回の計算のものが引き継がれます）
 <code>
-$ rm *.broyd*
+$ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -NI
-$ run_lapw -cc 0.0001 -ec 0.00001 -i 100
 </code>
@@ Line 242: / Line 249: @@
 </code>
-再び計算を実行するには、フラグ-NIをつけて run_lapw を実行します。
+計算の再実行は上で述べたようブロイデン法の履歴を削除してから行います。
-<code>
+最初は少ないk点数（10くらい）から初めて、徐々に増やしていきましょう。
-$ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -NI
-</code>
-最初は少ないk点数（それこそ10くらい）から初めて、徐々に増やしていきましょう。
 k点数に対して全エネルギーやフェルミ準位（半導体の場合はエネルギーギャップの大きさ）をプロットしてみて、収束しているかどうか必ず確認しましょう。
@@ Line 255: / Line 258: @@
-====== SCF計算結果の保存 ======
+====== 計算結果の保存 ======
 計算の入出力ファイルを残しておきたい場合は save_lapw を使用します。

雑多な記録

User Tools

Site Tools

Differences

Page Tools