Application of MR2T2 algorithm in statistical mechanics Canonical ensemble: Distribution function: p(rN)=exp(-E(rN)/KT)/Q Q= JdrN exp(-E(rN)/kT) Different states of the Markov chain: In the application of Metropolis algorithm to a simulation of molecular system, the states of the Markov chain correspond to different configurations of the molecular system