Suitability of the ReaxFF potential for MD modeling of lithium across low and high temperatures

Modeling lithium's atomic-scale behavior is critical for its roles in plasma-facing fusion components and lithium-ion batteries yet remains challenging across phase regimes. This study benchmarks ReaxFF, 2NN-MEAM, and SNAP potentials from 100-1000 K using molecular dynamics, with ensemble and cooling protocols carefully controlled. Compared to reliable experimental data, ReaxFF diverges above 800 K, underestimating density by ∼10% at 1000 K and overestimating diffusivity-yielding anomalously low-density liquid behavior. Radial distribution functions and coordination profiles further reveal possible excessive disorder above 800 K. While ReaxFF offers qualitative insight into disordered lithium, its quantitative reliability diminishes at high temperatures, requiring validation against ab initio or experimental benchmarks. These findings inform potential selection for fusion- and battery-relevant simulations and underscore the sensitivity of glassy phase modeling to potential choice and thermal history.