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:31] koudai [SCF計算の実行] |
wien2k:scf計算 [2020/03/20 16:16] koudai [計算の設定の初期化] |
||
---|---|---|---|
Line 114: | Line 114: | ||
</ | </ | ||
- | init.structというファイルができます(この他にもいろいろできますが、今は無視)ので、適当な名前に変更しておきます。 | + | init.structというファイルを (プロジェクト名).struct |
< | < | ||
Line 149: | Line 149: | ||
SCF計算に入る前に、計算の設定の初期化をおこないます。 | SCF計算に入る前に、計算の設定の初期化をおこないます。 | ||
- | ほとんどの場合はデフォルトのままでよいですが、ここではk点数だけを設定しましょう。 | + | ここではk点数とRKmaxを設定しましょう。 |
- | k点数を増やせば計算の精度はあがりますが、計算時間もかかるようになります。 | + | k点数やRKmaxを増やせば計算の精度はあがりますが、計算時間もかかるようになります。 |
- | 最終的にほしい物理量がk点数に対して収束しているか、確認する必要があります。 | + | 最終的にほしい物理量がこれらのユーザー側が設定した量に対して収束しているか確認する必要があります。 |
- | ここではk点数を1000にしてみましょう(ちなみにデフォルトも1000です)。 | + | なお、RKmaxは波動関数のカットオフを決めます。原子の種類によって必要な大きさが異なるので、選び方は http:// |
+ | |||
+ | ここではk点数を1000に、RKmaxを7.5にしてみましょう。 | ||
< | < | ||
- | $ init_lapw -b -numk 2000 | + | $ init_lapw -b -numk 1000 -rkmax 7.5 |
</ | </ | ||
- | ここでオプション | + | ここでフラグ |
- | 他の設定はオプション -h により参照できます。 | + | |
+ | オプションの意味は次のとおりです | ||
+ | |||
+ | ^オプション^初期値^説明^ | ||
+ | |-rkmax|7.0|RKmax| | ||
+ | |-numk|1000|k点数| | ||
+ | |||
+ | |||
+ | 他の設定はフラグ-hにより参照できます。 | ||
< | < | ||
Line 173: | Line 183: | ||
< | < | ||
- | $ run_lapw -ec 0.0001 -cc 0.0001 -i 100 | + | $ run_lapw -cc 0.00001 -ec 0.00001 -i 100 |
</ | </ | ||
- | ここでオプション -ec はエネルギーの収束で、オプション -cc は電荷の収束です。 | + | ^オプション^初期値^説明^ |
- | デフォルトはエネルギーの収束しか見ませんが、電荷の収束も見ておいたほうが無難です。 | + | |-cc|なし|電荷の収束。単位はe。指定しなければ電荷の収束はチェックされないが、見ておいたほうが無難| |
- | また、オプション | + | |-ec|0.0001|エネルギーの収束。単位はRy。デフォルトだと少し荒い| |
+ | |-i|100|SCF計算を行うときの、繰り返し回数の最大値。| | ||
- | 他のオプションはオプション-hで確認できます。 | + | 他のオプションはフラグ-hで確認できます。 |
< | < | ||
Line 198: | Line 209: | ||
計算の収束を見るために、k点数を変更したいとします。 | 計算の収束を見るために、k点数を変更したいとします。 | ||
- | 再び init_lapw を走らせると、すべての設定が初期化されますので(今はすべてデフォルトなのでさほど問題ないですが)、k点数だけ変更するために x_lapw kgen を使います。 | + | 再び init_lapw を走らせると全ての設定が初期化されますので、k点数だけ変更するために x_lapw kgen を使います。 |
< | < | ||
Line 214: | Line 225: | ||
< | < | ||
$ rm *.broyd* | $ rm *.broyd* | ||
- | $ run_lapw -ec 0.0001 -cc 0.0001 | + | $ run_lapw -ec 0.00001 -cc 0.0001 |
</ | </ | ||
case.broydX を削除せずに run_lapw を実行した場合は、60秒後に自動で削除してSCF計算がスタートします(その間にジョブのキャンセルが可能です)。 | case.broydX を削除せずに run_lapw を実行した場合は、60秒後に自動で削除してSCF計算がスタートします(その間にジョブのキャンセルが可能です)。 | ||
+ | 全エネルギーやフェルミ準位(半導体の場合はギャップ)が、k点数に対して収束しているか必ず確認しましょう。 | ||
====== ジョブの一時停止 ====== | ====== ジョブの一時停止 ====== | ||
Line 228: | Line 239: | ||
</ | </ | ||
+ | |||
+ | ====== プログラム名について ====== | ||
+ | |||
+ | プログラム名の最後に _lapw をつけていましたが、これを省略しても良いようになっています。 | ||
+ | 例えば x_lapw は単に x だけで実行可能です。 |