Processing math: 100%

 

 

 

The Metropolis algorithm

Having chosen the acceptance probabilities, we have guaranteed that if the P(n)i has equilibrated, that is if it is equal to pi, it will remain equilibrated. Next we need to find the circumstances for convergence to equilibrium.

The dynamical equation can be written as P(n)i=jMijP(n1)j with the matrix M given by Mij=δij[1kTikAik]+TjiAji. Summing over i shows that iMij=1, and since kTik=1, and Aik1, the elements of the matrix satisfy Mij0. The matrix M is therefore a stochastic matrix.