Differences

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

--- wien2k:scf計算 [2019/02/12 20:47]
koudai [k点数の変更]
+++ wien2k:scf計算 [2020/03/28 18:50]
koudai
@@ Line 14: / Line 14: @@
-====== 結晶構造データの作成 ======
+====== SCF計算 ======
+===== 結晶構造データの作成 =====
 WIEN2kでは結晶構造はcase.structというファイルに保存されます。
@@ Line 36: / Line 39: @@
-===== makestruct_lapwを使う方法 =====
+==== makestruct_lapwを使う方法 ====
 基本的には質問に答えていくだけで完了します。
@@ Line 114: / Line 117: @@
 </code>
-init.structというファイルができます（この他にもいろいろできますが、今は無視）ので、適当な名前に変更しておきます。
+init.structというファイルを (プロジェクト名).struct に変更しておきます。
 <code>
@@ Line 121: / Line 124: @@
-===== CIFファイルからstructファイルを生成する方法 =====
+==== CIFファイルからstructファイルを生成する方法 ====
-実際にWIEN2kを使った第一原理計算を行うときは、対話形式でstructファイルを作成するよりもCIFファイルから生成することが多くなると思います。
+実際にWIEN2kを使った第一原理計算を行うときは、対話形式でstructファイルを作成するよりも、結晶構造データの統一フォーマットであるCIFファイルから生成することが多くなると思います。
 上記のサイトからCIFファイルをダウンロードし、ファイル名をTiC.cifに変更して作業ディレクトリに保存します。
@@ Line 146: / Line 149: @@
-====== 計算の設定の初期化 ======
+===== 計算の設定の初期化 =====
 SCF計算に入る前に、計算の設定の初期化をおこないます。
@@ Line 161: / Line 164: @@
 </code>
-ここでオプション -b はバッチモードで、これをつけることで指定のもの以外はすべてまとめてデフォルトの値に設定されます。
+ここでフラグ -b はバッチモードで、これをつけることで指定のもの以外はすべてまとめてデフォルトの値に設定されます。
-他の設定はオプション -h により参照できます。
+オプションの意味は次のとおりです
+^オプション^初期値^説明^
+|-rkmax|7.0|RKmax|
+|-numk|1000|k点数|
+他の設定はフラグ-hにより参照できます。
 <code>
@@ Line 169: / Line 180: @@
-====== SCF計算の実行 ======
+===== SCF計算の実行 =====
 これで準備が整いましたので計算を実行します。
@@ Line 175: / Line 186: @@
 <code>
-$ run_lapw -ec 0.00001 -cc 0.0001 -i 100
+$ run_lapw -cc 0.0001 -ec 0.00001 -i 100
 </code>
-ここでオプション -ec はエネルギーの収束で、オプション -cc は電荷の収束です。
+^オプション^初期値^説明^
-デフォルトはエネルギーの収束しか見ませんが、電荷の収束も見ておいたほうが無難です。
+|-cc|なし|電荷の収束。単位はe。指定しなければ電荷の収束はチェックされないが、見ておいたほうが無難|
-また、オプション -i で最大何回までSCF計算を行うかを指定できます。
+|-ec|0.0001|エネルギーの収束。単位はRy。デフォルトだと少し荒い|
+|-i|40|SCF計算を行うときの、繰り返し回数の最大値。40回だと足りないこともあるので少し多めに。|
-他のオプションはオプション-hで確認できます。
+他のオプションはフラグ-hで確認できます。
 <code>
@@ Line 212: / Line 224: @@
 </code>
-再び計算を実行するには TiC.broydX を削除してから run_lapw を実行します。
+再び計算を実行するには、フラグ-NIをつけて run_lapw を実行します。
 <code>
-$ rm *.broyd*
+$ run_lapw -cc 0.0001 -ec 0.00001 -NI
-$ run_lapw -ec 0.00001 -cc 0.0001
 </code>
-case.broydX を削除せずに run_lapw を実行した場合は、60秒後に自動で削除してSCF計算がスタートします（その間にジョブのキャンセルが可能です）。
 全エネルギーやフェルミ準位（半導体の場合はギャップ）が、k点数に対して収束しているか必ず確認しましょう。
+====== SCF計算結果の保存 ======
+構造最適化前のstructファイルや計算結果のファイルを残しておきたい場合は save_lapw を使用します。
+例えばk点を増やして計算する際に、k点数を1000で計算したときの入力ファイルや結果を k1000 という名前のディレクトリに保存しておきたい場合は次のようにします。
+<code>
+$ save_lapw -d k1000
+</code>
+また、この保存したファイルを復活させる場合は次のようにします。-fは同じ名前のファイルがあった時に、強制的に上書きするオプションです。
+<code>
+$ restore_lapw -d k1000 -f
+</code>
 ====== ジョブの一時停止 ======
@@ Line 230: / Line 258: @@
 </code>
+この空のファイルは、ジョブが止まった後に自動的に削除されます。
 ====== プログラム名について ======
 プログラム名の最後に _lapw をつけていましたが、これを省略しても良いようになっています。
 例えば x_lapw は単に x だけで実行可能です。

雑多な記録

User Tools

Site Tools

Differences

Page Tools