Maxwell–Stefan Diffusivities in Binary Mixtures of Ionic Liquids with Dimethyl Sulfoxide (DMSO) and H2O

Influence of alkali metals on water dynamics inside imidazolium-based ionic liquid nano-domains

The global need to expand the design of energy-storage devices led to the investigation of alkali metal - Ionic Liquid (IL) mixtures as a possible class of electrolytes. In this study, 1D and 2D Nuclear Magnetic Resonance (NMR) and Electrochemical Impedance Spectroscopy (EIS) as well as Molecular Dynamics (MD) simulations were used to study the intermolecular interactions in imidazolium-based IL - water - alkali halide ternary mixtures. The 1H and 23Na 1D and 1H DOSY NMR spectra revealed that the presence of small quantities of NaCl does not influence the aggregation of water molecules in the IL nano-domains. The order of adding ionic compounds to water, as well as the certain water and NaCl molecular ratios, lead to the formation of isolated water clusters. Two ternary solutions representing different orders of compounds mixing (H2O+ IL + NaCl or H2O+ NaCl + IL) showed a strong dependence of the initial solvation shell of Na+ and the self-clustering of water. Furthermore, the behaviour of water was found to be independent from the conditions applied during the solution preparation, such as temperature and/or duration of stirring and aging. These findings could be confirmed by large differences in the amount of ionic species, observed in the ternary solutions and depending on the order of mixing/solute preparation.

Fick Diffusion Coefficient in Binary Mixtures of [HMIM][NTf2] and Carbon Dioxide by Dynamic Light Scattering and Molecular Dynamics Simulations

Dynamic light scattering (DLS) experiments and equilibrium molecular dynamics (EMD) simulations were performed in the saturated liquid phase of the binary mixture of 1-hexyl-3-methylimidazolium bis(trifluormethylsulfonyl)imide ([HMIM][NTf₂]) and carbon dioxide (CO₂) to access the Fick diffusion coefficient (D₁₁). The investigations were performed within or close to saturation conditions at temperatures between (298.15 and 348.15) K and CO₂ mole fractions (x_CO2) up to 0.81. The DLS experiments were combined with polarization-difference Raman spectroscopy (PDRS) to simultaneously access the composition of the liquid phase. For the first time in an electrolyte-based system, D₁₁ was directly calculated from EMD simulations by accessing the Maxwell-Stefan (MS) diffusion coefficient and the thermodynamic factor. Agreement within combined uncertainties was found between D₁₁ from DLS and EMD simulations for CO₂ mole fractions up to 0.5. In general, an increasing D₁₁ with increasing x_CO2 could be observed, with a local maximum present at a CO₂ mole fraction of about 0.75. The local maximum could be explained by an increasing MS diffusion coefficient with increasing x_CO2 over the entire studied composition range and a decreasing thermodynamic factor at x_CO2 above 0.7. Finally, PDRS and EMD simulations were combined to investigate the influence of the fluid structure on the diffusive process.

Machine Learning-Based Upscaling of Finite-Size Molecular Dynamics Diffusion Simulations for Binary Fluids

Molecular diffusion coefficients calculated using molecular dynamics (MD) simulations suffer from finite-size (i.e., finite box size and finite particle number) effects. Results from finite-sized MD simulations can be upscaled to infinite simulation size by applying a correction factor. For self-diffusion of single-component fluids, this correction has been well-studied by many researchers including Yeh and Hummer (YH); for binary fluid mixtures, a modified YH correction was recently proposed for correcting MD-predicted Maxwell-Stephan (MS) diffusion rates. Here we use both empirical and machine learning methods to identify improvements to the finite-size correction factors for both self-diffusion and MS diffusion of binary Lennard-Jones (LJ) fluid mixtures. Using artificial neural networks (ANNs), the error in the corrected LJ fluid diffusion is reduced by an order of magnitude versus existing YH corrections, and the ANN models perform well for mixtures with large dissimilarities in size and interaction energies where the YH correction proves insufficient.

Properties of Ion Complexes and Their Impact on Charge Transport in Organic Solvent-Based Electrolyte Solutions for Lithium Batteries: Insights from a Theoretical Perspective

Electrolyte formulations in standard lithium ion and lithium metal batteries are complex mixtures of various components. In this article, we review molecular key principles of ion complexes in multicomponent electrolyte solutions in regards of their influence on charge transport mechanisms. We outline basic concepts for the description of ion–solvent and ion–ion interactions, which can be used to rationalize recent experimental and numerical findings concerning modern electrolyte formulations. Furthermore, we discuss benefits and drawbacks of empirical concepts in comparison to molecular theories of solution for a more refined understanding of ion behavior in organic solvents. The outcomes of our discussion provide a rational for beneficial properties of ions, solvent, co-solvent and additive molecules, and highlight possible routes for further improvement of novel electrolyte solutions.

Shear Viscosity Computed from the Finite-Size Effects of Self-Diffusivity in Equilibrium Molecular Dynamics

A method is proposed for calculating the shear viscosity of a liquid from finite-size effects of self-diffusion coefficients in Molecular Dynamics simulations. This method uses the difference in the self-diffusivities, computed from at least two system sizes, and an analytic equation to calculate the shear viscosity. To enable the efficient use of this method, a set of guidelines is developed. The most efficient number of system sizes is two and the large system is at least four times the small system. The number of independent simulations for each system size should be assigned in such a way that 50%–70% of the total available computational resources are allocated to the large system. We verified the method for 250 binary and 26 ternary Lennard-Jones systems, pure water, and an ionic liquid ([Bmim][Tf₂N]). The computed shear viscosities are in good agreement with viscosities obtained from equilibrium Molecular Dynamics simulations for all liquid systems far from the critical point. Our results indicate that the proposed method is suitable for multicomponent mixtures and highly viscous liquids. This may enable the systematic screening of the viscosities of ionic liquids and deep eutectic solvents.

Finite-Size Effects of Binary Mutual Diffusion Coefficients from Molecular Dynamics

Molecular dynamics simulations were performed for the prediction of the finite-size effects of Maxwell-Stefan diffusion coefficients of molecular mixtures and a wide variety of binary Lennard-Jones systems. A strong dependency of computed diffusivities on the system size was observed. Computed diffusivities were found to increase with the number of molecules. We propose a correction for the extrapolation of Maxwell-Stefan diffusion coefficients to the thermodynamic limit, based on the study by Yeh and Hummer ( J. Phys. Chem. B , 2004 , 108 , 15873 - 15879 ). The proposed correction is a function of the viscosity of the system, the size of the simulation box, and the thermodynamic factor, which is a measure for the nonideality of the mixture. Verification is carried out for more than 200 distinct binary Lennard-Jones systems, as well as 9 binary systems of methanol, water, ethanol, acetone, methylamine, and carbon tetrachloride. Significant deviations between finite-size Maxwell-Stefan diffusivities and the corresponding diffusivities at the thermodynamic limit were found for mixtures close to demixing. In these cases, the finite-size correction can be even larger than the simulated (finite-size) Maxwell-Stefan diffusivity. Our results show that considering these finite-size effects is crucial and that the suggested correction allows for reliable computations.

Evaluation of thermophysical properties of ionic liquids with polar solvent: a comparable study of two families of ionic liquids with various ions

In this work, we explore and compare the role of the ion effect on the thermophysical properties of two families of ionic liquids (ILs), namely, tetra-alkyl ammonium cation [R4N](+) with hydroxide [OH](-) anion and 1-alkyl-3-methyl imidazolium cation [amim](+) with different anions (chloride, methyl sulfate, and tetrafluoroborate), with polar solvent such as dimethylsulfoxide (DMSO) in the temperature range from 25 to 40 °C and over the whole concentration range of ILs. Two families of ILs, namely, tetramethyl ammonium hydroxide [(CH3)4N][OH] (TMAH), tetraethyl ammonium hydroxide [(C2H5)4N][OH] (TEAH), tetrapropyl ammonium hydroxide [(C3H7)4N][OH] (TPAH), and tetrabutyl ammonium hydroxide [(C4H9)4N][OH] (TBAH) from ammonium-based ILs and 1-ethyl-3-methylimidazolium chloride [Emim][Cl], 1-ethyl-3-methylimidazolium methylsulfate [Emim][MeSO4], 1-butyl-3-methylimidazolium tetrafluoroborate [Bmim][BF4], and 1-butyl-3-methylimidazolium chloride ([Bmim][Cl]) from imidazolium family of ILs, are used in the present study. To address the molecular interactions of ILs with DMSO, densities (ρ), ultrasonic sound velocities (u), and viscosities (η) have been measured over the entire composition range and at four temperatures, 25, 30, 35, and 40 °C, under atmospheric pressure. From these experimental data, the excess molar volume (V(E)), the deviation in isentropic compressibility (Δκs), and the deviation in viscosity (Δη) were calculated and were adequately correlated by using the Redlich-Kister polynomial equation. The measured and predicted data were interpreted on the basis of intermolecular interactions and structural effects between like and unlike molecules upon mixing. The hydrogen-bonding features between ammonium-based ILs and DMSO were analyzed using molecular modeling program by HyperChem 7.