Density functional theory for differential capacitance of planar electric double layers in ionic liquids

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.

Integrative studies of ionic liquid interface layers: bridging experiments, theoretical models and simulations

Ionic liquids (ILs) are a class of salts existing in the liquid state below 100 °C, possessing low volatility, high thermal stability as well as many highly attractive solvent and electrochemical capabilities, etc., making them highly tunable for a great variety of applications, such as lubricants, electrolytes, and soft functional materials. In many applications, ILs are first either physi- or chemisorbed on a solid surface to successively create more functional materials. The functions of ILs at solid surfaces can differ considerably from those of bulk ILs, mainly due to distinct interfacial layers with tunable structures resulting in new ionic liquid interface layer properties and enhanced performance. Due to an almost infinite number of possible combinations among the cations and anions to form ILs, the diversity of various solid surfaces, as well as different external conditions and stimuli, a detailed molecular-level understanding of their structure-property relationship is of utmost significance for a judicious design of IL-solid interfaces with appropriate properties for task-specific applications. Many experimental techniques, such as atomic force microscopy, surface force apparatus, and so on, have been used for studying the ion structuring of the IL interface layer. Molecular Dynamics simulations have been widely used to investigate the microscopic behavior of the IL interface layer. To interpret and clarify the IL structure and dynamics as well as to predict their properties, it is always beneficial to combine both experiments and simulations as close as possible. In another theoretical model development to bridge the structure and properties of the IL interface layer with performance, thermodynamic prediction & property modeling has been demonstrated as an effective tool to add the properties and function of the studied nanomaterials. Herein, we present recent findings from applying the multiscale triangle "experiment-simulation-thermodynamic modeling" in the studies of ion structuring of ILs in the vicinity of solid surfaces, as well as how it qualitatively and quantitatively correlates to the overall ILs properties, performance, and function. We introduce the most common techniques behind "experiment-simulation-thermodynamic modeling" and how they are applied for studying the IL interface layer structuring, and we highlight the possibilities of the IL interface layer structuring in applications such as lubrication and energy storage.

Differential Capacitance of Ionic Liquid Confined between Metallic Interfaces

We present here a detailed analysis of the electric double layer at the gold electrode/[BMIM][BF4] interface using a polarizable model for the electrode, based on our recent approach to include image charges [Geada et al. Nat. Commun. 2018, 9, 716]. A double bell (camel) shape is obtained for the differential capacitance, where the inclusion of metal polarization allows for a higher density of ions in the double layer, particularly around the maxima, thereby increasing the capacitance. The charging mechanism differs for the positive and negative electrodes, with counterion adsorption prevailing at the anode and co-ion desorption prevailing at the cathode. The charging mechanism is predominantly governed by the BF4 anions, serving as counterions and co-ions at the anode and cathode, respectively. Within the considered range of potentials, only minor changes are observed in the dynamical properties, specifically in the diffusion coefficients. Notably, it is interesting to observe that bulk properties are restored at a shorter distance from the gold surface in the case of the anode compared to the cathode.

Unraveling Flow Effect on Capacitive Energy Extraction from Salinity Gradients

The harvesting of salinity gradient energy through a capacitive double-layer expansion (CDLE) technique is directly associated with ion adsorption and desorption in electrodes. Herein, we show that energy extraction can be modulated by regulating ion adsorption/desorption through water flow. The flow effects on the output energy, capacitance, and energy density under practical conditions are systematically investigated from a theoretical perspective, upon which the optimal operating condition is identified for energy extraction. We demonstrate that the net charge accumulation displays a negative correlation with the water flow velocity and so does the surface charge density, and this causes a nontrivial variation in the magnitude of output energy when water flows are introduced. When high water flows are introduced in both the charging and discharging processes, the energy extraction can be significantly reduced by 47.69-49.32%. However, when a high flow is solely exerted in the discharging process, the energy extraction can be enhanced by 12.94-14.49% even at low operation voltages. This study not only offers a comprehensive understanding of the microscopic mechanisms of surface-engineered energy extraction with water flows but also provides a novel direction for energy extraction enhancement.

Microscopic Simulations of Electrochemical Double-Layer Capacitors

Electrochemical double-layer capacitors (EDLCs) are devices allowing the storage or production of electricity. They function through the adsorption of ions from an electrolyte on high-surface-area electrodes and are characterized by short charging/discharging times and long cycle-life compared to batteries. Microscopic simulations are now widely used to characterize the structural, dynamical, and adsorption properties of these devices, complementing electrochemical experiments and in situ spectroscopic analyses. In this review, we discuss the main families of simulation methods that have been developed and their application to the main family of EDLCs, which include nanoporous carbon electrodes. We focus on the adsorption of organic ions for electricity storage applications as well as aqueous systems in the context of blue energy harvesting and desalination. We finally provide perspectives for further improvement of the predictive power of simulations, in particular for future devices with complex electrode compositions.

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.

Structural Forces in Ionic Liquids: The Role of Ionic Size Asymmetry

Ionic liquids (ILs) are charged fluids composed of anions and cations of different size and shape. The ordering of charge and density in ILs confined between charged interfaces underlies numerous applications of IL electrolytes. Here, we analyze the screening behavior and the resulting structural forces of a representative IL confined between two charge-varied plates. Using both molecular dynamics simulations and a continuum theory, we contrast the screening features of a more-realistic asymmetric system and a less-realistic symmetric one. The ionic size asymmetry plays a nontrivial role in charge screening, affecting both the ionic density profiles and the disjoining pressure distance dependence. Ionic systems with size asymmetry are stronger coupled systems, and this manifests itself both in their response to the electrode polarization and spontaneous structure formation at the interface. Analytical expressions for decay lengths of the disjoining pressure are obtained in agreement with the pressure profiles computed from molecular dynamics simulations.

Double-Edged Sword of Ion-Size Asymmetry in Energy Storage of Supercapacitors

The advanced supercapacitor is of great significance for renewable energy storage. Achieving its high energy and high power densities remains a huge challenge. Herein, the contribution of ion-size asymmetry to the charging behavior of a supercapacitor is systematically studied using time-dependent density functional theory (TDDFT). We track the time evolution of the ionic microstructure inside the porous electrode and its reservoir and reveal a kinetic charge inversion in the asymmetrical ion-size cases. Compared with the symmetrical ion-size case, we find that the ion-size asymmetry has a double-edged sword effect on the energy storage of a supercapacitor: it accelerates the charging process yet reduces the differential capacitance. Additionally, the energy density and power density can simultaneously increase in the asymmetrical cases, which provides important insights toward the experimental design of supercapacitors with high energy and high power densities.

Electrochemistry meets polymer physics: polymerized ionic liquids on an electrified electrode

Polymeric ionic liquids are emerging polyelectrolyte materials for modern electrochemical applications. In this paper, we propose a self-consistent field theory of a polymeric ionic liquid on a charged conductive electrode. Taking into account the conformational entropy of rather long polymerized cations within the Lifshitz theory and electrostatic and excluded volume interactions of ionic species within the mean-field approximation, we obtain a system of self-consistent field equations for the local electrostatic potential and average concentrations of monomeric units and counterions. We solve these equations in the linear approximation for the cases of a point-like charge and a flat infinite uniformly charged electrode immersed in a polymeric ionic liquid and derive analytical expressions for local ionic concentrations and electrostatic potential, and derive an analytical expression for the linear differential capacitance of the electric double layer. We also find a numerical solution to the self-consistent field equations for two types of boundary conditions for the local polymer concentration on the electrode, corresponding to the cases of the specific adsorption absence (indifferent surface) and strong short-range repulsion of the monomeric units near the charged surface (hard wall case). For both cases, we investigate the behavior of differential capacitance as a function of applied voltage for a pure polymeric ionic liquid and a polymeric ionic liquid dissolved in a polar organic solvent. We observe that the differential capacitance profile shape is strongly sensitive to the adopted boundary condition for the local polymer concentration on the electrode.

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.

Capacitance and Structure of Electric Double Layers: Comparing Brownian Dynamics and Classical Density Functional Theory

We present a study of the structure and differential capacitance of electric double layers of aqueous electrolytes. We consider electric double layer capacitors (EDLC) composed of spherical cations and anions in a dielectric continuum confined between a planar cathode and anode. The model system includes steric as well as Coulombic ion-ion and ion-electrode interactions. We compare results of computationally expensive, but “exact” , Brownian Dynamics (BD) simulations with approximate, but cheap, calculations based on classical Density Functional Theory (DFT). Excellent overall agreement is found for a large set of system parameters, including variations in concentration, ionic size- and valency-asymmetries, applied voltages and electrode separation, provided the differences between the canonical ensemble of the BD simulations and the grand-canonical ensemble of DFT are properly taken into account. In particular, a careful distinction is made between the differential capacitance $$C_N$$ C N at fixed number of ions and $$C_\mu $$ C μ at fixed ionic chemical potential. Furthermore, we derive and exploit their thermodynamic relations. In the future these relations will also be useful for comparing and contrasting experimental data with theories for supercapactitors and other systems. The quantitative agreement between simulation and theory indicates that the presented DFT is capable of accounting accurately for coupled Coulombic and packing effects. Hence it is a promising candidate to cheaply study room temperature ionic liquids at much lower dielectric constants than that of water.

Surface Charge Density in Electrical Double Layer Capacitors with Nanoscale Cathode-Anode Separation

Using a dynamic density functional theory, we study the charging dynamics, the final equilibrium structure, and the energy storage in an electrical double layer capacitor with nanoscale cathode-anode separation in a slit geometry. We derive a simple expression for the surface charge density that naturally separates the effects of the charge polarization due to the ions from those due to the polarization of the dielectric medium and allows a more intuitive understanding of how the ion distribution within the cell affects the surface charge density. We find that charge neutrality in the half-cell does not hold during the dynamic charging process for any cathode-anode separation, and also does not hold at the final equilibrium state for small separations. Therefore, the charge accumulation in the half-cell in general does not equal the surface charge density. The relationships between the surface charge density and the charge accumulation within the half-cell are systematically investigated by tuning the electrolyte concentration, cathode-anode separation, and applied voltage. For high electrolyte concentrations, we observe charge inversion at which the charge accumulation exceeds the surface charge at special values of the separation. In addition, we find that the energy density has a maximum at intermediate electrolyte concentrations for a high applied voltage.

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.

Dynamic Adsorption of Ions into Like-Charged Nanospace: A Dynamic Density Functional Theory Study

The adsorption processes of ions into charged nanospace are associated with many practical applications. Whereas a large number of microporous materials have been prepared toward efficient adsorption of ions from solutions, theoretical models that allow for capturing the characteristics of ion dynamic adsorption into like-charged nanopores are still few. The difficulty originates from the overlapping of electric potentials inside the pores. Herein, a theoretical model is proposed by incorporating dynamic density functional theory with modified Poisson equation for investigating the dynamic adsorption of ions into like-charged nanoslits. This model is rationalized by comparing the theoretical predictions with corresponding simulation results. Afterward, by analyzing the adsorption dynamics, we show that the overlapping effect is associated with the pore size, ion bulk concentration, and surface charge density, and it plays a dominant role in the coupling between the total adsorption amount of ions and total adsorption time. Specifically, with weak overlapping effect, the total adsorption amount is intuitively proportional to the total adsorption time; however, when the overlapping effect is strong, the total adsorption amount may be inversely proportional to the total adsorption time, indicating that both high adsorption amount and short adsorption time can be achieved simultaneously. This work provides a meaningful insight toward the rational design and optimization of microporous materials for efficient ion adsorption.

Does capillary evaporation limit the accessibility of nonaqueous electrolytes to the ultrasmall pores of carbon electrodes?

Porous carbons have been widely utilized as electrode materials for capacitive energy storage. Whereas the importance of pore size and geometry on the device performance has been well recognized, little guidance is available for identification of carbon materials with ideal porous structures. In this work, we study the phase behavior of ionic fluids in slit pores using the classical density functional theory. Within the framework of the restricted primitive model for nonaqueous electrolytes, we demonstrate that the accessibility of micropores depends not only on the ionic diameters (or desolvation) but also on their wetting behavior intrinsically related to the vapor-liquid or liquid-liquid phase separation of the bulk ionic systems. Narrowing the pore size from several tens of nanometers to subnanometers may lead to a drastic reduction in the capacitance due to capillary evaporation. The wettability of micropores deteriorates as the pore size is reduced but can be noticeably improved by raising the surface electrical potential. The theoretical results provide fresh insights into the properties of confined ionic systems beyond electric double layer models commonly employed for rational design/selection of electrolytes and electrode materials.

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.

A multi-scale approach to describe electrical impulses propagating along actin filaments in both intracellular and in vitro conditions

An accurate and efficient characterization of the polyelectrolyte properties for cytoskeleton filaments are key to the molecular understanding of electrical signal propagation, bundle and network formation, as well as their potential nanotechnological applications. In this article, we introduce an innovative multi-scale approach able to account for the atomistic details of a protein molecular structure, its biological environment, and their impact on electrical impulses propagating along wild type actin filaments. The formulation includes non-trivial contributions to the ionic electrical conductivity and capacitance coming from the diffuse part of the electrical double layer of G-actins. We utilize this monomer characterization in a non-linear inhomogeneous transmission line prototype model to account for the monomer–monomer interactions, dissipation and damping perturbations along the filament length. A novel, simple, accurate, approximate analytic expression has been obtained for the transmission line model. Our results reveal the propagation of electrical signal impulses in the form of solitons for the range of voltage stimulus and electrolyte solutions typically present for intracellular and in vitro conditions. The approach predicts a lower electrical conductivity with higher linear capacitance and non-linear accumulation of charge for intracellular conditions. Our results show a significant influence of the voltage input on the electrical impulse shape, attenuation and kern propagation velocity. The filament is able to sustain the soliton propagation at almost constant kern velocity for the in vitro condition, whereas the intracellular condition displays a remarkable deceleration. Additionally, the solitons are narrower and travel faster at higher voltage input. As a unique feature, this multi-scale theory is able to account for molecular structure conformation (mutation) and biological environment (protonations/deprotonations) changes often present in pathological conditions. It is also applicable to other highly charged rod-like polyelectrolytes with relevance in biomedicine and biophysics.

An innovative analytic solution accounting for the molecular structure, its biological environment, and their impact on electrical impulses along microfilaments.

Including diverging electrostatic potential in 3D-RISM theory: The charged wall case

Although three-dimensional site-site molecular integral equations of liquids are a powerful tool of the modern theoretical chemistry, their applications to the problem of characterizing the electrical double layer originating at the solid-liquid interface with a macroscopic substrate are severely limited by the fact that an infinitely extended charged plane generates a divergent electrostatic potential. Such potentials cannot be treated within the standard 3D-Reference Interaction Site Model equation solution framework since it leads to functions that are not Fourier transformable. In this paper, we apply a renormalization procedure to overcome this obstacle. We then check the validity and numerical accuracy of the proposed computational scheme on the prototypical gold (111) surface in contact with water/alkali chloride solution. We observe that despite the proposed method requires, to achieve converged charge densities, a higher spatial resolution than that suited to the estimation of biomolecular solvation with either 3D-RISM or continuum electrostatics approaches, it still is computationally efficient. Introducing the electrostatic potential of an infinite wall, which is periodic in 2 dimensions, we avoid edge effects, permit a robust integration of Poisson's equation, and obtain the 3D electrostatic potential profile for the first time in such calculations. We show that the potential within the electrical double layer presents oscillations which are not grasped by the Debye-Hückel and Gouy-Chapman theories. This electrostatic potential deviates from its average of up to 1-2 V at small distances from the substrate along the lateral directions. Applications of this theoretical development are relevant, for example, for liquid scanning tunneling microscopy imaging.

Theoretical study of the effect of π^{+}-π^{+} association in imidazolium ionic liquids at charged interfaces

We develop an extended classical density-functional theory to describe clustering of imidazolium-based cations into linear chains, driven by π-π stacking. We find that the associating system displays a similar short-ranged structure to the completely dissociated fluid. We also construct a restricted primitive model for associating ionic species in an RTIL+solvent mixture. The double-layer formed in these systems exhibits strong overscreening by the cation chains, as expected. Significantly enhanced capacitance is also observed for the case where counterions are the associating species. The established density-functional method can be also used to describe polydisperse polyelectrolyte models.

Structure of electric double layers in capacitive systems and to what extent (classical) density functional theory describes it

Ongoing scientific interest is aimed at the properties and structure of electric double layers (EDLs), which are crucial for capacitive energy storage, water treatment, and energy harvesting technologies like supercapacitors, desalination devices, blue engines, and thermocapacitive heat-to-current converters. A promising tool to describe their physics on a microscopic level is (classical) density functional theory (DFT), which can be applied in order to analyze pair correlations and charge ordering in the primitive model of charged hard spheres. This simple model captures the main properties of ionic liquids and solutions and it predicts many of the phenomena that occur in EDLs. The latter often lead to anomalous response in the differential capacitance of EDLs. This work constructively reviews the powerful theoretical framework of DFT and its recent developments regarding the description of EDLs. It explains to what extent current approaches in DFT describe structural ordering and in-plane transitions in EDLs, which occur when the corresponding electrodes are charged. Further, the review briefly summarizes the history of modeling EDLs, presents applications, and points out limitations and strengths in present theoretical approaches. It concludes that DFT as a sophisticated microscopic theory for ionic systems is expecting a challenging but promising future in both fundamental research and applications in supercapacitive technologies.

Lattice Model of an Ionic Liquid at an Electrified Interface

We study ionic liquids interacting with electrified interfaces. The ionic fluid is modeled as a Coulomb lattice gas. We compare the ionic density profiles calculated using a popular modified Poisson-Boltzmann equation with the explicit Monte Carlo simulations. The modified Poisson-Boltzmann theory fails to capture the structural features of the double layer and is also unable to correctly predict the ionic density at the electrified interface. The lattice Monte Carlo simulations qualitatively capture the coarse-grained structure of the double layer in the continuum. We propose a convolution relation that semiquantitatively relates the ionic density profiles of a continuum ionic liquid and its lattice counterpart near an electrified interface.

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.

Control of Selective Ion Transfer across Liquid-Liquid Interfaces: A Rectifying Heterojunction Based on Immiscible Electrolytes

The current rectification displayed by solid-state p-n semiconductor diodes relies on the abundance of electrons and holes near the interface between the p-n junction. In analogy to this electronic device, we propose here the construction of a purely ionic liquid-state electric rectifying heterojunction displaying an excess of monovalent cations and anions near the interface between two immiscible solvents with different dielectric properties. This system does not need any physical membrane or material barrier to show preferential ion transfer but relies on the ionic solvation energy between the two immiscible solvents. We construct a simple device, based on an oil/water interface, displaying an asymmetric behavior of the electric current as a function of the polarity of an applied electric field. This device also exhibits a region of negative differential conductivity, analogous to that observed in brain and heart cells via voltage clamp techniques. Computer simulations and mean field theory calculations for a model of this system show that the application of an external electric field is able to control the bulk concentrations of the ionic species in the immiscible liquids in a manner that is asymmetric with respect to the polarity or direction of the applied electric field. These properties make possible to enhance or suppress selective ion transport at liquid-liquid interfaces with the application of an external electric field or electrostatic potential, mimicking the function of biological ion channels, thus creating opportunities for varied applications.