This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
wien2k:軌道ポテンシャル [2020/07/03 15:51] koudai [計算方法] |
wien2k:軌道ポテンシャル [2020/07/05 21:33] koudai [計算方法] |
||
---|---|---|---|
Line 5: | Line 5: | ||
* 特にd電子系やf電子系といった局在が強い系に対して有効です。 | * 特にd電子系やf電子系といった局在が強い系に対して有効です。 | ||
* 組み合わせる手法によってLDA+U法、あるいはGGA+U法とも呼ばれます | * 組み合わせる手法によってLDA+U法、あるいはGGA+U法とも呼ばれます | ||
+ | * 非磁性体の場合を解説します。反強磁性体の場合は[[WIEN2k: | ||
====== 計算方法 | ====== 計算方法 | ||
- | - 事前に通常の方法でSCF計算を行います。収束は粗めで構いません。このときスピン分極を許した計算にします< | + | - 事前にSCF計算を行います。収束は粗めで構いません。このときスピン分極を許した計算にします< |
$ instgen -ask # 全部 n (non-magnetic) にする | $ instgen -ask # 全部 n (non-magnetic) にする | ||
$ init_lapw -b -numk 1000 -rkmax 8.5 -sp | $ init_lapw -b -numk 1000 -rkmax 8.5 -sp | ||
Line 15: | Line 15: | ||
</ | </ | ||
* runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します | * runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します | ||
+ | * プログラム上は軌道ポテンシャルを入れない runsp_c_lapw 計算はスキップできますが、変な状態に収束してしまうことがあります。必ず軌道ポテンシャルなしの計算を最初にしてください | ||
- 軌道ポテンシャルの設定ファイルを作成します< | - 軌道ポテンシャルの設定ファイルを作成します< | ||
$ init_orb_lapw -orb | $ init_orb_lapw -orb | ||
Line 63: | Line 64: | ||
</ | </ | ||
* Uを変更するときは NiO.inorb を直接編集してください。 | * Uを変更するときは NiO.inorb を直接編集してください。 | ||
- | * Uをいれる軌道は途中で変更できないので、変更したい場合は最初からやり直してください。 | + | * Uをいれる軌道は途中で変更できないので、変更したい場合は最初から(init_lapwから)やり直してください。 |
スピン分極を許しているので x_lapw を使う場合は -up が必要です。 | スピン分極を許しているので x_lapw を使う場合は -up が必要です。 |