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計算 [2017/10/28 01:10] koudai [計算の設定の初期化] |
wien2k:scf計算 [2017/10/28 01:51] koudai [SCF計算の実行] |
||
---|---|---|---|
Line 163: | Line 163: | ||
< | < | ||
- | $ run_lapw -ec 0.0001 -cc 0.0001 | + | $ run_lapw -ec 0.0001 -cc 0.0001 |
</ | </ | ||
ここでオプション-ecはエネルギーの収束で、オプション-ccは電荷の収束です。 | ここでオプション-ecはエネルギーの収束で、オプション-ccは電荷の収束です。 | ||
デフォルトはエネルギーの収束しか見ませんが、電荷の収束も見ておいたほうが無難です。 | デフォルトはエネルギーの収束しか見ませんが、電荷の収束も見ておいたほうが無難です。 | ||
+ | また、オプション-iで最大何回までSCF計算を行うかを指定できます。 | ||
+ | 他のオプションはオプション-hで確認できます。 | ||
+ | |||
+ | < | ||
+ | $ run_lapw -h | ||
+ | </ | ||
実行後、いろいろファイルができますが、重要なのはTiC.scfです。 | 実行後、いろいろファイルができますが、重要なのはTiC.scfです。 | ||
Line 175: | Line 181: | ||
これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。 | これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。 | ||
+ | |||
+ | ====== k点数の変更 ====== | ||
+ | |||
+ | 計算の収束を見るために、k点数を変更したいとします。 | ||
+ | 再びプログラムinit_lapwを走らせると、すべての設定が初期化されますので(今はすべてデフォルトなのでさほど問題ないですが)、k点数だけ変更するためにプログラムx_lapwを使います。 | ||
+ | |||
+ | < | ||
+ | $ x_lapw kgen | ||
+ | NUMBER OF K-POINTS IN WHOLE CELL: (0 allows to specify 3 divisions of G) | ||
+ | 2000 | ||
+ | | ||
+ | 72 k-points generated, ndiv= 12 12 12 | ||
+ | KGEN ENDS | ||
+ | 0.0u 0.0s 0:05.99 0.6% 0+0k 0+304io 0pf+0w | ||
+ | </ | ||
+ | |||
+ | 再び計算を実行するには、TiC.broydXを削除してからrun_lapwを実行します。 | ||
+ | |||
+ | < | ||
+ | $ rm *.broyd* | ||
+ | $ run_lapw -ec 0.0001 -cc 0.000 | ||
+ | </ | ||
+ | |||
+ | なお、case.broydXを削除しなかった場合は、60秒後に自動で削除してSCF計算がスタートします(その間にジョブのキャンセルが可能です)。 | ||