Observing nucleation and crystallization of rock salt LiF from molten state through molecular dynamics simulations with refined machine-learned force field

Lithium fluoride is a critical component for stabilizing lithium metal anodes and high-voltage cathodes toward the next-generation high-energy-density lithium batteries. A recent modeling study reported the formation of wurtzite LiF below ∼550 K (Hu et al., J. Am. Chem. Soc. 2023, 145, 1327-1333), in contrast to the experimental observation of rock salt LiF under ambient conditions. To address this discrepancy, we employ molecular dynamics simulations with a refined machine-learned force field (MLFF) and demonstrate the nucleation and crystallization of rock salt LiF from the molten phase at temperatures below ∼800 K. The rock salt phase remains stable in LiF nanoparticles. Complementary density functional theory calculations show that dispersion interactions are essential for correctly predicting the thermodynamic stability of rock salt LiF over the wurtzite phase on top of the commonly used PBE functional. Furthermore, we show that inclusion of virial stresses, alongside energies and forces, in the training of MLFFs is crucial for capturing phase nucleation and crystallization of rock salt LiF under the isothermal-isobaric ensemble. These findings underscore the critical role of dispersion interactions in atomistic simulations of battery materials, where such effects are often non-negligible, and highlight the necessity of incorporating virial stresses during the training of MLFFs to enable accurate modeling of solid-state systems.