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:25] koudai [計算方法] |
wien2k:軌道ポテンシャル [2020/07/03 14:46] koudai [計算方法] |
||
---|---|---|---|
Line 9: | Line 9: | ||
====== 計算方法 | ====== 計算方法 | ||
- | - 事前に通常の方法でSCF計算を行います。このときスピン分極を許した計算にします< | + | - 事前に通常の方法でSCF計算を行います。収束は粗めで構いません。このときスピン分極を許した計算にします< |
$ init_lapw -b -numk 1000 -rkmax 8.5 -sp | $ init_lapw -b -numk 1000 -rkmax 8.5 -sp | ||
- | $ runsp_c_lapw | + | $ runsp_c_lapw -ec 0.0001 -i 100 -p |
</ | </ | ||
* runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します | * runsp_c_lapw はSCF計算において各原子の磁気モーメントがゼロになるように制限します | ||
Line 53: | Line 53: | ||
</ | </ | ||
* https:// | * https:// | ||
- | | + | |
$ rm *.broyd* | $ rm *.broyd* | ||
$ runsp_c_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p | $ runsp_c_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p | ||
</ | </ | ||
- | * k点数やUをかえて再度実行するときは、ブロイデン法の履歴を削除しなくてもオプション | + | * 計算の収束の確認などでk点数を変更する場合は次のようにします< |
+ | $ x_lapw kgen | ||
+ | $ runsp_c_lapw -cc 0.0001 -ec 0.00001 -i 100 -orb -p -NI | ||
+ | </ | ||
* Uを変更するときは NiO.inorb を直接編集してください。 | * Uを変更するときは NiO.inorb を直接編集してください。 | ||
+ | * Uをいれる軌道は途中で変更できないので、変更したい場合は最初からやり直してください。 | ||
+ | 状態密度をかくときなどで x_lapw lapw1 を実行する際は、必ずオプション -orb をつけて実行してください。 | ||
+ | < | ||
+ | $ $ x_lapw lapw1 -up -orb -p | ||
+ | </ | ||