Communication: Nanoscale structure of tetradecyltrihexylphosphonium based ionic liquids

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.

Liquid Structure of Ionic Liquids with [NTf2]- Anions, Derived from Neutron Scattering

The liquid structure of three common ionic liquids (ILs) was investigated by neutron scattering for the first time. The ILs were based on the bis(trifluoromethanesulfonyl)imide anion, abbreviated in the literature as [NTf2]- or [TFSI]-, and on the following cations: 1-ethyl-3-methylimidazolium, [C2mim]+; 1-decyl-3-methylimidazolium, [C10mim]+; and trihexyl(tetradecyl)phosphonium, [P666,14]+. Comparative analysis of the three ILs confirmed increased size of nonpolar nanodomains with increasing bulk of alkyl chains. It also sheds light on the cation-anion interactions, providing experimental insight into strength, directionality, and angle of hydrogen bonds between protons on the imidazolium ring, as well as H-C-P protons in [P666,14]+, to oxygen and nitrogen atoms in the [NTf2]-. The new Dissolve data analysis package enabled, for the first time, the analysis of neutron scattering data of ILs with long alkyl chains, in particular, of [P666,14][NTf2]. Results generated with Dissolve were validated by comparing outputs from three different models, starting from three different sets of cation charges, for each of the three ILs, which gave convergent outcomes. Finally, a modified method for the synthesis of perdeuterated [P666,14][NTf2] has been reported, with the aim of reporting a complete set of synthetic and data processing approaches, laying robust foundations that enable the study of the phosphonium ILs family by neutron scattering.

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.

Coarse-grained simulations of ionic liquid materials: from monomeric ionic liquids to ionic liquid crystals and polymeric ionic liquids

Ionic liquid (IL) materials are promising electrolytes with striking physicochemical properties for energy and environmental applications. Heterogeneous structures and transport quantities of monomeric and polymeric ILs are intrinsically intercorrelated and span multiple spatiotemporal scales, which is more feasible for coarse-grained (CG) simulations than atomistic modelling. Herein we constructed a novel CG model for ethyl-imidazolium tetrafluoroborate ILs with varied cation alkyl chains ranging from C₂ to C₂₀, and the interaction parameters were validated against representative static and dynamic properties that were obtained from atomistic reference simulations and experimental characterizations at relevant thermodynamic states. This CG model was extended to study thermotropic phase behaviors of monomeric ILs and to explore ion association structures and ion transport quantities in polymeric ILs with different architectures. A systematic analysis of structural and dynamical quantities identifies an evolution of liquid morphology from homogeneous to nanosegregated structures and then a smectic mesomorphism via a gradual lengthening of cation alkyl chains, and thereafter a distinct structural transition characterized by a monotonic decrease in orientational and translational order parameters in a sequential heating cascade. Backbone and pendant polymeric ILs exhibit evident anion association structures with cation monomers and polymer chains, and striking intra- and interchain coordinations between cation monomers owing to an intrinsic polymer architecture effect. Such a peculiar ion pairing association leads to a progressive increase in anion intrachain hopping probabilities, and a concomitant decrease in anion interchain hopping events with a gradual lengthening of polymeric ILs. The anion diffusivities in polymeric ILs are intrinsically correlated with ion pairing association lifetimes and ion structural relaxation times via a universal power law correlation D ∼ τ^-1, irrespective of polymer architectures.

A Brief Guide to the Structure of High-Temperature Molten Salts and Key Aspects Making Them Different from Their Low-Temperature Relatives, the Ionic Liquids

High-temperature molten salt research is undergoing somewhat of a renaissance these days due to the apparent advantage of these systems in areas related to clean and sustainable energy harvesting and transfer. In many ways, this is a mature field with decades if not already a century of outstanding work devoted to it. Yet, much of this work was done with pioneering experimental and computational setups that lack the current day capabilities of synchrotrons and high-performance-computing systems resulting in deeply entrenched results in the literature that when carefully inspected may require revision. Yet, in other cases, access to isotopically substituted ions make those pioneering studies very unique and prohibitively expensive to carry out nowadays. There are many review articles on molten salts, some of them cited in this perspective, that are simply outstanding and we dare not try to outdo those. Instead, having worked for almost a couple of decades already on their low-temperature relatives, the ionic liquids, this is the perspective article that some of the authors would have wanted to read when embarking on their research journey on high-temperature molten salts. We hope that this will serve as a simple guide to those expanding from research on ionic liquids to molten salts and vice versa, particularly, when looking into their bulk structural features. The article does not aim at being comprehensive but instead focuses on selected topics such as short- and intermediate-range order, the constraints on force field requirements, and other details that make the high- and low-temperature ionic melts in some ways similar but in others diametrically opposite.

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.

Photodetachment and Electron Dynamics in 1-Butyl-1-methyl-pyrrolidinium Dicyanamide

The ultrafast transient absorption spectrum of 1-butyl-1-methyl-pyrrolidinium dicyanamide, [Pyr_1,4⁺][DCA^-], was measured in the visible and near-infrared (IR) spectral regions. Excitation of the liquid at 4.6 eV created initially delocalized and highly reactive electrons that either geminately recombined (69%) or localized onto a cavity with a time constant of ∼300 fs. Electron localization was reflected in the evolution of the TA spectrum and the time-dependent loss of reactivity with a dichloromethane quencher. The delocalized initial state and spectrum of the free electrons were consistent with computational predictions by Xu and Margulis [ J. Phys. Chem. B, 2015, 119, 532-542] on excess electrons in [Pyr_1,4⁺][DCA^-]. The computational study considered two possible localization mechanisms for excess electrons, localization on ions, and localization on cavities. In the case of photogenerated electron-hole pairs, the results presented here demonstrate localization to cavities as the dominant channel. Following localization onto a cavity, the free electrons underwent solvation and loss of reactivity with the quencher with rates that slowed in time. The dynamics were similar to an analogous prior study on the related liquid [Pyr_1,x⁺][NTf₂^-]. One significant difference was the larger yield of free electrons from photoexcitation of [Pyr_1,4⁺][DCA^-]. This was found to primarily reflect more efficient localization onto cavities rather than a slower geminate recombination rate.

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.

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.

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.

Surface Activities of a Lipid Analogue Room-Temperature Ionic Liquid and Its Effects on Phospholipid Membrane

There are great efforts of synthesizing imidazolium-based ionic liquids (ILs) for developing new antibiotics as these molecules have shown strong antibacterial activities. Compared to a single-hydrocarbon-chained IL, the lipid analogues (LAs) with two chains are more effective. In the present study, the LA molecule ^MeIm(COOH)^Me(Oleylamine)Iodide has been synthesized and its surface activities along with the effectiveness in restructuring of a model cellular membrane have been quantified. The molecule is found to be highly surface active as estimated from the area-pressure isotherm of a monolayer of the molecules formed at the air-water interface. The X-ray reflectivity (XRR) studies of a monolayer dip-coated on a hydrophilic substrate have shown the structural properties of the layer which resembles to those of unsaturated phospholipids. The LA molecules are observed to fluidize a phospholipid bilayer formed by the saturated lipid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). At a lower surface pressure, the lipid monolayer of DPPC has exhibited a thickening effect at a low concentration of added LA and a thinning effect at higher concentration. However, at a high surface pressure of the monolayer, the thickness is found to decrease monotonically. The in-plane pressure-dependent interaction of LA molecules with model cellular membrane and the corresponding perturbation in the structure and physical properties of the membrane may be linked to the strong lysing effect of these types of molecules.

Structure and dynamics of the molten alkali-chloride salts from an X-ray, simulation, and rate theory perspective

Molten salts are of great interest as alternative solvents, electrolytes, and heat transfer fluids in many emerging technologies.

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.

Investigation of structural ordering in network forming ionic liquids: A molecular dynamics study

Molecular dynamics simulations reveal anomalous short- and medium-range ordering with increasing temperature in network-forming ionic liquids (NIL) consisting of alkyl-diammonium cations with long side chains of 6 carbon atoms and citrate anions (NIL 5-6). This effect is weaker, and only a short-range order is observed in equivalent systems with side chains shortened to 3 C atoms (NIL 5-3). The short-range ordering can be attributed to volume expansion during heating, but the intermediate range order requires volume expansion as well as an increase in temperature. We find that the cross (cation-anion) interactions are the major contributors to the observed trend and the development of complex 3D correlations in the two-particle correlation functions. The simulations suggest that the above phenomenon can be correlated to local trapping of cation molecules in a variety of configurations at lower temperatures where molecular shape distributions show great variability; as temperature increases, the distribution of molecular radii of gyration becomes narrower, enabling the increased ordering.

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.

Molecular dynamics simulations of ammonium/phosphonium-based protic ionic liquids: influence of alkyl to aryl group

The variation of the center atom in the cation from an N to a P-atom leads to improved physiochemical properties of protic ionic liquids (PILs) which are suitable for electrolyte applications. We present an atomistic simulations study to compare the effect of an alkyl or aryl group on trioctylammonium triflate ([HN(Oct)3][TFO]) and triphenylammonium triflate ([HN(Ph)3][TFO]) with their phosphonium analogues. We have computed the binding energy from quantum chemical calculations and physical properties such as the viscosity and the electrical conductivity of PILs from molecular dynamics simulations. The influence of the aromatic character in PILs is found to be significant to the physical properties. Gas phase quantum chemical calculations on clusters of ion pairs have revealed the presence of C-H/π interactions in aromatic PILs along with hydrogen bonding. The variation in strength of the ion-pair affinities is examined using electric-current correlation and velocity autocorrelation functions. The qualitative differences observed are due to the aromatic rings and change in the central atom of the quaternary cation from an N to a P-atom, substantiated quantitatively by diffusion coefficients and electrical conductivities. The relatively weaker ion-pair interactions and low binding energy (-73.34 kcal mol-1) lead to the highest electrical conductivity in [HP(Ph)3][TFO].

Structures of Ionic Liquids Having Both Anionic and Cationic Octyl Tails: Lamellar Vacuum Interface vs Sponge-Like Bulk Order

Numerous experimental and computational studies have shown that the structure of ionic liquids is significantly influenced by confinement and by interactions with interfaces. The nature of the interface can affect the immediate ordering of cations and anions, changing important rheological characteristics relevant to lubrication. Most studies suggest that such changes are local or short-ranged and that bulk properties are reestablished on a length scale of a few nanometers. The current study focuses on the 1-methyl-3-octylimidazolium octylsulfate ionic liquid for which both the cation and anion have moderate length linear alkyl tails. For this system, we find that the bulk phase is dominated by the very common sponge-like morphology characteristic of many ionic liquids. However, at the vacuum interface, a lamellar structure is observed that is not restricted to the vicinity of the surface but instead extends across the full 9 nm slab of our simulation. We suspect that in reality it could extend significantly beyond this.