This shows you the differences between two versions of the page.
Next revision | Previous revision Next revision Both sides next revision | ||
wien2k:tb-mbj交換ポテンシャル [2020/07/02 18:52] koudai 作成 |
wien2k:tb-mbj交換ポテンシャル [2020/07/02 19:35] koudai [計算方法] |
||
---|---|---|---|
Line 20: | Line 20: | ||
$ run_lapw -p -i 1 -NI | $ run_lapw -p -i 1 -NI | ||
</ | </ | ||
- | - ブロイデン法の履歴を削除し、再度mBJの初期設定を行う< | + | - ブロイデン法の履歴を削除し、mBJの設定を完了させる< |
$ rm *.broyd* | $ rm *.broyd* | ||
$ init_mbj_lapw | $ init_mbj_lapw | ||
Line 30: | Line 30: | ||
0.5 | 0.5 | ||
</ | </ | ||
- | * $\bar{g}$の定義や$A$, | ||
- | * 特に設定がなければcは元のTB-mBJの論文に出てくる値[[http:// | ||
- | - | ||
- | |||
- | - [Utils. >> init_mbj_lapw] 二回クリックしたらダメになるので注意 | ||
- | - [Execution >> run SCF] -> [Remove the files case.broyd[1|2]] | ||
- | - [Execution >> run SCF] It-numberを1にして[Start SCF cycle] | ||
- | - [Utils. >> init_mbj_lapw] -> [Remove the files case.broyd[1|2]] | ||
- | - [Utils. >> init_mbj_lapw]を再びクリック | ||
- | - [Files >> input_files]から// | ||
- | * 例えば$A=−0.012$, | ||
- | -0.012 | ||
- | 1.023 | ||
- | 0.5 | ||
- | </ | ||
* $\bar{g}$の定義や$A$, | * $\bar{g}$の定義や$A$, | ||
* 特に設定がなければcは元のTB-mBJの論文に出てくる値[[http:// | * 特に設定がなければcは元のTB-mBJの論文に出てくる値[[http:// | ||
Line 50: | Line 35: | ||
$ run_lapw -p -i 1000 -ec 0.00001 -cc 0.0001 | $ run_lapw -p -i 1000 -ec 0.00001 -cc 0.0001 | ||
</ | </ | ||
- | * なかなか収束しないので通常の10倍くらいにする。 | + | * なかなか収束しないので、SCF計算の回数は通常の10倍くらいにする。 |
これでTB-mBJポテンシャルを使った計算は終わり。 | これでTB-mBJポテンシャルを使った計算は終わり。 | ||
あとは状態密度やバンドを描いて、電子状態の変化を確認する。 | あとは状態密度やバンドを描いて、電子状態の変化を確認する。 | ||
+ | x_lapwの実行時に特になにかフラグをつけなくても、自動的にmBJポテンシャルを考慮した計算になる。 | ||
適切なTB-mBJポテンシャルの選び方は、上で書いた論文の値をそのまま使うのも一つであるが、光学伝導度測定などの実験で得られたバンドギャップを再現する大きさに合わせるのが手っ取り早い。 | 適切なTB-mBJポテンシャルの選び方は、上で書いた論文の値をそのまま使うのも一つであるが、光学伝導度測定などの実験で得られたバンドギャップを再現する大きさに合わせるのが手っ取り早い。 |