Simulations of ionic liquids confined by metal electrodes using periodic Green functions

Electrolytes in conducting nanopores: Revisiting constant charge and constant potential simulations

Simulating electrolyte-electrode systems poses challenges due to the need to account for the electrode's response to ion movements in order to maintain a constant electrode potential, which slows down the simulations. To circumvent this, computationally more efficient constant charge (CC) simulations are sometimes employed. However, the accuracy of CC simulations in capturing the behavior of electrolyte-electrode systems remains unclear, especially for microporous electrodes. Herein, we consider electrolyte-filled slit nanopores and systematically analyze the in-pore ion structure and diffusivity using CC and constant potential simulations. Our results indicate that CC simulations provide comparable pore occupancies at high bulk ion densities and for highly charged pores, but they fail to accurately describe the ion structure and dynamics, particularly in quasi-2D (single-layer) pores and at low ion densities. We attribute these results to the superionic state emerging in conducting nanoconfinement and its interplay with excluded volume interactions.

Charge regulation of nanoparticles in the presence of multivalent electrolytes

We explore the charge regulation (CR) of spherical nanoparticles immersed in an asymmetric electrolyte of a specified pH. Using a recently developed reactive canonical Monte Carlo (MC) simulation method, titration isotherms are obtained for suspensions containing monovalent, divalent, and trivalent coions. A theory based on the modified Poisson-Boltzmann approximation, which incorporates the electrostatic ion solvation free energy and discrete surface charge effects, is used to compare with the simulation results. A remarkably good agreement is found without any fitting parameters, both for the ion distributions and titration curves, suggesting that ionic correlations between coions and hydronium ions at the nanoparticle surface play only a minor role in determining the association equilibrium between hydroniums and the functional sites on the nanoparticle surface. On the other hand, if suspension contains multivalent counterions, we observe a large deviation between theory and simulations, showing that the electrostatic correlations between counterions and hydronium ions at the nanoparticle surface are very significant and must be properly taken into account to correctly describe CR for such solutions.

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.

Theory and Simulations of Ionic Liquids in Nanoconfinement

Room-temperature ionic liquids (RTILs) have exciting properties such as nonvolatility, large electrochemical windows, and remarkable variety, drawing much interest in energy storage, gating, electrocatalysis, tunable lubrication, and other applications. Confined RTILs appear in various situations, for instance, in pores of nanostructured electrodes of supercapacitors and batteries, as such electrodes increase the contact area with RTILs and enhance the total capacitance and stored energy, between crossed cylinders in surface force balance experiments, between a tip and a sample in atomic force microscopy, and between sliding surfaces in tribology experiments, where RTILs act as lubricants. The properties and functioning of RTILs in confinement, especially nanoconfinement, result in fascinating structural and dynamic phenomena, including layering, overscreening and crowding, nanoscale capillary freezing, quantized and electrotunable friction, and superionic state. This review offers a comprehensive analysis of the fundamental physical phenomena controlling the properties of such systems and the current state-of-the-art theoretical and simulation approaches developed for their description. We discuss these approaches sequentially by increasing atomistic complexity, paying particular attention to new physical phenomena emerging in nanoscale confinement. This review covers theoretical models, most of which are based on mapping the problems on pertinent statistical mechanics models with exact analytical solutions, allowing systematic analysis and new physical insights to develop more easily. We also describe a classical density functional theory, which offers a reliable and computationally inexpensive tool to account for some microscopic details and correlations that simplified models often fail to consider. Molecular simulations play a vital role in studying confined ionic liquids, enabling deep microscopic insights otherwise unavailable to researchers. We describe the basics of various simulation approaches and discuss their challenges and applicability to specific problems, focusing on RTIL structure in cylindrical and slit confinement and how it relates to friction and capacitive and dynamic properties of confined ions.

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.

Phase Transitions and Electrochemical Properties of Ionic Liquids and Ionic Liquid-Solvent Mixtures

Recent advances in studies of ionic liquids (IL) and ionic liquid-solvent mixtures are reviewed. Selected experimental, simulation, and theoretical results for electrochemical, thermodynamical, and structural properties of IL and IL-solvent mixtures are described. Special attention is paid to phenomena that are not predicted by the classical theories of the electrical double layer or disagree strongly with these theories. We focus on structural properties, especially on distribution of ions near electrodes, on electrical double layer capacitance, on effects of confinement, including decay length of a dissjoining pressure between confinig plates, and on demixing phase transition. In particular, effects of the demixing phase transition on electrochemical properties of ionic liquid-solvent mixtures for different degrees of confinement are presented.

Molecular Simulation of Electrode-Solution Interfaces

Many key industrial processes, from electricity production, conversion, and storage to electrocatalysis or electrochemistry in general, rely on physical mechanisms occurring at the interface between a metallic electrode and an electrolyte solution, summarized by the concept of an electric double layer, with the accumulation/depletion of electrons on the metal side and of ions on the liquid side. While electrostatic interactions play an essential role in the structure, thermodynamics, dynamics, and reactivity of electrode-electrolyte interfaces, these properties also crucially depend on the nature of the ions and solvent, as well as that of the metal itself. Such interfaces pose many challenges for modeling because they are a place where quantum chemistry meets statistical physics. In the present review, we explore the recent advances in the description and understanding of electrode-electrolyte interfaces with classical molecular simulations, with a focus on planar interfaces and solvent-based liquids, from pure solvent to water-in-salt electrolytes.

Superionic Liquids in Conducting Nanoslits: Insights from Theory and Simulations

Mapping the theory of charging supercapacitors with nanostructured electrodes on known lattice models of statistical physics is an interesting task, aimed at revealing generic features of capacitive energy storage in such systems. The main advantage of this approach is the possibility to obtain analytical solutions that allow new physical insights to be more easily developed. But how general the predictions of such theories could be? How sensitive are they to the choice of the lattice? Herein, we address these questions in relation to our previous description of such systems using the Bethe-lattice approach and Monte Carlo simulations. Remarkably, we find a surprisingly good agreement between the analytical theory and simulations. In addition, we reveal a striking correlation between the ability to store energy and ion ordering inside a pore, suggesting that such ordering can be beneficial for energy storage.

Charge fluctuations from molecular simulations in the constant-potential ensemble

Statistical mechanics of constant-potential molecular simulations yields a new fluctuation–dissipation relation for the differential capacitance.

Feeling Your Neighbors across the Walls: How Interpore Ionic Interactions Affect Capacitive Energy Storage

Progress in low-dimensional carbon materials has intensified research on supercapacitors with nanostructured/nanoporous electrodes. The theoretical and simulation work so far has focused on charging single nanopores or nanoporous networks and the effects due to ionic interactions inside the pores, while the effect of interpore ion-ion correlations has received less attention. Herein, we study how the interactions between the ions in the neighboring pores across the pore walls affect capacitive energy storage. We develop a simple lattice model for the ions in a stack of parallel-aligned nanotubes, solve it by using the perturbation and "semi-mean-field" theories, and test the results by Monte Carlo simulations. We demonstrate that the interpore ionic interactions can have a profound effect on charge storage; in particular, such interactions can enhance or diminish the stored energy density, depending on the sign of like-charge interactions. We also find that charging can proceed either continuously or via a phase transition. Our results call for more detailed investigations of the properties of carbon pore walls and suggest that tuning their electrostatic response may be promising for the rational design of an optimal supercapacitor.

Like-Charge Attraction between Metal Nanoparticles in a 1∶1 Electrolyte Solution

We calculate the force between two spherical metal nanoparticles of charge Q_{1} and Q_{2} in a dilute 1∶1 electrolyte solution. Numerically solving the nonlinear Poisson-Boltzmann equation, we find that metal nanoparticles with the same sign of charge can attract one another. This is fundamentally different from what is found for like-charged, nonpolarizable, colloidal particles, the two-body interaction potential for which is always repulsive inside a dilute 1∶1 electrolyte. Furthermore, the existence of like-charge attraction between spherical metal nanoparticles is even more surprising in view of the result that such attraction is impossible between parallel metal slabs, showing the fundamental importance of curvature. To overcome a slow convergence of the numerical solution of the full nonlinear Poisson-Boltzmann equation, we developed a modified Derjaguin approximation which allows us to accurately and rapidly calculate the interaction potential between two metal nanoparticles or between a metal nanoparticle and a phospholipid membrane.

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.

Lattice model of ionic liquid confined by metal electrodes

We study, using Monte Carlo simulations, the density profiles and differential capacitance of ionic liquids confined by metal electrodes. To compute the electrostatic energy, we use the recently developed approach based on periodic Green's functions. The method also allows us to easily calculate the induced charge on the electrodes permitting an efficient implementation of simulations in a constant electrostatic potential ensemble. To speed up the simulations further, we model the ionic liquid as a lattice Coulomb gas and precalculate the interaction potential between the ions. We show that the lattice model captures the transition between camel-shaped and bell-shaped capacitance curves-the latter characteristic of ionic liquids (strong coupling limit) and the former of electrolytes (weak coupling). We observe the appearance of a second peak in the differential capacitance at ≈0.5 V for 2:1 ionic liquids, as the packing fraction is increased. Finally, we show that ionic size asymmetry decreases substantially the capacitance maximum, when all other parameters are kept fixed.

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.

Simulations of Coulomb systems confined by polarizable surfaces using periodic Green functions

We present an efficient approach for simulating Coulomb systems confined by planar polarizable surfaces. The method is based on the solution of the Poisson equation using periodic Green functions. It is shown that the electrostatic energy arising from the surface polarization can be decoupled from the energy due to the direct Coulomb interaction between the ions. This allows us to combine an efficient Ewald summation method, or any other fast method for summing over the replicas, with the polarization contribution calculated using Green function techniques. We apply the method to calculate density profiles of ions confined between the charged dielectric and metal surfaces.

An experimental/theoretical method to measure the capacitive compactness of an aqueous electrolyte surrounding a spherical charged colloid

The capacitive compactness has been introduced very recently [G. I. Guerrero-García et al., Phys. Chem. Chem. Phys. 20, 262-275 (2018)] as a robust and accurate measure to quantify the thickness, or spatial extension, of the electrical double layer next to either an infinite charged electrode or a spherical macroion. We propose here an experimental/theoretical scheme to determine the capacitive compactness of a spherical electrical double layer that relies on the calculation of the electrokinetic charge and the associated mean electrostatic potential at the macroparticle's surface. This is achieved by numerically solving the non-linear Poisson-Boltzmann equation of point ions around a colloidal sphere and matching the corresponding theoretical mobility, predicted by the O'Brien and White theory [J. Chem. Soc., Faraday Trans. 2 74, 1607-1626 (1978)], with experimental measurements of the electrophoretic mobility under the same conditions. This novel method is used to calculate the capacitive compactness of NaCl and CaCl₂ electrolytes surrounding a negatively charged polystyrene particle as a function of the salt concentration.

Quantifying the thickness of the electrical double layer neutralizing a planar electrode: the capacitive compactness

The spatial extension of the ionic cloud neutralizing a charged colloid or an electrode is usually characterized by the Debye length associated with the supporting charged fluid in the bulk. This spatial length arises naturally in the linear Poisson-Boltzmann theory of point charges, which is the cornerstone of the widely used Derjaguin-Landau-Verwey-Overbeek formalism describing the colloidal stability of electrified macroparticles. By definition, the Debye length is independent of important physical features of charged solutions such as the colloidal charge, electrostatic ion correlations, ionic excluded volume effects, or specific short-range interactions, just to mention a few. In order to include consistently these features to describe more accurately the thickness of the electrical double layer of an inhomogeneous charged fluid in planar geometry, we propose here the use of the capacitive compactness concept as a generalization of the compactness of the spherical electrical double layer around a small macroion (González-Tovar et al., J. Chem. Phys. 2004, 120, 9782). To exemplify the usefulness of the capacitive compactness to characterize strongly coupled charged fluids in external electric fields, we use integral equations theory and Monte Carlo simulations to analyze the electrical properties of a model molten salt near a planar electrode. In particular, we study the electrode's charge neutralization, and the maximum inversion of the net charge per unit area of the electrode-molten salt system as a function of the ionic concentration, and the electrode's charge. The behaviour of the associated capacitive compactness is interpreted in terms of the charge neutralization capacity of the highly correlated charged fluid, which evidences a shrinking/expansion of the electrical double layer at a microscopic level. The capacitive compactness and its first two derivatives are expressed in terms of experimentally measurable macroscopic properties such as the differential and integral capacity, the electrode's surface charge density, and the mean electrostatic potential at the electrode's surface.