  Electrical double layers (EDLs) arise when an electrolyte is in contact with
a charged surface, and are encountered in several application areas including
batteries, supercapacitors, electrocatalytic reactors, and colloids. In the
modeling of EDLs, a prominent knowledge gap has been the exclusion of van der
Waals (vdW) and soft repulsive interactions in modified Poisson-Boltzmann (PB)
theories. Although more short-ranged as compared to electrostatic interactions,
we show here that vdW interactions can play an important role in determining
the structure of the EDL via the formation of a Stern layer and in modulating
the differential capacitance of an electrode in solution. To this end, we
incorporate ion-ion and wall-ion vdW attraction and soft repulsion via a 12-6
Lennard-Jones (LJ) potential, resulting in a modified PB-LJ approach. The
wall-ion LJ interactions were found to have a significant effect on the
electrical potential and concentration profiles, especially close to the wall.
However, ion-ion LJ interactions do not affect the EDL structure at low bulk
ion concentrations (< 1 M). We also derive dimensionless numbers to quantify
the impact of ion-ion and wall-ion LJ interactions on the EDL. Furthermore, in
the pursuit of capturing ion-specific effects, we apply our model by
considering various combinations of ions. We observe how varying parameters
such as the electrolyte concentration and electrode potential affect the
structure of the EDL due to the competition between ion-specific LJ and
electrostatic interactions. Lastly, we show that the inclusion of vdW and soft
repulsion interactions as well as hydration effects lead to a better
qualitative agreement of the PB models with experimental double-layer
differential capacitance data. Overall, the modified PB-LJ approach presented
herein will lead to more accurate theoretical descriptions of EDLs in various
application areas.
