wien2k:軌道ポテンシャル
Differences
This shows you the differences between two versions of the page.
| Next revision | Previous revision Next revision Both sides next revision | ||
|
wien2k:軌道ポテンシャル [2020/07/03 11:55] koudai 作成 |
wien2k:軌道ポテンシャル [2020/07/03 15:51] koudai [計算方法] |
||
|---|---|---|---|
| Line 3: | Line 3: | ||
| * バンド計算では電子相関の効果が無視されるため、一般にバンドギャップを過小評価してしまいます | * バンド計算では電子相関の効果が無視されるため、一般にバンドギャップを過小評価してしまいます | ||
| * 相関の効果を取り入れることができ、計算コストがもっとも安いのが、軌道ポテンシャルを手で入れる方法です | * 相関の効果を取り入れることができ、計算コストがもっとも安いのが、軌道ポテンシャルを手で入れる方法です | ||
| + | * 特にd電子系やf電子系といった局在が強い系に対して有効です。 | ||
| * 組み合わせる手法によってLDA+U法、あるいはGGA+U法とも呼ばれます | * 組み合わせる手法によってLDA+U法、あるいはGGA+U法とも呼ばれます | ||
| Line 8: | Line 9: | ||
| ====== 計算方法 | ====== 計算方法 | ||
| - | | + | |
| + | $ instgen -ask # 全部 n (non-magnetic) にする | ||
| + | $ init_lapw -b -numk 1000 -rkmax 8.5 -sp | ||
| + | $ runsp_c_lapw -ec 0.0001 -i 100 -p | ||
| + | </ | ||
| + | * runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します | ||
| - 軌道ポテンシャルの設定ファイルを作成します< | - 軌道ポテンシャルの設定ファイルを作成します< | ||
| $ init_orb_lapw -orb | $ init_orb_lapw -orb | ||
| Line 23: | Line 28: | ||
| * NiO.indm | * NiO.indm | ||
| * 2行目: +Uを入れる原子の数 | * 2行目: +Uを入れる原子の数 | ||
| - | * 3,4行: +Uを入れる原子のインデックス（NiO_AFM.structで確認。Ni1の場合は1、Ni2の場合は2になっているはず）、+Uをいれる軌道の数（今回はd軌道の1つ）、+Uを入れる軌道の種類（d軌道ならl=2） | + | * 3,4行: +Uを入れる原子のインデックス（NiO.structで確認。Ni1の場合は1、Ni2の場合は2になっているはず）、+Uをいれる軌道の数（今回はd軌道の1つ）、+Uを入れる軌道の種類（d軌道ならl=2） |
| * 5行目: LDA+Uは 0 0 で固定。詳しくはUser' | * 5行目: LDA+Uは 0 0 で固定。詳しくはUser' | ||
| -12.0 Emin cutoff | -12.0 Emin cutoff | ||
| Line 31: | Line 36: | ||
| 0 0 r-index, | 0 0 r-index, | ||
| </ | </ | ||
| - | * NiO_AFM.inorb | + | * NiO.inorb |
| * 1行目: nmodはLDA+U計算なら1、natorbは+Uする原子の数、 iprは1にすれば出力がより詳細になる | * 1行目: nmodはLDA+U計算なら1、natorbは+Uする原子の数、 iprは1にすれば出力がより詳細になる | ||
| * 2行目: 変更不可 | * 2行目: 変更不可 | ||
| Line 45: | Line 50: | ||
| 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) |
| </ | </ | ||
| * https:// | * https:// | ||
| - | | + | |
| $ rm *.broyd* | $ rm *.broyd* | ||
| - | $ run_lapw | + | $ runsp_c_lapw |
| </ | </ | ||
| - | * k点数やUをかえて再度実行するときは、ブロイデン法の履歴を削除しなくてもオプション | + | * 計算の収束の確認などでk点数を変更する場合は次のようにします< |
| - | * Uを変更するときは | + | $ x_lapw kgen |
| + | $ runsp_c_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p -NI | ||
| + | </ | ||
| + | * Uを変更するときは | ||
| + | * Uをいれる軌道は途中で変更できないので、変更したい場合は最初からやり直してください。 | ||
| + | スピン分極を許しているので x_lapw を使う場合は -up が必要です。 | ||
| + | また x_lapw lapw1 を実行する際は、必ずオプション -orb をつけて実行してください。 | ||
| + | < | ||
| + | $ x_lapw lapw1 -up -orb -p | ||
| + | </ | ||
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
