Bridging Atomistic and Mesoscale Lithium Transport via Machine-Learned Force Fields and Markov State Models

Abstract

Lithium diffusion in solid-state battery anodes occurs through thermally activated hops between metastable sites often separated by large energy barriers, making such events rare on ab initio molecular dynamics (AIMD) timescales. Here, we present a bottom-up multiscale workflow that integrates AIMD, machine-learned force fields (MLFFs), and Markov state models (MSMs) to establish a quantitatively consistent link between atomistic hopping mechanisms and mesoscale transport. MLFFs fine-tuned on AIMD reference data retain near-DFT accuracy while enabling large-scale molecular dynamics simulations extending to tens of nanoseconds. These extended trajectories remove the strong finite-size bias present in AIMD and yield diffusion coefficients in excellent agreement with experiment. Furthermore, from these long MLFF trajectories, we obtain statistically converged lithium jump networks and construct MSMs that remain Markovian across more than two orders of magnitude in the lag times used for their construction. The resulting MSMs faithfully reproduce mean-square displacements and recover rare diffusion processes that do not occur on AIMD timescales. In addition to propagating lithium distributions, the MSM transition matrices provide mechanistic insight: their eigenvalues and eigenvectors encode characteristic relaxation timescales and dominant transport pathways. Although demonstrated for defect-free crystalline Li$_x$Si$_y$ phases, the AIMD$\rightarrow$MLFF$\rightarrow$MSM framework is general and provides a transferable approach for describing lithium transport in amorphous materials, defect-mediated diffusion, and next-generation solid-state anodes.