Temperature Dependence of Short and Intermediate Range Order in Molten MgCl2 and Its Mixture with KCl

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.

Deep Learning Potential Assisted Prediction of Local Structure and Thermophysical Properties of the SrCl2-KCl-MgCl2 Melt

The local structure and thermophysical properties of SrCl2-KCl-MgCl2 melt were revealed by deep potential molecular dynamicsdriven by machine learning to facilitate the development of molten salt electrolytic Mg-Sr alloys. The short- and intermediate-range order of the SrCl2-KCl-MgCl2 melts was explored through radial distribution functions and structure factors, respectively, and their component and temperature dependence were discussed comprehensively. In the MgCl2-rich system, the intermediate-range order is more pronounced, and its evolution with temperature exhibits a non-Debye-Waller behavior. Mg-Cl is dominated by 4,5 coordination and Sr-Cl by 6,7 coordination, and their coordination geometries exhibit distorted octahedra and distorted pentagonal bipyramids, respectively. A database of thermophysical properties of SrCl2-KCl-MgCl2 melts, including density, self-diffusion coefficient, viscosity, and ionic conductivity, was thus developed, covering the temperature range from 873 to 1173 K.

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.

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.

Tracing mechanistic pathways and reaction kinetics toward equilibrium in reactive molten salts

In the dynamic environment of multi-component reactive molten salts, speciation unfolds as a complex process, involving multiple competing reaction pathways that are likely to face free energy barriers before reaching the reaction equilibria. Herein, we unravel intricate speciation in the AlCl3-KCl melt compositions with rate theory and ab initio molecular dynamics simulations. We find that the compositions with 100 and 50 mol% AlCl3 exclusively comprise neutral Al2Cl6 dimers and charged AlCl4- monomers, respectively. In intermediate AlCl3-KCl compositions, the chemical speciation proves to be a very complex process, requiring over 0.5 nanosecond to reach an equilibrium distribution of multiple species. It is a consequence of the competitive formation and dissociation of additional species, including charged Al dimers, trimers, and tetramers. Here, the species formation occurs through ion exchange events, which we explain by computing free energy landscapes and employing a Marcus-like rate theory. We show that both interspecies and intraspecies ion exchanges are probable and are dictated by the local structural reorganization reflected in the change of local coulombic fields. The species distributions are validated by comparing computed Raman spectra and neutron structure factors with the available experimental data. We find an excellent simulation-experiment agreement in both cases. Nevertheless, Raman spectroscopy turns out to be particularly advantageous for distinguishing between unique species distributions because of the distinct vibrational signatures of different species. The mechanistic insight into reaction dynamics gained in this study will be essential for the advancement of molten salts as reactive media in high-temperature energy applications.

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.

Are High-Temperature Molten Salts Reactive with Excess Electrons? Case of ZnCl2

New and exciting frontiers for the generation of safe and renewable energy have brought attention to molten inorganic salts of fluorides and chlorides. This is because high-temperature molten salts can act both as coolants and liquid fuel in next-generation nuclear reactors. Whereas research from a few decades ago suggests that salts are mostly unreactive to radiation, recent experiments hint at the fact that electrons generated in such extreme environments can react with the melt and form new species including nanoparticles. Our study probes the fate of an excess electron in molten ZnCl2 using first-principles molecular dynamics calculations. We find that on the time scale accessible to our study, an excess electron can be found in one of three states; the lowest-energy state can be characterized as a covalent Zn2Cl5•2- radical ion, the other two states are a solvated Zn•+ species (ZnCl3•2-) and a more delocalized species that still has some ZnCl3•2- character. Since for each of these, the singly occupied molecular orbital (SOMO) where the excess charge resides has a distinct and well-separated energy, the different species can in principle be characterized by their own electronic spectra. The study also sheds light onto what is commonly understood as the spectrum of a transient radical species which can be from the SOMO onto higher energy states or from the melt to pair with the excess electron leaving a hole in the liquid.

Probing the local structure of FLiBe melts and solidified salts by in situ high-temperature NMR

The 2LiF-BeF2 (FLiBe) salt melt is considered the primary choice for a coolant and fuel carrier for the generation IV molten salt reactor (MSR). However, the basics of ionic coordination and short-range ordered structures have been rarely reported due to the toxicity and volatility of beryllium fluorides, as well as the lack of suitable high-temperature in situ probe methods. In this work, the local structure of FLiBe melts was investigated in detail using the newly designed HT-NMR method. It was found that the local structure was comprised of a series of tetrahedral coordinated ionic clusters (e.g., BeF42-, Be2F73-, Be3F104-, and polymeric intermediate-range units). Li+ ions were coordinated by BeF42- ions and the polymeric Be-F network through the analysis of the NMR chemical shifts. Using solid-state NMR, the structure of solid FLiBe solidified mixed salts was confirmed to form a 3D network structure, significantly similar to those of silicates. The above results provide new insights into the local structure of FLiBe salts, which verifies the strong covalent interactions of Be-F coordination and the specific structural transformation to the polymeric ions above 25% BeF2 concentration.

Development of Deep Potentials of Molten MgCl2-NaCl and MgCl2-KCl Salts Driven by Machine Learning

Molten MgCl2-based chlorides have emerged as potential thermal storage and heat transfer materials due to high thermal stabilities and lower costs. In this work, deep potential molecular dynamics (DPMD) simulations by a method combination of the first principle, classical molecular dynamics, and machine learning are performed to systemically study the relationships of structures and thermophysical properties of molten MgCl2-NaCl (MN) and MgCl2-KCl (MK) eutectic salts at the temperature range of 800-1000 K. The densities, radial distribution functions, coordination numbers, potential mean forces, specific heat capacities, viscosities, and thermal conductivities of these two chlorides are successfully reproduced under extended temperatures by DPMD with a larger size (5.2 nm) and longer timescale (5 ns). It is concluded that the higher specific heat capacity of molten MK is originated from the strong potential mean force of Mg-Cl bonds, whereas the molten MN performs better in heat transfer due to the larger thermal conductivity and lower viscosity, attributed to the weak interaction between Mg and Cl ions. Innovatively, the plausibility and reliability of microscopic structures and macroscopic properties for molten MN and MK verify the extensibilities of these two deep potentials in temperatures, and these DPMD results also provide detailed technical parameters to the simulations of other formulated MN and MK salts.

First-principles molecular dynamics simulations of UCln-MgCl2 (n = 3, 4) molten salts

Molten chlorides are a preferred choice for fast-spectrum molten salt reactors. Molten MgCl₂ forms eutectic mixtures with NaCl and is considered as a promising dilutant to dissolve fuel salts such as UCl₃ and UCl₄. However, the structure and chemical properties of UCl_n (n = 3, 4) in molten MgCl₂ are not well understood. Here we use first-principles molecular dynamics to investigate the molten salt system UCl_n-MgCl₂ (n = 3, 4) at various concentrations of U³⁺ and U⁴⁺. It is found that the coordination environment of Cl^- around U³⁺, especially in the first coordination shell, varies only slightly with the uranium concentration and that both the 7-fold coordinate (UCl₇^4-) and 6-fold coordinate (UCl₆^3-) structures dominate at ∼40%, leading to an average coordination number of 6.6-6.7. A network or polymeric structure of U³⁺ cations sharing Cl^- ions is extensively formed when the mole fraction of UCl₃ is greater than 0.2. In contrast, the average coordination number of Cl^- around U⁴⁺ is about 6.4 for a mole fraction of UCl₄, x(UCl₄), of 0.1 but decreases to 6.0 for x(UCl₄) = 0.2 and then stays at about 6.0-6.2 with the uranium concentration. The 6-fold coordinate structure (UCl₆^2-) is the most populous in UCl₄-MgCl₂, at about 60%. U-Cl network formation becomes dominant (>50%) only when x(UCl₄) > 0.5. Unlike Na⁺, Mg²⁺ forms a network structure with Cl^- ions and when x(UCl₃) or x(UCl₄) < 0.5, over 90% of Mg²⁺ ions are part of a network structure, implying the complex influences from Mg²⁺ on the coordination of Cl around U. The present work reveals the impact of MgCl₂ as a solvent for UCl_n (n = 3, 4) on the U-Cl coordination and structure, and motivates further studies of their transport properties and the tertiary systems containing MgCl₂-UCl_n.

Adding MgCl2 to Molten NaCl-UCl n (n=3, 4): Insights from First Principles Molecular Dynamics

Molten chlorides are proposed for fast-spectrum molten salt reactors. Molten MgCl 2 with NaCl forms eutectic mixtures and is considered as a promising dilutant to dissolve fuel salts such as UCl 3 and UCl 4 . Previous study suggests the formation of U-Cl network at U:Na=1:1 binary salt. However, it is unclear how the structure of UCl n (n = 3, 4) in NaCl will change after adding MgCl 2 in the salt. Here we use first-principles molecular dynamics to investigate the molten ternary salts NaCl-MgCl 2 -UCl n (n=3, 4) at various concentrations of Mg 2+ in NaCl-UCl n with a fixed ratio of Na:U at 1:1. It is found that the addition of Mg 2+ in NaCl-UCl 3 leads to higher coordination number (from 6.5 to 6.7) of Cl around U while the U-Cl network structure slightly decreases with the Mg concentration. Adding MgCl 2 to NaCl-UCl 4 , however, breaks down the U-Cl network more completely. We attribute the different behavior of adding Mg 2+ into NaCl-UCl 3 and NaCl-UCl 4 to the difference between U(III) and U(IV) in attracting Cl - ions to form the first coordination shell. The present work reveals the impact of MgCl 2 as a dilutant solvent to the NaCl-UCl n fuel salts, which will be helpful in further studies and understanding of the thermophysical and transport properties of the ternary systems.

A high temperature cell for investigating interfacial structure on the molecular scale in molten salt/alloy systems

In this work, we describe the design and development of an in situ neutron reflectometry cell for high temperature investigations of structural changes occurring at the interface between inorganic salts, in their molten state up to 800 °C, and corrosion resistant alloys or other surfaces. In the cell, a molten salt is confined by an annular ring of single crystal sapphire constrained between the sample substrate and a sapphire plate using two gold O-rings, enclosing a liquid salt volume of 20 ml, along with a dynamic cell volume to accommodate expansion of the liquid with heating. As a test case for the cell, we report on an in situ neutron reflectometry measurement of the interface between a eutectic salt mixture of MgCl₂-KCl (32:68 molar ratio) and a single crystal sapphire substrate at 450 °C, resulting in the formation of a 60 Å layer having a scattering length density of 1.72 × 10^-6 Å^-2. While the origin of this layer is uncertain, it is likely to have resulted from the salt reacting with an existing impurity layer on the sapphire substrate.

Computing the Liquidus of Binary Monatomic Salt Mixtures with Direct Simulation and Alchemical Free Energy Methods

We describe and validate a free-energy-based method for computing the liquidus for binary solid-liquid phase diagrams in molecular simulations of monatomic salts. The method is demonstrated by calculating the liquidus for LiCl-KCl and MgCl₂-KCl salt mixtures with the polarizable ion model (PIM). The free-energy-based method is cross-validated with direct coexistence simulations. Both techniques show excellent agreement with one another. Though the predictions of the PIM disagree with experiments, we use our free-energy-based approach to decouple the contributions of liquid mixture nonidealities and pure component solid-liquid equilibrium to the phase diagram. In both mixtures, the PIM accurately reproduces the liquid phase nonidealities but fails to predict the liquidus because it does not accurately predict the pure component melting temperature of LiCl or MgCl₂.

A Holistic Approach for Elucidating Local Structure, Dynamics, and Speciation in Molten Salts with High Structural Disorder

To examine ion solvation, exchange, and speciation for minority components in molten salts (MS) typically found as corrosion products, we propose a multimodal approach combining extended X-ray absorption fine structure (EXAFS) spectroscopy, optical spectroscopy, ab initio molecular dynamics (AIMD) simulations, and rate theory of ion exchange. Going beyond conventional EXAFS analysis, our method can accurately quantify populations of different coordination states of ions with highly disordered coordination environments via linear combination fitting of the EXAFS spectra of these coordination states computed from AIMD to the experimental EXAFS spectrum. In a case study of dilute Ni(II) dissolved in the ZnCl₂+KCl melts, our method reveals heterogeneous distributions of coordination states of Ni(II) that are sensitive to variations in temperature and melt composition. These results are fully explained by the difference in the chloride exchange dynamics at varied temperatures and melt compositions. This insight will enable a better understanding and control of ion solubility and transport in MS.

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.

Unraveling Local Structure of Molten Salts via X-ray Scattering, Raman Spectroscopy, and Ab Initio Molecular Dynamics

In this work, we resolve a long-standing issue concerning the local structure of molten MgCl₂ by employing a multimodal approach, including X-ray scattering and Raman spectroscopy, along with the theoretical modeling of the experimental spectra based on ab initio molecular dynamics (AIMD) simulations utilizing several density functional theory (DFT) methods. We demonstrate the reliability of AIMD simulations in achieving excellent agreement between the experimental and simulated spectra for MgCl₂ and 50 mol % MgCl₂ + 50 mol % KCl, and ZnCl₂, thus allowing structural insights not directly available from experiment alone. A thorough computational analysis using five DFT methods provides a convergent view that octahedrally coordinated magnesium in pure MgCl₂ upon melting preferentially coordinates with five chloride anions to form distorted square pyramidal polyhedra that are connected via corners and to a lesser degree via edges. This is contrasted with the results for ZnCl₂, which does not change its tetrahedral coordination on melting. Although the five-coordinate MgCl₅^3- complex was not considered in the early literature, together with an increasing tendency to form a tetrahedrally coordinated complex with decreasing the MgCl₂ content in the mixture with alkali metal chloride systems, current work reconciles the results of most previous seemingly contradictory experimental studies.

X-ray scattering reveals ion clustering of dilute chromium species in molten chloride medium

Enhancing the solar energy storage and power delivery afforded by emerging molten salt-based technologies requires a fundamental understanding of the complex interplay between structure and dynamics of the ions in the high-temperature media. Here we report results from a comprehensive study integrating synchrotron X-ray scattering experiments, ab initio molecular dynamics simulations and rate theory concepts to investigate the behavior of dilute Cr3+ metal ions in a molten KCl-MgCl2 salt. Our analysis of experimental results assisted by a hybrid transition state-Marcus theory model reveals unexpected clustering of chromium species leading to the formation of persistent octahedral Cr-Cr dimers in the high-temperature low Cr3+ concentration melt. Furthermore, our integrated approach shows that dynamical processes in the molten salt system are primarily governed by the charge density of the constituent ions, with Cr3+ exhibiting the slowest short-time dynamics. These findings challenge several assumptions regarding specific ionic interactions and transport in molten salts, where aggregation of dilute species is not statistically expected, particularly at high temperature.

Machine-Learning-Driven Simulations on Microstructure and Thermophysical Properties of MgCl2-KCl Eutectic

Theoretical studies on the MgCl₂-KCl eutectic heavily rely on ab initio calculations based on density functional theory (DFT). However, neither large-scale nor long-time calculations are feasible in the framework of the ab initio method, which makes it challenging to accurately predict some properties. To address this issue, a scheme based on ab initio calculation, deep neural networks, and machine learning is introduced. By training on high-quality data sets generated by ab initio calculations, a deep potential (DP) is constructed to describe the interaction between atoms. This work shows that the DP enables higher efficiency and similar accuracy relative to DFT. By performing molecular dynamics simulations with DP, the microstructure and thermophysical properties of the MgCl₂-KCl eutectic (32:68 mol %) are investigated. The structural evolution with temperature is analyzed through partial radial distribution functions, coordination numbers, angular distribution functions, and structural factors. Meanwhile, the estimated thermophysical properties are discussed, including density, thermal expansion coefficient, shear viscosity, self-diffusion coefficient, and specific heat capacity. It reveals that the Mg²⁺ ions in this system have a distorted tetrahedral geometry rather than an octahedral one (with vacancies). The microstructure of the MgCl₂-KCl eutectic shows the feature of medium-range order, and this feature will be enhanced at a higher temperature. All predicted thermophysical properties are in good agreement with the experimental results. The hydrodynamic radius determined from the shear viscosity and self-diffusion coefficient shows that the Mg²⁺ ions have a strong local structure and diffuse as if with an intact coordination shell. Overall, this work provides a thorough understanding of the microstructure and enriches the data of the thermophysical properties of the MgCl₂-KCl eutectic.

SEM-Drude Model for the Accurate and Efficient Simulation of MgCl2-KCl Mixtures in the Condensed Phase

There is a long history of models that to different extents reproduce structural and dynamical properties of high-temperature molten salts. Whereas rigid ion models can work fairly well for some of the monovalent salts, polarizability is fundamentally important when small divalent or multivalent cations are combined with significantly polarizable anions such as Cl^- to form networked liquids that display a first sharp diffraction peak. There are excellent polarizable ion models (PIMs) for these systems, but there has been little success with the less expensive Core-Shell type models, which are often described as unwieldy or difficult to fit. In this article, we present the Sharma-Emerson-Margulis (SEM)-Drude model for MgCl₂/KCl mixtures that with the same ingredients used in the latest and most accurate PIM models overcome the aforementioned obstacles at significantly less computational cost; structural and dynamical properties are for all practical purposes very similar to what we obtain from the PIM but typical simulations can be more than 30 times faster. This has allowed us not only to expand our recent studies on the temperature and composition dependence of intermediate range order in MgCl₂/KCl mixtures but also to access transport properties that were simply too costly to properly sample in our recently published studies.

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

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