「Kerkerの方法」の版間の差分
削除された内容 追加された内容
編集の要約なし |
加筆（表式及び説明などを追加） |
||
1行目:
'''Kerkerの方法'''(Kerker method)：[[セルフコンシステント]]な電子状態計算おいて、[[電荷密度]]を混合する方法の一つ[1]。パラメーターを適宜設定することにより、単純混合による方法より計算の収束を速く出来る（速くならない場合もある）。
この方法での電荷密度混合の表式は、電子状態の計算をセルフコンシステントに進める上で、そのための繰り返し計算におけるn回目の繰り返しでの電荷密度をρ(n)とすると、
:<math> \rho (n+1) = \rho_{\mathrm{in}} (n) + \alpha { | \mathbf{G} |^2 \over { | \mathbf{G} |^2 + \lambda^2 } } \{ \rho_{\mathrm{out}} (n) - \rho_{\mathrm{in}} (n) \} </math>
となる。<b>G</b>は[[逆格子ベクトル]]、α、λは適当な定数である。電荷密度ρ(n)に付いている添え字in、outは、それぞれn回目の繰り返しでの入力と出力の電荷密度であることを意味している。λ=0の時、上式は、
:<math> \rho (n+1) = \rho_{\mathrm{in}} (n) + \alpha \{ \rho_{\mathrm{out}} (n) - \rho_{\mathrm{in}} (n) \} = (1 - \alpha) \rho_{\mathrm{in}} (n) + \alpha \rho_{\mathrm{in}} (n) </math>
となり。これは<b>単純混合</b>の表式に帰着する。定数、α、λは計算が速く収束するような値を選ぶが、それは扱う系や計算条件に依存する。
==参考文献==
| |||