Temperature-dependent structure of methyltributylammonium bis(trifluoromethylsulfonyl)amide: X ray scattering and simulations

Are NaTFSI and NaFSI Salt-Based Water-in-Salt Electrolytes Structurally Similar or Different?

Water-in-salt electrolytes (WiSEs) are a promising class of electrolytes due to their wide electrochemical stability window and nonflammability. In this study, we explore the structural organization of sodium bis(trifluoromethylsulfonyl)imide (NaTFSI) and sodium bis(fluorosulfonyl)imide (NaFSI) salt-based aqueous electrolytes, covering dilute to highly concentrated regions, by employing an all-atom molecular dynamics simulation. For the NaTFSI-based electrolyte, we observe that Na+ ions are mostly surrounded by water molecules at all the salt concentrations due to the very strong interaction between them. While TFSI anions weakly coordinate with Na+ ions and other TFSI anions, they also mostly prefer to be surrounded by water molecules. These interactions were found to have moderate dependence on the concentration of the NaTFSI salt. For the NaFSI-based electrolyte, while the Na+-water interaction is stronger at lower salt concentrations, the number of nearest neighbor FSI anions is found to be more than that of water at higher concentrations (≥20 m). This is because the increase in the salt concentration leads to expulsion of water molecules from the solvation shell of Na+ ions and enhances the interaction between Na+ ions and oxygen atoms of FSI. At the highest salt concentration (solubility limit), the bulk-like water structure is completely disrupted and dominated by an anionic network in the FSI-based electrolyte. In contrast, water-water hydrogen bonding network is still present even in the highly concentrated TFSI-based electrolyte. The simulated X-ray scattering pattern displays a low-q peak, revealing the presence of an intermediate range ordering due to alternating anion-rich and water/Na+-rich regions in both the electrolytes. However, the characteristic length scale corresponding to the low-q peak decreases with increasing the salt content in both the electrolytes.

A Computational and Spectroscopic Analysis of Solvate Ionic Liquids Containing Anions with Long and Short Perfluorinated Alkyl Chains

Anion-driven, nanoscale polar-apolar structural organization is investigated in a solvate ionic liquid (SIL) setting by comparing sulfonate-based anions with long and short perfluorinated alkyl chains. Representative SILs are created from 1,2-bis(2-methoxyethoxy)ethane ("triglyme" or "G3"), lithium nonafluoro-1-butanesulfonate, and lithium trifluoromethanesulfonate. Molecular dynamics simulations, density functional theory computations, and vibrational spectroscopy provide insight into the overall liquid structure, cation-solvent interactions, and cation-anion association. Significant competition between G3 and anions for cation-binding sites characterizes the G3-LiC4F9SO3 mixtures. Only 50% of coordinating G3 molecules form tetradentate complexes with Li+ in [(G3)1Li][C4F9SO3]. Moreover, the SIL is characterized by extensive amounts of ion pairing. Based on these observations, [(G3)1Li][C4F9SO3] is classified as a "poor" SIL, similar to the analogous [(G3)1Li][CF3SO3] system. Even though the comparable basicity of the CF3SO3- and C4F9SO3- anions leads to similar SIL classifications, the hydrophobic fluorobutyl groups support extensive apolar domain formation. These apolar moieties permeate throughout [(G3)1Li][C4F9SO3] and persist even at relatively low dilution ratios of [(G3)10Li][C4F9SO3]. By way of comparison, the CF3 group is far too short to sustain polar-apolar segregation. This demonstrates how chemically modifying the anions to include hydrophobic groups can impart unique nanoscale organization to a SIL. Moreover, tuning these nano-segregated fluorinated domains could, in principle, control the presence of dimensionally ordered states in these mixtures without changing the coordination of the lithium ions.

Effect of ionic groups on the morphology and transport properties in a novel perfluorinated ionomer containing sulfonic and phosphonic acid groups: a molecular dynamics study

Combining the phosphonic acid group with the sulfonic acid group in PEMs has been shown to be an effective strategy for improving the fuel cell performance. However, the interplay of two different ionic groups and the resulting effect on the membrane properties have not been fully elucidated. Here, we used classical molecular dynamics simulation to investigate the morphologies, transport properties and effects of ionic groups in a novel perfluorinated PEM containing two ionic groups (PFSA-PFPA) in comparison to the corresponding homopolymers. Phase separations between hydrophilic and hydrophobic domains are confirmed in these PEMs and result from the evolution of water clusters formed around the ionic groups. The combination of both ionic groups brings a complicated morphological feature in PFSA-PFPA, with near-cylindrical aqueous domains of large length scales interconnected by tortuous domains of small sizes. And we found that the self-diffusion coefficients of water molecules are strongly related to morphologies, with the water transport in PFSA-PFPA lying between two analogous homopolymers. At the molecular level, we found that the sulfonic and phosphonic acid groups have distinct effects on the coordination behaviors and the dynamics of water molecules and hydronium ions. Strong electrostatic interactions lead to compact coordination structures and sluggish dynamics of hydronium ions around phosphonic acid groups, which determine the morphological evolution and transport properties in PFSA-PFPA. Our study affords insights into the relationship between molecular characteristics and transport properties bridged by phase-separated morphologies in a novel PEM containing both sulfonic acid and phosphonic acid groups, which deepens the understanding of the interplay between two ionic groups and may inspire the rational design of high-performance PEMs.

The influence of charge ordering in the microscopic structure of monohydroxy alcohols

While radiation scattering data provides insight inside the microstructure of liquids, the Debye relation relating the scattering intensityI(k) to the atom-atom structure factorsSab(k)shows that, ultimately, it is these individual structure correlation functions which contain the relevant information about the micro-structure. However, these quantities are not observables, except in few cases where one can invert the Debye relation to obtain the structure functions. In the majority of other cases, the need for model dependent computer simulations is unavoidable. The resulting calculations reveal that the scattering pre-peak is the result of cancellations between positive pre-peaks and negative anti-peaks contributions from the atom-atom structure factors. What of systems where this cancellation is such that it entirely suppresses the scattering pre-peak? One would be tempted to falsely conclude that there is no underlying micro-heterogeneity. Hence, the structure functions appear as hidden variables, and it is important to understand the relation between their features and the micro-structure of the system. Through the computer simulation study of various mono-ols, ranging from methanol to 1-nonanol, as well as the branched octanols, we show how the features of the atom-atom pair correlation functiongab(r)affect that of the structure factorsSab(k), and reveal that the micro-structure is ultimately the result of the charge ordering between different atoms in the system.

Do Ionic Liquids Slow Down in Stages?

High impact recent articles have reported on the existence of a liquid-liquid (L-L) phase transition as a function of both pressure and temperature in ionic liquids (ILs) containing the popular trihexyltetradecylphosphonium cation (P666,14+), sometimes referred to as the "universal liquifier". The work presented here reports on the structural-dynamic pathway from liquid to glass of the most well-studied IL comprising the P666,14+ cation. We present experimental and computational evidence that, on cooling, the path from the room-temperature liquid to the glass state is one of separate structural-dynamic changes. The first stage involves the slowdown of the charge network, while the apolar subcomponent is fully mobile. A second, separate stage entails the slowdown of the apolar domain. Whereas it is possible that these processes may be related to the liquid-liquid and glass transitions, more research is needed to establish this conclusively.

Self-Assembled Nanostructures in Aprotic Ionic Liquids Facilitate Charge Transport at Elevated Pressure

Ionic liquids (ILs), revealing a tendency to form self-assembled nanostructures, have emerged as promising materials in various applications, especially in energy storage and conversion. Despite multiple reports discussing the effect of structural factors and external thermodynamic variables on ion organization in a liquid state, little is known about the charge-transport mechanism through the self-assembled nanostructures and how it changes at elevated pressure. To address these issues, we chose three amphiphilic ionic liquids containing the same tetra(alkyl)phosphonium cation and anions differing in size and shape, i.e., thiocyanate [SCN]-, dicyanamide [DCA]-, and tricyanomethanide [TCM]-. From ambient pressure dielectric and mechanical experiments, we found that charge transport of all three examined ILs is viscosity-controlled at high temperatures. On the other hand, ion diffusion is much faster than structural dynamics in a nanostructured supercooled liquid (at T < 210 ± 3 K), which constitutes the first example of conductivity independent from viscosity in neat aprotic ILs. High-pressure measurements and MD simulations reveal that the created nanostructures depend on the anion size and can be modified by compression. For small anions, increasing pressure shapes immobile alkyl chains into lamellar-type phases, leading to increased anisotropic diffusivity of anions through channels. Bulky anions drive the formation of interconnected phases with continuous 3D curvature, which render ion transport independent of pressure. This work offers insight into the design of high-density electrolytes with percolating conductive phases providing efficient ion flow.

Role of Charge Ordering in the Dynamics of Cluster Formation in Associated Liquids

Liquids are archetypes of disordered systems, yet liquids of polar molecules are locally more ordered than nonpolar molecules, due to the Coulomb interaction based charge ordering phenomenon. Hydrogen bonded liquids, such as water or alcohols, for example, represent a special type of polar liquids, in that they form labile clustered local structures. For water, in particular, hydrogen bonding and the related local tetrahedrality, play an important role in the various attempts to understand this liquid. However, labile structures imply dynamics, and it is not clear how it affects the understanding of this type of liquids from purely static point of view. Herein, we propose to reconsider hydrogen bonding as a charge ordering process. This concept allows us to demonstrate the insufficiency of the analysis of the microscopic structure based solely on static pair correlation functions, and the need for dynamical correlation functions, both in real and reciprocal space. The subsequent analysis allows to recover several aspects of our understanding of hydrogen bonded liquids, but from the charge order perspective. For water, it confirms the jump rotation picture found recently, and it allows to rationalize the contradicting pictures that arise when following the interpretations based on hydrogen bonding. For alcohols, it allows to understand the dynamical origin of the scattering prepeak, which does not exist for water, despite the fact that both these liquids have very similar hydroxyl group chain clusters. The concept of charge ordering complemented by the analysis of dynamical correlation functions appear as a promising way to understand microheterogeneity in complex liquids and mixtures from kinetics point of view.

How NaFTA salt affects the structural landscape and transport properties of Pyrr1,3FTA ionic liquid

Recently, it has been demonstrated that ionic liquids (ILs) with an asymmetric anion render a wider operational temperature range and can be used as a solvent in sodium ion batteries. In the present study, we examine the microscopic structure and dynamics of pure 1-methyl-1-propylpyrrolidinium fluorosulfonyl(trifluoromethylsulfonyl)amide (Pyrr_1,3FTA) IL using atomistic molecular dynamics simulations. How the addition of the sodium salt (NaFTA) having the same anion changes the structural landscape and transport properties of the pure IL has also been explored. The simulated x-ray scattering structure functions reveal that the gradual addition of NaFTA salt (up to 1.2 molal) suppresses the charge alternating feature of the pure IL because of the replacement of the Pyrr⁺ cations with the Na⁺ ions. The Na⁺ ions are majorly found near the oxygen atoms of the anions, but the probability of finding the Na⁺ ions near these atoms slightly decreases with increasing salt concentration. As expected, the Na⁺ ions stay away from the Pyrr⁺ cations. However, the probability of finding the anions around anions increases with increasing salt concentration. The simulated self-diffusion coefficients of the ions in the pure IL reveal slightly faster diffusion of the Pyrr⁺ cations as compared to the FTA^- anions. Interestingly, in the salt solution, despite having smaller size, the diffusion of the Na⁺ ions is found to be lesser than the Pyrr⁺ cations and the FTA^- anions. The analysis of the ionic conductivity and transport numbers reveals that the fractional contribution of the FTA^- anion to the overall conductivity remains nearly constant with increasing salt concentration, but the contribution of Pyrr⁺ cation decreases and Na⁺ ion increases.

Collective Effects in Ionic Liquid [emim][Tf2N] and Ionic Paramagnetic Nitrate Solutions without Long-Range Structuring

Mesoscopic shear elasticity has been revealed in ordinary liquids both experimentally by reinforcing the liquid/surface interfacial energy and theoretically by nonextensive models. The elastic effects are here examined in the frame of small molecules with strong electrostatic interactions such as room temperature ionic liquids [emim][Tf2N] and nitrate solutions exhibiting paramagnetic properties. We first show that these charged fluids also exhibit a nonzero low-frequency shear elasticity at the submillimeter scale, highlighting their resistance to shear stress. A neutron scattering study completes the dynamic mechanical analysis of the paramagnetic nitrate solution, evidencing that the magnetic properties do not induce the formation of a structure in the solution. We conclude that the elastic correlations contained in liquids usually considered as viscous away from any phase transition contribute in an effective way to collective effects under external stress whether mechanical or magnetic fields.

Microheterogeneity-Induced Vibrational Spectral Dynamics of Aqueous 1-Alkyl-3-methylimidazolium Tetrafluoroborate Ionic Liquids of Different Cationic Chain Lengths

We have monitored the impacts of an increment in the alkyl chain length of the imidazolium-based tetrafluoroborate ionic liquids on the local deuteroxyl probe modes of interest. For this study, we have taken 1-ethyl-3-methylimidazolium tetrafluoroborate [EMIm][BF₄], 1-butyl-3-methylimidazolium tetrafluoroborate [BMIm][BF₄], 1-octyl-3-methylimidazolium tetrafluoroborate [OMIm][BF₄], and 1-decyl-3-methylimidazolium tetrafluoroborate [DMIm][BF₄] ionic liquid solutions with 5% HOD in H₂O as the vibrational reporter of the associated ultrafast system dynamics. Classical molecular dynamics (MD) simulations were employed to determine molecular structure and dynamic properties, while the spectral profiles were derived by applying the wavelet analysis of classical trajectories. Spatial distribution functions reveal the heterogeneity within the molecular structures of the ionic liquids (ILs) with varying alkyl chain lengths. The intense position of the spectral peak, the frequency corresponding to the shoulder peak, and the spectral linewidth of the O-D stretch distribution are not influenced by the increment in the cationic chain length. In addition, the ionic liquid (IL) [BMIm][BF₄] exhibits a notable trend; the dynamic timescales are longer than the other studied systems. Therefore, we have performed the Voronoi decomposition analysis of the ionic and the polar-apolar domains, symmetrically increasing the length of alkyl chains on the IL cations. Domain analysis reveals structural microheterogeneity; the anions form discrete domains, and the ionic liquid constituting cations form continuous domains irrespective of the alkyl chain length on the imidazolium cations. Therefore, this computational ultrafast spectroscopy study aids in forming a molecular-level picture of the ionic liquid cations and anions in the liquid phase, providing a detailed interpretation of the spectral properties of the probe stretching vibrations.

Heterogeneity and Nanostructure of Superconcentrated LiTFSI-EmimTFSI Hybrid Aqueous Electrolytes: Beyond the 21 m Limit of Water-in-Salt Electrolyte

Ionic liquids such as EmimTFSI (1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide) have been found to improve the solubility of LiTFSI salt in water-in-salt electrolyte (WiSE) from 21 to 60 m. However, the molecular origin of such enhancement in the solubility is still unknown. In the present work, we elucidate the microscopic structures of LiTFSI-EmimTFSI-based hybrid aqueous electrolytes and compare them with the structure of LiTFSI-based WiSE using molecular dynamics simulations. Our analysis reveals the presence of alternating water-rich clusters and TFSI-rich extended domains in the WiSE. In these clusters and domains, the Li⁺ ions reside such that the total number of oxygen atoms around them is conserved to four, where water contributes about three oxygen atoms. The addition of EmimTFSI in the WiSE results in removal of water from the nearest-neighbor solvation shell of TFSI^- ions but not from the Li⁺ ions. Significant structural changes are observed when LiTFSI salt is further added to LiTFSI-EmimTFSI aqueous solution. In both the hybrid electrolytes, water and Emim⁺ cations are found to avoid each other. The simulated X-ray scattering structure factor reveals the presence of larger length-scale heterogeneity in the most concentrated solution of the hybrid aqueous electrolyte. We observe that this nanoscale heterogeneity originates from a water-TFSI-Emim-TFSI-water-TFSI-Emim-TFSI-like arrangement in which Li⁺ ions are dispersed such that the coordination number of oxygen atoms around them is enhanced to five, wherein the major contribution comes from the TFSI^- ions. We envision that the enhanced LiTFSI solubility originates from the replacement of water molecules with TFSI^- ions in the first solvation shell of Li⁺ ions.

Long-Range Organization Study of Piperidinium-Based Ionic Liquids by Polarizable Molecular Dynamics Simulations

The nanoscale organization of some classes of ionic liquids is responsible for their singular properties. In this paper, we use polarizable molecular dynamics simulations and small-angle X-ray scattering to probe the structure of two piperidinium- and (trifluoromethylsulfonyl)imide-based ionic liquids ([EBPip⁺][NTf₂^-] and [EOPip⁺][NTf₂^-]) that differ in the alkyl chain length of their cation. The X-ray scattering intensities calculated numerically, from the radial distribution functions, are in excellent agreement with the experimental data. The analysis of the different contributions of the X-ray scattering data allowed us to highlight the correlations responsible for the low q peak observed for the long-chain alkyl cations. New angular analyses showed that anions were more likely to align with alkyl chains as their size increased, inducing angular correlation between anions at larger distances. They also showed that the long alkyl chains of the cations aligned more with each other than the short ones. These more aligned alkyl chains induce a smaller volume of the apolar microdomains compared to the well-studied imidazolium-based ionic liquids, leading to the smaller correlation distance for piperidinium-based ionic liquids.

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.

Modelling biocompatible ionic liquids based on organic acids and amino acids: challenges for computational models and future perspectives

In this short review I shall highlight the basic principle and the difficulties that arise in attempting the computational modeling of seemingly simple systems which hide an unexpected complexity. Biocompatible ionic liquids which are based on the coupling of organic or amino acid anions with metabolic cations such as cholinium are the target of this review. These substances have been the subject of intense research activities in the last few years and have attracted the attention of computational chemists. I shall show that the computational description of these substances is far from trivial and requires the use of sophisticated techniques in order to account for a surprisingly rich chemistry that is due to several phenomena such as polarization, charge transfer, proton transfer equilibria and tautomerization reactions.

Unique and generic structural features of cholinium amino acid-based biocompatible ionic liquids

Cholinium amino acid-based (Ch-AA) biocompatible ionic liquids (bio-ILs) are synthesized from renewable components and are efficiently used for biomass processing. However, their microscopic structural features that lead to their application as biomass solvents remain undetermined. Herein, we use atomistic simulations to investigate the structures of six different Ch-AA bio-ILs up to the nanometer length scale and demonstrate that, depending on the anion side chain structure, the respective IL exhibits structural ordering at different length scales. All the six Ch-AA bio-ILs investigated here show a generic feature of having a strong hydrogen bonding network between the hydroxyl group of the cholinium cation and the carboxyl group of the amino acid anions. We illustrate that each of these bio-ILs also displays a unique feature. Distinctive intermediate range structural ordering leads to heterogeneity in methioninate- and phenylalaninate-based ILs caused by the anion side chain segregation. Intermediate range ordering is not observed in glutaminate- and glutamate-based ILs because significant anion side chain and backbone interactions hinder the formation of side chain clusters. Interestingly, for the cysteinate-based IL, the side chains do not interact with the backbones and the intermediate range ordering is not observed because of a shorter anionic side chain.

On the X-ray Scattering Pre-peak of Linear Mono-ols and the Related Microstructure from Computer Simulations

The X-ray scattering intensities (I(k)) of linear alkanols OH(CH₂)_n-1CH₃ obtained from experiments (methanol to 1-undecanol) and computer simulations (methanol to 1-nonanol) of different force field models are comparatively studied particularly in order to explain the origin and the properties of the scattering pre-peak in the k-vector range 0.3-1 Å^-1. The experimental I(k) values show two apparent features: the pre-peak position k_P decreases with increasing n, and more intriguingly, the amplitude A_P goes through a maximum at 1-butanol (n = 4). The first feature is well reproduced by all force-field models, while the second shows strong model dependence. The simulations reveal various shapes of clusters of the hydroxyl head-group from n>2. k_P is directly related to the size of the meta-objects corresponding to such clusters surrounded by their alkyl tails. The explanation of the A_P turnover at n = 4 is more involved in terms of cancellations of atom-atom structure factor S(k) contributions related to domain ordering. The flexibility of the alkyl tails tends to reduce the cross contributions, thus revealing the crucial importance of this parameter in the models. Force fields with all-atom representation are less successful in reproducing the pre-peak features for smaller alkanols, n<6, possibly because they blur the charge ordering process since all atoms bear partial charges. The analysis clearly shows that it is not possible to obtain a model-free explanation of the features of I(k).

SEM-Drude Model for the Accurate and Efficient Simulation of MgCl2-KCl Mixtures in the Condensed Phase

There is a long history of models that to different extents reproduce structural and dynamical properties of high-temperature molten salts. Whereas rigid ion models can work fairly well for some of the monovalent salts, polarizability is fundamentally important when small divalent or multivalent cations are combined with significantly polarizable anions such as Cl^- to form networked liquids that display a first sharp diffraction peak. There are excellent polarizable ion models (PIMs) for these systems, but there has been little success with the less expensive Core-Shell type models, which are often described as unwieldy or difficult to fit. In this article, we present the Sharma-Emerson-Margulis (SEM)-Drude model for MgCl₂/KCl mixtures that with the same ingredients used in the latest and most accurate PIM models overcome the aforementioned obstacles at significantly less computational cost; structural and dynamical properties are for all practical purposes very similar to what we obtain from the PIM but typical simulations can be more than 30 times faster. This has allowed us not only to expand our recent studies on the temperature and composition dependence of intermediate range order in MgCl₂/KCl mixtures but also to access transport properties that were simply too costly to properly sample in our recently published studies.

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.

Heterogeneity in the microstructure and dynamics of tetraalkylammonium hydroxide ionic liquids: insight from classical molecular dynamics simulations and Voronoi tessellation analysis

Microscopic structural and dynamic heterogeneities were investigated for three ionic liquids (ILs), tetraethylammonium hydroxide, tetrapropylammonium hydroxide, and tetrabutylammonium hydroxide employing classical molecular dynamics (MD) simulations.

Ether functionalisation, ion conformation and the optimisation of macroscopic properties in ionic liquids

Ionic liquids are an attractive material class due to their wide liquid range, intrinsic ionic conductivity, and high chemical as well as electrochemical stability. However, the widespread use of ionic liquids is hindered by significantly higher viscosities compared to conventional molecular solvents. In this work, we show how the transport properties of ionic liquids can be altered significantly, even for isostructural ions that have the same backbone. To this end, structure-property relationships have been determined for a set of 16 systematically varied representative ionic liquids. Variations in molecular structure include ammonium vs. phosphonium, ether vs. alkyl side chains, and rigid vs. flexible anions. Ab initio calculations are used to relate molecular structures to the thermal, structural and transport properties of the ionic liquids. We find that the differences in properties of ether and alkyl functionalised ionic liquids are primarily dependent on minimum energy geometries, with the conformational flexibility of ether side chains appearing to be of secondary importance. We also show unprecedented correlations between anion conformational flexibility and transport properties. Critically, increasing fluidity upon consecutive introduction of ether side chains and phosphonium centres into the cation is found to be dependent on whether the anion is flexible or rigid. We demonstrate that targeted design of functional groups based on structure-property relationships can yield ionic liquids of exceptionally high fluidity.

Elucidating Ionic Correlations Beyond Simple Charge Alternation in Molten MgCl2-KCl Mixtures

The development of technologies for nuclear reactors based on molten salts has seen a big resurgence. The success of thermodynamic models for these hinges in part on our ability to predict at the atomistic level the behavior of pure salts and their mixtures under a range of conditions. In this letter, we present high-energy X-ray scattering experiments and molecular dynamics simulations that describe the molten structure of mixtures of MgCl₂ and KCl. As one would expect, KCl is a prototypical salt in which structure is governed by simple charge alternation. In contrast, MgCl₂ and its mixtures with KCl display more complex correlations including intermediate-range order and the formation of Cl^--decorated Mg²⁺ chains. A thorough computational analysis suggests that intermediate-range order beyond charge alternation may be traced to correlations between these chains. An analysis of the coordination structure for Mg²⁺ ions paints a more complex picture than previously understood, with multiple accessible states of distinct geometries.

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.

Structural analysis of ionic liquids with symmetric and asymmetric fluorinated anions

Ionic liquids (ILs) with relatively low viscosities and broad windows of electrochemical stability are often constructed by pairing asymmetric cations with bisfluorosulfonylimide (FSI^-) or bistriflimide (NTf₂ ^-) anions. In this work, we systematically studied the structures of ILs with these anions and related perfluorobis-sulfonylimide anions with asymmetry and/or longer chains: (fluorosulfonyl)(trifluoromethylsulfonyl)imide (BSI_0,1 ^-), bis(pentafluoroethylsulfonyl)imide (BETI^-), and (trifluoromethylsulfonyl) (nonafluorobutylsulfonyl)imide (BSI_1,4 ^-) using high energy X-ray scattering and molecular dynamics simulation methods. 1-alkyl-3-methylimidazolium cations with shorter (ethyl, Im_2,1 ⁺) and longer (octyl, Im_1,8 ⁺) hydrocarbon chains were selected to examine how the sizes of nonpolar hydrocarbon and fluorous chains affect IL structures and properties. In comparison with these, we also computationally explored the structure of ionic liquids with anions having longer fluorinated tails.

Temperature Dependence of Static Structure Factor Peak Intensities for a Pyrrolidinium-Based Ionic Liquid

Static structure factors ( S( q)) for many ionic liquids show low-wavenumber peaks whose intensities increase with increasing temperature. The greater peak intensities might seem to imply increasing intermediate-range order with increasing temperature. Molecular dynamics (MD) simulations for a representative ionic liquid, 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide (C₄C₁pyrrTFSI), were used to calculate S( q) and partial S( q) (cation-cation, anion-anion, and cation-anion) at 298, 363, and 500 K. S( q) and partial S( q) were further decomposed into positive and negative components (which each indicate structural ordering) by separately summing positive and negative Fourier transform summands. Increasing temperature causes the negative components of each partial S( q) to decrease in magnitude more than the positive components, causing the total S( q) to increase in magnitude. Thus, structural ordering with periodicities corresponding to observed peaks in S( q) does not increase but instead decoheres with increasing temperature, even though S( q) peak heights increase. Fourier transform summands also show where in real space the positive and negative component contributions to S( q) change when the temperature increases. This new, detailed analysis based on Fourier transform summands comprising S( q) argues for great caution when interpreting S( q) intensities and highlights the value of simulations as a complement to X-ray (or neutron) scattering experiments.

Cation influence on heterocyclic ammonium ionic liquids: a molecular dynamics study

Four different ionic liquids (ILs) consisting of the bis(trifluoromethanesulfonyl)imide ([NTf₂]⁻) anion, with structurally similar systematically varying cations, are investigated herein through classical molecular dynamics.

Neat ionic liquids versus ionic liquid mixtures: a combination of experimental data and molecular simulation

Simple mixtures of ionic liquids (IL–IL mixtures) can become a promising approach for the substitution of task-specific ILs.

Structural analysis of zwitterionic liquids vs. homologous ionic liquids

Zwitterionic liquids (Zw-ILs) have been developed that are homologous to monovalent ionic liquids (ILs) and show great promise for controlled dissolution of cellulosic biomass. Using both high energy X-ray scattering and atomistic molecular simulations, this article compares the bulk liquid structural properties for novel Zw-ILs with their homologous ILs. It is shown that the significant localization of the charges on Zw-ILs leads to charge ordering similar to that observed for conventional ionic liquids with monovalent anions and cations. A low-intensity first sharp diffraction peak in the liquid structure factor S(q) is observed for both the Zw-IL and the IL. This is unexpected since both the Zw-IL and IL have a 2-(2-methoxyethoxy)ethyl (diether) functional group on the cationic imidazolium ring and ether functional groups are known to suppress this peak. Detailed analyses show that this intermediate range order in the liquid structure arises for slightly different reasons in the Zw-IL vs. the IL. For the Zw-IL, the ether tails in the liquid are shown to aggregate into nanoscale domains.

Finding the best density functional approximation to describe interaction energies and structures of ionic liquids in molecular dynamics studies

Ionic liquids raise interesting but complicated questions for theoretical investigations due to the fact that a number of different inter-molecular interactions, e.g., hydrogen bonding, long-range Coulomb interactions, and dispersion interactions, need to be described properly. Here, we present a detailed study on the ionic liquids ethylammonium nitrate and 1-ethyl-3-methylimidazolium acetate, in which we compare different dispersion corrected density functional approximations to accurate local coupled cluster data in static calculations on ionic liquid clusters. The efficient new composite method B97-3c is tested and has been implemented in CP2K for future studies. Furthermore, tight-binding based approaches which may be used in large scale simulations are assessed. Subsequently, ab initio as well as classical molecular dynamics simulations are conducted and structural analyses are presented in order to shed light on the different short- and long-range structural patterns depending on the method and the system size considered in the simulation. Our results indicate the presence of strong hydrogen bonds in ionic liquids as well as the aggregation of alkyl side chains due to dispersion interactions.

Direct experimental observation of mesoscopic fluorous domains in fluorinated room temperature ionic liquids

Fluorinated room temperature ionic liquids (FRTILs) represent a class of solvent media that are attracting great attention due to their IL-specific properties as well as features stemming from their fluorous nature. Medium-to-long fluorous tails constitute a well-defined apolar moiety in the otherwise polar environment. Similarly to the case of alkyl tails, such chains are expected to result in the formation of self-assembled fluorous domains. So far, however, no direct experimental observation has been made of the existence of such structural heterogeneities on the nm scale. We report here the first experimental evidence of the existence of mesoscopic spatial segregation of fluorinated domains, on the basis of highly complementary X-ray and neutron scattering data sets (highlighting the importance of the latter probe) and NMR spectroscopy. Data are interpreted using atomistic molecular dynamics simulations, emphasizing the existence of a self-assembly mechanism that delivers segregated fluorous domains, where preferential solubilisation of fluorinated compounds can occur, thus paving the way for several smart applications.