One-dimensional conduction through supporting electrolytes: two-scale cathodic Debye layer.

Supporting-electrolyte solutions comprise chemically inert cations and anions, produced by salt dissolution, together with a reactive ionic species that may be consumed and generated on bounding ion-selective surfaces (e.g., electrodes or membranes). Upon application of an external voltage, a Faraday current is thereby established. It is natural to analyze this ternary-system process through a one-dimensional transport problem, employing the thin Debye-layer limit. Using a simple model of ideal ion-selective membranes, we have recently addressed this problem for moderate voltages [Yariv and Almog, Phys. Rev. Lett. 105, 176101 (2010)], predicting currents that scale as a fractional power of Debye thickness. We address herein the complementary problem of moderate currents. We employ matched asymptotic expansions, separately analyzing the two inner thin Debye layers adjacent to the ion-selective surfaces and the outer electroneutral region outside them. A straightforward calculation following comparable singular-perturbation analyses of binary systems is frustrated by the prediction of negative ionic concentrations near the cathode. Accompanying numerical simulations, performed for small values of Debye thickness, indicate a number unconventional features occurring at that region, such as inert-cation concentration amplification and electric-field intensification. The current-voltage correlation data of the electrochemical cell, obtained from compilation of these simulations, does not approach a limit as the Debye thickness vanishes. Resolution of these puzzles reveals a transformation of the asymptotic structure of the cathodic Debye layer. This reflects the emergence of an internal boundary layer, adjacent to the cathode, wherein field and concentration scaling differs from those of the Gouy-Chapman theory. The two-scale feature of the cathodic Debye layer is manifested through a logarithmic voltage scaling with Debye thickness. Accounting for this scaling, the complied current-voltage data collapses upon a single curve. This curve practically coincides with an asymptotically calculated universal current-voltage relation.