A new computational methodology for the characterization of complex molecular environments using IR spectroscopy: bridging the gap between experiments and computations

The molecular interactions and dynamics of complex liquid solutions are now routinely measured using IR and 2DIR spectroscopy. In particular, the use of the latter allows the determination of the frequency fluctuation correlation function (FFCF), while the former provides us with the average frequency. In turn, the FFCF can be used to quantify the vibrational dynamics of a molecule in a solution, and the center frequency provides details about the chemical environment, solvatochromism, of the vibrational mode. In simple solutions, the IR methodology can be used to unambiguously assign the interactions and dynamics observed by a molecule in solution. However, in complex environments with molecular heterogeneities, this assignment is not simple. Therefore, a method that allows for such an assignment is essential. Here, a parametrization free method, called Instantaneous Frequencies of Molecules or IFM, is presented. The IFM method, when coupled to classical molecular simulations, can predict the FFCF of a molecule in solutions. Here, N-methylacetamide (NMA) in seven different chemical environments, both simple and complex, is used to test this new method. The results show good agreement with experiments for the NMA solvatochromism and FFCF dynamics, including characteristic times and amplitudes of fluctuations. In addition, the new method shows equivalent or improved results when compared to conventional frequency maps. Overall, the use of the new method in conjunction with molecular dynamics simulations allows unlocking the full potential of IR spectroscopy to generate molecular maps from vibrational observables, capable of describing the interaction landscape of complex molecular systems.

Crystallization Kinetics of Phosphonium Ionic Liquids: Effect of Cation Alkyl Chain Length and Thermal History

The effects of alkyl chain length on the crystallization kinetics and ion mobility of tetraalkylphosphonium, [P666,n][TFSI], (n = 2, 6, 8, and 12) ionic liquids were studied by differential scanning calorimetry (DSC) and broadband dielectric spectroscopy (BDS) over a wide temperature range. The liquid-glass transition temperature (Tg) and ion dynamics examined over a broad T range were almost insensitive to structural modifications of the phosphonium cation. In contrast, the crystallization kinetics were strongly affected by the length of the fourth alkyl chain. Furthermore, the thermal history of the sample (cold vs melt crystallization) significantly impacted the crystallization rate. It has been found that the nature of crystallization phenomena is the same across the homologous series, while the kinetic aspect differs. Finally, electric conductivity in supercooled liquid and crystalline solid phases was measured for all samples, revealing significant ionic conductivity, largely independent of the cation structure.

Ordering of ionic liquids at a charged sapphire interface: Evolution with cationic chain length

HYPOTHESIS

Room Temperature Ionic Liquids (RTILs) bulk's molecular layering dominates their structure also at the RTIL/sapphire interface, increasing the layer spacing with the cationic alkyl chain length n. However, the negatively-charged sapphire surface compresses the layers, increases the layering range, and affects the intra-layer structure in yet unknown ways.

EXPERIMENTS

X-ray reflectivity (XR) off the RTIL/sapphire interface, for a broad homologous RTIL series 1-alkyl-3-methylimidazolium bis(trifluoromethansulfonyl)imide, hitherto unavailable for any RTIL.

FINDINGS

RTIL layers against the sapphire, exhibit two spacings: da and db. da is n-varying, follows the behavior of the bulk spacing but exhibits a downshift, thus showing significant layer compression, and over twofold polar slab thinning. The latter suggests exclusion of anions from the interfacial region due to the negative sapphire charging by x-ray-released electrons. The layering range is larger than the bulk's. db is short and near n-independent, suggesting polar moieties' layering, the coexistence mode of which with the da-spaced layering is unclear. Comparing the present layering with the bulk's and the RTIL/air interface's provides insight into the Coulomb and dispersion interaction balance dominating the RTIL's structure and the impact thereon of the presence of a charged solid interface.

Bulk Nanostructure of Mixtures of Choline Arginate, Choline Lysinate, and Water

Neutron diffraction with empirical potential structure refinement was used to investigate the bulk liquid nanostructure of mixtures of choline arginate (Ch[Arg]), choline lysinate (Ch[Lys]), and water at mole ratios of 1Ch[Arg]:1Ch[Lys]:6H2O (balanced), 1Ch[Arg]:1Ch[Lys]:20H2O (balanced dilute), 3Ch[Arg]:1Ch[Lys]:12H2O (Arg- rich), and 1Ch[Arg]:3Ch[Lys]:12H2O (Lys- rich). The Arg- and Lys- anions tend not to associate due to electrostatic repulsion between charge groups and weak anion-anion attractions. This means that the local ion structures around the anions in these mixtures resemble the parent single-component systems. The bulk liquid nanostructure varies with the Arg-:Lys- ratio. In the Lys--rich mixture (1Ch[Arg]:3Ch[Lys]:12H2O), Lys- side chains cluster into a continuous apolar domain separated from a charged domain of polar groups. In the balanced mixture (1Ch[Arg]:1Ch[Lys]:6H2O), Lys- side chains form discrete apolar aggregates within a continuous polar domain of Arg-, Ch+, and water, and in the Arg--rich mixture (3Ch[Arg]:1Ch[Lys]:12H2O), the distribution of Lys- and Arg- is nearly homogeneous. Finally, in the balance dilute system (1Ch[Arg]:1Ch[Lys]:20H2O), a percolating water domain forms.

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.

Small-Angle X-ray Scattering Study of the Amphiphilic Bulk Nanostructure of Tetraalkylammonium Deep Eutectic Solvents

Deep eutectic solvents (DESs) are low-melting mixtures, often prepared from a salt and a molecular hydrogen bond donor. Like ionic liquids, DESs that contain at least one sufficiently amphiphilic component can form bicontinuous nanostructures consisting of polar and nonpolar domains, although this has not been widely explored for many DES combinations. Here, the bulk nanostructures of DESs comprising tetraalkylammonium bromide salts (tetrabutylammonium bromide, tetraoctylammonium bromide, and methyltrioctylammonium bromide) with alkanols and alkanoic acids of systematically varied chain lengths (C2, C6, C8, and C10) as hydrogen bond donors have been studied. Small-angle X-ray scattering techniques were used to identify the relationship between the alkyl chain length and functionality of the hydrogen bond donor on the nature of the amphiphilic nanostructures formed. These findings demonstrated that the amphiphilic nanostructures of the DESs were not affected by the functional group on the hydrogen bond donor, with these nanostructures influenced primarily by both the absolute and relative alkyl chain lengths of the salt and hydrogen bond donor.

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.

Attraction between Like Charged Ions in Ionic Liquids: Unveiling the Enigma of Tetracyanoborate Anions

Intermolecular interactions in ionic liquids are mainly governed by Coulombic forces. Attraction between cations has been previously observed and was attributed to dispersion interactions between nonpolar moieties, hydrogen bonding, or π stacking. In this study, we present the intriguing behavior of tetracyanoborate anions in ionic liquids that, unlike their dicyanamide and tricyanomethanide counterparts, form dimers in both solid and liquid phases. A joint simulation and experimental study uncovers the origin of such anion-anion attraction: stabilization by induction and dispersion forces between several cyano groups, which is strong enough to overcome electrostatic repulsion. These findings open up new opportunities in the rational design of ionic liquids, where interactions between ions of the same charge can be controlled and fine-tuned by the presence of cyano groups.

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.

Nanostructure in Amphiphile-Based Deep Eutectic Solvents

Deep eutectic solvents (DESs) are an emerging class of modern, often "green" solvents with unique properties. Recently, a deep eutectic system based on amphiphilic surfactant N-alkyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (C12 & C14 sulfobetaine) and (1S)-(+)-10-camphor-sulfonic acid in the molar ratio 1:1.5 has been reported. Nanostructuring can be expected in this DES due to the nature of the components. In this work, we have investigated the native nanostructure in the DES comprising C12-C18 alkyl chain sulfobetaines with camphor sulfonic acid and how it interacts with polar and nonpolar species, water and dodecane, respectively, using small angle neutron scattering. By using contrast variation to highlight the relative position of the solvent components and additives, we can resolve the structure of the solvent and how it changes upon interaction with water and dodecane. Scattering from the neat DES shows structures corresponding to the self-assembly of sulfobetaines; the size of the structure increases as the alkyl chain length of the sulfobetaines increases. Water and dodecane interact, respectively, with the hydrophilic and hydrophobic moieties in the DES structure, primarily the sulfobetaine, thereby swelling and solvating the entire structure. The extent of the shift of the peak position, and the swelling, depend on concentration of the additive. The solution phase organization and the interaction of polar and nonpolar species as observed here, have the potential to affect the ordering of inorganic or polymeric materials grown in such solvents, paving new avenues for templating applications.

Molecular Architecture Effects on Bulk Nanostructure in Bis(Orthoborate) Ionic Liquids

A series of 19 ionic liquids (ILs) based on phosphonium and imidazolium cations of varying alkyl-chain lengths with the orthoborate anions bis(oxalato)borate [BOB]- , bis(mandelato)borate, [BMB]- and bis(salicylato)borate, [BScB]- , are synthesized and studied using small-angle neutron scattering (SANS). All measured systems display nanostructuring, with 1-methyl-3-n-alkyl imidazolium-orthoborates forming clearly bicontinuous L3 spongelike phases when the alkyl chains are longer than C6 (hexyl). L3 phases are fitted using the Teubner and Strey model, and diffusely-nanostructured systems are primarily fitted using the Ornstein-Zernicke correlation length model. Strongly-nanostructured systems have a strong dependence on the cation, with molecular architecture variation explored to determine the driving forces for self-assembly. The ability to form well-defined complex phases is effectively extinguished in several ways: methylation of the most acidic imidazolium ring proton, replacing the imidazolium 3-methyl group with a longer hydrocarbon chain, substitution of [BOB]- by [BMB]- , or exchanging the imidazolium for phosphonium systems, irrespective of phosphonium architecture. The results suggest there is only a small window of opportunity, in terms of molecular amphiphilicity and cation:anion volume matching, for the formation of stable extensive bicontinuous domains in pure bulk orthoborate-based ILs. Particularly important for self-assembly processes appear to be the ability to form H-bonding networks, which offer additional versatility in imidazolium systems.

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.

Relating the Induced Free Charge Density Gradient in a Room-Temperature Ionic Liquid to Molecular-Scale Organization

We report on dilution-dependent changes in the local environments of chromophores incorporated into room-temperature ionic liquid (RTIL)-molecular solvent binary systems where the ionic liquid cation and molecular solvent possess the same alkyl chain length. We have used the RTIL 1-decyl-1-methylpyrrolidinium bis(trifluoromethanesulfonyl)imide (DMPyrr⁺TFSI^-) and the molecular solvent 1-decanol. Perylene was used as a non-polar probe, and cresyl violet (CV⁺) was used as a polar probe chromophore. We observe that in both regions there is a change in the chromophore local environments with increasing 1-decanol content. The changes in the nonpolar regions of the binary RTIL-molecular solvent system occur at a lower 1-decanol concentration than changes in the polar regions. Both chromophores reorient as oblate rotors in this binary system, allowing detailed information on the relative values of the Cartesian components of the rotational diffusion constants to be extracted from the experimental data. The induced free charge density gradient, ρ_f, known to exist in RTILs, persists to high 1-decanol content (1-decanol mole fraction of 0.75), with the structural details of the gradient being reflected in depth-dependent changes in the Cartesian components of the rotational diffusion constants of CV⁺. This is the first time that changes in molecular organization have been correlated with ρ_f.

The Nanostructure of Alkyl-Sulfonate Ionic Liquids: Two 1-Alkyl-3-methylimidazolium Alkyl-Sulfonate Homologous Series

The functionalization of polymers with sulfonate groups has many important uses, ranging from biomedical applications to detergency properties used in oil-recovery processes. In this work, several ionic liquids (ILs) combining 1-alkyl-3-methylimidazolium cations [C_nC₁im]⁺ (4 ≤ n ≤ 8) with alkyl-sulfonate anions [C_mSO₃]^- (4 ≤ m ≤ 8) have been studied using molecular dynamics simulations, totalizing nine ionic liquids belonging to two homologous series. The radial distribution functions, structure factors, aggregation analyses, and spatial distribution functions reveal that the increase in aliphatic chain length induces no significant change in the structure of the polar network of the ILs. However, for imidazolium cations and sulfonate anions with shorter alkyl chains, the nonpolar organization is conditioned by the forces acting on the polar domains, namely, electrostatic interactions and hydrogen bonding.

The structuring effect of the alkyl domains on the polar network of ionic liquid mixtures: a molecular dynamics study

By using molecular dynamics simulations, we investigate the structural and dynamic properties of mixtures of 1,3-dimethylimidazolium bis(trifluoromethanesulfonyl)imide, [C₁C₁im][Tf₂N] and 1-dodecyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide, [C₁₂C₁im][Tf₂N] (also C₁ and C₁₂ in short). Such mixtures feature an imidazolium bistriflimide salt with a very short alkyl chain, not giving rise to any nano-segregation as a pure component, with another one with a longer alkyl chain that exhibits a substantial nano-segregation as a pure liquid. As the mole fraction of the long-chain component C₁₂ is increased, the so-called pre-peak of the structure factor S(q), occurring in the region 1-3 nm^-1, shows a shift to higher values of the wavevector q, mirroring a decrease of the corresponding correlation length. Moreover, the intensity of the pre-peak strongly increases with the C₁₂ concentration. These results are in very good agreement with experimental X-ray scattering data in the literature. On the other hand, the diffusion of the ions is found to exhibit a simple behaviour consistent with the increased viscosity of the mixture, and these results are also in good agreement with NMR experimental data from the literature. The simulation results are rationalized as caused by a structuring effect, similar to the hydrophobic effect, of the alkyl domains of the C₁₂ component dissolved in the "solvent" represented by the C₁ cation, the Tf₂N^- anion and the C₁₂ cation head. In short, the exclusion of the alkyl chains from the polar network, a process mostly governed by electrostatic interactions, favours the formation of hydrophobic domains, which in turn exert a structuring effect on the ions of the polar domains, favouring a stronger ionic interaction. This is finally reflected in a shorter correlation length and a higher intensity of the pre-peak of the structure factor S(q) as the C₁₂ mole fraction is increased. At variance with the microscopic structure, the diffusion of all three types of ions is not strongly influenced by the nano-segregation and is essentially dependent on the viscosity of the mixture.

An Overview of Structure and Dynamics Associated with Hydrophobic Deep Eutectic Solvents and Their Applications in Extraction Processes

Recent development of novel water-immiscible green solvents known as hydrophobic deep eutectic solvents (HDESs) has opened the gates for applications requiring media where presence of water is undesirable. Ever since they were prepared, researchers have used HDESs in diverse fields such as extraction processes, CO 2 sequestration, membrane formation, and catalysis. The microstructure and dynamics associated with the species comprising HDESs guide their suitability for specific applications. For example, varying the alkyl tail length of HDES components significantly affects the dynamics of the components and thus helps in tuning the efficiency of extraction processes. The development of HDESs is still in infancy and very few theoretical studies are available in the literature that help in understanding the structure and dynamics of HDESs. This review highlights the recent work focused on the microscopic structure and dynamics of HDESs and their potential applications, particularly in extraction processes. We have also provided a glimpse of how the integration of experiments and computational techniques can help understand the mechanism of extraction processes.

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.

Polycation radius of gyration in a polymeric ionic liquid (PIL): the PIL melt is not a theta solvent

The conformation of the polycation in the prototypical polymeric ionic liquid (PIL) poly(3-methyl-1-aminopropylimidazolylacrylamide) bis(trifluoromethylsulfonyl)imide (poly(3MAPIm)TFSI) was probed using small-angle neutron scattering (SANS) and ultra-small-angle neutron scattering (USANS) at 25 °C and 80 °C. Poly(3MAPIm)TFSI contains microvoids which lead to intense low q scattering that can be mitigated using mixtures of hydrogen- and deuterium-rich materials, allowing determination of the polycation conformation and radius of gyration (R_g). In the pure PIL, the polycation adopts a random coil conformation with R_g = 52 ± 0.5 Å. In contrast to conventional polymer melts, the pure PIL is not a theta solvent for the polycation. The TFSI^- anions, which comprise 48% v/v of the PIL, are strongly attracted to the polycation and act like small solvent molecules which leads to chain swelling analogous to an entangled, semi-dilute, or concentrated polymer solution in a good solvent.

Nanostructure, electrochemistry and potential-dependent lubricity of the catanionic surface-active ionic liquid [P6,6,6,14] [AOT]

HYPOTHESIS

A catanionic surface-active ionic liquid (SAIL) trihexyltetradecylphosphonium 1,4-bis(2-ethylhexoxy)-1,4-dioxobutane-2-sulfonate ([P_6,6,6,14] [AOT]) is nanostructured in the bulk and at the interface. The interfacial nanostructure and lubricity may be changed by applying a potential.

EXPERIMENTS

The bulk structure and viscosity have been investigated using small angle X-ray scattering (SAXS) and rheometry. The interfacial structure and lubricity as a function of potential have been investigated using atomic force microscopy (AFM). The electrochemistry has been investigated using cyclic voltammetry.

FINDINGS

[P_6,6,6,14] [AOT] shows sponge-like bulk nanostructure with distinct interdigitation of cation-anion alkyl chains. Shear-thinning occurs at 293 K and below, but becomes less obvious on heating up to 313 K. Voltammetric analysis reveals that the electrochemical window of [P_6,6,6,14] [AOT] on a gold micro disk electrode exceeds the potential range of the AFM experiments and that negligible redox activity occurs in this range. The interfacial layered structure of [P_6,6,6,14] [AOT] is weaker than conventional ILs and SAILs, whereas lubricity is better, confirming the inverse correlation between the near-surface structure and lubricity. The adhesive forces of [P_6,6,6,14] [AOT] are lower at -1.0 V than at open circuit potential and +1.0 V, likely due to reduced electrostatic interactions caused by shielding of charge centres via long alkyl chains.

Ionic Liquids and Deep Eutectics as a Transformative Platform for the Synthesis of Nanomaterials

Ionic liquids (ILs) are becoming a revolutionary synthesis medium for inorganic nanomaterials, permitting more efficient, safer and environmentally benign preparation of high quality products. A smart combination of ILs and...

Structure and dynamics of highly concentrated LiTFSI/acetonitrile electrolytes

High salt concentration has been shown to induce increased electrochemical stability in organic solvent-based electrolytes. Accompanying the change in bulk properties is a structural ordering on mesoscopic length scales and changes in the ion transport mechanism have also been suggested. Here we investigate the local structure and dynamics in highly concentrated acetonitrile electrolytes as a function of salt concentration. Already at low concentrations ordering on microscopic length scales in the electrolytes is revealed by small angle X-ray scattering, as a result of correlations of Li+ coordinating clusters. For higher salt concentrations a charge alternation-like ordering is found as anions start to take part in the solvation. Results from quasi-elastic neutron spectroscopy reveal a jump diffusion dynamical process with jump lengths virtually independent of both temperature and Li-salt concentration. The jump can be envisaged as dissociation of a solvent molecule or anion from a particular Li+ solvation structure. The residence time, 50-800 ps, between the jumps is found to be highly temperature and Li-salt concentration dependent, with shorter residence times for higher temperature and lower concentrations. The increased residence time at high Li-salt concentration can be attributed to changes in the interaction of the solvation shell as a larger fraction of TFSI anions take part in the solvation, forming more stable solvation shells.

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.

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.

Evidence of a liquid-liquid transition in a glass-forming ionic liquid

A liquid-liquid transition (LLT) is a transformation from one liquid to another through a first-order transition. The LLT is fundamental to the understanding of the liquid state and has been reported in a few materials such as silicon, phosphorus, triphenyl phosphite, and water. Furthermore, it has been suggested that the unique properties of materials such as water, which is critical for life on the planet, are linked to the existence of the LLT. However, the experimental evidence for the existence of an LLT in many molecular liquids remains controversial, due to the prevalence and high propensity of the materials to crystallize. Here, we show evidence of an LLT in a glass-forming trihexyltetradecylphosphonium borohydride ionic liquid that shows no tendency to crystallize under normal laboratory conditions. We observe a step-like increase in the static dielectric permittivity at the transition. Furthermore, the sizes of nonpolar local domains and ion-coordination numbers deduced from wide-angle X-ray scattering also change abruptly at the LLT. We independently corroborate these changes in local organization using Raman spectroscopy. The experimental access to the evolution of local order and structural dynamics across a liquid-liquid transition opens up unprecedented possibilities to understand the nature of the liquid state.

Interplay of Structure and Dynamics in Lithium/Ionic Liquid Electrolytes: Experiment and Molecular Simulation

Despite their promising use in electrochemical and electrokinetic devices, ionic-liquid-based electrolytes often exhibit complex behavior arising from a subtle interplay of their structure and dynamics. Here, we report a joint experimental and molecular simulation study of such electrolytes obtained by mixing 1-butyl 3-methylimidazolium tetrafluoroborate with lithium tetrafluoroborate. More in detail, experiments consisting of X-ray scattering, pulsed field gradient NMR, and complex impedance spectroscopy are analyzed in the light of molecular dynamics simulations to probe the structural, dynamical, and electrochemical properties of this ionic-liquid-based electrolyte. Lithium addition promotes the nanostructuration of the liquid as evidenced from the appearance of a scattering prepeak that becomes more pronounced. Microscopically, using the partial structure factors determined from molecular dynamics, this prepeak is shown to correspond to the formation of well-ordered positive/negative charge series and also large aggregates (Li_n(BF₄)_4-m)^(4-m+n)-, which develop upon lithium addition. Such nanoscale ordering entails a drastic decrease in both the molecular mobility and ionic conductivity. In particular, the marked association of Li⁺ cations with four BF₄^- anions and long ion pairing times, which are promoted upon lithium addition, are found to severely hinder the Li⁺ transport properties.

Microscopic structural features of water in aqueous-reline mixtures of varying compositions

Reline, a mixture of urea and choline chloride in a 2 : 1 molar ratio, is one of the most frequently used deep eutectic solvents. Pure reline and its aqueous solution have large scale industrial use. Owing to the presence of active hydrogen bond formation sites, urea and choline cations can disrupt the hydrogen-bonded network in water. However, a quantitative understanding of the microscopic structural features of water in the presence of reline is still lacking. We carry out extensive all-atom molecular dynamics simulations to elucidate the effect of the gradual addition of co-solvents on the microscopic arrangements of water molecules. We consider four aqueous solutions of reline, between 26.3 and 91.4 wt%. A disruption of the local hydrogen-bonded structure in water is observed upon inclusion of urea and choline chloride. The extent of deviation of the water structure from tetrahedrality is quantified using the tetrahedral order parameter (q_tet). Our analyses show a monotonic increase in the structural disorder as the co-solvents are added. Increase in the q_tet values are observed when highly electro-negative hetero-atoms like nitrogen, oxygen of urea and choline cations are counted as partners of the central water molecules. Further insights are drawn from the characterization of the hydrogen-bonded network in water and we observe the gradual rupturing of water-water hydrogen bonds and their subsequent replacement by the water-urea hydrogen bonds. A negligible contribution from the hydrogen bonds between water and bulky choline cations has also been found. Considering all the constituents as the hydrogen bond partners we calculate the possibility of a successful hydrogen bond formation with a central water molecule. This gives a clear picture of the underlying mechanism of water replacement by urea.

Heterogeneity in hydrophobic deep eutectic solvents: SAXS prepeak and local environments

The observation of the prepeak in the simulated total X-ray scattering structure function (S(q)) reveals the presence of intermediate-range structural heterogeneity in hydrophobic deep eutectic solvents.

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.

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.

Capacitance of thin films containing polymerized ionic liquids

Electrode-polymer interfaces dictate many of the properties of thin films such as capacitance, the electric field experienced by polymers, and charge transport. However, structure and dynamics of charged polymers near electrodes remain poorly understood, especially in the high concentration limit representative of the melts. To develop an understanding of electric field-induced transformations of electrode-polymer interfaces, we have studied electrified interfaces of an imidazolium-based polymerized ionic liquid (PolyIL) using combinations of broadband dielectric spectroscopy, specular neutron reflectivity, and simulations based on the Rayleigh's dissipation function formalism. Overall, we obtained the camel-shaped dependence of the capacitance on applied voltage, which originated from the responses of an adsorbed polymer layer to applied voltages. This work provides additional insights related to the effects of molecular weight in affecting structure and properties of electrode-polymer interfaces, which are essential for designing next-generation energy storage and harvesting devices.

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.

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.