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/04/03 14:44] koudai [計算の設定の初期化] |
wien2k:scf計算 [2020/07/07 15:33] koudai [k点数の変更] |
||
---|---|---|---|
Line 137: | Line 137: | ||
これで入力ファイルTiC.structができます。 | これで入力ファイルTiC.structができます。 | ||
ただし、このままだとマフィンチン半径$R_{\rm MT}$が正しく設定されていないので計算します。 | ただし、このままだとマフィンチン半径$R_{\rm MT}$が正しく設定されていないので計算します。 | ||
+ | 拡張子の.structは不要です。 | ||
< | < | ||
Line 212: | Line 213: | ||
計算の収束を見るために、k点数を変更したいとします。 | 計算の収束を見るために、k点数を変更したいとします。 | ||
- | 再び init_lapw を走らせると全ての設定が初期化されますので、k点数だけ変更するために x_lapw kgen を使います。 | + | 再び init_lapw を走らせると全ての設定や計算結果が初期化されますので、k点数だけ変更するために x_lapw kgen を使います。 |
< | < | ||
Line 238: | Line 239: | ||
最初は少ないk点数(それこそ10くらい)から初めて、徐々に増やしていきましょう。 | 最初は少ないk点数(それこそ10くらい)から初めて、徐々に増やしていきましょう。 | ||
k点数に対して全エネルギーやフェルミ準位(半導体の場合はエネルギーギャップの大きさ)をプロットしてみて、収束しているかどうか必ず確認しましょう。 | k点数に対して全エネルギーやフェルミ準位(半導体の場合はエネルギーギャップの大きさ)をプロットしてみて、収束しているかどうか必ず確認しましょう。 | ||
+ | |||
+ | * run_lapw に新しいフラグをつけるときは、基本的に rm *.broyd* を使うと安全です | ||