Ionic Layering and Overcharging in Electrical Double Layers in a Poisson-Boltzmann Model

Electrical double layers (EDLs) play a significant role in a broad range of physical phenomena related to colloidal stability, diffuse-charge dynamics, electrokinetics, and energy storage applications. Recently, it has been suggested that for large ion sizes or multivalent electrolytes, ions can arrange in a layered structure inside the EDLs. However, the widely used Poisson--Boltzmann models for EDLs are unable to capture the details of ion concentration oscillations and the effect of electrolyte valence on such oscillations. Here, by treating a pair of ions as hard spheres below the distance of closest approach and as point charges otherwise, we are able to predict ionic layering without any additional parameters or boundary conditions while still being compatible with the Poisson--Boltzmann framework. Depending on the combination of ion valence, size, and concentration, our model reveals a structured EDL with spatially oscillating ion concentrations. We report the dependence of critical ion concentration, i.e., the ion concentration above which the oscillations are observed, on the counter-ion valence and the ion size. More importantly, our model displays quantitative agreement with the results of computationally intensive models of the EDL. Finally, we analyze the nonequilibrium problem of EDL charging and demonstrate that ionic layering increases the total charge storage capacity and the charging timescale.