Differences

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

--- wien2k:scf計算 [2020/06/24 10:56]
koudai [k点数の変更]
+++ wien2k:scf計算 [2021/07/14 21:45]
koudai
@@ Line 165: / Line 165: @@
 </code>
   * ここでフラグ -b はバッチモードで、これをつけることで指定のもの以外はすべてまとめてデフォルトの値に設定されます。
-  * 最後に init_lapw finished ok と出たら正孔です。このメッセージが出てこなかった場合、どこかで問題が起こっているので修正します。
+  * 最後に init_lapw finished ok と出たらうまく行っています。このメッセージが出てこなかった場合、どこかで問題が起こっているので修正します。
 オプションの意味は次のとおりです
@@ 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>
 これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。
+===== SCF計算の再実行 =====
+SCF計算を再度実行したい場合は、ブロイデン法の履歴を削除してから再実行します。
+（電子密度の分布は前回の計算のものが引き継がれます）
+<code>
+$ rm *.broyd*
+$ run_lapw -cc 0.0001 -ec 0.00001 -i 100
+</code>
+フラグ-NIをつけて再実行することも可能ですが、計算条件を変更する場合は使えません。
+<code>
+$ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -NI
+</code>
@@ Line 229: / Line 253: @@
 $ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -NI
 </code>
-あるいはブロイデン法を使ってSCF計算を収束させているので、その計算履歴を削除して
-<code>
-$ rm *.broyd*
-$ run_lapw -cc 0.0001 -ec 0.00001 -i 100
-</code>
-でも可能です（電子密度の分布は前回の計算のものが引き継がれます）。
 最初は少ないk点数（それこそ10くらい）から初めて、徐々に増やしていきましょう。
@@ Line 245: / Line 262: @@
-====== SCF計算結果の保存 ======
+====== 計算結果の保存 ======
 計算の入出力ファイルを残しておきたい場合は save_lapw を使用します。

雑多な記録

User Tools

Site Tools

Differences

Page Tools