Differences

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

--- wien2k:scf計算 [2019/01/26 01:34]
koudai [k点数の変更]
+++ wien2k:scf計算 [2019/02/12 20:47]
koudai [k点数の変更]
@@ Line 149: / Line 149: @@
 SCF計算に入る前に、計算の設定の初期化をおこないます。
-ほとんどの場合はデフォルトのままでよいですが、ここではk点数だけを設定しましょう。
+ここではk点数とRKmaxを設定しましょう。
-k点数を増やせば計算の精度はあがりますが、計算時間もかかるようになります。
+k点数やRKmaxを増やせば計算の精度はあがりますが、計算時間もかかるようになります。
-最終的にほしい物理量がk点数に対して収束しているか、確認する必要があります。
+最終的にほしい物理量がこれらのユーザー側が設定した量に対して収束しているか確認する必要があります。
-ここではk点数を1000にしてみましょう（ちなみにデフォルトも1000です）。
+なお、RKmaxは波動関数のカットオフを決めます。原子の種類によって必要な大きさが異なるので、選び方は http://susi.theochem.tuwien.ac.at/reg_user/faq/rkmax.html を参照してください。
+ここではk点数を1000に、RKmaxを7.5にしてみましょう。
 <code>
-$ init_lapw -b -numk 2000
+$ init_lapw -b -numk 1000 -rkmax 7.5
 </code>
@@ Line 173: / Line 175: @@
 <code>
-$ run_lapw -ec 0.0001 -cc 0.0001 -i 100
+$ run_lapw -ec 0.00001 -cc 0.0001 -i 100
 </code>
@@ Line 198: / Line 200: @@
 計算の収束を見るために、k点数を変更したいとします。
-再び init_lapw を走らせると全ての設定が初期化されますので（今は全てデフォルトなのでさほど問題ないですが）、k点数だけ変更するために x_lapw kgen を使います。
+再び init_lapw を走らせると全ての設定が初期化されますので、k点数だけ変更するために x_lapw kgen を使います。
 <code>
@@ Line 214: / Line 216: @@
 <code>
 $ rm *.broyd*
-$ run_lapw -ec 0.0001 -cc 0.0001
+$ run_lapw -ec 0.00001 -cc 0.0001
 </code>
 case.broydX を削除せずに run_lapw を実行した場合は、60秒後に自動で削除してSCF計算がスタートします（その間にジョブのキャンセルが可能です）。
+全エネルギーやフェルミ準位（半導体の場合はギャップ）が、k点数に対して収束しているか必ず確認しましょう。
 ====== ジョブの一時停止 ======
@@ Line 228: / Line 230: @@
 </code>
+====== プログラム名について ======
+プログラム名の最後に _lapw をつけていましたが、これを省略しても良いようになっています。
+例えば x_lapw は単に x だけで実行可能です。

雑多な記録

User Tools

Site Tools

Differences

Page Tools