Molecular Dynamics Simulations of Molten Magnesium Chloride Using Machine‐Learning‐Based Deep Potential

Transferability and Accuracy of Ionic Liquid Simulations with Equivariant Machine Learning Interatomic Potentials

Ionic liquids (ILs) are an exciting class of electrolytes finding applications in many areas from energy storage to solvents, where they have been touted as "designer solvents" as they can be mixed to precisely tailor the physiochemical properties. As using machine learning interatomic potentials (MLIPs) to simulate ILs is still relatively unexplored, several questions need to be answered to see if MLIPs can be transformative for ILs. Since ILs are often not pure, but are either mixed together or contain additives, we first demonstrate that a MLIP can be trained to be compositionally transferable; i.e., the MLIP can be applied to mixtures of ions not directly trained on, while only being trained on a few mixtures of the same ions. We also investigated the accuracy of MLIPs for a novel IL, which we experimentally synthesize and characterize. Our MLIP trained on ∼200 DFT frames is in reasonable agreement with our experiments and DFT.

Machine-Learning-Based Interatomic Potentials for Group IIB to VIA Semiconductors: Toward a Universal Model

Rapid advancements in machine-learning methods have led to the emergence of machine-learning-based interatomic potentials as a new cutting-edge tool for simulating large systems with ab initio accuracy. Still, the community awaits universal interatomic models that can be applied to a wide range of materials without tuning neural network parameters. We develop a unified deep-learning interatomic potential (the DPA-Semi model) for 19 semiconductors ranging from group IIB to VIA, including Si, Ge, SiC, BAs, BN, AlN, AlP, AlAs, InP, InAs, InSb, GaN, GaP, GaAs, CdTe, InTe, CdSe, ZnS, and CdS. In addition, independent deep potential models for each semiconductor are prepared for detailed comparison. The training data are obtained by performing density functional theory calculations with numerical atomic orbitals basis sets to reduce the computational costs. We systematically compare various properties of the solid and liquid phases of semiconductors between different machine-learning models. We conclude that the DPA-Semi model achieves GGA exchange-correlation functional quality accuracy and can be regarded as a pretrained model toward a universal model to study group IIB to VIA semiconductors.

Development of NaCl-MgCl2-CaCl2 Ternary Salt for High-Temperature Thermal Energy Storage Using Machine Learning

NaCl-MgCl2-CaCl2 eutectic ternary chloride salts are potential heat transfer and storage materials for high-temperature thermal energy storage. In this study, first-principles molecular dynamics simulation results were used as a data set to develop an interatomic potential for ternary chloride salts using a neural network machine learning method. Deep potential molecular dynamics (DPMD) simulations were performed to predict the microstructure and thermophysical properties of the NaCl-MgCl2-CaCl2 ternary salt. This work reveals that DPMD simulations can accurately calculate the microstructure and thermophysical properties of ternary chloride salts. The association strength of chloride ions and cations follows the order of Mg2+ > Ca2+ > Na+, and the coordination number decreases gradually with increasing temperature, indicating a progressively looser and more disordered molten structure. Furthermore, thermophysical properties, such as density, specific heat capacity, thermal conductivity, and viscosity, are in good agreement with the experimental measurements. Machine learning molecular dynamics will provide a feasible multivariate molten salt exploration method for the design of next-generation solar power plants and thermal energy storage systems.

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.

Computational methods to simulate molten salt thermophysical properties

Molten salts are important thermal conductors used in molten salt reactors and solar applications. To use molten salts safely, accurate knowledge of their thermophysical properties is necessary. However, it is experimentally challenging to measure these properties and a comprehensive evaluation of the full chemical space is unfeasible. Computational methods provide an alternative route to access these properties. Here, we summarize the developments in methods over the last 70 years and cluster them into three relevant eras. We review the main advances and limitations of each era and conclude with an optimistic perspective for the next decade, which will likely be dominated by emerging machine learning techniques. This article is aimed to help researchers in peripheral scientific domains understand the current challenges of molten salt simulation and identify opportunities to contribute.

Exploring the Effects of Ionic Defects on the Stability of CsPbI3 with a Deep Learning Potential

Inorganic metal halide perovskites, such as CsPbI₃ , have recently drawn extensive attention due to their excellent optical properties and high photoelectric efficiencies. However, the structural instability originating from inherent ionic defects leads to a sharp drop in the photoelectric efficiency, which significantly limits their applications in solar cells. The instability induced by ionic defects remains unresolved due to its complicated reaction process. Herein, to explore the effects of ionic defects on stability, we develop a deep learning potential for a CsPbI₃ ternary system based upon density functional theory (DFT) calculated data for large-scale molecular dynamics (MD) simulations. By exploring 2.4 million configurations, of which 7,730 structures are used for the training set, the deep learning potential shows an accuracy approaching DFT-level. Furthermore, MD simulations with a 5,000-atom system and a one nanosecond timeframe are performed to explore the effects of bulk and surface defects on the stability of CsPbI₃ . This deep learning potential based MD simulation provides solid evidence together with the derived radial distribution functions, simulated diffraction of X-rays, instability temperature, molecular trajectory, and coordination number for revealing the instability mechanism of CsPbI₃ . Among bulk defects, Cs defects have the most significant influence on the stability of CsPbI₃ with a defect tolerance concentration of 0.32 %, followed by Pb and I defects. With regards to surface defects, Cs defects have the largest impact on the stability of CsPbI₃ when the defect concentration is less than 15 %, whereas Pb defects act play a dominant role for defect concentrations exceeding 20 %. Most importantly, this machine-learning-based MD simulation strategy provides a new avenue to explore the ionic defect effects on the stability of perovskite-like materials, laying a theoretical foundation for the design of stable perovskite materials.

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.

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.