wien2k:軌道ポテンシャル
Differences
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/06 14:20] koudai [計算方法] |
||
|---|---|---|---|
| Line 5: | Line 5: | ||
| * 特にd電子系やf電子系といった局在が強い系に対して有効です。 | * 特にd電子系やf電子系といった局在が強い系に対して有効です。 | ||
| * 組み合わせる手法によってLDA+U法、あるいはGGA+U法とも呼ばれます | * 組み合わせる手法によってLDA+U法、あるいはGGA+U法とも呼ばれます | ||
| + | * 非磁性体の場合を解説します。反強磁性体の場合は[[WIEN2k: | ||
| ====== 計算方法 | ====== 計算方法 | ||
| - | - 事前に通常の方法でSCF計算を行います。収束は粗めで構いません。このときスピン分極を許した計算にします< | + | - 事前にSCF計算を行います。収束は粗めで構いません。このときスピン分極を許した計算にします< |
| - | $ instgen | + | $ instgen_lapw |
| $ 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 | ||
| Line 62: | 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 | ||
| </ | </ | ||
| - | * Uを変更するときは NiO.inorb を直接編集してください。 | + | * 軌道ポテンシャルを変更するときは NiO.inorb を直接編集してください。 |
| - | * Uをいれる軌道は途中で変更できないので、変更したい場合は最初からやり直してください。 | + | * 軌道ポテンシャルをいれる軌道の種類は途中で変更できないので、変更したい場合は最初から（init_lapwから）やり直してください。 |
| スピン分極を許しているので x_lapw を使う場合は -up が必要です。 | スピン分極を許しているので x_lapw を使う場合は -up が必要です。 | ||
wien2k/軌道ポテンシャル.txt · Last modified: 2021/06/27 22:04 (external edit)
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
