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

Impedance of nanocapacitors from molecular simulations to understand the dynamics of confined electrolytes

Nanoelectrochemical devices have become a promising candidate technology across various applications, including sensing and energy storage, and provide new platforms for studying fundamental properties of electrode/electrolyte interfaces. In this work, we employ constant-potential molecular dynamics simulations to investigate the impedance of gold-aqueous electrolyte nanocapacitors, exploiting a recently introduced fluctuation-dissipation relation. In particular, we relate the frequency-dependent impedance of these nanocapacitors to the complex conductivity of the bulk electrolyte in different regimes, and use this connection to design simple but accurate equivalent circuit models. We show that the electrode/electrolyte interfacial contribution is essentially capacitive and that the electrolyte response is bulk-like even when the interelectrode distance is only a few nanometers, provided that the latter is sufficiently large compared to the Debye screening length. We extensively compare our simulation results with spectroscopy experiments and predictions from analytical theories. In contrast to experiments, direct access in simulations to the ionic and solvent contributions to the polarization allows us to highlight their significant and persistent anticorrelation and to investigate the microscopic origin of the timescales observed in the impedance spectrum. This work opens avenues for the molecular interpretation of impedance measurements, and offers valuable contributions for future developments of accurate coarse-grained representations of confined electrolytes.

Charging and discharging a supercapacitor in molecular simulations

As the world moves more toward unpredictable renewable energy sources, better energy storage devices are required. Supercapacitors are a promising technology to meet the demand for short-term, high-power energy storage. Clearly, understanding their charging and discharging behaviors is essential to improving the technology. Molecular Dynamics (MD) simulations provide microscopic insights into the complex interplay between the dynamics of the ions in the electrolyte and the evolution of the charge distributions on the electrodes. Traditional MD simulations of (dis)charging supercapacitors impose a pre-determined evolving voltage difference between the electrodes, using the Constant Potential Method (CPM). Here, we present an alternative method that explicitly simulates the charge flow to and from the electrodes. For a disconnected capacitor, i.e., an open circuit, the charges are allowed to redistribute within each electrode while the sum charges on both electrodes remain constant. We demonstrate, for a model capacitor containing an aqueous salt solution, that this method recovers the charge-potential curve of CPM simulations. The equilibrium voltage fluctuations are related to the differential capacitance. We next simulate a closed circuit by introducing equations of motion for the sum charges, by explicitly accounting for the external circuit element(s). Charging and discharging of the model supercapacitor via a resistance proceed by double exponential processes, supplementing the usual time scale set by the electrolyte dynamics with a novel time scale set by the external circuit. Finally, we propose a simple equivalent circuit that reproduces the main characteristics of this supercapacitor.

A simple efficient algorithm for molecular simulations of constant potential electrodes

Increasingly, society requires high power, high energy storage devices for applications ranging from electric vehicles to buffers on the electric grid. Supercapacitors are a promising contribution to meeting these demands, though there still remain unsolved practical problems. Molecular dynamics simulations can shed light on the relevant molecular level processes in electric double layer capacitors, but these simulations are computationally very demanding. Our focus here is on the algorithmic complexity of the constant potential method (CPM), which uses dedicated electrostatics solvers to maintain a fixed potential difference between two conducting electrodes. We show how any standard electrostatics solver-capable of calculating the energies and forces on all atoms-can be used to implement CPM with a minimum of coding. As an example, we compare our generalized implementation of CPM, based on invocations of the particle-particle-particle-mesh routine of the Large-scale Atomic/Molecular Massively Parallel Simulator, with a traditional implementation based on a dedicated re-implementation of Ewald summation. Both methods yield comparable results on four test systems, with the former achieving a substantial gain in speed and improved scalability. The step from dedicated electrostatic solvers to generic routines is made possible by noting that CPM's traditional narrow Gaussian point-spread of atomic charges on the electrodes effectively endows point-like atoms with chemical hardness, i.e., an intra-atomic energy quadratic in the charge.

Nonmonotonic electrophoretic mobility of rodlike polyelectrolytes by multivalent coions in added salt

It is well established that when multivalent counterions or salts are added to a solution of highly charged polyelectrolytes (PEs), correlation effects can cause charge inversion of the PE, leading to electrophoretic mobility (EM) reversal. In this work, we use coarse-grained molecular-dynamics simulations to unravel the less understood effect of coion valency on EM reversal for rigid DNA-like PEs. We find that EM reversal induced by multivalent counterions is suppressed with increasing coion valency in the salt added and eventually vanishes. Further, we find that EM is enhanced at fixed low salt concentrations for salts with monovalent counterions when multivalent coions with increasing valency are introduced. However, increasing the salt concentration causes a crossover that leads to EM reversal which is enhanced by increasing coion valency at high salt concentration. Remarkably, this multivalent coion-induced EM reversal persists even for low values of PE linear charge densities where multivalent counterions alone cannot induce EM reversal. These results facilitate tuning PE-PE interactions and self-assembly with both coion and counterion valencies.

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.

Electrical noise in electrolytes: a theoretical perspective

Seemingly unrelated experiments such as electrolyte transport through nanotubes, nano-scale electrochemistry, NMR relaxometry and surface force balance measurements, all probe electrical fluctuations: of the electric current, the charge and polarization, the field gradient (for quadrupolar nuclei) and the coupled mass/charge densities. The fluctuations of such various observables arise from the same underlying microscopic dynamics of the ions and solvent molecules. In principle, the relevant length and time scales of these dynamics are encoded in the dynamic structure factors. However, modelling the latter for frequencies and wavevectors spanning many orders of magnitude remains a great challenge to interpret the experiments in terms of physical processes such as solvation dynamics, diffusion, electrostatic and hydrodynamic interactions between ions, interactions with solid surfaces, etc. Here, we highlight the central role of the charge-charge dynamic structure factor in the fluctuations of electrical observables in electrolytes and offer a unifying perspective over a variety of complementary experiments. We further analyze this quantity in the special case of an aqueous NaCl electrolyte, using simulations with explicit ions and an explicit or implicit solvent. We discuss the ability of the standard Poisson-Nernst-Planck theory to capture the simulation results, and how the predictions can be improved. We finally discuss the contributions of ions and water to the total charge fluctuations. This work illustrates an ongoing effort towards a comprehensive understanding of electrical fluctuations in bulk and confined electrolytes, in order to enable experimentalists to decipher the microscopic properties encoded in the measured electrical noise.

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.

Interactions between Rigid Polyelectrolytes Mediated by Ordering and Orientation of Multivalent Nonspherical Ions in Salt Solutions

Multivalent ions in solutions with polyelectrolytes (PEs) induce electrostatic correlations that can drastically change ion distributions around the PEs and their mutual interactions. Using coarse-grained molecular dynamics simulations, we show how in addition to valency, ion shape and concentration can be harnessed as tools to control rigid like-charged PE-PE interactions. We demonstrate a correlation between the orientational ordering of aspherical ions and how they mediate the effective PE-PE attraction induced by multivalency. The interaction type, strength, and range can thus be externally controlled in ionic solutions. Our results can be used as generic guidelines to tune the self-assembly of like-charged polyelectrolytes by variation of the characteristics of the ions.

The electrostatic effect and its role in promoting electrocatalytic reactions by specifically adsorbed anions

The adsorption of anions and its impact on electrocatalytic reactions are fundamental topics in electrocatalysis. Previous studies revealed that adsorbed anions display an overall poisoning effect in most cases. However, for a few reactions such as the hydrogen evolution reaction (HER), oxidation of small organic molecules (SOMs), and reduction of CO₂ and O₂, some specifically adsorbed anions can promote their reaction kinetics under certain conditions. The promotion effect is frequently attributed to the adsorbate induced modification of the nature of the active sites, the change of the adsorption configuration and free energy of the key reactive intermediate which consequently change the activation energy, the pre-exponential factor of the rate determining step etc. In this paper, we will give a mini review of the indispensable role of the classical double layer effect in enhancing the kinetics of electrocatalytic reactions by anion adsorption. The ubiquitous electrostatic interactions change both the potential distribution and the concentration distribution of ionic species across the electric double layer (EDL), which alters the electrochemical driving force and effective concentration of the reactants. The contribution to the overall kinetics is highlighted by taking HER, oxidation of SOMs, reduction of CO₂ and O₂, as examples.

Frequency and field-dependent response of confined electrolytes from Brownian dynamics simulations

Using Brownian dynamics simulations, we investigate the effects of confinement, adsorption on surfaces, and ion-ion interactions on the response of confined electrolyte solutions to oscillating electric fields in the direction perpendicular to the confining walls. Nonequilibrium simulations allows to characterize the transitions between linear and nonlinear regimes when varying the magnitude and frequency of the applied field, but the linear response, characterized by the frequency-dependent conductivity, is more efficiently predicted from the equilibrium current fluctuations. To that end, we (rederive and) use the Green-Kubo relation appropriate for overdamped dynamics, which differs from the standard one for Newtonian or underdamped Langevin dynamics. This expression highlights the contributions of the underlying Brownian fluctuations and of the interactions of the particles between them and with external potentials. Although already known in the literature, this relation has rarely been used to date, beyond the static limit to determine the effective diffusion coefficient or the DC conductivity. The frequency-dependent conductivity always decays from a bulk-like behavior at high frequency to a vanishing conductivity at low frequency due to the confinement of the charge carriers by the walls. We discuss the characteristic features of the crossover between the two regimes, most importantly how the crossover frequency depends on the confining distance and the salt concentration, and the fact that adsorption on the walls may lead to significant changes both at high and low frequencies. Conversely, our results illustrate the possibility to obtain information on diffusion between walls, charge relaxation, and adsorption by analyzing the frequency-dependent conductivity.

Frequency-Dependent Impedance of Nanocapacitors from Electrode Charge Fluctuations as a Probe of Electrolyte Dynamics

The frequency-dependent impedance is a fundamental property of electrical components. We show that it can be determined from the equilibrium dynamical fluctuations of the electrode charge in constant-potential molecular simulations, extending in particular a fluctuation-dissipation relation for the capacitance recovered in the low-frequency limit and provide an illustration on water-gold nanocapacitors. This Letter opens the way to the interpretation of electrochemical impedance measurements in terms of microscopic mechanisms, directly from the dynamics of the electrolyte, or indirectly via equivalent circuit models as in experiments.

Absence of anomalous underscreening in highly concentrated aqueous electrolytes confined between smooth silica surfaces

Recent surface forces apparatus experiments that measured the forces between two mica surfaces and a series of subsequent theoretical studies suggest the occurrence of universal underscreening in highly concentrated electrolyte solutions. We performed a set of systematic Atomic Force Spectroscopy measurements for aqueous salt solutions in a concentration range from 1 mM to 5 M using chloride salts of various alkali metals as well as mixed concentrated salt solutions (involving both mono- and divalent cations and anions), that mimic concentrated brines typically encountered in geological formations. Experiments were carried out using flat substrates and submicrometer-sized colloidal probes made of smooth oxidized silicon immersed in salt solutions at pH values of 6 and 9 and temperatures of 25 °C and 45 °C. While strong repulsive forces were observed for the smallest tip-sample separations, none of the conditions explored displayed any indication of anomalous long range electrostatic forces as reported for mica surfaces. Instead, forces are universally dominated by attractive van der Waals interactions at tip-sample separations of ≈2 nm and beyond for salt concentrations of 1 M and higher. Complementary calculations based on classical density functional theory for the primitive model support these experimental observations and display a consistent decrease in screening length with increasing ion concentration.

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.

A molecular perspective on induced charges on a metallic surface

Understanding the response of the surface of metallic solids to external electric field sources is crucial to characterize electrode-electrolyte interfaces. Continuum electrostatics offer a simple description of the induced charge density at the electrode surface. However, such a simple description does not take into account features related to the atomic structure of the solid and to the molecular nature of the solvent and of the dissolved ions. In order to illustrate such effects and assess the ability of continuum electrostatics to describe the induced charge distribution, we investigate the behavior of a gold electrode interacting with sodium or chloride ions fixed at various positions, in a vacuum or in water, using all-atom constant-potential classical molecular dynamics simulations. Our analysis highlights important similarities between the two approaches, especially under vacuum conditions and when the ion is sufficiently far from the surface, as well as some limitations of the continuum description, namely, neglecting the charges induced by the adsorbed solvent molecules and the screening effect of the solvent when the ion is close to the surface. While the detailed features of the charge distribution are system-specific, we expect some of our generic conclusions on the induced charge density to hold for other ions, solvents, and electrode surfaces. Beyond this particular case, the present study also illustrates the relevance of such molecular simulations to serve as a reference for the design of improved implicit solvent models of electrode-electrolyte interfaces.

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.