Structure of 1-Alkyl-1-methylpyrrolidinium Bis(trifluoromethylsulfonyl)amide Ionic Liquids with Linear, Branched, and Cyclic Alkyl Groups

What is the origin of slow relaxation modes in highly viscous ionic liquids?

Room temperature ionic liquids (RTILs) are molten salts consisting entirely of ions and have over the past decades gained increased interest due to their high potential in applications. These structurally complex systems often display multiple relaxation modes in the response functions at lower frequencies, hinting to complex underlying mechanisms. While the existence of these multimodal spectra in the shear mechanical, dielectric, and light scattering response of RTILs has been confirmed multiple times, controversy still surrounds the origin. This paper, therefore, aims to provide additional insights into the multimodal spectra seen in RTILs by presenting new shear mechanical results on seven different RTILs: Pyr1n-TFSI with n = 4, 6, and 8; Pyr18-TFSI mixed with Li-TFSI in two high concentrations; and Cn-mim-BF4 with n = 3 and 8. Dynamic depolarized light scattering was also measured on one of the Pyr18-TFSI Li-salt mixtures. These specific cases were analyzed in detail and put into a bigger perspective together with an overview of the literature. Recent literature offers two specific explanations for the origin of the multimodal shear mechanical spectra: (1) cation-anion time scale separation or (2) combined cation-anion relaxation in addition to a dynamic signal from mesoscale aggregates at lower frequencies. However, neither of these two pictures can consistently explain all the results on different ionic liquids. Instead, we conclude that the origin of the multimodal spectrum is system specific. This underlines the complexity of this class of liquids and shows that great care must be taken when making general conclusions based on specific cases.

Heterogeneous Structure, Mechanisms of Counterion Exchange, and the Spacer Salt Effect in Complex Molten Salt Mixtures Including LaCl3

Complex molten chloride salt mixtures of uranium, magnesium, and sodium are top candidates for promising nuclear energy technologies to produce electricity based on molten salt reactors. From a local structural perspective, LaCl3 is similar to UCl3 and hence a good proxy to study these complex salt mixtures. As fission products, lanthanide salts and their mixtures are also very important in their own right. This article describes from an experimental and theory perspective how very different the structural roles of MgCl2 and NaCl are in mixtures with LaCl3. We find that, whereas MgCl2 becomes an integral part of multivalent ionic networks, NaCl separates them. In a recent article (J. Am. Chem. Soc. 2022, 144, 21751-21762) we have called the disruptive behavior of NaCl "the spacer salt effect". Because of the heterogeneous nature of these salt mixtures, there are multiple structural motifs in the melt, each with its particular free energetics. Our work identifies and quantifies these; it also elucidates the mechanisms through which Cl- ions exchange between Mg2+-rich and La3+-rich environments.

Nanostructure and Dynamics of Aprotic Ionic Liquids at Graphite Electrodes as a Function of Potential

Atomic force microscope (AFM) videos reveal the near-surface nanostructure and dynamics of the ionic liquids (ILs) 1-butyl-3-methylimidazolium dicyanamide (BMIM DCA) and 1-hexyl-3-methylimidazolium dicyanamide (HMIM DCA) above highly oriented pyrolytic graphite (HOPG) electrodes as a function of surface potential. Molecular dynamics (MD) simulations reveal the molecular-level composition of the nanostructures. In combination, AFM and MD show that the near-surface aggregates form via solvophobic association of the cation alkyl chains at the electrode interface. The diffusion coefficients of interfacial nanostructures are ≈0.01 nm2 s-1 and vary with the cation alkyl chain length and the surface potential. For each IL, the nanostructure diffusion coefficients are similar at open-circuit potential (OCP) and OCP + 1V, but BMIM DCA moves about twice as fast as HMIM DCA. At negative potentials, the diffusion coefficient decreases for BMIM DCA and increases for HMIM DCA. When the surface potential is switched from negative to positive, a sudden change in the direction of the nanostructure motion is observed for both BMIM DCA and HMIM DCA. No transient dynamics are noted following other potential jumps. This study provides a new fundamental understanding regarding the dynamics of electrochemically stable ILs at electrodes vital for the rational development of IL-based electrochemical devices.

Solvent effects on extractant conformational energetics in liquid-liquid extraction: a simulation study of molecular solvents and ionic liquids

Extractant design in liquid-liquid extraction (LLE) is a research frontier of metal ion separations that typically focuses on the direct extractant-metal interactions. However, a more detailed understanding of energetic drivers of separations beyond primary metal coordination is often lacking, including the role of solvent in the extractant phase. In this work, we propose a new mechanism for enhancing metal-complexant energetics with nanostructured solvents. Using molecular dynamics simulations with umbrella sampling, we find that the organic solvent can reshape the energetics of the extractant's intramolecular conformational landscape. We calculate free energy profiles of different conformations of a representative bidentate extractant, n-octyl(phenyl)-N,N-diisobutyl carbamoyl methyl phosphinoxide (CMPO), in four different solvents: dodecane, tributyl phosphate (TBP), and dry and wet ionic liquid (IL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf2N]). By promoting reorganization of the extractant molecule into its binding conformation, our findings reveal how particular solvents can ameliorate this unfavorable step of the metal separation process. In particular, the charge alternating nanodomains formed in ILs substantially reduce the free energy penalty associated with extractant reorganization. Importantly, using alchemical free energy calculations, we find that this stabilization persists even when we explicitly include the extracted cation. These findings provide insight into the energetic drivers of metal ion separations and potentially suggest a new approach to designing effective separations using a molecular-level understanding of solvent effects.

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.

Microstructures of Choline Amino Acid based Biocompatible Ionic Liquids

Bio-compatible ionic liquids (Bio-ILs) represent a class of solvents with peculiar properties and exhibit huge potential for their applications in different fields of chemistry. Ever since they were discovered, researchers have used bio-ILs in diverse fields such as biomass dissolution, CO2 sequestration, and biodegradation of pesticides. This review highlights the ongoing research studies focused on elucidating the microscopic structure of bio-ILs based on cholinium cation ([Ch]+ ) and amino acid ([AA]- ) anions using the state-of-the-arta b i n i t i o ${ab\hskip0.25eminitio}$ and classical molecular dynamics (MD) simulations. The microscopic structure associated with these green ILs guides their suitability for specific applications. ILs of this class differ in the side chain of the amino acid anions, and varying the side chain significantly affects the structure of these ILs and thus helps in tuning the efficiency of biomass dissolution. This review demonstrates the central role of the side chain on the morphology of choline amino acid ([Ch][AA]) bio-ILs. The seemingly matured field of bio-ILs and their employment in various applications still holds significant potential, and the insights on their microscopic structure would steer the field of target specific application of these green ILs.

Dielectric relaxation and dielectric decrement in ionic acetamide deep eutectic solvents: Spectral decomposition and comparison with experiments

Frequency-dependent dielectric relaxation in three deep eutectic solvents (DESs), (acetamide+LiClO4/NO3/Br), was investigated in the temperature range, 329 ≤ T/K ≤ 358, via molecular dynamics simulations. Subsequently, decomposition of the real and the imaginary components of the simulated dielectric spectra was carried out to separate the rotational (dipole-dipole), translational (ion-ion), and ro-translational (dipole-ion) contributions. The dipolar contribution, as expected, was found to dominate all the frequency-dependent dielectric spectra over the entire frequency regime, while the other two components together made tiny contributions only. The translational (ion-ion) and the cross ro-translational contributions appeared in the THz regime in contrast to the viscosity-dependent dipolar relaxations that dominated the MHz-GHz frequency window. Our simulations predicted, in agreement with experiments, anion-dependent decrement of the static dielectric constant (ɛs ∼ 20 to 30) for acetamide (ɛs ∼ 66) in these ionic DESs. Simulated dipole-correlations (Kirkwood g factor) indicated significant orientational frustrations. The frustrated orientational structure was found to be associated with the anion-dependent damage of the acetamide H-bond network. Single dipole reorientation time distributions suggested slowed down acetamide rotations but did not indicate presence of any "rotationally frozen" molecule. The dielectric decrement is, therefore, largely static in origin. This provides a new insight into the ion dependence of the dielectric behavior of these ionic DESs. A good agreement between the simulated and the experimental timescales was also noticed.

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.

Physical properties and nanostructuring of long-chained homobaric imidazolium ionic liquids

Understanding the structure-property relationship and nanoscopic behaviour of ionic liquids is of utmost importance for their potential applications. Focusing these studies on sets of homobaric ionic liquids could provide important insight into the effects of specific chemical groups on the overall interaction profile, bringing researchers one step closer to succesfully designing ionic liquids which are tailor-made for specific applications. This work focuses on ionic liquids with 12 total carbons on their side chains, studying both their bulk physical properties (such as densities and viscosities) and their nanostructuring. The results reveal that by keeping the total number of carbons constant, but arranging them differently around the imidazolium ring, either in a linear or in a branched-chain formation, can result in compounds with quite distinct properties. Some of those (such as diffusivity) appear to be more sensitive to symmetry variations, while others (such as density) are not significantly affected. X-ray scattering is used in order to get a clearer understanding of the nanostructuring of the studied compounds and to investigate to what extent the observed macroscopic properties are directly linked to the nanoscale ordering.

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...

In situ nanoscale evaluation of pressure-induced changes in structural morphology of phosphonium phosphate ionic liquid at single-asperity contacts

In this work, we perform atomic force microscopy (AFM) experiments to evaluate in situ the dependence of the structural morphology of trihexyltetradecylphosphonium bis(2-ethylhexyl) phosphate ([P_6,6,6,14][DEHP]) ionic liquid (IL) on applied pressure. The experimental results obtained upon sliding a diamond-like-carbon-coated silicon AFM tip on mechanically polished steel at an applied pressure up to 5.5 ± 0.3 GPa indicate a structural transition of confined [P_6,6,6,14][DEHP] molecules. This pressure-induced morphological change of [P_6,6,6,14][DEHP] IL leads to the generation of a lubricious, solid-like interfacial layer, whose growth rate increases with applied pressure and temperature. The structural variation of [P_6,6,6,14][DEHP] IL is proposed to derive from the well-ordered layering of the polar groups of ions separated by the apolar tails. These results not only shed new light on the structural organization of phosphonium-based ILs under elevated pressure, but also provide novel insights into the normal pressure-dependent lubrication mechanisms of ILs in general.

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.

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.

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.

Molecular simulation of osmometry in aqueous solutions of the BMIMCl ionic liquid: a potential route to force field parameterization of liquid mixtures

Despite widespread development and use of ionic liquids (ILs) in both academic and industrial research, computational force fields (FFs) for most of those are not available for a precise description of inter-species interactions in aqueous environments. In the scope of this study, by means of molecular simulations, the osmotic coefficient of an aqueous solution of an IL is calculated and used as a basis to reparameterize popular IL-FFs existing in the literature. We first calculate the osmotic coefficients (at 298.15 K and 1 atm pressure) of aqueous solutions of 1-butyl-3-methylimidazolium chloride (BMIMCl), a generic IL, popularly used in biomass processing and the subsequent conversion to value-added intermediates. The performance of two popular atomic, nonpolarizable FFs developed for BMIMCl, one by Lopes, Pádua, and coworkers (FF-LP) and the other by Sambasivarao, Acevedo, and coworkers (FF-SA), when mixed with the SPC/E water model, is tested with respect to their ability to reproduce the experimental osmotic coefficient data. Interestingly, the osmotic coefficient is found to be increasing with a gradual increase in IL molality within the concentration range of our investigation, which is contrary to the experimental trend reported in the literature for the same IL-water mixture. Henceforth, necessary corrections to the nonbonded ion-ion and ion-water interactions are made to match the experimental osmotic coefficient. To further assess the reliability of the new FF, we extensively explore the thermodynamic (density, isothermal compressibility, and thermal expansion coefficient), dynamic (diffusivity and viscosity), and association/dissociation properties (rationalized with the help of radial distribution functions) with both the original and reparameterized FF for a wider range of concentrations up to a molality of 18.50 mol kg-1. The calculated quantities are compared against experimental data wherever available. The modified FF parameters exhibit significant improvements in terms of its ability to match experimental solution properties, such as density, viscosity, association/dissociation, etc. We report that excessive dissociation of BMIMCl in water is responsible for the shortcomings observed in the original FFs and improved prediction of physicochemical properties could be achieved using the modified FFs.

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.

Self-Stabilized Giant Aggregates in Water from Room-Temperature Ionic Liquids with an Asymmetric Polar-Apolar-Polar Architecture

We report the assembly of four imidazolium bromides, each of which bears a naphthyl on one side of the imidazolium cation and a branched alkyl chain on the other. This design creates a new type of amphiphilic ionic liquid with an apolar-polar-apolar structure and a low melting point (m_p, <-20 °C), which has not been achieved by reported counterparts bearing linear alkyl chains. In solvent-free states, microphase segregation occurs where polar and apolar domains arrange bicontinuously as proved by molecular dynamics (MD) simulations. When dispersed in water, self-stabilized giant aggregates formed with ultrahigh colloidal stability (up to years). MD simulations provide clues of discrete bicontinuous phases within the giant aggregates. These newly discovered self-assemblies provide a heterogeneous reservoir that can accommodate guest molecules including the highly apolar fullerene C₆₀, paving the way for a wide range of potential applications.

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.

Photon Upconversion in TTA-Inducing Ionic Liquids: Pinpointing the Role of IL Nanostructured Media Using MD Simulations

The use of task-specific chromophoric ionic liquids as energy transfer media in triplet-triplet annihilation photon upconversion (TTA-UC) processes has produced several examples of systems with signifficantly enhanced performances. In this work, we use molecular dynamics simulations to probe the relation between the nanostructure of chromophoric ionic liquids and their ability to achieve high TTA-UC quantum yields. The existing atomistic and systematic force fields commonly used to model different ionic liquids are extended to include substituted anthracene moieties, thus allowing the modeling of several chromophoric ionic liquids. The simulation results show that the polar network of the ionic liquids can orient the anthracene moieties within the nonpolar domains preventing direct contacts between them but allowing orientations at the optimal distance for triplet energy migration.

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.

Influence of Hydration on the Structure of Reline Deep Eutectic Solvent: A Molecular Dynamics Study

In this article, we have performed an all-atom molecular dynamics simulation study to investigate the influence of water on the molecular level arrangement of reline deep eutectic solvent for different hydration levels ranging from 3.4 to 58.1 wt % of water and complemented the observations of recently measured neutron scattering experimental data. This study is particularly important because water is being introduced as a second hydrogen bond donor/acceptor in reline, wherein the structure is primarily governed by hydrogen bonding and electrostatic interactions. We have analyzed the simulated X-ray scattering structure functions, their partial components, and hydrogen bonding interactions to understand the effects of water on various intermolecular interactions in reline-water mixtures. It is observed that at lower hydration level, reline structure is qualitatively retained. At higher hydration level, most water molecules preferentially solvate chloride anions and ammonium group of choline cations mostly impacting choline-choline, choline-chloride, and chloride-chloride interactions. The present study reveals that at and above 41 wt % of water, the molecular arrangement of reline drastically changes and set to transition from reline to an aqueous solution of reline components with further increase in the hydration level. Hydrogen bond analysis reveals the presence of strong chloride-water H-bonding interaction, which gradually replaces choline-chloride and urea-chloride hydrogen bondings as the hydration level in the mixture increases.

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.

Effect of molecular solvents of varying polarity on the self-assembly of 1-n-dodecyl-3-methylimidazolium octylsulfate ionic liquid

Large-scale molecular dynamics simulations consisting of more than 88,000–106,000 atoms for approximately 250 ns (including equilibration and production) were conducted to assess the effect of polar, nonpolar and amphiphilic molecular solvents on the nanoscale structuring of 1-[Formula: see text]-dodecyl-3-methylimidazolium [C[Formula: see text]mim] octylsulfate [C8SO4] ionic liquid (IL). Water [H2O], [Formula: see text]-octane [C8H[Formula: see text]] and 1-octanol [C8H[Formula: see text]OH] are employed as examples of polar, nonpolar, and amphiphilic molecules, respectively. The results indicate that each of these molecular solvents modify the nanosegregation behavior of the ionic liquid in a unique way. Water induces a high order of structuring of the ionic liquid as indicated by extremely high nematic order parameter for the system. In addition, the morphology of the neat ionic liquid is transformed from layer-like to that of bilayer-like in which the polar and nonpolar domains alternate. The presence of water also causes the stretching of the nonpolar domain, thus, increasing its size. At the concentration examined in this work, [Formula: see text]-octane is found to be only partially miscible with the ionic liquid. The polar network is maintained; however, the continuous cationic nonpolar domain is split into multiple domains. [Formula: see text]-octane is accommodated in the ionic liquid nonpolar domain. Similarly, the amphiphilicity of 1-octanol leads to an increase in the number of cationic as well as anionic domains. The overall nonpolar domain length, however, remains nearly identical to that found for the pure ionic liquid. Additional characterization of structural features of the three systems is discussed in terms of one-dimensional number densities, nematic order parameters for the overall systems and their components and structure factors.

Ionic Liquids Treated within the Grand Canonical Adaptive Resolution Molecular Dynamics Technique

We use the Grand Canonical Adaptive Resolution Molecular Dynamics Technique (GC-AdResS) to examine the essential degrees of freedom necessary for reproducing the structural properties of the imidazolium class of ionic liquids. In this technique, the atomistic details are treated as an open sub-region of the system while the surrounding environment is modelled as a generic coarse-grained model. We systematically characterize the spatial quantities such as intramolecular, intermolecular radial distribution functions, other structural and orientational properties of ILs. The spatial quantities computed in an open sub-region of the system are in excellent agreement with the equivalent quantities calculated in a full atomistic simulation, suggesting that the atomistic degrees of freedom outside the sub-region are negligible. The size of the sub-region considered in this study is 2 nm, which is essentially the size of a few ions. Insight from the study suggests that a higher degree of spatial locality seems to play a crucial role in characterizing the properties of imidazolium based ionic liquids.