Transport Properties of the Ionic Liquid 1-Ethyl-3-Methylimidazolium Chloride from Equilibrium Molecular Dynamics Simulation. The Effect of Temperature

Transfer learning for accurate description of atomic transport in Al-Cu melts

Machine learning interatomic potentials (MLIPs) provide an optimal balance between accuracy and computational efficiency and allow studying problems that are hardly solvable by traditional methods. For metallic alloys, MLIPs are typically developed based on density functional theory with generalized gradient approximation (GGA) for the exchange-correlation functional. However, recent studies have shown that this standard protocol can be inaccurate for calculating the transport properties or phase diagrams of some metallic alloys. Thus, optimization of the choice of exchange-correlation functional and specific calculation parameters is needed. In this study, we address this issue for Al-Cu alloys, in which standard Perdew-Burke-Ernzerhof (PBE)-based MLIPs cannot accurately calculate the viscosity and melting temperatures at Cu-rich compositions. We have built MLIPs based on different exchange-correlation functionals, including meta-GGA, using a transfer learning strategy, which allows us to reduce the amount of training data by an order of magnitude compared to a standard approach. We show that r2SCAN- and PBEsol-based MLIPs provide much better accuracy in describing thermodynamic and transport properties of Al-Cu alloys. In particular, r2SCAN-based deep machine learning potential allows us to quantitatively reproduce the concentration dependence of dynamic viscosity. Our findings contribute to the development of MLIPs that provide quantum chemical accuracy, which is one of the most challenging problems in modern computational materials science.

Electrokinetic properties of NaCl solution via molecular dynamics simulations with scaled-charge electrolytes

Scaling ionic charges has become an alternative to polarizable force fields for representing indirect charge transfer effects in molecular simulations. In our work, we apply molecular dynamics simulations to investigate the properties of NaCl aqueous solutions in homogeneous and confined media. We compare classical integer- and scaled-charge force fields for the ions. In the bulk, we validate the force fields by computing equilibrium and transport properties and comparing them with experimental data. Integer-charge ions overestimate dielectric saturation and ionic association. Both force fields present an excess in ion-ion correlation, which leads to a deviation in the ionic conductivity at higher ionic strengths. Negatively charged quartz is used to simulate the confinement effect. Electrostatic interactions dominate counter-ion adsorption. Full-charge ions have stronger and more defined adsorption planes. We obtain the electroosmotic mobility of the solution by combining the shear plane location from non-equilibrium simulations with the ionic distribution from equilibrium simulations. From the Helmholtz-Smoluchowski equation, the zeta potential and the streaming potential coupling coefficient are computed. From an atomic-scale perspective, our molecular dynamics simulations corroborate the hypothesis of maximum packing of the Stern layer, which results in a stable and non-zero zeta potential at high salinity. The scaled-charge model representation of both properties is in excellent qualitative and quantitative agreement with experimental data. With our work, we demonstrate how useful and precise simple scaled-charge models for electrolytes can be to represent complex systems, such as the electrical double layer.

Efficient single-run implementation of generalized Einstein relation to compute transport coefficients: A binary-based time sampling

A robust and simple implementation of the generalized Einstein formulation using single equilibrium molecular dynamics simulation is introduced to compute diffusion and shear viscosity. The unique features underlying this framework are as follows: (1) The use of a simple binary-based method to sample time-dependent transport coefficients results in a uniform distribution of data on a logarithmic time scale. Although we sample "on-the-fly," the algorithm is readily applicable for post-processing analysis. Overlapping same-length segments are not sampled as they indicate strong correlations. (2) Transport coefficients are estimated using a power law fitting function, a generalization of the standard linear relation, that accurately describes the long-time plateau. (3) The use of a generalized least squares (GLS) fitting estimator to explicitly consider correlations between fitted data points results in a reliable estimate of the statistical uncertainties in a single run. (4) The covariance matrix for the GLS method is estimated analytically using the Wiener process statistics and computed variances. (5) We provide a Python script to perform the fits and automate the procedure to determine the optimal fitting domain. The framework is applied to two fluids, binary hard sphere and a Lennard-Jones near the triple point, and the validity of the single-run estimates is verified against multiple independent runs. The approach should be applicable to other transport coefficients since the diffusive limit is universal to all of them. Given its rigor and simplicity, this methodology can be readily incorporated into standard molecular dynamics packages using on-the-fly or post-processing analysis.

Efficient Molecular Dynamics Simulations of Deep Eutectic Solvents with First-Principles Accuracy Using Machine Learning Interatomic Potentials

In recent years, deep eutectic solvents emerged as highly tunable and ecofriendly alternatives to common organic solvents and liquid electrolytes. In the present work, the ability of machine learning (ML) interatomic potentials for molecular dynamics (MD) simulations of these liquids is explored, showcasing a trained neural network potential for a 1:2 ratio mixture of choline chloride and urea (reline). Using the ML potentials trained on density functional theory data, MD simulations for large systems of thousands of atoms and nanosecond-long time scales are feasible at a fraction of the computational cost of the target first-principles simulations. The obtained structural and dynamical properties of reline from MD simulations using our machine learning models are in good agreement with the first-principles MD simulations and experimental results. Running a single MD simulation is highlighted as a general shortcoming of typical first-principles studies if the dynamic properties are investigated. Furthermore, velocity cross-correlation functions are employed to study the collective dynamics of the molecular components in reline.

Effect of Ion Mass on Dynamic Correlations in Ionic Liquids

Ionic liquids (ILs) are a class of liquid salts with distinct properties such as high ionic conductivity, low volatility, and a broad electrochemical window, making them appealing for use in energy storage applications. The ion-ion correlations are some of the key factors that play a critical role in the ionic conductivity of ILs. In this work, we present the study of the impact of ion mass on ion-ion correlations in ILs, applying a combination of broadband dielectric spectroscopy measurements and molecular dynamics simulations. We examined three ILs with the same cation but different anions to consider three different cases of cation-anion masses: M+ > M-, M+ ≈ M-, and M+ < M-. We applied the momentum conservation approach to estimate the contribution of distinct ion-ion correlations from experimental data and obtained good agreement with direct calculations of distinct ion-ion correlations from molecular dynamics simulations. Our findings reveal that relative ion mass has a strong effect on the distinct ion-ion correlations, leading to swapping of the relative amplitude of distinct cation-cation and anion-anion correlations.

Computation of Shear Viscosity by a Consistent Method in Equilibrium Molecular Dynamics Simulations: Applications to 1-Decene Oligomers

Accurate computation of shear viscosity is fundamental for describing fluid flow and designing and developing new processes. Poly-α-olefins (PAO's), particularly from 1-decene, have been applied to a variety of industrial processes. Recently, these molecules have been applied as carbon dioxide thickeners, enhancing carbon dioxide viscosity, which is important in carbon dioxide injection, either for enhanced oil recovery or sequestration in geological formations. For these applications, knowledge of the pure oligomer viscosity is crucial to design and operate the oligomer upstream pipelines before mixing them with carbon dioxide. Using Green-Kubo formalism with equilibrium molecular dynamics simulations, two methods are presented in the literature to generate the traceless, symmetric pressure tensor. In this work, we show that these two methods provide different values of shear viscosity, from the analysis of how the diagonal components of the traceless, symmetric pressure tensor are computed in each method. Then, we examine the consistency and correctness of each method: one is found to be consistent. The other is corrected by scaling the fluctuations of the diagonal components. Shear viscosities of supercritical carbon dioxide, vapor and liquid n-pentane, and liquid n-decane are computed to illustrate the analysis. We also apply the consistent method to compute the viscosity of 1-decene oligomers, including for the first time larger-than-dimer oligomers (trimer, tetramer, hexamer, and decamer).

Viscosity Prediction of High-Concentration Antibody Solutions with Atomistic Simulations

The computational prediction of the viscosity of dense protein solutions is highly desirable, for example, in the early development phase of high-concentration biopharmaceutical formulations where the material needed for experimental determination is typically limited. Here, we use large-scale atomistic molecular dynamics (MD) simulations with explicit solvation to de novo predict the dynamic viscosities of solutions of a monoclonal IgG1 antibody (mAb) from the pressure fluctuations using a Green-Kubo approach. The viscosities at simulated mAb concentrations of 200 and 250 mg/mL are compared to the experimental values, which we measured with rotational rheometry. The computational viscosity of 24 mPa·s at the mAb concentration of 250 mg/mL matches the experimental value of 23 mPa·s obtained at a concentration of 213 mg/mL, indicating slightly different effective concentrations (or activities) in the MD simulations and in the experiments. This difference is assigned to a slight underestimation of the effective mAb-mAb interactions in the simulations, leading to a too loose dynamic mAb network that governs the viscosity. Taken together, this study demonstrates the feasibility of all-atom MD simulations for predicting the properties of dense mAb solutions and provides detailed microscopic insights into the underlying molecular interactions. At the same time, it also shows that there is room for further improvements and highlights challenges, such as the massive sampling required for computing collective properties of dense biomolecular solutions in the high-viscosity regime with reasonable statistical precision.

Probing the Origin of Viscosity of Liquid Electrolytes for Lithium Batteries

Viscosity is an extremely important property for ion transport and wettability of electrolytes. Easy access to viscosity values and a deep understanding of this property remain challenging yet critical to evaluating the electrolyte performance and tailoring electrolyte recipes with targeted properties. We proposed a screened overlapping method to efficiently compute the viscosity of lithium battery electrolytes by molecular dynamics simulations. The origin of electrolyte viscosity was further comprehensively probed. The viscosity of solvents exhibits a positive correlation with the binding energy between molecules, indicating viscosity is directly correlated to intermolecular interactions. Salts in electrolytes enlarge the viscosity significantly with increasing concentrations while diluents serve as the viscosity reducer, which is attributed to the varied binding strength from cation-anion and cation-solvent associations. This work develops an accurate and efficient method for computing the electrolyte viscosity and affords deep insight into viscosity at the molecular level, which exhibits the huge potential to accelerate advanced electrolyte design for next-generation rechargeable batteries.

Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property

In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.

Tracing the Effect of Replacing [Gly]- with [Ala]- and Hydroxylation of [emim]+ on the Fine-Tuning of the Transport Properties of the Corresponding Amino Acid-Based Ionic Liquids Using MD Simulation

Natural amino acid-based ionic liquids (AAILs) composed of deprotonated amino acids, [AA]^-, as anions and hydroxylated imidazolium cations provide an eco-friendly nontoxic IL family with the growing number of chemical and biochemical revolutionary applications. In this paper, the transport properties of four AAILs composed of 1-(2-hydroxyethyl)-3-methylimidazolium ([HOemim]⁺) and 1-ethyl-3-methylimidazolium ([emim]⁺) cations with alaninate and glycinate anions were studied by molecular dynamics (MD) simulations. A nonpolarizable all-atom force field with the scaled charge (±0.8e) on each of the ions was applied and compared with the unit charge model in some cases. The tunable effects of the presence of the hydroxyl group in the side chain of the imidazolium cation, the type of amino acid anion, and the varied temperature on the dynamical behavior of AAILs were investigated in detail. The experimentally compatible trends of the simulated ionic self-diffusion coefficients, ionic conductivity, and ionicity were found to be inverse to the viscosity and ionic association of these ILs as [emim][Gly] > [emim][Ala] > [HOemim][Gly] > [HOemim][Ala]. The main reason behind these trends is the higher ability of the hydroxylated cation for the hydrogen-bond formation with [AA]^-. The mean square displacement (MSD), self-diffusion, and transference number of imidazolium cations are larger than those of [AA]^- anions, except in the case of [HOemim][Gly]. It was found that the activation energy for diffusion of [AA]^- is lower than that of [HOemim]⁺ but higher than that of [emim]⁺ in [HOemim][AA] and [emim][AA] ILs, respectively. The computed velocity autocorrelation function (VACF) showed that [Gly]^-, as the lightest ion, has the shortest mean collision time and velocity randomization time among the ions, especially in the [HOemim][Gly] IL. Replacing [emim]⁺ with [HOemim]⁺, similar to the effect of decreasing temperature, causes significant decreasing of the ionic self-diffusion and increasing of the well depth of the first minimum of the ionic VACFs. Current findings show that introducing suitable functional groups in the side chain of imidazolium cations can be a viable approach for efficient engineering design and fine-tuning of the transport properties of these AAILs.

Collectivity in ionic liquids: a temperature dependent, polarizable molecular dynamics study

We use polarizable molecular dynamics simulations to study the thermal dependence of both structural and dynamic properties of two ionic liquids sharing the same cation (1-ethyl-3-methylimidazolium). The linear temperature trend in the structure is accompanied by an exponential Arrhenius-like behavior of the dynamics. Our parameter-free Voronoi tessellation analysis directly casts doubt on common concepts such as the alternating shells of cations and anions and the ionicity. The latter tries to explain the physico-chemical properties of the ionic liquids based on the association and dissociation of an ion pair. However, cations are in the majority of both ion cages, around cations and around anions. There is no preference of a cation for a single anion. Collectivity is a key factor in the dynamic properties of ionic liquids. Consequently, collective rotation relaxes faster than single-particle rotations, and the activation energies for collective translation and rotation are lower than those of the single molecules.

Fast evaluation technique for the shear viscosity and ionic conductivity of electrolyte solutions

With the growing need to obtain ideal materials for various applications, there is an increasing interest in computational methods to rapidly and accurately search for materials. Molecular dynamics simulation is one of the successful methods used to investigate liquid electrolytes with high transport properties applied in lithium-ion batteries. However, further reduction in computational cost is required to find a novel material with the desired properties from a large number of combinations. In this study, we demonstrate an effective fast evaluation technique for shear viscosity and ionic conductivity by molecular dynamics simulation for an exhaustive search of electrolyte materials with high transport properties. The proposed model was combined with a short-time correlation function of the stress tensor and empirical relationships to address the issues of inefficient and uncertain evaluation by conventional molecular dynamics methods. Because we focus on liquid electrolytes consisting of organic solvents and lithium salts, our model requires dissociation ratio and effective diffusion size of lithium salts. Our method is applied to search for the compositional combinations of electrolytes with superior transport properties even at low temperatures. These results correlate well with experimental results.

Ion Correlation in Choline Chloride-Urea Deep Eutectic Solvent (Reline) from Polarizable Molecular Dynamics Simulations

In recent years, deep eutectic solvents (DESs) emerged as highly tunable and environmentally friendly alternatives to common ionic liquids and organic solvents. In this work, a polarizable model based on the CHARMM Drude polarizable force field is developed for a 1:2 ratio mixture of choline chloride/urea (reline) DES. To successfully reproduce the structure of the liquid as compared to first-principles molecular dynamics simulations, a damping factor was introduced to correct the observed over-binding between the chloride and the hydrogen bonding site of choline. Investigated radial distributions reveal the formation of hydrogen bonds between all the constituents of reline and similar interactions for chloride and urea's oxygen atoms, which could contribute to the melting point depression of the mixture. Predicted dynamic properties from our polarizable force field are in good agreement with experiments, showing significant improvements over nonpolarizable models. Similar to some ionic liquids, an oscillatory behavior in the velocity autocorrelation function of the anion is visible, which can be interpreted as a rattling motion of the lighter anion surrounded by the heavier cations. The obtained results for ionic conductivity of reline show some degree of correlated ion motion in this DES. However, a joint diffusion of ion pairs cannot be observed during the simulations.

Ion Transport in the EMITFSI/PVDF System at Different Temperatures: A Molecular Dynamics Simulation

We used all-atom molecular dynamics simulations to study the ion transport in the 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide/poly(vinylidene fluoride) (EMITFSI/PVDF) system with 40.05 wt % EMITFSI at different temperatures. The glass-transition temperature (T g = 204 K) of this system shows a good agreement with the experimental value (200 K). With the increase of temperature, the peaks of the pair correlation function show an increasing trend. Interestingly, the coordination numbers of ion pairs and the degree of independent ion motion are mainly affected by the binding energy between ion pairs as the temperature increases. In addition, the ion transport properties with increasing temperature can be studied by the ion-pair relaxation times, ion-pair lifetimes, and diffusion coefficients. The simulation results illustrate that the ion transport is intensified. Especially, the cations can always diffuse faster than the anions. The power law shows that mobilities of anions and cations are seen to exhibit a "superionic" behavior. With the increase of temperature, transference numbers of anions decrease first and then increase and transference numbers of cations show the opposite changes; ionic conductivity increases gradually; and viscosity decreases gradually, indicating that the diffusion resistance of ions decreases. In general, after adding PVDF into the EMITFSI system, the glass-transition temperature and viscosity increase, the ionic conductivity and degree of independent ion motion decrease, and diffusion coefficients of cations decrease faster than those of the anions.

Properties of aqueous 1,4-dioxane solution via molecular dynamics

Polyethers are promising compounds for the creation of electrochemical energy storage systems. The molecular dynamics method can facilitate the search of compounds that have the most potential. However, the application of this method requires verification of the force fields. We perform molecular dynamics calculations of the physical properties of the aqueous 1,4-dioxane solution (density, enthalpy of mixing, and viscosity) and compare them to the available experimental data. In addition, we confirm the idea that the solution structure depends on the dioxane molar fraction, proposed in the experiment of Takamuku et al. [J. Mol. Liq. 83(1-3), 163-177 (1999)]. The hydrogen bonds between dioxane and water are analyzed. The correlation between the excess viscosity and enthalpy of mixing is demonstrated.

Molecular Dynamics and Emerging Network Graphs of Interactions in Dinitrile-Based Li-Ion Battery Electrolytes

Advancements in battery research have shown interesting formulations of battery electrolytes that have helped improve the efficiency of Li-ion batteries over the decades. However, the quest for a safer and affordable battery electrolyte still proceeds with more unique formulations reported in the literature regularly. The dinitriles, especially adiponitrile and glutaronitrile, have caught the attention of the research community as part of this quest. In this work, we performed molecular dynamics simulations of dinitrile electrolytes with lithium bistrifluorosulfonimide (LiTFSI) as the electrolyte salt at varying concentrations and temperatures. On analysis of our simulations, we find that the densities of the mixtures follow the same trend as that of experimental values. The solvation properties were explored using the radial distribution functions. The connectivity of the Li⁺ with the dinitrile molecules and anions is established for all of the electrolyte concentrations using network graphs. We observe that the electrolytes form highly networked structures as the concentration increases without being affected by the rise in temperature. The networking of ionic interactions was quantified by calculating the average degree of each graph. Ionic conductivity calculations were computed using three methods: Nernst-Einstein relation, correlated method, and current autocorrelation function. We report the importance of accounting for the correlated motion of ions while estimating the ionic conductivity. The correlated conductivity and current autocorrelation function calculations provide a satisfactory estimation of the ionic conductivity compared to the experimental values.

Correlated Ion Transport and the Gel Phase in Room Temperature Ionic Liquids

Here we present a theory of ion aggregation and gelation of room temperature ionic liquids (RTILs). Based on it, we investigate the effect of ion aggregation on correlated ion transport-ionic conductivity and transference numbers-obtaining closed-form expressions for these quantities. The theory depends on the maximum number of associations a cation and anion can form and the strength of their association. To validate the presented theory, we perform molecular dynamics simulations on several RTILs and a range of temperatures for one RTIL. The simulations indicate the formation of large clusters, even percolating through the system under certain circumstances, thus forming a gel, with the theory accurately describing the obtained cluster distributions in all cases. However, based on the strength and lifetime of associations in the simulated RTILs, we expect free ions to dominate ionic conductivity despite the presence of clusters, and we do not expect the percolating cluster to trigger structural arrest in the RTIL.

Interfacial Polarization and Ionic Structure at the Ionic Liquid-Metal Interface Studied by Vibrational Spectroscopy and Molecular Dynamics Simulations

Ionic liquids (ILs) have both fundamental and practical value in interfacial science and electrochemistry. However, understanding their behavior near a surface is challenging because of strong Coulomb interactions and large and irregular ionic sizes, which affect both their structure and energetics. To understand this problem, we present a combined experimental and computational study using a vibrational probe molecule, 4-mercaptobenzonitrile, inserted at the junction between a metal and a variety of ILs. The vibrational frequency of the nitrile in the probe molecule reports on the local solvation environment and the electrostatic field at this junction. Within the ethylmethyl imidazolium (EMIM⁺) cation family of ILs, we varied the anions over a range of sizes and types. Complementing our surface spectroscopy, we also ran molecular dynamics simulations of these interfaces to better understand the ionic structures that produced the measured fields. The magnitude of the frequency shifts, and thereby fields, shows a general correlation with the size of anions, with larger anions corresponding to smaller fields. We find that the source of this correlation is partial intercalation of smaller anions into the probe monolayer, resulting in tighter packing of ionic layers near the surface. Larger anions reduce the overall lateral ion packing density near the surface, which reduces the net charge per unit area and explains the smaller observed fields. The insight from this work is important for developing a fundamental picture of concentrated electrolytes near interfaces and can help with designing ILs to create tailored electric fields near an electrode.

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.

Extension of transferable coarse-grained models to dicationic ionic liquids

In this study, we extended the previously developed coarse-grained (CG) models of mono-cationic ionic liquids (MILs) to di-cationic ILs (DILs). To achieve this purpose, the MD simulations in three different mapping schemes of CG were done and the results of RDF (as a structural property), density (as a volumetric property) and the diffusion coefficient (as a dynamical property) were compared with the corresponding results of the all-atom (AA) simulations for [C₅(mim)₂][BF₄]₂. The previously developed CG models for MILs with the least refinement in parameters were used to extend the CG models for DILs. Since, the first mapping scheme of the CG model showed the best agreement with the results of the AA simulations for the three mentioned studied properties, this scheme was selected to simulate DILs using the CG model. The transferability of the selected CG model to DILs was investigated by comparing the different volumetric, structural and dynamical properties of [C_n(mim)₂][BF₄]₂ (with n = 3, 6, 9, and 12) obtained from the CG model with those obtained using the corresponding atomistic simulations at different thermodynamic state points. The average deviation for the densities of the CG model with respect to the AA results is less than 2%. Furthermore, in both CG and AA models, the densities and isobaric expansion coefficients decrease with increasing temperature and linkage alkyl chain. The structural properties of the studied DILs, i.e. RDFs, nano segregation of domains, heterogeneity order parameters (HOPs) and angle distributions showed good agreements between the results of the CG and AA models. The CG-based calculated diffusion coefficients of the studied DILs at different temperatures showed that this model leads to faster dynamics with respect to the AA model due to the sacrifice of some degrees of freedom in this model. However, the trend of increasing diffusion coefficients with increasing temperature and linkage alkyl chain length is the same in both CG and AA models. Also, there are good agreements between the results of these two models for other dynamical properties, i.e. electrical conductivity, transference numbers and non-Gaussian parameter with increasing linkage alkyl chain and at various temperatures.

A molecular dynamics study on magnetic imidazolium-based ionic liquids: the effect of an external magnetic field

In this paper, we have reported a molecular dynamics (MD) study on the properties of three different magnetic imidazolium-based ionic liquids in the absence and presence of an external magnetic field. In this regard, the volumetric properties such as density and isobaric expansion coefficient, dynamical properties, namely, viscosity, mean square displacement of ions, diffusion coefficients, transport numbers of cations and anions, and electrical conductivity, and structural properties such as radial distribution function (RDF) and spatial distribution function (SDF) of [emim][FeCl4], [bmim][FeCl4] and [hmim][FeCl4] have been studied at different temperatures using molecular dynamics simulations. After studying the different volumetric, structural, and dynamical properties of the above-mentioned magnetic ILs in the absence of a magnetic field, we investigated the effect of an external magnetic field on the structural properties of one of these systems, i.e., [bmim][FeCl4]. In this regard, we established different contributions in the interactions between the external magnetic field and the studied magnetic ionic liquid (MIL). The number density profiles of the studied MIL before and after imposing an external magnetic field of 1.5 T showed a significant variation in the molecular distribution. The results indicated that the external magnetic field reduced the intensity of RDFs due to the reduction in the interactions between different ion sites as a result of changes in their orientations. After applying the external magnetic field, it was observed that due to the oppositely directed forces on the cations and anions, they moved in opposite directions. The snapshots showed that the static motion of the anion was smaller because of its small size. In the presence of an external magnetic field, the ions distributed more homogeneously compared to that observed in the absence of this field. The results of this study can be used in the rational and accurate design of viscomagnetic fluids and reaction systems in the presence and absence of magnetic fields.

The Effect of Concentration of Lithium Salt on the Structural and Transport Properties of Ionic Liquid-Based Electrolytes

Ionic liquids (ILs) are used as electrolytes in high-performance lithium-ion batteries, which can effectively improve battery safety and energy storage capacity. All atom molecular dynamics simulation and experiment were combined to investigate the effect of the concentration of lithium salt on the performance of electrolytes of four IL solvents ([C_nmim][TFSI] and [C_nmim][FSI], n = 2, 4). The IL electrolytes exhibit higher density and viscosity; meanwhile, larger lithium ion transfer numbers as the concentration of LiTFSI increases. Furthermore, in order to explore the effect of the concentration of lithium salt on the ionic associations of Li⁺ and anion of IL, the microstructures of the lithium salt in various IL electrolytes at different concentrations were investigated. The structural analysis indicated that strong bidentate and monodentate coordination was found between Li⁺ and anion of all IL electrolytes. Both cis and trans isomerism of [FSI]⁻ were observed in [FSI]⁻-type IL electrolyte systems. Furthermore, the existence of the ion cluster [Li[anion]_x]^(x−1)− in IL electrolytes and the cluster became more closed and compact as the concentration of LiTFSI increases.

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.

Contributions of force field interaction forms to Green-Kubo viscosity integral in n-alkane case

Decades of molecular simulation history proved that the Green-Kubo method for shear viscosity converges without any problems in atomic and simple molecular liquids, unlike liquids with high values of viscosity. In the case of highly viscous liquids, the time decomposition method was developed in 2015 by Maginn and co-authors [Y. Zhang, A. Otani, and E. J. Maginn, J. Chem. Theory Comput. 11, 3537-3546 (2015)] which allows us to improve the convergence of the Green-Kubo integral. In this paper, the contributions of intramolecular and intermolecular force field parts to the viscosity integral are discovered to gain the understanding of the Green-Kubo method. The n-alkanes from n-ethane to n-pentane at 330 K in the optimized potentials for liquid simulations-all atom force field are used as reference models. The dependencies of these contributions and decay times of the corresponding correlation functions on the chain length are observed. The nonequilibrium simulations are carried out to verify the Green-Kubo results. The obtained values of viscosity are compared with experimental data.

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields

Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.

Thermodynamic, structural and dynamic properties of ionic liquids [C4mim][CF3COO], [C4mim][Br] in the condensed phase, using molecular simulations

In this work a series of thermodynamic, structural, and dynamical properties for the 1-butyl-3-methylimidazolium trifluoroacetate ([C4mim][CF3COO]) and 1-butyl-3-methylimidazolium bromide, ([C4mim][Br]) ionic liquids (ILs) were calculated using Non-polarizable Force Fields (FF), parameterized using a methodology developed previously within the research group, for condensed phase applications. Properties such as the Vapor-Liquid Equilibrium (VLE) curve, critical points (ρ c, T c), Radial, Spatial and Combined Distribution Functions and self-diffusion coefficients were calculated using Equilibrium Molecular Dynamics simulations (EMD); other properties such as shear viscosities and thermal conductivities were calculated using Non-Equilibrium Molecular Dynamics simulations (NEMD). The results obtained in this work indicated that the calculated critical points are comparable with those available in the literature. The calculated structural information for these two ILs indicated that the anions interact mainly with hydrogen atoms from both the imidazolium ring and the methyl chain; the bromide anion displays twice the hydrogen coordination number than the oxygen atoms from the trifluoroacetate anion. Furthermore, Non-Covalent interactions (NCI index), determined by DFT calculations, revealed that some hydrogen bonds in the [C4mim][Br] IL displayed similar strength to those in the [C4mim][CF3COO] IL, in spite of the shorter O--H distances found in the latter IL. The majority of the calculated transport properties presented reasonable agreement with the experimental available data. Nonetheless, the self-diffusion coefficients determined in this work are under-estimated with respect to experimental values; however, by escalating the electrostatic atomic charges for the anion and cation to ±0.8e, only for this property, a remarkable improvement was obtained. Experimental evidence was recovered for most of the calculated properties and to the best of our knowledge, some new predictions were done mainly in thermodynamic states where data are not available. To validate the FF, developed previously within the research group, dynamic properties were also evaluated for a series of ILs such as [C4mim][PF6], [C4mim][BF4], [C4mim][OMs], and [C4mim][NTf2] ILs.

Scalable Screening of Soft Matter: A Case Study of Mixtures of Ionic Liquids and Organic Solvents

Room-temperature ionic liquids (RTILs) are a class of organic salts that are liquid at room temperature. Their physiochemical properties, including low vapor pressure and wide electrochemical stability window, have driven their use as electrolytes in many electrochemical applications; however, the slow transport properties of many RTILs have limited their utility in some applications. This issue is often mitigated by solvating ionic liquids in neutral organic solvents. To date, however, solvent interactions have only been explored for a small number of solvents, particularly acetonitrile and propylene carbonate, at only a few compositions. In this work, we use molecular dynamics simulations in the context of a computational screening approach to study mixtures of ionic liquids in many different solvents at a range of concentrations. Building on prior work, we again find that ionic liquid diffusivity increases monotonically with greater solvent concentration. In contrast to prior work, we find that pure solvent diffusivity, not polarity, is the most influential solvent property on mixture behavior. We also explore the concentration dependence of ionic conductivity and find maxima at intermediate concentrations. Experimental conductivity measurements, inspired by the computational screening study, support this observation with qualitatively consistent results. These results can further guide the selection of solvents for electrochemical applications of RTILs.

Insights into the influence of the molecular structures of fluorinated ionic liquids on their thermophysical properties. A soft-SAFT based approach

Development of predictive models for FILs.

Insights into High Conductivity of the Two-Dimensional Iodine-Oxidized sp2-c-COF

A recent experiment [ Jin , E. ; Science 2017 , 357 , 673 - 676 ] shows that the conductivity of a two-dimensional (2D) sp²-carbon-hybridized π-conjugated covalent organic framework (sp²-c-COF) can be enhanced by as much as 12 orders of magnitude after iodine oxidation processing. To understand the physical mechanism underlying such a huge increase in the conductivity, we perform multiscale computations and find that the high conductivity of the iodine-oxidized 2D COF can be attributed to both hole transfer and ion transfer within the 2D COF. The computed dominant charge distribution corresponding to the valence band maximum (VBM) suggests that the delocalized π electrons occur mostly at the active reaction sites. The computed low ionization energy at the active reaction sites further supports that the 2D COF tends to lose electrons during iodine oxidation and to yield cationic COF and anionic triiodide I₃^-. Complementary classical molecular dynamics simulation shows a relatively high anion conductivity of 13.63 × 10^-2 S m^-1, consistent with the high conductivity measured from the experiment (7.1 × 10^-2 S m^-1). Meanwhile, we find that the cations in 2D COF can also induce a shift of the Fermi level to cross the valence band, thereby enhancing the hole mobility to 86.75 cm² V^-1 s^-1. For proposing a potential application of the highly conductive iodine-oxidized 2D sp²-c-COF, we design a prototypical model of the 2D spirally wound lithium-ion battery and find that it exhibits enhanced stability than a typical electrolyte material.

Rheology of phosphonium ionic liquids: a molecular dynamics and experimental study

Comparison between the theoretical and experimental viscosity of an ionic liquid.