Spiers Memorial Lecture: From cold to hot, the structure and structural dynamics of dense ionic fluids

The structure of ionic liquids (ILs), which a decade or two ago was the subject of polite but heated debate, is now much better understood. This has opened opportunities to ask more sophisticated questions about the role of structure in transport, the structure of systems with ions that are not prototypical, and the similarity between ILs and other dense ionic fluids such as the high-temperature inorganic molten salts; let alone the fact that new areas of research have emerged that sprung from our collective understanding of the structure of ILs such as the deep eutectic solvents, the polymerized ionic liquids, and the zwitterionic liquids. Yet, our understanding of the structure of prototypical ILs may not be as complete as we think it to be, given that recent experiments appear to show that in some cases there may be more than one liquid phase resulting in liquid-liquid (L-L) phase transitions. This article presents a perspective on what we think are key topics related to the structure and structural dynamics of ILs and to some extent high-temperature molten salts.

Structural effects of water clusters on viscosity at high shear rates

In this study, we use molecular dynamics simulations of liquid water to investigate how shear thinning affects the viscosity of liquid water by structural changes of the hydrogen bond network. The effect of shear on viscosity can be divided into two parts: shear-induced destruction of the hydrogen bond network and the influence of the water structure on shear viscosity. First, strong shear destroys tetrahedral structures and thus reduces the connectivity of the hydrogen bond network. It is mainly because shear deformation, characterized by compression and expansion axes, respectively, triggers the destruction and formation of hydrogen bonds, resulting in anisotropic effects on water structures. At the same time, shear destroys large clusters and enhances the formation of small ones, resulting in a decrease in average cluster sizes. Second, the change of viscosity obeys a power law relationship with the change of hydrogen bond structures, highlighting a one-to-one correspondence between structure and property. Meanwhile, in order to explain why the structure affects viscosity, we define hydrogen-bond viscosity and find that the cooperative motion of the water structures can promote momentum transfer in the form of aggregations. Hydrogen-bond viscosity accounts for 5%-50% of the total viscosity. Our results elucidate that water structures are the important structural units to explain the change of water properties.

Dynamics in Quaternary Ionic Liquids with Non-Flexible Anions: Insights from Mechanical Spectroscopy

The present work investigates how mechanical properties and ion dynamics in ionic liquids (ILs) can be affected by ILs' design while considering possible relationships between different mechanical and transport properties. Specifically, we study mechanical properties of quaternary ionic liquids with rigid anions by means of Dynamical Mechanical Analysis (DMA). We are able to relate the DMA results to the rheological and transport properties provided by viscosity, conductivity, and diffusion coefficient measurements. A good agreement is found in the temperature dependence of different variables described by the Vogel-Fulcher-Tammann model. In particular, the mechanical spectra of all the measured liquids showed the occurrence of a relaxation, for which the analysis suggested its attribution to a diffusive process, which becomes evident when the ion dynamics are not affected by the fast structural reorganization of flexible anions on a local level.

Elucidating the interplay of local and mesoscale ion dynamics and transport properties in aprotic ionic liquids

Ion dynamics and charge transport in 1-methyl-3-octylimidazolium ionic liquids with chloride, bromide, tetrafluoroborate, tricyanomethanide, hexafluorophosphate, triflate, tetrachloroaluminate, bis(trifluoromethylsulfonyl)imide, and heptachlorodialuminate anions are investigated by broadband dielectric spectroscopy, rheology, viscometry, and differential scanning calorimetry. A detailed analysis reveals an anion and temperature-dependent separation of characteristic molecular relaxation rates extracted from various representations of the dielectric spectra. The separation in rates extracted from the electric modulus and conductivity formalisms is interpreted as an experimental signature of significant heterogeneity in the local ion dynamics associated with the structural glass transition, viscosity, and dc ion conductivity. It is further found that the degree of dynamic heterogeneity correlates with the strengths of slow dielectric and mechanical relaxations previously attributed to the dynamics of mesoscale solvophobic aggregates. Increasing local dynamic heterogeneity correlates with an increase in the strength of the slow, aggregate dielectric relaxation and a decrease in the strength of the slow, aggregate mechanical relaxation. Accordingly, increasing local dynamic heterogeneity, brought about by change in temperature and/or cation/anion chemical structure, correlates with an increase in the static dielectric permittivities and a decrease in the contribution of aggregate dynamics to the zero-shear viscosities. The established correlation provides a new ability to distinguish between the influence of mesoscale aggregate shape/morphology versus local and mesoscale ion dynamics on the transport properties of ionic liquids.

Flexibility is the key to tuning the transport properties of fluorinated imide-based ionic liquids

Ionic liquids are becoming increasingly popular for practical applications such as biomass processing and lithium-ion batteries. However, identifying ionic liquids with optimal properties for specific applications by trial and error is extremely inefficient since there are a vast number of potential candidate ions. Here we combine experimental and computational techniques to determine how the interplay of fluorination, flexibility and mass affects the transport properties of ionic liquids with the popular imide anion. We observe that fluorination and flexibility have a large impact on properties such as viscosity, whereas the influence of mass is negligible. Using targeted modifications, we show that conformational flexibility provides a significant contribution to the success of fluorination as a design element. Contrary to conventional wisdom, fluorination by itself is thus not a guarantor for beneficial properties such as low viscosity.

Volume viscosity and ultrasonic relaxation of ethanol-water mixtures studied by molecular dynamics simulations

The volume viscosity of ethanol-water mixtures at various compositions was calculated by means of equilibrium molecular dynamics simulations, and the results were compared with acoustic experiments. The volume viscosity exhibited a strong maximum around the ethanol mole fraction of 0.27, which is in agreement with experiments. The relaxation in the 100 ps regime, which had been revealed by ultrasonic spectroscopy, was also reproduced. Analysis of the two-body density and the adiabatic pressure fluctuation demonstrated that the large volume viscosity of the mixture originates from the long-range concentration fluctuation.

Coupling between Translational Diffusion of a Solute and Dynamics of the Heterogeneous Structure: Higher Alcohols and Ionic Liquids

Translational diffusion of nonpolar monoatomic solutes in a room-temperature ionic liquid and 1-octanol was studied by molecular dynamics simulation. The diffusion coefficient was evaluated in two different ways: (1) from the mean-square displacement of a freely diffusing solute and (2) from the time correlation function of force acting on a fixed solute. The diffusion of the free solute is much greater than the prediction of the Stokes-Einstein (SE) relation when the size of the solute is small, as has been reported by many experimental works. In contrast, the friction on fixed small solutes follows the SE relation. The mechanism of the solute diffusion in both solvents was then analyzed based on the coupling between the translational motion of the solute and the collective dynamics of the heterogeneous intermediate-range structure characteristic to these solvents. Analysis revealed that the coupling is present in all systems, but the relaxation is fast in the cases of free and small solutes. This suggests that the coupling can relax through the motion of the solute when the solute is free and small, while the relaxation of the heterogeneous structure itself is required for large or fixed solutes. The difference in the relaxation dynamics of the friction on the solute and the shear viscosity is explained as the coupling with different dynamic modes of the solvent. Therefore, the validity of the SE relation may not be a good criterion to judge whether the mechanisms of the diffusion and the viscosity are the same or not.

Reorganization Energy of Electron Transfer in Ionic Liquids

Bandshape analysis of charge-transfer optical bands in room-temperature ionic liquids (ILs) was performed to extract the reorganization energy of electron transfer. Remarkably, the reorganization energies in ILs are close to those in cyclohexane. This result runs against common wisdom in the field since conducting ILs, which are characterized by an infinite static dielectric constant, and nonpolar cyclohexane fall to the opposite ends of the polarity scale based on their dielectric constants. Theoretical calculations employing structure factors of ILs from molecular dynamics simulations support the low values of the reorganization energy. Standard dielectric arguments do not apply to solvation in ILs, and nonergodic reorganization energies are required for a quantitative analysis.

Pressing matter: why are ionic liquids so viscous?

Room temperature ionic liquids are considered to have huge potential for practical applications such as batteries. However, their high viscosity presents a significant challenge to their use changing from niche to ubiquitous. The modelling and prediction of viscosity in ionic liquids is the subject of an ongoing debate involving two competing hypotheses: molecular and local mechanisms versus collective and long-range mechanisms. To distinguish between these two theories, we compared an ionic liquid with its uncharged, isoelectronic, isostructural molecular mimic. We measured the viscosity of the molecular mimic at high pressure to emulate the high densities in ionic liquids, which result from the Coulomb interactions in the latter. We were thus able to reveal that the relative contributions of coulombic compaction and the charge network interactions are of similar magnitude. We therefore suggest that the optimisation of the viscosity in room temperature ionic liquids must follow a dual approach.

Relationship between the Relaxation of Ionic Liquid Structural Motifs and That of the Shear Viscosity

In a set of recent articles, we have highlighted that friction is highly inhomogeneous in a typical ionic liquid (IL) with charge networks that are stiff and charge-depleted regions that are soft. This has consequences not only for the dynamics of ILs but also for the transport properties of solutes dissolved in them. In this article, we explore whether the family of alkylimidazolium ILs coupled with bis(trifluoromethylsulfonyl)imide (with similar Coulombic interactions but different alkyl tails), when dynamically “equalized” by having a similar shear viscosity, display q-dependent structural relaxation time scales that are the same across the family. Our results show that this is not the case, and in fact, the relaxation of in-network charge alternation appears to be significantly affected by the presence of separate polar and apolar domains. However, we also find that if one was to assign weight factors to the relaxation of the structural motifs, charge alternation always contributes about the same amount (between 62.1 and 66.3%) across systems to the running integral of the stress tensor correlation function from which the shear viscosity is derived. Adjacency correlations between positive and negative moieties also contribute an identical amount if a prepeak is not present (about 38%) and a slightly smaller amount (about 28%) when intermediate range order exists. The prepeak only contributes about 6% to viscoelastic relaxation, highlighting that the dynamics of the smaller scale motifs is the most important.

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.

A Pictorial View of Viscosity in Ionic Liquids and the Link to Nanostructural Heterogeneity

Prototypical ionic liquids (ILs) are characterized by three structural motifs associated with (1) vicinal interactions, (2) the formation of positive-negative charge-alternating chains or networks, and (3) the alternation of these networks with apolar domains. In recent articles, we highlighted that the friction and mobility in these systems are nowhere close to being spatially homogeneous. This results in what one could call mechanical heterogeneity, where charge networks are intrinsically stiff and charge-depleted regions are softer, flexible, and mobile. This Letter attempts to provide a clear and visual connection between friction-associated with the dynamics of the structural motifs (in particular, the charge network)-and recent theoretical work by Yamaguchi linking the time-dependent viscosity of ILs to the decay of the charge alternation peak in the dynamic structure function. We propose that charge blurring associated with the loss of memory of where positive and negative charges are within networks is the key mechanism associated with viscosity in ILs. An IL will have low viscosity if a characteristic charge-blurring decorrelation time is low. With this in mind, engineering new low-viscosity ILs is reduced to understanding how to minimize this quantity.

Mesoscale Organization and Dynamics in Binary Ionic Liquid Mixtures

The impact of mesoscale organization on dynamics and ion transport in binary ionic liquid mixtures is investigated by broad-band dielectric spectroscopy, dynamic-mechanical spectroscopy, X-ray scattering, and molecular dynamics simulations. The mixtures are found to form distinct liquids with macroscopic properties that significantly deviate from weighted contributions of the neat components. For instance, it is shown that the mesoscale morphologies in ionic liquids can be tuned by mixing to enhance the static dielectric permittivity of the resulting liquid by as high as 100% relative to the neat ionic liquid components. This enhancement is attributed to the intricate role of interfacial dynamics associated with the changes in the mesoscopic aggregate morphologies in these systems. These results demonstrate the potential to design the physicochemical properties of ionic liquids through control of solvophobic aggregation.

Shear Thinning and Nonlinear Structural Deformation of Ionic Liquids with Long Alkyl Chains Studied by Molecular Dynamics Simulation

Nonequilibrium molecular dynamics (MD) simulation was performed on imidazolium-based ionic liquids of two different alkyl chain lengths, and shear-rate-dependent viscosity was evaluated. Compared with the frequency-dependent linear shear viscosity determined by equilibrium MD simulation, shear thinning occurs at the shear rate several times lower than that predicted by the Cox-Merz rule. The deformation of the structure factor was also evaluated as the function of shear rate. The onset of shear thinning corresponds to that of the nonlinearity in the deformation of the charge-alternation mode in both ionic liquids, which is in harmony with the result of our previous work that the shear relaxation of ionic liquids is mainly coupled to the structural relaxation of the charge-alternation mode. In the presence of the polar-nonpolar domain structure characteristic to ionic liquids with a long alkyl chain, in particular, the nonlinearity in the domain mode begins within the Newtonian regime of shear viscosity.

Mesoscale Aggregates and Dynamic Asymmetry in Ionic Liquids: Evidence from Depolarized Dynamic Light Scattering

Nanoscale structures in ionic liquids (ILs) are usually identified by X-ray or neutron scattering techniques and occur when the alkyl chains of the cations are long enough to show the tendency to segregate into apolar domains. In search of dynamic evidence for these nanostructures, different experimental techniques recently reported bimodal dynamic susceptibility spectra. In all cases, the faster process observed was ascribed to the structural α-relaxation and the slower one to the relaxation of long-lived aggregates. By contrast, we show by depolarized dynamic light scattering (DDLS) experiments on a systematic series of imidazolium-based ILs that the dynamics of the cation and anion are clearly separated for long alkyl chains. Therefore, the observation of a bimodal behavior is not related to any nanostructure but reflects the two-component nature of ILs. Thus, a consistent picture is obtained across different experimental methods, like dielectric and shear mechanical relaxation. Finally, the actual dynamic signature of nanostructures is identified for the first time as a weak feature in some of the DDLS spectra at even lower frequencies.

Translation-orientation coupling and Cox-Merz rule of liquid hexane

Equilibrium and non-equilibrium molecular dynamics simulations are performed on liquid hexane in order to clarify the origin of the Cox-Merz rule of liquids composed of chain-like molecules. The relation between the frequency-dependent complex shear viscosity and the shear-rate dependent nonlinear viscosity follows the Cox-Merz rule as expected. The slowest viscoelastic relaxation mode is explained by the translation-orientation coupling mechanism, and the saturation of the shear-induced orientational order is observed in the non-equilibrium simulation at the onset of the shear thinning. The origin of the Cox-Merz rule is discussed in terms of the translation-orientation coupling.

In an ionic liquid, high local friction is determined by the proximity to the charge network

Structural heterogeneity in Ionic Liquids (ILs) is to a large extent defined by nanoscale apolar pockets that act as spacers between strings of positive and negative charges that alternate. In contrast to this, recent work from our group and that of others appear to indicate that dynamic, energetic, and mechanical heterogeneities are governed by the charged part of the liquid. In this article, we study the dynamics of methane, a small apolar solute, in the family of ILs 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ( Im 1 , n + /NTf₂ ^-), with n = 2, 4, 8 at temperatures that make the viscosity for each liquid similar and around 8 cP. We do this in an attempt to equalize the effect of the solvent on the dynamics of the solute. In all cases, we find that solute proximity to charge-enhanced regions coincides with translationally caged regimes (high local friction) whereas the opposite is true in charge-depleted regions. In a way, these ILs behave like a liquid within a liquid where the charge network is the high friction component.