Origin of heterogeneous dynamics in local molecular structures of ionic liquids

Modulation of Diffusion Mechanism and Its Correlation with Complexation in Aqueous Deep Eutectic Solvents

Aqueous mixtures of deep eutectic solvents (DESs) have gained traction recently as an effective template to tailor their physicochemical properties. But detailed microscopic insights into the effects of water on the molecular relaxation phenomenon in DESs are not entirely understood. DESs are strong network-forming liquids due to the extensive hydrogen bonding and complex formation between their species, and therefore, water can behave as a controlled disruptor altering the microscopic structure and dynamics in DESs. In this study, the role of water in the diffusion mechanism of acetamide in the aqueous mixtures of DESs synthesized using acetamide and lithium perchlorate is investigated using molecular dynamics (MD) simulation and quasielastic neutron scattering (QENS). The acetamide dynamics comprises localized diffusion within transient cages and a jump diffusion process across cages. The jump diffusion process is observed to be strongly enhanced by about a factor of 10 as the water content in the system is increased. Meanwhile, the geometry of the localized dynamics is unaltered by addition of water, but the localized diffusion becomes significantly faster and more heterogeneous with increasing water concentration. The accelerating effects of water on localized diffusion are also substantiated by QENS experiments. The water concentration in the DES is observed to control the solvation structure of lithium ions, with the ions becoming significantly hydrated at 20 wt % water. The formation of interwater and water-acetamide hydrogen bonds is observed. The increase in water concentration is found to increase the number of H-bonds; however, their lifetimes are found to decrease substantially. Similarly, the lifetimes of acetamide-lithium complexes are also found to be diminished by increasing water concentration. A power-law scaling relationship between lifetimes and diffusion constants is established, elucidating the extent of coupling between diffusive processes and hydrogen bonding and microscopic complexation. This study demonstrates the ability to use water as an agent to probe the role of structural relaxation and complex lifetimes of diffusive processes at different time and length scales.

Orientational wetting and dynamical correlations toward glass transition on the surface of imidazolium-based ionic liquids

Surface induces many fascinating physical phenomena, such as dynamic acceleration, surface anchoring, and orientational wetting, and, thus, is of great interest to study. Here, we report classic molecular dynamics simulations on the free-standing surface of imidazolium-based ionic liquids (ILs) [C4mim][PF6] and [C10mim][PF6]. On [C10mim][PF6] surface, a significant orientational wetting is observed, with the wetting strength showing a diverging tendency. Depth of the wetting was captured from the density and orientational order profile by a static length, which remarkably increases below the temperature Tstat upon cooling down. The dynamical correlation length that measures the distance of surface-dynamics acceleration into the bulk was characterized via the spatial-dependent mobility. The translational correlation exhibits a similar drastic increment at Tstat, while the rotational correlation drastically increases at a lower temperature Trot. We connect these results to the dynamics in bulk liquids, by finding Tstat and Trot that correspond to the onset temperatures where the liquids become cooperative for translational and rotational relaxation, respectively. This signifies the importance of collective dynamics in the bulk on the orientational wetting and surface dynamics in the ILs.

Confinement Effects in Dynamics of Ionic Liquids with Polymer-Grafted Nanoparticles

Ionic liquid mixed with poly(methyl methacrylate)-grafted nanoparticle aggregates at low particle concentrations was shown to exhibit different dynamics and ionic conductivity than that of pure ionic liquid in our previous studies. In this work, we report on the quasi-elastic neutron scattering results on ionic liquid containing polymer-grafted nanoparticles at the higher particle concentration. The diffusivity of imidazolium (HMIM + ) cations of 1-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (HMIM-TFSI) in the presence of poly(methyl methacrylate)-grafted iron oxide nanoparticles and the ionic conductivity of solutions were discussed through the confinement. Analysis of the elastic incoherent structure factor suggested the confinement radius decreased with the addition of grafted particles in HMIM-TFSI/solvent mixture, indicating the confinement that is induced by the high concentration of grafted particles, shrinks the HMIM-TFSI restricted volume. We further conjecture that this enhanced diffusivity occurs as a result of the local ordering of cations within aggregates of poly(methyl methacrylate)-grafted particles.

Choline chloride and ethylene glycol based deep eutectic solvent (DES) versus hydroxyl functionalized room temperature ionic liquids (RTILs): assessing the differences in microscopic behaviour between the DES and RTILs

With the aim of understanding the differences in the behavior of deep eutectic solvents (DESs) and room temperature ionic liquids (RTILs) in terms of their structure, dynamics, and intra- and intermolecular interactions, three different ILs and one DES having similar functionalities (hydroxyl) have been investigated by using both ensembled average and single-molecule spectroscopic techniques. Specifically, for this purpose, a choline chloride based DES (ethaline) and three hydroxyl functionalized ILs (1-(2-hydroxyethyl)-3-imidazolium bis(trifluoromethanesulfonyl)imide ([OHEMIM][NTF₂]), N-(2-hydroxyl ethyl)-N-methylpyrrolidinium bis(trifluoromethanesulfonyl)imide ([OHEMPy][NTf₂]), and N-(2-hydroxyethyl)-N,N-dimethylpropan-1-aminium bis(trifluoromethanesulfonyl)imide ([OHC3CH][NTf₂])) are employed and investigated by EPR, time-resolved fluorescence, NMR and FCS studies. Estimation of polarity through EPR spectroscopy has revealed that the hydroxyl ILs employed in these studies are hyper-polar (close to water) in nature, whereas the polarity of the DES is found to be close to those of aliphatic polyhydroxy-alcohols. Interestingly, both time-resolved fluorescence anisotropy and FCS studies on these systems have suggested that the hydroxyl ILs are more dynamically heterogeneous than the DES. More interestingly, PFG-NMR measurements have indicated that the fluid structure of ethaline is relatively more associated as compared to those of the ILs despite the fact that all the cations have the same hydroxyl functionalities. All these investigations have essentially demonstrated that, despite having similar functionalities, both the DES and hydroxyl ILs employed in the present study exhibit microscopic behaviours that are significantly different from each other, indicating the interplay of various intermolecular interactions within the constituent species in governing the behaviours of these solvent systems.

Validation of Structural Grounds for Anomalous Molecular Mobility in Ionic Liquid Glasses

Ionic liquid (IL) glasses have recently drawn much interest as unusual media with unique physicochemical properties. In particular, anomalous suppression of molecular mobility in imidazolium IL glasses vs. increasing temperature was evidenced by pulse Electron Paramagnetic Resonance (EPR) spectroscopy. Although such behavior has been proven to originate from dynamics of alkyl chains of IL cations, the role of electron spin relaxation induced by surrounding protons still remains unclear. In this work we synthesized two deuterated imidazolium-based ILs to reduce electron-nuclear couplings between radical probe and alkyl chains of IL, and investigated molecular mobility in these glasses. The obtained trends were found closely similar for deuterated and protonated analogs, thus excluding the relaxation-induced artifacts and reliably demonstrating structural grounds of the observed anomalies in heterogeneous IL glasses.

Nonlinear measurements of kinetics and generalized dynamical modes. II. Application to a simulation of solvation dynamics in an ionic liquid

Solvation dynamics in ionic liquids show features that are often associated with supercooled liquids, including "stretched" nonexponential relaxation. To better understand the mechanism behind the stretching, the nonlinear mode-correlation methods proposed in Paper I [S. R. Hodge and M. A. Berg, J. Chem. Phys. 155, 024122 (2021)] are applied to a simulation of a prototypical ionic liquid. A full Green's function is recovered. In addition, specific tests for non-Gaussian dynamics are made. No deviations from Gaussian dynamics are found. This finding is incompatible with rate heterogeneity as a cause of the nonexponential relaxation and appears to be in conflict with an earlier multidimensional analysis of the same data. Although this conflict is not resolved here, this work does demonstrate the practicality of mode-correlation analysis in the face of finite datasets and calculations.

Phase Transitions and Electrochemical Properties of Ionic Liquids and Ionic Liquid-Solvent Mixtures

Recent advances in studies of ionic liquids (IL) and ionic liquid-solvent mixtures are reviewed. Selected experimental, simulation, and theoretical results for electrochemical, thermodynamical, and structural properties of IL and IL-solvent mixtures are described. Special attention is paid to phenomena that are not predicted by the classical theories of the electrical double layer or disagree strongly with these theories. We focus on structural properties, especially on distribution of ions near electrodes, on electrical double layer capacitance, on effects of confinement, including decay length of a dissjoining pressure between confinig plates, and on demixing phase transition. In particular, effects of the demixing phase transition on electrochemical properties of ionic liquid-solvent mixtures for different degrees of confinement are presented.

Probing the Effect of Side Alkyl Chain Length on the Structural and Dynamical Micro-heterogeneities in Dicationic Ionic Liquids

The molecular dynamics simulations and Voronoi tessellation analysis of two dicationic ionic liquids (DILs) including [C₅(mim)₂][NTf₂]₂ and [C₅(mim)₂C₄][NTf₂]₂ have been carried out to investigate the effects of side alkyl chain length on the structural and dynamical micro-heterogeneity of these DILs. Radial distribution functions (RDFs), spatial distribution functions (SDFs), and also neighborhood analysis of ions have been calculated to determine the arrangement of the nearest neighboring ions. To better understand the hydrogen-bonding network, microstructures, inter- and intramolecular orientations of ions in the studied DILs, different kinds of combined distribution functions (CDFs) were computed and analyzed. Also, qualitative and quantitative analyses of the structural heterogeneity were explored through total/partial structure factors, heterogeneity order parameters (HOPs), and domain analysis from Voronoi tessellation. The results showed that the side alkyl chains in DILs have significant effects on their micro-organizations in such a way that [C₅(mim)₂C₄][NTf₂]₂ with longer side chains has more microstructural heterogeneity than [C₅(mim)₂][NTf₂]₂ where the linkage alkyl chain is the same in both of them. Furthermore, to shed light on the dynamical heterogeneity, ion pair, ion cage, and hydrogen-bond stabilities and also the reorientation dynamics of ions have been investigated. Results demonstrated that local dynamics differences originate from local structural heterogeneity.

Ion Dissociation in Ionic Liquids and Ionic Liquid Solutions

The extent to which cations and anions in ionic liquids (ILs) and ionic liquid solutions are dissociated is of both fundamental scientific interest and practical importance because ion dissociation has been shown to impact viscosity, density, surface tension, volatility, solubility, chemical reactivity, and many other important chemical and physical properties. When mixed with solvents, ionic liquids provide the unique opportunity to investigate ion dissociation from infinite dilution in the solvent to a completely solvent-free state, even at ambient conditions. The most common way to estimate ion dissociation in ILs and IL solutions is by comparing the molar conductivity determined from ionic conductivity measurements such as electrochemical impedance spectroscopy (EIS) (which measure the movement of only the charged, i.e., dissociated, ions) with the molar conductivity calculated from ion diffusivities measured by pulse field gradient nuclear magnetic resonance spectroscopy (PFG-NMR, which gives movement of all of the ions). Because the NMR measurements are time-consuming, the number of ILs and IL solutions investigated by this method is relatively limited. We have shown that use of the Stokes-Einstein equation with estimates of the effective ion Stokes radii allows ion dissociation to be calculated from easily measured density, viscosity, and ionic conductivity data (ρ, η, λ), which is readily available in the literature for a much larger number of pure ILs and IL solutions. Therefore, in this review, we present values of ion dissociation for ILs and IL solutions (aqueous and nonaqueous) determined by both the traditional molar conductivity/PFG-NMR method and the ρ, η, λ method. We explore the effect of cation and anion alkyl chain length, structure, and interaction motifs of the cation and anion, temperature, and the strength of the solvent in IL solutions.

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.

The Effect of Phenyl Substitutions on Microstructures and Dynamics of Tetraalkylphosphonium Bis(trifluoro- methylsulfonyl)imide Ionic Liquids

Extensive atomistic simulations demonstrated that a gradual substitution of hexyl chains with phenyl groups in tetraalkylphosphonium cations results in remarkable changes in hydrogen bonding interactions, liquid structures and scattering structural functions, and rotational dynamics of hexyl chains and phenyl groups in tetraalkylphosphonium bis(trifluoromethylsulfonyl)imide ionic liquids. Hydrogen donor sites in hexyl chains present competitive characteristics with those in phenyl groups in coordinating anions, as well as their continuous and intermittent hydrogen bonding dynamics. Cation-cation and anion-anion spatial correlations show concomitant shift to short distances with decreased peak intensities with variations of cation structures, whereas cation-anion correlations have a distinct shift to large radial distances due to decreased associations of anions with neighboring cations. These microstructural changes are qualitatively manifested in shifts of prominent peaks for prevalent charge alternations and adjacency correlations between ion species in scattering structural functions. Meanwhile, rotational dynamics of hexyl chains speed up, which, in turn, slow down rotations of phenyl groups, whereas anions exhibit imperceptible changes in their rotational dynamics. These computational results are intrinsically correlated with conformational flexibilities, molecular sizes, and steric hindrance effects of phenyl groups in comparison with hexyl chains, and constrained distributions of anions around cations in heterogeneous ionic environments.

Magnetic Resonance Imaging and Molecular Dynamics Characterization of Ionic Liquid in Poly(ethylene oxide)-Based Polymer Electrolytes

Ternary systems consisting of polymers, lithium salts, and ionic liquids (ILs) are promising materials for the development of next-generation lithium batteries. The ternary systems combine the advantages of polymer-salt and IL-salt systems, thus providing media with high ionic conductivity and solid-like mechanical properties. In this work, we apply nuclear magnetic resonance 1H microimaging [magnetic resonance imaging (MRI)] techniques and molecular dynamics (MD) simulations to study the translational and rotational dynamics of the N-butyl-N-methylpyrrolidinium (PYR14) cation in poly(ethylene oxide) (PEO) matrices containing the lithium bis(trifluoromethanesulfonyl) imide salt (LiTFSI) and the PYR14TFSI IL. The analysis of diffusion-weighted images in PEO/LiTFSI/PYR14TFSI samples with varying mole ratios (10:1:x, with x = 1, 2, 3, and 4) shows, in a wide range of temperatures, a spatially heterogeneous distribution of PYR14 diffusion coefficients. Their weight-averaged values increase with IL content but remain well below the values estimated for the neat IL. The analysis of T2 (spin-spin relaxation) parametric images shows that the PEO matrix significantly hinders PYR14 rotational freedom, which is only partially restored by increasing the IL content. The MD simulations, performed on IL-filled cavities within the PEO matrix, reveal that PYR14 diffusion is mainly affected by Li/TFSI coordination within the IL phase. In agreement with MRI experiments, increasing the IL content increases the PYR14 diffusion coefficients. Finally, the analysis of MD trajectories suggests that Li diffusion mostly develops within the IL phase, although a fraction of Li cations is strongly coordinated by PEO oxygen atoms.

How Molecular Chiralities of Bis(mandelato)borate Anions Affect Their Binding Structures With Alkali Metal Ions and Microstructural Properties in Tetraalkylphosphonium Ionic Liquids

Spiroborate anion-based inorganic electrolytes and ionic liquids (ILs) have fascinating electrochemical and tribological properties and have received widespread attention in industrial applications. The molecular chiralities of spiroborate anions have a significant effect on the microstructures and macroscopic functionalities of these ionic materials in application and thus deserve fundamental consideration. In the current work, we performed quantum chemistry calculations to address the binding strength and coordination structures of chiral bis(mandelato)borate ([BMB]) anions with representative alkali metal ions, as well as the electronic properties of alkali metal ion-[BMB] ion pair complexes. The optimized [BMB] conformers are categorized into V-shaped, bent, and twisted structures with varied electrostatic potential contours and conformational energies and distinct alkali metal ion-[BMB] binding structures. Alkali metal ions have additional associations with phenyl groups in V-shaped [BMB] conformers owing to preferential cation-π interactions. Furthermore, the effects of the molecular chiralities of [BMB] anions on the thermodynamics and microstructural properties of tetraalkylphosphonium [BMB] ILs were studied by performing extensive atomistic interactions. Oxygen atoms in [BMB] anions have competitive hydrogen bonding interactions with hydrogen atoms in cations depending on the molecular chiralities and steric hindrance effects of [BMB] anions. However, the molecular chiralities of [BMB] anions have a negligible effect on the liquid densities of tetraalkylphosphonium [BMB] ILs and the spatial distributions of boron atoms in anions around phosphorous atoms in cations. Enlarging tetraalkylphosphonium cation sizes leads to enhanced cation-anion intermolecular hydrogen bonding and Coulombic interactions due to enhanced segregation of polar groups in apolar networks in heterogeneous IL matrices, as verified by scattering structural functions.

Dynamical properties of a room temperature ionic liquid: Using molecular dynamics simulations to implement a dynamic ion cage model

The transport behavior of ionic liquids (ILs) is pivotal for a variety of applications, especially when ILs are used as electrolytes. Many aspects of the transport dynamics of ILs remain to be understood. Here, a common ionic liquid, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (BmimNTf₂), was studied with molecular dynamics simulations. The results show that BmimNTf₂ displays typical structural relaxation, subdiffusive behavior, and a breakdown of the Stokes-Einstein diffusion relation as in glass-forming liquids. In addition, the simulations show that the translational dynamics, reorientation dynamics, and structural relaxation dynamics are well described by the Vogel-Fulcher-Tammann equation like fragile glass forming liquids. Building on previous work that employed ion cage models, it was found that the diffusion dynamics of the cations and anions were well described by a hopping process random walk where the step time is the ion cage lifetime obtained from the cage correlation function. Detailed analysis of the ion cage structures indicated that the electrostatic potential energy of the ion cage dominates the diffusion dynamics of the caged ion. The ion orientational relaxation dynamics showed that ion reorientation is a necessary step for ion cage restructuring. The dynamic ion cage model description of ion diffusion presented here may have implications for designing ILs to control their transport behavior.

Neutron scattering studies on short- and long-range layer structures and related dynamics in imidazolium-based ionic liquids

Alkyl-methyl-imidazolium ionic liquids CnmimX (n: alkyl-carbon number, X: anion) have short-range layer structures consisting of ionic and neutral (alkylchain) domains. To investigate the temperature dependences of the interlayer, interionic group, and inter-alkylchain correlations, we have measured the neutron diffraction (ND) of C16mimPF₆, C9.5mimPF₆, and C8mimPF₆ in the temperature region from 4 K to 470 K. The quasielastic neutron scattering (QENS) of C16mimPF₆ was also measured to study the dynamics of each correlation. C16mimPF₆ shows a first-order transition between the liquid (L) and liquid crystalline (LC) phases at T_c = 394 K. C8mimPF₆ exhibits a glass transition at T_g = 200 K. C9.5mimPF₆, which is a 1:3 mixture between C8mimPF₆ and C10mimPF₆, has both transitions at T_c = 225 K and T_g = 203 K. In the ND experiments, all samples exhibit three peaks corresponding to the correlations mentioned above. The widths of the interlayer peak at ca. 0.2 Å^-1 changed drastically at the L-LC transitions, while the interionic peaks at ca. 1 Å^-1 exhibited a small jump at T_c. The peak position and area of the three peaks did not change much at the transition. The structural changes were minimal at T_g. The QENS experiments demonstrated that the relaxation time of the interlayer motion increased tenfold at T_c, while those of other motions were monotonous in the whole temperature region. The structural and dynamical changes mentioned above are characteristic of the L-LC transition in imidazolium-based ionic liquids.

Electrostatic interactions in soft particle systems: mesoscale simulations of ionic liquids

Computer simulations provide a unique insight into the microscopic details, molecular interactions and dynamic behavior responsible for many distinct physicochemical properties of ionic liquids. Due to the sluggish and heterogeneous dynamics and the long-ranged nanostructured nature of ionic liquids, coarse-grained meso-scale simulations provide an indispensable complement to detailed first-principles calculations and atomistic simulations allowing studies over extended length and time scales with a modest computational cost. Here, we present extensive coarse-grained simulations on a series of ionic liquids of the 1-alkyl-3-methylimidazolium (alkyl = butyl, heptyl-, and decyl-) family with Cl, [BF4], and [PF6] counterions. Liquid densities, microstructures, translational diffusion coefficients, and re-orientational motion of these model ionic liquid systems have been systematically studied over a wide temperature range. The addition of neutral beads in cationic models leads to a transition of liquid morphologies from dispersed apolar beads in a polar framework to that characterized by bi-continuous sponge-like interpenetrating networks in liquid matrices. Translational diffusion coefficients of both cations and anions decrease upon lengthening of the neutral chains in the cationic models and by enlarging molecular sizes of the anionic groups. Similar features are observed in re-orientational motion and time scales of different cationic models within the studied temperature range. The comparison of the liquid properties of the ionic systems with their neutral counterparts indicates that the distinctive microstructures and dynamical quantities of the model ionic liquid systems are intrinsically related to Coulombic interactions. Finally, we compared the computational efficiencies of three linearly scaling O(N log N) Ewald summation methods, the particle-particle particle-mesh method, the particle-mesh Ewald summation method, and the Ewald summation method based on a non-uniform fast Fourier transform technique, to calculate electrostatic interactions. Coarse-grained simulations were performed using the GALAMOST and the GROMACS packages and hardware efficiently utilizing graphics processing units on a set of extended [1-decyl-3-methylimidazolium][BF4] ionic liquid systems of up to 131 072 ion pairs.

From Nanoscale to Microscale: Crossover in the Diffusion Dynamics within Two Pyrrolidinium-Based Ionic Liquids

Knowledge of the ion motion in room temperature ionic liquids (RTILs) is critical for their applications in a number of fields, from lithium batteries to dye-sensitized solar cells. Experiments on a limited number of RTILs have shown that on macroscopic time scales the ions typically undergo conventional, Gaussian diffusion. On shorter time scales, however, non-Gaussian behavior has been observed, similar to supercooled fluids, concentrated colloidal suspensions, and more complex systems. Here we characterize the diffusive motion of ionic liquids based on the N-butyl-N-methylpyrrolidinium (PYR₁₄) cation and bis(trifluoro methanesulfonyl)imide (TFSI) or bis(fluorosulfonyl)imide (FSI) anions. A combination of pulsed gradient spin-echo (PGSE) NMR experiments and molecular dynamics (MD) simulations demonstrates a crossover from subdiffusive behavior to conventional Gaussian diffusion at ∼10 ns. The deconvolution of molecular displacements into a continuous spectrum of diffusivities shows that the short-time behavior is related to the effects of molecular caging. For PYR₁₄FSI, we identify the change of short-range ion-counterion associations as one possible mechanism triggering long-range displacements.