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計算 [2019/01/26 01:22] koudai [計算の設定の初期化] |
wien2k:scf計算 [2019/01/26 01:47] koudai [計算の設定の初期化] |
||
---|---|---|---|
Line 153: | Line 153: | ||
最終的にほしい物理量がk点数に対して収束しているか、確認する必要があります。 | 最終的にほしい物理量がk点数に対して収束しているか、確認する必要があります。 | ||
- | ここではk点数を1000にしてみましょう(ちなみにデフォルトも1000です)。 | + | ここではk点数を2000に、RKmaxを7.5にしてみましょう。 |
+ | (RKmaxは波動関数のカットオフを決めます。原子の種類によって必要な大きさが異なるので、選び方は http:// | ||
< | < | ||
- | $ init_lapw -b -numk 2000 | + | $ init_lapw -b -numk 2000 -rkmax 7.5 |
</ | </ | ||
Line 187: | Line 188: | ||
実行後、いろいろなファイルができますが、計算の経過や重要な情報は TiC.scf で確認できます。 | 実行後、いろいろなファイルができますが、計算の経過や重要な情報は TiC.scf で確認できます。 | ||
- | 今回は金属なのでフェルミエネルギーが出力されています。 | + | |
- | 半導体であればエネルギーギャップがわかります。 | + | $ grep ": |
+ | $ grep ": | ||
+ | $ grep ": | ||
これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。 | これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。 | ||
Line 196: | Line 199: | ||
計算の収束を見るために、k点数を変更したいとします。 | 計算の収束を見るために、k点数を変更したいとします。 | ||
- | 再び init_lapw を走らせると、すべての設定が初期化されますので(今はすべてデフォルトなのでさほど問題ないですが)、k点数だけ変更するために x_lapw kgen を使います。 | + | 再び init_lapw を走らせると全ての設定が初期化されますので(今は全てデフォルトなのでさほど問題ないですが)、k点数だけ変更するために x_lapw kgen を使います。 |
< | < |