[論文レビュー] Ion-ion correlations: an improved one-component plasma correction

Márcia C. Barbosa, Markus Deserno+1 more|arXiv (Cornell University)|Oct 22, 1999

Electrostatics and Colloid Interactions被引用数 40

ひとこと要約

本稿では、中心イオンの周囲の均一な背景電荷を除くDebye-Hückel-Hole-Cavity (DHHC)手法を導入することで、Poisson-Boltzmann理論に対する熱力学的に安定な補正を提案する。これにより、標準的なDebye-Hückel＋ホール(DHH)理論の『構造崩壊』問題が解消される。得られた自由エネルギー汎関数は、凸な過剰自由エネルギー密度から導出されており、帯電したロッド系におけるイオン相関を正確に捉えており、平均場理論のPB理論よりも分子動力学シミュレーションと一致する。

ABSTRACT

Based on a Debye-Hueckel approach to the one-component plasma we propose a new free energy for incorporating ionic correlations into Poisson-Boltzmann like theories. Its derivation employs the exclusion of the charged background in the vicinity of the central ion, thereby yielding a thermodynamically stable free energy density, applicable within a local density approximation. This is an improvement over the existing Debye-Hueckel plus hole theory, which in this situation suffers from a "structuring catastrophe". For the simple example of a strongly charged stiff rod surrounded by its counterions we demonstrate that the Poisson-Boltzmann free energy functional augmented by our new correction accounts for the correlations present in this system when compared to molecular dynamics simulations.

研究の動機と目的

• 一成分プラズマ(OCP)モデルの局所密度近似(LDA)における熱力学的不安定性、特に相関補正付きPoisson-Boltzmann理論における『構造崩壊』を解消すること。
• LDAフレームワークで物理的に整合性があり、凸な自由エネルギー汎関数を提供し、非物理的崩壊を回避するイオン相関のための自由エネルギー汎関数を構築すること。
• 中心イオン周囲の背景電荷のキャビティ除去に基づく、OCP自由エネルギーを修正し、その結果を用いてPoisson-Boltzmann理論にイオン間相関を組み込むこと。
• 帯電した剛直なロッドと反イオンを含むモデル系において、新しい補正を分子動力学シミュレーションと照合すること。

提案手法

• 中心イオンの周囲に半径hのキャビティを設け、その内部の均一な背景電荷を除去することで、熱力学的に安定なDebye-Hückel-Hole-Cavity (DHHC)モデルを導入する。
• キャビティ領域での背景電荷の除外により、標準的なDHH手法が示す凹型自由エネルギーを避ける、凸で熱力学的に安定なOCPの過剰自由エネルギー密度を導出する。
• 得られた自由エネルギーを、Poisson-Boltzmann汎関数の局所密度近似(LDA)補正として用い、過剰自由エネルギー項を修正する。
• 全自由エネルギー汎関数を、平均場PB寄与項と新しいDHHC補正項の和として定式化し、局所イオン密度n(r)に依存する形にする。
• 剛直なロッドのシリンダーモデルにこの汎関数を適用し、モンテカルロ法を用いてイオン分布を解く。
• 修正された汎関数から予測される反イオン分布P(r)を、分子動力学シミュレーションの結果と比較する。

実験結果

リサーチクエスチョン

• RQ1 修正された一成分プラズマモデルは、局所密度近似における熱力学的不安定性（構造崩壊）を、標準的なDebye-Hückel＋ホール(DHH)理論から解消できるか？
• RQ2 提案されたDHHCアプローチは、非一様なイオン系で非物理的崩壊を防ぐ凸な自由エネルギー汎関数を提供するか？
• RQ3 DHHC補正付きPoisson-Boltzmann汎関数は、強帯電ロッド系におけるイオン分布の予測を、平均場PB理論に比べてどの程度改善するか？
• RQ4 新しい相関補正は、多価反イオンのシミュレーションで観察されるイオン凝縮効果を、分子動力学シミュレーションとどれほどよく再現するか？

主な発見

• DHHC手法により、OCPの過剰自由エネルギーが凸になるようにすることで『構造崩壊』が完全に解消され、安定な局所密度近似が可能になる。
• 新しい自由エネルギー汎関数は、特にマクロイオン周囲での反イオン凝縮を高精度で捉えており、帯電ロッド系におけるイオン間相関を的確に表現している。
• 一価ロッド（v=1, ℓB/r₀=3）の場合、DHHC補正PB理論は、径方向分布関数に関して分子動力学シミュレーション結果とよく一致する。
• 三価ロッド（v=3, ℓB/r₀=1）の場合、理論はわずかにイオン凝縮を過剰に予測しているが、これはモデルに含まれないシミュレーション上の硬い反発力の寄与による可能性が大きい。
• 複雑な積分方程式法を用いずに、単純かつ物理的に整合性のある方法で、平均場理論にイオン相関を組み込む手法を提供している。
• DHH理論への最小限の修正で効果が得られるため、既存のPoisson-Boltzmannフレームワークへの実用的かつ容易な拡張が可能である。