Differential capacitance of the electric double layer: The interplay between ion finite size and dielectric decrement

Differential capacitance of curved electrodes: role of hydration interactions and charge regulation

The functioning of supercapacitors relies on establishing electrostatic double-layer capacitance across a larger surface area, offering numerous advantages over conventional batteries, such as an extended lifespan and elevated safety standards. The differential capacitance is a fundamental property within the electrical double layer, playing a pivotal role in the advancement of electrical double-layer supercapacitors. In addition to electrostatic interactions, multiple theoretical and experimental studies have indicated that the differential capacitance is influenced by factors such as the physical structure of the electrode, solvent-mediated hydration interactions, and the specific type of electrolyte utilized. In this work, we incorporate hydration interactions into the Poisson-Boltzmann theory to explore curved electrodes whose surfaces can be covered by either acidic or basic groups. We examine how the electrostatic interaction, charge regulation, hydration effects, and the finite size of ions collectively modify the differential capacitance. Furthermore, we explore different scenarios of electrode curvature and how it may be used to achieve larger capacitance depending on the electrolyte type and pH.

Diffusioosmotic flow reversals due to ion-ion electrostatic correlations

Existing theories of diffusioosmosis have neglected ion-ion electrostatic correlations, which are important in concentrated electrolytes. Here, we develop a mathematical model to numerically compute the diffusioosmotic mobilities of binary symmetric electrolytes across low to high concentrations in a charged parallel-plate channel. We use the modified Poisson equation to model the ion-ion electrostatic correlations and the Bikerman model to account for the finite size of ions. We report two key findings. First, ion-ion electrostatic correlations can cause a unique reversal in the direction of diffusioosmosis. Such a reversal is not captured by existing theories, occurs at ≈ 0.4 M for a monovalent electrolyte, and at a much lower concentration of ≈ 0.003 M for a divalent electrolyte in a channel with the same surface charge. This highlights that diffusioosmosis of a concentrated electrolyte can be qualitatively different from that of a dilute electrolyte, not just in its magnitude but also its direction. Second, we predict a separate diffusioosmotic flow reversal, which is not due to electrostatic correlations but the competition between the underlying chemiosmosis and electroosmosis. This reversal can be achieved by varying the magnitude of the channel surface charge without changing its sign. However, electrostatic correlations can radically change how this flow reversal depends on the channel surface charge and ion diffusivity between a concentrated and a dilute electrolyte. The mathematical model developed here can be used to design diffusioosmosis of dilute and concentrated electrolytes, which is central to applications such as species mixing and separation, enhanced oil recovery, and reverse electrodialysis.

Analytical solutions for viscoelectric effects in electrokinetic nanochannels

Understanding electrokinetic transport in nanochannels and nanopores is essential for emerging biological and electrochemical applications. The viscoelectric effect is an important mechanism implicated in the increase of local viscosity due to the polarization of a solvent under a strong electric field. However, most analyses of the viscoelectric effect have been limited to numerical analyses. In this work, we present a set of analytical solutions applicable to the physical description of viscoelectric effects in nanochannel electrokinetic systems. To achieve such closed-form solutions, we employ the Debye-Hückel approximation of small diffuse charge layer potentials compared to the thermal potential. We analyze critical parameters, including electroosmotic flow profiles, electroosmotic mobility, flow rate, and channel conductance. We compare and benchmark our analytical solutions with published predictions from numerical models. Importantly, we leverage these analytical solutions to identify essential thermophysical and nondimensional parameters that govern the behavior of these systems. We identify scaling parameters and relations among surface charge density, ionic strength, and nanochannel height.

Surface tension between liquids containing antagonistic and regular salts

We investigate theoretically the surface tension between two immiscible liquids containing antagonistic and regular salts. The mean-field model includes the electrostatic energy and the Gibbs transfer energy of the ions. We find the Donnan potential difference between the liquids and solve the Poisson-Boltzmann equation to calculate the potential and ion density profiles. When the liquids contain only a regular salt, the surface tension increases when compared to the no-salt case; in contrast, the surface tension is reduced when they contain only an antagonistic salt. When both salts are present, the surface tension can either further decrease or increase depending on the preferential solvation of the ions, in agreement with published experimental observations.

Neural network predicts ion concentration profiles under nanoconfinement

Modeling the ion concentration profile in nanochannel plays an important role in understanding the electrical double layer and electro-osmotic flow. Due to the non-negligible surface interaction and the effect of discrete solvent molecules, molecular dynamics (MD) simulation is often used as an essential tool to study the behavior of ions under nanoconfinement. Despite the accuracy of MD simulation in modeling nanoconfinement systems, it is computationally expensive. In this work, we propose neural network to predict ion concentration profiles in nanochannels with different configurations, including channel widths, ion molarity, and ion types. By modeling the ion concentration profile as a probability distribution, our neural network can serve as a much faster surrogate model for MD simulation with high accuracy. We further demonstrate the superior prediction accuracy of neural network over XGBoost. Finally, we demonstrated that neural network is flexible in predicting ion concentration profiles with different bin sizes. Overall, our deep learning model is a fast, flexible, and accurate surrogate model to predict ion concentration profiles in nanoconfinement.

Incorporating ion-specific van der Waals and soft repulsive interactions in the Poisson-Boltzmann theory of electrical double layers

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. Over the last century, the development of Poisson-Boltzmann (PB) models and their modified versions have provided significant physical insight into the structure and dynamics of the EDL. Incorporation of physics such as finite-ion-size effects, dielectric decrement, and ion-ion correlations has made such models increasingly accurate when compared to more computationally expensive approaches such as molecular simulations and classical density functional theory. However, a prominent knowledge gap has been the exclusion of van der Waals (vdW) and soft repulsive interactions in modified PB models. Although short-ranged as compared to electrostatic interactions, we show here that vdW and soft repulsive 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 an electrolyte. 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 ions such as Na, K+, Mg2+, Cl-, and SO42-. 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, leads 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.

Oxide- and Silicate-Water Interfaces and Their Roles in Technology and the Environment

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by "anomalies" or "non-idealities" such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide- and silicate-water interfaces. This critical review discusses how science progresses from understanding ideal solid-water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. We anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

Nonlinear dielectric decrement of electrolyte solutions: An effective medium approach

HYPOTHESIS

The dielectric constant of an electrolyte solution, which determines electrostatic interactions between colloids and interfaces, depends nonlinearly on the salinity and also on the type of salt. The linear decrement at dilute solutions is due to the reduced polarizability in the hydration shell around an ion. However, the full hydration volume cannot explain the experimental solubility, which indicates the hydration volume should decrease at high salinity. Volume reduction of the hydration shell is supposed to weaken dielectric decrement and thus should be relevant to the nonlinear decrement.

SIMULATIONS

According to the effective medium theory for the permittivity of heterogeneous media, we derive an equation which relates the dielectric constant with the dielectric cavities created by the hydrated cations and anions, and the effect of partial dehydration at high salinity is taken into account.

FINDINGS

Analysis of experiments on monovalent electrolytes suggests that weakened dielectric decrement at high salinity originates primarily from the partial dehydration. Furthermore, the onset volume fraction of the partial dehydration is found to be salt-specific, and is correlated with the solvation free energy. Our results suggest that while the reduced polarizability of the hydration shell determines the linear dielectric decrement at low salinity, ion-specific tendency of dehydration is responsible for nonlinear dielectric decrement at high salinity.

Electrostatics of soft (bio)interfaces: Corrections of mean-field Poisson-Boltzmann theory for ion size, dielectric decrement and ion-ion correlation

HYPOTHESIS

Electrostatics of soft (ion-permeable) (bio)particles (e.g. microorganisms, core/shell colloids) in aqueous electrolytes is commonly formulated by the mean-field Poisson-Boltzmann theory and integration of the charge contributions from electrolyte ions and soft material. However, the effects connected to the size of the electrolyte ions and that of the structural charges carried by the particle, to dielectric decrement and ion-ion correlations on soft interface electrostatics have been so far considered at the margin, despite the limits of the Gouy theory for condensed and/or multivalent electrolytes.

EXPERIMENTS

Accordingly, we modify herein the Poisson-Boltzmann theory for core/shell (bio)interfaces to include the aforementioned molecular effects considered separately or concomitantly. The formalism is applicable for poorly to highly charged particles in the thin electric double layer regime and to unsymmetrical multivalent electrolytes.

FINDINGS

Computational examples of practical interests are discussed with emphasis on how each considered molecular effect or combination thereof affects the interfacial potential distribution depending on size and valence of cations and anions, size of particle charges, length scale of ionic correlations and shell-to-Debye layer thickness ratio. The origins of here-evidenced pseudo-harmonic potential profile and ion size-dependent screening of core/shell particle charges are detailed. In addition, the existence and magnitude of the Donnan potential when reached in the shell layer are shown to depend on the excluded volumes of the electrolyte ions.

Thermodynamics of Charge Regulation during Ion Transport through Silica Nanochannels

Ion-surface interactions can alter the properties of nanopores and dictate nanofluidic transport in engineered and biological systems central to the water-energy nexus. The ion adsorption process, known as "charge regulation", is ion-specific and is dependent on the extent of confinement when the electric double layers (EDLs) between two charged surfaces overlap. A fundamental understanding of the mechanisms behind charge regulation remains lacking. Herein, we study the thermodynamics of charge regulation reactions in 20 nm SiO₂ channels via conductance measurements at various concentrations and temperatures. The effective activation energies (E_a) for ion conductance at low concentrations (strong EDL overlap) are ∼2-fold higher than at high concentrations (no EDL overlap) for the electrolytes studied here: LiCl, NaCl, KCl, and CsCl. We find that E_a values measured at high concentrations result from the temperature dependence of viscosity and its influence on ion mobility, whereas E_a values measured at low concentrations result from the combined effects of ion mobility and the enthalpy of cation adsorption to the charged surface. Notably, the E_a for surface reactions increases from 7.03 kJ mol^-1 for NaCl to 16.72 ± 0.48 kJ mol^-1 for KCl, corresponding to a difference in surface charge of -8.2 to -0.8 mC m^-2, respectively. We construct a charge regulation model to rationalize the cation-specific charge regulation behavior based on an adsorption equilibrium. Our findings show that temperature- and concentration-dependent conductance measurements can help indirectly probe the ion-surface interactions that govern transport and colloidal interactions at the nanoscale─representing a critical step forward in our understanding of charge regulation and adsorption phenomena under nanoconfinement.

Thermodynamics of Electrolyte Solutions Near Charged Surfaces: Constant Surface Charge vs. Constant Surface Potential

Electric double layers are ubiquitous in science and engineering and are of current interest, owing to their applications in the stabilization of colloidal suspensions and as supercapacitors. While the structure and properties of electric double layers in electrolyte solutions near a charged surface are well characterized, there are subtleties in calculating thermodynamic properties from the free energy of a system with charged surfaces. These subtleties arise from the difference in the free energy between systems with constant surface charge and constant surface potential. In this work, we present a systematic, pedagogical framework to properly account for the different specifications on charged bodies in electrolyte solutions. Our approach is fully variational---that is, all free energies, boundary conditions, relevant electrostatic equations, and thermodynamic quantities are systematically derived using variational principles of thermodynamics. We illustrate our approach by considering a simple electrolyte solution between two charged surfaces using the Poisson--Boltzmann theory. Our results highlight the importance of using the proper thermodynamic potential and provide a general framework for calculating thermodynamic properties of electrolyte solutions near charged surfaces. Specifically, we present the calculation of the pressure and the surface tension between two charged surfaces for different boundary conditions, including mixed boundary conditions.

On the importance of the electric double layer structure in aqueous electrocatalysis

To design electrochemical interfaces for efficient electric-chemical energy interconversion, it is critical to reveal the electric double layer (EDL) structure and relate it with electrochemical activity; nonetheless, this has been a long-standing challenge. Of particular, no molecular-level theories have fully explained the characteristic two peaks arising in the potential-dependence of the EDL capacitance, which is sensitively dependent on the EDL structure. We herein demonstrate that our first-principles-based molecular simulation reproduces the experimental capacitance peaks. The origin of two peaks emerging at anodic and cathodic potentials is unveiled to be an electrosorption of ions and a structural phase transition, respectively. We further find a cation complexation gradually modifies the EDL structure and the field strength, which linearly scales the carbon dioxide reduction activity. This study deciphers the complex structural response of the EDL and highlights its catalytic importance, which bridges the mechanistic gap between the EDL structure and electrocatalysis.

Implicit Solvation Methods for Catalysis at Electrified Interfaces

Implicit solvation is an effective, highly coarse-grained approach in atomic-scale simulations to account for a surrounding liquid electrolyte on the level of a continuous polarizable medium. Originating in molecular chemistry with finite solutes, implicit solvation techniques are now increasingly used in the context of first-principles modeling of electrochemistry and electrocatalysis at extended (often metallic) electrodes. The prevalent ansatz to model the latter electrodes and the reactive surface chemistry at them through slabs in periodic boundary condition supercells brings its specific challenges. Foremost this concerns the difficulty of describing the entire double layer forming at the electrified solid–liquid interface (SLI) within supercell sizes tractable by commonly employed density functional theory (DFT). We review liquid solvation methodology from this specific application angle, highlighting in particular its use in the widespread ab initio thermodynamics approach to surface catalysis. Notably, implicit solvation can be employed to mimic a polarization of the electrode’s electronic density under the applied potential and the concomitant capacitive charging of the entire double layer beyond the limitations of the employed DFT supercell. Most critical for continuing advances of this effective methodology for the SLI context is the lack of pertinent (experimental or high-level theoretical) reference data needed for parametrization.

Cation Overcrowding Effect on the Oxygen Evolution Reaction

The influence of electrolyte ions on the catalytic activity of electrode/electrolyte interfaces is a controversial topic for many electrocatalytic reactions. Herein, we focus on an effect that is usually neglected, namely, how the local reaction conditions are shaped by nonspecifically adsorbed cations. We scrutinize the oxygen evolution reaction (OER) at nickel (oxy)hydroxide catalysts, using a physicochemical model that integrates density functional theory calculations, a microkinetic submodel, and a mean-field submodel of the electric double layer. The aptness of the model is verified by comparison with experiments. The robustness of model-based insights against uncertainties and variations in model parameters is examined, with a sensitivity analysis using Monto Carlo simulations. We interpret the decrease in OER activity with the increasing effective size of electrolyte cations as a consequence of cation overcrowding near the negatively charged electrode surface. The same reasoning could explain why the OER activity increases with solution pH on the RHE scale and why the OER activity decreases in the presence of bivalent cations. Overall, this work stresses the importance of correctly accounting for local reaction conditions in electrocatalytic reactions to obtain an accurate picture of factors that determine the electrode activity.

Monte Carlo simulations and mean-field modeling of electric double layers at weakly and moderately charged spherical macroions

Monte Carlo simulations are employed to determine the differential capacitance of an electric double layer formed by small size-symmetric anions and cations in the vicinity of weakly to moderately charged macroions. The influence of interfacial curvature is deduced by investigating spherical macroions, ranging from flat to moderately curved. We also calculate the differential capacitance using a previously developed mean-field model where, in addition to electrostatic interactions, the excluded volumes of the ions are taken into account using either the lattice-gas or the Carnahan-Starling equation of state. For both equations of state, we compare the mean-field model for arbitrary curvature with a recently developed second-order curvature expansion. Our Monte Carlo simulations predict an increase in the differential capacitance with growing macroion curvature if the surface charge density is small, whereas for moderately charged macroions the differential capacitance passes through a local minimum. Both mean-field models tend to somewhat overestimate the differential capacitance as compared with Monte Carlo simulations. At the same time, they do reproduce the curvature dependence of the differential capacitance, especially for small surface charge density. Our study suggests that the quality of mean-field modeling does not worsen when weakly or moderately charged macroions exhibit spherical curvature.

The differential capacitance as a probe for the electric double layer structure and the electrolyte bulk composition

In this work, we theoretically study the differential capacitance of an aqueous electrolyte in contact with a planar electrode, using classical density functional theory, and show how this measurable quantity can be used as a probe to better understand the structure and composition of the electric double layer at play. Specifically, we show how small trace amounts of divalent ions can influence the differential capacitance greatly and also how small ions dominate its behavior for high electrode potentials. In this study, we consider primitive model electrolytes and not only use the standard definition of the differential capacitance but also derive a new expression from mechanical equilibrium in a planar geometry. This expression reveals explicitly that the first layer of ions near the charged surface is key to its understanding. Our insights might be used as a guide in experiments to better understand the electrolyte-electrode interface as well as the (composition of the) bulk electrolyte.

Effects of Confinement and Ion Adsorption in Ionic Liquid Supercapacitors with Nanoporous Electrodes

We investigate the effects of pore size and ion adsorption on the room-temperature ionic liquid capacitor with nanoporous electrodes, with a focus on optimizing the capacitance and energy storage. Using a recently developed modified BSK model accounting for both ion correlations and nonelectrostatic interactions, we find that ion crowding proximate to the electrode surface induced by the spontaneous charge separation due to strong ion correlations is responsible for the anomalous increase in the capacitance with decreasing pore sizes observed in experiments. Reducing the strength of ion correlations increases the capacitance and suppresses the anomalous size dependence. For a given pore size, the capacitance peak diverges when the ion correlation strength α reaches a critical value, α_sc,L. The capacitance peak shifts to smaller pore size as α decreases because of rapid decrease of α_sc,L with decreasing pore size. Asymmetric preferential ion adsorption is shown to lead to significantly enhanced energy storage close to the transition point for any pore sizes. For a given correlation strength, the energy storage is optimal at a pore size where α = α_sc,L.

Electrokinetic actuation of an uncharged polarizable dielectric droplet in charged hydrogel medium

Electrokinetic transport of an uncharged nonconducting microsized liquid droplet in a charged hydrogel medium is studied. Dielectric polarization of the liquid drop under the action of an externally imposed electric field induces a non-homogeneous charge density at the droplet surface. The interactions of the induced surface charge of the droplet with the immobile charges of the hydrogel medium generates an electric force to the droplet, which actuates the drop through the charged hydrogel medium. A numerical study based on the first principle of electrokinetics is adopted. Dependence of the droplet velocity on its dielectric permittivity, bulk ionic concentration, and immobile charge density of the gel is analyzed. The surface conduction is significant in presence of charged gel, which creates a concentration polarization. The impact of the counterion saturation in the Debye layer due to the dielectric decrement of the medium is addressed. The modified Nernst-Planck equation for ion transport and the Poisson equation for the electric field is considered to take into account the dielectric polarization. A quadrupolar vortex around the uncharged droplet is observed when the gel medium is considered to be uncharged, which is similar to the induced charge electroosmosis around an uncharged dielectric colloid in free-solution. We find that the induced charge electrokinetic mechanism creates a strong recirculation of liquid within the droplet and the translational velocity of the droplet strongly depends on its size for the dielectric droplet embedded in a charged gel medium.

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.

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.

Charging Dynamics of Overlapping Double Layers in a Cylindrical Nanopore

The charging of electrical double layers inside a cylindrical pore has applications to supercapacitors, batteries, desalination and biosensors. The charging dynamics in the limit of thin double layers, i.e., when the double layer thickness is much smaller than the pore radius, is commonly described using an effective RC transmission line circuit. Here, we perform direct numerical simulations (DNS) of the Poisson-Nernst-Planck equations to study the double layer charging for the scenario of overlapping double layers, i.e., when the double layer thickness is comparable to the pore radius. We develop an analytical model that accurately predicts the results of DNS. Also, we construct a modified effective circuit for the overlapping double layer limit, and find that the modified circuit is identical to the RC transmission line but with different values and physical interpretation of the capacitive and resistive elements. In particular, the effective surface potential is reduced, the capacitor represents a volumetric current source, and the charging timescale is weakly dependent on the ratio of the pore radius and the double layer thickness.

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.

Capacitance of thin films containing polymerized ionic liquids

Electrode-polymer interfaces dictate many of the properties of thin films such as capacitance, the electric field experienced by polymers, and charge transport. However, structure and dynamics of charged polymers near electrodes remain poorly understood, especially in the high concentration limit representative of the melts. To develop an understanding of electric field-induced transformations of electrode-polymer interfaces, we have studied electrified interfaces of an imidazolium-based polymerized ionic liquid (PolyIL) using combinations of broadband dielectric spectroscopy, specular neutron reflectivity, and simulations based on the Rayleigh's dissipation function formalism. Overall, we obtained the camel-shaped dependence of the capacitance on applied voltage, which originated from the responses of an adsorbed polymer layer to applied voltages. This work provides additional insights related to the effects of molecular weight in affecting structure and properties of electrode-polymer interfaces, which are essential for designing next-generation energy storage and harvesting devices.

Effects of Kinetic Dielectric Decrement on Ion Diffusion and Capacitance in Electrochemical Systems

Diffusion of ionic components in electrolytes not only eliminates the gradients of ionic concentrations but also alters the local dielectric environment, and the coupling effect between kinetic dielectric decrement and ionic concentration gradient on the diffusion dynamics is not well understood. Herein, taking the charging process in electrical double layer systems as a case study, we conduct a multiscale investigation of ion diffusions in aqueous electrolytes by combining the dynamic density functional theory and an ion-concentration-dependent dielectric constant model. By properly considering the time evolutions of local dielectric constant coupled with ion density, we report an interesting phenomenon on the suppression of surface charge density that is not captured by conventional models. In addition, we show that the usage of aqueous electrolyte with small dielectric decrement coefficients promotes the capacitance, in quantitative agreement with experimental measurements.

Effects of Surface Transition and Adsorption on Ionic Liquid Capacitors

Room-temperature ionic liquids (RTILs) are synthetic electrolytes with electrochemical stability superior to that of conventional aqueous-based electrolytes, allowing a significantly enlarged electrochemical window for application as capacitors. In this study, we propose a variant of an existing RTIL model for solvent-free RTILs, accounting for both ion-ion correlations and nonelectrostatic interactions. Using this model, we explore the phenomenon of spontaneous surface charge separation in RTIL capacitors and find that this transition is a common feature for realistic choices of the model parameters in most RTILs. In addition, we investigate the effects of asymmetric preferential ion adsorption on this charge separation transition and find that proximity of the transition in this case can result in greatly enhanced energy storage. Our work suggests that differential chemical treatment of electrodes can be a simple and useful means for optimizing energy storage in RTIL capacitors.

Continuum models of the electrochemical diffuse layer in electronic-structure calculations

Continuum electrolyte models represent a practical tool to account for the presence of the diffuse layer at electrochemical interfaces. However, despite the increasing popularity of these in the field of materials science, it remains unclear which features are necessary in order to accurately describe interface-related observables such as the differential capacitance (DC) of metal electrode surfaces. We present here a critical comparison of continuum diffuse-layer models that can be coupled to an atomistic first-principles description of the charged metal surface in order to account for the electrolyte screening at electrified interfaces. By comparing computed DC values for the prototypical Ag(100) surface in an aqueous solution to experimental data, we validate the accuracy of the models considered. Results suggest that a size-modified Poisson-Boltzmann description of the electrolyte solution is sufficient to qualitatively reproduce the main experimental trends. Our findings also highlight the large effect that the dielectric cavity parameterization has on the computed DC values.

Grand-canonical approach to density functional theory of electrocatalytic systems: Thermodynamics of solid-liquid interfaces at constant ion and electrode potentials

Properties of solid-liquid interfaces are of immense importance for electrocatalytic and electrochemical systems, but modeling such interfaces at the atomic level presents a serious challenge and approaches beyond standard methodologies are needed. An atomistic computational scheme needs to treat at least part of the system quantum mechanically to describe adsorption and reactions, while the entire system is in thermal equilibrium. The experimentally relevant macroscopic control variables are temperature, electrode potential, and the choice of the solvent and ions, and these need to be explicitly included in the computational model as well; this calls for a thermodynamic ensemble with fixed ion and electrode potentials. In this work, a general framework within density functional theory (DFT) with fixed electron and ion chemical potentials in the grand canonical (GC) ensemble is established for modeling electrocatalytic and electrochemical interfaces. Starting from a fully quantum mechanical description of multi-component GC-DFT for nuclei and electrons, a systematic coarse-graining is employed to establish various computational schemes including (i) the combination of classical and electronic DFTs within the GC ensemble and (ii) on the simplest level a chemically and physically sound way to obtain various (modified) Poisson-Boltzmann (mPB) implicit solvent models. The detailed and rigorous derivation clearly establishes which approximations are needed for coarse-graining as well as highlights which details and interactions are omitted in vein of computational feasibility. The transparent approximations also allow removing some of the constraints and coarse-graining if needed. We implement various mPB models within a linear dielectric continuum in the GPAW code and test their capabilities to model capacitance of electrochemical interfaces as well as study different approaches for modeling partly periodic charged systems. Our rigorous and well-defined DFT coarse-graining scheme to continuum electrolytes highlights the inadequacy of current linear dielectric models for treating properties of the electrochemical interface.

Electrical Double Layers: Effects of Asymmetry in Electrolyte Valence on Steric Effects, Dielectric Decrement, and Ion-Ion Correlations

We study the effects of asymmetry in electrolyte valence (i.e., non z: z electrolytes) on mean field theory of the electrical double layer. Specifically, we study the effect of valence asymmetry on finite ion-size effects, the dielectric decrement, and ion-ion correlations. For a model configuration of an electrolyte near a charged surface in equilibrium, we present comprehensive analytical and numerical results for the potential distribution, electrode charge density, capacitance, and dimensionless salt uptake. We emphasize that the asymmetry in electrolyte valence significantly influences the diffuse-charge relations, and prior results reported in the literature are readily extended to non z: z electrolytes. We develop scaling relations and invoke physical arguments to examine the importance of asymmetry in electrolyte valence on the aforementioned effects. We conclude by providing implications of our findings on diffuse-charge dynamics and other electrokinetic phenomena.

Theory of the Double Layer in Water-in-Salt Electrolytes

One challenge in developing the next generation of lithium-ion batteries is the replacement of organic electrolytes, which are flammable and most often contain toxic and thermally unstable lithium salts, with safer, environmentally friendly alternatives. Recently developed water-in-salt electrolytes (WiSEs), which are nonflammable, nontoxic, and also have enhanced electrochemical stability, are promising alternatives. In this work, we develop a simple modified Poisson-Fermi theory for WiSEs, which demonstrates the fine interplay between electrosorption, solvation, and ion correlations. The phenomenological parameters are extracted from molecular dynamics simulations, also performed here. The theory reproduces the WiSEs' electrical double-layer structure with remarkable accuracy.

Modeling the camel-to-bell shape transition of the differential capacitance using mean-field theory and Monte Carlo simulations

Mean-field electrostatics is used to calculate the differential capacitance of an electric double layer formed at a planar electrode in a symmetric 1:1 electrolyte. Assuming the electrolyte is also ion-size symmetric, we derive analytic expressions for the differential capacitance valid up to fourth order in the surface charge density or surface potential. Our mean-field model accounts exclusively for electrostatic interactions but includes an arbitrary non-ideality in the mixing entropy of the mobile ions. The ensuing criterion for the camel-to-bell shape transition of the differential capacitance is analyzed using commonly used mixing models (one based on a lattice gas and the other based on the Carnahan-Starling equation of state) and compared with Monte Carlo simulations. We observe a reasonable agreement between all our mean-field models and the simulation data for the camel-to-bell shape transition. The absolute value of the differential capacitance for an uncharged (or weakly charged) electrode is, however, not reproduced by our mean-field approaches, not even upon introducing a Stern layer with a thickness equal of the ion radius. We show that, if a Stern layer is introduced, its thickness dependence on the ion size is non-monotonic or, depending on the salt concentration, even inversely proportional.

General theory of charge regulation and surface differential capacitance

A generalization of the mean-field approach will be derived that will take into account the ion-ion as well as ion-surface non-electrostatic effects on an equal footing, being based on the bulk and surface equations of state in the absence of electrostatic interactions. This approach will be applied to the analysis of a single planar surface with dissociable sites with several models of the specific ion-surface non-electrostatic interactions, providing a general thermodynamic insight into the characteristics of the surface differential capacitance. The ion-surface interactions and ion-ion packing considerations at the surface will be shown to be more relevant than the bulk packing constraints for ions vicinal to the surface, as well as to set in prior to the conditions where the bulk packing constraints would become relevant.

Diffuse double-layer structure in mixed electrolytes considering ions as dielectric spheres

The structure of the diffuse part of the electric double layer at solid-electrolyte solution interfaces is examined using a theoretical model that takes into account the finite ion size by modeling the solution as a suspension of polarizable insulating spheres in water. This formalism is applied to mixed electrolyte solutions using the "Boublik-Mansoori-Carnahan-Starling-Leland" (BMCSL) theory for the steric interactions among ions. It is shown that the ionic size differences have a strong bearing on the diffuse part of the electric double-layer structure of these systems. Moreover, for strong potential values, the different size-related effects become important even for binary electrolyte solutions due to the presence of H⁺ and OH^- ions that are substantially smaller than hydrated ions originated from salt dissociation. The obtained results display some of the qualitative features observed in experiments on aqueous systems that are generally interpreted in terms of totally different mechanisms.

Analytical Interfacial Layer Model for the Capacitance and Electrokinetics of Charged Aqueous Interfaces

We construct an analytical model to account for the influence of the subnanometer-wide interfacial layer on the differential capacitance and the electro-osmotic mobility of solid-electrolyte interfaces. The interfacial layer is incorporated into the Poisson-Boltzmann and Stokes equations using a box model for the dielectric properties, the viscosity, and the ionic potential of mean force. We calculate the differential capacitance and the electro-osmotic mobility as a function of the surface charge density and the salt concentration, both with and without steric interactions between the ions. We compare the results from our theoretical model with experimental data on a variety of systems (graphite and metallic silver for capacitance and titanium oxide and silver iodide for electro-osmotic data). The differential capacitance of silver as a function of salinity and surface charge density is well reproduced by our theory, using either the width of the interfacial layer or the ionic potential of mean force as the only fitting parameter. The differential capacitance of graphite, however, needs an additional carbon capacitance to explain the experimental data. Our theory yields a power-law dependence of the electro-osmotic mobility on the surface charge density for high surface charges, reproducing the experimental data using both the interfacial parameters extracted from molecular dynamics simulations and fitted interfacial parameters. Finally, we examine different types of hydrodynamic boundary conditions for the power-law behavior of the electro-osmotic mobility, showing that a finite-viscosity layer explains the experimental data better than the usual hydrodynamic slip boundary condition. Our analytical model thus allows us to extract the properties of the subnanometer-wide interfacial layer by fitting to macroscopic experimental data.

The effects of ion adsorption on the potential of zero charge and the differential capacitance of charged aqueous interfaces

Using a box profile approximation for the non-electrostatic surface adsorption potentials of anions and cations, we calculate the differential capacitance of aqueous electrolyte interfaces from a numerical solution of the Poisson-Boltzmann equation, including steric interactions between the ions and an inhomogeneous dielectric profile. Preferential adsorption of the positive (negative) ion shifts the minimum of the differential capacitance to positive (negative) surface potential values. The trends are similar for the potential of zero charge; however, the potential of zero charge does not correspond to the minimum of the differential capacitance in the case of asymmetric ion adsorption, contrary to the assumption commonly used to determine the potential of zero charge. Our model can be used to obtain more accurate estimates of ion adsorption properties from differential capacitance or electrocapillary measurements. Asymmetric ion adsorption also affects the relative heights of the characteristic maxima in the differential capacitance curves as a function of the surface potential, but even for strong adsorption potentials the effect is small, making it difficult to reliably determine the adsorption properties from the peak heights.

On the dielectric decrement of electrolyte solutions: a dressed-ion theory analysis

Based on the dressed-ion theory and a simple physical argument regarding the conductivity of the solution, we derive a relation between the ionic strength and dielectric constant of an electrolyte solution. At its simplest, this model gives the dielectric constant at low ionic strength I as ε_r(I) = ε_r(0)(1 + αI)^-1, where α (the excess polarization) is directly related to the dressed-ion charge. One contribution to the origin of the dielectric decrement is thus seen to stem from the electrostatic screening of the ions in solution, with no solvent contributions necessary.

Incorporation of ion and solvent structure into mean-field modeling of the electric double layer

An electric double layer forms when the small mobile ions of an electrolyte interact with an extended charged object, a macroion. The competition between electrostatic attraction and translational entropy loss of the small ions results in a diffuse layer of partially immobilized ions in the vicinity of the macroion. Modeling structure and energy of the electric double layer has a long history that has lead to the classical Poisson-Boltzmann theory and numerous extensions that account for ion-ion correlations and structural ion and solvent properties. The present review focuses on approaches that instead of going beyond the mean-field character of Poisson-Boltzmann theory introduce structural details of the ions and the solvent into the Poisson-Boltzmann modeling framework. The former include not only excluded volume effects but also the presence of charge distributions on individual ions, spatially extended ions, and internal ionic degrees of freedom. The latter treat the solvent either explicitly as interacting Langevin dipoles or in the form of effective non-electrostatic interactions, in particular Yukawa interactions, that are added to the Coulomb potential. We discuss how various theoretical models predict structural properties of the electric double layer such as the differential capacitance and compare some of these predictions with computer simulations.

Differential capacitance of the diffuse double layer at electrode-electrolyte interfaces considering ions as dielectric spheres: Part I. Binary electrolyte solutions

A full theoretical account of the differential capacitance of the diffuse part of the electric double layer at electrode-electrolyte solution interfaces is presented. It builds upon the standard electrokinetic model adding all the additional effects related to the finite ionic size. This includes steric interactions among ions by means of either the Bikerman or Carnahan-Starling expressions and all the permittivity related effects that arise when ions are represented as dielectric spheres. These include the solution permittivity dependence on the local ionic concentration, calculated by means of the Maxwell mixture formula, and two additional forces acting on the ions, namely the Born and the dielectrophoretic forces that depend on the permittivity and the electric field gradients, respectively. The obtained results show that the diffuse double layer behavior is sufficient to qualitatively account for the observed differential capacitance dependence on the electrode voltage. Moreover, when combined with an inner layer capacitance and using the Carnahan-Starling expression, a remarkably good quantitative agreement is achieved.

Evaluating continuum solvation models for the electrode-electrolyte interface: Challenges and strategies for improvement

Ab initio modeling of electrochemical systems is becoming a key tool for understanding and predicting electrochemical behavior. Development and careful benchmarking of computational electrochemical methods are essential to ensure their accuracy. Here, using charging curves for an electrode in the presence of an inert aqueous electrolyte, we demonstrate that most continuum models, which are parameterized and benchmarked for molecules, anions, and cations in solution, undersolvate metal surfaces, and underestimate the surface charge as a function of applied potential. We examine features of the electrolyte and interface that are captured by these models and identify improvements necessary for realistic electrochemical calculations of metal surfaces. Finally, we reparameterize popular solvation models using the surface charge of Ag(100) as a function of voltage to find improved accuracy for metal surfaces without significant change in utility for molecular and ionic solvation.

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.