Differences

This shows you the differences between two versions of the page.

--- wien2k:scf計算 [2019/01/26 01:22]
koudai [計算の設定の初期化]
+++ wien2k:scf計算 [2019/01/26 01:47]
koudai [計算の設定の初期化]
@@ Line 153: / Line 153: @@
 最終的にほしい物理量がk点数に対して収束しているか、確認する必要があります。
-ここではk点数を1000にしてみましょう（ちなみにデフォルトも1000です）。
+ここではk点数を1000に、RKmaxを7.5にしてみましょう。
+（RKmaxは波動関数のカットオフを決めます。原子の種類によって必要な大きさが異なるので、選び方は http://susi.theochem.tuwien.ac.at/reg_user/faq/rkmax.html を参照してください）
 <code>
-$ init_lapw -b -numk 2000
+$ init_lapw -b -numk 1000 -rkmax 7.5
 </code>
-ここでオプション-bはバッチモードで、これをつけることで指定のもの以外はすべてまとめてデフォルトの値に設定されます。
+ここでオプション -b はバッチモードで、これをつけることで指定のもの以外はすべてまとめてデフォルトの値に設定されます。
-他の設定はオプション-hにより参照できます。
+他の設定はオプション -h により参照できます。
 <code>
@@ Line 187: / Line 188: @@
 実行後、いろいろなファイルができますが、計算の経過や重要な情報は TiC.scf で確認できます。
-今回は金属なのでフェルミエネルギーが出力されています。
-半導体であればエネルギーギャップがわかります。
+  $ grep ":ENE" TiC.scf     # 全エネルギーを確認します。
+  $ grep ":FER" TiC.scf     # フェルミ準位を確認します。単位は Ry です。
+  $ grep ":GAP" TiC.scf     # エネルギーギャップを確認します。金属の場合は0になります。
 これで電子状態が得られましたので、これをもとに状態密度やバンド分散などを計算することになります。
@@ Line 196: / Line 199: @@
 計算の収束を見るために、k点数を変更したいとします。
-再び init_lapw を走らせると、すべての設定が初期化されますので（今はすべてデフォルトなのでさほど問題ないですが）、k点数だけ変更するために x_lapw kgen を使います。
+再び init_lapw を走らせると全ての設定が初期化されますので（今は全てデフォルトなのでさほど問題ないですが）、k点数だけ変更するために x_lapw kgen を使います。
 <code>

雑多な記録

User Tools

Site Tools

Differences

Page Tools