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計算 [2020/03/30 15:47] koudai [k点数の変更] |
wien2k:scf計算 [2020/03/31 14:48] koudai [k点数の変更] |
||
---|---|---|---|
Line 202: | Line 202: | ||
実行後、いろいろなファイルができますが、計算の経過や重要な情報は TiC.scf で確認できます。 | 実行後、いろいろなファイルができますが、計算の経過や重要な情報は TiC.scf で確認できます。 | ||
- | $ grep ": | + | $ grep ": |
- | $ grep ": | + | $ grep ": |
- | $ grep ": | + | $ grep ": |
これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。 | これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。 | ||
Line 225: | Line 225: | ||
再び計算を実行するには、フラグ-NIをつけて run_lapw を実行します。 | 再び計算を実行するには、フラグ-NIをつけて run_lapw を実行します。 | ||
+ | < | ||
+ | $ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -NI | ||
+ | </ | ||
+ | あるいは計算履歴を削除して | ||
< | < | ||
- | $ run_lapw -cc 0.0001 -ec 0.00001 -NI | + | $ rm *.broyd* |
+ | $ run_lapw -cc 0.0001 -ec 0.00001 | ||
</ | </ | ||
+ | でも可能です。 | ||
全エネルギーやフェルミ準位(半導体の場合はエネルギーギャップの大きさ)が、k点数に対して収束しているか必ず確認しましょう。 | 全エネルギーやフェルミ準位(半導体の場合はエネルギーギャップの大きさ)が、k点数に対して収束しているか必ず確認しましょう。 | ||
Line 239: | Line 245: | ||
====== SCF計算結果の保存 ====== | ====== SCF計算結果の保存 ====== | ||
- | 構造最適化前のstructファイルや計算結果のファイルを残しておきたい場合は save_lapw を使用します。 | + | 計算の入出力ファイルを残しておきたい場合は save_lapw を使用します。 |
例えばk点を増やして計算する際に、k点数を1000で計算したときの入力ファイルや結果を k1000 という名前のディレクトリに保存しておきたい場合は次のようにします。 | 例えばk点を増やして計算する際に、k点数を1000で計算したときの入力ファイルや結果を k1000 という名前のディレクトリに保存しておきたい場合は次のようにします。 | ||