Effects of Alkyl Chain Length on Interfacial Structure and Differential Capacitance in Graphene Supercapacitors: A Molecular Dynamics Simulation Study

Structural Properties of [N1888][TFSI] Ionic Liquid: A Small Angle Neutron Scattering and Polarizable Molecular Dynamics Study

In this study, we investigate the quaternary ammonium-based ionic liquid (QAIL), methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N1888][TFSI], utilizing small angle neutron scattering (SANS) measurements and polarizable molecular dynamics (MD) simulations to characterize the short- and long-range liquid structure. Scattering structure factors show signatures of three length scales in reciprocal space indicative of alternating polarity (k ∼ 0.44 Å-1), charge (k ∼ 0.75 Å-1), and neighboring or adjacent (k ∼ 1.46 Å-1) domains. Excellent agreement between simulation and experimental scattering structure factors validates various simulation analyses that provide detailed atomistic characterization of the different length scale correlations. The first solvation shell structure is illustrated by obtaining radial, angular, dihedral, and combined distribution functions, where two dominant spatial motifs, N+···N- and N+···O-, compete for optimal packing around the polar head of the [N1888]+ cation. Intermediate and long-range structures are governed by the balance between local electroneutrality and octyl chain networking, respectively. By computing the charge-correlation structure factor, SZZ, and the spatial extent of the octyl chain network using graph theory, the bulk-phase structure of [N1888][TFSI] is characterized in terms of electrostatic screening and apolar domain formation length scales.

Nanostructure and Dynamics of the Locally Concentrated Ionic Liquid 2:1 (wt:wt) HMIM FAP:TFTFE and HMIM FAP on Graphite and Gold Electrodes as a Function of Potential

Video-rate atomic force microscopy (AFM) is used to record the near-surface nanostructure and dynamics of one pure ionic liquid (IL), 1-hexyl-3-methylimidazolium tris(pentafluoroethyl)trifluorophosphate (HMIM FAP), and a locally-concentrated IL comprising HMIM FAP with the low viscosity diluent 1,1,2,2-tetrafluoroethyl 2,2,2-trifluoroethyl ether (TFTFE), on highly oriented pyrolytic graphite (HOPG) and Au(111) electrodes as a function of potential. Over the potential range measured (open-circuit potential ± 1 V), different near-surface nanostructures are observed. For pure HMIM FAP, globular aggregates align in rows on HOPG, whereas elongated and worm-like nanostructures form on Au(111). For 2:1 (wt:wt) HMIM FAP:TFTFE, larger and less defined diluent swollen IL aggregates are present on both electrodes. Long-lived near-surface nanostructures for HMIM FAP and the 2:1 (wt:wt) HMIM FAP:TFTFE persist on both electrodes. 2:1 (wt:wt) HMIM FAP:TFTFE mixture diffuses more rapidly than pure HMIM FAP on both electrodes with obviously higher diffusion coefficients on HOPG than on Au(111) due to weaker electrostatic and solvophobic interactions between near-surface aggregates and Stern layer ions. These outcomes provide valuable insights for a wide range of IL applications in interface sciences, including electrolytes, catalysts, lubricants, and sensors.

Effect of Surface Chemistry on the Electrical Double Layer in a Long-Chain Ionic Liquid

Room temperature ionic liquids (ILs) can create a strong accumulation of charges at solid interfaces by forming a very thin and dense electrical double layer (EDL). The structure of this EDL has important consequences in numerous applications involving ILs, for example, in supercapacitors, sensors, and lubricants, by impacting the interfacial capacitance, the charge carrier density of semiconductors, as well as the frictional properties of the interfaces. We have studied the interfacial structure of a long chain imidazolium-based IL (1-octyl-3-methylimidazolium dicyanamide) on several substrates: mica, silica, silicon, and molybdenum disulfide (MoS2), using atomic force microscopy (AFM) experiments and molecular dynamics (MD) simulations. We have observed 3 types of interfacial structures for the same IL, depending on the chemistry of the substrate and the water content, showing that the EDL structure is not an intrinsic property of the IL. We evidenced that at a low water content, neutral and apolar (thus hydrophobic) substrates promote a thin layer structure, where the ions are oriented parallel to the substrate and cations and anions are mixed in each layer. In contrast, a strongly charged (thus hydrophilic) substrate yields an extended structuration into several bilayers, while a heterogeneous layering with loose bilayer regions was observed on an intermediate polar and weakly charged substrate and on an apolar one at a high bulk water content. In the latter case, water contamination favors the formation of bilayer patches by promoting the segregation of the long chain IL into polar and apolar domains.

Molecular dynamics study of fluorosulfonyl ionic liquids as electrolyte for electrical double layer capacitors

The development of high-performance supercapacitors is an important goal in the field of energy storage. Ionic liquids (ILs) are promising electrolyte materials for efficient energy storage in supercapacitors, because of the high stability, low volatility, and wider electrochemical stability window than traditional electrolytes. However, ILs-based supercapacitors usually show a relatively lower power density owing to the inherent viscosity-induced low electrical conductivity. Fluorosulfonyl ILs have aroused much attention in energy storage devices due to its low toxicity and excellent stability. Here, we propose that structural modification is an effective way to improve the energy storage performance of fluorosulfonyl ILs through the classical molecular dynamics (MD) method. Four fluorosulfonyl ILs with different sizes and symmetries were considered. Series of properties including conductivity, interface structure, and double-layer capacitance curves were systematically investigated. The results show that smaller size and more asymmetric structure can enhance self-diffusion coefficient and conductivity, and improve the electrochemical performance. Appropriate modification of the electrodes can further enhance the capacitive performance. Our work provides an opportunity to further understand and develop the fluorosulfonyl ILs electrolyte in supercapacitors.

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.

Theory-augmented informatics of ionic liquid electrolytes for co-design with nanoporous electrode materials

Ionic liquids possess compelling properties and vast chemical diversity, promising unprecedented performance and tunability for advanced electrochemical applications in catalysis, sensing, and energy storage. However, with broad tunability comes intractable, multidimensional parameter spaces not easily traversed by empirical approaches, limiting both scientific understanding and technological breakthroughs with these novel materials. In this Communication, we propose an extensible figure of merit that co-optimizes key ionic liquid properties, including electrochemical stability window, viscosity, and molecular ion size with respect to pore sizes of nanoporous electrodes typically utilized in electrochemical technologies. We coupled density functional theory (DFT) with informatics to augment physiochemical property databases to screen for high-performance room-temperature ionic liquid (RTIL) candidate compounds. This co-design framework revealed a number of promising RTILs that are underrepresented in the literature and thus warrant future follow-up investigations.

Molecular dynamics investigation of charging process in polyelectrolyte-based supercapacitors

Supercapacitors are one of the technologically impressive types of energy storage devices that are supposed to fill the gap between chemical batteries and dielectric capacitors in terms of power and energy density. Many kinds of materials have been investigated to be used as supercapacitors' electrolytes to overcome the known limitations of them. The properties of polymer-based electrolytes show a promising way to defeat some of these limitations. In this paper, a simplified model of polymer-based electrolytes between two electrodes is numerically investigated using the Molecular Dynamics simulation. The simulations are conducted for three different Bjerrum lengths and a typical range of applied voltages. The results showed a higher differential capacitance compared to the cases using ionic-liquid electrolytes. Our investigations indicate a rich domain in molecular behaviors of polymer-based electrolytes that should be considered in future supercapacitors.

Imidazolium Triflate Ionic Liquids' Capacitance-Potential Relationships and Transport Properties Affected by Cation Chain Lengths

In this paper we report the effects of five imidazolium cations with varying alkyl chain lengths to study the effects of cation size on capacitance versus voltage behavior. The cations include ethyl-, butyl-, hexyl-, octyl-, and decyl-3-methylimidazolium, all paired with a triflate anion. We analyze the capacitance with respect to the cation alkyl chain length qualitatively and quantitatively by analyzing changes in the capacitance-potential curvature shape and magnitude across several standard scanning protocols and electrochemical techniques. Further, three transport properties (viscosity, diffusion coefficient, and electrical conductivity) are experimentally determined and integrated into the outcomes. Ultimately, we find higher viscosities, lower diffusion coefficients, and lower electrical conductivities when the alkyl chain length is increased. Also, capacitance values increase with cation size, except 1-octyl-3-methylimidazolium, which does not follow an otherwise linear trend. This capacitive increase is most pronounced when sweeping the potential in the cathodic direction. These findings challenge the conventional hypothesis that increasing the length of the alkyl chain of imidazolium cations diminishes the capacitance and ionic liquid performance in charge storage.

The Influence of Anion Structure on the Ionic Liquids/Au (100) Interface by Molecular Dynamics Simulations

The microstructure of electrical double layers (EDLs) of 1-butyl-3-methylimidazolium tetrafluoroborate ([Bmim]BF₄), 1-butyl-3-methylimidazoliumhexafluorophosphate ([Bmim]PF₆), and 1-butyl-3-methylimidazoliumbis (trifluoromethylsulfonyl) imide ([Bmim]TFSI) were studied by quantum chemical calculation and molecular dynamics simulation. For the set of ionic liquids investigated here, we found some interesting universal laws due to differences geometry and electronic structure of anions. We show that the morphology of the three anions on the electrode surface is different due to the different geometric structure. The plane formed by the bottom three atoms of the symmetrically tetrahedral BF₄^- and the bottom atom of the symmetrically octahedral PF₆^- face the electrode whether the electrode is charged or not, while the conformation of twisted V-shaped TFSI^- changes with different surface charges on the electrode. Meanwhile, we also demonstrate that the energy of highest occupied molecular orbital (E_HOMO), the energy of lowest unoccupied molecular orbital (E_LUMO) and their energies gap (ΔE) are very interesting due to different electronic structure of anions. Specially, the E_HOMO, E_LUMO, and ΔE were related to the electronegativity of the central atom in the case of the same symmetry on the neutral surface. The more electronegative the central atom is, the lower E_HOMO, E_LUMO and higher ΔE values are. However, on the charged surface, the interaction between anion and electrode is opposite to ΔE. Moreover, different arrangements of anion and cation are related to the interaction between particles. The stronger interaction leads a double-row structure and the weak interaction lead worm-like and island patterns on Au (100) surface. In general, we observed that the higher ΔE cause stronger interaction, which lead to different patterns on Au (100) surface. Meanwhile, we also confirmed that the stronger interaction between particles and electrode lead to the thinner effective EDL and a large differential capacitance value. These results provide a new perspective for double-layer structure in atomic and molecular level. This is helpful to deepen the understanding of the interface phenomena and characteristics of [Bmim]BF₄, [Bmim]PF₆, and [Bmim]TFSI on Au (100) system and provide theoretical basis for the application of these kind of systems.

Effects of Layering and Supporting Substrate on Liquid Slip at the Single-Layer Graphene Interface

Understanding modulation of liquid molecule slippage along graphene surfaces is crucial for many promising applications of two-dimensional materials, such as in sensors, nanofluidic devices, and biological systems. Here, we use force measurements by atomic force microscopy (AFM) to directly measure hydrodynamic, solvation, and frictional forces along the graphene plane in seven liquids. The results show that the greater slip lengths correlate with the interfacial ordering of the liquid molecules, which suggests that the ordering of the liquid forming multiple layers promotes slip. This phenomenon appears to be more relevant than solely the wetting behavior of graphene or the solid-liquid interaction energy, as traditionally assumed. Furthermore, the slip boundary condition of the liquids along the graphene plane is sensitive to the substrate underneath graphene, indicating that the underlying substrate affects graphene's interaction with the liquid molecules. Because interfacial slip can have prominent consequences on the pressure drop, on electrical and diffusive transport through nanochannels, and on lubrication, this work can inspire innovation in many applications through the modulation of the substrate underneath graphene and of the interfacial ordering of the liquid.

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

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

Investigation of the Ionic Liquid Graphene Electric Double Layer in Supercapacitors Using Constant Potential Simulations

In this work, we investigate the effect of the cation structure on the structure and dynamics of the electrode–electrolyte interface using molecular dynamics simulations. A constant potential method is used to capture the behaviour of 1-ethyl-3-methylimidazolium bis (trifluoromethane)sulfonimide ([C2mim][NTf2]) and butyltrimethylammonium bis(trifluoromethane) sulfonimide ([N4,1,1,1][NTf2]) ionic liquids at varying potential differences applied across the supercapacitor. We find that the details of the structure in the electric double layer and the dynamics differ significantly, yet the charge profile and capacitance do not vary greatly. For the systems considered, charging results in the rearrangement and reorientation of ions within ∼1 nm of the electrode rather than the diffusion of ions to/from the bulk region. This occurs on timescales of O(10 ns) for the ionic liquids considered, and depends on the viscosity of the fluid.

Interfacial structuring of non-halogenated imidazolium ionic liquids at charged surfaces: effect of alkyl chain length

Control of the interfacial structures of ionic liquids (ILs) at charged interfaces is important to many of their applications, including in energy storage solutions, sensors and advanced lubrication technologies utilising electric fields.

Many-chain effects on the co-nonsolvency of polymer brushes in a good solvent mixture

Polymer brushes normally swell in a good solvent and collapse in a poor solvent. An abnormal response of polymer brushes, so-called co-nonsolvency, is the phenomenon where the brush counter-intuitively collapses in a good solvent mixture. In this work, we employed molecular dynamics simulations to investigate the structural properties of the grafted polymers in the occurrence of co-nonsolvency. Brushes with various grafting densities were considered to study the effect of topologically excluded volumes on the co-nonsolvency. We found that the brush height follows a novel scaling behavior of the grafting density h ∼ σg0.71 in the co-nonsolvent mixture. Using the scaling exponent and Alexander-de Gennes theory, an analytic function that predicts the monomer density was obtained. The many-chain effects in the co-nonsolvent lead to the formation of both intermolecular and intramolecular bridging structures. Increasing the grafting density entails lower looping events occuring because of the intermolcular bridging, causing diverse structural properties. We report how the average thickness, the polymer orientation, and the looping probability vary as the grafting density increases. Based on these observations, we constructed a phase diagram of the polymer brush system using the average thickness and orientation as order parameters. Our simulations and analytical results reveal the nature of co-nonsolvency in polymer brushes in an explicit way and will help to provide practical guidelines for applications such as drug delivery and sensor devices.

Delfos: deep learning model for prediction of solvation free energies in generic organic solvents

Prediction of aqueous solubilities or hydration free energies is an extensively studied area in machine learning applications in chemistry since water is the sole solvent in the living system. However, for non-aqueous solutions, few machine learning studies have been undertaken so far despite the fact that the solvation mechanism plays an important role in various chemical reactions. Here, we introduce Delfos (deep learning model for solvation free energies in generic organic solvents), which is a novel, machine-learning-based QSPR method which predicts solvation free energies for various organic solute and solvent systems. A novelty of Delfos involves two separate solvent and solute encoder networks that can quantify structural features of given compounds via word embedding and recurrent layers, augmented with the attention mechanism which extracts important substructures from outputs of recurrent neural networks. As a result, the predictor network calculates the solvation free energy of a given solvent-solute pair using features from encoders. With the results obtained from extensive calculations using 2495 solute-solvent pairs, we demonstrate that Delfos not only has great potential in showing accuracy comparable to that of the state-of-the-art computational chemistry methods, but also offers information about which substructures play a dominant role in the solvation process.

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.

Graphene Meets Ionic Liquids: Fermi Level Engineering via Electrostatic Forces

Graphene-based two-dimensional (2D) materials are promising candidates for a number of different energy applications. A particularly interesting one is in next generation supercapacitors, where graphene is being explored as an electrode material in combination with room temperature ionic liquids (ILs) as electrolytes. Because the amount of energy that can be stored in such supercapacitors critically depends on the electrode-electrolyte interface, there is considerable interest in understanding the structure and properties of the graphene/IL interface. Here, we report the changes in the properties of graphene upon adsorption of a homologous series of alkyl imidazolium tetrafluoroborate ILs using a combination of experimental and theoretical tools. Raman spectroscopy reveals that these ILs cause n-type doping of graphene, and the magnitude of doping increases with increasing cation chain length despite the expected decrease in the density of surface-adsorbed ions. Molecular modeling simulations show that doping originates from the changes in the electrostatic potential at the graphene/IL interface. The findings described here represent an important step in developing a comprehensive understanding of the graphene/IL interface.

Understanding Electric Double-Layer Gating Based on Ionic Liquids: from Nanoscale to Macroscale

In electric double-layer transistors (EDLTs), it is well known that the EDL formed by ionic liquids (ILs) can induce an ultrahigh carrier density at the semiconductor surface, compared to solid dielectric. However, the mechanism of device performance is still not fully understood, especially at a molecular level. Here, we evaluate the gating performance of amorphous indium gallium zinc oxide (a-IGZO) transistor coupled with a series of imidazolium-based ILs, using an approach combining of molecular dynamics simulation and finite element modeling. Results reveal that the EDL with different ion structures could produce inhomogeneous electric fields at the solid-electrolyte interface, and the heterogeneity of electric field-induced charge distributions at semiconductor surface could reduce the electrical conductance of a-IGZO during gating process. Meanwhile, a resistance network analysis was adopted to bridge the nanoscopic data with the macroscopic transfer characteristics of IL-gated transistor, and showed that our theoretical results could well estimate the gating performance of practical devices. Thereby, our findings could provide both new concepts and modeling techniques for IL-gated transistors.

Design of Supercapacitor Electrodes Using Molecular Dynamics Simulations

Electric double-layer capacitors (EDLCs) are advanced electrochemical devices for energy storage and have attracted strong interest due to their outstanding properties. Rational optimization of electrode-electrolyte interactions is of vital importance to enhance device performance for practical applications. Molecular dynamics (MD) simulations could provide theoretical guidelines for the optimal design of electrodes and the improvement of capacitive performances, e.g., energy density and power density. Here we discuss recent MD simulation studies on energy storage performance of electrode materials containing porous to nanostructures. The energy storage properties are related to the electrode structures, including electrode geometry and electrode modifications. Altering electrode geometry, i.e., pore size and surface topography, can influence EDL capacitance. We critically examine different types of electrode modifications, such as altering the arrangement of carbon atoms, doping heteroatoms and defects, which can change the quantum capacitance. The enhancement of power density can be achieved by the intensified ion dynamics and shortened ion pathway. Rational control of the electrode morphology helps improve the ion dynamics by decreasing the ion diffusion pathway. Tuning the surface properties (e.g., the affinity between the electrode and the ions) can affect the ion-packing phenomena. Our critical analysis helps enhance the energy and power densities of EDLCs by modulating the corresponding electrode structures and surface properties.