Do Electrostatics Control the Diffusive Dynamics of Solitary Water? NMR and MD Studies of Water Translation and Rotation in Dipolar and Ionic Solvents

NMR-based measurements of the diffusion coefficients and rotation times of solitary water and benzene at 300 K are reported in a diverse collection of 13 conventional organic solvents and 10 imidazolium ionic liquids. Proton chemical shifts of water are found to be correlated to water OH-stretching frequencies, confirming the importance of electrostatic interactions in these shifts. However, the influence of magnetic interactions in aromatic solvents renders chemical shifts a less reliable indicator of electrostatics. Diffusion coefficients (DB) and rotational correlation times (τB) of benzene in the solvents examined are accurately described as functions of viscosity (η) by DB ∝ η-0.81 and τB ∝ η0.64. Literature values of DB and τB in alkane and normal alcohols, which were not included among the solvents studied here, are systematically faster than predicted by these correlations, indicating that factors beyond solvent viscosity play a role in determining the friction on benzene. In contrast to benzene, water diffusion and rotation are poorly described in terms of viscosity alone, even in the dipolar and ionic solvents measured here. The present data and the substantial literature data already available on dilute water diffusion show a systematic dependence of DW on solvent polarity among isoviscous solvents. The aspect of solvent polarity most relevant to water dynamics is the ability of a solvent to accept hydrogen bonds from water, as conveniently quantified by the frequency of water's OH stretching band, ΔνOH. The friction on translation, ζtr = kBT/DW, and rotation, ζrot = kBTτW, are both well correlated by functions of the form ζ(η, ΔνOH) = a1ηa2 exp (a3ΔνOH), where the ai are adjustable parameters. Molecular dynamics simulations reveal a strong coupling between electrostatic and nonelectrostatic water-solvent interactions, which makes it impossible to dissect the friction on water into additive dielectric and hydrodynamic components. Simulations also provide a tentative explanation for the unusual form of the correlating function ζ(η, ΔνOH), at least in the case of ζrot.

Influence of Residual Water Traces on the Electrochemical Performance of Hydrophobic Ionic Liquids for Magnesium-Containing Electrolytes

A trace amount of water is typically unavoidable as an impurity in ionic liquids, which is a huge challenge for their application in Mg-ion batteries. Here, we employed molecular sieves of different pore diameters (3, 4, and 5 Å), to effectively remove the trace amounts of water from 1-methyl-1-propylpiperidinium bis(trifluoromethylsulfonyl)imide (MPPip-TFSI) and 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide (BMP-TFSI). Notably, after sieving (water content <1 mg ⋅ L-1 ), new anodic peaks arise that are attributed to the formation of different anion-cation structures induced by minimizing the influence of hydrogen bonds. Furthermore, electrochemical impedance spectroscopy (EIS) reveals that the electrolyte resistance decreases by ∼10 % for MPPip-TFSI and by ∼28 % for BMP-TFSI after sieving. The electrochemical Mg deposition/dissolution is investigated in MPPip-TFSI/tetraglyme (1 : 1)+100 mM Mg(TFSI)2 +10 mM Mg(BH4 )2 using Ag/AgCl and Mg reference electrodes. The presence of a trace amount of water leads to a considerable shift of 0.9 V vs. Mg2+/ Mg in the overpotential of Mg deposition. In contrast, drying of MPPip-TFSI enhances the reversibility of Mg deposition/dissolution and suppresses the passivation of the Mg electrode.

MDSuite: comprehensive post-processing tool for particle simulations

Particle-Based (PB) simulations, including Molecular Dynamics (MD), provide access to system observables that are not easily available experimentally. However, in most cases, PB data needs to be processed after a simulation to extract these observables. One of the main challenges in post-processing PB simulations is managing the large amounts of data typically generated without incurring memory or computational capacity limitations. In this work, we introduce the post-processing tool: MDSuite. This software, developed in Python, combines state-of-the-art computing technologies such as TensorFlow, with modern data management tools such as HDF5 and SQL for a fast, scalable, and accurate PB data processing engine. This package, built around the principles of FAIR data, provides a memory safe, parallelized, and GPU accelerated environment for the analysis of particle simulations. The software currently offers 17 calculators for the computation of properties including diffusion coefficients, thermal conductivity, viscosity, radial distribution functions, coordination numbers, and more. Further, the object-oriented framework allows for the rapid implementation of new calculators or file-readers for different simulation software. The Python front-end provides a familiar interface for many users in the scientific community and a mild learning curve for the inexperienced. Future developments will include the introduction of more analysis associated with ab-initio methods, colloidal/macroscopic particle methods, and extension to experimental data.

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.

Specific Ion Effects in Different Media: Current Status and Future Challenges

We discuss the current state of research as well as the future challenges for a deeper understanding of specific ion effects in protic and aprotic solvents as well as various additional media. Despite recent interest in solute or interfacial effects, we focus exclusively on the specific properties of ions in bulk electrolyte solutions. Corresponding results show that many mechanisms remain unknown for these simple media, although theoretical, computational, and experimental studies have provided some insights into explaining individual observations. In particular, the importance of local interactions and electronic properties is emphasized, which enabled a more consistent interpretation of specific ion effects over the past years. Despite current insufficient knowledge, we also discuss future challenges in relation to dynamic properties as well as the influence of different concentrations, different solvents, and solute contributions to gain a deeper understanding of specific ion effects for technological applications.

Cyclodextrin solubilization in hydrated reline: Resolving the unique stabilization mechanism in a deep eutectic solvent

By complexing with hydrophobic compounds, cyclodextrins afford increased solubility and thermodynamic stability to hardly soluble compounds, thereby underlining their invaluable applications in pharmaceutical and other industries. However, common cyclodextrins such as β-cyclodextrin, suffer from limited solubility in water, which often leads to precipitation and formation of unfavorable aggregates, driving the search for better solvents. Here, we study the solvation of cyclodextrin in deep eutectic solvents (DESs), environmentally friendly media that possess unique properties. We focus on reline, the DES formed from choline chloride and urea, and resolve the mechanism through which its constituents elevate β-cyclodextrin solubility in hydrated solutions compared to pure water or dry reline. Combining experiments and simulations, we determine that the remarkable solubilization of β-cyclodextrin in hydrated reline is mostly due to the inclusion of urea inside β-cyclodextrin's cavity and at its exterior surfaces. The role of choline chloride in further increasing solvation is twofold. First, it increases urea's solubility beyond the saturation limit in water, ultimately leading to much higher β-cyclodextrin solubility in hydrated reline in comparison to aqueous urea solutions. Second, choline chloride increases urea's accumulation in β-cyclodextrin's vicinity. Specifically, we find that the accumulation of urea becomes stronger at high reline concentrations, as the solution transitions from reline-in-water to water-in-reline, where water alone cannot be regarded as the solvent. Simulations further suggest that in dry DES, the mechanism of β-cyclodextrin solvation changes so that reline acts as a quasi-single component solvent that lacks preference for the accumulation of urea or choline chloride around β-cyclodextrin.

How sensitive are physical properties of choline chloride-urea mixtures to composition changes: Molecular dynamics simulations and Kirkwood-Buff theory

Deep eutectic solvents (DESs) have emerged as a cheaper and greener alternative to conventional organic solvents. Choline chloride (ChCl) mixed with urea at a molar ratio of 1:2 is one of the most common DESs for a wide range of applications such as electrochemistry, material science, and biochemistry. In this study, molecular dynamics simulations are performed to study the effect of urea content on the thermodynamic and transport properties of ChCl and urea mixtures. With increased mole fraction of urea, the number of hydrogen bonds (HBs) between cation-anion and ion-urea decreases, while the number of HBs between urea-urea increases. Radial distribution functions (RDFs) for ChCl-urea and ChCl-ChCl pairs shows a significant decrease as the mole fraction of urea increases. Using the computed RDFs, Kirkwood-Buff Integrals (KBIs) are computed. KBIs show that interactions of urea-urea become stronger, while interactions of urea-ChCl and ChCl-ChCl pairs become slightly weaker with increasing mole fraction of urea. All thermodynamic factors are found larger than one, indicating a non-ideal mixture. Our results also show that self- and collective diffusivities increase, while viscosities decrease with increasing urea content. This is mainly due to the weaker interactions between ions and urea, resulting in enhanced mobilities. Ionic conductivities exhibit a non-monotonic behavior. Up to a mole fraction of 0.5, the ionic conductivities increase with increasing urea content and then reach a plateau.

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.

Use of Ionic Liquids in Protein and DNA Chemistry

Ionic liquids (ILs) have been receiving much attention as solvents in various areas of biochemistry because of their various beneficial properties over the volatile solvents and ILs availability in myriad variants (perhaps as many as 108) owing to the possibility of paring one cation with several anions and vice-versa as well as formulations as zwitterions. Their potential as solvents lies in their tendency to offer both directional and non-directional forces toward a solute molecule. Because of these forces, ionic liquids easily undergo intermolecular interactions with a range of polar/non-polar solutes, including biomolecules such as proteins and DNA. The interaction of genomic species in aqueous/non-aqueous states assists in unraveling their structure and functioning, which have implications in various biomedical applications. The charge density of ionic liquids renders them hydrophilic and hydrophobic, which retain intact over long-range of temperatures. Their ability in stabilizing or destabilizing the 3D-structure of a protein or the double-helical structure of DNA has been assessed superior to the water and volatile organic solvents. The aptitude of an ion in influencing the structure and stability of a native protein depends on their ranking in the Hofmeister series. However, at several instances, a reverse Hofmeister ordering of ions and specific ion-solute interaction has been observed. The capability of an ionic liquid in terms of the tendency to promote the coiling/uncoiling of DNA structure is noted to rely on the basicity, electrostatic interaction, and hydrophobicity of the ionic liquid in question. Any change in the DNA's double-helical structure reflects a change in its melting temperature (T m), compared to a standard buffer solution. These changes in DNA structure have implications in biosensor design and targeted drug-delivery in biomedical applications. In the current review, we have attempted to highlight various aspects of ionic liquids that influence the structure and properties of proteins and DNA. In short, the review will address the issues related to the origin and strength of intermolecular interactions, the effect of structural components, their nature, and the influence of temperature, pH, and additives on them.

Microstructural and Dynamical Heterogeneities in Ionic Liquids

Ionic liquids (ILs) are a special category of molten salts solely composed of ions with varied molecular symmetry and charge delocalization. The versatility in combining varied cation-anion moieties and in functionalizing ions with different atoms and molecular groups contributes to their peculiar interactions ranging from weak isotropic associations to strong, specific, and anisotropic forces. A delicate interplay among intra- and intermolecular interactions facilitates the formation of heterogeneous microstructures and liquid morphologies, which further contributes to their striking dynamical properties. Microstructural and dynamical heterogeneities of ILs lead to their multifaceted properties described by an inherent designer feature, which makes ILs important candidates for novel solvents, electrolytes, and functional materials in academia and industrial applications. Due to a massive number of combinations of ion pairs with ion species having distinct molecular structures and IL mixtures containing varied molecular solvents, a comprehensive understanding of their hierarchical structural and dynamical quantities is of great significance for a rational selection of ILs with appropriate properties and thereafter advancing their macroscopic functionalities in applications. In this review, we comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.

Structural investigations of molecular solutes within nanostructured ionic liquids

Ionic liquids (ILs) containing sufficiently long alkyl chains form amphiphilic nanostructures with well-defined polar and non-polar domains. Here we have explored the robustness of these amphiphilic nanostructures to added solutes and gained insight into how the nature of the solute and IL ions affect the partitioning of these solutes within the nanostructured domains of ILs. To achieve this, small angle X-ray scattering (SAXS) investigations were performed and discussed for mixtures of 9 different molecular compounds with 6 different ILs containing imidazolium cations. The amphiphilic nanostructure of ILs persisted to high solute concentrations, over 50 mol% of added solute for most 1-butyl-3-methylimidazolium ILs and above 80 mol% for most 1-decyl-3-methylimidazolium ILs. Solute partitioning within these domains was found to be controlled by the inherent polarity and size of the solute, as well as specific interactions between the solute and IL ions, with SAXS results corroborated with IR spectroscopy and molecular dynamics simulations. Molecular dynamics simulations also revealed the ability to induce π⁺-π⁺ stacking between imidazolium cations through the use of these added molecular compounds. Collectively, these results provide scope for the selection of IL ions to rationally influence and control the partitioning behaviour of given solutes within the amphiphilic nanostructure of ILs.

Thermodynamic Analysis of n-Hexane-Ethanol Binary Mixtures Using the Kirkwood-Buff Theory

A complete thermodynamic analysis of mixtures consisting of molecules with complex chemical constitution can be rather demanding. The Kirkwood-Buff theory of solutions allows the estimation of thermodynamic properties, which cannot be directly extracted from atomistic simulations, such as the Gibbs energy of mixing (Δ_mix G). In this work, we perform molecular dynamics simulations of n-hexane-ethanol binary mixtures in the liquid state under two temperature-pressure conditions and at various mole fractions. On the basis of the recently published methodology of Galata [ Fluid Phase Equilib. 2018 , 470 , 25 - 37 ] , we first calculate the Kirkwood-Buff integrals in the isothermal-isobaric ( NpT) ensemble, identifying how system size affects their estimation. We then extract the activity coefficients, excess Gibbs energy, excess enthalpy, and excess entropy for the n-hexane-ethanol binary mixtures we simulate. We employ two approaches for quantifying composition fluctuations: one based on counting molecular centers of mass and a second one based on counting molecular segments. Results from the two approaches are practically indistinguishable. We compare our results against predictions of vapor-liquid equilibria obtained in a previous simulation work using the same force field, as well as with experimental data, and find very good agreement. In addition, we develop a simple methodology to identify the hydrogen bonds between ethanol molecules and analyze their effects on mixing properties.

3D structure of the electric double layer of ionic liquid-alcohol mixtures at the electrochemical interface

Mixtures of the ionic liquid 1-butyl-3-methylimidazolium tetrafluoroborate with amphiphilic cosolvents, such as methanol and ethanol, nanoconfined between graphene walls are studied by means of molecular dynamics simulations and the results are compared with those of the pure ionic liquid and its mixtures with water confined in the same conditions. We investigate the adsorption of cosolvent molecules at the graphene walls as well as their distribution across the system. The results show that, due to a higher affinity of the polar groups to be close to the anions in combination with the electrostatic and excluded volume interactions, there exists a high tendency of the OH groups to lie close to the anode, inducing small changes in the first cation layer. The orientation of cosolvent molecules is found to be closely related to the alignment of the molecular dipole moment. We also investigate the lateral ionic distribution in the layers close to the electrodes, which shows a structural transition from liquid-like lamellar ordering to solid-like hexagonal patterns as the size of the cosolvent molecules increases leading to smaller position fluctuations of the ions. The dependence of the specific patterns on the nature of the electrodes is also studied. This study strongly suggests that the ionic patterns formed in the first ionic layers next to the charged interfaces are universal since their existence does not crucially depend on the atomic composition of the interfacial material, but only on the net charge density of the considered ionic layer, which significantly changes the ionic mobility in this region.

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.

Aqueous ionic liquids in comparison with standard co-solutes : Differences and common principles in their interaction with protein and DNA structures

Ionic liquids (ILs) are versatile solvents for a broad range of biotechnological applications. Recent experimental and simulation results highlight the potential benefits of dilute ILs in aqueous solution (aqueous ILs) in order to modify protein and DNA structures systematically. In contrast to a limited number of standard co-solutes like urea, ectoine, trimethylamine-N-oxide (TMAO), or guanidinium chloride, the large amount of possible cation and anion combinations in aqueous ILs can be used to develop tailor-made stabilizers or destabilizers for specific purposes. In this review article, we highlight common principles and differences between aqueous ILs and standard co-solutes with a specific focus on their underlying macromolecular stabilization or destabilization behavior. In combination with statistical thermodynamics theories, we present an efficient framework, which is used to classify structure modification effects consistently. The crucial importance of enthalpic and entropic contributions to the free energy change upon IL-assisted macromolecular unfolding in combination with a complex destabilization mechanism is described in detail. A special focus is also set on aqueous IL-DNA interactions, for which experimental and simulation outcomes are summarized and discussed in the context of previous findings.