A numerical technique that propagates the positions and momenta of interacting particles by integrating equations of motion (classical or ab initio) to model the time evolution of many-body systems and to sample dynamical and thermodynamic properties.