User Tools

Site Tools


wien2k:tb-mbj交換ポテンシャル

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
wien2k:tb-mbj交換ポテンシャル [2020/07/02 19:28]
koudai [計算方法]
wien2k:tb-mbj交換ポテンシャル [2021/06/27 22:04] (current)
Line 17: Line 17:
 $ init_mbj_lapw $ init_mbj_lapw
 </code> </code>
-  - 一度だけSCF計算を再実行する<code>+  - 一度だけSCF計算を再実行し、ブロイデン法の履歴を削除する<code>
 $ run_lapw -p -i 1 -NI $ run_lapw -p -i 1 -NI
-</code> 
-  - ブロイデン法の履歴を削除し、再度mBJの設定を完了させる<code> 
 $ rm *.broyd* $ rm *.broyd*
 +</code>
 +  - mBJの設定を完了させる<code>
 $ init_mbj_lapw $ init_mbj_lapw
 </code> </code>
Line 39: Line 39:
 これでTB-mBJポテンシャルを使った計算は終わり。 これでTB-mBJポテンシャルを使った計算は終わり。
 あとは状態密度やバンドを描いて、電子状態の変化を確認する。 あとは状態密度やバンドを描いて、電子状態の変化を確認する。
 +x_lapwの実行時に特になにかフラグをつけなくても、自動的にmBJポテンシャルを考慮した計算になる。
  
 適切なTB-mBJポテンシャルの選び方は、上で書いた論文の値をそのまま使うのも一つであるが、光学伝導度測定などの実験で得られたバンドギャップを再現する大きさに合わせるのが手っ取り早い。 適切なTB-mBJポテンシャルの選び方は、上で書いた論文の値をそのまま使うのも一つであるが、光学伝導度測定などの実験で得られたバンドギャップを再現する大きさに合わせるのが手っ取り早い。
Line 48: Line 49:
 ===== 収束しない場合  ===== ===== 収束しない場合  =====
  
-  - 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) PRATT   0.0   YES  (BROYD/PRATT, extra charge (+1 for additional e), norm)
 0.10            mixing FACTOR for BROYD/PRATT scheme 0.10            mixing FACTOR for BROYD/PRATT scheme
wien2k/tb-mbj交換ポテンシャル.1593685735.txt.gz · Last modified: 2021/06/27 22:01 (external edit)