Evaluation of the constant potential method in simulating electric double-layer capacitors

The influence of water polarization on slip friction at charged interfaces

The present study employs equilibrium molecular dynamics simulations to explore the potential mechanism for controlling friction by applying electrostatic fields in nanoconfined aqueous electrolytes. The slip friction coefficient demonstrates a gradual increase corresponding to the surface charge density for pure water and aqueous electrolytes, exhibiting a similar trend across both nanochannel walls. An expression is formulated to rationalize the observed slip friction behavior, describing the effect of the electric field on the slip friction coefficient. According to this formulation, the slip friction coefficient increases proportionally to the square of the uniform electric field emanating from the charged electrode. This increase in slip friction results from the energy change due to the orientation polarization of interfacial water dipoles. The minimal variations in the empirically determined proportionality constant for pure water and aqueous electrolytes indicate that water polarization primarily governs slip friction at charged interfaces. These findings offer insights into the electrical effects on nanoscale lubrication of aqueous electrolytes, highlighting the significant role of water polarization in determining slip.

Li-ion solvation structure at electrified solid-liquid interface: Understanding solvation structure dynamics and its role in electrochemical energy storage through binary ethylene carbonate and dimethyl carbonate solvent

Binary solvent electrolytes can provide interpretations for designing advanced electrolytes of next generation batteries. This study investigates the adsorption mechanisms of solvated lithium ions in binary solvents near charged electrodes. Molecular dynamic simulations are performed for lithium hexafluorophosphate (LiPF6) in ethylene carbonate and dimethyl carbonate (EC:DMC) solvent sandwiched between two electrodes. Results show that lithium ions form a tetrahedral solvation structure with two EC and two DMC molecules. The solvated lithium ion shows anti-electrostatic interaction with electrodes. This can be attributed to the electrostatic attraction of the polar end of the DMC molecule, which keeps the cation anchored to the positive electrode. Meanwhile, the solvation structure adopts a fix orientation at the negative electrode, which leads to unchanged electrostatic interaction at high charge density. Finally, EC molecules are swapped by DMC molecules near the negative electrode at high charge density. This leads to a decrease in local relative permittivity and, therefore, a decrease in differential capacitance. The differential capacitance of the positive electrode continuously decreases with increasing charge density. This is caused by the partial anchoring of solvent molecules holding the cations, which cancels the adsorption of anions near the positive electrode. This study provides insights into designing better electrolytes for efficient battery performance.

Symmetric effect on electrical double-layer characteristics and molecular assembly interplay in imidazolium-based Ionic liquid electrolytes in supercapacitor models

Studies on the ion-layer formation of imidazolium-based ionic liquids have extensively explored how to improve in-depth knowledge of electrical double-layer (EDL) properties. In this computational study, 1-alkyl-methylimidazolium bis(trifluoromethylsulfonyl)imide ([Cnmim][NTf2]), namely, [C1mim][NTf2] and [C2mim][NTf2], inside a simulated supercapacitor were investigated to expose an symmetric alkyl chain effect. Molecular dynamic simulations of a supercapacitor model with graphite electrodes were conducted. Changes in charging dynamics and EDL structures at different voltages were studied. Although [C1mim][NTf2] equilibrated much quicker than [C2mim][NTf2], surface charge development on the symmetrical imidazolium ionic liquid was slower than that on the asymmetrical counterpart. Physical EDL structural analysis showed that [C1mim][NTf2] could not rearrange in a rigid co-ion layer, whereby the [C1mim]+ cation stayed adsorbed on the positive electrode throughout all the tested voltages. The strongly attached [C1mim]+ on the electrode surface contributed to low responsiveness in symmetrical [C1mim][NTf2], which was supported by lower overall differential capacitance (CD) magnitude and less sharp CD wings at high voltage when compared to [C2mim][NTf2].

Exploring the Impact of In Situ-Formed Solid-Electrolyte Interphase on the Cycling Performance of Aluminum Metal Anodes

Unwanted processes in metal anode batteries, e.g., non-uniform metal electrodeposition, electrolyte decomposition, and/or short-circuiting, are not fully captured by the electrolyte bulk solvation structure but rather defined by the electrode-electrolyte interface and its changes induced by cycling conditions. Specifically, for aluminum-ion batteries (AIBs), the role of the solid-electrolyte interphase (SEI) on the Al0 electrodeposition mechanism and associated changes during resting or cycling remain unclear. Here, we investigated the current-dependent changes at the electrified aluminum anode/ionic liquid electrolyte interface to reveal the conditions of the SEI formation leading to irreversible cycling in the AIBs. We identified that the mechanism of anode failure depends on the nature of the counter electrode, where the areal capacity and cycling current for Al0 electrodeposition dictates the number of successful cycles. Notwithstanding the differences behind unstable aluminum anode cycling in symmetrical cells and AIBs, the uniform removal of electrochemically inactive SEI components, e.g., oxide-rich or solvent-derived organic-rich interphases, leads to more efficient cycling behavior. These understandings raise the importance of using specific conditioning protocols for efficient cycling of the aluminum anode in conjugation with different cathode materials.

Selective adsorption of divalent and trivalent cations in porous electrodes

The capacitive deionization technology uses the electrochemical adsorption of ions in porous electrodes to desalinate seawater or brackish water. Recently, capacitive deionization has gained significant attention as a technology for selective adsorption of ionic species from multicomponent aqueous electrolytes. To investigate the mechanism of selective adsorption at the molecular level, we performed molecular dynamics simulations of aqueous electrolytes and porous electrodes with different divalent or trivalent ions, electrode pore sizes, and applied voltages. We calculated the free energy barriers preventing ions from entering the pores of the electrode and the structure of the water molecules near the ions and the electrode surface under various conditions. Our results suggest that, when the pore and ion sizes are comparable, the steric and electrostatic interactions between the hydrated ions and electrode pores are comparable in magnitude. Moreover, the relative importance of the two interactions can be reversed by slight changes in the external conditions, such as the ion size, valence of the ions, electrode pore size, and applied voltage. Thus, by finely tuning the electrode pore size and the applied voltage, it may be possible to selectively adsorb a particular ionic species from a multicomponent electrolyte through capacitive deionization using a porous electrode.

Modeling of Nanomaterials for Supercapacitors: Beyond Carbon Electrodes

Capacitive storage devices allow for fast charge and discharge cycles, making them the perfect complements to batteries for high power applications. Many materials display interesting capacitive properties when they are put in contact with ionic solutions despite their very different structures and (surface) reactivity. Among them, nanocarbons are the most important for practical applications, but many nanomaterials have recently emerged, such as conductive metal-organic frameworks, 2D materials, and a wide variety of metal oxides. These heterogeneous and complex electrode materials are difficult to model with conventional approaches. However, the development of computational methods, the incorporation of machine learning techniques, and the increasing power in high performance computing now allow us to tackle these types of systems. In this Review, we summarize the current efforts in this direction. We show that depending on the nature of the materials and of the charging mechanisms, different methods, or combinations of them, can provide desirable atomic-scale insight on the interactions at play. We mainly focus on two important aspects: (i) the study of ion adsorption in complex nanoporous materials, which require the extension of constant potential molecular dynamics to multicomponent systems, and (ii) the characterization of Faradaic processes in pseudocapacitors, that involves the use of electronic structure-based methods. We also discuss how recently developed simulation methods will allow bridges to be made between double-layer capacitors and pseudocapacitors for future high power electricity storage devices.

Improving Aqueous Zinc Ion Batteries with Alkali Metal Ions

Aqueous zinc (Zn) ion batteries have received broad attention recently. However, their practical application is limited by severe Zn dendrite growth and the hydrogen evolution reaction. In this study, three alkali metal ions (Li+, Na+, and K+) are added in ZnSO4 electrolytes, which are subjected to electrochemical measurements and molecular dynamics simulations. The studies show that since K+ has the highest mobility and self-diffusion coefficient among the four ions (Li+, Na+, K+, and Zn2+), it enables K+ to preferentially approach a zinc dendrite at an earlier time, driven by a negative electric field during a cathodic process. The electric double layer, with K+ around the negatively charged Zn dendrite, inhibits dendrite growth and mitigates the hydrogen evolution reaction on the Zn anode. Under this kinetic effect, the Zn-Zn symmetric cell with K+ exhibits a long cycling stability of 1000 h at 1 mA·cm-2 of 1 mAh·cm-2 and 190 h at 30 mA·cm-2 of 2 mAh·cm-2. Such a kinetic effect is also observed with additives Na+ and Li+, though less profound than that of K+.

Molecular insights into the nanoconfinement effect on the structure and dynamics of ionic liquids in carbon nanotubes

The properties and applications of ionic liquids (ILs) have been widely investigated when they are confined within nanochannels such as carbon nanotubes (CNTs). The confined ILs exhibit very different properties from their bulk state due to a nanoconfinement effect, which plays an important role in the performances of devices with ILs. In this work, we studied the effect of the charge carried by CNTs on confined ILs inside CNTs using molecular dynamics simulations. In charged CNTs, cations and anions are distributed separately along the radial directions, and the transition of orientations of the cations between parallel and vertical to CNTs occurs by changing the charge state of CNTs. The number of hydrogen bonds (HBs) formed by the confined ILs can be reduced by switching the surface charge of CNTs from positive to negative due to the contact modes between cations and anions as well as the distributions of cations in CNTs. The diffusivities along and vertical to the axial direction of CNTs were found to be non-monotonic owing to the "trade-off" effect from both ion pair interlocking and anchoring ILs on the CNT walls. Additionally, the region-dependent dynamics of ILs were also related to the intermolecular interactions in different regions of CNTs. Furthermore, the vibrational modes of ILs were obviously influenced in highly charged CNTs as determined by calculating the density of vibrational states, which demonstrated the transitions in the structure and interactions. The density distributions changed from single layer to double layers when increasing the pore size of neutral CNTs while the hydrogen bonds exhibited a non-monotonic tendency versus the pore sizes. Our results might help to understand the structure and dynamics of confined ILs as well as aid optimizing the performance of devices with ILs.

Understanding electrochemical interfaces through comparing experimental and computational charge density-potential curves

Electrode-electrolyte interfaces play a decisive role in electrochemical charge accumulation and transfer processes. Theoretical modelling of these interfaces is critical to decipher the microscopic details of such phenomena. Different force field-based molecular dynamics protocols are compared here in a view to connect calculated and experimental charge density-potential relationships. Platinum-aqueous electrolyte interfaces are taken as a model. The potential of using experimental charge density-potential curves to transform cell voltage into electrode potential in force-field molecular dynamics simulations, and the need for that purpose of developing simulation protocols that can accurately calculate the double-layer capacitance, are discussed.

Reactive molecular dynamics simulations of lithium-ion battery electrolyte degradation

The development of reliable computational methods for novel battery materials has become essential due to the recently intensified research efforts on more sustainable energy storage materials. Here, we use a recently developed framework allowing to consistently incorporate quantum-mechanical activation barriers to classical molecular dynamics simulations to study the reductive solvent decomposition and formation of the solid electrolyte interphase for a graphite/carbonate electrolyte interface. We focus on deriving condensed-phase effective rates based on the elementary gas-phase reduction and decomposition energy barriers. After a short initial transient limited by the elementary barriers, we observe that the effective rate shows a transition to a kinetically slow regime influenced by the changing coordination environment and the ionic fluxes between the bulk electrolyte and the interface. We also discuss the impact of the decomposition on the ionic mobility. Thus, our work shows how elementary first-principles properties can be mechanistically leveraged to provide fundamental insights into electrochemical stability of battery electrolytes.

Structure and self-diffusivity of mixed-cation electrolytes between neutral and charged graphene sheets

Graphene-based applications, such as supercapacitors or capacitive deionization, take place in an aqueous environment, and they benefit from molecular-level insights into the behavior of aqueous electrolyte solutions in single-digit graphene nanopores with a size comparable to a few molecular diameters. Under single-digit graphene nanoconfinement (smallest dimension <2 nm), water and ions behave drastically different than in the bulk. Most aqueous electrolytes in the graphene-based applications as well as in nature contain a mix of electrolytes. We study several prototypical aqueous mixed alkali-chloride electrolytes containing an equimolar fraction of Li/Na, Li/K, or Na/K cations confined between neutral and positively or negatively charged parallel graphene sheets. The strong hydration shell of small Li+ vs a larger Na+ or large K+ with weaker or weak hydration shells affects the interplay between the ions's propensity to hydrate or dehydrate under the graphene nanoconfinement and the strength of the ion-graphene interactions mediated by confinement-induced layered water. We perform molecular dynamics simulations of the confined mixed-cation electrolytes using the effectively polarizable force field for electrolyte-graphene systems and focused on a relation between the electrochemical adsorption and structural properties of the water molecules and ions and their diffusion behavior. The simulations show that the one-layer nanoslits have the biggest impact on the ions' adsorption and the water and ions' diffusion. The positively charged one-layer nanoslits only allow for Cl- adsorption and strengthen the intermolecular bonding, which along with the ultrathin confinement substantially reduces the water and Cl- diffusion. In contrast, the negatively charged one-layer nanoslits only allow for adsorption of weakly hydrated Na+ or K+ and substantially break up the non-covalent bond network, which leads to the enhancement of the water and Na+ or K+ diffusion up to or even above the bulk diffusion. In wider nanoslits, cations adsorb closer to the graphene surfaces than Cl-'s with preferential adsorption of a weakly hydrated cation over a strongly hydrated cation. The positive graphene charge has an intuitive effect on the adsorption of weakly hydrated Na+'s or K+'s and Cl-'s and a counterintuitive effect on the adsorption of strongly hydrated Li+'s. On the other hand, the negative surface charge has an intuitive effect on the adsorption of both types of cations and only mild intuitive or counterintuitive effects on the Cl- adsorption. The diffusion of water molecules and ions confined in the wider nanoslits is reduced with respect to the bulk diffusion, more for the positive graphene charge, which strengthened the intermolecular bonding, and less for the negative surface charge, which weakened the non-covalent bond network.

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.

Suppressing Zn pulverization with three-dimensional inert-cation diversion dam for long-life Zn metal batteries

Tremendous attention has been paid to the water-associated side reactions and zinc (Zn) dendrite growth on the electrode-electrolyte interface. However, the Zn pulverization that can cause continuous depletion of active Zn metal and exacerbate hydrogen evolution is severely neglected. Here, we disclose that the excessive Zn feeding that causes incomplete crystallization is responsible for Zn pulverization formation through analyzing the thermodynamic and kinetics process of Zn deposition. On the basis, we introduce 1-ethyl-3-methylimidazolium cations (EMIm+) into the electrolyte to form a Galton-board-like three-dimensional inert-cation (3DIC) region. Modeling test shows that the 3DIC EMIm+ can induce the Zn2+ flux to follow in a Gauss distribution, thus acting as elastic sites to buffer the perpendicular diffusion of Zn2+ and direct the lateral diffusion, thus effectively avoiding the local Zn2+ accumulation and irreversible crystal formation. Consequently, anti-pulverized Zn metal deposition behavior is achieved with an average Coulombic efficiency of 99.6% at 5 mA cm-2 over 2,000 cycles and superb stability in symmetric cell over 1,200 h at -30 °C. Furthermore, the Zn||KVOH pouch cell can stably cycle over 1,200 cycles at 2 A g-1 and maintain a capacity of up to 12 mAh.

Water dynamics and sum-frequency generation spectra at electrode/aqueous electrolyte interfaces

The dynamics of water at interfaces between an electrode and an electrolyte is essential for the transport of redox species and for the kinetics of charge transfer reactions next to the electrode. However, while the effects of electrode potential and ion concentration on the electric double layer structure have been extensively studied, a comparable understanding of dynamical aspects is missing. Interfacial water dynamics presents challenges since it is expected to result from the complex combination of water-water, water-electrode and water-ion interactions. Here we perform molecular dynamics simulations of aqueous NaCl solutions at the interface with graphene electrodes, and examine the impact of both ion concentration and electrode potential on interfacial water reorientational dynamics. We show that for all salt concentrations water dynamics exhibits strongly asymmetric behavior: it slows down at increasingly positively charged electrodes but it accelerates at increasingly negatively charged electrodes. At negative potentials water dynamics is determined mostly by the electrode potential value, but in contrast at positive potentials it is governed both by ion-water and electrode-water interactions. We show how these strikingly different behaviors are determined by the interfacial hydrogen-bond network structure and by the ions' surface affinity. Finally, we indicate how the structural rearrangements impacting water dynamics can be probed via vibrational sum-frequency generation spectroscopy.

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.

Effect of Interlayer Spaces and Interfacial Structures on High-Performance MXene/Ionic Liquid Supercapacitors: A Molecular Dynamics Simulation

The combination of high-capacitance MXenes and wide-electrochemical-window ionic liquids (ILs) has exhibited bright prospects in supercapacitors. Several strategies, such as surficial functionalization and interlayer spacing tuning, have been used to enhance the electrochemical performance of supercapacitors. However, the lack of theoretical guidance on these strategies, including the effects of the microenvironment in the interlayer of confined ILs, hindered the further exploration of such devices. Herein, we performed molecular dynamics simulations to comprehensively investigate the effects of the interlayer space and surface terminations of MXene electrodes on capacity. The results show that the electrical double layer (EDL) structure was found to form on the interface between the MXene electrode and ILs electrolyte by analyzing the ion number density and charge density in the nanometer confined spaces. Under the same potential, the -OH terminations significantly impact the ion orientation in the EDL, particularly near the electrode surface, where cations tend to align vertically, allowing the retention of more cations at the electrode surfaces. Interestingly, such an orientation distribution was decisively from the hydrogen bonds expressed by O-H···O between the -OH termination of MXene and -OH groups of ILs. The differential capacitances of the supercapacitors were calculated by the surficial electron density, and it showed that the capacitance is a nearly one-quarter increase in the 14 Å interlayer spacing compared with that of 10 Å under an applied potential of 2 V. At the same time, the Ti3C2(OH)2 electrode had a higher differential capacitance than the Ti3C2O2 electrode, which possibly originates from the stronger hydrogen bonds to contribute to the vertical aggregation of the cations. Our results highlighted the roles of the interlayer spacing distance and surface terminations of the MXene on the performance of the type of supercapacitor.

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.

Development of Heteroatomic Constant Potential Method with Application to MXene-Based Supercapacitors

We describe a method for modeling constant-potential charges in heteroatomic electrodes, keeping pace with the increasing complexity of electrode composition and nanostructure in electrochemical research. The proposed "heteroatomic constant potential method" (HCPM) uses minimal added parameters to handle differing electronegativities and chemical hardnesses of different elements, which we fit to density functional theory (DFT) partial charge predictions in this paper by using derivative-free optimization. To demonstrate the model, we performed molecular dynamics simulations using both HCPM and conventional constant potential method (CPM) for MXene electrodes with Li-TFSI/AN (lithium bis(trifluoromethane sulfonyl)imide/acetonitrile)-based solvent-in-salt electrolytes. Although the two methods show similar accumulated charge storage on the electrodes, the results indicated that HCPM provides a more reliable depiction of electrode atom charge distribution and charge response compared with CPM, accompanied by increased cationic attraction to the MXene surface. These results highlight the influence of elemental composition on electrode performance, and the flexibility of our HCPM opens up new avenues for studying the performance of diverse heteroatomic electrodes including other types of MXenes, two-dimensional materials, metal-organic frameworks (MOFs), and doped carbonaceous electrodes.

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.

Strain induced electrochemical behaviors of ionic liquid electrolytes in an electrochemical double layer capacitor: Insights from molecular dynamics simulations

Electrochemical Double Layer Capacitors (EDLCs) with ionic liquid electrolytes outperform conventional ones using aqueous and organic electrolytes in energy density and safety. However, understanding the electrochemical behaviors of ionic liquid electrolytes under compressive/tensile strain is essential for the design of flexible EDLCs as well as normal EDLCs, which are subject to external forces during assembly. Despite many experimental studies, the compression/stretching effects on the performance of ionic liquid EDLCs remain inconclusive and controversial. In addition, there is hardly any evidence of prior theoretical work done in this area, which makes the literature on this topic scarce. Herein, for the first time, we developed an atomistic model to study the processes underlying the electrochemical behaviors of ionic liquids in an EDLC under strain. Constant potential non-equilibrium molecular dynamics simulations are conducted for EMIM BF4 placed between two graphene walls as electrodes. Compared to zero strain, low compression of the EDLC resulted in compromised performance as the electrode charge density dropped by 29%, and the performance reduction deteriorated significantly with a further increase in compression. In contrast, stretching is found to enhance the performance by increasing the charge storage in the electrodes by 7%. The performance changes with compression and stretching are due to changes in the double-layer structure. In addition, an increase in the value of the applied potential during the application of strain leads to capacity retention with compression revealed by the newly performed simulations.

The concerted proton-electron transfer mechanism of proton migration in the electrochemical interface

The proton migration in the electrochemical interface is a fundamental electrochemical processes in proton involved reactions. We find fractional electron transfer, which is inversely proportional to the distance between the proton and electrode, during the proton migration under constant potential. The electrical energy carried by the transferred charge facilitates the proton to overcome the chemical barrier in the migration pathway, which is accounting for more than half electrical energy in the proton involved reactions. Consequently, less charge transfer and energy exchange take place in the reduction process. Therefore, the proton migration in the electrochemical interface is an essential component of the electrochemical reaction in terms of electron transfer and energy conversation, and are worthy of more attention in the rational design and optimization of electrochemical systems.

Enabling High-Capacitance Supercapacitors by Polyelectrolyte Brushes

Polyelectrolyte brushes (PEBs) hold excellent potential for designing high-capacitance electrical double-layer capacitors (EDLCs), a crucial component of supercapacitors. Both experiments and computational simulations have shown their energy-storage advantage. However, the effect of PEBs on the energy storage of EDLCs is not yet fully understood. Herein, we systematically study the energy-storage effects of polyanionic (PA) and polycationic (PC) brushes using polymer density functional theory (DFT). First, the application of polymer DFT in polyelectrolyte-grafted EDLCs is successfully validated using molecular dynamics simulations. With the help of polymer DFT, an interfacial adhesion microstructure of the PA/PC brushes is observed. Most importantly, the results show that polyelectrolyte-grafted EDLCs achieve a significant increase in capacitance at low salt concentrations and surface voltages, offering an excellent energy-storage advantage over traditional EDLCs. However, this advantage is considerably diminished at high salt concentrations or surface voltages, showing unusual salt- and voltage-dependent behaviors of energy-storage capacity. Nonetheless, the PC-grafted EDLCs maintain their outstanding energy-storage performance, even at relatively high salt concentrations and surface voltages. These findings deepen our comprehension of PEBs at the molecular level and provide insights for the molecular design of high-capacitance supercapacitors.

Molecular dynamics simulations of electrified interfaces including the metal polarisation

Understanding electrified interfaces requires an accurate description of the electric double layer which also takes into account the metal polarisation. Here we present a simple approach to the molecular dynamics simulation of electrified interfaces which combines fixed charges and a core-shell model for the description of the polarisable electron density on the metal electrode. The approach has been applied to the Au(111) surface in contact with a NaCl aqueous electrolyte solution in order to calculate the differential capacitance and to gain a detailed picture of the charging mechanism. Metal polarisation enhances the interfacial capacitance with a difference between the cathode and anode. In particular, we find that the influence of the metal polarisation on the electric double layer depends on the ion's solvation shell structure and, for the investigated interface, is more important at the cathode, where it modifies the sodium ion distribution.

Structure and self-diffusivity of alkali-halide electrolytes in neutral and charged graphene nanochannels

Understanding the microscopic behaviour of aqueous electrolyte solutions in graphene-based ultrathin nanochannels is important in nanofluidic applications such as water purification, fuel cells, and molecular sensing. Under extreme confinement (<2 nm), the properties of water and ions differ drastically from those in the bulk phase. We studied the structural and diffusion behaviour of prototypical aqueous solutions of electrolytes (LiCl, NaCl, and KCl) confined in both neutral and positively-, and negatively-charged graphene nanochannels. We performed molecular dynamics simulations of the solutions in the nanochannels with either one, two- or three-layer water structures using the effectively polarisable force field for graphene. We analysed the structure and intermolecular bond network of the confined solutions along with their relation to the self-diffusivity of water and ions. The simulations show that Na and K cations can more easily rearrange their solvation shells under the graphene nanoconfinement and adsorb on the graphene surfaces or dissolve in the confinement-induced layered water than the Li cation. The negative surface charge together with the presence of ions orient water molecules with hydrogens towards the graphene surfaces, which in turn weakens the intermolecular bond network. The one-layer nanochannels have the biggest effect on the water structure and intermolecular bonding as well as on the adsorption of ions with only co-ions entering these nanochannels. The self-diffusivity of confined water is strongly reduced with respect to the bulk water and decreases with diminishing nanochannel heights except for the negatively-charged one-layer nanochannel. The self-diffusivity of ions also decreases with the reducing the nanochannel heights except for the self-diffusivity of cations in the negatively-charged one-layer nanochannel, evidencing cooperative diffusion of confined water and ions. Due to the significant break-up of the intermolecular bond network in the negatively-charged one-layer nanochannel, self-diffusion coefficients of water and cations exceed those for the two- and three-layer nanochannels and become comparable to the bulk values.

Charge equilibration model with shielded long-range Coulomb for reactive molecular dynamics simulations

Atomic description of electrochemical systems requires reactive interaction potential to explicitly describe the chemistry between atoms and molecules and the evolving charge distribution and polarization effects. Calculating Coulomb electrostatic interactions and polarization effects requires a better estimate of the partial charge distribution in molecular systems. However, models such as reactive force fields and charge equilibration (QEq) include Coulomb interactions up to a short-distance cutoff for better computational speeds. Ignoring long-distance electrostatic interaction affects the ability to describe electrochemistry in large systems. We studied the long-range Coulomb effects among charged particles and extended the QEq method to include long-range effects. By this extension, we anticipate a proper account of Coulomb interactions in reactive molecular dynamics simulations. We validate the approach by computing charges on a series of metal-organic frameworks and some simple systems. Results are compared to regular QEq and quantum mechanics calculations. The study shows slightly overestimated charge values in the regular QEq approach. Moreover, our method was combined with Ewald summation to compute forces and evaluate the long-range effects of simple capacitor configurations. There were noticeable differences between the calculated charges with/without long-range Coulomb interactions. The difference, which may have originated from the long-range influence on the capacitor ions, makes the Ewald method a better descriptor of Coulomb electrostatics for charged electrodes. The approach explored in this study enabled the atomic description of electrochemical systems with realistic electrolyte thickness while accounting for the electrostatic effects of charged electrodes throughout the dielectric layer in devices like batteries and emerging solid-state memory.

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.

Rational Design of Solid Polymer Electrolyte Based on Ionic Liquid Monomer for Supercapacitor Applications via Molecular Dynamics Study

The safety concern arising from flammable liquid electrolytes used in batteries and supercapacitors drives technological advances in solid polymer electrolytes (SPEs) in which flammable organic solvents are absent. However, there is always a trade-off between the ionic conductivity and mechanical properties of SPEs due to the lack of interaction between the ionic liquid and polymer resin. The inadequate understanding of SPEs also limits their future exploitation and applications. Herein, we provide a complete approach to develop a new SPE, consisting of a cation (monomer), anion and hardener from ions-monomers using molecular dynamics (MD) simulations. The results show that the strong solid-liquid interactions between the SPE and graphene electrode lead to a very small gap of ∼5.5 Å between the components of SPE and electrode, resulting in a structured solid-to-liquid interface, which can potentially improve energy storage performance. The results also indicated the critical role of the mobility of free-standing anions in the SPE network to achieve high ionic conductivity for applications requiring fast charge/discharge. In addition, the formations of hardener-depleted regions and cation-anion-poor/rich regions near the uncharged/charged electrode surfaces were observed at the molecular level, providing insights for rationally designing the SPEs to overcome the boundaries for further breakthroughs in energy storage technology.

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

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

Constant-Potential Molecular Dynamics Simulations of Molten-Salt Double Layers for FLiBe and FLiNaK

We report the results of constant-potential molecular dynamics simulations of the double- layer interface between molten FLiBe and FLiNaK fluoride mixtures and idealized solid electrodes. Employing methods similar to those used in studies of chloride double layers, we compute the structure and differential capacitance of molten fluoride electric double layers as a function of applied voltage. The role of molten salt structure is probed through comparisons between FLiBe and FLiNaK, which serve as models for strong and weak associate- forming salts, respectively. In FLiBe, screening involves changes in Be-F-Be angles and alignment of the oligomers parallel to the electrode, while in FLiNaK the electric field is screened mainly by rearrangement of individual ions, predominantly the polarizable potassium cation.

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.

Applying Classical, Ab Initio, and Machine-Learning Molecular Dynamics Simulations to the Liquid Electrolyte for Rechargeable Batteries

Rechargeable batteries have become indispensable implements in our daily life and are considered a promising technology to construct sustainable energy systems in the future. The liquid electrolyte is one of the most important parts of a battery and is extremely critical in stabilizing the electrode-electrolyte interfaces and constructing safe and long-life-span batteries. Tremendous efforts have been devoted to developing new electrolyte solvents, salts, additives, and recipes, where molecular dynamics (MD) simulations play an increasingly important role in exploring electrolyte structures, physicochemical properties such as ionic conductivity, and interfacial reaction mechanisms. This review affords an overview of applying MD simulations in the study of liquid electrolytes for rechargeable batteries. First, the fundamentals and recent theoretical progress in three-class MD simulations are summarized, including classical, ab initio, and machine-learning MD simulations (section 2). Next, the application of MD simulations to the exploration of liquid electrolytes, including probing bulk and interfacial structures (section 3), deriving macroscopic properties such as ionic conductivity and dielectric constant of electrolytes (section 4), and revealing the electrode-electrolyte interfacial reaction mechanisms (section 5), are sequentially presented. Finally, a general conclusion and an insightful perspective on current challenges and future directions in applying MD simulations to liquid electrolytes are provided. Machine-learning technologies are highlighted to figure out these challenging issues facing MD simulations and electrolyte research and promote the rational design of advanced electrolytes for next-generation rechargeable batteries.

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.

Engineering a passivating electric double layer for high performance lithium metal batteries

In electrochemical devices, such as batteries, traditional electric double layer (EDL) theory holds that cations in the cathode/electrolyte interface will be repelled during charging, leaving a large amount of free solvents. This promotes the continuous anodic decomposition of the electrolyte, leading to a limited operation voltage and cycle life of the devices. In this work, we design a new EDL structure with adaptive and passivating properties. It is enabled by adding functional anionic additives in the electrolyte, which can selectively bind with cations and free solvents, forming unique cation-rich and branch-chain like supramolecular polymer structures with high electrochemical stability in the EDL inner layer. Due to this design, the anodic decomposition of ether-based electrolytes is significantly suppressed in the high voltage cathodes and the battery shows outstanding performances such as super-fast charging/discharging and ultra-low temperature applications, which is extremely hard in conventional electrolyte design principle. This unconventional EDL structure breaks the inherent perception of the classical EDL rearrangement mechanism and greatly improve electrochemical performances of the device.

EDL structure of ionic liquid-MXene-based supercapacitor and hydrogen bond role on the interface: a molecular dynamics simulation investigation

As a new class of electrodes, MXenes have shown excellent performance in supercapacitors. At the same time, ionic liquid (IL) electrolytes with wider electrochemical windows are expected to substantially increase the supercapacitor capacitance. The combination of MXenes and ILs is promising for energy storage devices with a high energy density and power density. The studies have indicated that the surface terminations of MXenes and the functional groups of ILs, can both strongly influence the supercapacitor's performance. However, studies at the molecular level are still lacking. In this work, we performed molecular dynamics simulations to investigate the interfacial structures and their influence on the energy storage mechanism. The results show that the two ILs exhibit very different charging rates, though the charge densities are similar after charging equilibrium. The interfacial analysis reveals different electrical double-layer (EDL) structures, in which most cations stay perpendicular to the Ti₃C₂(OH)₂ electrode when some cations shift to a vertical arrangement near the Ti₃C₂O₂ electrode. Such structures have led to the higher capacitance of the Ti₃C₂(OH)₂ electrode, even more than 2 times that of the Ti₃C₂O₂ electrode as the potential difference ranges from 0 to 2 V. It was also found that hydrogen bonds between the -OH groups of HEMIm⁺ cations and terminations of the MXene play an important role in improving the capacitances by aggregating more HEMIm⁺ cations on the surface of the Ti₃C₂(OH)₂ electrode. Our work provides clear mechanistic evidence that both terminations of the MXene electrodes and functional groups of the IL electrolytes affect the interfacial structures and the EDL formation, further leading to the different supercapacitor performance, which will be helpful in designing highly efficient energy-storage devices.

Electromagnetic bioeffects: a multiscale molecular simulation perspective

Electromagnetic bioeffects remain an enigma from both the experimental and theoretical perspectives despite the ubiquitous presence of related technologies in contemporary life. Multiscale computational modelling can provide valuable insights into biochemical systems and predict how they will be perturbed by external stimuli. At a microscopic level, it can be used to determine what (sub)molecular scale reactions various stimuli might induce; at a macroscopic level, it can be used to examine how these changes affect dynamic behaviour of essential molecules within the crowded biomolecular milieu in living tissues. In this review, we summarise and evaluate recent computational studies that examined the impact of externally applied electric and electromagnetic fields on biologically relevant molecular systems. First, we briefly outline the various methodological approaches that have been employed to study static and oscillating field effects across different time and length scales. The practical value of such modelling is then illustrated through representative case-studies that showcase the diverse effects of electric and electromagnetic field on the main physiological solvent - water, and the essential biomolecules - DNA, proteins, lipids, as well as some novel biomedically relevant nanomaterials. The implications and relevance of the theoretical multiscale modelling to practical applications in therapeutic medicine are also discussed. Finally, we summarise ongoing challenges and potential opportunities for theoretical modelling to advance the current understanding of electromagnetic bioeffects for their modulation and/or beneficial exploitation in biomedicine and industry.

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.

Understanding of Bulk and Interfacial Structures Ternary and Binary Deep Eutectic Solvents with a Constant Potential Method: A Molecular Dynamics Study

In the last decade, deep eutectic solvents (DESs) emerge as promising electrolytes in supercapacitors and rechargeable batteries due to their unique properties, wide electrochemical window, low viscosity, and high ionic...

Salt-in-water and water-in-salt electrolytes: the effects of the asymmetry in cation and anion valence on their properties

We investigated the structural, dynamic, energetic, and electrostatic properties of electrolytes based on the ion pairs LiCl and Li₂SO₄. Atomistic molecular dynamics simulations were used to simulate these aqueous electrolytic solutions at two different concentrations 2 M (normal) and 21 M (superconcentrated, WiSE). The effects of the valence asymmetry of the Li₂SO₄ electrolyte were also discussed for both salt concentrations. Our results differ in the physical aspect of pure electrolytes, showing the drastic effect of high concentration, in particular on the viscosity, which is dramatically increased in WiSE. This is a consequence of their reduced ionic mobility and has a direct effect on ionic conductivity. Also, our results for graphene-based supercapacitors, as indicated by some experimental work, do not indicate any better performance of WiSEs over normal electrolytes. In fact, the differences in the total capacitance, due to the concentration of ions, presented by both electrolytes are negligible. The valence asymmetry can be clearly observed in some properties but for most of them its effects could not be quantified or isolated.

Probing the thermodynamics and kinetics of ethylene carbonate reduction at the electrode-electrolyte interface with molecular simulations

Understanding the formation of the solid-electrolyte interphase (SEI) in lithium-ion batteries is an ongoing area of research due to its high degree of complexity and the difficulties encountered by experimental studies. Herein, we investigate the initial stage of SEI growth, the reduction reaction of ethylene carbonate (EC), from both a thermodynamic and a kinetic approach with theory and molecular simulations. We employed both the potential distribution theorem and the Solvation Method based on Density (SMD) to EC solvation for the estimation of reduction potentials of Li⁺, EC, and Li⁺-solvating EC (s-EC) as well as reduction rate constants of EC and s-EC. We find that solvation effects greatly influence these quantities of interest, particularly the Li⁺/Li reference electrode potential in EC solvent. Furthermore, we also compute the inner- and outer-sphere reorganization energies for both EC and s-EC at the interface of liquid EC and a hydroxyl-terminated graphite surface, where total reorganization energies are predicted to be 76.6 and 88.9 kcal/mol, respectively. With the computed reorganization energies, we estimate reduction rate constants across a range of overpotentials and show that EC has a larger electron transfer rate constant than s-EC at equilibrium, despite s-EC being more thermodynamically favorable. Overall, this manuscript demonstrates how ion solvation effects largely govern the prediction of reduction potentials and electron transfer rate constants at the electrode-electrolyte interface.

Modeling galvanostatic charge-discharge of nanoporous supercapacitors

Molecular modeling has been considered indispensable in studying the energy storage of supercapacitors at the atomistic level. The constant potential method (CPM) allows the electric potential to be kept uniform in the electrode, which is essential for a realistic description of the charge repartition and dynamics process in supercapacitors. However, previous CPM studies have been limited to the potentiostatic mode. Although widely adopted in experiments, the galvanostatic mode has rarely been investigated in CPM simulations because of a lack of effective methods. Here we develop a modeling approach to simulating the galvanostatic charge-discharge process of supercapacitors under constant potential. We show that, for nanoporous electrodes, this modeling approach can capture experimentally consistent dynamics in supercapacitors. It can also delineate, at the molecular scale, the hysteresis in ion adsorption-desorption dynamics during charging and discharging. This approach thus enables the further accurate modeling of the physics and electrochemistry in supercapacitor dynamics.

Cation- and pH-Dependent Hydrogen Evolution and Oxidation Reaction Kinetics

The production of molecular hydrogen by catalyzing water splitting is central to achieving the decarbonization of sustainable fuels and chemical transformations. In this work, a series of structure-making/breaking cations in the electrolyte were investigated as spectator cations in hydrogen evolution and oxidation reactions (HER/HOR) in the pH range of 1 to 14, whose kinetics was found to be altered by up to 2 orders of magnitude by these cations. The exchange current density of HER/HOR was shown to increase with greater structure-making tendency of cations in the order of Cs⁺ < Rb⁺ < K⁺ < Na⁺ < Li⁺, which was accompanied by decreasing reorganization energy from the Marcus-Hush-Chidsey formalism and increasing reaction entropy. Invoking the Born model of reorganization energy and reaction entropy, the static dielectric constant of the electrolyte at the electrified interface was found to be significantly lower than that of bulk, decreasing with the structure-making tendency of cations at the negatively charged Pt surface. The physical origin of cation-dependent HER/HOR kinetics can be rationalized by an increase in concentration of cations on the negatively charged Pt surface, altering the interfacial water structure and the H-bonding network, which is supported by classical molecular dynamics simulation and surface-enhanced infrared absorption spectroscopy. This work highlights immense opportunities to control the reaction rates by tuning interfacial structures of cation and solvents.