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 00:42] koudai |
wien2k:scf計算 [2017/10/28 01:11] koudai [SCF計算の実行] |
||
---|---|---|---|
Line 41: | Line 41: | ||
< | < | ||
- | $ makestruct | + | $ makestruct_lapw |
************************************************ | ************************************************ | ||
Line 153: | Line 153: | ||
< | < | ||
- | $ init -h | + | $ init_lapw |
</ | </ | ||
Line 161: | Line 161: | ||
これで準備が整いましたので、計算を実行します。 | これで準備が整いましたので、計算を実行します。 | ||
プログラムrun_lapwを使用します。 | プログラムrun_lapwを使用します。 | ||
+ | |||
+ | < | ||
+ | $ run_lapw -ec 0.0001 -cc 0.0001 | ||
+ | </ | ||
+ | |||
+ | ここでオプション-ecはエネルギーの収束で、オプション-ccは電荷の収束です。 | ||
+ | デフォルトはエネルギーの収束しか見ませんが、電荷の収束も見ておいたほうが無難です。 | ||
+ | 他のオプションはオプション-hで確認できます。 | ||
+ | |||
+ | < | ||
+ | $ run_lapw -h | ||
+ | </ | ||
+ | |||
+ | 実行後、いろいろファイルができますが、重要なのはTiC.scfです。 | ||
+ | 今回は金属なのでフェルミエネルギーが出力されています。 | ||
+ | 半導体であればギャップサイズが出力されます。 | ||
+ | |||
+ | これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。 | ||
+ | |||