Differences

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

--- wien2k:軌道ポテンシャル [2020/07/03 11:57]
koudai [計算方法]
+++ wien2k:軌道ポテンシャル [2020/07/03 14:39]
koudai [計算方法]
@@ Line 3: / Line 3: @@
   * バンド計算では電子相関の効果が無視されるため、一般にバンドギャップを過小評価してしまいます
   * 相関の効果を取り入れることができ、計算コストがもっとも安いのが、軌道ポテンシャルを手で入れる方法です
+  * 特にd電子系やf電子系といった局在が強い系に対して有効です。
   * 組み合わせる手法によってLDA+U法、あるいはGGA+U法とも呼ばれます
@@ Line 8: / Line 9: @@
 ====== 計算方法  ======
-  * 事前に通常の方法でSCF計算を行います
+  - 事前に通常の方法でSCF計算を行います。k点数や収束は粗めで構いません。このときスピン分極を許した計算にします<code>
+$ init_lapw -b -numk 100 -rkmax 8.5 -sp
+$ runsp_c_lapw -ec 0.0001 -i 100 -p
+</code>
+    * runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します
   - 軌道ポテンシャルの設定ファイルを作成します<code>
 $ init_orb_lapw -orb
@@ Line 23: / Line 27: @@
     * NiO.indm
       * 2行目: +Uを入れる原子の数
-      * 3,4行: +Uを入れる原子のインデックス（NiO_AFM.structで確認。Ni1の場合は1、Ni2の場合は2になっているはず）、+Uをいれる軌道の数（今回はd軌道の1つ）、+Uを入れる軌道の種類（d軌道ならl=2）
+      * 3,4行: +Uを入れる原子のインデックス（NiO.structで確認。Ni1の場合は1、Ni2の場合は2になっているはず）、+Uをいれる軌道の数（今回はd軌道の1つ）、+Uを入れる軌道の種類（d軌道ならl=2）
       * 5行目: LDA+Uは 0 0 で固定。詳しくはUser's Guideの7.10節<file - NiO_AFM.indm>
 -12.0         Emin cutoff
@@ 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>
 $ 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を変更するときは NiO_AFM.inorb を直接編集してください。
+    * Uをいれる軌道は途中で変更できないので、変更したい場合は最初からやり直してください。

雑多な記録

User Tools

Site Tools

Differences

Page Tools