Differences

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

--- wien2k:軌道ポテンシャル [2020/07/03 12:08]
koudai [概要]
+++ wien2k:軌道ポテンシャル [2020/07/05 21:28]
koudai [計算方法]
@@ Line 5: / Line 5: @@
   * 特にd電子系やf電子系といった局在が強い系に対して有効です。
   * 組み合わせる手法によってLDA+U法、あるいはGGA+U法とも呼ばれます
+  * 非磁性体の場合を解説します。反強磁性体の場合は[[WIEN2k:反強磁性]]を参照してください
 ====== 計算方法  ======
-  * 事前に通常の方法でSCF計算を行います
+  - 事前にSCF計算を行います。収束は粗めで構いません。このときスピン分極を許した計算にします<code>
+$ instgen -ask  # 全部 n (non-magnetic) にする
+$ init_lapw -b -numk 1000 -rkmax 8.5 -sp
+$ runsp_c_lapw -ec 0.0001 -i 100 -p
+</code>
+    * runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します
+    * プログラム上は軌道ポテンシャルを入れない runsp_c_lapw 計算はスキップできますが、変な状態に収束してしまうことがあります。必ず軌道ポテンシャルなしの計算を最初にしてください
   - 軌道ポテンシャルの設定ファイルを作成します<code>
 $ init_orb_lapw -orb
@@ Line 46: / Line 51: @@
   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点数を変更する場合は次のようにします<code>
+$ x_lapw kgen
+$ runsp_c_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p -NI
 </code>
-    * k点数やUをかえて再度実行するときは、ブロイデン法の履歴を削除しなくてもオプション -NI をつけるだけで十分です。
     * Uを変更するときは NiO.inorb を直接編集してください。
+    * Uをいれる軌道は途中で変更できないので、変更したい場合は最初からやり直してください。
+スピン分極を許しているので x_lapw を使う場合は -up が必要です。
+また x_lapw lapw1 を実行する際は、必ずオプション -orb をつけて実行してください。
+<code>
+$ x_lapw lapw1 -up -orb -p
+</code>

雑多な記録

User Tools

Site Tools

Differences

Page Tools