Differences

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

--- wien2k:scf計算 [2021/07/14 21:45]
koudai
+++ wien2k:scf計算 [2023/04/23 19:47] (current)
koudai [SCF計算の実行]
@@ Line 127: / Line 127: @@
 実際にWIEN2kを使った第一原理計算を行うときは、対話形式でstructファイルを作成するよりも、結晶構造データの統一フォーマットであるCIFファイルから生成することが多くなると思います。
-上記のサイトからCIFファイルをダウンロードし、ファイル名をTiC.cifに変更して作業ディレクトリに保存します。
+上記のCODのサイトからCIFファイルをダウンロードし、ファイル名をTiC.cifに変更して作業ディレクトリに保存します。
 用意ができたら、プログラムcif2structを使ってstructファイルに変換します。
@@ Line 209: / Line 209: @@
 （略）
 :FER  : F E R M I - ENERGY(TETRAH.M.)=   0.7492912035
-$ grep ":GAP" TiC.scf     # エネルギーギャップ。金属の場合は0になります。
+$ grep ":GAP" TiC.scf     # バンドギャップの大きさ。金属の場合は0になります。
 （略）
 :GAP (global)   :    0.0    Ry =     0.0   eV (metal)
@@ Line 227: / Line 227: @@
 </code>
-フラグ-NIをつけて再実行することも可能ですが、計算条件を変更する場合は使えません。
+フラグ-NIをつけて再実行することも可能ですが、計算条件を変更する場合は使えません（エラーが出ます）。
 <code>
@@ Line 236: / Line 236: @@
 ====== k点数の変更 ======
-計算の収束を見るために、k点数を変更したいとします。
+k点数が十分であったかを見るために、k点数を変更したいとします。
 再び init_lapw を走らせると全ての設定や計算結果が初期化されますので、k点数だけ変更するために x_lapw kgen を使います。
@@ Line 249: / Line 249: @@
 </code>
-再び計算を実行するには、フラグ-NIをつけて run_lapw を実行します。
+計算の再実行は上で述べたようブロイデン法の履歴を削除してから行います。
-<code>
+最初は少ないk点数（10くらい）から初めて、徐々に増やしていきましょう。
-$ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -NI
-</code>
-最初は少ないk点数（それこそ10くらい）から初めて、徐々に増やしていきましょう。
 k点数に対して全エネルギーやフェルミ準位（半導体の場合はエネルギーギャップの大きさ）をプロットしてみて、収束しているかどうか必ず確認しましょう。

雑多な記録

User Tools

Site Tools

Differences

Page Tools