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 14:47] koudai [計算方法] |
wien2k:軌道ポテンシャル [2020/07/05 21:28] 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) にする | ||
$ 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 -ec 0.0001 -i 100 -p | ||
</ | </ | ||
* runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します | * runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します | ||
+ | * プログラム上は軌道ポテンシャルを入れない runsp_c_lapw 計算はスキップできますが、変な状態に収束してしまうことがあります。必ず軌道ポテンシャルなしの計算を最初にしてください | ||
- 軌道ポテンシャルの設定ファイルを作成します< | - 軌道ポテンシャルの設定ファイルを作成します< | ||
$ init_orb_lapw -orb | $ init_orb_lapw -orb |