Electrolytes confined between polarizable surfaces in slit pores with anisotropic permittivity tensor

技术摘要：各向异性介电张量作用下受限在极化表面间的电解质溶液

问题陈述
纳米级受限下的电解质溶液表现出与其体相截然不同的物理性质。近期的实验和模拟证据表明，对于小于约 1 nm 的受限空间，介质的介电常数呈现高度各向异性。具体而言，垂直于限制表面的介电常数分量（$\varepsilon_\perp$）与平行分量（$\varepsilon_\parallel$）相比显著降低。虽然各向同性介电边界对极化表面的影响已得到充分理解，但现有的连续介质和粗粒化理论描述未能充分纳入这种介电各向异性。目前缺乏高效的计算方法来将微观观察与中观建模联系起来，特别是在必须使用张量而非标量来描述介电常数且需考虑极化（介电或金属）边界的系统中。

研究方法
本文提出了一种计算框架，用于模拟受限在具有各向异性介电张量（$\varepsilon = \text{diag}(\varepsilon_\parallel, \varepsilon_\parallel, \varepsilon_\perp)$）的窄缝孔隙中的粗粒化电解质溶液。该方法结合了两种主要技术来求解静电方程：

1. 坐标拉伸变换：为了处理各向异性泊松方程，作者沿限制轴应用了坐标变换：$\tilde{z} = \sqrt{\varepsilon_\parallel/\varepsilon_\perp} \, z$。这将其映射为一个受有效介电常数 $\varepsilon_{\text{eff}} = \sqrt{\varepsilon_\parallel \varepsilon_\perp}$ 控制的各向同性框架。
2. 混合静电求解器：
   • 表面极化：利用 2D 周期性格林函数法求解极化表面的静电方程。该方法利用坐标变换推导出了涉及镜像电荷的极化能（$U_p$）的解析表达式。
   • 直接离子-离子相互作用：直接库仑相互作用能（$U_{ES}$）通过经过板层修正的 3D Ewald 求和进行计算。模拟盒在三个维度上进行了人工复制，并施加了特定约束（$L_z > 2L_{xy}$ 且 $L_z > 2L$），以在变换后的坐标空间内准确应用板层修正。

该方法应用于正则蒙特卡洛（MC）模拟，研究受限在间距为 $L = 1$ nm 的狭缝孔隙中的 1:1 尺寸不对称电解质（阳离子直径 $d_+ = 0.3$ nm，阴离子直径 $d_- = 0.6$ nm）。研究对比了极化介电表面（$\varepsilon_o = 4$）和理想金属表面（$\varepsilon_o = \infty$）的情景。

主要结果
模拟表明，介电各向异性剧烈地重塑了双电层结构，其驱动机制由两个截然不同的机制组成：

• 离子-壁相互作用：解析推导显示，孤立离子靠近平面边界的自镜像势能随距离的标度关系为 $U_{\text{self}}(d) \propto \gamma q^2 / (4\varepsilon_\parallel d)$，其中 $\gamma$ 是介电对比因子。至关重要的是，垂直介电常数 $\varepsilon_\perp$ 在空间标度项中被抵消了。对于金属边界（$\gamma = -1$），单个离子-镜像相互作用变得与 $\varepsilon_\perp$ 完全无关，与均匀体相介质的情况一致。
• 离子-离子相关性：与壁相互作用不同，直接的离子-离子库仑相互作用由于 $\varepsilon_\perp$ 的受抑而显著增强。在变换后的框架中，处于相同 $z$ 坐标处的离子间的面内相互作用势标度为 $1/\sqrt{\varepsilon_\parallel \varepsilon_\perp}$。随着 $\varepsilon_\perp$ 的减小（例如从 78.54 降至 4），相反电荷离子之间的横向相关性显著增强。
• 结构重组：
  • 介电表面：由于镜像排斥作用，大尺寸阴离子仍被固定在其立体接触平面处；然而，增强的面内吸引力迫使较小的阳离子脱离其原本偏好的体相位置，直接锁定在与阴离子相同的 $z$ 平面上。这优化了二维层内的局部电中性。
  • 金属表面：尽管镜像力具有吸引性，但低 $\varepsilon_\perp$ 诱导的强烈横向相关性与之竞争，并最终在结构堆积中占据主导地位。
  • 收敛性：在强各向异性（$\varepsilon_\perp = 4$）条件下，金属边界与介电边界的区别消失了。两种边界条件下的离子密度分布图几乎完全一致，因为直接的面内离子-离子相关性主导了热力学过程，使得边界的具体性质（镜像电荷效应）成为了次要的修正项。

意义与主张
本文主张，垂直介电分量的受抑从根本上改变了受限电解质的静电景观。作者证明，使用简单的各向同性体相介电常数来处理纳米受限水环境，无法捕捉到控制双电层组装的基本物理过程。

其核心意义在于识别了一种由于各向异性介电受抑而产生低维相关区域的机制。这项工作确立了这样一个事实：虽然 $\varepsilon_\perp$ 会影响边界反射系数，但其对结构的影响主要是通过放大直接的面内离子-离子相关性，而非通过改变壁面势能来实现。所提出的方法提供了一个稳健且高效的计算框架，用于探索此类现象，而无需承担显式溶剂计算的沉重负担，为纳米流体通道、超级电容器和电化学储能器件中普遍存在的各向异性行为提供了预测性见解。