Multiscale and Multi‐Timestep Switching of Multiple Machine Learning Force Fields for Artificial Intelligence‐Driven Materials Simulations

Molecular dynamics (MD) is essential for investigating atomic‐scale processes in materials and molecular systems, but the cost of high‐accuracy machine learning force field simulations still limits accessible system sizes and timescales. Here, we propose a practical model‐switching strategy for Deep Potential (DP)‐based MD simulations that alternates between independently trained DP models with different cutoff radii: a standard 6 Å model for higher accuracy and a reduced‐cutoff 4 Å model for faster inference. The method was implemented in LAMMPS/DeePMD and evaluated using solid‐phase anatase TiO ₂ and liquid‐phase polyethylene glycol (PEG). For anatase TiO ₂ , the 1:3 4–6 Å switching scheme preserved radial distribution function (RDF) correlations of 0.996 or higher relative to the 6 Å baseline while achieving a 1.24‐fold speedup. For PEG, the switching scheme maintained RDF correlations of 0.996 or higher with a 1.18‐fold speedup. Additional optimization using network‐size reduction and mixed‐precision inference achieved a 2.53‐fold speedup with RDF correlations of 0.975–0.988. Constant particle‐number, pressure, and temperature (NPT) simulations remained stable, whereas constant particle‐number, volume, and energy (NVE) simulations revealed system‐dependent energy‐drift behavior, particularly for aggressively optimized models. These results demonstrate that DP model switching provides a simple and practical route for accelerating structural MD simulations while highlighting the need for validation when strict energy conservation is required.