This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
wien2k:scf計算 [2019/01/26 01:46] koudai [計算の設定の初期化] |
wien2k:scf計算 [2019/02/08 19:08] koudai [SCF計算の実行] |
||
---|---|---|---|
Line 149: | Line 149: | ||
SCF計算に入る前に、計算の設定の初期化をおこないます。 | SCF計算に入る前に、計算の設定の初期化をおこないます。 | ||
- | ほとんどの場合はデフォルトのままでよいですが、ここではk点数だけを設定しましょう。 | + | ここではk点数とRKmaxを設定しましょう。 |
- | k点数を増やせば計算の精度はあがりますが、計算時間もかかるようになります。 | + | k点数やRKmaxを増やせば計算の精度はあがりますが、計算時間もかかるようになります。 |
- | 最終的にほしい物理量がk点数に対して収束しているか、確認する必要があります。 | + | 最終的にほしい物理量がこれらのユーザー側が設定した量に対して収束しているか確認する必要があります。 |
- | ここではk点数を2000に、RKmaxを8にしてみましょう。 | + | なお、RKmaxは波動関数のカットオフを決めます。原子の種類によって必要な大きさが異なるので、選び方は http:// |
- | (RKmaxは波動関数のカットオフを決めます。原子の種類によって必要な大きさが異なるので、選び方は http:// | + | |
+ | ここではk点数を1000に、RKmaxを7.5にしてみましょう。 | ||
< | < | ||
- | $ init_lapw -b -numk 2000 -rkmax | + | $ init_lapw -b -numk 1000 -rkmax |
</ | </ | ||
Line 174: | Line 175: | ||
< | < | ||
- | $ run_lapw -ec 0.0001 -cc 0.0001 -i 100 | + | $ run_lapw -ec 0.00001 -cc 0.0001 -i 100 |
</ | </ | ||
Line 199: | Line 200: | ||
計算の収束を見るために、k点数を変更したいとします。 | 計算の収束を見るために、k点数を変更したいとします。 | ||
- | 再び init_lapw を走らせると全ての設定が初期化されますので(今は全てデフォルトなのでさほど問題ないですが)、k点数だけ変更するために x_lapw kgen を使います。 | + | 再び init_lapw を走らせると全ての設定が初期化されますので、k点数だけ変更するために x_lapw kgen を使います。 |
< | < | ||
Line 220: | Line 221: | ||
case.broydX を削除せずに run_lapw を実行した場合は、60秒後に自動で削除してSCF計算がスタートします(その間にジョブのキャンセルが可能です)。 | case.broydX を削除せずに run_lapw を実行した場合は、60秒後に自動で削除してSCF計算がスタートします(その間にジョブのキャンセルが可能です)。 | ||
+ | 全エネルギーやフェルミ準位(半導体の場合はギャップ)が、k点数に対して収束しているか必ず確認しましょう。 | ||
====== ジョブの一時停止 ====== | ====== ジョブの一時停止 ====== | ||
Line 229: | Line 230: | ||
</ | </ | ||
+ | |||
+ | ====== プログラム名について ====== | ||
+ | |||
+ | プログラム名の最後に _lapw をつけていましたが、これを省略しても良いようになっています。 | ||
+ | 例えば x_lapw は単に x だけで実行可能です。 |