Theoretical-Computational Modeling of Gas-State Thermodynamics in Flexible Molecular Systems: Ionic Liquids in the Gas Phase as a Case Study

A theoretical-computational procedure based on the quasi-Gaussian entropy (QGE) theory and molecular dynamics (MD) simulations is proposed for the calculation of thermodynamic properties for molecular and supra-molecular species in the gas phase. The peculiarity of the methodology reported in this study is its ability to construct an analytical model of all the most relevant thermodynamic properties, even within a wide temperature range, based on a practically automatic sampling of the entire conformational repertoire of highly flexible systems, thereby bypassing the need for an explicit search for all possible conformers/rotamers deemed relevant. In this respect, the reliability of the presented method mainly depends on the quality of the force field used in the MD simulations and on the ability to discriminate in a physically coherent way between semi-classical and quantum degrees of freedom. The method was tested on six model systems (n-butane, n-butane, n-octanol, octadecane, 1-butyl-3-methylimidazolium hexafluorophosphate and 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic pairs), which, being experimentally characterized and already addressed by other theoretical-computational methods, were considered as particularly suitable to allow us to evaluate the method’s accuracy and efficiency, bringing out advantages and possible drawbacks. The results demonstrate that such a physically coherent yet relatively simple method can represent a further valid computational tool that is alternative and complementary to other extremely efficient computational methods, as it is particularly suited for addressing the thermodynamics of gaseous systems with a high conformational complexity over a large range of temperature.     

Aggregation Behavior of Trisiloxane-Tailed Surface Active Ionic Liquids in Aqueous Solution: Coarse-Grained Molecular Dynamics Study

The aggregation behaviors of two trisiloxane-tailed surface active ionic liquids in water have been investigated by coarse-grained (CG) molecular dynamics simulation on the basis of MARTINI force field. The new CG model is developed from the optimized molecule computed by using density functional theory. Direct comparison of angles and bonds obtained from all-atom (AA) simulations with those calculated from the CG model has been conducted to validate the latter model. Excellent agreement between AA and CG demonstrates that the potential of the new CG model can represent the complex system well. The long time CG simulation has been performed to understand the formation process of micelles when dissolving ionic liquids in water. Vesicles were observed at the final stage of the simulation and their partially truncated views and density profiles were obtained to describe the structure in detail.     

Effects of Cationic Structure on Cellulose Dissolution in Ionic Liquids: A Molecular Dynamics Study

In recent years, great progress has been made in the dissolution of cellulose with ionic liquids (ILs). However, the mechanism of cellulose dissolution, especially the role the IL cation played in the dissolution process, has not been clearly understood. Herein, the mixtures of cellulose with a series of imidazolium‐based chloride ionic liquids and 1‐butyl‐3‐methyl pyridinium chloride ([C₄mpy]Cl) were simulated to study the effect that varying the heterocyclic structure and alkyl chain length of the IL cation has on the dissolution of cellulose. It was shown that the dissolution of cellulose in [C₄mpy]Cl is better than that in [C₄mim]Cl. For imidazolium‐based ILs, the shorter the alkyl chain is, the higher the solubility will be. In addition, an all‐atom force field for 1‐allyl‐3‐methyl imidazolium cation ([Amim]⁺) was developed, for the first time, to investigate the effect the electron‐withdrawing group within the alkyl chain of the IL cation has on the dissolution of cellulose. It was found that the interaction energy between [Amim]⁺ and cellulose was greater than that between [C₃mim]⁺ and cellulose, indicating that the presence of electron‐withdrawing group in alkyl chain of the cation enhanced the interaction between the cation and cellulose due to the increase of electronegativity of the cations. These findings are used to assess the cationic effect on the dissolution of cellulose in ILs. They are also expected to be important for rational design of novel ILs for efficient dissolution of cellulose.     

Broadband NMR Relaxometry as a Powerful Technique to Study Molecular Dynamics of Ionic Liquids

Fast field cycling nuclear magnetic resonance (FFC NMR) relaxometry technique has been demonstrated to be a useful analytical tool to investigate molecular dynamics in very diverse systems during the last decades. Of particular importance has been its application in studying ionic liquids, upon which this review article is based. Some of the research carried out on ionic liquids during the last ten years using this technique is highlighted in this article with the aim of promoting the favorable features of FFC NMR applied toward understanding dynamics of complex systems.     

What Far‐Infrared Spectra Can Contribute to the Development of Force Fields for Ionic Liquids Used in Molecular Dynamics Simulations

Symbiosis: Far‐infrared spectra can be used to check the quality of force fields for molecular dynamics simulations of ionic liquids. On the other hand, MD simulations can explain the molecular basis of measured properties for this new liquid material (see picture).

     

Role of Hydrogen Bond Defects for Cluster Formation and Distribution in Ionic Liquids by Means of Neutron Diffraction and Molecular Dynamics Simulations

Defects fundamentally govern the properties of all real materials. Correlating molecular defects to macroscopic quantities remains a challenge, particularly in the liquid phase. Herein, we report the influence of hydrogen bonds (HB) acting as defects in mixtures of non-hydroxyl-functionalized ionic liquids (ILs) with an increasing concentration of hydroxyl-functionalized ILs. We observed two types of HB defects: The conventional HBs between cation and anion (c–a), and the elusive HBs between cations (c–c) despite the repulsive Coulomb forces. We use neutron diffraction with isotopic substitution in combination with molecular dynamics simulations for measuring the geometry, strength, and distribution of mobile OH defects in the IL mixtures. In principle, this procedure allows relating the number and stability of defects to macroscopic properties such as diffusion, viscosity, and conductivity, which are of utmost importance for the performance of electrolytes in batteries and other electrical devices.     

Ether-Functionalized Pyrrolidinium-Based Room Temperature Ionic Liquids: Physicochemical Properties,Molecular Dynamics,and the Lithium Ion Coordination Environment

The physicochemical properties of room temperature ionic liquids (RTILs) consisting of bis(trifluoromethanesulfonyl)amide (TFSA⁻) combined with 1-hexyl-1-methylpyrrolidinium (Pyr_1,6⁺), 1-(butoxymethyl)-1-methylpyrrolidinium (Pyr_1,1O4⁺), 1-(4-methoxybutyl)-1-methyl pyrrolidinium (Pyr_1,4O1⁺), and 1-((2-methoxyethoxy)methyl)-1-methylpyrrolidinium (Pyr_1,1O2O1⁺) were investigated using both experimental and computational approaches. Pyr_1,1O2O1TFSA, which contains two ether oxygen atoms, showed the lowest viscosity, and the relationship between its physicochemical properties and the position and number of the ether oxygen atoms was discussed by a careful comparison with Pyr_1,1O4TFSA and Pyr_1,4O1TFSA. Ab initio calculations revealed the conformational flexibility of the side chain containing the ether oxygen atoms. In addition, molecular dynamics (MD) calculations suggested that the ion distributions have a significant impact on the transport properties. Furthermore, the coordination environments of the Li ions in the RTILs were evaluated using Raman spectroscopy, which was supported by MD calculations using 1000 ion pairs. The presented results will be valuable for the design of functionalized RTILs for various applications.     

Ionic Liquids: Dissecting the Enthalpies of Vaporization

We calculate the heats of vaporisation for imidazolium‐based ionic liquids [C_nmim][NTf₂] with n=1, 2, 4, 6, 8 by means of molecular dynamics (MD) simulations and discuss their behavior with respect to temperature and the alkyl chain length. We use a force field developed recently. The different cohesive energies contributing to the overall heats of vaporisations are discussed in detail. With increasing alkyl chain length, the Coulomb contribution to the heat of vaporisation remains constant at around 80 kJ mol^?1, whereas the van der Waals interaction increases continuously. The calculated increase of about 4.7 kJ mol^?1 per CH₂‐group of the van der Waals contribution in the ionic liquid exactly coincides with the increase in the heats of vaporisation for n‐alcohols and n‐alkanes, respectively. The results support the importance of van der Waals interactions even in systems completely composed of ions.     

The Double‐Faced Nature of Hydrogen Bonding in Hydroxy‐Functionalized Ionic Liquids Shown by Neutron Diffraction and Molecular Dynamics Simulations

We characterize the double‐faced nature of hydrogen bonding in hydroxy‐functionalized ionic liquids by means of neutron diffraction with isotopic substitution (NDIS), molecular dynamics (MD) simulations, and quantum chemical calculations. NDIS data are fit using the empirical potential structure refinement technique (EPSR) to elucidate the nearest neighbor H???O and O???O pair distribution functions for hydrogen bonds between ions of opposite charge and the same charge. Despite the presence of repulsive Coulomb forces, the cation–cation interaction is stronger than the cation–anion interaction. We compare the hydrogen‐bond geometries of both “doubly charged hydrogen bonds” with those reported for molecular liquids, such as water and alcohols. In combination, the NDIS measurements and MD simulations reveal the subtle balance between the two types of hydrogen bonds: The small transition enthalpy suggests that the elusive like‐charge attraction is almost competitive with conventional ion‐pair formation.     

Effects of Ionic Liquids on the Thermodynamics of Hydrogen Activation by Frustrated Lewis Pairs: A Density Functional Theory Study**

Nowadays, hydrogen activation by frustrated Lewis pairs (FLPs) and their applications are one of the emerging research topics in the field of catalysis. Previous studies have shown that the thermodynamics of this reaction is determined by electronic structures of FLPs and solvents. Herein, we investigated systems consisting of typical FLPs and ionic liquids (ILs), which are well known by their large number of types and excellent solvent effects. The density functional theory (DFT) calculations were performed to study the thermodynamics for H₂ activation by both inter- and intra-molecular FLPs, as well as the individual components. The results show that the computed overall Gibbs free energies in ILs are more negative than that computed in toluene. Through the thermodynamics partitioning, we find that ILs favor the H−H cleavage elemental step over the elemental steps of proton attachment, hydride attachment and zwitterionic stabilization. Moreover, the results show that these effects are strongly dependent on the type of FLPs, where intra-molecular FLPs are more affected compared to the inter-molecular FLPs.     

TiCl4 Dissolved in Ionic Liquid Mixtures from Аb Initio Molecular Dynamics Simulations

To gain a deeper understanding of the TiCl4 solvation effects in multi-component ionic liquids, we performed ab initio molecular dynamics simulations of 1-butyl-3-methylimidazolium [C4C1Im]+, tetrafluoroborate [BF4]−, chloride [Cl]− both with and without water and titanium tetrachloride TiCl4. Complex interactions between cations and anions are observed in all investigated systems. By further addition of water and TiCl4 this complex interaction network is extended. Observations of the radial distribution functions and number integrals show that water and TiCl4 not only compete with each other to interact mainly with [Cl]−, which strongly influences the cation-[BF4]− interaction, but also interact with each other, which leads to the fact that in certain systems the cation-anion interaction is enhanced. Further investigations of the Voronoi polyhedra analysis have demonstrated that water has a greater impact on the nanosegregated system than TiCl4 which is also due to the fact of the shear amount of water relative to all other components and its higher mobility compared to TiCl4. Overall, the polar network of the IL mixture collapses by including water and TiCl4. In the case of [Cl]− chloride enters the water continuum, while [BF4]− remains largely unaffected, which deeply affects the interaction of the ionic liquid (IL) network.     

Correlation between Dynamic Heterogeneity and Local Structure in a Room‐Temperature Ionic Liquid: A Molecular Dynamics Study of [bmim][PF6]

The complex dynamics of a room‐temperature ionic liquid, 1‐n‐butyl‐3‐methylimidazolium hexafluorophosphate ([bmim][PF₆]), is studied using equilibrium classical molecular dynamics simulations in the temperature range of 250–450 K. The activation energies for the self‐diffusion of ions are around 30–34 kJ mol^?1, with that of the anion a little higher than that for the cation. The electrical conductivity of the liquid is calculated and good agreement with experiments is obtained. Structural relaxation is studied through the decay of coherent (total density–density correlation) and incoherent (self part of density–density correlation) intermediate scattering functions over a range of temperatures and wave vectors relevant to the system. The relaxation data are used to identify and characterize two processes, α and β. The dependence of the two relaxation times on temperature and wave vector is obtained. The dynamical heterogeneity of the ions determined through the non‐Gaussian parameter indicates the motion of the cation to be more heterogeneous than that of the anion. The faster ones among the cations are coordinated to faster anions, while slower cations are surrounded predominantly by slower anions. Thus, the dynamical heterogeneity in this ionic liquid is shown to have structural signatures.     

n-Butane,iso-Butane and 1-Butene Adsorption on Imidazolium-Based Ionic Liquids Studied with Molecular Beam Techniques

The interaction of molecules, especially hydrocarbons, at the gas/ionic liquid (IL) surface plays a crucial role in supported IL catalysis. The dynamics of this process is investigated by measuring the trapping probabilities of n-butane, iso-butane and 1-butene on a set of frozen 1-alkyl-3-methylimidazolium-based ILs [C_nC₁Im]X, where n=4, 8 and X⁻=Cl⁻, Br⁻, [PF₆]⁻ and [Tf₂N]⁻. The decrease of the initial trapping probability with increasing surface temperature is used to determine the desorption energy of the hydrocarbons at the IL surfaces. It increases with increasing alkyl chain length n and decreasing anion size for the ILs studied. We attribute these effects to different degrees of alkyl chain surface enrichment, while interactions between the adsorbate and the anion do not play a significant role. The adsorption energy also depends on the adsorbing molecule: It decreases in the order n-butane>1-butene>iso-butane, which can be explained by different dispersion interactions.     

Understanding Structural and Transport Properties of Dissolved Li2S8 in Ionic Liquid Electrolytes through Molecular Dynamics Simulations

Lithium-sulfur batteries with high energy density are considered as one of the most promising future energy storage devices. However, the parasitic lithium polysulfides shuttle phenomenon severely hinders the commercialization of such batteries. Ionic liquids have been found to suppress the lithium polysulfides solubility, diminishing the shuttle effect effectively. Herein, we performed classical molecular dynamics simulations to explore the microscopic mechanism and transport behaviors of typical Li₂S₈ species in ionic liquids and ionic liquid-based electrolyte systems. We found that the trifluoromethanesulfonate anions ([OTf]⁻) exhibit higher coordination strength with lithium ions compared with bis(trifluoromethanesulfonyl)imide anions ([TFSI]⁻) in static microstructures. However, the dynamical characteristics indicate that the presence of the [OTf]⁻ anions in ionic liquid electrolytes bring faster Li⁺ exchange rate and easier dissociation of Li⁺ solvation structures. Our simulation models offer a significant guidance to future studies on designing ionic liquid electrolytes for lithium-sulfur batteries.     

In the footsteps of August Michaelis: Syntheses and Thermodynamics of Extremely Low-Volatile Ionic Liquids

A series of nine different known ionic liquids or low melting salts was synthesised and purified. They are composed of the [NTf₂]^– (bis(trifluoromethane)sulfonimide), [OTf]^– (trifluoro-methane-sulfonate), or [B(CN)₄]^– (tetracyanidoborate) anion and [Ph₄P]⁺ (tetraphenylphosphonium), [Ph₃BzP]⁺ (triphenylbenzyl phosphonium), [ⁿBu₄P]⁺ (tetra-ⁿbutylphosphonium), [ⁿBuPh₃P]⁺ (tri-phenyl-ⁿbutylphosphonium), [ⁿBu₄N]⁺ (tetra-ⁿbutylammonium), or the [PPN]⁺ (bis(triphenylphosphine)iminium) cation. Precise vapour pressure data and enthalpies of vaporisation were measured using the Quartz Crystal Microbalance (QCM) method and evaluated. Structure-property relations are established using the obtained data as well as literature known data of ILs with alkyl-substituted imidazolium cations. It turns out that ILs with the tetracyanidoborate anion have even higher values of the enthalpy of vaporisation than those with the common [NTf₂]⁻ or [OTf]⁻ anion and therefore are even less volatile.     

Ionic Liquids and Water: Hydrophobicity vs. Hydrophilicity

Many chemical processes rely extensively on organic solvents posing safety and environmental concerns. For a successful transfer of some of those chemical processes and reactions to aqueous media, agents acting as solubilizers, or phase-modifiers, are of central importance. In the present work, the structure of aqueous solutions of several ionic liquid systems capable of forming multiple solubilizing environments were modeled by molecular dynamics simulations. The effect of small aliphatic chains on solutions of hydrophobic 1-alkyl-3-methylimidazolium bis(trifluoromethyl)sulfonylimide ionic liquids (with alkyl = propyl [C₃C₁im][NTf₂], butyl [C₄C₁im][NTf₂] and isobutyl [iC₄C₁im][NTf₂]) are covered first. Next, we focus on the interactions of sulphonate- and carboxylate-based anions with different hydrogenated and perfluorinated alkyl side chains in solutions of [C₂C₁im][C_nF_2n+1SO₃], [C₂C₁im][C_nH_2n+1SO₃], [C₂C₁im][CF₃CO₂] and [C₂C₁im][CH₃CO₂] (n = 1, 4, 8). The last system considered is an ionic liquid completely miscible with water that combines the cation N-methyl-N,N,N-tris(2-hydroxyethyl)ammonium [N_{1 2OH 2OH 2OH}]⁺, with high hydrogen-bonding capability, and the hydrophobic anion [NTf₂]^–. The interplay between short- and long-range interactions, clustering of alkyl and perfluoroalkyl tails, and hydrogen bonding enables a wealth of possibilities in tailoring an ionic liquid solution according to the needs.     

Spectroscopic and Chemical Properties of Ionic Liquids: Computational Study

A brief account is given of highlights of our computational efforts – often in collaboration with experimental groups – to understand spectroscopic and chemical properties of ionic liquids (ILs). Molecular dynamics, including their inhomogeneous character, responsible for key spectral features observed in dielectric absorption, infra-red (IR) and fluorescence correlation spectroscopy (FCS) measurements are elucidated. Mechanisms of chemical processes involving imidazolium-based ILs are illustrated for CO₂ capture and related reactions, transesterification of cellulose, and Au nanocluster-catalyzed Suzuki cross-coupling reaction with attention paid to differing roles of IL ions. A comparison with experiments is also made.     

Aqueous Solutions of Ionic Liquids: Microscopic Assembly

Aqueous solutions of ionic liquids are of special interest, due to the distinctive properties of ionic liquids, in particular, their amphiphilic character. A better understanding of the structure–property relationships of such systems is hence desirable. One of the crucial molecular‐level interactions that influences the macroscopic behavior is hydrogen bonding. In this work, we conduct molecular dynamics simulations to investigate the effects of ionic liquids on the hydrogen‐bond network of water in dilute aqueous solutions of ionic liquids with various combinations of cations and anions. Calculations are performed for imidazolium‐based cations with alkyl chains of different lengths and for a variety of anions, namely, [Br]^?, [NO₃]^?, [SCN]^?, [BF₄]^?, [PF₆]^?, and [Tf₂N]^?. The structure of water and the water–ionic liquid interactions involved in the formation of a heterogeneous network are analyzed by using radial distribution functions and hydrogen‐bond statistics. To this end, we employ the geometric criterion of the hydrogen‐bond definition and it is shown that the structure of water is sensitive to the amount of ionic liquid and to the anion type. In particular, [SCN]^? and [Tf₂N]^? were found to be the most hydrophilic and hydrophobic anions, respectively. Conversely, the cation chain length did not influence the results.