[論文レビュー] Density profiles in a classical Coulomb fluid near a dielectric wall. I. Mean-field scheme

Jean‐Noël Aqua, F. Cornu|arXiv (Cornell University)|May 29, 2001

Electrostatics and Colloid Interactions参考文献 9被引用数 7

ひとこと要約

本稿は、相対誘電率 $\epsilon_{\rm w}$ を持つ誘電体壁に近接する古典的多成分クーロン流体の密度分布を計算するための平均場解析的枠組みを構築する。第一のボーン＝グリーン＝ヨン階層から導かれる非一様デバイ＝フッケル方程式を解くことで、幾何的反発力、誘電応答（$\epsilon_{\rm w} > 1$ または $< 1$）、自己無撞着な電場層形成が共同で電荷分布を形作り、$\epsilon_{\rm w}$ の変化に伴い二重層構造から三重層構造への遷移を明らかにする。

ABSTRACT

The equilibrium density profiles in a classical multicomponent plasma near a hard wall made with a dielectric material characterized by a relative dielectric constant $\\ew$ are studied from the first Born-Green-Yvon equation combined with Poisson equation in a regime where Coulomb coupling is weak inside the fluid. In order to prevent the collapse between charges with opposite signs or between each charge and its dielectric image inside the wall when $\\ew &gt;1$, hard-core repulsions are added to the Coulomb pair interaction. The charge-image interaction cannot be treated perturbatively and the density profiles vary very fast in the vicinity of the wall when $\\ew \ eq 1$. The formal solution of the associated inhomogeneous Debye-H\\"uckel equations will be given in Paper II, together with a systematic fugacity expansion which allows to retrieve the results obtained from the truncated \\bgy hierarchy. In the present paper the exact density profiles are calculated analytically up to first order in the coupling parameter. The expressions show the interplay between three effects~: the geometric repulsion from the impenetrable wall; the electrostatic effective attraction ($\\ew &gt;1$) or repulsion ($\\ew &lt;1$) due to its dielectric response; and the Coulomb interaction between each charge and the potential drop created by the electric layer which appears as soon as the system is not symmetric. We exhibit how the charge density profile evolves between a structure with two oppositely-charged layers and a three-layer organization when $\\ew$ varies. (The case of two ideally conducting walls will be displayed elsewhere)

研究の動機と目的

• 誘電率 $\epsilon_{\rm w}$ を持つ誘電体壁に近接する古典的多成分プラズマの解析的密度分布を導出すること。
• 弱い結合条件下で、幾何的ハードコア反発力、誘電体像力、クーロン相関の相互作用を平均場近似で扱うこと。
• $\epsilon_{\rm w}$ が二重層構造と三重層構造の間の遷移をどのように制御するかを特定すること。
• スクリーニングと像電荷による自己無撞着な電場層が密度分布をどのように形作るかの役割を確立すること。

提案手法

• 弱い結合下で、第一のボーン＝グリーン＝ヨン（bgy）階層から導かれる非一様デバイ＝フッケル方程式を解く。
• $\Delta_{\rm el} = (\epsilon_{\rm w}-1)/(\epsilon_{\rm w}+1)$ を用いたクーロン対相互作用 $v_{\rm w}({\mathbf{r}};{\mathbf{r}}') = \frac{1}{|{\mathbf{r}}-{\mathbf{r}}'|} - \Delta_{\rm el}\frac{1}{|{\mathbf{r}}-{\mathbf{r}}'^*|}$ を使用する。
• $\epsilon_{\rm w} > 1$ の場合に電荷とその誘電体像の収縮を防ぐために、範囲 $b$ のハードコア反発を導入する。
• 結合パラメータ $\rho \varepsilon_D^2 \ln \varepsilon_D$ のべき級数展開を一次まで行い、主要な補正項を保持する。
• 自己エネルギー $V_{\rm self}(x) = -\Delta_{\rm el}/(4x)$ を有効ポテンシャルへの非摂動的寄与として取り扱う。
• ポアソン方程式と線形応答を用いて、電気的ポテンシャル $\Phi(x)$ と電荷密度分布 $\rho_\gamma(x)$ を導出する。

実験結果

リサーチクエスチョン

• RQ1 壁の誘電率 $\epsilon_{\rm w}$ は、古典的クーロン流体における電荷密度層の形成にどのように影響を与えるか？
• RQ2 ハードウォールからの幾何的反発力、誘電体像力、クーロンスクリーニングの相互作用が、密度分布をどのように決定づけるか？
• RQ3 $\epsilon_{\rm w}$ の増加に伴い、系が二重層構造から三重層構造にどのように遷移するか？
• RQ4 像電荷相互作用から生じる自己無撞着な電場層が、壁近傍の密度分布をどのように修正するか？

主な発見

• $\epsilon_{\rm w} > 1$ の場合、誘電体像は有効な引力を生成し、幾何的反発力と競合することで、壁に近い位置に正電荷密度のピークを持つ三重層構造を形成する。
• $\epsilon_{\rm w} < 1$ の場合、幾何的および電気的反発力が壁近傍の密度を抑制し、分布に単一の最小値が現れる。
• $\epsilon_{\rm w} = 1$ の場合、幾何的反発力のみが作用し、標準的な二重層構造に還元され、既知の結果と一致する。
• 自己エネルギー寄与 $e^2 V_{\rm self}(x)$ は短距離で支配的であり、$\epsilon_{\rm w} > 1$ の場合に密度分布の非単調性を引き起こす要因である。
• 二重層と三重層構造の遷移は、$\Delta_{\rm el}$ とデバイスクリーニング長 $\kappa_D^{-1}$ の競合によって支配され、$\kappa_D b \ll 1$ のとき三重層構造が出現する。
• 密度分布の解析的解は、結合パラメータの一次まで正確であり、$\rho \varepsilon_D^2 \ln \varepsilon_D$ のオーダーの補正を含む。