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

Structure of Novel Phosphonium-Based Ionic Liquids with S and O Substitutions from Experiments and a Mixed Quantum-Classical Approach

This article presents experimental characterization information and synchrotron X-ray scattering measurements on a set of novel O- and S-substituted phosphonium-based ionic liquids (ILs) all coupled with the bis(fluorosulfonyl)imide (FSI-) anion. The ILs include the ethoxyethyltriethylphosphonium (P222(2O2)+) and triethyl[2-(ethylthio)ethyl]phosphonium (P222(2S2)+) cations, and we contrast results on these with those for unsubstituted triethylpentylphosphonium (P2225+). The article also introduces a physics-based protocol that combines classical force field studies on larger simulation boxes with classical and first-principles studies on smaller boxes. The method produces significantly improved S(q) functions in the regime which in prior publications we have associated with inter- and intraionic adjacency correlations. By understanding which shorter-range structural changes improve S(q) in the q-regime of interest, we are also able to pinpoint specific deficiencies in the classical force field model. The approach we take should be quite general and could help study other complex liquids on different length scales.

Chlorine gas and anion radical reactivity in molten salts and the link to chlorobasicity

Next-generation nuclear power plants may include exciting novel designs in which molten salts are the coolant or a combination of the coolant and fuel. Whereas it is straightforward to see why having a low volatility coolant can be advantageous for safety, much is not understood about the production of volatile halogen gases as a result of radiation and even less is known about the distribution of these species at and away from interfaces. Using first principles molecular dynamics simulations, we investigate the product of the disproportionation reaction between chlorine anion radicals (nominally Cl2˙-) in the bulk and slab configurations. We find that the product depends on the chlorobasicity of the medium. For example, in ZnCl2, Cl2 forms, but in a eutectic mixture of LiCl and KCl, Cl3- is formed as a product. We also find that Cl3- prefers to form at the vapor interface and this may have implications for corrosion and reactivity. Furthermore, the mechanisms of the mobility of Cl2 and Cl3- are radically different, the first one being vehicular and the second Grotthus-like. Chlorobasicity is linked to the electronic structure of the host melt; ZnCl2 forms extended networks along which metal ions and anionic counterions have significant electronic orbital overlap forming long, linear, molecular-like constructs; the opposite is true for the alkali metal eutectic salt.

Synthesis and applications of ionic liquids for chromatographic analysis

Ionic liquids (ILs) have emerged as more desirable liquids than conventional solvents for chemistry, material science, engineering and environmental science. The scientific literature reveals an exponential increase in the number of research projects aimed at exploring the chromatographic features of ionic liquids. The review provides sound scientific data to examine the structural characteristics of ionic liquids that make them ideal for use in chromatography. This contribution is distinctive since it integrates the synthesis, benefits, drawbacks, and possible uses of ionic liquids in several chromatographic separation processes. Keeping the cation the same, the introduction of different anions is also possible, and this strategy leads to the synthesis of a series of different ionic liquids with varying properties. A detailed probe is given on the influence of ionic liquid structure and properties on their chromatographic behavior, both as stationary phase and mobile phase and/or mobile phase additives. Ionic liquid based immobilized stationary phases and their analyte retention mechanisms (hydrogen bonding, electrostatic forces of attraction, π-π stacking, ion exchange, and hydrophilic interactions, etc.) are critically discussed. Finally, a thorough analysis of the literature suggests that IL-based stationary phases may undergo multi-mode and more flexible retention mechanisms. Their dual polarity can facilitate interaction with both polar and non-polar compounds. Similarly, using IL as a mobile phase can offer more pragmatic and sustainable options for enantiomer separation.

Short-time collective dynamics of an ionic liquid: A computer simulation study with non-polarizable and polarizable models, and ab initio molecular dynamics

Molecular dynamics (MD) simulation is used to study the intermolecular dynamics in the THz frequency range of the ionic liquid 1-ethyl-3-methylimidazolium bis(fluorosulfonyl)imide, [C2C1im][FSI]. Non-polarizable and polarizable models for classical MD simulation are compared using as quality criteria ab initio molecular dynamics (AIMD) and experimental data from far-infrared (FIR) spectroscopy and previously published data of inelastic x-ray scattering (IXS). According to data from IXS spectroscopy, incorporating polarization in the classical MD simulation has relatively little effect on the dispersion curve (excitation frequency vs wavevector) for longitudinal acoustic modes. When the AIMD simulation is used as a reference, the polarizable model leads the time correlation functions of velocity, mass, and charge currents to relax abnormally quickly. The charge current spectra from the AIMD simulation and the non-polarizable model agree with the experimental FIR spectrum, while the polarizable model gives an excessively broad band. When compared to the non-polarizable model, the polarizable model does improve the calculation of transport coefficients (diffusion coefficient, viscosity, and conductivity); however, it yields overdamped short-time collective dynamics.

Vibrational spectroscopy of dynamic interactions and relaxation processes in ionic systems

In recent years, there has been a boom in the publication of research results on ionic (electrolyte) systems used in the creation of electrochemical energy systems (batteries, supercapacitors, etc.). One of the alternative ways to obtain dynamic information in ionic systems is spectroscopic experiments, and the most valuable information on the behavior of condensed systems at picosecond and adjacent time intervals is provided by vibrational spectroscopy. The source of this information is the analysis of the shape and width of the lines in the vibrational spectra of ion systems. This review aims to systematize the theoretical and methodological aspects of the study of the dynamics of ion systems by vibrational spectroscopy.

Spiers Memorial Lecture: From cold to hot, the structure and structural dynamics of dense ionic fluids

The structure of ionic liquids (ILs), which a decade or two ago was the subject of polite but heated debate, is now much better understood. This has opened opportunities to ask more sophisticated questions about the role of structure in transport, the structure of systems with ions that are not prototypical, and the similarity between ILs and other dense ionic fluids such as the high-temperature inorganic molten salts; let alone the fact that new areas of research have emerged that sprung from our collective understanding of the structure of ILs such as the deep eutectic solvents, the polymerized ionic liquids, and the zwitterionic liquids. Yet, our understanding of the structure of prototypical ILs may not be as complete as we think it to be, given that recent experiments appear to show that in some cases there may be more than one liquid phase resulting in liquid-liquid (L-L) phase transitions. This article presents a perspective on what we think are key topics related to the structure and structural dynamics of ILs and to some extent high-temperature molten salts.

Unusual dynamics of tetrahedral liquids caused by the competition between dynamic heterogeneity and structural heterogeneity

Tetrahedral liquids exhibit intriguing thermodynamic and transport properties because of the various ways tetrahedra can be packed and connected. Recently, an unusual temperature dependence of the stretching exponent β in a model tetrahedral liquid ZnCl2 from Tm + 85 K to Tm + 35 K has been reported using neutron-spin echo spectroscopy. This discovery stands in sharp contrast to other glass-forming liquids. In this study, we conducted neural network force field driven molecular dynamic simulations of ZnCl2. We found a non-monotonic temperature dependence of β from liquid to supercooled liquid temperatures. Further structural decomposition and dynamic analysis suggest that this unusual dynamic behavior is a result of the competition between the decrease in the diversity of tetrahedra motifs (structural heterogeneity) and the increase in glassy dynamic heterogeneity. This result may contribute to new understandings of the structural relaxation of other network liquids.

Structure and dynamics of Li1.24K0.76CO3 molten carbonate electrolyte from molecular simulations with explicit polarization

Molten carbonate electrolysis cells represent a key technology for harnessing surplus energy from renewable sources and converting it into gaseous energy carriers. To optimize their efficiency, a comprehensive understanding of each step in the operational process is essential. Here, we focus on the electrolyte of choice in molten carbonate cells: the Li1.24K0.76CO3 melt. Utilizing molecular dynamics with explicit polarization, we demonstrate that the structure of this molten mixture is characterized by a dense network of lithium-carbonate complexes, with K+ ions loosely embedded within this network. This structural insight enables us to rationalize from an atomistic perspective the conductivity trends observed experimentally in molten carbonates. Moreover, our work highlights the importance of including polarization for the simulations of dense liquid carbonates. It also acts as a foundational step towards more advanced theoretical studies for elucidating the role of the electrolyte in these devices.

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.

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.

Phase behaviour and aggregate structures of the surface-active ionic liquid [BMIm][AOT] in water

HYPOTHESIS

The surface-active ionic liquid, 1-butyl-3-methylimidazolium 1,4-bis-2-ethylhexylsulfosuccinate ([BMIm][AOT]), has a sponge-like bulk nanostructure consisting of percolating polar and apolar domains formed by the ion charge groups and alkyl chains, respectively. We hypothesise that added water will swell the polar domains and change the liquid nanostructure.

EXPERIMENTS

Small angle X-ray scattering (SAXS), small angle neutron scattering (SANS) and polarizing optical microscopy (POM) were used to investigate the nanostructure of [BMIm][AOT] as a function of water content. Differential scanning calorimetry (DSC) was employed to probe the thermal transitions of [BMIm][AOT]-water mixtures and the mobility of water molecules.

FINDINGS

SAXS, SANS and POM show that at lower water contents, [BMIm][AOT]-water mixtures have a sponge-like nanostructure similar to the pure SAIL, at medium water contents a lamellar phase forms, and at high water contents vesicles form. DSC results reveal that water molecules are supercooled in the lamellar phase. For the first time, results reveal a series of transitions from inverse sponge, to lamellar then to vesicles, for [BMIm][AOT] upon dilution with water.

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.

Effect of ion structure on the physicochemical properties and gas absorption of surface active ionic liquids

Surface active ionic liquids (SAILs) combine useful characteristics of both ionic liquids (ILs) and surfactants, hence are promising candidates for a wide range of applications. However, the effect of SAIL ionic structures on their physicochemical properties remains unclear, which limits their uptake. To address this knowledge gap, in this work we investigated the density, viscosity, surface tension, and corresponding critical micelle concentration in water, as well as gas absorption of SAILs with a variety of cation and anion structures. SAILs containing anions with linear alkyl chains have smaller molar volumes than those with branched alkyl chains, because linear alkyl chains are interdigitated to a greater extent, leading to more compact packing. This interdigitation also results in SAILs being about two orders of magnitude more viscous than comparable conventional ILs. SAILs at the liquid-air interface orient alkyl chains towards the air, leading to low surface tensions closer to n-alkanes than conventional ILs. Critical temperatures of about 900 K could be estimated for all SAILs from their surface tensions. When dissolved in water, SAILs adsorb at the liquid-air interface and lower the surface tension, like conventional surfactants in water, after which micelles form. Molecular simulations show that the micelles are spherical and that lower critical micelle concentrations correspond to the formation of aggregates with a larger number of ion pairs. CO₂ and N₂ absorption capacities are examined and we conclude that ionic liquids with larger non-polar domains absorb larger quantities of both gases.

Molecular Modelling of Ionic Liquids: Situations When Charge Scaling Seems Insufficient

Charge scaling as an effective solution to the experiment-computation disagreement in molecular modelling of ionic liquids (ILs) could bring the computational results close to the experimental reference for various thermodynamic properties. According to the large-scale benchmark calculations of mass density, solvation, and water-ILs transfer-free energies in our series of papers, the charge-scaling factor of 0.8 serves as a near-optimal option generally applicable to most ILs, although a system-dependent parameter adjustment could be attempted for further improved performance. However, there are situations in which such a charge-scaling treatment would fail. Namely, charge scaling cannot really affect the simulation outcome, or minimally perturbs the results that are still far from the experimental value. In such situations, the vdW radius as an additional adjustable parameter is commonly tuned to minimize the experiment-calculation deviation. In the current work, considering two ILs from the quinuclidinium family, we investigate the impacts of this vdW-scaling treatment on the mass density and the solvation/partition thermodynamics in a fashion similar to our previous charge-scaling works, i.e., scanning the vdW-scaling factor and computing physical properties under these parameter sets. It is observed that the mass density exhibits a linear response to the vdW-scaling factor with slopes close to -1.8 g/mL. By further investigating a set of physiochemically relevant temperatures between 288 K and 348 K, we confirm the robustness of the vdW-scaling treatment in the estimation of bulk properties. The best vdW-scaling parameter for mass density would worsen the computation of solvation/partition thermodynamics, and a marginal decrease in the vdW-scaling factor is considered as an intermediate option balancing the reproductions of bulk properties and solvation thermodynamics. These observations could be understood in a way similar to the charge-scaling situation. i.e., overfitting some properties (e.g., mass density) would degrade the accuracy of the other properties (e.g., solvation free energies). Following this principle, the general guideline for applying this vdW-tuning protocol is by using values between the density-derived choice and the solvation/partition-derived solution. The charge and current vdW scaling treatments cover commonly encountered ILs, completing the protocol for accurate modelling of ILs with fixed-charge force fields.

Atomistic simulation framework for molten salt vapor–liquid equilibrium prediction and its application to NaCl

Knowledge of the vapor–liquid equilibrium (VLE) properties of molten salts is important in the design of thermal energy storage systems for solar power and nuclear energy production applications. The high temperatures involved make their experimental determination problematic, and the development of both macroscopic thermodynamic correlations and predictive molecular-based methodologies are complicated by the requirement to appropriately incorporate the chemically reacting vapor-phase species. We derive a general thermodynamic-based atomistic simulation framework for molten salt VLE prediction and show its application to NaCl. Its input quantities are temperature-dependent ideal-gas free energy data for the vapor phase reactions and density and residual chemical potential data for the liquid. If these are not available experimentally, the former may be predicted using standard electronic structure software, and the latter may be predicted by means of classical atomistic simulation methodology. The framework predicts the temperature dependence of vapor pressure, coexisting phase densities, vapor phase composition, and vaporization enthalpy. It also predicts the concentrations of vapor phase species present in minor amounts (such as the free ions), quantities that are extremely difficult to measure experimentally. We furthermore use the results to obtain an approximation to the complete VLE binodal dome and the critical properties. We verify the framework for molten NaCl, for which experimentally based density and chemical potential data are available in the literature. We then apply it to the analysis of NaCl simulation data for two commonly used atomistic force fields. The framework can be readily extended to molten salt mixtures and to ionic liquids.

Chemistry Dissolved in Ionic Liquids. A Theoretical Perspective

The theoretical treatment of ionic liquids must focus now on more realistic models while at the same time keeping an accurate methodology when following recent ionic liquids research trends or allowing predictability to come to the foreground. In this Perspective, we summarize in three cases of advanced ionic liquid research what methodological progress has been made and point out difficulties that need to be overcome. As particular examples to discuss we choose reactions, chirality, and radicals in ionic liquids. All these topics have in common that an explicit or accurate treatment of the electronic structure and/or intermolecular interactions is required (accurate methodology), while at the same time system size and complexity as well as simulation time (realistic model) play an important role and must be covered as well.

Fine-Tuning the Polarizable CL&Pol Force Field for the Deep Eutectic Solvent Ethaline

Polarizable force fields are gradually becoming a common choice for ionic soft matter, in particular, for molecular dynamics (MD) simulations of ionic liquids (ILs) and deep eutectic solvents (DESs). The CL&Pol force field introduced in 2019 is the first general, transferable, and polarizable force field for MD simulations of different types of DESs. The original formulation contains, however, some problems that appear in simulations of ethaline and may also have a broader impact. First, the originally proposed atomic diameter parameters are unbalanced, resulting in too weak interactions between the chlorides and the hydroxyl groups of the ethylene glycol molecules. This, in turn, causes an artificial phase separation in long simulations. Second, there is an overpolarization of chlorides due to strong induced dipoles that give rise to the presence of peaks and antipeaks at very low q-vector values (2.4 nm^-1) in the partial components of the structure factors. In physical terms, this is equivalent to overestimated spatial nanoscale heterogeneity. To correct these problems, we adjusted the chloride-hydroxyl radial distribution functions against ab initio data and then extended the use of the Tang-Toennis damping function for the chlorides' induced dipoles. These adjustments correct the problems without losing the robustness of the CL&Pol force field. The results were also compared with the nonpolarizable version, the CL&P force field. We expect that the corrections will facilitate reliable use of the CL&Pol force field for other types of DESs.