Influence of nonelectrostatic ion-ion interactions on double-layer capacitance

How charge regulation and ion-surface affinity affect the differential capacitance of an electrical double layer

The differential capacitance of an electrical double layer is a topic of great importance to develop more efficient and environment-friendly energy storage devices: electric double layer supercapacitors. In addition to the bare electrostatic interactions, recent experimental and computational studies suggest that electrodes covered by ionizable groups do interact selectively with specific ion types, an effect that can increase the maximal conductivity and voltage of a supercapacitor. Inspired by this, in the present work we investigate how ion-specific non-electrostatic interactions modify the differential capacitance of a flat electrode whose surface is covered by ionizable groups subject to a charge regulation process. The incorporation of hydration interactions by means of ion-specific Yukawa potential into the Poisson-Boltzmann theory allows our model to describe different scenarios of ion-surface affinity and, hence, the selective depletion or accumulation of specific ion types close to a charged surface. We obtained larger capacitance values when considering electrodes that favor the accumulation of cations and the depletion of anions.

Review and Modification of Entropy Modeling for Steric Effects in the Poisson-Boltzmann Equation

The classical Poisson-Boltzmann model can only work when ion concentrations are very dilute, which often does not match the experimental conditions. Researchers have been working on the modification of the model to include the steric effect of ions, which is non-negligible when the ion concentrations are not dilute. Generally the steric effect was modeled to correct the Helmholtz free energy either through its internal energy or entropy, and an overview is given here. The Bikerman model, based on adding solvent entropy to the free energy through the concept of volume exclusion, is a rather popular steric-effect model nowadays. However, ion sizes are treated as identical in the Bikerman model, making an extension of the Bikerman model to include specific ion sizes desirable. Directly replacing the ions of non-specific size by specific ones in the model seems natural and has been accepted by many researchers in this field. However, this straightforward modification does not have a free energy formula to support it. Here modifications of the Bikerman model to include specific ion sizes have been developed iteratively, and such a model is achieved with a guarantee that: (1) it can approach Boltzmann distribution at diluteness; (2) it can reach saturation limit as the reciprocal of specific ion size under extreme electrostatic conditions; (3) its entropy can be derived by mean-field lattice gas model.

Debye-Hückel theory of weakly curved macroions: Implementing ion specificity through a composite Coulomb-Yukawa interaction potential

The free energy of a weakly curved, isolated macroion embedded in a symmetric 1:1 electrolyte solution is calculated on the basis of linear Debye-Hückel theory, thereby accounting for nonelectrostatic Yukawa pair interactions between the mobile ions and of the mobile ions with the macroion surface, present in addition to the electrostatic Coulomb potential. The Yukawa interactions between anion-anion, cation-cation, and anion-cation pairs are independent from each other and serve as a model for solvent-mediated ion-specific effects. We derive expressions for the free energy of a planar surface, the spontaneous curvature, the bending stiffness, and the Gaussian modulus. It is shown that a perturbation expansion, valid if the Yukawa interactions make a small contribution to the overall free energy, yields simple analytic results that exhibit good agreement with the general free energy over the range of experimentally relevant interaction parameters.

Surface tension of a Yukawa fluid according to mean-field theory

Yukawa fluids consist of particles that interact through a repulsive or attractive Yukawa potential. A surface tension arises at the walls of the container that encloses the fluid or at the interface between two coexisting phases. We calculate that surface tension on the level of mean-field theory, thereby either ignoring the particle size (ideal Yukawa fluid) or accounting for a non-vanishing particle size through a nonideal contribution to the free energy, exemplified either on the level of a lattice gas (lattice Yukawa fluid) or based on the Carnahan-Starling equation of state (Carnahan-Starling Yukawa fluid). Our mean-field results, which do not rely on assuming small gradients of the particle concentrations, become exact in the limit of large temperature and large screening length. They are calculated numerically in the general case and analytically in the two limits of small particle concentration and close to the critical point for a phase-separating system. For a sufficiently small particle concentration, our predicted surface tension is accurate whereas for a phase boundary, we expect good agreement with exact calculations in the limit of a large screening length and if the mean-field model employs the Carnahan-Starling equation of state.

Solvent-mediated nonelectrostatic ion-ion interactions predicting anomalies in electrophoresis

We study the effects of solvent-mediated nonelectrostatic ion-ion interactions on electrophoretic mobility of a charged spherical particle. To this end, we consider the case of low surface electrostatic potential resulting in the linearization of the governing equations, which enables us to deduce a closed-form analytical solution to the electrophoretic mobility. We subsequently compare our results to the standard model using Henry's approach and report the changes brought about by the nonelectrostatic potential. The classical approach to determine the electrophoretic mobility underpredicts the particle velocity when compared with experiments. We show that this issue can be resolved by taking into account nonelectrostatic interactions. Our analysis further reveals the phenomenon of electrophoretic mobility reversal that has been experimentally observed in numerous previous studies.

Differential capacitance of an electric double layer with asymmetric solvent-mediated interactions: mean-field theory and Monte Carlo simulations

The differential capacitance of an electrical double layer is directly affected by properties of the electrolyte solution such as temperature, salt concentration, ionic size, and solvent structure.

Effects of solvent-mediated nonelectrostatic ion-ion interactions on a streaming potential in microchannels and nanochannels

Here, we capture the consequences of solvent-mediated nonelectrostatic ion-ion interactions, coupled with the considerations of finite-sized effects of the ionic species, on electrokinetic transport in narrow fluidic confinements. We consider pressure-driven flow in microchannels and nanochannels in the presence of electrical double layer effects and analyze the establishment of a streaming potential as mediated by a Yukawa-like pair potential that integrates the ion specificity with the governing nonelectrostatic interactions. We bring out the influences of these interactions on electroviscous effects manifested due to the establishment of the streaming potential. Our considerations provide a plausible explanation for the gross overestimation of electrokinetic energy conversion efficiencies as predicted by classical electrical double layer theories that ignore nonelectrostatic interactions between the ionic species.

The Influence of Dielectric Decrement on Electrokinetics

We treat the dielectric decrement induced by excess ion polarization as a source of ion specificity and explore its impact on electrokinetics. We employ a modified Poisson-Nernst-Planck (PNP) equations accounting for the dielectric decrement. The dielectric decrement is determined by the excess ion polarization parameter α and when α = 0 the standard PNP model is recovered. Our model shows that ions saturate at large zeta potentials (ζ). Because of ion saturation, a condensed counterion layer forms adjacent to the charged surface, introducing a new length scale, the thickness of the condensed layer (l_c ). For the electro-osmotic mobility, the dielectric decrement weakens the electro-osmotic flow owing to the decrease of the dielectric permittivity. At large ζ, when α ≠ 0, the electro-osmotic mobility is found to be proportional to ζ/2, in contrast to ζ predicted by the standard PNP model. This is attributed to ion saturation at large ζ. In terms of the electrophoretic mobility M_e , we carry out both an asymptotic analysis in the thin-double-layer limit and solve the full modified PNP model to compute M_e . Our analysis reveals that the impact of the dielectric decrement is intriguing. At small and moderate ζ, the dielectric decrement decreases M_e with an increasing α. At large ζ, it is well known that the surface conduction becomes significant and plays an important role in determining M_e . It is observed that the dielectric decrement effectively reduces the surface conduction. Hence in stark contrast, M_e increases as α increases. Our predictions of the contrast dependence of the mobility on α at different zeta potentials qualitatively agree with experimental results on the dependence of the mobility among ions and provide a possible explanation for such ion specificity. Finally, the comparisons between the thin-double-layer asymptotic analysis and the full simulations of the modified PNP model suggest that at large ζ the validity of the thin-double-layer approximation is determined by l_c rather than the traditional Debye length.