Molecular dynamics (MD) simulations have been widely applied to investigate various physical and chemical processes in aqueous and interfacial environments, which are crucial for the design of energy materials and for understanding several chemical processes at the heart of life itself. However, the applicability of MD simulations has been constrained by several inherent challenges, including the accuracy of force fields, limitations in simulation size and time scales. One promising solution to
