The Kramers equation simulation algorithm I. Operator analysis

Using an operatorial formalism, we study the Kramers equation and its applications to numerical simulations. We obtain classes of algorithms which may be made precise at every desired order in the time step $\epsilon$ and with a set of free parameters which can be used to reduce autocorrelations. We show that it is possible to use a global Metropolis test to restore Detailed Balance.