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/03 15:55] 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 |