In molecular dynamics simulations of systems with very hot particles in contact with a bath of cold particles one must choose a very small time step dictated by the fast (hot) particles. This requires the recalculation of all the forces after each of these small time steps. We show how reference system propagator algorithms (RESPA) can be used to speed up these simulations, and also comment on some interesting physical phenomena involved in these nonequilibrium systems.