Influence of solvation on the structural and capacitive properties of electrical double layer capacitors

Organic Solvent Boosts Charge Storage and Charging Dynamics of Conductive MOF Supercapacitors

Conductive metal-organic frameworks (c-MOFs) and ionic liquids (ILs) have emerged as auspicious combinations for high-performance supercapacitors. However, the nanoconfinement from c-MOFs and high viscosity of ILs slow down the charging process. This hindrance can, however, be resolved by adding solvent. Here, constant-potential molecular simulations are performed to scrutinize the solvent impact on charge storage and charging dynamics of MOF-IL-based supercapacitors. Conditions for >100% enhancement in capacity and ≈6 times increase in charging speed are found. These improvements are confirmed by synthesizing near-ideal c-MOFs and developing multiscale models linking molecular simulations to electrochemical measurements. Fundamentally, the findings elucidate that the solvent acts as an "ionophobic agent" to induce a substantial enhancement in charge storage, and as an "ion traffic police" to eliminate convoluted counterion and co-ion motion paths and create two distinct ion transport highways to accelerate charging dynamics. This work paves the way for the optimal design of MOF 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.

Predictive Molecular Models for Charged Materials Systems: From Energy Materials to Biomacromolecules

Electrostatic interactions play a dominant role in charged materials systems. Understanding the complex correlation between macroscopic properties with microscopic structures is of critical importance to develop rational design strategies for advanced materials. But the complexity of this challenging task is augmented by interfaces present in the charged materials systems, such as electrode-electrolyte interfaces or biological membranes. Over the last decades, predictive molecular simulations that are founded in fundamental physics and optimized for charged interfacial systems have proven their value in providing molecular understanding of physicochemical properties and functional mechanisms for diverse materials. Novel design strategies utilizing predictive models have been suggested as promising route for the rational design of materials with tailored properties. Here, an overview of recent advances in the understanding of charged interfacial systems aided by predictive molecular simulations is presented. Focusing on three types of charged interfaces found in energy materials and biomacromolecules, how the molecular models characterize ion structure, charge transport, morphology relation to the environment, and the thermodynamics/kinetics of molecular binding at the interfaces is discussed. The critical analysis brings two prominent field of energy materials and biological science under common perspective, to stimulate crossover in both research field that have been largely separated.

Assessing the impact of valence asymmetry in ionic solutions and its consequences on the performance of supercapacitors

Molecular dynamics simulations were performed to describe the properties of hypothetical salt electrolytic solutions. The main focus of this work is the valence asymmetry, which in recent years has been considered an important aspect in the physical chemistry of aqueous electrolytes. In general, our results show that the structural, energetic, and dynamic properties respond differently to the asymmetry of ionic solutions, but in all cases, appreciable changes were observed. Graphene supercapacitors based on the investigated electrolytes were studied in light of their electrostatic properties. We observed that the electrode capacitances, positive and negative, were greatly influenced by the presence of cations in the electrical double layer of the negative electrode and by the absence of these cations, in the double layer of the positive electrode. In general, we assess that quantitative variations due to valence asymmetry may indeed be an important factor for the development of new and more efficient electrolytes.

Microscopic Simulations of Electrochemical Double-Layer Capacitors

Electrochemical double-layer capacitors (EDLCs) are devices allowing the storage or production of electricity. They function through the adsorption of ions from an electrolyte on high-surface-area electrodes and are characterized by short charging/discharging times and long cycle-life compared to batteries. Microscopic simulations are now widely used to characterize the structural, dynamical, and adsorption properties of these devices, complementing electrochemical experiments and in situ spectroscopic analyses. In this review, we discuss the main families of simulation methods that have been developed and their application to the main family of EDLCs, which include nanoporous carbon electrodes. We focus on the adsorption of organic ions for electricity storage applications as well as aqueous systems in the context of blue energy harvesting and desalination. We finally provide perspectives for further improvement of the predictive power of simulations, in particular for future devices with complex electrode compositions.

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.

Origin of Enhanced Performance in Nanoporous Electrical Double Layer Capacitors: Insights on Micropore Structure and Electrolyte Composition from Molecular Simulations

We explore the effect of solvation and micropore structure on the energy storage performance of electrical double layer capacitors using constant potential molecular dynamics simulations of realistically modeled nanoporous carbon electrodes and ionic liquid/organic solvent mixtures. We show that the time-dependent charging profiles of electrodes with larger pores reach the plateau regime faster, while the charging time has a nonmonotonic dependence on ion concentration, mirroring the composition dependence of bulk electrolyte conductivity. When the average pore size of the electrode is similar to or slightly larger than the size of a solvated ion, the solvation enhances ion electrosorption into nanopores by disrupting anion-cation coordination and decreasing the barrier to counterion penetration while blocking the co-ions. In these systems, areal capacitance exhibits a significant nonmonotonic dependence on ion concentration, in which capacitance increases with the introduction of solvent in the concentrated regime followed by a decrease with further dilution. This gives rise to a maximum in capacitance at intermediate dilution levels. When pores are significantly larger than solvated ions, capacitance maximum weakens and eventually disappears. These findings provide novel insights on the combined effect of electrolyte composition and electrode pore size on the charging kinetics and equilibrium behavior of realistically modeled electrical double layer capacitors. Generalization of the approach developed here can facilitate the rational optimization of material properties for electrical double layer capacitor applications.

Multi-scale simulation of the adsorption of lithium ion on graphite surface: From quantum Monte Carlo to molecular density functional theory

The structure of the double-layer formed at the surface of carbon electrodes is governed by the interactions between the electrode and the electrolyte species. However, carbon is notoriously difficult to simulate accurately, even with well-established methods such as electronic density functional theory and molecular dynamics. Here, we focus on the important case of a lithium ion in contact with the surface of graphite, and we perform a series of reference quantum Monte Carlo calculations that allow us to benchmark various electronic density functional theory functionals. We then fit an accurate carbon–lithium pair potential, which is used in molecular density functional theory calculations to determine the free energy of the adsorption of the ion on the surface in the presence of water. The adsorption profile in aqueous solution differs markedly from the gas phase results, which emphasize the role of the solvent on the properties of the double-layer.

Elucidating Curvature-Capacitance Relationships in Carbon-Based Supercapacitors

Nanoscale surface curvatures, either convex or concave, strongly influence the charging behavior of supercapacitors. Rationalizing individual influences of electrode atoms to the capacitance is possible by interpreting distinct elements of the charge-charge covariance matrix derived from individual charge variations of the electrode atoms. An ionic liquid solvated in acetonitrile and confined between two electrodes, each consisting of three undulated graphene layers, serves as a demonstrator to illustrate pronounced and nontrivial features of the capacitance with respect to the electrode curvature. In addition, the applied voltage determines whether a convex or concave surface contributes to increased capacitance. While at lower voltages capacitance variations are in general correlated with ion number density variations in the double layer formed in the concave region of the electrode, for certain electrode designs a surprisingly strong contribution of the convex part to the differential capacitance is found both at higher and lower voltages.

Structure and Capacitance of Electrical Double Layers in Tricationic Ionic Liquids with Organic Solvents

Tricationic ionic liquid (TIL) electrolytes have been successfully employed in supercapacitors with graphene electrodes, but the low power density of the TILs-based supercapacitors caused by strong cations-anions associations requires enhancement by adding organic solvents to the liquid. In this paper, the role of the solvents acetonitrile (ACN) and ethylene carbonate (EC) on the ion diffusion, the conductivity of the TIL [C₆(mim)₃](Tf₂N)₃, and the structures and the capacitances of the electrical double layers (EDLs) in TIL/ACN and TIL/EC electrolytes were probed by molecular dynamics (MD) simulations. The results indicate that adding organic solvents to the liquid significantly reduces interactions between ions, thereby greatly improving the ion diffusion coefficients and the conductivity of the TIL, and the maximum conductivity is found at the 0.55 M TIL/ACN electrolyte concentration. Moreover, the reduced packing of counterions and the strong expulsion of coions near charged electrodes are observed in the organic electrolytes, especially in the TIL/EC electrolyte. Further analyses on EDLs affirm that the asymmetric camel-shaped differential capacitance-voltage (C-V) curve in the pure TIL electrolyte is weakly changed by the solvent ACN or EC. Besides, the EDL capacitance in the TIL-based hybrid electrolytes is improved slightly by the organic solvents. Comparing the integral capacitances in TIL/ACN and TIL/EC with different solvent contents, it is found that reducing the solvent polarity may be more beneficial to promote the EDL capacitance. Comprehensively, in this work, the 0.55 M TIL/ACN electrolyte is the optimal choice for the high-performance supercapacitor. Hence, solvating TIL electrolytes in supercapacitors by suitable solvents can effectively enhance the power density without compromising energy density.

Mesoscopic simulations of the in situ NMR spectra of porous carbon based supercapacitors: electronic structure and adsorbent reorganisation effects

In situ NMR spectroscopy is a powerful technique to investigate charge storage mechanisms in carbon-based supercapacitors thanks to its ability to distinguish ionic and molecular species adsorbed in the porous electrodes from those in the bulk electrolyte. The NMR peak corresponding to the adsorbed species shows a clear change of chemical shift as the applied potential difference is varied. This variation in chemical shift is thought to originate from a combination of ion reorganisation in the pores and changes in ring current shifts due to the changes of electronic density in the carbon. While previous Density Functional Theory calculations suggested that the electronic density has a large effect, the relative contributions of these two effects is challenging to untangle. Here, we use mesoscopic simulations to simulate NMR spectra and investigate the relative importance of ion reorganisation and ring currents on the resulting chemical shift. The model is able to predict chemical shifts in good agreement with NMR experiments and indicates that the ring currents are the dominant contribution. A thorough analysis of a specific electrode/electrolyte combination for which detailed NMR experiments have been reported allows us to confirm that local ion reorganisation has a very limited effect but the relative quantities of ions in pores of different sizes, which can change upon charging/discharging, can lead to a significant effect. Our findings suggest that in situ NMR spectra of supercapacitors may provide insights into the electronic structure of carbon materials in the future.

Structure of cholinium glycinate biocompatible ionic liquid at graphite electrode interface

We use constant potential molecular dynamics simulations to investigate the interfacial structure of the cholinium glycinate biocompatible ionic liquid (bio-IL) sandwiched between graphite electrodes with varying potential differences. Through number density profiles, we observe that the cation and anion densities oscillate up to ∼1.5 nm from the nearest electrode. The range of these oscillations does not change significantly with increasing electrode potential. However, the amplitudes of the cation (anion) density oscillations show a notable increase with increasing potential at the negative (positive) electrode. At higher potential differences, the bulkier N(CH₃)₃CH₂ group of cholinium cations ([Ch]⁺) overcomes the steric barrier and comes closer to the negative electrode as compared to oxygen atom (O_[Ch]⁺ ). We observe an increase in the interaction between O_[Ch]⁺ and the positive electrode with a decrease in the distance between them on increasing the potential difference. We also observe hydrogen bonding between the hydroxyl group of [Ch]⁺ cations and oxygens of glycinate anions through the simulated tangential radial distribution function. Orientational order parameter analysis shows that the cation (anion) prefers to align parallel to the negative (positive) electrode at higher applied potential differences. Charge density profiles show a positive charge density peak near the positive electrode at all the potential differences because of the presence of partially positive charged hydrogen atoms of cations and anions. The differential capacitance (C_d) of the bio-IL shows two constant regimes, one for each electrode. The magnitude of these C_d values clearly suggests potential application of such bio-ILs as promising battery electrolytes.

From Water Solutions to Ionic Liquids with Solid State Nanopores as a Perspective to Study Transport and Translocation Phenomena

Solid state nanopores are single-molecular devices governed by nanoscale physics with a broad potential for technological applications. However, the control of translocation speed in these systems is still limited. Ionic liquids are molten salts which are commonly used as alternate solvents enabling the regulation of the chemical and physical interactions on solid-liquid interfaces. While their combination can be challenging to the understanding of nanoscopic processes, there has been limited attempts on bringing these two together. While summarizing the state of the art and open questions in these fields, several major advances are presented with a perspective on the next steps in the investigations of ionic-liquid filled nanopores, both from a theoretical and experimental standpoint. By analogy to aqueous solutions, it is argued that ionic liquids and nanopores can be combined to provide new nanofluidic functionalities, as well as to help resolve some of the pertinent problems in understanding transport phenomena in confined ionic liquids and providing better control of the speed of translocating analytes.

Molecular Simulation of Electrode-Solution Interfaces

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

Atomic structure and electrical property of ionic liquids at the MoS2 electrode with varying interlayer spacing

Understanding the structure and properties at the electrolyte-electrode interface is vital for the rational design of the supercapacitors or other electrochemical devices. In this work, we explored the influence of interlayer spacing of the MoS₂ electrode on the interfacial structure and electrical properties of sodium-ionic liquids (ILs) electrolytes via performing the all-atom molecular dynamics simulations. From the number density, charge density, and electrical potential distribution near the surface, the Mo- and S-terminal edges possess positive and negative features when the interlayer spacing is less than 8.5 Å. Meanwhile, the strength of the first density layer of ILs increases with the increase of the interlayer spacing of MoS₂ for both Mo- and S- terminal surfaces in the neutral or charging state. Furthermore, the coordination number of sodium ion at the electrode surface was analyzed, and it was shown that the S-terminal surface has a larger coordination number than that on the Mo-terminal surface. Interestingly, the coordination number of MoS₂ with the interlayer spacing of 8.0 Å is the lowest in the ranges of 6.5~8.5 Å. The electrolyte's charge screening factor also reflects the opposite electrical state of Mo- and S-terminal surfaces and weakens with increasing the interlayer spacing and surface charge density. The obtained understanding of ILs at electrode interfaces with different layer spacings in this work will provide insight into the molecular mechanisms of ILs-based sodium supercapacitors or other electrochemical devices in critical chemical engineering processes.

Exploring doped or vacancy-modified graphene-based electrodes for applications in asymmetric supercapacitors

We report here density functional theory calculations and molecular dynamics atomistic simulations to determine the total capacitance of graphene-modified supercapacitors.

Charge fluctuations from molecular simulations in the constant-potential ensemble

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

Carbons with Regular Pore Geometry Yield Fundamental Insights into Supercapacitor Charge Storage

We conduct molecular dynamics simulations of electrical double-layer capacitors (EDLCs) using a library of ordered, porous carbon electrode materials called zeolite templated carbons (ZTCs). The well-defined pore shapes of the ZTCs enable us to determine the influence of pore geometry on both charging dynamics and charge storage mechanisms in EDLCs, also referred to as supercapacitors. We show that charging dynamics are negatively correlated with the pore-limiting diameter of the electrode material and display signatures of both progressive charging and ion trapping. However, the equilibrium capacitance, unlike charging dynamics, is not strongly correlated to commonly used, purely geometric descriptors such as pore size. Instead, we find a strong correlation of capacitance to the charge compensation per carbon (CCpC), a descriptor we define in this work as the average charge of the electrode atoms within the coordination shell of a counterion. A high CCpC indicates efficient charge storage, as the strong partial charges of the electrode are able to screen counterion charge, enabling higher ion loading and thus more charge storage within the electrode at a fixed applied voltage. We determine that adsorption sites with a high CCpC tend to be found within pockets with a smaller radius of curvature, where the counterions are able to minimize their distance with multiple points on the electrode surface, and therefore induce stronger local partial charges.

Structural Forces in Mixtures of Ionic Liquids with Organic Solvents

Using molecular dynamics simulations, we study the impact of electrode charging and addition of solvent (acetonitrile, ACN) on structural forces of the BMIM PF₆ ionic liquid (IL) confined by surfaces at nanometer separations. We establish relationships between the structural forces and the microscopic structure of the confined liquid. Depending on the structural arrangements of cations and anions across the nanofilm, the load-induced squeeze-out of liquid layers occurs via one-layer or bilayer steps. The cations confined between charged plates orient with their aliphatic chain perpendicular to the surface planes and link two adjacent IL layers. These structures facilitate the squeeze-out of single layers. For both pure IL and IL-ACN mixtures, we observe a strong dependence of nanofilm structure on the surface charge density, which affects the simulated pressure-displacement curves. Addition of solvent to the IL modifies the layering in the confined film. At high electrode charges and high dilution of IL (below 10% molar fraction), the layered structure of the nanofilm is less well defined. We predict a change in the squeeze-out mechanism under pressure, from a discontinuous squeeze-out (for high IL concentrations) to an almost continuous one (for low IL concentrations). Importantly, our simulations show that charged electrodes are coated with ions even at low IL concentrations. These ion-rich layers adjacent to the charged plate surfaces are not squeezed out even under very high normal pressures of ∼5 GPa. Hence, we demonstrate the high performance of IL-solvent mixtures to protect surfaces from wear and to provide lubrication at high loads.

Electrode polarization effects on interfacial kinetics of ionic liquid at graphite surface: An extended lagrangian-based constant potential molecular dynamics simulation study

Computational models including electrode polarization can be essential to study electrode/electrolyte interfacial phenomena more realistically. We present here a constant-potential classical molecular dynamics simulation method based on the extended Lagrangian formulation where the fluctuating electrode atomic charges are treated as independent dynamical variables. The method is applied to a graphite/ionic liquid system for the validation and the interfacial kinetics study. While the correct adiabatic dynamics is achieved with a sufficiently small fictitious mass of charge, static properties have been shown to be almost insensitive to the fictitious mass. As for the kinetics study, electrical double layer (EDL) relaxation and ion desorption from the electrode surface are considered. We found that the polarization slows EDL relaxation greatly whereas it has little impact on the ion desorption kinetics. The findings suggest that the polarization is essential to estimate the kinetics in nonequilibrium processes, not in equilibrium. © 2019 Wiley Periodicals, Inc.

Adding Solvent into Ionic Liquid-Gated Transistor: The Anatomy of Enhanced Gating Performance

Most studies of ionic liquid (IL)-gated field effect transistors (FETs) focus on the extremely large electric field and capacitance induced in liquid/solid interfaces and correspondingly the significantly enhanced carrier density in semiconductors, which can appreciably improve the gating performance. However, how to boost the switching speed, another key property of gating performance of FETs, has been rarely explored. In this work, the gating performance of molybdenum disulfide (MoS₂) FETs, gated by the mixtures of IL/organic solvent (1-butyl-3-methylimidazolium tetrafluoroborate/acetonitrile, [Bmim][BF₄]/ACN) at different ion concentrations, is investigated for both dynamic and static properties by a combination of molecular dynamics simulation and resistance network analysis. Results reveal that organic solvent can speed up the IL response time by a factor of about 40 times at the optimal ion concentration of 1.94 M, which is mainly attributed to the increased ionic conductivity of IL via the addition of organic solvent. Meanwhile, the surface charge distribution of MoS₂ becomes more homogenous after the addition of organic solvent, which increases the conductivity of MoS₂ by up to 2.4 times. Surprisingly, the optimal ion concentration for increased switching speed is nearly the same as that for achieving the highest MoS₂ conductivity. Thus, our findings provide a strategy to simultaneously improve the dynamic and static gating performance of IL-gated FETs as well as a modeling technique to screen out the ideal ion concentration.

Interfacial structure and structural forces in mixtures of ionic liquid with a polar solvent

Oscillatory and monotonic decay in mixtures of salt and solvent at interfaces with varying surface charge.

A coarse-grained polarizable force field for the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate

We present a coarse-grained polarizable molecular dynamics force field for the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIm][PF₆]). For the treatment of electronic polarizability, we employ the Drude model. Our results show that the new explicitly polarizable force field reproduces important static and dynamic properties such as mass density, enthalpy of vaporization, diffusion coefficients, or electrical conductivity in the relevant temperature range. In situations where an explicit treatment of electronic polarizability might be crucial, we expect the force field to be an improvement over non-polarizable models, while still profiting from the reduction of computational cost due to the coarse-grained representation.

Structure of electric double layers in capacitive systems and to what extent (classical) density functional theory describes it

Ongoing scientific interest is aimed at the properties and structure of electric double layers (EDLs), which are crucial for capacitive energy storage, water treatment, and energy harvesting technologies like supercapacitors, desalination devices, blue engines, and thermocapacitive heat-to-current converters. A promising tool to describe their physics on a microscopic level is (classical) density functional theory (DFT), which can be applied in order to analyze pair correlations and charge ordering in the primitive model of charged hard spheres. This simple model captures the main properties of ionic liquids and solutions and it predicts many of the phenomena that occur in EDLs. The latter often lead to anomalous response in the differential capacitance of EDLs. This work constructively reviews the powerful theoretical framework of DFT and its recent developments regarding the description of EDLs. It explains to what extent current approaches in DFT describe structural ordering and in-plane transitions in EDLs, which occur when the corresponding electrodes are charged. Further, the review briefly summarizes the history of modeling EDLs, presents applications, and points out limitations and strengths in present theoretical approaches. It concludes that DFT as a sophisticated microscopic theory for ionic systems is expecting a challenging but promising future in both fundamental research and applications in supercapacitive technologies.

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.

Computational Insights into Materials and Interfaces for Capacitive Energy Storage

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

Molecular Origin of Electric Double-Layer Capacitance at Multilayer Graphene Edges

Multilayer graphenes have been widely used as active materials for electric double-layer capacitors (EDLCs), where their numerous edges are demonstrated to play a crucial role in charge storage. In this work, the interfacial structure and capacitive behaviors of multilayer graphene edges with representative interlayer spacing are studied via molecular dynamics (MD) simulations. Compared with planar graphite surfaces, edges can achieve a 2-fold increase in the specific capacitance at a wider interlayer spacing of ∼5.0 Å. Unusually, the molecular origins for achieved charge storage are predominantly attributed to the structural evolutions of solvents occurring in the double layer, going beyond the traditional views of regulating the capacitance by ion adsorption/separation. Specifically, diverse ionic distributions exhibit similar screening ability and EDLC thickness, while water molecules can counterbalance the interfacial electric fields more effectively at edge site. The as-obtained findings will be instructive in designing graphene-based EDLCs for advanced capacitive performances.

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.

A comparative study of room temperature ionic liquids and their organic solvent mixtures near charged electrodes

The structural properties of electrolytes consisting of solutions of ionic liquids in a polar solvent at charged electrode surfaces are investigated using classical atomistic simulations. The studied electrolytes consisted of tetraethylammonium tetrafluoroborate (NEt4-BF4), 1-ethyl-3-methylimidazolium tetrafluoroborate (c2mim-BF4) and 1-octyl-3-methylimidazolium tetrafluoroborate (c8mim-BF4) salts dissolved in acetonitrile solvent. We discuss the influence of electrolyte concentration, chemical structure of the ionic salt, temperature, conducting versus semiconducting nature of the electrode, electrode geometry and surface roughness on the electric double layer structure and capacitance and compare these properties with those obtained for pure room temperature ionic liquids. We show that electrolytes consisting of solutions of ions can behave quite differently from pure ionic liquid electrolytes.

Concentration Fluctuations and Capacitive Response in Dense Ionic Solutions

We use molecular dynamics simulations in a constant potential ensemble to study the effects of solution composition on the electrochemical response of a double layer capacitor. We find that the capacitance first increases with ion concentration following its expected ideal solution behavior but decreases upon approaching a pure ionic liquid in agreement with recent experimental observations. The nonmonotonic behavior of the capacitance as a function of ion concentration results from the competition between the independent motion of solvated ions in the dilute regime and solvation fluctuations in the concentrated regime. Mirroring the capacitance, we find that the characteristic decay length of charge density correlations away from the electrode is also nonmonotonic. The correlation length first decreases with ion concentration as a result of better electrostatic screening but increases with ion concentration as a result of enhanced steric interactions. When charge fluctuations induced by correlated ion-solvent fluctuations are large relative to those induced by the pure ionic liquid, such capacitive behavior is expected to be generic.

Dense ionic fluids confined in planar capacitors: in- and out-of-plane structure from classical density functional theory

The ongoing scientific interest in the properties and structure of electric double layers (EDLs) stems from their pivotal role in (super)capacitive energy storage, energy harvesting, and water treatment technologies. Classical density functional theory (DFT) is a promising framework for the study of the in- and out-of-plane structural properties of double layers. Supported by molecular dynamics simulations, we demonstrate the adequate performance of DFT for analyzing charge layering in the EDL perpendicular to the electrodes. We discuss charge storage and capacitance of the EDL and the impact of screening due to dielectric solvents. We further calculate, for the first time, the in-plane structure of the EDL within the framework of DFT. While our out-of-plane results already hint at structural in-plane transitions inside the EDL, which have been observed recently in simulations and experiments, our DFT approach performs poorly in predicting in-plane structure in comparison to simulations. However, our findings isolate fundamental issues in the theoretical description of the EDL within the primitive model and point towards limitations in the performance of DFT in describing the out-of-plane structure of the EDL at high concentrations and potentials.

Understanding Capacitance Variation in Sub-nanometer Pores by in Situ Tuning of Interlayer Constrictions

The contribution of subnanometer pores in carbon electrodes to the charge-storage mechanism in supercapacitors has been the subject of intense debate for over a decade. Here, we provide a model system based on graphene oxide, which employs interlayer constrictions as a model for pore sizes that can be both controllably tuned and studied in situ during supercapacitor device use. Correlating electrochemical performance and in situ tuning of interlayer constrictions, we observe a peak in specific capacitance when interlayer constriction size reaches the diameters of unsolvated ions, supporting the hypothesized link between loss of ion solvation shell and anomalous capacitance increase for subnanometer pores.

Molecular dynamics simulation of amino acid ionic liquids near a graphene electrode: effects of alkyl side-chain length

The nanoscopic structure of amino acid ionic liquids (AAILs) as biodegradable electrolytes near a neutral graphene surface was studied by molecular dynamics (MD) simulation.

Tailoring graphene-based electrodes from semiconducting to metallic to increase the energy density in supercapacitors

The semiconducting character of graphene and some carbon-based electrodes can lead to noticeably lower total capacitances and stored energy densities in electric double layer (EDL)capacitors. This paper discusses the chemical and electronic structure modifications that enhance the available energy bands, density of states and quantum capacitance of graphene substrates near the Fermi level, therefore restoring the conducting character of these materials. The doping of graphene with p or n dopants, such as boron and nitrogen atoms, or the introduction of vacancy defects that introduce zigzag edges, can significantly increase the quantum capacitance within the potential range of interest for the energy storage applications by either shifting the Dirac point away from the Fermi level or by eliminating the Dirac point. We show that a combination of doping and vacancies at realistic concentrations is sufficient to increase the capacitance of a graphene-based electrode to within 1 μF cm(−2) from that of a metallic surface.Using a combination of ab initio calculations and classical molecular dynamics simulations we estimate how the changes in the quantum capacitance of these electrode materials affect the total capacitance stored by the open structure EDL capacitors containing room temperature ionic liquid electrolytes.