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/03 12:08]
koudai [概要]
wien2k:軌道ポテンシャル [2020/07/03 14:39]
koudai [計算方法]
Line 9: 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>   - 軌道ポテンシャルの設定ファイルを作成します<code>
 $ init_orb_lapw -orb $ init_orb_lapw -orb
Line 46: Line 49:
 2  1 2      index of atom, number of l, l 2  1 2      index of atom, number of l, l
   1               nsic 0..AMF, 1..SIC, 2..HFM   1               nsic 0..AMF, 1..SIC, 2..HFM
- 0.88   0      U J (Ry) + 0.44   0      U J (Ry) 
- 0.88   0      U J (Ry)+ 0.44   0      U J (Ry)
 </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>   * ブロイデン法の履歴を削除して、オプション-orbをつけて実行します<code>
 $ rm *.broyd* $ 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> </code>
-    * k点数やUをかえて再度実行するときは、ブロイデン法の履歴を削除しなくてもオプション -NI をつけるだけで十分です。 
     * Uを変更するときは NiO.inorb を直接編集してください。     * Uを変更するときは NiO.inorb を直接編集してください。
 +    * Uをいれる軌道は途中で変更できないので、変更したい場合は最初からやり直してください。
  
  
  
wien2k/軌道ポテンシャル.txt · Last modified: 2021/06/27 22:04 (external edit)