Evaluation of molecular dynamics simulation methods for ionic liquid electric double layers

Temperature-dependent differential capacitance of an ionic liquid-graphene-based supercapacitor

One of the critical factors affecting the performance of supercapacitors is thermal management. The design of supercapacitors that operate across a broad temperature range and at high charge/discharge rates necessitates understanding the correlation of the molecular characteristics of the device (such as interfacial structure and inter-ionic and ion-electrode interactions) with its macroscopic properties. In this study, we use molecular dynamics (MD) simulations to investigate the influence of Joule heating on the structure and dynamics of the ionic liquid (IL)/graphite-based supercapacitors. The temperature-dependent electrical double layer (EDL) and differential capacitance-potential (CD-V) curves of two different ([Bmim][BF4] and [Bmim][PF6]) IL-graphene pairs were studied under various thermal gradients. For the [Bmim][BF4] system, the differential capacitance curves transition from 'U' to bell shape under an applied thermal gradient (∇T) in the range from 3.3 K nm-1 to 16.7 K nm-1. Whereas in [Bmim][PF6], we find a positive dependence of differential capacitance with ∇T with a U-shaped CD-V curve. We examine changes in the EDL structure and screening potential (ϕ(z)) as a function of ∇T and correlate them with the trends observed in the CD-V curve. The identified correlation between the interfacial charge density and differential capacitance with thermal gradient would be helpful for the molecular design of the IL-electrode interface in supercapacitors or other chemical engineering applications.

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.

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.

Innermost Ion Association Configuration Is a Key Structural Descriptor of Ionic Liquids at Electrified Interfaces

The structure of electric double layers (EDLs) is crucial for all types of electrochemical processes. While in dilute solutions EDL structure can be approximately treated within the Gouy-Chapman-Stern regime, in highly ionic electrolytes the description of EDL has been largely elusive. Here we study the EDL structure of an ionic liquid on a series of crystalline electrodes. Through molecular dynamics (MD) simulations, we observe strong intermolecular interaction among cations and anions and propose that the cation-anion association structure at the innermost layer is a key descriptor of the EDL. Using our recently developed electrochemical 3D atomic force microscopy (EC-3D-AFM) technique, we confirm the theoretical prediction and further find that the width of the first EDL is an experimental gauge of the ion association structure in that layer. We expect such ion association descriptors to be broadly applicable to a large range of highly ionic electrolytes on various electrode surfaces.

ELECTRODE: An electrochemistry package for atomistic simulations

Constant potential methods (CPM) enable computationally efficient simulations of the solid-liquid interface at conducting electrodes in molecular dynamics (MD). They have been successfully used, for example, to realistically model the behavior of ionic liquids or water-in-salt electrolytes in supercapacitors and batteries. The CPM models conductive electrodes by updating charges of individual electrode atoms according to the applied electric potential and the (time-dependent) local electrolyte structure. Here we present a feature-rich CPM implementation, called ELECTRODE, for the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), which includes a constrained charge method and a thermo-potentiostat. The ELECTRODE package also contains a finite-field approach, multiple corrections for non-periodic boundary conditions of the particle-particle particle-mesh solver, and a Thomas-Fermi model for using non-ideal metals as electrodes. We demonstrate the capabilities of this implementation for a parallel-plate electrical double-layer capacitor, for which we have investigated the charging times with the different implemented methods and found an interesting relationship between water and ionic dipole relaxations. To prove the validity of the one-dimensional correction for the long-range electrostatics, we estimated the vacuum capacitance of two co-axial carbon nanotubes and compared it to structureless cylinders, for which an analytical expression exists. In summary, the ELECTRODE package enables efficient electrochemical simulations using state-of-the-art methods, allowing one to simulate even heterogeneous electrodes. Moreover, it allows unveiling more rigorously how electrode curvature affects the capacitance with the one-dimensional correction.

Unified polarizable electrode models for open and closed circuits: Revisiting the effects of electrode polarization and different circuit conditions on electrode-electrolyte interfaces

A precise understanding of the interfacial structure and dynamics is essential for the optimal design of various electrochemical devices. Herein, we propose a method for classical molecular dynamics simulations to deal with electrochemical interfaces with polarizable electrodes under the open circuit condition. Less attention has been paid to electrochemical circuit conditions in computation despite being often essential for a proper assessment, especially comparison between different models. The present method is based on the chemical potential equalization principle, as is a method developed previously to deal with systems under the closed circuit condition. These two methods can be interconverted through the Legendre transformation, so that the difference in the circuit conditions can be compared on the same footing. Furthermore, the electrode polarization effect can be correctly studied by comparing the present method with the conventional simulations with the electrodes represented by fixed charges, since both of the methods describe systems under the open circuit condition. The method is applied to a parallel-plate capacitor composed of platinum electrodes and an aqueous electrolyte solution. The electrode polarization effects have an impact on the interfacial structure of the electrolyte solution. We found that the difference in the circuit conditions significantly affects the dynamics of the electrolyte solution. The electric field at the charged electrode surface is poorly screened by the nonequilibrium solution structure in the open circuit condition, which accelerates the motion of the electrolyte solution.

Fully periodic, computationally efficient constant potential molecular dynamics simulations of ionic liquid supercapacitors

Molecular dynamics (MD) simulations of complex electrochemical systems, such as ionic liquid supercapacitors, are increasingly including the constant potential method (CPM) to model conductive electrodes at a specified potential difference, but the inclusion of CPM can be computationally expensive. We demonstrate the computational savings available in CPM MD simulations of ionic liquid supercapacitors when the usual non-periodic slab geometry is replaced with fully periodic boundary conditions. We show how a doubled cell approach, previously used in non-CPM MD simulations of charged interfaces, can be used to enable fully periodic CPM MD simulations. Using either a doubled cell approach or a finite field approach previously reported by others, fully periodic CPM MD simulations produce comparable results to the traditional slab geometry simulations with a nearly double speedup in computational time. Indeed, these savings can offset the additional cost of the CPM algorithm, resulting in periodic CPM MD simulations that are computationally competitive with the non-periodic, fixed charge equivalent simulations for the ionic liquid supercapacitors studied here.

Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface

Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. We discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.

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.

An evaluation of the capacitive behavior of supercapacitors as a function of the radius of cations using simulations with a constant potential method

We report on the atomistic molecular dynamics, applying the constant potential method to determine the structural and electrostatic interactions at the electrode-electrolyte interface of electrochemical supercapacitors as a function of the cation radius (Cs⁺, Rb⁺, K⁺, Na⁺, Li⁺). We find that the electrical double layer is susceptible to the size, hydration layer volume, and cations' mobility and analyzed them. Besides, the transient potential shows an increase in magnitude and length as a function of the monocation size, i.e., Cs⁺ > Rb⁺ > K⁺ > Na⁺ > Li⁺. On the other hand, the charge distribution along the electrode surface is less uniform for large monocations. Nonetheless, the difference is not observed as a function of the radius of the cation for the integral capacitance. Our results are comparable to studies that employed the fixed charge method for treating such systems.

Interfacial water asymmetry at ideal electrochemical interfaces

Controlling electrochemical reactivity requires a detailed understanding of the charging behavior and thermodynamics of the electrochemical interface. Experiments can independently probe the overall charge response of the electrochemical double layer by capacitance measurements and the thermodynamics of the inner layer with potential of maximum entropy measurements. Relating these properties by computational modeling of the electrochemical interface has so far been challenging due to the low accuracy of classical molecular dynamics (MD) for capacitance and the limited time and length scales of ab initio MD. Here, we combine large ensembles of long-time-scale classical MD simulations with charge response from electronic density functional theory to predict the potential-dependent capacitance of a family of ideal aqueous electrochemical interfaces with different peak capacitances. We show that while the potential of maximum capacitance varies, this entire family exhibits an electrode charge of maximum capacitance (CMC) between -2.9 and -2.2 μC/cm², regardless of the details in the electronic response. Simulated heating of the same interfaces reveals that the entropy peaks at a charge of maximum entropy (CME) of -5.1 ± 0.6 μC/cm², in agreement with experimental findings for metallic electrodes. The CME and CMC both indicate asymmetric response of interfacial water that is stronger for negatively charged electrodes, while the difference between CME and CMC illustrates the richness in behavior of even the ideal electrochemical interface.

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.

Structural effects in nanotribology of nanoscale films of ionic liquids confined between metallic surfaces

Room Temperature Ionic Liquids (RTILs) attract significant interest in nanotribology. However, their microscopic lubrication mechanism is still under debate. Here, using non-equilibrium molecular dynamics simulations, we investigate the lubrication performance of ultra-thin (<2 nm) films of [C₂MIM]⁺ [NTf₂]^- confined between plane-parallel neutral surfaces of Au(111) or Au(100). We find that films consisting of tri-layers or bilayers, form ordered structures with a flat orientation of the imidazolium rings with respect to the gold surface plane. Tri-layers are unstable against loads >0.5 GPa, while bi-layers sustain pressures in the 1-2 GPa range. The compression of these films results in monolayers that can sustain loads of several GPa without significant loss in their lubrication performance. Surprisingly, in such ultra-thin films the imidazolium rings show higher orientational in-plane disorder, with and the rings adopting a tilted orientation with respect to the gold surface. The friction force and friction coefficient of the monolayers depends strongly on the structure of the gold plates, with the friction coefficient being four times higher for monolayers confined between Au(100) surfaces than for more compact Au(111) surfaces. We show that the general behaviour described here is independent of whether the metallic surfaces are modelled as polarizable or non-polarizable surfaces and speculate on the nature of this unexpected conclusion.

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.

Interfacial Polarization and Ionic Structure at the Ionic Liquid-Metal Interface Studied by Vibrational Spectroscopy and Molecular Dynamics Simulations

Ionic liquids (ILs) have both fundamental and practical value in interfacial science and electrochemistry. However, understanding their behavior near a surface is challenging because of strong Coulomb interactions and large and irregular ionic sizes, which affect both their structure and energetics. To understand this problem, we present a combined experimental and computational study using a vibrational probe molecule, 4-mercaptobenzonitrile, inserted at the junction between a metal and a variety of ILs. The vibrational frequency of the nitrile in the probe molecule reports on the local solvation environment and the electrostatic field at this junction. Within the ethylmethyl imidazolium (EMIM⁺) cation family of ILs, we varied the anions over a range of sizes and types. Complementing our surface spectroscopy, we also ran molecular dynamics simulations of these interfaces to better understand the ionic structures that produced the measured fields. The magnitude of the frequency shifts, and thereby fields, shows a general correlation with the size of anions, with larger anions corresponding to smaller fields. We find that the source of this correlation is partial intercalation of smaller anions into the probe monolayer, resulting in tighter packing of ionic layers near the surface. Larger anions reduce the overall lateral ion packing density near the surface, which reduces the net charge per unit area and explains the smaller observed fields. The insight from this work is important for developing a fundamental picture of concentrated electrolytes near interfaces and can help with designing ILs to create tailored electric fields near an electrode.

Effect of side chain modifications in imidazolium ionic liquids on the properties of the electrical double layer at a molybdenum disulfide electrode

Knowledge of how the molecular structures of ionic liquids (ILs) affect their properties at electrified interfaces is key to the rational design of ILs for electric applications. Polarizable molecular dynamics simulations were performed to investigate the structural, electrical, and dynamic properties of electric double layers (EDLs) formed by imidazolium dicyanamide ([ImX1][DCA]) at the interface with the molybdenum disulfide electrode. The effect of side chain of imidazolium on the properties of EDLs was analyzed by using 1-ethyl-3-methylimidazolium ([Im21]), 1-octyl-3-methylimidazolium ([Im81]), 1-benzyl-3-methylimidazolium ([ImB1]), and 1-(2-hydroxyethyl)-3-methylimidazolium ([ImO1]) as cations. Using [Im21] as reference, we find that the introduction of octyl or benzyl groups significantly alters the interfacial structures near the cathode because of the reorientation of cations. For [Im81], the positive charge on the cathode induces pronounced polar and non-polar domain separation. In contrast, the hydroxyl group has a minor effect on the interfacial structures. [ImB1] is shown to deliver slightly larger capacitance than other ILs even though it has larger molecular volume than [Im21]. This is attributed to the limiting factor for capacitance being the strong association between counter-ions, instead of the free space available to ions at the interface. For [Im81], the charging mechanism is mainly the exchange between anions and octyl tails, while for the other ILs, the mechanism is mainly the exchange of counter-ions. Analysis on the charging process shows that the charging speed does not correlate strongly with macroscopic bulk dynamics like viscosity. Instead, it is dominated by local displacement and reorientation of ions.

Three-Dimensional Molecular Mapping of Ionic Liquids at Electrified Interfaces

Electric double layers (EDLs), occurring ubiquitously at solid-liquid interfaces, are critical for electrochemical energy conversion and storage processes such as capacitive charging and redox reactions. However, to date the molecular-scale structure of EDLs remains elusive. Here we report an advanced technique, electrochemical three-dimensional atomic force microscopy (EC-3D-AFM), and use it to directly image the molecular-scale EDL structure of an ionic liquid under different electrode potentials. We observe not only multiple discrete ionic layers in the EDL on a graphite electrode but also a quasi-periodic molecular density distribution within each layer. Furthermore, we find pronounced 3D reconfiguration of the EDL at different voltages, especially in the first layer. Combining the experimental results with molecular dynamics simulations, we find potential-dependent molecular redistribution and reorientation in the innermost EDL layer, both of which are critical to EDL capacitive charging. We expect this mechanistic understanding to have profound impacts on the rational design of electrode-electrolyte interfaces for energy conversion and storage.

Evolution and Reversible Polarity of Multilayering at the Ionic Liquid/Water Interface

Highly correlated positioning of ions underlies Coulomb interactions between ions and electrified interfaces within dense ionic fluids such as biological cells and ionic liquids. Recent work has shown that highly correlated ionic systems behave differently than dilute electrolyte solutions, and interest is focused upon characterizing the electrical and structural properties of the dense electrical double layers (EDLs) formed at internal interfaces. It has been a challenge for experiments to characterize the progressive development of the EDL on the nanoscale as the interfacial electric potential is varied over a range of positive and negative values. Here we address this challenge by measuring X-ray reflectivity from the interface between an ionic liquid (IL) and a dilute aqueous electrolyte solution over a range of interfacial potentials from -450 to 350 mV. The growth of alternately charged cation-rich and anion-rich layers was observed along with a polarity reversal of the layers as the potential changed sign. These data show that the structural development of an ionic multilayer-like EDL with increasing potential is similar to that suggested by phenomenological theories and MD simulations, although our data also reveal that the excess charge beyond the first ionic layer decays more rapidly than predicted.

Mass-zero constrained molecular dynamics for electrode charges in simulations of electrochemical systems

Classical molecular dynamics simulations have recently become a standard tool for the study of electrochemical systems. State-of-the-art approaches represent the electrodes as perfect conductors, modeling their responses to the charge distribution of electrolytes via the so-called fluctuating charge model. These fluctuating charges are additional degrees of freedom that, in a Born-Oppenheimer spirit, adapt instantaneously to changes in the environment to keep each electrode at a constant potential. Here, we show that this model can be treated in the framework of constrained molecular dynamics, leading to a symplectic and time-reversible algorithm for the evolution of all the degrees of freedom of the system. The computational cost and the accuracy of the new method are similar to current alternative implementations of the model. The advantage lies in the accuracy and long term stability guaranteed by the formal properties of the algorithm and in the possibility to systematically introduce additional kinematic conditions of arbitrary number and form. We illustrate the performance of the constrained dynamics approach by enforcing the electroneutrality of the electrodes in a simple capacitor consisting of two graphite electrodes separated by a slab of liquid water.

Fast multipole method for three-dimensional systems with periodic boundary condition in two directions

We derived a new expression for the electrostatic interaction of three-dimensional charge-neutral systems with two-dimensional periodic boundary conditions (slab geometry) using a fast multipole method (FMM). Contributions from all the image cells are expressed as a sum of real and reciprocal space terms, and a self-interaction term. The reciprocal space contribution consists of two parts: zero and nonzero terms of the absolute value of the reciprocal lattice vector. To test the new expressions, electrostatic interactions were calculated for a randomly placed charge distribution in a cubic box and liquid water produced by molecular dynamics calculation. The accuracy could be controlled by the degree of expansion of the FMM. In the present expression, the computational complexity of the electrostatic interaction of N-particle systems is order N, which is superior to that of the conventional two-dimensional periodic Ewald method for a slab geometry and the particle mesh Ewald method with a large empty space at an interface of the unit cell. © 2020 Wiley Periodicals, Inc.

An electric double layer structure and differential capacitance at the electrode interface of tributylmethylammonium bis(trifluoromethanesulfonyl)amide studied using a molecular dynamics simulation

A molecular dynamics simulation at the electrode interface of a quaternary ammonium ionic liquid, tributylmethylammonium bis(trifluoromethanesulfonyl)amide ([N₁₄₄₄⁺][TFSA^-]), has been performed. Unlike the commonly used cations, such as 1-alkyl-3-methylimidazolium and 1,1-alkylmethylpyrrolidinium cations, N₁₄₄₄⁺ has multiple long-alkyl groups (three butyl groups). The behavior of ions at the electrode interface, especially these butyl groups, has been investigated. N₁₄₄₄⁺ at the first layer mainly has two types of orientations, lying and standing. The lying orientation is dominant at moderately negative potentials. However, the standing one becomes dominant at the more negative potentials. Due to this orientational change, the number of N₁₄₄₄⁺ increases at the first layer as the potential becomes negative even at the potentials where the anions are completely depleted there. The change in orientation results in the upward deviation of the differential capacitance from the theoretical prediction at the negative potentials. The results suggest that the orientational preference caused by the steric constraint between alkyl groups plays an important role in the behavior of the electric double layer of the ionic liquids.

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.

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields

Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.

Probing Electric Double-Layer Composition via in Situ Vibrational Spectroscopy and Molecular Simulations

At an electrode, ions and solvent accumulate to screen charge, leading to a nanometer-scale electric double layer (EDL). The EDL guides electrode passivation in batteries, while in (super)capacitors, it determines charge storage capacity. Despite its importance, quantification of the nanometer-scale and potential-dependent EDL remains a challenging problem. Here, we directly probe changes in the EDL composition with potential using in situ vibrational spectroscopy and molecular dynamics simulations for a Li-ion battery electrolyte (LiClO₄ in dimethyl carbonate). The accumulation rate of Li⁺ ions at the negative surface and ClO₄^- ions at the positive surface from vibrational spectroscopy compares well to that predicted by simulations using a polarizable APPLE&P force field. The ion solvation shell structure and ion-pairing within the EDL differs significantly from the bulk, especially at the negative electrode, suggesting that the common rationalization of interfacial electrochemical processes in terms of bulk ion solvation should be applied with caution.

Storing Energy in Biodegradable Electrochemical Supercapacitors

The development of green and biodegradable electrical components is one of the main fronts of research to overcome the growing ecological problem related to the issue of electronic waste. At the same time, such devices are highly desirable in biomedical applications such as integrated bioelectronics, for which biocompatibility is also required. Supercapacitors for storage of electrochemical energy, designed only with biodegradable organic matter would contemplate both aspects, that is, they would be ecologically harmless after their service lifetime and would be an important component for applications in biomedical engineering. By means of atomistic simulations of molecular dynamics, we propose a supercapacitor whose electrodes are formed exclusively by self-organizing peptides and whose electrolyte is a green amino acid ionic liquid. Our results indicate that this supercapacitor has a high potential for energy storage with superior performance than conventional supercapacitors. In particular its capacity to store energy was estimated to be almost 20 times greater than an analogue one of planar metallic electrodes.

Graphene/ionic liquid ultracapacitors: does ionic size correlate with energy storage performance?

Electrolytes formed by ions of similar volumes show better performance in supercapacitor applications.

Computational Insights into Materials and Interfaces for Capacitive Energy Storage

Supercapacitors such as electric double-layer capacitors (EDLCs) and pseudocapacitors are becoming increasingly important in the field of electrical energy storage. Theoretical study of energy storage in EDLCs focuses on solving for the electric double-layer structure in different electrode geometries and electrolyte components, which can be achieved by molecular simulations such as classical molecular dynamics (MD), classical density functional theory (classical DFT), and Monte-Carlo (MC) methods. In recent years, combining first-principles and classical simulations to investigate the carbon-based EDLCs has shed light on the importance of quantum capacitance in graphene-like 2D systems. More recently, the development of joint density functional theory (JDFT) enables self-consistent electronic-structure calculation for an electrode being solvated by an electrolyte. In contrast with the large amount of theoretical and computational effort on EDLCs, theoretical understanding of pseudocapacitance is very limited. In this review, we first introduce popular modeling methods and then focus on several important aspects of EDLCs including nanoconfinement, quantum capacitance, dielectric screening, and novel 2D electrode design; we also briefly touch upon pseudocapactive mechanism in RuO2. We summarize and conclude with an outlook for the future of materials simulation and design for capacitive energy storage.

Computational and Experimental Study of Li-Doped Ionic Liquids at Electrified Interfaces

We evaluate the influence of Li-salt doping on the dynamics, capacitance, and structure of three ionic liquid electrolytes, [pyr14][TFSI], [pyr13][FSI], and [EMIM][BF4], using molecular dynamics and polarizable force fields. In this respect, our focus is on the properties of the electric double layer (EDL) formed by the electrolytes at the electrode surface as a function of surface potential (Ψ). The rates of EDL formation are found to be on the order of hundreds of picoseconds and only slightly influenced by the addition of Li-salt. The EDLs of three electrolytes are shown to have different energy storage capacities, which we relate to the EDL formation free energy. The differential capacitance obtained from our computations exhibits asymmetry about the potential of zero charge and is consistent with the camel-like profiles noted from mean field theories and experiments on metallic electrodes. The introduction of Li-salt reduces the noted asymmetry in the differential capacitance profile. Complementary experimental capacitance measurements have been made on our three electrolytes in their neat forms and with Li-salt. The measurements, performed on glassy carbon electrodes, produce U-like profiles, and Li-salt doping is shown to strongly affect capacitance at high magnitudes of Ψ. Differences in the theoretical and experimental shapes and magnitudes of capacitance are rationalized in terms of the electrode surface and pseudocapacitive effects. In both neat and Li-doped liquids, the details of the computational capacitance profile are well described by Ψ-induced changes in the density and molecular orientation of ions in the molecular layer closest to the electrode. Our results suggest that the addition of Li+ induces disorder in the EDL, which originates from the strong binding of anions to Li+. An in-depth analysis of the distribution of Li+ in the EDL reveals that it does not readily enter the molecular layer at the electrode surface, preferring instead to be localized farther away from the surface in the second molecular layer. This behavior is validated through an analysis of the free energy of Li+ solvation as a function of distance from the electrode. Free energy wells are found to coincide with localized concentrations of Li+, the depths of which increase with Ψ and suggest a source of impedance for Li+ to reach the electrode. Finally, we make predictions of the specific energy at ideal graphite utilizing the computed capacitance and previously derived electrochemical windows of the liquids.