双电层三大模型辨析：Helmholtz、Gouy-Chapman 与 Stern 的假设、电容表达式及实测电容曲线的适用边界

金属电极浸入电解液后，金属一侧的过剩电荷只分布在表面几个原子层内，面电荷密度 σ就是这部分电荷按几何面积折算的结果。

溶液一侧要把这份电荷补偿掉，可动的只有离子和水分子：离子靠改变局部浓度贡献净电荷，水分子靠偶极取向贡献极化。两侧电荷符号相反、总量相等，合起来称为双电层。

补偿过程发生在多大的空间范围内，决定了电位从电极表面降到本体溶液的方式。把电极相对参比电极的电位记作 E，溶液本体取作零点，σ 随 E 变化的斜率 C=∂σ/∂E 定义为微分电容，银、汞等理想极化电极在稀电解液中的实测值多在20 μF·cm⁻²上下。

σ 恰好为零的那个电位称为零电荷电位，稀溶液中它对应电容曲线上的极小值。

微分电容测得出来，电荷在空间上如何排列却测不出来。C 的数值只给出 σ 与 E 的换算关系，并不指明反离子停在离电极多远的位置，也不指明有多少离子仍被完整的水合壳层隔开。

双电层模型要补的就是这段空间信息：给定 σ，溶液一侧的净电荷密度沿垂直于电极的坐标怎样分布。

空间分布一旦写定，σ–E 曲线和 C–E 曲线随之确定，模型是否成立因而能用实测电容检验。Helmholtz、Gouy-Chapman 和 Stern 给出的正是三套这样的分布假设；三者描述的是同一个界面，区别在于依次放松了哪一条限制。

放松的对象各不相同：Helmholtz 把反离子固定在一层平面上，Gouy-Chapman 允许它们在热运动中散开，Stern 同时保留有限的离子尺寸和热运动。每放松一项，电容表达式里就多出一个随电位或浓度改变的量，模型能对上的实测区段也随之移动。

Helmholtz 在 1853 年提出的模型最简单：反离子被表面电荷吸引，紧贴电极排成一层，溶液侧的电荷全部集中在这一层上，层外保持电中性。电极表面和这层离子构成两块相距一个离子半径的带电平面，界面因此等价于一只平行板电容器。

平行板电容器的电容只由极板间距和介质决定，Helmholtz 模型给出的表达式因此是C_H=ε_rε₀/d。间距 d 取水合离子半径量级的 0.3 nm，紧密层内水分子取向受限，有效相对介电常数远低于本体的 78，取 6–10 代入后得到 18–30 μF·cm⁻²。这个量级与实测电容相符，估计紧密层电容时至今仍沿用这个式子。

C_H 不随电位改变，也不随电解液浓度改变，表达式里既没有 E 也没有离子浓度。实测曲线在这两个方向上都有变化：同一电极在 1 mM 和 0.1 M 溶液中测得的电容可以相差数倍，电位扫过零电荷电位时曲线还会出现极小。把反离子固定在一层平面上的假设，同时去掉了这两种变化。

Gouy 和 Chapman 在 1910 年和 1913 年换掉了这条假设：离子被当作没有体积的点电荷，在静电吸引和热运动的竞争中散布开，局部浓度按 Boltzmann 分布随电位变化。把这个分布代入 Poisson 方程得到 Poisson–Boltzmann 方程，解出的电位沿远离电极的方向按指数衰减。

衰减到 1/e 所经过的距离称为德拜长度 κ⁻¹，25 ℃ 的 1:1 水溶液中它约等于 0.30 nm 除以浓度的平方根，1 mM 对应 9.6 nm，0.1 M 只剩 0.96 nm。

电荷散布在一段有厚度的区域里，电容由这段厚度决定。Gouy-Chapman 模型给出C_GC=εε₀κ·cosh(zeψ₀/2kT)，零电荷电位处 cosh 项等于 1，电容退回 εε₀/κ⁻¹：1 mM 时约7 μF·cm⁻²，0.1 M 时升到约 72 μF·cm⁻²。浓度和电位这两个变量就此进入表达式。

cosh 项随电位升高增长很快，25 ℃ 时 kT/e 等于 25.7 mV，界面电位降达到 200 mV 时 cosh 项已接近 25，1 mM 溶液算出的电容超过 170 μF·cm⁻²。实测曲线在同样条件下只有几十 μF·cm⁻²，数值随电位增大后反而趋于平缓。

这种发散来自点电荷假设：离子没有体积，靠近电极的一侧就能无限制地堆积，模型里没有约束最近接近距离的条件。Helmholtz 给出了这个距离却丢掉了热运动，Gouy-Chapman 保留了热运动却丢掉了尺寸。两者失效的位置正好互补，前者算不出浓度依赖，后者在高电位和高浓度下算出过大的电容。

Stern 在 1924 年把两条限制同时写进模型：离子有确定的体积，中心最近只能到达距电极一个离子半径的位置，这个平面称为 Stern 面。Stern 面以内没有可动电荷，电位沿距离线性下降；Stern 面以外的离子仍服从 Boltzmann 分布，电位继续按指数衰减。

界面因此由两段串联，紧密层贡献 C_H，扩散层贡献 C_GC，总电容按1/C=1/C_H+1/C_GC给出。串联电容由较小的一项主导：取 C_H 为 20 μF·cm⁻²，1 mM 溶液中 C_GC 只有 7 μF·cm⁻²，总电容降到 5.3 μF·cm⁻²；0.1 M 时 C_GC 升到 72 μF·cm⁻²，总电容回到 15.7 μF·cm⁻²，接近紧密层自身的数值。

浓度改变的只有 κ⁻¹，紧密层的厚度和介电常数不随之改变。稀溶液中扩散层厚、C_GC 小，总电容跟着浓度走；浓溶液中扩散层收缩到零点几纳米，C_GC 远大于 C_H，总电容停在紧密层给出的常数上。Helmholtz 模型只对应后一种极限，Gouy-Chapman 模型只对应前一种极限。

电位的作用同样落在 cosh 项上。零电荷电位处 cosh 项取最小值，串联后的电容在这里出现极小，电位偏离零电荷电位时 C_GC 快速增大，总电容随之升高并逐渐由 C_H 封顶。这个极小在稀溶液中清楚可见，浓度升高后由 C_H 主导而变平。测定零电荷电位因此常在几毫摩尔量级的电解液中进行。

串联表达式对所有离子使用同一个最近接近距离，实测中阴离子的行为却与阳离子不同。Grahame 在 1947 年把紧密层再分成两个平面：内 Helmholtz 面穿过脱去部分水合壳层、直接吸附在电极上的离子中心，外 Helmholtz 面穿过仍保持完整水合壳层的离子中心。

氟离子、氯离子一类阴离子在正电位下常发生特性吸附，位置由外 Helmholtz 面移到内 Helmholtz 面，紧密层的有效厚度随电位和离子种类改变。

串联模型给出的 C–E 曲线只有一个极小，两侧单调上升。银单晶在氟化钠溶液中的实测曲线在稀溶液里却有两个峰，浓度升高时两峰互相靠拢并合成一个峰，这一变化称为驼峰形到钟形的转变。三个经典模型都给不出峰的位置，也给不出合并所需的浓度。

全原子模拟把两个峰分别对应到界面上的两类结构改变。

阴极极化区里，界面水的取向在某个电位附近集体翻转，由一个氢指向溶液的状态转到一个氢指向电极的状态，翻转以一级相变的方式发生，面电荷密度—电位曲线出现 S 形区段，电容在相变电位处升高成峰，峰位几乎不随浓度移动。

阳极极化区里，氟离子部分脱水后在内 Helmholtz 层凝聚，峰位随浓度升高向低电位移动，浓度超过约 0.24 M 后与阴极峰合并。

同一批模拟给出的几何量也与串联模型的设定不同。阳离子在距电极 5.1 Å 和 7.4 Å 处形成两个外 Helmholtz 层，扩散层的德拜长度拟合为 5.6 Å 并且不随本体浓度改变；两层之间的水按氢键网络关联取向，拟合出的有效介电常数取负值。Ag(100) 面的零电荷电位在模拟和实验中都落在 −0.63 V（相对标准氢电极）。

峰的位置还能通过改变离子的水合状态移动。把钠离子与 15-冠-5 络合后，阳离子无法在内 Helmholtz 层保持紧贴，电荷分离距离增大到原来的约 1.5 倍，阴极峰随之移向更负的电位。

三个经典模型把这个距离当作固定值，随组分和电位改变的界面结构因此落在它们的适用范围之外；同一体系中二氧化碳还原的一氧化碳分电流密度随冠醚浓度呈对数下降。

三个模型仍有明确的可用区间：电解液稀、离子不发生特性吸附、电位停留在零电荷电位附近的窄窗口内时，串联表达式给出的电容与实测值在同一量级，浓度依赖的方向也对得上。

窗口之外，界面水的取向、阴离子的特性吸附和离子之间的关联都要单独写进模型，电容曲线上的峰才有对应的结构来源。