Computationally efficient canonical molecular dynamics simulations by using a multiple time-step integrator algorithm combined with the particle mesh Ewald method and with the fast multipole method.