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

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.

In-plane structure of the electric double layer in the primitive model using classical density functional theory

The electric double layer (EDL) has a pivotal role in screening charges on surfaces as in supercapacitor electrodes or colloidal and polymer solutions. Its structure is determined by correlations between the finite-sized ionic charge carriers of the underlying electrolyte, and, this way, these correlations affect the properties of the EDL and of applications utilizing EDLs. We study the structure of EDLs within classical density functional theory (DFT) in order to uncover whether a structural transition in the first layer of the EDL that is driven by changes in the surface potential depends on specific particle interactions or has a general footing. This transition has been found in full-atom simulations. Thus far, investigating the in-plane structure of the EDL for the primitive model (PM) using DFT has proved a challenge. We show here that the use of an appropriate functional predicts the in-plane structure of EDLs in excellent agreement with molecular dynamics simulations. This provides the playground to investigate how the structure factor within a layer parallel to a charged surface changes as a function of both the applied surface potential and its separation from the surface. We discuss pitfalls in properly defining an in-plane structure factor and fully map out the structure of the EDL within the PM for a wide range of electrostatic electrode potentials. However, we do not find any signature of a structural crossover and conclude that the previously reported effect is not fundamental but rather occurs due to the specific force field of ions used in the simulations.

Local measures of fluctuations in inhomogeneous liquids: statistical mechanics and illustrative applications

We show in detail how three one-body fluctuation profiles, namely the local compressibility, the local thermal susceptibility, and the reduced density, can be obtained from a statistical mechanical many-body description of classical particle-based systems. We present several different and equivalent routes to the definition of each fluctuation profile, facilitating their explicit numerical calculation in inhomogeneous equilibrium systems. This underlying framework is used for the derivation of further properties such as hard wall contact theorems and novel types of inhomogeneous one-body Ornstein-Zernike equations. The practical accessibility of all three fluctuation profiles is exemplified by grand canonical Monte Carlo simulations that we present for hard sphere, Gaussian core and Lennard-Jones fluids in confinement.

A Gaussian field approach to the planar electric double layer structures in electrolyte solutions

In this work, the planar, electric, double-layer structures of non-polarizable electrodes in electrolyte solutions are studied with Gaussian field theory. A response function with two Yukawa functions is used to capture the electrostatic response of the electrolyte solution, from which the modified response function in the planar symmetry is derived analytically. The modified response function is further used to evaluate the induced charge density and the electrostatic potential near an electrode. The Gaussian field theory, combined with a two-Yukawa response function, can reproduce the oscillatory decay behavior of the electric potentials in concentrated electrolyte solutions. When the exact sum rules for the bulk electrolyte solutions and the electric double layers are used as constraints to determine the parameters of the response function, the Gaussian field theory could at least partly capture the nonlinear response effect of the surface charge density. Comparison with results for a planar electrode with fixed surface charge densities from molecular simulations demonstrates the validity of Gaussian field theory.

Molecular dynamics of preferential adsorption in mixed alkali–halide electrolytes at graphene electrodes

Understanding the microscopic behavior of aqueous electrolyte solutions in contact with graphene and related carbon surfaces is important in electrochemical technologies, such as capacitive deionization or supercapacitors. In this work, we focus on preferential adsorption of ions in mixed alkali–halide electrolytes containing different fractions of Li+/Na+ or Li+/K+ and/or Na+/K+ cations with Cl− anions dissolved in water. We performed molecular dynamics simulations of the solutions in contact with both neutral and positively and negatively charged graphene surfaces under ambient conditions, using the effectively polarizable force field. The simulations show that large ions are often intuitively attracted to oppositely charged electrodes. In contrast, the adsorption behavior of small ions tends to be counterintuitive. In mixed-cation solutions, one of the cations always supports the adsorption of the other cation, while the other cation weakens the adsorption of the first cation. In mixed-cation solutions containing large and small cations simultaneously, adsorption of the larger cations varies dramatically with the electrode charge in an intuitive way, while adsorption of the smaller cations changes oppositely, i.e., in a counterintuitive way. For (Li/K)Cl mixed-cation solutions, these effects allow the control of Li+ adsorption by varying the electrode charge, whereas, for LiCl single-salt solutions, Li+ adsorption is nearly independent of the electrode charge. We rationalize this cation–cation lever effect as a result of a competition between three driving forces: (i) direct graphene–ion interactions, (ii) the strong tendency of the solutions to saturate the network of non-covalent intermolecular bonds, and (iii) the tendency to suppress local charge accumulation in any region larger than typical interparticle distances. We analyze the driving forces in detail using a general method for intermolecular bonding based on spatial distribution functions and different contributions to the total charge density profiles. The analysis helps to predict whether an ion is more affected by each of the three driving forces, depending on the strength of the ion solvation shells and the compatibility between the contributions of the charge density profiles due to the ion and water molecules. This approach is general and can also be applied to other solutions under different thermodynamic conditions.

Force density functional theory in- and out-of-equilibrium

When a fluid is subject to an external field, as is the case near an interface or under spatial confinement, then the density becomes spatially inhomogeneous. Although the one-body density provides much useful information, a higher level of resolution is provided by the two-body correlations. These give a statistical description of the internal microstructure of the fluid and enable calculation of the average interparticle force, which plays an essential role in determining both the equilibrium and dynamic properties of interacting fluids. We present a theoretical framework for the description of inhomogeneous (classical) many-body systems, based explicitly on the two-body correlation functions. By consideration of local Noether-invariance against spatial distortion of the system we demonstrate the fundamental status of the Yvon-Born-Green (YBG) equation as a local force-balance within the fluid. Using the inhomogeneous Ornstein-Zernike equation we show that the two-body correlations are density functionals and, thus, that the average interparticle force entering the YBG equation is also a functional of the one-body density. The force-based theory we develop provides an alternative to standard density functional theory for the study of inhomogeneous systems both in- and out-of-equilibrium. We compare force-based density profiles to the results of the standard potential-based (dynamical) density functional theory. In-equilibrium, we confirm both analytically and numerically that the standard approach yields profiles that are consistent with the compressibility pressure, whereas the force-density functional gives profiles consistent with the virial pressure. For both approaches we explicitly prove the hard-wall contact theorem that connects the value of the density profile at the hard-wall with the bulk pressure. The structure of the theory offers deep insights into the nature of correlation in dense and inhomogeneous systems.

Understanding the Electric Double-Layer Structure, Capacitance, and Charging Dynamics

Significant progress has been made in recent years in theoretical modeling of the electric double layer (EDL), a key concept in electrochemistry important for energy storage, electrocatalysis, and multitudes of other technological applications. However, major challenges remain in understanding the microscopic details of the electrochemical interface and charging mechanisms under realistic conditions. This review delves into theoretical methods to describe the equilibrium and dynamic responses of the EDL structure and capacitance for electrochemical systems commonly deployed for capacitive energy storage. Special emphasis is given to recent advances that intend to capture the nonclassical EDL behavior such as oscillatory ion distributions, polarization of nonmetallic electrodes, charge transfer, and various forms of phase transitions in the micropores of electrodes interfacing with an organic electrolyte or ionic liquid. This comprehensive analysis highlights theoretical insights into predictable relationships between materials characteristics and electrochemical performance and offers a perspective on opportunities for further development toward rational design and optimization of electrochemical systems.

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.