User Tools

Site Tools


wien2k:軌道ポテンシャル

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
wien2k:軌道ポテンシャル [2020/07/05 21:28]
koudai [計算方法]
wien2k:軌道ポテンシャル [2020/07/31 11:26]
koudai [計算方法]
Line 10: Line 10:
  
   - 事前にSCF計算を行います。収束は粗めで構いません。このときスピン分極を許した計算にします<code>   - 事前にSCF計算を行います。収束は粗めで構いません。このときスピン分極を許した計算にします<code>
-instgen -ask  # 全部 n (non-magnetic) にする+instgen_lapw -ask  # 全部 n (non-magnetic) にする
 $ init_lapw -b -numk 1000 -rkmax 8.5 -sp $ init_lapw -b -numk 1000 -rkmax 8.5 -sp
-$ runsp_c_lapw -ec 0.0001 -i 100 -p+$ runsp_c_lapw -i 100 -p
 </code> </code>
     * runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します     * runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します
-    * プログラム上は軌道ポテンシャルを入れない runsp_c_lapw 計算はスキップできますが、変な状態に収束してしまうことがあります。必ず軌道ポテンシャルなしの計算を最初にしてください+    * プログラム上は軌道ポテンシャルを入れない runsp_c_lapw 計算はスキップできますが、予期せぬ状態に収束してしまうことがあります。必ず軌道ポテンシャルなしの計算を最初にしてください
   - 軌道ポテンシャルの設定ファイルを作成します<code>   - 軌道ポテンシャルの設定ファイルを作成します<code>
 $ init_orb_lapw -orb $ init_orb_lapw -orb
Line 26: Line 26:
 Enter the name, l, U(Ry) and J of the Atoms (eg. Fe 2 0.3 0.0; exit with RETURN): Enter the name, l, U(Ry) and J of the Atoms (eg. Fe 2 0.3 0.0; exit with RETURN):
 </code> </code>
-  - 入力ファイルが表示されるので、正しく設定されたか確認します。割と間違って出力されることが多いで+  - 入力ファイルが表示されるので、正しく設定されたか確認します
     * NiO.indm     * NiO.indm
       * 2行目: +Uを入れる原子の数       * 2行目: +Uを入れる原子の数
Line 55: Line 55:
 </file> </file>
     * https://wien.zeus.theochem.tuwien.ac.narkive.com/dV0E1fCH/how-to-include-d-and-f-states-in-calculation     * https://wien.zeus.theochem.tuwien.ac.narkive.com/dV0E1fCH/how-to-include-d-and-f-states-in-calculation
-  - ブロイデン法の履歴を削除して、オプション-orbをつけて実行します<code>+  - ブロイデン法の履歴を削除して(削除しないままSCF計算を実行してしまうとファイルが壊れてしまうので注意)、オプション-orbをつけて実行します<code>
 $ rm *.broyd* $ rm *.broyd*
 $ runsp_c_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p $ runsp_c_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p
Line 63: Line 63:
 $ runsp_c_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p -NI $ runsp_c_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p -NI
 </code> </code>
-    * Uを変更するときは NiO.inorb を直接編集してください。 +    * 軌道ポテンシャルの大きさを変更するときは NiO.inorb を直接編集してください。 
-    * Uをいれる軌道は途中で変更できないので、変更したい場合は最初からやり直してください。+    * 軌道ポテンシャルが入っている状態でSCFを再計算するときは、ブロイデン法の履歴を削除しないでください。 
 +    * 軌道ポテンシャルをいれる軌道の種類は途中で変更できないので、変更したい場合は最初から(init_lapwから)やり直してください。
  
 スピン分極を許しているので x_lapw を使う場合は -up が必要です。 スピン分極を許しているので x_lapw を使う場合は -up が必要です。
-また x_lapw lapw1 を実行する際は、必ずオプション -orb をつけて実行してください。+また x_lapw lapw1 および x_lapw lapwso を実行する際は、必ずオプション -orb をつけて実行してください。
 <code> <code>
 $ x_lapw lapw1 -up -orb -p $ x_lapw lapw1 -up -orb -p
wien2k/軌道ポテンシャル.txt · Last modified: 2021/06/27 22:04 (external edit)