Kerkerの方法のソースを表示

{{出典の明記|date=2021年4月6日 (火) 11:03 (UTC)}}
'''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>

となる。'''G'''は[[逆格子ベクトル]]、α、λは適当な定数である。電荷密度ρ(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{out}} (n) </math>

となり。これは'''単純混合'''の表式となっている。定数、α、λは計算が速く収束するような値を選ぶが、それは扱う系や計算条件に依存する。

== 参考文献 ==
[1] G. P. Kerker, Phys. Rev. B'''23''' (1981) 3082.

== 関連項目 ==
*[[電荷密度の混合の仕方]]

[[Category:バンド計算]]