Differences

This shows you the differences between two versions of the page.

--- wien2k:tb-mbj交換ポテンシャル [2020/07/02 19:29] – [計算方法] koudai
+++ wien2k:tb-mbj交換ポテンシャル [2024/11/19 22:22] (current) – [計算方法] koudai
@@ Line 3: / Line 3: @@
 バンド計算では一般にバンドギャップを過小評価してしまう、あるいは半導体なのに金属と見積もってしまうという問題がある。
 これを回避するためにいくつかの手法が考案されているが（例えばLDA+U法やハイブリッド汎関数法など）、現在もっとも精度よく値が出るのがGW近似と呼ばれるものである。
-しかし、この手法はかなりの計算時間が必要で、極めて単純な構造の物質を除いてあまり実用的ではない。
+しかし、この手法は多くの計算コストが必要で、単純な構造の物質を除いてあまり実用的ではない。
 そこで、GW近似の結果をお手軽に再現できるのが、手でポテンシャルを入れてバンドギャップを適正な値にする[[http://dx.doi.org/10.1103/PhysRevLett.102.226401|Tran-Blaha Modified Becke-Johnson (TB-mBJ)交換ポテンシャル]]を使った近似手法である。
-WIEN2kではTB-mBJポテンシャルを公式にサポートしているので、簡単に計算をすることができる。
@@ Line 17: / Line 16: @@
 $ init_mbj_lapw
 </code>
-  - 一度だけSCF計算を再実行する<code>
+  - 一度だけSCF計算を再実行し、ブロイデン法の履歴を削除する<code>
 $ run_lapw -p -i 1 -NI
-</code>
-  - ブロイデン法の履歴を削除し、mBJの設定を完了させる<code>
 $ rm *.broyd*
+</code>
+  - mBJの設定を完了させる<code>
 $ init_mbj_lapw
 </code>
@@ Line 39: / Line 38: @@
 これでTB-mBJポテンシャルを使った計算は終わり。
 あとは状態密度やバンドを描いて、電子状態の変化を確認する。
+x_lapwの実行時に特になにかフラグをつけなくても、自動的にmBJポテンシャルを考慮した計算になる。
-適切なTB-mBJポテンシャルの選び方は、上で書いた論文の値をそのまま使うのも一つであるが、光学伝導度測定などの実験で得られたバンドギャップを再現する大きさに合わせるのが手っ取り早い。
+適切なTB-mBJポテンシャルの選び方は、上で書いた論文の値をそのまま使うのも一つであるが、光学測定などの実験で得られたバンドギャップを再現する大きさに合わせるのが手っ取り早い。
 <code>
 $ grep ":GAP" case.scf
@@ Line 48: / Line 48: @@
 ===== 収束しない場合  =====
-  - case.inmを以下のように書き換える。<file - case.inm>
+  - case.inm を開いて MSR1 -> PRATT にして、mixing factorを0.1にする<file - case.inm>
 PRATT   0.0   YES  (BROYD/PRATT, extra charge (+1 for additional e), norm)
 .10            mixing FACTOR for BROYD/PRATT scheme