This shows you the differences between two versions of the page.
wien2k:tb-mbj交換ポテンシャル [2020/07/02 19:28] koudai [計算方法] |
wien2k:tb-mbj交換ポテンシャル [2021/06/27 22:04] |
||
---|---|---|---|
Line 1: | Line 1: | ||
- | ====== 概要 | ||
- | |||
- | バンド計算では一般にバンドギャップを過小評価してしまう、あるいは半導体なのに金属と見積もってしまうという問題がある。 | ||
- | これを回避するためにいくつかの手法が考案されているが(例えばLDA+U法やハイブリッド汎関数法など)、現在もっとも精度よく値が出るのがGW近似と呼ばれるものである。 | ||
- | しかし、この手法はかなりの計算時間が必要で、極めて単純な構造の物質を除いてあまり実用的ではない。 | ||
- | そこで、GW近似の結果をお手軽に再現できるのが、手でポテンシャルを入れてバンドギャップを適正な値にする[[http:// | ||
- | WIEN2kではTB-mBJポテンシャルを公式にサポートしているので、簡単に計算をすることができる。 | ||
- | |||
- | |||
- | ====== 計算方法 | ||
- | |||
- | * 事前に通常の方法でSCF計算を行う。 | ||
- | * 磁気秩序を入れた系やスピン軌道相互作用のある系でもできるが、先にそちらをとりいれたSCF計算を済ませておくこと。 | ||
- | * LDA+Uと重複した使用もできるが、電子相関の効果を複数回取り入れることになり物理的によくわからない結果となってしまうので避けること。 | ||
- | |||
- | - mBJの初期設定を行う< | ||
- | $ init_mbj_lapw | ||
- | </ | ||
- | - 一度だけSCF計算を再実行する< | ||
- | $ run_lapw -p -i 1 -NI | ||
- | </ | ||
- | - ブロイデン法の履歴を削除し、再度mBJの設定を完了させる< | ||
- | $ rm *.broyd* | ||
- | $ init_mbj_lapw | ||
- | </ | ||
- | - case.in0abpにTB-mBJポテンシャル$c=A+B\bar{g}^e$を入力 | ||
- | * 例えば$A=−0.012$, | ||
- | -0.012 | ||
- | 1.023 | ||
- | 0.5 | ||
- | </ | ||
- | * $\bar{g}$の定義や$A$, | ||
- | * 特に設定がなければcは元のTB-mBJの論文に出てくる値[[http:// | ||
- | - SCF計算を実行する< | ||
- | $ run_lapw -p -i 1000 -ec 0.00001 -cc 0.0001 | ||
- | </ | ||
- | * なかなか収束しないので、SCF計算の回数は通常の10倍くらいにする。 | ||
- | |||
- | これでTB-mBJポテンシャルを使った計算は終わり。 | ||
- | あとは状態密度やバンドを描いて、電子状態の変化を確認する。 | ||
- | |||
- | 適切なTB-mBJポテンシャルの選び方は、上で書いた論文の値をそのまま使うのも一つであるが、光学伝導度測定などの実験で得られたバンドギャップを再現する大きさに合わせるのが手っ取り早い。 | ||
- | < | ||
- | $ grep ": | ||
- | </ | ||
- | |||
- | |||
- | ===== 収束しない場合 | ||
- | |||
- | - case.inmを以下のように書き換える。< | ||
- | PRATT | ||
- | 0.10 mixing FACTOR for BROYD/PRATT scheme | ||
- | 1.00 1.00 PW and CLM-scaling factors | ||
- | 9999 8 idum, HISTORY | ||
- | </ | ||
- | - SCF計算を実行する< | ||
- | $ run_lapw -p -i 1000 -ec 0.00001 -cc 0.0001 -NI | ||
- | </ | ||
- | - case.inmで書き換えたものをもとに戻す。< | ||
- | MSR1 | ||
- | 0.20 mixing FACTOR for BROYD/PRATT scheme | ||
- | 1.00 1.00 PW and CLM-scaling factors | ||
- | 9999 8 idum, HISTORY | ||
- | </ | ||
- | - SCF計算を再実行する< | ||
- | $ run_lapw -p -i 1000 -ec 0.00001 -cc 0.0001 -NI | ||
- | </ | ||
- | |||
- | もしこれでも収束しない場合、mixing factorを0.1から徐々に上げて0.2に近づけていくとよい。 | ||
- | |||
- | 参考:[[http:// | ||
- | |||
- | |||
- | |||
- | ====== 参考 TB-mBJポテンシャルcの値 | ||
- | |||
- | ===== オリジナル論文 ===== | ||
- | |||
- | < | ||
- | -0.012 | ||
- | 1.023 | ||
- | 0.5 | ||
- | </ | ||
- | [[http:// | ||
- | |||
- | ===== 絶縁体 ===== | ||
- | < | ||
- | 0.488 | ||
- | 0.500 | ||
- | 1 | ||
- | </ | ||
- | [[http:// | ||
- | |||
- | 実験で観測されているバンドギャップが7eV以上のものを絶縁体としています。 | ||
- | |||
- | ===== 半導体 ===== | ||
- | < | ||
- | 0.267 | ||
- | 0.656 | ||
- | 1 | ||
- | </ | ||
- | [[http:// | ||
- | |||
- | |||