Differences

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

--- wien2k:軌道ポテンシャル [2020/07/03 11:58]
koudai [計算方法]
+++ wien2k:軌道ポテンシャル [2020/07/03 14:44]
koudai [計算方法]
@@ Line 3: / Line 3: @@
   * バンド計算では電子相関の効果が無視されるため、一般にバンドギャップを過小評価してしまいます
   * 相関の効果を取り入れることができ、計算コストがもっとも安いのが、軌道ポテンシャルを手で入れる方法です
+  * 特にd電子系やf電子系といった局在が強い系に対して有効です。
   * 組み合わせる手法によってLDA+U法、あるいはGGA+U法とも呼ばれます
@@ Line 8: / Line 9: @@
 ====== 計算方法  ======
-  * 事前に通常の方法でSCF計算を行います
+  - 事前に通常の方法でSCF計算を行います。収束は粗めで構いません。このときスピン分極を許した計算にします<code>
+$ init_lapw -b -numk 1000 -rkmax 8.5 -sp
+$ runsp_c_lapw -ec 0.0001 -i 100 -p
+</code>
+    * runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します
   - 軌道ポテンシャルの設定ファイルを作成します<code>
 $ init_orb_lapw -orb
@@ Line 45: / Line 49: @@
   1 2      index of atom, number of l, l
                nsic 0..AMF, 1..SIC, 2..HFM
-.88   0      U J (Ry)
+.44   0      U J (Ry)
-.88   0      U J (Ry)
+.44   0      U J (Ry)
 </file>
     * https://wien.zeus.theochem.tuwien.ac.narkive.com/dV0E1fCH/how-to-include-d-and-f-states-in-calculation
-  * ブロイデン法の履歴を削除して、オプション-orbをつけて実行します<code>
+  - ブロイデン法の履歴を削除して、オプション-orbをつけて実行します<code>
 $ rm *.broyd*
-$ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p
+$ runsp_c_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p
 </code>
-    * k点数やUをかえて再度実行するときは、ブロイデン法の履歴を削除しなくてもオプション -NI をつけるだけで十分です。
+    * 計算の収束の確認などでk点数を変更する場合は次のようにします<code>
-    * Uを変更するときは NiO_AFM.inorb を直接編集してください。
+$ x_lapw kgen
+$ runsp_c_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p -NI
+</code>
+    * Uを変更するときは NiO.inorb を直接編集してください。
+    * Uをいれる軌道は途中で変更できないので、変更したい場合は最初からやり直してください。
+状態密度をかくときなどで x_lapw lapw1 を実行する際は、必ずオプション -orb をつけて実行してください。

雑多な記録

User Tools

Site Tools

Differences

Page Tools