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計算 [2021/07/14 21:45] koudai |
wien2k:scf計算 [2023/04/23 19:45] koudai [k点数の変更] |
||
---|---|---|---|
Line 227: | Line 227: | ||
</ | </ | ||
- | フラグ-NIをつけて再実行することも可能ですが、計算条件を変更する場合は使えません。 | + | フラグ-NIをつけて再実行することも可能ですが、計算条件を変更する場合は使えません(エラーが出ます)。 |
< | < | ||
Line 236: | Line 236: | ||
====== k点数の変更 ====== | ====== k点数の変更 ====== | ||
- | 計算の収束を見るために、k点数を変更したいとします。 | + | k点数が十分であったかを見るために、k点数を変更したいとします。 |
再び init_lapw を走らせると全ての設定や計算結果が初期化されますので、k点数だけ変更するために x_lapw kgen を使います。 | 再び init_lapw を走らせると全ての設定や計算結果が初期化されますので、k点数だけ変更するために x_lapw kgen を使います。 | ||
Line 249: | Line 249: | ||
</ | </ | ||
- | 再び計算を実行するには、フラグ-NIをつけて run_lapw を実行します。 | + | 計算の再実行は上で述べたようブロイデン法の履歴を削除してから行います。 |
- | < | + | 最初は少ないk点数(10くらい)から初めて、徐々に増やしていきましょう。 |
- | $ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -NI | + | |
- | </ | + | |
- | + | ||
- | 最初は少ないk点数(それこそ10くらい)から初めて、徐々に増やしていきましょう。 | + | |
k点数に対して全エネルギーやフェルミ準位(半導体の場合はエネルギーギャップの大きさ)をプロットしてみて、収束しているかどうか必ず確認しましょう。 | k点数に対して全エネルギーやフェルミ準位(半導体の場合はエネルギーギャップの大きさ)をプロットしてみて、収束しているかどうか必ず確認しましょう。 | ||