Polarizable Force Field Development and Molecular Dynamics Simulations of Ionic Liquids

Relaxation in a Prototype Ionic Liquid: Influence of Water on the Dynamics

The influence of water on the relaxation of a prototype ionic liquid (IL) C₈mimBF₄ is examined in the IL-rich regime combining quasi-elastic neutron scattering (QENS) and molecular dynamics (MD) simulations. The QENS and MD simulations results for relaxation of IL and the equimolar mixture with water probed by the dynamics of the C₈mim hydrogen atoms in the time range of 2 ps to 1 ns are in excellent agreement. The QENS data show that translational relaxation increases by a factor of 7 on the addition of water, while rotational relaxation involving multiple processes fitted by a KWW function with low β values is speeded up by a factor of 3 on the time scale of QENS measurements. The MD simulations show that the cation diffusion coefficient, inverse viscosity, and ionic conductivity increase on the addition of water, consistent with the very small change in ionicity. The difficulties in obtaining rotational and translational diffusion coefficients from fits to QENS experiments of pure ILs and IL-water mixtures are discussed.

Quantum Chemical Methods for the Prediction of Energetic, Physical, and Spectroscopic Properties of Ionic Liquids

The accurate prediction of physicochemical properties of condensed systems is a longstanding goal of theoretical (quantum) chemistry. Ionic liquids comprising entirely of ions provide a unique challenge in this respect due to the diverse chemical nature of available ions and the complex interplay of intermolecular interactions among them, thus resulting in the wide variability of physicochemical properties, such as thermodynamic, transport, and spectroscopic properties. It is well understood that intermolecular forces are directly linked to physicochemical properties of condensed systems, and therefore, an understanding of this relationship would greatly aid in the design and synthesis of functionalized materials with tailored properties for an application at hand. This review aims to give an overview of how electronic structure properties obtained from quantum chemical methods such as interaction/binding energy and its fundamental components, dipole moment, polarizability, and orbital energies, can help shed light on the energetic, physical, and spectroscopic properties of semi-Coulomb systems such as ionic liquids. Particular emphasis is given to the prediction of their thermodynamic, transport, spectroscopic, and solubilizing properties.

Molecular Structure and Dynamics of Ionic Liquids in a Rigid-Rod Polyanion-Based Ion Gel

The recent fabrication of liquid crystalline ion gels featuring rigid-rod polyanions aligned within room-temperature ionic liquids (RTILs) opens up exciting new avenues for engineering ion conducting materials. These gels exhibit an unusual combination of properties including high ionic conductivity, distinct transport anisotropy, and widely tunable elastic modulus. Using molecular simulations, we study the structure and dynamics of the ions in an ion gel consisting of rigid-rod polyanions and [C₂mim][TfO] RTILs. We show that the ion distribution in the interstitial space between polymer rods exhibits the hallmarks of the RTIL structure near charged surfaces; i.e., cations (C₂mim⁺) and anions (TfO^-) form alternating layers around the polymer rods and the charge on the rod is overscreened by the ionic layer surrounding it. The distinct ordering of ions suggests the formation of a long-range "electrostatic network" in the ion gel, which may contribute to its mechanical cohesion and high modulus. The dynamics of both C₂mim⁺ and TfO^- ions slow down due to the fact that some C₂mim⁺ ions become associated with the sulfonate groups of the polymer rod on nanosecond time scales, which hinders the dynamics of all ions in the gel. C₂mim⁺ and TfO^- ion diffusion in the gel are only 2-10 times slower than in bulk RTILs, which is still much faster than, e.g., Li ions in typical ion conducting polymers. This fast ion transport combined with strong mechanical cohesion open up exciting opportunities for application of these gels in electrochemical devices including Li-metal batteries.

Vibrational Spectroscopy of Ionic Liquids

Vibrational spectroscopy has continued use as a powerful tool to characterize ionic liquids since the literature on room temperature molten salts experienced the rapid increase in number of publications in the 1990's. In the past years, infrared (IR) and Raman spectroscopies have provided insights on ionic interactions and the resulting liquid structure in ionic liquids. A large body of information is now available concerning vibrational spectra of ionic liquids made of many different combinations of anions and cations, but reviews on this literature are scarce. This review is an attempt at filling this gap. Some basic care needed while recording IR or Raman spectra of ionic liquids is explained. We have reviewed the conceptual basis of theoretical frameworks which have been used to interpret vibrational spectra of ionic liquids, helping the reader to distinguish the scope of application of different methods of calculation. Vibrational frequencies observed in IR and Raman spectra of ionic liquids based on different anions and cations are discussed and eventual disagreements between different sources are critically reviewed. The aim is that the reader can use this information while assigning vibrational spectra of an ionic liquid containing another particular combination of anions and cations. Different applications of IR and Raman spectroscopies are given for both pure ionic liquids and solutions. Further issues addressed in this review are the intermolecular vibrations that are more directly probed by the low-frequency range of IR and Raman spectra and the applications of vibrational spectroscopy in studying phase transitions of ionic liquids.

π+–π+ stacking of imidazolium cations enhances molecular layering of room temperature ionic liquids at their interfaces

The interfacial structure of room temperature ionic liquids (RTILs) is governed by the competing effects of the randomization due to the molecular polarizability and the ordered structure stabilized by π⁺–π⁺ interactions between the cationic molecules of RTILs.

Induced polarization restricts the conformational distribution of a light-harvesting molecular triad in the ground state

Free energy surface of the light-harvesting triad employing a non-polarizable force field (NFF) and a polarizable force field (PFF) shows that induced polarization limits the motion of rotation about chemical bonds as well as bending at the porphyrin, which are prominent using the NFF, thus limiting the conformational space of the triad.

Capacitive performance of amino acid ionic liquid electrolyte-based supercapacitors by molecular dynamics simulation

Capacitance–electric potential curves of amino acid ionic liquid electrolyte-based supercapacitors.

F. Gold nanoparticles optical properties induced by water and an ionic liquid (bmimBF4) inside cationic reverse micelles

Polar solvent interactions at the interface are a key factor affecting the final properties of nanoparticles, including size, shape and crystallinity.

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.

Charged, dipolar soft matter systems from a combined microscopic-mesoscopic viewpoint

As an example of charged, dipolar soft matter, the ionic liquid 1-ethyl-3-methyl-imidazolium dicyanamide is studied by coarse-grained molecular dynamics simulations. We focus on the link between microscopic and mesoscopic properties for both structure and dynamics. Thereby, the generalized Kirkwood g K-factor plays a central role in establishing this link which is not possible on the basis of molecular hydrodynamics. The decoupling between translational and rotational motion is indicative of the dynamical heterogeneity in ionic liquids.

Computational and Experimental Study of Li-Doped Ionic Liquids at Electrified Interfaces

We evaluate the influence of Li-salt doping on the dynamics, capacitance, and structure of three ionic liquid electrolytes, [pyr14][TFSI], [pyr13][FSI], and [EMIM][BF4], using molecular dynamics and polarizable force fields. In this respect, our focus is on the properties of the electric double layer (EDL) formed by the electrolytes at the electrode surface as a function of surface potential (Ψ). The rates of EDL formation are found to be on the order of hundreds of picoseconds and only slightly influenced by the addition of Li-salt. The EDLs of three electrolytes are shown to have different energy storage capacities, which we relate to the EDL formation free energy. The differential capacitance obtained from our computations exhibits asymmetry about the potential of zero charge and is consistent with the camel-like profiles noted from mean field theories and experiments on metallic electrodes. The introduction of Li-salt reduces the noted asymmetry in the differential capacitance profile. Complementary experimental capacitance measurements have been made on our three electrolytes in their neat forms and with Li-salt. The measurements, performed on glassy carbon electrodes, produce U-like profiles, and Li-salt doping is shown to strongly affect capacitance at high magnitudes of Ψ. Differences in the theoretical and experimental shapes and magnitudes of capacitance are rationalized in terms of the electrode surface and pseudocapacitive effects. In both neat and Li-doped liquids, the details of the computational capacitance profile are well described by Ψ-induced changes in the density and molecular orientation of ions in the molecular layer closest to the electrode. Our results suggest that the addition of Li+ induces disorder in the EDL, which originates from the strong binding of anions to Li+. An in-depth analysis of the distribution of Li+ in the EDL reveals that it does not readily enter the molecular layer at the electrode surface, preferring instead to be localized farther away from the surface in the second molecular layer. This behavior is validated through an analysis of the free energy of Li+ solvation as a function of distance from the electrode. Free energy wells are found to coincide with localized concentrations of Li+, the depths of which increase with Ψ and suggest a source of impedance for Li+ to reach the electrode. Finally, we make predictions of the specific energy at ideal graphite utilizing the computed capacitance and previously derived electrochemical windows of the liquids.

Photoinduced and Thermal Relaxation in Surface-Grafted Azobenzene-Based Monolayers: A Molecular Dynamics Simulation Study

Extensive atomistic molecular dynamics simulations have been employed to study the structure and molecular orientational relaxation of azobenzene-based monolayers grafted to a solid substrate. Systems with surface coverage of 0.6 nm(2)/molecule were investigated over a wide temperature range ranging from 298 K, where the mesogens show local ordering and the monolayer dynamics was found to be glassy, up to 700 K, where the azobenzene groups have a nearly isotropic orientational distribution, with a subnanosecond characteristic orientational relaxation time scale. Biased simulations that model single-molecule thermal excitation and conformational isomerization have been conducted to obtain insight into the mechanisms for photoinduced athermal fluidization and monolayer reorganization observed experimentally in this system. Our simulations clearly indicate that trans-cis conformational isomerization transitions of azobenzene units can lead to reorientation of mesogens and to the formation of a monolayer with strong macroscopic in-plane nematic order. While local heating created by excitation process can facilitate this process, thermal excitation alone is not sufficient to induce ordering in the monolayer. Instead, the work done by a molecule undergoing cis-trans isomerization on the cage of neighboring molecules is the key mechanism for photofluidization and orientational ordering in dMR monolayers exposed to linearly polarized light leading to relaxation dynamics that can be described in terms of higher effective temperature. The obtained simulation results are discussed in light of recent experimental data reported for these systems.

Computer Simulations of Ion Transport in Polymer Electrolyte Membranes

Understanding the mechanisms and optimizing ion transport in polymer membranes have been the subject of active research for more than three decades. We present an overview of the progress and challenges involved with the modeling and simulation aspects of the ion transport properties of polymer membranes. We are concerned mainly with atomistic and coarser level simulation studies and discuss some salient work in the context of pure binary and single ion conducting polymer electrolytes, polymer nanocomposites, block copolymers, and ionic liquid-based hybrid electrolytes. We conclude with an outlook highlighting future directions.

Additive polarizabilities in ionic liquids

An extended designed regression analysis of experimental data on density and refractive indices of several classes of ionic liquids yielded statistically averaged atomic volumes and polarizabilities of the constituting atoms. These values can be used to predict the molecular volume and polarizability of an unknown ionic liquid as well as its mass density and refractive index. Our approach does not need information on the molecular structure of the ionic liquid, but it turned out that the discrimination of the hybridization state of the carbons improved the overall result. Our results are not only compared to experimental data but also to quantum-chemical calculations. Furthermore, fractional charges of ionic liquid ions and their relation to polarizability are discussed.

Self-Consistent Determination of Atomic Charges of Ionic Liquid through a Combination of Molecular Dynamics Simulation and Density Functional Theory

A self-consistent scheme is developed to determine the atomic partial charges of ionic liquid. Molecular dynamics (MD) simulation was conducted to sample a set of ion configurations, and these configurations were subject to density functional theory (DFT) calculations to determine the partial charges. The charges were then averaged and used as inputs for the subsequent MD simulation, and MD and DFT calculations were repeated until the MD results are not altered any more. We applied this scheme to 1,3-dimethylimidazolium bis(trifluoromethylsulfonyl) imide ([C1mim][NTf2]) and investigated its structure and dynamics as a function of temperature. At convergence, the average ionic charges were ±0.84 e at 350 K due to charge transfer among ions, where e is the elementary charge, while the reduced ionic charges do not affect strongly the density of [C1mim][NTf2] and radial distribution function. Instead, major effects are found on the energetics and dynamics, with improvements of the overestimated heat of vaporization and the too slow motions of ions observed in MD simulations using commonly used force fields.

Fused coarse-grained model of aromatic ionic liquids and their behaviour at electrodes

A fused coarse-grained model of aromatic ionic liquids 1-alkyl-3-methylimidazoliums tetrafluoroborate ([C_nMIM⁺][BF₄⁻]) has been constructed.

Effect of water on the structure of a prototype ionic liquid

The influence of water on the structure of a prototype ionic liquid (IL) 1-octyl-3-methylimidazolium tetrafluoroborate (C₈mimBF₄) is examined in the IL-rich regime using high-energy X-ray diffraction (HEXRD) and molecular dynamics (MD) simulations.

Molecular simulations of imidazolium-based tricyanomethanide ionic liquids using an optimized classical force field

An atomistic force field is optimized to accurately predict the equilibrium and transport properties of technologically important imidazolium-based tricyanomethanide ionic liquids.

The influence of cations on lithium ion coordination and transport in ionic liquid electrolytes: a MD simulation study

The dynamical and structural properties in two ionic liquid electrolytes (ILEs) based on [emim][TFSI] and [pyr₁₃][TFSI] were compared as a function of LiTFSI salt concentrations using atomistic molecular dynamics (MD) simulations.

Thermophysical properties of imidazolium tricyanomethanide ionic liquids: experiments and molecular simulation

The low-viscous tricyanomethanide ([TCM]⁻)-based ionic liquids (ILs) are gaining increasing interest as attractive fluids for a variety of industrial applications.

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.

Ab Initio Simulations and Electronic Structure of Lithium-Doped Ionic Liquids: Structure, Transport, and Electrochemical Stability

Density functional theory (DFT), density functional theory molecular dynamics (DFT-MD), and classical molecular dynamics using polarizable force fields (PFF-MD) are employed to evaluate the influence of Li(+) on the structure, transport, and electrochemical stability of three potential ionic liquid electrolytes: N-methyl-N-butylpyrrolidinium bis(trifluoromethanesulfonyl)imide ([pyr14][TFSI]), N-methyl-N-propylpyrrolidinium bis(fluorosulfonyl)imide ([pyr13][FSI]), and 1-ethyl-3-methylimidazolium boron tetrafluoride ([EMIM][BF4]). We characterize the Li(+) solvation shell through DFT computations of [Li(Anion)n]((n-1)-) clusters, DFT-MD simulations of isolated Li(+) in small ionic liquid systems, and PFF-MD simulations with high Li-doping levels in large ionic liquid systems. At low levels of Li-salt doping, highly stable solvation shells having two to three anions are seen in both [pyr14][TFSI] and [pyr13][FSI], whereas solvation shells with four anions dominate in [EMIM][BF4]. At higher levels of doping, we find the formation of complex Li-network structures that increase the frequency of four anion-coordinated solvation shells. A comparison of computational and experimental Raman spectra for a wide range of [Li(Anion)n]((n-1)-) clusters shows that our proposed structures are consistent with experiment. We then compute the ion diffusion coefficients and find measures from small-cell DFT-MD simulations to be the correct order of magnitude, but influenced by small system size and short simulation length. Correcting for these errors with complementary PFF-MD simulations, we find DFT-MD measures to be in close agreement with experiment. Finally, we compute electrochemical windows from DFT computations on isolated ions, interacting cation/anion pairs, and liquid-phase systems with Li-doping. For the molecular-level computations, we generally find the difference between ionization energy and electron affinity from isolated ions and interacting cation/anion pairs to provide upper and lower bounds, respectively, to experiment. In the liquid phase, we find the difference between the lowest unoccupied and highest occupied electronic levels in pure and hybrid functionals to provide lower and upper bounds, respectively, to experiment. Li-doping in the liquid-phase systems results in electrochemical windows little changed from the neat systems.

Do TFSA Anions Slither? Pressure Exposes the Role of TFSA Conformational Exchange in Self-Diffusion

Multinuclear ((1)H, (2)H, and (19)F) magnetic resonance spectroscopy techniques as functions of temperature and pressure were applied to the study of selectively deuterated 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)amide (EMIM TFSA) ionic liquid isotopologues and related ionic liquids. For EMIM TFSA, temperature-dependent (2)H T1 data indicate stronger electric field gradients in the alkyl chain region compared to the imidazolium ring. Most significantly, the pressure dependences of the EMIM and TFSA self-diffusion coefficients revealed that the displacements of the cations and anions are independent, with diffusion of the TFSA anions being slowed much more by increasing pressure than for the EMIM cations, as shown by their respective activation volumes (28.8 ± 2.5 cm(3)/mol for TFSA vs 14.6 ± 1.3 cm(3)/mol for EMIM). Increasing pressure may lower the mobility of the TFSA anion by hindering its interconversion between trans and cis conformers, a process that is coupled to diffusion according to published molecular dynamics simulations. Measured activation volumes (ΔV(‡)) for ion self-diffusion in EMIM bis(fluoromethylsulfonyl)amide and EMIM tetrafluoroborate support this hypothesis. In addition, (2)H T1 data suggest increased ordering with increasing pressure, with two T1 regimes observed for the MD3 and D2 isotopologues between 0.1-100 and 100-250 MPa, respectively. The activation volumes for T1 were 21 and 25 cm(3)/mol (0-100 MPa) and 11 and 12 cm(3)/mol (100-250 MPa) for the MD3 and D2 isotopologues, respectively.

Capacitive Energy Storage: Current and Future Challenges

Capacitive energy storage devices are receiving increasing experimental and theoretical attention due to their enormous potential for energy applications. Current research in this field is focused on the improvement of both the energy and the power density of supercapacitors by optimizing the nanostructure of porous electrodes and the chemical structure/composition of the electrolytes. However, the understanding of the underlying correlations and the mechanisms of electric double layer formation near charged surfaces and inside nanoporous electrodes is complicated by the complex interplay of several molecular scale phenomena. This Perspective presents several aspects regarding the experimental and theoretical research in the field, discusses the current atomistic and molecular scale understanding of the mechanisms of energy and charge storage, and provides a brief outlook to the future developments and applications of these devices.