Diffusion and structure of propylene carbonate-metal salt electrolyte solutions for post-lithium-ion batteries: From experiment to simulation

The diffusion of cations in organic solvent solutions is important for the performance of metal-ion batteries. In this article, pulsed field gradient nuclear magnetic resonance experiments and fully atomistic molecular dynamic simulations were employed to study the temperature-dependent diffusive behavior of various liquid electrolytes representing 1M propylene carbonate solutions of metal salts with bis(trifluoromethylsulfonyl)imide (TFSI-) or hexafluorophosphate (PF6-) anions commonly used in lithium-ion batteries and beyond. The experimental studies revealed the temperature dependence of the diffusion coefficients for the propylene carbonate (PC) solvent and for the anions following an Arrhenius type of behavior. It was observed that the PC molecules are the faster species. For the monovalent cations (Li+, Na+, K+), the PC solvent diffusion was enhanced as the cation size increased, while for the divalent cations (Mg2+, Ca2+, Sr2+, Ba2+), the opposite trend was observed, i.e., the diffusion coefficients decreased as the cation size increased. The anion diffusion in LiTFSI and NaTFSI solutions was found to be similar, while in electrolytes with divalent cations, a decrease in anion diffusion with increasing cation size was observed. It was shown that non-polarizable charge-scaled force fields could correspond perfectly to the experimental values of the anion and PC solvent diffusion coefficients in salt solutions of both monovalent (Li+, Na+, K+) and divalent (Mg2+, Ca2+, Sr2+, Ba2+) cations at a range of operational temperatures. Finally, after calculating the radial distribution functions between cations, anions, and solvent molecules, the increase in the PC diffusion coefficient established with the increase in cation size for monovalent cations was clearly explained by the large hydration shell of small Li+ cations, due to their strong interaction with the PC solvent. In solutions with larger monovalent cations, such as Na+, and with a smaller solvation shell of PC, the PC diffusion is faster due to more liberated solvent molecules. In the salt solutions with divalent cations, both the anion and the PC diffusion coefficients decreased as the cation size increased due to an enhanced cation-anion coordination, which was accompanied by an increase in the amount of PC in the cation solvation shell due to the presence of anions.

Temperature-dependent differential capacitance of an ionic liquid-graphene-based supercapacitor

One of the critical factors affecting the performance of supercapacitors is thermal management. The design of supercapacitors that operate across a broad temperature range and at high charge/discharge rates necessitates understanding the correlation of the molecular characteristics of the device (such as interfacial structure and inter-ionic and ion-electrode interactions) with its macroscopic properties. In this study, we use molecular dynamics (MD) simulations to investigate the influence of Joule heating on the structure and dynamics of the ionic liquid (IL)/graphite-based supercapacitors. The temperature-dependent electrical double layer (EDL) and differential capacitance-potential (CD-V) curves of two different ([Bmim][BF4] and [Bmim][PF6]) IL-graphene pairs were studied under various thermal gradients. For the [Bmim][BF4] system, the differential capacitance curves transition from 'U' to bell shape under an applied thermal gradient (∇T) in the range from 3.3 K nm-1 to 16.7 K nm-1. Whereas in [Bmim][PF6], we find a positive dependence of differential capacitance with ∇T with a U-shaped CD-V curve. We examine changes in the EDL structure and screening potential (ϕ(z)) as a function of ∇T and correlate them with the trends observed in the CD-V curve. The identified correlation between the interfacial charge density and differential capacitance with thermal gradient would be helpful for the molecular design of the IL-electrode interface in supercapacitors or other chemical engineering applications.

Size Matters: A Computational Study of Hydrogen Absorption in Ionic Liquids

We combined both density functional theory and classical molecular dynamics simulations to investigate the molecular mechanisms governing hydrogen solvation in a total of 12 ionic liquids. Overall, the analysis of the structural properties under high temperature and pressure conditions revealed weak interactions between hydrogen and the ionic liquids, with a slight preference of this gas to be placed at the apolar domains. Interestingly, those ionic liquids comprising nitrate anions allow the accommodation of hydrogen molecules also in the polar areas. The study of the hydrogen velocity autocorrelation functions supports this observation. In addition, the structure of all of the tested ionic liquids was almost insensitive to the addition of hydrogen, so the available free volume and cavity formation are presumably the most important factors affecting solubility.

Polypropylene carbonate-based electrolytes as model for a different approach towards improved ion transport properties for novel electrolytes

Linear poly(alkylene carbonates) such as polyethylene carbonate (PEC) and polypropylene carbonate (PPC) have gained increasing interest due to their remarkable ion transport properties such as high Li⁺ transference numbers. The cause of these properties is not yet fully understood which makes it challenging to replicate them in other polymer electrolytes. Therefore, it is critical to understand the underlying mechanisms in polycarbonate electrolytes such as PPC. In this work we present insights from impedance spectroscopy, transference number measurements, PFG-NMR, IR and Raman spectroscopy as well as molecular dynamics simulations to address this issue. We find that in addition to plasticization, the lithium ion coordination by the carbonate groups of the polymer is weakened upon gelation, leading to a rapid exhange of the lithium ion solvation shell and consequently a strong increase of the conductivity. Moreover, we study the impact of the anions by employing different conducting salts. Interestingly, while the total conductivity decreases with increasing anion size, the reverse trend can be observed for the lithium ion transference numbers. Via our holistic approach, we demonstrate that this behavior can be attributed to differences in the collective ion dynamics.

Treating Polarization Effects in Charged and Polar Bio-Molecules Through Variable Electrostatic Parameters

Polarization plays important roles in charged and hydrogen bonding containing systems. Much effort ranging from the construction of physics-based models to quantum mechanism (QM)-based and machine learning (ML)-assisted models have been devoted to incorporating the polarization effect into the conventional force fields at different levels, such as atomic and coarse grained (CG). The application of polarizable force fields or polarization models was limited by two aspects, namely, computational cost and transferability. Different from physics-based models, no predetermining parameters were required in the QM-based approaches. Taking advantage of both the accuracy of QM calculations and efficiency of molecular mechanism (MM) and ML, polarization effects could be treated more efficiently while maintaining the QM accuracy. The computational cost could be reduced with variable electrostatic parameters, such as the charge, dipole, and electronic dielectric constant with the help of linear scaling fragmentation-based QM calculations and ML models. Polarization and entropy effects on the prediction of partition coefficient of druglike molecules are demonstrated by using both explicit or implicit all-atom molecular dynamics simulations and machine learning-assisted models. Directions and challenges for future development are also envisioned.

Impact of Dispersion Force Schemes on Liquid Systems: Comparing Efficiency and Drawbacks for Well-Targeted Test Cases

First-principles molecular dynamics (FPMD) calculations were performed on liquid GeSe₄ with the aim of inferring the impact of dispersion (van der Waals, vdW) forces on the structural properties. Different expressions for the dispersion forces were employed, allowing us to draw conclusions on their performances in a comparative fashion. These results supersede previous FPMD calculations obtained in smaller systems and shorter time trajectories by providing data of unprecedented accuracy. We obtained a substantial agreement with experiments for the structure factor regardless of the vdW scheme employed. This objective was achieved by using (in addition to FPMD with no dispersion forces) a selection of vdW schemes available within density functional theory. The first two are due to Grimme, D2 and D3, and the third one is devised within the so-called maximally localized Wannier functions approach (MLWF). D3 results feature a sizeable disagreement in real space with D2 and MLWF in terms of the partial and total pair correlation functions as well as the coordination numbers. More strikingly, total and partial structure factors calculated with D3 exhibit an unexpected sharp increase at low k. This peculiarity goes along with large void regions within the network, standing for a phase separation of indecipherable physical meaning. In view of these findings, further evidence of unconventional structural properties found by employing D3 is presented by relying on results obtained for a complex ionic liquid supported on a solid surface. The novelty of our study is multifold: new, reliable FPMD data for a prototypical disordered network system, convincing agreement with experimental data and assessment of the impact of dispersion forces, with emphasis on the intriguing behavior of one specific recipe and the discovery of common structural features shared by drastically dissimilar physical systems when the D3 vdW scheme is employed.

Rational Design of Solid Polymer Electrolyte Based on Ionic Liquid Monomer for Supercapacitor Applications via Molecular Dynamics Study

The safety concern arising from flammable liquid electrolytes used in batteries and supercapacitors drives technological advances in solid polymer electrolytes (SPEs) in which flammable organic solvents are absent. However, there is always a trade-off between the ionic conductivity and mechanical properties of SPEs due to the lack of interaction between the ionic liquid and polymer resin. The inadequate understanding of SPEs also limits their future exploitation and applications. Herein, we provide a complete approach to develop a new SPE, consisting of a cation (monomer), anion and hardener from ions-monomers using molecular dynamics (MD) simulations. The results show that the strong solid-liquid interactions between the SPE and graphene electrode lead to a very small gap of ∼5.5 Å between the components of SPE and electrode, resulting in a structured solid-to-liquid interface, which can potentially improve energy storage performance. The results also indicated the critical role of the mobility of free-standing anions in the SPE network to achieve high ionic conductivity for applications requiring fast charge/discharge. In addition, the formations of hardener-depleted regions and cation-anion-poor/rich regions near the uncharged/charged electrode surfaces were observed at the molecular level, providing insights for rationally designing the SPEs to overcome the boundaries for further breakthroughs in energy storage technology.

Probing the distribution of ionic liquid mixtures at charged and neutral interfaces via simulations and lattice-gas theory

Room temperature ionic liquid solutions confined between neutral and charged surfaces are investigated by means of atomistic Molecular Dynamics simulations. We study 1-ethyl-3-methylimidazolium dicyanamide ([EMIm]⁺[DCA]^-) in water or dimethylsulfoxide (DMSO) mixtures in confinement between two interfaces. The analysis is based on the comparison of the molecular species involved and the charged state of the surfaces. Focus is given on the influence of different water/DMSO concentrations on the microstructuring and accumulation of each species. Thermodynamic aspects, such as the entropic contributions in the observed trends are obtained from the simulations using a lattice-gas theory. The results clearly underline the differences in these properties for the water and DMSO mixtures and unravel the underlying mechanisms and inherent details. We were able to pinpoint the importance of the size and the relative permittivity of the molecules in guiding their microstructuring in the vicinity of the surfaces, as well as their interactions with the latter, i.e. the solute-surface interactions. The influence of water and DMSO on the overscreening at charged interfaces is also discussed. The analysis on the molecular accumulation at the interfaces allows us to predict whether the accumulation is entropy or enthalpy driven, which has an impact in the removal of the molecular species from the surfaces. We discuss the impact of this work in providing an essential understanding towards a careful design of electrochemical elements.

Recent Developments in Polarizable Molecular Dynamics Simulations of Electrolyte Solutions

Polarizable molecular dynamics simulations are a fast progressing field in the scientific research of ionic liquids. The fundamentals of polarizable simulations, as well as their application to ionic liquids, were summarized in a review [Bedrov, D.; Piquemal, J.-P.; Borodin, O.; MacKerell, Jr., A. D.; Roux, B.; Schröder, C. Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields. Chem. Rev. 2019, 119, 7940–7995] in 2019. Since then, new methods to treat intermolecular interaction of induced dipoles in these highly charged systems were developed. This concerns the damping of these interactions and additional charge transfer as well as the prediction of ionic materials with ultrahigh refractive indices. In addition to the progress of the polarizable force fields, also thermostats and barostats for polarizable simulations evolved recently.

On the Role of Water in the Formation of a Deep Eutectic Solvent Based on NiCl2·6H2O and Urea

The metal-based deep eutectic solvent (MDES) formed by NiCl₂·6H₂O and urea in 1:3.5 molar ratio has been prepared for the first time and characterized from a structural point of view. Particular accent has been put on the role of water in the MDES formation, since the eutectic could not be obtained with the anhydrous form of the metal salt. To this end, mixtures at different water/MDES molar ratios (W) have been studied with a combined approach exploiting molecular dynamics and ab initio simulations, UV–vis and near-infra-red spectroscopies, small- and wide-angle X-ray scattering, and X-ray absorption spectroscopy measurements. In the pure MDES, a close packing of Ni²⁺ ion clusters forming oligomeric agglomerates is present thanks to the mediation of bridging chloride anions and water molecules. Conversely, urea poorly coordinates the metal ion and is mostly found in the interstitial regions among the Ni²⁺ ion oligomers. This nanostructure is disrupted upon the introduction of additional water, which enlarges the Ni–Ni distances and dilutes the system up to an aqueous solution of the MDES constituents. In the NiCl₂·6H₂O 1:3.5 MDES, the Ni²⁺ ion is coordinated on average by one chloride anion and five water molecules, while water easily saturates the metal solvation sphere to provide a hexa-aquo coordination for increasing W values. This multidisciplinary study allowed us to reconstruct the structural arrangement of the MDES and its aqueous mixtures on both short- and intermediate-scale levels, clarifying the fundamental role of water in the eutectic formation and challenging the definition at the base of these complex systems.

The metal-based deep eutectic solvent formed by NiCl₂·6H₂O and urea in 1:3.5 a molar ratio was prepared for the first time, and its aqueous mixtures were characterized from a structural point of view, highlighting the fundamental role of water in the eutectic formation.

Comparison between Ab Initio Molecular Dynamics and OPLS-Based Force Fields for Ionic Liquid Solvent Organization

OPLS-based force fields (FFs) have been shown to provide accurate bulk-phase properties for a wide variety of imidazolium-based ionic liquids (ILs). However, the ability of OPLS to reproduce an IL solvent structure is not as well-validated given the relative lack of high-level theoretical or experimental data available for comparison. In this study, ab initio molecular dynamics (AIMD) simulations were performed for three widely used ILs: the 1-butyl-3-methylimidazolium cation with chloride, tetrafluoroborate, or hexafluorophosphate anions, that is, [BMIM][Cl], [BMIM][BF₄], and [BMIM][PF₆], respectively, as a basis for further assessment of two unique IL FFs: the ±0.8 charge-scaled OPLS-2009IL FF and the OPLS-VSIL FF. The OPLS-2009IL FF employs a traditional all-atom functional form, whereas the OPLS-VSIL FF was developed using a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to reproducing hydrogen bonding. Detailed comparisons between AIMD and the OPLS FFs were made based on radial distribution functions (RDFs), combined distribution functions (CDFs), and spatial distribution functions (SDFs) to examine cation-anion interactions and π⁺-π⁺ stacking between the imidazolium rings. While both FFs were able to correctly capture the general solvent structure of these popular ILs, the OPLS-VSIL FF quantitatively reproduced interaction distances more accurately. In addition, this work provides further insights into the different short- and long-range structure patterns of these popular ILs.

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.

Ionic screening in bulk and under confinement

Recent experiments have shown that the repulsive force between atomically flat, like-charged surfaces confining room-temperature ionic liquids or concentrated electrolytes exhibits an anomalously large decay length. In our previous publication [J. Zeman, S. Kondrat, and C. Holm, Chem. Commun. 56, 15635 (2020)], we showed by means of extremely large-scale molecular dynamics simulations that this so-called underscreening effect might not be a feature of bulk electrolytes. Herein, we corroborate these findings by providing additional results with more detailed analyses and expand our investigations to ionic liquids under confinement. Unlike in bulk systems, where screening lengths are computed from the decay of interionic potentials of mean force, we extract such data in confined systems from cumulative charge distributions. At high concentrations, our simulations show increasing screening lengths with increasing electrolyte concentration, consistent with classical liquid state theories. However, our analyses demonstrate that-also for confined systems-there is no anomalously large screening length. As expected, the screening lengths determined for ionic liquids under confinement are in good quantitative agreement with the screening lengths of the same ionic systems in bulk. In addition, we show that some theoretical models used in the literature to relate the measured screening lengths to other observables are inapplicable to highly concentrated electrolytes.

The effect of hydration on the stability of ionic liquid crystals: MD simulations of [C14C1im]Cl and [C14C1im]Cl·H2O

The thermal range of the stability of Ionic Liquid Crystal (ILC) phases of imidazolium ILCs, and the type of the mesophase itself are affected by several molecular structural features, the two prominent ones being the alkyl chain length and the counter-anion. Hydration is also very important: monohydrate samples of 1-alkyl-3-methylimidazolium halides have a higher clearing point and a wider thermal range of the stability of the ionic smectic phase, compared with the analogous anhydrous sample. To understand the reasons, at a microscopic level, for such increased stability due to hydration, we run classical Molecular Dynamics (MD) simulations of a typical ionic liquid crystal, 1-tetradecyl-3-methylimidazolium chloride, and of its monohydrate form. We tested a full-charge non-polarizable force field and a scaled-charge version having the total charge of the ions scaled by a factor of 0.80. Comparison of the structural and dynamic properties with available experimental data reveals that the scaling of the charge by a factor of 0.80 results in a good agreement between simulated and experimental data and it sheds light on the microscopic mechanism responsible for the increased stability of the monohydrated phase. A hydrogen-bond network between water and the chloride anion is established in the ionic layer which increases the stability of the ionic layer; this in turn increases the nano-segregation between the ionic and hydrophobic layers which eventually produce an increased order of the alkylic layer as well.

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM+][OAc-]

Reactive molecular dynamics simulations enable a detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the empirical valence bond framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NNs). The physics-based terms ensure the correct asymptotic behavior of electrostatic, polarization, and dispersion interactions and are compatible with existing solvent force fields. NNs are utilized for a versatile description of short-range orbital interactions within the transition state region and accurate rendering of vibrational motion of the reacting complex. We demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium cation with acetate to form 1-ethyl-3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ∼1 kJ mol^-1 mean absolute error accuracy within the transition state region for the gas-phase complex. To characterize the solvent modulation of the reaction profile, we compute potentials of mean force for the gas-phase reaction as well as the reaction within a four-ion cluster and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes the formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

Xenon Diffusion in Ionic Liquids with Blurred Nanodomain Separation

The dynamics of xenon gas, loaded in a series of 1-alkyl-3-methylimidazolium based ionic liquids, probes the formation of increasingly blurred polar/apolar nanodomains as a function of the anion type and the cation chain length. Exploiting ¹²⁹ Xe NMR spectroscopy techniques, like Pulse Gradient Spin Echo (PGSE) and inversion recovery (IR), the diffusion motion and relaxation times are determined for 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [C_n C₁ im][TFSI]. A correlation between the ILs nano-structure and both xenon diffusivity and relaxation times, as well as chemical shifts, is outlined. Interestingly, comparison with previous results of the same properties in the homologous imidazolium chlorides and hexafluorophospate shows an opposite trend with the alkyl chain length. Classical molecular dynamics (MD) simulations are used to calculate the xenon and cation and anion diffusion coefficients in the same systems, including imidazolium cations with longer chains (n=4, 6, 8 … 20). An almost quantitative agreement with the experiments validates the MD simulations and, at the same time, provides the necessary structural and dynamic microscopic insights on the nano-segregation and diffusion of xenon in bistriflimide, chloride and hexafluorphosphate salts allowing to observe and rationalize the shaping effect of the cation in the nanostructure.

Efficient Parametrization of Force Field for the Quantitative Prediction of the Physical Properties of Ionic Liquid Electrolytes

The prediction of transport properties of room-temperature ionic liquids from nonpolarizable force field-based simulations has long been a challenge. The uniform charge scaling method has been widely used to improve the agreement with the experiment by incorporating the polarizability and charge transfer effects in an effective manner. While this method improves the performance of the force fields, this prescription is ad hoc in character; further, a quantitative prediction is still not guaranteed. In such cases, the nonbonded interaction parameters too need to be refined, which requires significant effort. In this work, we propose a three-step semiautomated refinement procedure based on (1) atomic site charges obtained from quantum calculations of the bulk condensed phase; (2) quenched Monte Carlo optimizer to shortlist suitable force field candidates, which are then tested using pilot simulations; and (3) manual refinement to further improve the accuracy of the force field. The strategy is designed in a sequential manner with each step improving the accuracy over the previous step, allowing the users to invest the effort commensurate with the desired accuracy of the refined force field. The refinement procedure is applied on N,N-diethyl-N-methyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide (DEME-TFSI), a front-runner as an electrolyte for electric double-layer capacitors and single-molecule-based devices. The transferability of the refined force field is tested on N,N-dimethyl-N-ethyl-N-methoxyethoxyethylammonium bis(trifluoromethanesulfonyl)imide (N_112,2O2O1-TFSI). The refined force field is found to be better at predicting both structural and transport properties compared to the uniform charge scaling procedure, which showed a discrepancy in the X-ray structure factor. The refined force field showed quantitative agreement with structural (density and X-ray structure factor) and transport properties-diffusion coefficients, ionic conductivity, and shear viscosity over a wide temperature range, building a case for the wide adoption of the procedure.

Modelling biocompatible ionic liquids based on organic acids and amino acids: challenges for computational models and future perspectives

In this short review I shall highlight the basic principle and the difficulties that arise in attempting the computational modeling of seemingly simple systems which hide an unexpected complexity. Biocompatible ionic liquids which are based on the coupling of organic or amino acid anions with metabolic cations such as cholinium are the target of this review. These substances have been the subject of intense research activities in the last few years and have attracted the attention of computational chemists. I shall show that the computational description of these substances is far from trivial and requires the use of sophisticated techniques in order to account for a surprisingly rich chemistry that is due to several phenomena such as polarization, charge transfer, proton transfer equilibria and tautomerization reactions.

A review on machine learning algorithms for the ionic liquid chemical space

There are thousands of papers published every year investigating the properties and possible applications of ionic liquids. Industrial use of these exceptional fluids requires adequate understanding of their physical properties, in order to create the ionic liquid that will optimally suit the application. Computational property prediction arose from the urgent need to minimise the time and cost that would be required to experimentally test different combinations of ions. This review discusses the use of machine learning algorithms as property prediction tools for ionic liquids (either as standalone methods or in conjunction with molecular dynamics simulations), presents common problems of training datasets and proposes ways that could lead to more accurate and efficient models.

Spatial-Decomposition Analysis of Electrical Conductivity in Mixtures of Ionic Liquid and Sodium Salt for Sodium-Ion Battery Electrolytes

Ionic liquid (IL)-based electrolytes are a promising material for the development of sodium-ion batteries, and their performance can be quantified by electrical conductivity. In this highly concentrated ionic system, the correlated motions of ion pairs are influential on the ionic transport properties. Herein, all-atom analyses are conducted through molecular dynamics simulations to bridge the macroscopically observable electrical conductivity with the molecular pictures of correlated motion of ion pairs. The analysis is applied to three mixtures of IL with sodium salt that are relevant to the electrolyte for a sodium-ion battery: [1-ethyl-3-methylimidazolium, Na][bis(fluorosulfonyl)amide] ([C₂C₁im, Na][FSA]), [N-methyl-N-propylpyrrolidinium, Na][FSA] ([C₃C₁pyrr, Na][FSA]), and [K, Na][FSA]. The computational results on electrical conductivities are in agreement with the experimental reports, and their dependency on temperature and sodium-ion composition is reproduced well. The overall contributions from cross-correlated motions are found to be negative in all the IL mixtures; thus, the total conductivities are less than their Nernst-Einstein estimates. The spatial view of cross-correlated motions is further obtained by decomposing the time correlation functions of velocities according to the distances between ion pairs. It is observed that ion pairs are moving in the same direction for ∼0.3 ps when they were initially within the first coordination shell, followed by motions toward opposite directions. The cross-correlation terms are also dissected into local and nonlocal components, and it is commonly seen for all the ion pairs that the local component is negative for cation-anion pairs and is positive for cation-cation and anion-anion pairs. The motions of ion pairs are accompanied by a "backflow" that manifests in the form of the nonlocal component whose sign is opposite to the corresponding local component. In fact, the contributions of the correlated motions of ions to the electrical conductivity are not localized to contact pairs and extend spatially beyond the first coordination shell of the cation-anion pairs.

Interplay of Structure and Dynamics in Lithium/Ionic Liquid Electrolytes: Experiment and Molecular Simulation

Despite their promising use in electrochemical and electrokinetic devices, ionic-liquid-based electrolytes often exhibit complex behavior arising from a subtle interplay of their structure and dynamics. Here, we report a joint experimental and molecular simulation study of such electrolytes obtained by mixing 1-butyl 3-methylimidazolium tetrafluoroborate with lithium tetrafluoroborate. More in detail, experiments consisting of X-ray scattering, pulsed field gradient NMR, and complex impedance spectroscopy are analyzed in the light of molecular dynamics simulations to probe the structural, dynamical, and electrochemical properties of this ionic-liquid-based electrolyte. Lithium addition promotes the nanostructuration of the liquid as evidenced from the appearance of a scattering prepeak that becomes more pronounced. Microscopically, using the partial structure factors determined from molecular dynamics, this prepeak is shown to correspond to the formation of well-ordered positive/negative charge series and also large aggregates (Li_n(BF₄)_4-m)^(4-m+n)-, which develop upon lithium addition. Such nanoscale ordering entails a drastic decrease in both the molecular mobility and ionic conductivity. In particular, the marked association of Li⁺ cations with four BF₄^- anions and long ion pairing times, which are promoted upon lithium addition, are found to severely hinder the Li⁺ transport properties.

A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization

Molecular simulations can elucidate atomistic-level mechanisms of key biological processes, which are often hardly accessible to experiment. However, the results of the simulations can only be as trustworthy as the underlying simulation model. In many of these processes, interactions between charged moieties play a critical role. Current empirical force fields tend to overestimate such interactions, often in a dramatic way, when polyvalent ions are involved. The source of this shortcoming is the missing electronic polarization in these models. Given the importance of such biomolecular systems, there is great interest in fixing this deficiency in a computationally inexpensive way without employing explicitly polarizable force fields. Here, we review the electronic continuum correction approach, which accounts for electronic polarization in a mean-field way, focusing on its charge scaling variant. We show that by pragmatically scaling only the charged molecular groups, we qualitatively improve the charge-charge interactions without extra computational costs and benefit from decades of force field development on biomolecular force fields.

Structure, Molecular Interactions, and Dynamics of Aqueous [BMIM][BF4] Mixtures: A Molecular Dynamics Study

Molecular dynamics simulations with many-body polarizable force fields were carried out to investigate the thermodynamic, structural, and dynamic properties of aqueous solutions of 1-butyl-3-methylimidazolium tetrafluoroborate ([bmim][BF₄]). The radial distribution functions exhibit well-defined features, revealing favored structural correlations between [bmim]⁺, [BF₄]^-, and H₂O. The addition of water is shown to alter ionic liquid structural organizations by replacing counterions in the coordination shells and disrupt the cation-anion network. At low water concentration, the majority of water molecules are isolated from each other and have lower average dipole moment than that in pure water. With increasing hydration level, while [bmim][BF₄] ionic network breaks up and becomes isolated ion pairs or free ions in the dilute limit, water begins to form clusters of increasing sizes and eventually forms a percolating network. As a result, the average water dipole moment increases and approaches its bulk value. Water is also observed to have a substantial influence on the dynamics of ionic liquids. At low water content, the cation and anion have similar diffusion coefficients due to the correlated ionic motion of long-lived ion pairs. As the water concentration increases, both ions exhibit greater mobility and faster rotations from the breakup of ionic network. Consequently, the ionic conductivity of [bmim][BF₄] aqueous solutions rises with increasing water composition.

Pressure in Molecular Simulations with Scaled Charges. 1. Ionic Systems

Charge scaling, rationalized as MDEC (molecular dynamics in electronic continuum) or ECC (electronic continuum correction), has become a widely used simple approach to how to avoid self-consistent induced dipoles yet approximately take into account the effects of electronic polarizability. It has been assumed that the continuum permittivity does not depend on density; in turn, pressure is calculated by standard formulas. In this work, we elaborate a complementary approximation of density-independent molecular polarizability and derive formulas for pressure corrections within the MDEC framework; real behavior lies between these two extremes. The pressure corrections for test ionic systems are huge and negative, leading to sizable densities in constant-pressure MDEC simulations. A comparison of MDEC results with equivalent polarizable systems gives a good pressure match for a crystal but very low MDEC pressures for ionic liquids. These results witness about the importance of a correct density dependence not only of continuum permittivity in MDEC simulations but also of polarizability in polarizable simulations.

Xenon Dynamics in Ionic Liquids: A Combined NMR and MD Simulation Study

The translational dynamics of xenon gas dissolved in room-temperature ionic liquids (RTILs) is revealed by ¹²⁹Xe NMR and molecular dynamics (MD) simulations. The dynamic behavior of xenon gas loaded in 1-alkyl-3-methylimidazolium chloride, [C_nC₁im]Cl (n = 6, 8, 10), and hexafluorophosphate, [C_nC₁im][PF₆] (n = 4, 6, 8, 10) has been determined by measuring the ¹²⁹Xe diffusion coefficients and NMR relaxation times. The analysis of the experimental NMR data demonstrates that, in these representative classes of ionic liquids, xenon motion is influenced by the length of the cation alkyl chain and anion type. ¹²⁹Xe spin–lattice relaxation times are well described with a monoexponential function, indicating that xenon gas in ILs effectively experiences a single average environment. These experimental results can be rationalized based on the analysis of classical MD trajectories. The mechanism described here can be particularly useful in understanding the separation and adsorption properties of RTILs.

Simulation of Wetting and Interfacial Behavior of Quaternary Ammonium and Phosphonium Ionic Liquid Nanodroplets Over Face-Centered Cubic Metal Surfaces

We studied the behavior of quaternary ammonium- and phosphonium-based ([N₂₂₂₅][NTf₂] and [P₂₂₂₅][NTf₂]) ionic liquid (IL) nanodroplets on copper (Cu) and platinum (Pt) metal surfaces using classical molecular dynamics simulations. Solid-liquid interactions were underpinned by studying the dynamic spreading, contact angle, time-dependent binding energies, mean-square displacements, number and charge density profiles, and orientational distribution of these two IL nanodroplets. In particular, the role and importance of different crystallographic facets of face-centered cubic metal surface [Cu(100), Cu(111), Pt(100), and Pt(111)] were investigated. The results indicate a vast variety of properties that are dictated highly by the structure of different crystallographic facets of the metal substrate. The nature of the facet (e.g., the atomic arrangement symmetry, the extent of closed packed, and the unit cell size) leads to an extended scale spreading of the ILs on the surface with the (111) index but not with the (100) index, while the nature of metals itself plays a role.

Accurate Diels-Alder Energies and Endo Selectivity in Ionic Liquids Using the OPLS-VSIL Force Field

Our recently developed optimized potentials for liquid simulations-virtual site ionic liquid (OPLS-VSIL) force field has been shown to provide accurate bulk phase properties and local ion-ion interactions for a wide variety of imidazolium-based ionic liquids. The force field features a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to hydrogen bonding interactions. In this study, the Diels-Alder reaction between cyclopentadiene and methyl acrylate was computationally investigated in the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate, [BMIM][PF₆], as a basis for the validation of the OPLS-VSIL to properly reproduce a reaction medium environment. Mixed ab initio quantum mechanics and molecular mechanics (QM/MM) calculations coupled to free energy perturbation and Monte Carlo sampling (FEP/MC) that utilized M06-2X/6-31G(d) and OPLS-VSIL gave activation free energy barriers of 14.9 and 16.0 kcal/mol for the endo-cis and exo-cis Diels-Alder reaction pathways, respectively (exptl. ΔH^‡ of 14.6 kcal/mol). The endo selectivity trend was correctly predicted with a calculated 73% endo preference. The rate and selectivity enhancements present in the endo conformation were found to arise from preferential hydrogen bonding with the exposed C4 ring hydrogen on the BMIM cation. Weaker electronic stabilization of the exo transition state was predicted. For comparison, our earlier ±0.8 charge-scaled OPLS-2009IL force field also yielded a ΔG^‡ of 14.9 kcal/mol for the favorable endo reaction pathway but did not adequately capture the highly organized solvent interactions present between the cation and Diels-Alder transition state.

Current Status of AMOEBA-IL: A Multipolar/Polarizable Force Field for Ionic Liquids

Computational simulations of ionic liquid solutions have become a useful tool to investigate various physical, chemical and catalytic properties of systems involving these solvents. Classical molecular dynamics and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations of IL systems have provided significant insights at the atomic level. Here, we present a review of the development and application of the multipolar and polarizable force field AMOEBA for ionic liquid systems, termed AMOEBA–IL. The parametrization approach for AMOEBA–IL relies on the reproduction of total quantum mechanical (QM) intermolecular interaction energies and QM energy decomposition analysis. This approach has been used to develop parameters for imidazolium– and pyrrolidinium–based ILs coupled with various inorganic anions. AMOEBA–IL has been used to investigate and predict the properties of a variety of systems including neat ILs and IL mixtures, water exchange reactions on lanthanide ions in IL mixtures, IL–based liquid–liquid extraction, and effects of ILs on an aniline protection reaction.

Influence of the ion size on the stability of the smectic phase of ionic liquid crystals

The thermotropic phase behavior of ionic liquids and ionic liquid crystals based on novel N-alkyl-3-methylpyridinium halides, trihalides and dichloroiodates was experimentally studied by polarized optical spectroscopy (POM) and differential scanning calorimetry (DSC) as well as by molecular dynamics (MD) simulation. In the experiments, the existence and thermal range of stability of the smectic phase of these ionic liquid crystals are found to strongly depend on the volume ratio between the cation and anion, that is their relative size. Only compounds with a relatively large volume ratio of the cation to anion, i.e., those with longer cationic alkyl chains and monoatomic halide anions, have a stable smectic A phase. Both melting points and clearing points increase with such a ratio. The MD simulation results qualitatively agree very well with the experimental data and provide molecular details which can explain the experimentally observed phenomena: the stronger van der Waals interactions from the longer alkyl chains and the stronger electrostatic interactions from the smaller anions with a higher charge density increase the stability of both the crystal phase and the smectic phase; this also prevents the ionic layers from easily mixing with the hydrophobic regions, a mechanism that ultimately leads to a nanosegregated isotropic liquid phase.

Solvate ionic liquids based on lithium bis(trifluoromethanesulfonyl)imide–glyme systems: coordination in MD simulations with scaled charges

Charge scaling in molecular dynamics simulations of lithium bis(trifluoromethanesulfonyl)imide–glyme solvate ionic liquids yields better agreement with experiments.

Predicting Melting Points of Biofriendly Choline-Based Ionic Liquids with Molecular Dynamics

In this work, we introduce a simulation-based method for predicting the melting point of ionic liquids without prior knowledge of their crystal structure. We run molecular dynamics simulations of biofriendly, choline cation-based ionic liquids and apply the method to predict their melting point. The root-mean-square error of the predicted values is below 24 K. We advocate that such precision is sufficient for designing ionic liquids with relatively low melting points. The workflow for simulations is available for everyone and can be adopted for any species from the wide chemical space of ionic liquids.

Transport Properties of Li-TFSI Water-in-Salt Electrolytes

Water-in-salts are a new family of electrolytes that may allow the development of aqueous Li-ion batteries. They have a structure that is reminiscent of ionic liquids, and they are characterized by a high concentration of ionic species. In this work, we study their transport properties and how they evolve with concentration by using molecular dynamic simulations. We first focus on the choice of the force field. By comparing the simulated viscosities and self-diffusion coefficients with experimental measurements, we select a set of parameters that reproduces well the transport properties. We then use the selected force field to study in detail the variations of the self and collective diffusivities of all the species as well as the transport number of the lithium ion. We show that correlations between ions and water play an important role over the whole concentration range. In the water-in-salt regime, the anions form a percolating network that reduces the cation-anion correlations and leads to rather large values for the transport number compared to other standard electrolytes.

Ion transport in backbone-embedded polymerized ionic liquids

We use atomistic computer simulations to examine ion-transport phenomena for backbone polymerized cationic liquids with bistrifluoromethylesulfonylimide (TFSI^-) counterions. We consider a system in which the polymerized cation moiety is the imidazolium ring and study the structural characteristics and ion mobilities for cases in which the cations are separated by four, six, and eight methylene units on the backbone. A pendant polymerized ionic liquid, 1-butyl-3-vinylimidazolium, is compared to the backbone series across ion coordination and hopping features. The anion diffusivity in backbone polymerized cationic liquids is found to decrease with increasing spacer length, which is shown to result from a decrease in intramolecular and intermolecular hopping frequencies due to an increasing distance separating imidazolium moieties. In comparison with pendant polymerized ionic liquids, we observe that the participation rates of intermolecular hopping events in the backbone polymers far exceed that of the pendant, and the intrapolymeric ionic coordination profile shows the TFSI^- of the pendant polymer with a high propensity for coordination by multiple imidazolium, compared with one monomer from a given polymer for the backbone series. Despite these differences, backbone polymerized ionic liquids are seen to possess correlated diffusivity and ion-association relaxation times, in a manner similar to the results observed in past studies for pendant variants.