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.

Adsorption of Tunable aryl alkyl ionic liquids (TILs) on the graphene and Defective graphene nanosheets: A DFT Study

Optical and electronic characteristics of the graphene nanosheets (GNS) could be altered by some structural defects such as double-vacancy and Stone-Wales ones. The physisorption manner of [MPI][BF4], [MPT1][BF4], [MPT2][BF4], and [MPTT][BF4] ionic liquids on intact and defective GNS surfaces were investigated using M06-2X/cc-pVDZ computational method. Capability for adsorption on the DV and SW graphene surfaces by TILs is increased by about 1.0-4.3 and 0.4-2.0 kcal/mol respectively. The electrostatic potential of the GNS-DV surface is more negative than the GNS-SW one which enables it to interact with cation parts of the adsorbed TILs so extensively. The highest adsorption energy belongs to the [MPI][BF4]/GNS-DV system. Adsorption of the TILs on the GNS surfaces leads to a decrease in the energy of the LUMO molecular orbital as well as their energy gap of them. Results revealed that the electrical conductivity, as well as absorption spectra of the GNS surfaces, are affected by TILs adsorption and defect nature.

Toward a New Generation of Fire-Safe Energy Storage Devices: Recent Progress on Fire-Retardant Materials and Strategies for Energy Storage Devices

Over the last few decades, tremendous progress has been achieved in the development of advanced materials for energy storage devices. These achievements have largely enabled the adoption and transition to key technologies such as mobile phones, electric vehicles, and internet of things. However, the recent surge in fire accidents and explosions emanating from energy storage devices have been closely associated with the highly flammable components that make up these devices which have often led to the loss of life and property. Therefore, replacing flammable materials with fire retardant materials has been recognized as the critical solution to the ever-growing fire problem in these devices. This review summarizes the progress achieved so far in the field of fire retardant materials for energy storage devices. Finally, a perspective on the current state of the art is provided, and a future outlook for these fire-retardant materials, strategies, and new characterization methods is discussed.

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...

Carbon-carbon supercapacitors: Beyond the average pore size or how electrolyte confinement and inaccessible pores affect the capacitance

Carbon-carbon supercapacitors are high power electrochemical energy storage systems, which store energy through reversible ion adsorption at the electrode-electrolyte interface. Due to the complex structure of the porous carbons used as electrodes, extracting structure-property relationships in these systems remains a challenge. In this work, we conduct molecular simulations of two model supercapacitors based on nanoporous electrodes with the same average pore size, a property often used when comparing porous materials, but different morphologies. We show that the carbon with the more ordered structure, and a well defined pore size, has a much higher capacitance than the carbon with the more disordered structure and a broader pore size distribution. We analyze the structure of the confined electrolyte and show that the ions adsorbed in the ordered carbon are present in larger quantities and are also more confined than for the disordered carbon. Both aspects favor a better charge separation and thus a larger capacitance. In addition, the disordered electrodes contain a significant amount of carbon atoms, which are never in contact with the electrolyte, carry a close to zero charge, and are thus not involved in the charge storage. The total quantities of adsorbed ions and degrees of confinement do not change much with the applied potential, and as such, this work opens the door to computationally tractable screening strategies.

Pimples reduce and dimples enhance flat dielectric surface image repulsion

In solid-liquid, or liquid-liquid, interfaces with dielectric contrast, charged particles interact with the induced polarization charge of the interface. These interactions contribute to an effective self-energy of the bulk ions and mediate ion-ion interactions. For flat interfaces, the self-energy and the mediated interactions are neatly constructed by the image charge method. For other geometries, explicit results are scarce and the problem must be treated via approximations or direct computation. The case of interfaces with roughness is of great practical importance. This article provides analytical results, valid to first-order in perturbation theory, for the self-energy of particles near rough substrates. Explicit formulas are provided for the case of a sinusoidal deformation of a flat surface. Generic deformations can be treated by superposition. In addition to results for the self-energy, the surface polarization charge is presented as a quadrature. The interaction between an ion and the deformed surface is modified by the change in relative distance as well as by the local curvature of the surface. Solid walls, with a lower dielectric constant than the liquid, repel all ions. We show that the repulsion is reduced by local convexity and enhanced by concavity; dimples are more repulsive than pimples.

Nanoconfinement between Graphene Walls Suppresses the Near-Wall Diffusion of the Ionic Liquid [BMIM][PF6]

We identify two distinct regimes for the diffusion of the ionic liquid [BMIM][PF₆] confined between parallel graphene walls using molecular dynamics simulations. Within 2 nm of the wall, the cations and anions form a well-defined layered structure. In this region, the in-plane diffusion coefficients are suppressed when compared to their bulk values and increase monotonically with the distance away from the wall. Beyond 2 nm from the wall, the density profile and in-plane diffusion coefficients recover their bulk values. The channel-averaged in-plane diffusion coefficients increase monotonically with wall separation and recover the bulk values at a separation of 15 nm. A simple semianalytical model is proposed that mirrors this trend. The results also highlight the importance of applying a finite-size correction to molecular dynamics-predicted diffusion coefficients of confined liquids, which may otherwise be unusually larger than their bulk values.

A QM/MD Coupling Method to Model the Ion-Induced Polarization of Graphene

We report a new Quantum Mechanical/Molecular Dynamics (QM/MD) simulation loop to model the coupling between the electron and atom dynamics in solid/liquid interfacial systems. The method can describe simultaneously both the quantum mechanical surface polarizability emerging from the proximity to the electrolyte and the electrolyte structure and dynamics. In the current setup, Density Functional Tight Binding calculations for the electronic structure calculations of the surface are coupled with classical molecular dynamics to simulate the electrolyte solution. The reduced computational cost of the QM part makes the coupling with a classical simulation engine computationally feasible and allows simulation of large systems for hundreds of nanoseconds. We tested the method by simulating both a noncharged graphene flake and a noncharged and charged infinite graphene sheet immersed in an NaCl electrolyte solution. We found that, when no bias is applied, ions preferentially remained in solution, and only cations are mildly attracted to the surface of the graphene. This preferential adsorption of cations vs anions seems to persist also when the surface is moderately charged and rules out any substantial ions/surface charge transfer.

Computational Study of the Properties of Acetonitrile/Water-in-Salt Hybrid Electrolytes as Electrolytes for Supercapacitors

Normal and water-in-salt Li-bis(trifluoromethane) sulfonimide anion-based electrolytes were modeled using atomistic molecular dynamics simulations. Their acetonitrile (ACN) mixtures, in various concentrations, were also studied to evaluate the impact of a cosolvent on the structural, dynamical, and electrical properties of the electrolytes using liquid electrolyte and supercapacitor models. Our simulations for pure and ACN-based electrolytes revealed a drastic difference that exists between normal electrolytes and water-in-salt electrolytes and a systematic reduction of the diffusion of species by approximately a factor of 2 because of the ACN impact. Electrolytic cells for each electrolyte were built with graphene as the electrode. Our results for capacitance reveal an asymmetry between the electrode capacitances, with negative electrode capacitance systematically higher than those of the positive electrode. The total capacitance of the electrode exhibited negligible variations regardless of the concentration and composition of the electrolyte.

Understanding the charging dynamics of an ionic liquid electric double layer capacitor via molecular dynamics simulations

We investigate the charging phenomena of an electric double layer capacitor (EDLC) by conducting both equilibrium and non-equilibrium molecular dynamics (MD) simulations. A graphene electrode and 1-ethyl-3-methylimidazolium thiocyanate ([EMIM]+[SCN]-) ionic liquid were used as a system for the EDLC. We clarify the ionic layer structure and show that an abrupt change of the ionic layers leads to a high differential capacitance of the EDLC. The charging simulations reveal that the charging dynamics of the EDLC is highly dependent on the rearrangement of the ionic layer structure. Particularly, the electrode charge during the charging process is consistent with the perpendicular displacement of ionic liquid molecules. From this property, we analyze the contribution of each molecular ion to the electrode charge stored during charging. Charging of the EDLC is largely dependent on the desorption of the co-ions from the electrode rather than the adsorption of the counter-ions. In addition, the contribution of bulk ions to the charge stored in the EDLC is as important as that of ions adjacent to the electrode surface contrary to the conventional viewpoint. From these results, we identify the charging mechanism of the EDLC and discuss the relevance to experimental results. Our findings in the present study are expected to play an important role in designing an efficient EDLC with a novel perspective on the charging of the EDLC.

Investigating the asymmetry in the EDL response of C60/graphene supercapacitors

Molecular dynamics simulations were employed to model C₆₀/graphene composite electrodes that can expand the effective area and performance of supercapacitors.

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.

Active Mechanism of the Interphase Film-Forming Process for an Electrolyte Based on a Sulfolane Solvent and a Chelato-Borate Complex

Electrolytes based on sulfolane (SL) solvents and lithium bis(oxalato)borate (LiBOB) chelato-borate complexes have been reported many times for use in advanced lithium-ion batteries due to their many advantages. This study aims to clarify the active mechanism of the interphase film-forming process to optimize the properties of these batteries by experimental analysis and theoretical calculations. The results indicate that the self-repairing film-forming process during the first cycle is divided into three stages: the initial film formation with an electric field force of ∼1.80 V, the further growth of the preformation solid electrolyte interphase (SEI) film at ∼1.73 V, and the final formation of a complete SEI film at a potential below 0.7 V. Additionally, we can deduce that the decomposition of LiBOB and SL occurs throughout nearly the entire process of the formation of the SEI film. The decomposition product of BOB^- anions tends to form films with an irregular structure, whereas the decomposition product of SL is in favor of the formation of a uniform SEI film.

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.

Kinetic-Dominated Charging Mechanism within Representative Aqueous Electrolyte-based Electric Double-Layer Capacitors

The chemical nature of electrolytes has been demonstrated to play a pivotal role in the charge storage of electric double-layer capacitors (EDLCs), whereas primary mechanisms are still partially resolved but controversial. In this work, a systematic exploration into EDL structures and kinetics of representative aqueous electrolytes is performed with numerical simulation and experimental research. Unusually, a novel charging mechanism exclusively predominated by kinetics is recognized, going beyond traditional views of manipulating capacitances preferentially via interfacial structural variations. Specifically, strikingly distinctive EDL structures stimulated by diverse ion sizes, valences, and mixtures manifest a virtually identical EDL capacitance, where the dielectric nature of solvents attenuates ionic effects on electrolyte redistributions, in stark contradiction with solvent-free counterpart and traditional Helmholtz theory. Meanwhile, corresponding kinetics evolve conspicuously with ionic species, intimately correlated with ion-solvent interactions. The achieved mechanisms are subsequently illuminated by electrochemical measurements, highlighting the crucial interplay between ions and solvents in regulating EDLC performances.

The effect of defect types on the electronic and optical properties of graphene nanoflakes physisorbed by ionic liquids

Defect engineering and non-covalent interaction strategies allow for dramatically tuning the optoelectronic properties of graphene. Using ab initio density functional theory (M06-2X/cc-pVDZ), we find that the nature of defects on the graphene nanoflakes (GNFs) and the size of defective GNF (DGNF) surfaces affect the binding energy (ΔE_b) of ionic liquids (ILs) and the UV-Vis absorption spectra of DGNFIL complexes. Further, our results indicate that increasing the size of DGNFs affects the geometrical structure of the surfaces and increases the binding energy of ILs by about 10%. Analysis based on AIM and EDA shows that the interactions between ILs and DGNFs are non-covalent in nature (dispersion energy being dominant) and associated with charge transfer between the IL and nanoflakes. A comparison between the ΔE_b values of ILs on DGNFs, GNFs, and h-BN nanoflakes (h-BNNF) shows that the presence of defects on the GNF surfaces increases the binding energy values as follows: DGNFIL > pristine GNFIL > h-BNNFIL. Our calculations indicate that increasing the size of DGNF surfaces leads to a decrease in the HOMO-LUMO energy gap (E_g) of the DGNF surfaces. Orbital energy and density of state calculations show that the E_g of DV(SW)-GNFs decreases upon IL adsorption and their Fermi energy level is shifted depending on the type of IL, thus enabling better conductivity. Reactivity descriptors generally indicate that the chemical potential (μ) and chemical hardness (η) of nanoflakes decrease upon IL adsorption, whereas the electrophilicity index (ω) increases. The UV-Vis absorption spectrum of DV-GNF and SW-GNF shows four bands in the visible spectrum which correspond to π → π* transitions with the absorption bands of SW-GNF appearing at higher wavelengths than those of DV-GNF. The most intense absorption bands in DV-GNF (λ = 348 nm) and SW-GNF (λ = 375 nm) are associated with electronic transitions HOMO-1 → LUMO+2 and HOMO → LUMO+1, respectively. In addition, these absorption bands undergo a red-shift by both increasing the size of the DV(SW)-GNF surfaces and IL adsorption. We also observe that the energy gaps and absorption spectra can be altered by varying the defect types and the type of IL adsorbate, where the defect types affect the spectral shapes of the bands and adsorbates at the first absorption peak, thus having potential application for light-emitting devices.

A review of the structure and dynamics of nanoconfined water and ionic liquids via molecular dynamics simulation

During the past decade, the research on fluids in nanoconfined geometries has received considerable attention as a consequence of their wide applications in different fields. Several nanoconfined systems such as water and ionic liquids, together with an equally impressive array of nanoconfining media such as carbon nanotube, graphene and graphene oxide have received increasingly growing interest in the past years. Water is the first system that has been reviewed in this article, due to its important role in transport phenomena in environmental sciences. Water is often considered as a highly nanoconfined system, due to its reduction to a few layers of water molecules between the extended surface of large macromolecules. The second system discussed here is ionic liquids, which have been widely studied in the modern green chemistry movement. Considering the great importance of ionic liquids in industry, and also their oil/water counterpart, nanoconfined ionic liquid system has become an important area of research with many fascinating applications. Furthermore, the method of molecular dynamics simulation is one of the major tools in the theoretical study of water and ionic liquids in nanoconfinement, which increasingly has been joined with experimental procedures. In this way, the choice of water and ionic liquids in nanoconfinement is justified by applying molecular dynamics simulation approaches in this review article.

The diffusion behavior and capacitance of tetraethylammonium/tetrafluoroborate ions in acetonitrile with different molar concentrations: a molecular dynamics study

A molecular dynamics (MD) simulation with the optimized potentials for liquid simulations-all atom (OPLS-AA) force field was carried out to investigate the dynamic behaviors of organic electrolyte molecules between a graphite cathode and anode.

Electric potential calculation in molecular simulation of electric double layer capacitors

For the molecular simulation of electric double layer capacitors (EDLCs), a number of methods have been proposed and implemented to determine the one-dimensional electric potential profile between the two electrodes at a fixed potential difference. In this work, we compare several of these methods for a model LiClO4-acetonitrile/graphite EDLC simulated using both the traditional fixed-charged method (FCM), in which a fixed charge is assigned a priori to the electrode atoms, or the recently developed constant potential method (CPM) (2007 J. Chem. Phys. 126 084704), where the electrode charges are allowed to fluctuate to keep the potential fixed. Based on an analysis of the full three-dimensional electric potential field, we suggest a method for determining the averaged one-dimensional electric potential profile that can be applied to both the FCM and CPM simulations. Compared to traditional methods based on numerically solving the one-dimensional Poisson's equation, this method yields better accuracy and no supplemental assumptions.

Computer Simulation Study of Graphene Oxide Supercapacitors: Charge Screening Mechanism

Graphene oxide supercapacitors in the parallel plate configuration are studied via molecular dynamics (MD) simulations. The full range of electrode oxidation from 0 to 100% is examined by oxidizing the graphene surface with hydroxyl groups. Two different electrolytes, 1-ethyl-3-methylimidazolium tetrafluoroborate (EMI(+)BF4(-)) as an ionic liquid and its 1.3 M solution in acetonitrile as an organic electrolyte, are considered. While the area-specific capacitance tends to decrease with increasing electrode oxidation for both electrolytes, its details show interesting differences between the organic electrolyte and ionic liquid, including the extent of decrease. For detailed insight into these differences, the screening mechanisms of electrode charges by electrolytes and their variations with electrode oxidation are analyzed with special attention paid to the aspects shared by and the contrasts between the organic electrolyte and ionic liquid.

Heterogeneous dynamics and its length scale in simple ionic liquid models: a computational study

We numerically investigate the dynamic heterogeneity and its length scale found in coarse-grained ionic liquid model systems.

Electrolytes between dielectric charged surfaces: Simulations and theory

We present a simulation method to study electrolyte solutions in a dielectric slab geometry using a modified 3D Ewald summation. The method is fast and easy to implement, allowing us to rapidly resum an infinite series of image charges. In the weak coupling limit, we also develop a mean-field theory which allows us to predict the ionic distribution between the dielectric charged plates. The agreement between both approaches, theoretical and simulational, is very good, validating both methods. Examples of ionic density profiles in the strong electrostatic coupling limit are also presented. Finally, we explore the confinement of charge asymmetric electrolytes between neutral surfaces.

Computational chemistry for graphene-based energy applications: progress and challenges

Research in graphene-based energy materials is a rapidly growing area. Many graphene-based energy applications involve interfacial processes. To enable advances in the design of these energy materials, such that their operation, economy, efficiency and durability is at least comparable with fossil-fuel based alternatives, connections between the molecular-scale structure and function of these interfaces are needed. While it is experimentally challenging to resolve this interfacial structure, molecular simulation and computational chemistry can help bridge these gaps. In this Review, we summarise recent progress in the application of computational chemistry to graphene-based materials for fuel cells, batteries, photovoltaics and supercapacitors. We also outline both the bright prospects and emerging challenges these techniques face for application to graphene-based energy materials in future.

Effects of Oxygen-Containing Functional Groups on Supercapacitor Performance

Molecular dynamics (MD) simulations of the interface between graphene and the ionic liquid 1-butyl-3-methylimidazolium trifluoromethanesulfonate (BMIM OTf) were carried out to gain molecular-level insights into the performance of graphene-based supercapacitors and, in particular, determine the effects of the presence of oxygen-containing defects at the graphene surface on their integral capacitance. The MD simulations predict that increasing the surface coverage of hydroxyl groups negatively affects the integral capacitance, whereas the effect of the presence of epoxy groups is much less significant. The calculated variations in capacitance are found to be directly correlated to the interfacial structure. Indeed, hydrogen bonding between hydroxyl groups and SO3 moieties prevents BMIM(+) and OTf(-) ions from interacting favorably in the interfacial layer and restrains the orientation and mobility of OTf(-) ions, thereby reducing the interfacial permittivity of the ionic liquid. The results of the simulations can facilitate the rational design of electrode materials for supercapacitors.

Unusual effects of solvent polarity on capacitance for organic electrolytes in a nanoporous electrode

The interplay between ions and solvent molecules inside the nanoporous electrodes of a supercapacitor has not been well understood but could be a fertile ground for new insights into the device's performance. By tuning the dipole moment of the solvent in an organic electrolyte, we find, from classical density functional theory calculations, pronounced oscillation of capacitance with the pore size for a moderately to weakly polar solvent. A quantitative analysis of the electric-double-layer (EDL) structure indicates that the capacitance oscillation shares a similar physical origin to that of an ionic liquid electrolyte: the oscillatory behavior arises from the formation of alternating layers of counterions and coions near strongly charged surfaces. More interestingly, we find that in the large-pore region, the capacitance versus the pore size has a volcano-shaped trend; in other words, there exists a solvent dipole moment that yields a maximal capacitance. These theoretical predictions can be validated with future experiments and highlight the great potential in tuning the organic solvent to achieve optimal performance of EDL capacitors.

Carbons and electrolytes for advanced supercapacitors

Electrical energy storage (EES) is one of the most critical areas of technological research around the world. Storing and efficiently using electricity generated by intermittent sources and the transition of our transportation fleet to electric drive depend fundamentally on the development of EES systems with high energy and power densities. Supercapacitors are promising devices for highly efficient energy storage and power management, yet they still suffer from moderate energy densities compared to batteries. To establish a detailed understanding of the science and technology of carbon/carbon supercapacitors, this review discusses the basic principles of the electrical double-layer (EDL), especially regarding the correlation between ion size/ion solvation and the pore size of porous carbon electrodes. We summarize the key aspects of various carbon materials synthesized for use in supercapacitors. With the objective of improving the energy density, the last two sections are dedicated to strategies to increase the capacitance by either introducing pseudocapacitive materials or by using novel electrolytes that allow to increasing the cell voltage. In particular, advances in ionic liquids, but also in the field of organic electrolytes, are discussed and electrode mass balancing is expanded because of its importance to create higher performance asymmetric electrochemical capacitors.