Abstract: Ionic compounds exhibit a variety of crystal structures that can critically affect their applications. Traditionally, relative sizes of cations and anions have been used to explain coordination of ions within the crystals. Such approaches assume atoms to be hard spheres and they cannot explain the observed structures of some crystals. Here we develop an atomistic method and use it to explore the structure-determining factors beyond the limitations of the hard sphere approach. Our approach is based upon a calibrated interatomic potential database that uses independent intrinsic bond lengths to measure atomic sizes. By carrying out extensive atomistic simulations, striking relationships among intrinsic bond lengths are discovered to determine the B1 (NaCl), B2 (CsCl), and B3 (zinc-blende) structure of binary ionic compounds.

Abstract: Unraveling the liquid structure of multicomponent molten salts is challenging due to the difficulty in conducting and interpreting high-temperature diffraction experiments. Motivated by this challenge, we developed composition-transferable Gaussian approximation potential (GAP) for molten LiCl-KCl. A DFT-SCAN accurate GAP is active-learned from only ~1100 training configurations drawn from 10 unique mixture compositions enriched with metadynamics. The GAP-computed structures show strong agreement across high-energy x-ray diffraction experiments, including for a eutectic not explicitly included in model training, thereby opening the possibility of composition discovery.

Abstract: The in silico modeling of molten salts is critical for emerging "carbon-free" energy applications but is inhibited by the cost of quantum mechanically treating the high polarizabilities of molten salts. Here, we integrate configurational sampling using classical force fields with active learning to automate and accelerate the generation of Gaussian approximation potentials (GAP) for molten salts. This methodology reduces the number of expensive ab initio evaluations required for training set generation to O(100), enabling the facile parametrization of a molten LiCl GAP model that exhibits a 19000-fold speedup relative to AIMD. The developed molten LiCl GAP model is applied to sample extended spatiotemporal scales, permitting new physical insights into molten LiCl's coordination structure as well as experimentally validated predictions of structures, densities, self-diffusion constants, and ionic conductivities. The developed methodology significantly lowers the barrier to the in silico understanding and design of molten salts across the periodic table.

Abstract: A computational approach combining dispersion-corrected density functional theory (DFT) and classical molecular dynamics is employed to characterize the geometrical and thermomechanical properties of a recently proposed 2D transition metal dihalide NiCl₂. The characterization is performed using a classical interatomic force field whose parameters are determined and verified through the comparison with the results of DFT calculations. The developed force field is used to study the mechanical response, thermal stability, melting and solidification of a NiCl₂ monolayer on the atomistic level of detail. The 2D NiCl₂ sheet is found to be thermally stable at temperatures below its melting point of ~695 K. At higher temperatures, several subsequent structural transformations of NiCl₂ are observed, namely a transition into a porous 2D sheet and a 1D nanowire. The MD simulations of NiCl₂ cooling show that the molten NiCl₂ system solidifies into an amorphous porous 2D structure at T ~450 K. The resulting structure has lower cohesive energy with respect to the initial 2D sheet. The computational methodology presented through the case study of NiCl₂ can also be utilized to study the properties of other novel 2D materials, including recently synthesized NiO₂, NiS₂, and NiSe₂.