Mesoscale Organization and Dynamics in Binary Ionic Liquid Mixtures

Chain-length dependent organisation in mixtures of hydrogenous and fluorous ionic liquids

As part of an ongoing study of the structure and properties of mixtures of ionic liquids in which one component has a hydrocarbon chain and the other a semiperfluorocarbon chain, we now report a study of the mixtures [C8MIM]1-x[C10MIM-F17]x[Tf2N], [C10MIM]1-x[C8MIM-F13]x[Tf2N] and [C10MIM]1-x[C10MIM-F17]x[Tf2N], where [C8MIM][Tf2N] is 1-methyl-3-octylimidazolium bis(trifluoromethylsulfonyl)imide, [C10MIM][Tf2N] is 1-decyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [C8MIM-F13][Tf2N] is 1-(1H,1H,2H,2H-perfluorooctyl)-3-methylimidizolium bis(trifluoromethylsulfonyl)imide and [C10MIM-F17][Tf2N] is 1-(1H,1H,2H,2H-perfluorodecyl)-3-methylimidizolium bis(trifluoromethylsulfonyl)imide. The mixtures were investigated using small-angle X-ray (SAXS) and neutron (SANS) scattering complemented by molecular dynamics simulations (with viscosity and surface tension measurements also possible for the mixtures [C10MIM]1-x[C8MIM-F13]x[Tf2N]). Unlike previous studies of [C8MIM]1-x[C8MIM-F13]x[Tf2N], where no strong evidence of alkyl/fluoroalkyl chain segregation or triphilic behaviour was seen (Elstone et al., J. Phys. Chem. B, 2023, 127, 7394-7407), these new mixtures show the formation of small aggregates of varying sizes of each component, even though all were co-miscible across the full range of compositions. Thus, while a clear polar non-polar peak (PNPP) was observed at large or small values of x, at intermediate compositions the small-angle neutron scattering at low q was dominated by scattering from these small aggregates, while at other compositions, there was little or no evidence of the PNPP. The origins of this behaviour are discussed in terms of inter-chain interactions.

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.

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.

Understanding the Liquid Structure in Mixtures of Ionic Liquids with Semiperfluoroalkyl or Alkyl Chains

By mixing ionic liquids (ILs), it is possible to fine-tune their bulk and interfacial structure. This alters their physical properties and solvation behavior and is a simple way to prepare a collection of ILs whose properties can be tuned to optimize a specific application. In this study, mixtures of perfluorinated and alkylated ILs have been prepared, and links between composition, properties, and nanostructure have been investigated. These different classes of ILs vary substantially in the flexibility and polarizability of their chains. Thus, a range of useful structural and physical property variations are accessible through mixing that will expand the library of IL mixtures available in an area that to this point has received relatively little attention. In the experiments presented herein, the physical properties and bulk structure of mixtures of 1-methyl-3-octylimidazolium bis(trifluoromethylsulfonyl)imide [C8MIM][Tf2N] and 1-(1H,1H,2H,2H-perfluorooctyl)-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [C8MIM-F13][Tf2N] have been prepared. The bulk liquid structure was investigated using a combination of small-angle X-ray and neutron scattering (SAXS and SANS, respectively) experiments in combination with atomistic molecular dynamics simulations and the measurement of density and viscosity. We observed that the addition of [C8MIM-F13][Tf2N] to [C8MIM][Tf2N] causes changes in the nanostructure of the IL mixtures that are dependent on composition so that variation in the characteristic short-range correlations is observed as a function of composition. Thus, while the length scales associated with the apolar regions (polar non-polar peak─PNPP) increase with the proportion of [C8MIM-F13][Tf2N] in the mixtures, perhaps surprisingly given the greater volume of the fluorocarbon chains, the length scale of the charge-ordering peak decreases. Interestingly, consideration of the contact peak shows that its origins are both in the direct anion···cation contact length scale and the nature (and hence volume) of the chains appended to the imidazolium cation.

Phase coexistence in [C22/C1MIm]+[NO3]- ionic-liquid mixtures and first-order phase transitions from homogeneous liquid to smectic B by varying the cation ratio

We perform molecular dynamics simulations to investigate the transition processes of [C22/C1MIm]+[NO3]- binary mixtures by varying the cation ratio of C22 to C1 at a fixed temperature of 400 K. The cation ratio is tuned by ranging C22 percentage from 0% to 100% with a fixed number of 4096 total simulated ion pairs. Our simulated-annealing results indicate that, at 400 K, pure C1 is a homogeneous liquid whilst pure C22 is an ionic liquid crystal (ILC) of smectic-B (SmB) type. With increasing C22 percentage, the system goes through a first-order phase transition from homogeneous liquid to nano-fragment liquid in the range from 15% to 17.5%, during which some of the individual cationic alkyl side chains locally aggregate to form small bundles "floating" in the polar "solvent" composed of anions and cationic head groups. Although the side chains in each bundle are parallelly aligned, the bundles distribute randomly without a global orientation. As the C22 percentage further increases, another first-order phase transition occurs to bring the system into the SmB ILC phase. Particularly, when the C22 percentage is in the range from 45% to 50%, the SmB phase coexists with the liquid phase containing both individual and bundled alkyl side chains.

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.

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

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

Mixing Ionic Liquids Affects the Kinetics and Thermodynamics of the Oxygen/Superoxide Redox Couple in the Context of Oxygen Sensing

The electrochemical oxygen reduction reaction is vital for applications such as fuel cells, metal air batteries and for oxygen gas sensing. Oxygen undergoes a 1-electron reduction process in dry ionic liquids (ILs) to form the electrogenerated superoxide ion that is solvated and stabilized by IL cations. In this work, the oxygen/superoxide (O₂/O₂ ^•-) redox couple has been used to understand the effect of mixing ILs with different cations in the context of developing designer electrolytes for oxygen sensing, by employing cyclic voltammetry at both gold and platinum electrodes. Different cations with a range of sizes, geometries and aromatic/aliphatic character were studied with a common bis(trifluoromethylsulfonyl)imide ([NTf₂]^-) anion. Diethylmethylsulfonium ([S_2,2,1]⁺), N-butyl-N-methylpyrrolidinum ([C₄mpyrr]⁺) and tetradecyltrihexylphosphonium ([P_14,6,6,6]⁺) cations were mixed with a common 1-butyl-3-methylimidazolium ([C₄mim]⁺) cation at mole fractions (x) of [C₄mim]⁺ of 0, 0.2, 0.4, 0.6, 0.8, and 1. Both the redox kinetics and thermodynamics were found to be highly dependent on the cation structure and the electrode material used. Large deviations from "ideal" mixtures were observed for mixtures of [C₄mim][NTf₂] with [C₄mpyrr][NTf₂] on gold electrodes, suggesting a much higher amount of [C₄mim]⁺ ions near the electrode surface despite the large excess of [C₄mpyrr]⁺ in the bulk. The electrical double layer structure was probed for a mixture of [C₄mim]_0.2[C₄mpyrr]_0.8[NTf₂] using atomic force microscopy measurements on Au, revealing that the first layer was more like [C₄mim][NTf₂] than [C₄mpyrr][NTf₂]. Unusually fast kinetics for O₂/O₂ ^•- in mixtures of [C₄mim]⁺ with [P_14,6,6,6]⁺ were also observed in the electrochemistry results, which warrants further follow-up studies to elucidate this promising behavior. Overall, it is important to understand the effect on the kinetic and thermodynamic properties of electrochemical reactions when mixing solvents, to aid in the creation of designer electrolytes with favorable properties for their intended application.

Small-angle neutron scattering from mixtures of long- and short-chain 3-alkyl-1-methyl imidazolium bistriflimides

The preparation of mixtures of ionic liquids (ILs) represents an attractive strategy to tune their properties, an important aspect of which is to understand how the structure of the bulk varies with composition. In this study, small-angle neutron scattering (SANS) was used to probe mixtures of methylimidazolium-based ionic liquids [C_nmim][Tf₂N] with [C₂mim][Tf₂N]) (n = 4, 6, 8 and 10) and of [C_mmim][Tf₂N] with [C₁₂mim][Tf₂N] (m = 2, 4, 6 and 8). Mixtures were prepared in both contrasts, which is to say that one component would be fully hydrogenated while the other was fully deuterated, and vice versa. Data were fitted using a range of appropriate models, of which the Teubner-Strey model provided most useful information and the pure materials showed a nascent Polar Non-polar Peak (PNPP) for n = 6, which became more evident as n increased. In the mixtures [C_nmim]_x[C₂mim]_1-x[Tf₂N], the PNPP was evident for n = 10 and 8, nascent for n = 6 and absent for n = 4, with percolation showing a very strong dependence on the chain length of the added IL, [C_nmim][Tf₂N]. In contrast, while the ability of [C₁₂mim][Tf₂N] to form percolated structures was damped when mixed with [C_mmim][Tf₂N], as m increased from 2 to 6, this effect was less strong. However, data obtained for mixtures of [C₁₂mim][Tf₂N] and [C₈mim][Tf₂N], both of which percolate as pure materials, did not fit easily in any of the models applied to the previous systems and gave results that depended on the contrast used. Complementary small-angle X-ray scattering (SAXS) data, however, showed the expected evolution and behaviour of the PNPP, COP and CP, revealing that the unexpected observations were due to an adventitious matching out of isotopic contrasts. As well as revealing details of the structures of these IL mixtures, the results also point to complementary strategies for generating bulk percolated structures as a function of cation chain length.

Computer Simulations of Deep Eutectic Solvents: Challenges, Solutions, and Perspectives

Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications, such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can predict and reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.

Hydrogen Bond Redistribution Effects in Mixtures of Protic Ionic Liquids Sharing the Same Cation: Non-ideal Mixing with Large Negative Mixing Enthalpies

We report a joint experimental and theoretical study characterising the hydrogen bond (HB) redistribution in mixtures of two different protic ionic liquids (PILs) sharing the same cation: triethylammonium-methanesulfonate ([TEA][OMs]) and...

Design of a homologous series of molecular glassformers

We design and synthesize a set of homologous organic molecules by taking advantage of facile and tailorable Suzuki cross coupling reactions to produce triarylbenzene derivatives. By adjusting the number and the arrangement of conjugated rings, the identity of heteroatoms, lengths of fluorinated alkyl chains, and other interaction parameters, we create a library of glassformers with a wide range of properties. Measurements of the glass transition temperature (T_g) show a power-law relationship between T_g and molecular weight (M_W), with of the molecules, with an exponent of 0.3 ± 0.1, for T_g values spanning a range of 300-450 K. The trends in indices of refraction and expansion coefficients indicate a general increase in the glass density with M_W, consistent with the trends observed in T_g variations. A notable exception to these trends was observed with the addition of alkyl and fluorinated alkyl groups, which significantly reduced T_g and increased the dynamical fragility (which is otherwise insensitive to M_W). This is an indication of reduced density and increased packing frustrations in these systems, which is also corroborated by the observations of the decreasing index of refraction with increasing length of these groups. These data were used to launch a new database for glassforming materials, glass.apps.sas.upenn.edu.

In Silico Elucidation of Molecular Picture of Water-Choline Chloride Mixture

Choline chloride (ChCl) is a component of several deep eutectic solvents (DESs) having numerous applications. Recent studies have reported manifold promising use of aqueous choline chloride solution as an alternative to DES, where water plays the role of the hydrogen-bond donor. The characteristic physical properties of the DESs and aqueous DES originate from the "inter-" and intraspecies hydrogen-bond network formed by the constituents. However, a detailed molecular-level picture of choline chloride and water mixture is largely lacking in the literature. This motivates us to carry out extensive all-atom molecular dynamics simulations of the ChCl-water mixture of varying compositions. Our analyses clearly show an overall increase in the interspecies association with an increase in ChCl concentration. At higher concentrations, the trimethylammonium groups of choline are stabilized by a nonpolar interaction, whereas the hydroxyl groups preferentially interact with water. Chloride ions are found to be involved in two types of interactions: one where chloride ions intercalate two or more choline cations, and the other one where they are surrounded by five to six water molecules forming solvated chloride ions. However, the relative fractions of these two types of associations depend on the concentration of ChCl in the mixture. Another important structural aspect is the disruption of the hydrogen-bonded water network due to the presence of both choline cations and chloride ions. However, chloride ions participate to partially restore the tetrahedral arrangement of partners around water molecules.

Binary mixtures of ionic liquids: Ideal, non-ideal, or quasi-ideal?

The mixing of ILs provides an opportunity for fine tuning the physiochemical properties of ILs for various applications. However, a suitable mixture having desired properties can only be designed when the physiochemical properties of the mixtures of ILs along with their spectroscopic properties are well understood. With an aim to achieve this objective, three different mixtures with a common anion, namely, [C₂C₁im][C₄C₁im][NTf₂], [C₃C₁pyr][C₄C₁pyr][NTf₂], and [C₃C₁im][C₃C₁pyr][NTf₂], have been investigated in the current study. Investigations have been carried out at the macroscopic level by observing the thermophysical properties, such as molar volume and thermal expansion coefficient, and at the microscopic level with time-resolved fluorescence measurements and the pulse field gradient nuclear magnetic resonance (NMR) technique. The results obtained from the thermophysical study have indicated that excess molar volume for imidazolium-based IL-IL mixtures may be linked to the free volume created by the alkyl chain of the imidazolium cation whereas for the mixture of pyrrolidinium ILs, lowering of density can give rise to free volume. Analysis of time-resolved fluorescence anisotropy data has provided clear evidence in favor of the presence of free volume in the binary mixture of ILs. NMR studies have also supported the fluorescence anisotropy data. The outcome of the present investigation reveals that the mixtures show appreciable deviation from ideal behavior and the deviation from the ideal behavior is caused due to the generation of free volume in the resultant mixture, describing these IL mixtures as quasi-ideal rather than ideal or non-ideal.

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.

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.

Thorough studies of tricyanomethanide-based ionic liquids - the influence of alkyl chain length of the cation

The glassy, supercooled, and normal liquid states of the 1-alkyl-3-methylimidazolium tricyanomethanide series [CnC1im][TCM] (n = 2, 4, 6, 8, and 16) were investigated by dielectric and mechanical (rheological) experiments supplemented by X-ray diffraction. The conductivity relaxation was found to be accompanied by a pronounced secondary relaxation. However, based on ambient and high-pressure results as well as the coupling model, we assumed that the latter one can not be classified as Johari-Goldstein relaxation. Moreover, the studies on the nanoscale organization of ionic liquids indicated that 1-alkyl-3-methylimidazolium tricyanomethanide ILs begin to form nanoscale aggregates when the alkyl chain of the cation has six carbon atoms.

Macroscopic Differentiators for Microscopic Structural Nonideality in Binary Ionic Liquid Mixtures

Combining two ionic liquids to form a binary ionic liquid mixture is a simple yet effective strategy to not only expand the number of ionic liquids but also precisely control various physicochemical properties of resultant ionic liquid mixtures. From a fundamental thermodynamic point of view, it is not entirely clear whether such mixtures can be classified as ideal solutions. Given a large number of binary ionic liquid mixtures that emerge, the ability to predict the presence of nonideality in such mixtures a priori without the need for experimentation or molecular simulation-based calculations is immensely valuable for their rational design. In this research report, we demonstrate that the difference in the molar volumes (ΔV) of the pure ionic liquids and the difference in the hydrogen-bonding ability of anions (Δβ) are the primary determinants of nonideal behavior of binary ionic liquid mixtures containing a common cation and two anions. Our conclusion is derived from a comparison of microscopic structural properties expressed in terms of radial, spatial, and angular distributions for binary mixtures and those of the corresponding pure ionic liquids. Molecular dynamics simulations of 16 binary ionic liquid mixtures, containing a common cation 1-n-butyl-3-methylimidazolium [C₄mim]⁺ and combinations of (less basic) fluorinated {trifluoromethylacetate [TFA]^-, trifluoromethanesulfonate [TFS]^-, bis(trifluoromethanesulfonyl)imide [NTf₂]^-, and tris (pentafluoroethyl) trifluorophosphate [eFAP]^-} versus (more basic) nonfluorinated {chloride Cl^-, acetate [OAC]^-, methylsulfate [MeSO₄]^-, and dimethylphosphate [Me₂PO₄]^-} anions, were conducted. The large number of binary ionic liquid mixtures examined here enabled us to span a broad range of ΔV and Δβ values. The results indicate that binary mixtures of two ionic liquids for which ΔV > 60 cm³/mol and Δβ > 0.4 are expected to be microscopically nonideal. On the other hand, ΔV < 60 cm³/mol and Δβ < 0.4 will lead to molecular structures that are not differentiated from those of their pure ionic liquid counterparts.

Elucidating the effect of the ionic liquid type and alkyl chain length on the stability of ionic liquid-iron porphyrin complexes

The present study is motivated by the long-term objective of understanding how ionic liquids are biodegraded by cytochrome P450, which contains iron porphyrin (FeP) serving as the catalytic center. To this end, the current study is designed to elucidate the impact of types and conformations of ionic liquids on the binding energy with FeP, the key interactions that stabilize the ionic liquid-FeP complex, and how the electron uptake ability of FeP is altered in the presence of ionic liquids. Four classes of ionic liquids are considered: 1-alkyl-3-methylimidazolium, 1-alkyl-pyridinium, 1-alkylsulfonium, and N-methyl-N-alkylpyrrolidinium. The influence of linear alkyl chains of ethyl, butyl, hexyl, octyl, and decyl is examined on the favorable binding modes with FeP, considering two widely different conformations: tail up and tail down with respect to FeP. Electronic structure calculations are performed at the M06 level of theory with the 6-31G(d,p) basis set for C, H, and N atoms, while the Lanl2DZ basis set is employed for Fe. Donor-acceptor interactions contributing to the binding of ionic liquids to FeP are unraveled through the natural bond orbital analysis. The results from this study indicate that the binding energies are dependent not only on the class of ionic liquids but also on the conformations presented to FeP. The propensity of FeP to acquire an electron is significantly enhanced in the presence of ionic liquid cations, irrespective of the type and the alkyl chain length.

Decelerated Liquid Dynamics Induced by Component-Dependent Supercooling in Hydrogen and Deuterium Quantum Mixtures

Isotopic mixtures of p-H₂ and o-D₂ molecules have been an attractive binary system because they include two kinds of purely isotopic molecules which possess the same electronic potential but the twice different mass inducing differently pronounced nuclear quantum effects (NQEs). Accessing details of structures and dynamics in such quantum mixtures combining complex molecular dynamics with NQEs of different strengths remains a challenging problem. Taking advantage of the nonempirical molecular dynamics method which describes p-H₂ and o-D₂ molecules, we found that the liquid dynamics slows down at a specific mixing ratio, which can be connected to the observed anomalous slowdown of crystallization in the quantum mixtures. We attributed the decelerated dynamics to the component-dependent supercooling of p-H₂ taking place in the mixtures, demonstrating that there is an optimal mixing ratio to hinder crystallization. The obtained physical insights will help in experimentally controlling and achieving unknown quantum mixtures including superfluid.

Anions as Dynamic Probes for Ionic Liquid Mixtures

Ionic liquid (IL) mixtures have been proposed as a viable alternative to rationally fine-tune the physicochemical properties of ILs for a variety of applications. The understanding of the effects of mixing ILs on the properties of the mixtures is however only in the very early stages. Two series of ionic liquid mixtures, based on the 1-ethyl-3-methylimidazolium and 1-dodecyl-3-methylimidazolium cations, and having a common anion (tetrafluoroborate or bis(trifluoromethylsulfonyl)imide), have been prepared and deeply characterized via multiple NMR techniques. Diffusion and relaxation methods combined with 2D ion-ion correlation (nuclear Overhauser enhancement) experiments have been used for a better understanding of the interplay between dynamics and structure of IL mixtures. A crucial role of the anion in driving the mixture's behavior emerged, making them important "dynamic probes" for gaining information of the polar and nonpolar regions of ionic liquids and their mixtures.