A method combining the Metropolis algorithm with molecular dynamics to sample phase space more efficiently.