[论文解读] 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}$）附近的密度分布。通过求解从第一Born-Green-Yvon (bgy) 层次推导出的非齐次德拜-休克尔方程，揭示了几何排斥、介电响应（$\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}$ 如何控制双电层与三层电荷密度结构之间的转变。
• 阐明由屏蔽电荷与像电荷共同形成的自洽电荷层在塑造密度分布中的作用。

提出的方法

• 在弱耦合条件下，求解从第一Born-Green-Yvon (bgy) 层次导出的非齐次德拜-休克尔方程。
• 采用库仑对相互作用 $v_{\rm w}({\mathbf{r}};{\mathbf{r}}') = \frac{1}{|{\mathbf{r}}-{\mathbf{r}}'|} - \Delta_{\rm el}\frac{1}{|{\mathbf{r}}-{\mathbf{r}}'^*|}$，其中 $\Delta_{\rm el} = (\epsilon_{\rm w}-1)/(\epsilon_{\rm w}+1)$。
• 引入范围为 $b$ 的硬核排斥势，以防止当 $\epsilon_{\rm w} > 1$ 时电荷与其介电像之间发生坍塌。
• 将密度分布按耦合参数 $\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$。