Differential capacitance of ionic liquid interface with graphite: the story of two double layers

Direct Measurement of the Differential Capacitance of Deep Eutectic Solvents on Platinum and Glassy Carbon Electrodes

Differential capacitance is a crucial parameter that connects the experimental observation of electrical double-layer behavior with theoretical models. However, the current number of reported differential capacitance values for deep eutectic solvents remains limited, making it challenging to verify or refute existing models. In this study, we systematically investigate the differential capacitance in deep eutectic solvents using chronoamperometry. By comparing metal and glassy carbon electrodes across various liquid combinations and ion concentrations, we observed a range of distinct capacitance characteristics. While some findings align with the existing mean-field model for ionic liquids, others clearly reflect the influence of electrode materials, with certain cases resisting full explanation by current theoretical models. These results underscore the importance of selecting appropriate electrode materials in experimental studies of such electrolytes and highlight the need for further theoretical advancements in understanding this complex liquid system.

Accounting for the Quantum Capacitance of Graphite in Constant Potential Molecular Dynamics Simulations

Molecular dynamics (MD) simulations at a constant electric potential are an essential tool to study electrochemical processes, providing microscopic information on the structural, thermodynamic, and dynamical properties. Despite the numerous advances in the simulation of electrodes, they fail to accurately represent the electronic structure of materials such as graphite. In this work, a simple parameterization method that allows to tune the metallicity of the electrode based on a quantum chemistry calculation of the density of states (DOS) is introduced. As a first illustration, the interface between graphite electrodes and two different liquid electrolytes, an aqueous solution of NaCl and a pure ionic liquid, at different applied potentials are studied. It is shown that the simulations reproduce qualitatively the experimentally-measured capacitance; in particular, they yield a minimum of capacitance at the point of zero charge (PZC), which is due to the quantum capacitance (QC) contribution. An analysis of the structure of the adsorbed liquids allows to understand why the ionic liquid displays a lower capacitance despite its large ionic concentration. In addition to its relevance for the important class of carbonaceous electrodes, this method can be applied to any electrode materials (e.g. 2D materials, conducting polymers, etc), thus enabling molecular simulation studies of complex electrochemical devices in the future.

Generalized Helmholtz model describes capacitance profiles of ionic liquids and concentrated aqueous electrolytes

In this work, we propose and validate a generalization of the Helmholtz model that can account for both "bell-shaped" and "camel-shaped" differential capacitance profiles of concentrated electrolytes, the latter being characteristic of ionic liquids. The generalization is based on introducing voltage dependence of both the dielectric constant "ϵr(V)" and thickness "L(V)" of the inner Helmholtz layer, as validated by molecular dynamics (MD) simulations. We utilize MD simulations to study the capacitance profiles of three different electrochemical interfaces: (1) graphite/[BMIm+][BF4-] ionic liquid interface; (2) Au(100)/[BMIm+][BF4-] ionic liquid interface; (3) Au(100)/1M [Na+][Cl-] aqueous interface. We compute the voltage dependence of ϵr(V) and L(V) and demonstrate that the generalized Helmholtz model qualitatively describes both camel-shaped and bell-shaped differential capacitance profiles of ionic liquids and concentrated aqueous electrolytes (in lieu of specific ion adsorption). In particular, the camel-shaped capacitance profile that is characteristic of ionic liquid electrolytes arises simply from combination of the voltage-dependent trends of ϵr(V) and L(V). Furthermore, explicit analysis of the inner layer charge density for both concentrated aqueous and ionic liquid double layers reveal similarities, with these charge distributions typically exhibiting a dipolar region closest to the electrode followed by a monopolar peak at larger distances. It is appealing that a generalized Helmholtz model can provide a unified description of the inner layer structure and capacitance profile for seemingly disparate aqueous and ionic liquid electrolytes.

Impedance Analysis of Capacitive and Faradaic Processes in the Pt/[Dema][TfO] Interface

The electrochemical reaction kinetics, especially the oxygen reduction reaction (ORR) at the cathode, is crucial for the performance of a fuel cell. In this study, the electrochemical processes on a polycrystalline Pt electrode in the presence of protic ionic liquid (PIL) electrolyte diethylmethylammonium triflate [Dema][TfO] are investigated by means of cyclic voltammetry and electrochemical impedance spectroscopy. Since water is continually produced during fuel cell operation, the effect of the water content in the PIL has been intensively analyzed. In order to reveal the dependence of the interfacial reaction characteristics on the electrode potential, the impedance spectra were simulated by an equivalent circuit whose parameters can be related to both Faradaic and capacitive processes. Two interfacial resistances were identified, which differ by about 3 orders of magnitude. The larger one is a charge transfer resistance that can be associated with slow Faradaic processes like the ORR and platinum oxidation/oxide reduction. The smaller resistance is probably linked with fast processes that involve water molecules, such as hydrogen deposition and oxidation. The high- and midfrequency capacitive processes are attributed to "classical" double layer and pseudocapacitive behavior, similar to those identified under nitrogen atmosphere.

Theoretical Insight into the Imidazolium-Based Ionic Liquid Interface Structure and Differential Capacitance on Au(111): Effects of the Cationic Substituent Group

Electric double layers (EDLs) play a key role in the electrochemical and energy storage of supercapacitors. It is important to understand the structure and properties of EDLs. In this work, quantum chemical calculations and molecular dynamics (MD) simulations are used to study the microstructure of EDLs of four different substituents of imidazolium-based bis(trifluoromethylsulfonyl)imide ionic liquids (ILs) on the Au(111) surface. It is shown that the particle interactions influence the different arrangements of the anion and cation. More alkyl substitutions and longer alkyl chains result in a higher ELUMO and thus a stronger interaction energy between cations and electrodes. Strong interactions produce linear patterns of anions/cations on the electrode and a maximum value of differential capacitance near PZC, whereas weak interactions generate worm-like patterns of anions/cations on Au(111) and a minimum value of differential capacitance near the PZC. We hold the opinion that the alkyl substitution has more effects on the EDLs. Our analysis provides a new perspective on EDLs structures at the atomic and molecular level. This study provides a good basis and guidance for further understanding the interface phenomena and characteristics of ionic liquids in electrochemical and energy device applications.

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.

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.

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.

Electronic screening using a virtual Thomas-Fermi fluid for predicting wetting and phase transitions of ionic liquids at metal surfaces

Of relevance to energy storage, electrochemistry and catalysis, ionic and dipolar liquids display unexpected behaviours-especially in confinement. Beyond adsorption, over-screening and crowding effects, experiments have highlighted novel phenomena, such as unconventional screening and the impact of the electronic nature-metallic versus insulating-of the confining surface. Such behaviours, which challenge existing frameworks, highlight the need for tools to fully embrace the properties of confined liquids. Here we introduce a novel approach that involves electronic screening while capturing molecular aspects of interfacial fluids. Although available strategies consider perfect metal or insulator surfaces, we build on the Thomas-Fermi formalism to develop an effective approach that deals with any imperfect metal between these asymptotes. Our approach describes electrostatic interactions within the metal through a 'virtual' Thomas-Fermi fluid of charged particles, whose Debye length sets the screening length λ. We show that this method captures the electrostatic interaction decay and electrochemical behaviour on varying λ. By applying this strategy to an ionic liquid, we unveil a wetting transition on switching from insulating to metallic conditions.

Interfacial water asymmetry at ideal electrochemical interfaces

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

Microscopic origin of the effect of substrate metallicity on interfacial free energies

We investigate the effect of the metallic character of solid substrates on solid-liquid interfacial thermodynamics using molecular simulations. Building on the recent development of a semiclassical Thomas-Fermi model to tune the metallicity in classical molecular dynamics simulations, we introduce a thermodynamic integration framework to compute the evolution of the interfacial free energy as a function of the Thomas-Fermi screening length. We validate this approach against analytical results for empty capacitors and by comparing the predictions in the presence of an electrolyte with values determined from the contact angle of droplets on the surface. The general expression derived in this work highlights the role of the charge distribution within the metal. We further propose a simple model to interpret the evolution of the interfacial free energy with voltage and Thomas-Fermi length, which allows us to identify the charge correlations within the metal as the microscopic origin of the evolution of the interfacial free energy with the metallic character of the substrate. This methodology opens the door to the molecular-scale study of the effect of the metallic character of the substrate on confinement-induced transitions in ionic systems, as reported in recent atomic force microscopy and surface force apparatus experiments.

Electrochemical Characterization of Single Layer Graphene/Electrolyte Interface: Effect of Solvent on the Interfacial Capacitance

The development of the basic understanding of the charge storage mechanisms in electrodes for energy storage applications needs deep characterization of the electrode/electrolyte interface. In this work, we studied the charge of the double layer capacitance at single layer graphene (SLG) electrode used as a model material, in neat (EMIm-TFSI) and solvated (with acetonitrile) ionic liquid electrodes. The combination of electrochemical impedance spectroscopy and gravimetric electrochemical quartz crystal microbalance (EQCM) measurements evidence that the presence of solvent drastically increases the charge carrier density at the SLG/ionic liquid interface. The capacitance is thus governed not only by the electronic properties of the graphene, but also by the specific organization of the electrolyte side at the SLG surface originating from the strong interactions existing between the EMIm+ cations and SLG surface. EQCM measurements also show that the carbon structure, with the presence of sp2 carbons, affects the charge storage mechanism by favoring counter-ion adsorption on SLG electrode versus ion exchange mechanism in amorphous porous carbons.

Effect of Surface Oxygen on the Wettability and Electrochemical Properties of Boron-Doped Nanocrystalline Diamond Electrodes in Room-Temperature Ionic Liquids

This paper reports on how the surface chemistry of boron-doped nanocrystalline diamond (BDD) thin-film electrodes (H vs O) affects the wettability and electrochemical properties in two room-temperature ionic liquids (RTILs): [BMIM][PF₆] and [HMIM][PF₆]. Comparative measurements were made in 0.5 mol L^-1 H₂SO₄. The BDD electrodes were modified by microwave or radio-frequency (RF) plasma treatment in H₂ (H-BDD), Ar (Ar-BDD), or O₂ (O-BDD). These modifications produced low-, medium-, and high-oxygen surface coverages. Atomic O/C ratios, as determined by X-ray photoelectron spectroscopy (XPS), were 0.01 for H-BDD, 0.08 for Ar-BDD, and 0.17 for O-BDD. The static contact angle of ultrapure water on the modified electrodes decreased from 110° (H-BDD) to 41° (O-BDD) with increasing surface oxygen coverage, as expected as the surface becomes more hydrophilic. Interestingly, the opposite trend was seen for both RTILs as the contact angle increased from 20° (H-BDD) to 50° (O-BDD) with increasing surface oxygen coverage. The cyclic voltammetric background current and potential-dependent capacitance in both RTILs were largest for BDD electrodes with the lowest O/C ratio (H-BDD) and smallest contact angle. Slightly larger voltammetric background currents and capacitance were observed in [HMIM][PF₆] than in [BMIM][PF₆]. Capacitance values ranged from 8 to 16 μF cm^-2 over the potential range for H-BDD and from 4 to 6 μF cm^-2 for O-BDD. The opposite trend was observed in H₂SO₄ as the voltammetric background current and capacitance were largest for BDD electrodes with the highest O/C ratio (O-BDD) and smallest contact angle. In summary, reducing the surface oxygen on BDD electrodes increases the wettability to two RTILs and this increases the voltammetric background current and capacitance.

Superionic liquids in conducting nanoslits: A variety of phase transitions and ensuing charging behavior

We develop a theory of charge storage in ultranarrow slitlike pores of nanostructured electrodes. Our analysis is based on the Blume-Capel model in an external field, which we solve analytically on a Bethe lattice. The obtained solutions allow us to explore the complete phase diagram of confined ionic liquids in terms of the key parameters characterizing the system, such as pore ionophilicity, interionic interaction energy, and voltage. The phase diagram includes the lines of first- and second-order, direct and re-entrant phase transitions, which are manifested by singularities in the corresponding capacitance-voltage plots. Testing our predictions experimentally requires monodisperse, conducting ultranarrow slit pores, to permit only one layer of ions, and thick pore walls, to prevent interionic interactions across the pore walls. However, some qualitative features, which distinguish the behavior of ionophilic and ionophobic pores and their underlying physics, may emerge in future experimental studies of more complex electrode structures.

Effects of Solvent Concentration on the Performance of Ionic-Liquid/Carbon Supercapacitors

We use molecular dynamics simulations to investigate the effects of solvent concentration on the bulk properties of an ion liquid electrolyte and the electrochemical performance on carbon-based electrodes, including pristine graphene, oxidized graphene, graphene armchair edge, graphene zigzag edge, onion-like carbon, and slit-pore carbon. We find that diluting the electrolyte reduces the number of ion pairs in the bulk and improves ion dynamics. The capacitance of the two-edge electrodes decreases monotonically as the solvent concentration increases, while the capacitance of other nonedge electrodes exhibits nonmonotonic behavior and a capacitance maximum is observed. Further analyses on the electric double layer reveals two competing factors: solvation reduces the charge overscreening effect, but it also causes the dilution of absorbed ion concentration. While the former increases the capacitance in the low dilution regime, the latter decreases the capacitance in the high dilution regime. In addition, the dilution also significantly improves the ion dynamics at the interface. Our simulation results demonstrate that diluting an ionic liquid electrolyte could potentially boost the power density while maintaining or even slightly increasing the energy density with a careful selection of solvent concentrations on a nonedge carbon electrode.

Carbon-Nanotube-Electrolyte Interface: Quantum and Electric Double Layer Capacitance

We present a comprehensive study of the electrochemical capacitance between a one-dimensional electronic material and an electrolyte. In contrast to a conventional, planar electrode, the nanoscale dimension of the electrode (with diameter smaller than the Debye length and approaching the size of the ions in solution) qualitatively changes the capacitance, which we measure and model herein. Furthermore, the finite density of states in these low dimensional electronic systems results in a quantum capacitance, which is comparable to the electrochemical capacitance. Using electrochemical impedance spectroscopy (EIS), we measure the ensemble average, complex, frequency dependent impedance (from 0.1 Hz to 1 MHz) between a purified (99.9%) semiconducting nanotube network and an aqueous electrolyte (KCl) at different concentrations between 10 mM and 1 M. The potential dependence of the capacitance is convoluted with the potential dependence of the in-plane conductance of the nanotube network, which we model using a transmission-line model to account for the frequency dependent in-plane impedance as well as the total interfacial impedance between the network and the electrolyte. The ionic strength dependence of the capacitance is expected to have a root cause from the double layer capacitance, which we model using a modified Poisson-Boltzmann equation. The relative contributions from those two capacitances can be quantitatively decoupled. We find a total capacitance per tube of 0.67-1.13 fF/μm according to liquid gate potential varying from -0.5 to -0.7 V.

Direct Measurement of the Differential Capacitance of Solvent-Free and Dilute Ionic Liquids

Differential capacitance is a key quantity in the understanding of electrical double-layer charging of electrolytes. However, experimental observations of ionic liquid systems are controversial, inconsistent, and often unable of confirming or refuting existing theories as well as highlighting discrepancies between the measurement techniques. We study the differential capacitance in both pure and dilute ionic liquids at room temperature. Using chronoamperometry to measure the differential capacitance of the liquids at a polycrystalline platinum electrode, we find good agreement between the measured capacitance curves and the extended mean-field model of Goodwin-Kornyshev [Goodwin, Z. A.; et al. Electrochim. Acta. 2017, 225, 190-197]. A crossover is found from the pure to the dilute regime, as shown by a transition from a camel-shape capacitance curve to a U-like one, together with a nonmonotonic dependence of capacitance with electrolyte concentration.

The influence of water content in a proton-conducting ionic liquid on the double layer properties of the Pt/PIL interface

The influence of the water content of 2-sulfoethylmethylammonium trifluoromethanesulfonate [2-Sema][TfO] on the double layer properties of the interface of platinum and the proton conducting ionic liquid (PIL) is investigated by means of impedance spectroscopy and cyclic voltammetry. By fitting the impedance spectra as complex capacitances, up to four differential double layer capacitances and corresponding time constants are obtained, depending on the potential (U = 0-1.6 V/RHE), water content (0.7-6.1 wt%) and temperature (T = 70-110 °C). Within the whole potential range investigated, a high frequency capacitance, C₁, and a low frequency capacitance, C₂, can be calculated. In the potential region of hydrogen underpotential deposition (H_UPD), C₁ can be separated into two parts, C_1a and C_1b. Whereas the high frequency capacitive processes can mainly be attributed to ion transport processes in the double layer, the low frequency process is ascribed to changes in the interfacial layer, including ad-/desorption and Faradaic processes. Alternative interpretations regarding the reorientation of ions, reconstruction of the metal surface and partial electron transfer between anions and Pt are considered.

Capacitance of Nanoporous Carbon-Based Supercapacitors Is a Trade-Off between the Concentration and the Separability of the Ions

Nanoporous carbon-based supercapacitors store electricity through adsorption of ions from the electrolyte at the surface of the electrodes. Room temperature ionic liquids, which show the largest ion concentrations among organic liquid electrolytes, should in principle yield larger capacitances. Here, we show by using electrochemical measurements that the capacitance is not significantly affected when switching from a pure ionic liquid to a conventional organic electrolyte using the same ionic species. By performing additional molecular dynamics simulations, we interpret this result as an increasing difficulty of separating ions of opposite charges when they are more concentrated, that is, in the absence of a solvent that screens the Coulombic interactions. The charging mechanism consistently changes with ion concentration, switching from counterion adsorption in the diluted organic electrolyte to ion exchange in the pure ionic liquid. Contrarily to the capacitance, in-pore diffusion coefficients largely depend on the composition, with a noticeable slowing of the dynamics in the pure ionic liquid.

Characterization of the Electric Double Layer Formation Dynamics of a Metal/Ionic Liquid/Metal Structure

Although ionic liquids (ILs) have been used extensively in recent years as a high-capacitance "dielectric" in electric double layer transistors, the dynamics of the double layer formation have remained relatively unexplored. Better understanding of the dynamics and relaxation processes involved in electric double layer formation will guide device optimization, particularly with regard to switching speed. In this paper, we explore the dynamical characteristics of an IL in a metal/ionic liquid/metal (M/IL/M) capacitor. In particular, we examine a Au/IL/Au structure where the IL is 1-butyl-1-methylpyrrolidinium tris(pentafluoroethyl)trifluorophosphate. The experiments consist of frequency-dependent impedance measurements and time-dependent current vs voltage measurements for applied linear voltage ramps and abrupt voltage steps. The parameters of an equivalent circuit model are determined by fits to the impedance vs frequency data and subsequently verified by calculating the current vs voltage characteristics for the applied potential profiles. The data analysis indicates that the dynamics of the structure are characterized by a wide distribution of relaxation times spanning the range of less than microseconds to longer than seconds. Possible causes for these time scales are discussed.

Pressing a spring: what does it take to maximize the energy storage in nanoporous supercapacitors?

We discuss the nonlinear effects and efficiency of charge storage in supercapacitors with nanoporous electrodes and ionic liquids, and demonstrate that to maximize the stored energy, it may be beneficial to create 'obstacles' or 'difficulties' in charging. This can be achieved by making thermodynamically unfavourable conditions for ions inside nanopores, or more favourable outside. We show by means of Monte Carlo simulations that such 'ionophobic' pores store energy more efficiently and can provide equivalent or even better energy capacity. Since the recent analysis predicts much faster charging of ionophobic nanopores, we conclude that such pores offer a better option for simultaneous energy/power optimization.

The interface between HOPG and 1-butyl-3-methyl-imidazolium hexafluorophosphate

The interface between highly oriented pyrolytic graphite (HOPG) and 1-butyl-3-metyl-imidazolium hexafluorophosphate (BMIPF₆) has been studied using cyclic voltammetry, electrochemical impedance spectroscopy, immersion charge measurements and in situ scanning tunneling microscopy (in situ STM).

Anomalous Capacitance Maximum of the Glassy Carbon-Ionic Liquid Interface through Dilution with Organic Solvents

We use electrochemical impedance spectroscopy to measure the effect of diluting a hydrophobic room temperature ionic liquid with miscible organic solvents on the differential capacitance of the glassy carbon-electrolyte interface. We show that the minimum differential capacitance increases with dilution and reaches a maximum value at ionic liquid contents near 5-10 mol% (i.e., ∼1 M). We provide evidence that mixtures with 1,2-dichloroethane, a low-dielectric constant solvent, yield the largest gains in capacitance near the open circuit potential when compared against two traditional solvents, acetonitrile and propylene carbonate. To provide a fundamental basis for these observations, we use a coarse-grained model to relate structural variations at the double layer to the occurrence of the maximum. Our results reveal the potential for the enhancement of double-layer capacitance through dilution.

Interfacial ionic 'liquids': connecting static and dynamic structures

It is well known that room temperature ionic liquids (RTILs) often adopt a charge-separated layered structure, i.e. with alternating cation- and anion-rich layers, at electrified interfaces. However, the dynamic response of the layered structure to temporal variations in applied potential is not well understood. We used in situ, real-time x-ray reflectivity to study the potential-dependent electric double layer (EDL) structure of an imidazolium-based RTIL on charged epitaxial graphene during potential cycling as a function of temperature. The results suggest that the graphene-RTIL interfacial structure is bistable in which the EDL structure at any intermediate potential can be described by the combination of two extreme-potential structures whose proportions vary depending on the polarity and magnitude of the applied potential. This picture is supported by the EDL structures obtained by fully atomistic molecular dynamics simulations at various static potentials. The potential-driven transition between the two structures is characterized by an increasing width but with an approximately fixed hysteresis magnitude as a function of temperature. The results are consistent with the coexistence of distinct anion- and cation-adsorbed structures separated by an energy barrier (∼0.15 eV).

Effects of repulsive interaction on the electric double layer of an imidazolium-based ionic liquid by molecular dynamics simulation

Effects of repulsion between BMIM⁺/PF₆⁻ and the graphite electrode on electric double layer was studied by molecular dynamics simulation.