Calculation of the Metropolis and the Glauber

Transition Probabilities for the Ising Model

and for the q-state Potts Modelby Peter MeyerAppendix I from the author's M.Phil. thesis

(i) Transition Probabilities in the Ising Model

Consider a spin model in a specific state, a particular spin S

_{i}and the set { S_{r}: S_{r}is a nearest neighbour of S_{i}}. The energy E_{i}contributed by the spin S_{i}is the sum of its interaction energies with the S_{r}and soE

_{i}= -J.Σ_{r}(S_{i}.S_{r}) = -J.S_{i}.Σ_{r}S_{r}In the Ising model there are only two possible spin values, +1 and -1, so the only possible change to the spin S

_{i}is for it to take the opposite value, namely, -S_{i}. Call this spin S_{j}. In this case the energy E_{j}associated with the new spin value is the sum of the interaction energies of S_{j}with the S_{r}and soE

_{j}= -J.Σ_{r}(S_{j}.S_{r}) = -J.Σ_{r}(-S_{i}.S_{r}) = -J.-S_{i}.Σ_{r}S_{r}= J.S_{i}.Σ_{r}S_{r}Thus the difference in the "before" and "after" energies is

ΔE = E _{j}- E_{i}= ( J.S_{i}.Σ_{r}S_{r}) - ( -J.S_{i}.Σ_{r}S_{r}) = 2.J.S_{i}.Σ_{r}S_{r}(1)Suppose that there is a maximum number of nearest neighbours which can be possessed by any spin, and let this number be denoted by

n. Then the maximum value for Σ_{r}S_{r}is +n and the minimum value is -n. Thus there are, at most, 2n+1 possible values for ΔE, corresponding to the possible values for Σ_{r}S_{r}, namely:-n, -n+1, -n+2, ..., -2, -1, 0, 1, 2, ..., n-1, n

although in a particular Ising spin model not all of these values may in fact be possible.

The transition probabilities for the Metropolis algorithm and for the Glauber algorithm were given in Section 1.7 as follows:

Metropolis algorithm: W(S _{i}→ S_{j})= 1 if ΔE _{ji}≤ 0= e ^{-ΔEji/(k}B^{T)}otherwise.Glauber algorithm: W(S _{i}→ S_{j})= 1 / [ 1 + e ^{ΔEji/(k}B^{T)}]where

Tis the temperature and k_{B}is Boltzmann's constant. From equation (1) we obtain:

Metropolis: W(S _{i}→ S_{j})= 1 if 2.J.S _{i}.Σ_{r}S_{r}≤ 0= e ^{-2JS}i^{.Σ}r^{S}r^{/k}B^{T}otherwise.Glauber: W(S _{i}→ S_{j})= 1 / [ 1 + e ^{2JS}i^{.Σ}r^{S}r^{/(k}B^{T)}]Taking J = k

_{B}= 1 the algorithms become:

Metropolis: W(S _{i}→ S_{j})= 1 if 2.S _{i}.Σ_{r}S_{r}≤ 0= e ^{-2.S}i^{.Σ}r^{S}r^{/T}otherwise.Glauber: W(S _{i}→ S_{j})= 1 / [ 1 + e ^{2.S}i^{.Σ}r^{S}r^{/T)}]Let S

_{2}= 2.S_{i}.Σ_{r}S_{r}then we have:

Metropolis: W(S _{i}→ S_{j})= 1 if S _{2}≤ 0= e ^{-S}2^{/T}otherwise.Glauber: W(S _{i}→ S_{j})= 1 / ( 1 + e ^{S}2^{/T})There are at most 2n+1 possible values for S

_{2}(assuming a maximum ofnnearest neighbours) as follows:-2n, -2n+2, -2n+4, ..., -2, 0, 2, ..., 2n-2, 2n

so the transition probabilities that we seek are:

Metropolis: S_{2}: -2n -2n+2 ... -2 0 2 4 ... 2n W: 1 1 ... 1 1 e^{-2/T}e^{-4/T }...^{ }e^{-2n/T}Glauber: S_{2}: -2n ... -2 0 2 ... 2n W: 1/(1+e^{-2n/T}) ... 1/(1+e^{-2/T}) 1/2 1/(1+e^{2/T}) ... 1/(1+e^{2n/T})

(ii) Transition Probabilities in the Potts Model, Version A.

Consider a spin model in a specific state, a particular spin S

_{i}(1 ≤ S_{i }≤ q) and the set { S_{r}: S_{r}is a nearest neighbour of S_{i}}. The energy E_{i}associated with the spin S_{i}is the sum of its interaction energies with the S_{r}and so E_{i}= -J.Σ_{r}[δ(S_{i},S_{r})].In the q-state Potts model there are

qpossible spin values, 1, ...,q, so the spin S_{i}may change to any of q-1 possible new values. We select one at random; call this spin S_{j}. In this case the energy E_{j}associated with the new spin value is the sum of the interaction energies of S_{j}with the S_{r}and so E_{j}= -J.Σ_{r}[δ(S_{j},S_{r})], where δ() is the Kronecker delta function: δ(x,y) = 1 if and only if x = y (otherwise 0). Thus the transition energy, ΔE, isE

_{j}- E_{i}= -J.Σ_{r}[δ(S_{j},S_{r})] - -J.Σ_{r}[δ(S_{i},S_{r})] = -J.{ Σ_{r}[δ(S_{j},S_{r})] - Σ_{r}[δ(S_{i},S_{r})] } soΔE = J.{ Σ _{r}[δ(S_{i},S_{r})] - Σ_{r}[δ(S_{j},S_{r})] } (2)Suppose that there is a maximum number of nearest neighbours which can be possessed by any spin, and let this number be denoted by

n. Then the maximum value for each of Σ_{r}[δ(S_{i},S_{r})] and Σ_{r}[δ(S_{j},S_{r})] is +n and the minimum value is 0. Thus the maximum value forS

_{r}[δ(S_{i},S_{r})] - Σ_{r}[δ(S_{j,}S_{r})] (3)is +n and the minimum is -n, so (as with the Ising model) there are, at most, 2n+1 possible values for ΔE (and thus for the transition probabilities) depending on the possible values for (3), namely:

-n, -n+1, -n+2, ..., -2, -1, 0, 1, 2, ..., n-1, n

although in a particular q-state Potts spin model not all of these values may in fact be possible.

The Metropolis and the Glauber algorithms for the q-state Potts model, version A, are obtained by substituting the RHS of (2) for ΔE into the definitions given in Section (i) above. As before we take J = k

_{B}= 1, so the algorithms become:Metropolis: W(S

_{i}→ S_{j}) = 1 if Σ_{r}[δ(S_{i},S_{r})] - Σ_{r}[δ(S_{j},S_{r})] ≤ 0,= e

^{-{S}r^{[δ(S}i^{,S}r^{)] - S}r^{[δ(S}j^{,S}r^{)]}/T}otherwise,that is, W(S

_{i}→ S_{j}) = 1 if Σ_{r}[δ(S_{i},S_{r})] ≤ Σ_{r}[δ(S_{j},S_{r})],= e

^{-{S}r^{[δ(S}i^{,S}r^{)] - S}r^{[δ(S}j^{,S}r^{)]}/T}otherwise.Glauber: W(S

_{i}→ S_{j}) = 1 / ( 1 + e^{{S}r^{[δ(S}i^{,S}r^{)] - S}r^{[δ(S}j^{,S}r^{)]}/T})Let S

_{3}= Σ_{r}[δ(S_{i},S_{r})] - Σ_{r}[δ(S_{j},S_{r})] then we have:Metropolis: W(S

_{i}→ S_{j}) = 1 if S_{3}≤ 0,= e

^{-S}3^{/T}otherwise.Glauber: W(S

_{i}→ S_{j}) = 1 / ( 1 + e^{S}3^{/T})There are at most 2n+1 possible values for S

_{3}(assuming a maximum ofnnearest neighbours) as follows:-n, -n+1, -n+2, ..., -2, -1, 0, 1, 2, ..., n-1, n

so the transition probabilities that we seek are:

Metropolis: S_{3}: -n -n+1 ... -1 0 1 2 ... n W: 1 1 ... 1 1 e^{-1/T }e^{-2/T}... e^{-n/T}Glauber: S_{3}: -n ... -1 0 1 ... n W: 1/(1+e^{-n/T}) ... 1/(1+e^{-1/T}) 1/2 1/(1+e^{1/T}) ... 1/(1+e^{n/T})It will be noted that these are not the same transition probabilities as in the Ising case.

(iii) Transition Probabilities in the Potts Model, Version B.

We follow the same reasoning as in the previous section.

Consider a spin model in a specific state, a particular spin S

_{i}(1 ≤ S_{i }≤ q) and the set { S_{r}: S_{r}is a nearest neighbour of S_{i}}. The energy E_{i}associated with the spin S_{i}is the sum of its interaction energies with the S_{r}and so E_{i}= -J.Σ_{r}[2.δ(S_{i},S_{r})-1].In the q-state Potts model there are

qpossible spin values, 1, ...,q, so the spin S_{i}may change to any of q-1 possible new values. We select one at random; call this spin S_{j}. In this case the energy E_{j}associated with the new spin value is the sum of the interaction energies of S_{j}with the S_{r}, so E_{j}= -J.Σ_{r}[2.δ(S_{j},S_{r})-1]. Thus the transition energy, ΔE, isE

_{j}- E_{i}= -J.Σ_{r}[2.δ(S_{j},S_{r})-1] - -J.Σ_{r}[2.δ(S_{i},S_{r})-1] = -J.{ Σ_{r}[2.δ(S_{j},S_{r})] - Σ_{r}[2.δ(S_{i},S_{r})] } soΔE = 2.J.{ Σ _{r}[δ(S_{i},S_{r})] - Σ_{r}[δ(S_{j},S_{r})] } (4)Suppose that there is a maximum number of nearest neighbours which can be possessed by any spin, and let this number be denoted by

n. Then the maximum value for each of Σ_{r}[δ(S_{i},S_{r})] and Σ_{r}[δ(S_{j},S_{r})] is +n and the minimum value is 0. Thus the maximum value forS

_{r}[δ(S_{i},S_{r})] - Σ_{r}[δ(S_{j,}S_{r})] (5)is +n and the minimum is -n, so (as with the Ising model and the Potts model, Version A) there are, at most, 2n+1 possible values for ΔE (and thus for the transition probabilities) depending on the possible values for (5), namely:

-n, -n+1, -n+2, ..., -2, -1, 0, 1, 2, ..., n-1, n

although in a particular q-state Potts spin model not all of these values may in fact be possible.

The Metropolis and the Glauber algorithms for the q-state Potts model, Version B, are obtained by substituting the RHS of (4) for ΔE into the definitions given in Section (i) above. As before we take J/k

_{B}= 1, so the algorithms become:Metropolis: W(S

_{i}→ S_{j}) = 1 if 2.{ Σ_{r}[δ(S_{i},S_{r})] - Σ_{r}[δ(S_{j},S_{r})] } ≤ 0,= e

^{-2.{S}r^{[δ(S}i^{,S}r^{)] - S}r^{[δ(S}j^{,S}r^{)]}/T}otherwise,that is, W(S

_{i}→ S_{j}) = 1 if Σ_{r}[δ(S_{i},S_{r})] ≤ Σ_{r}[δ(S_{j},S_{r})],= e

^{-2.{S}r^{[δ(S}i^{,S}r^{)] - S}r^{[δ(S}j^{,S}r^{)]}/T}otherwise.Glauber: W(S

_{i}→ S_{j}) = 1 / ( 1 + e^{2.{S}r^{[δ(S}i^{,S}r^{)] - S}r^{[δ(S}j^{,S}r^{)]}/T})Let S

_{4}= 2.{Σ_{r}[δ(S_{i},S_{r})] - Σ_{r}[δ(S_{j},S_{r})]} then we have:Metropolis: W(S

_{i}→ S_{j}) = 1 if S_{4}≤ 0,= e

^{-S}4^{/T}otherwise.Glauber: W(S

_{i}→ S_{j}) = 1 / ( 1 + e^{S}4^{/T})There are at most 2n+1 possible values for S

_{4}(assuming a maximum ofnnearest neighbours) as follows:-2n, -2n+2, -2n+4, ..., -2, 0, 2, ..., 2n-2, 2n

so the transition probabilities that we seek are:

Metropolis: S_{4}: -2n -2n+2 ... -2 0 2 4 ... 2n W: 1 1 ... 1 1 e^{-2/T }e^{-4/T }...^{ }e^{-2n/T}Glauber: S_{4}: -2n ... -2 0 2 ... 2n W: 1/(1+e^{-2n/T}) ... 1/(1+e^{-2/T}) 1/2 1/(1+e^{2/T}) ... 1/(1+e^{2n/T})It will be noted that these are not the same transition probabilities as for Version A of the Potts Model, but they are the same as those for the Ising model.

(iv) The Equivalence of the 2-state Potts Model and the Ising Model.

It is commonly said that the 2-state Potts model is equivalent to the Ising model. This is true if Version B of the Potts model is meant, because then the transition probabilities are the same and so the dynamics of the two models are the same. Usually, however, when the Potts model is discussed it is Version A which is meant, and in this case, even for q=2, the transition probabilities are not the same. Thus at the microscopic level there are differences between the Ising model and the 2-state Potts model, Version A, so we are led to ask whether the 2-state Potts model, Version A, is

strictlyequivalent to the Ising model.The transition probabilities (in both the Metropolis algorithm and the Glauber algorithm) are determined by (i) the transition energies, i.e., the energy difference between the initial state and a possible new state, (ii) the value of J/k

_{B}and (iii) the temperature. Since (ii) and (iii) are not properties of the spin model itself the question becomes whether (i) is the same regardless of whether the spins are described in terms of one model definition (Ising) or in terms of the other (2-state Potts, Version A). We shall first consider the 2-state Potts model, Version B.Version B:

Consider a spin S

_{i}and the setS= { S_{r}: S_{r}is a nearest neighbour of S_{i}}.Smay be partitioned intoS_{1}= { S_{r}: S_{r}= S_{i}} andS_{2}= { S_{r}: S_{r}≠ S_{i}}.Suppose S

_{i}and the S_{r}are viewed as Ising spins. Suppose S_{i}= +1, then (from section (i) above and the definition of the Ising model in Section 1.5):

ΔE

_{Ising}= 2.J_{Ising}.(+1).Σ_{r}S_{r}= 2.J_{Ising}.( |S_{1}| - |S_{2}| )Suppose S

_{i}= -1, thenΔE

_{Ising}= 2.J_{Ising}.(-1).Σ_{r}S_{r}= -2.J_{Ising}.( -|S_{1}| + |S_{2}| ) = 2.J_{Ising}.( |S_{1}| - |S_{2}| )Thus in each case ΔE

_{Ising }= 2.J_{Ising}.( |S_{1}| - |S_{2}| ).Now suppose S

_{i}and the S_{r}are viewed as spins in the 2-state Potts model, Version B. Then (from the definition of the Potts model, Version B, in Section 1.5)ΔE

_{Potts,B}= -J_{Potts,B}.2.{ Σ_{r}[δ(S_{j},S_{r})] - Σ_{r}[δ(S_{i},S_{r})] }where S

_{j}is some spin value other than S_{i}, soΔE

_{Potts,B}= -J_{Potts,B}.2.( |S_{2}| - |S_{1}| ) = J_{Potts,B}.2.( |S_{1}| - |S_{2}| )Since ΔE

_{Ising }= ΔE_{Potts,B}we may conclude that the Ising model is strictly equivalent to the 2-state Potts model, Version B.Version A:

As above, consider a spin S

_{i}and the setS= { S_{r}: S_{r}is a nearest neighbour of S_{i}}.Smay be partitioned intoS_{1}= { S_{r}: S_{r}= S_{i}} andS_{2}= { S_{r}: S_{r}≠ S_{i}}. As before, ΔE_{Ising }= 2.J_{Ising}.( |S_{1}| - |S_{2}| ).Now suppose S

_{i}and the S_{r}are viewed as spins in the 2-state Potts model, Version A. Then (from the definition of the Potts model, Version A, in Section 1.5)ΔE

_{Potts,A}= -J_{Potts,A}.{ Σ_{r}[δ(S_{j},S_{r})] - Σ_{r}[δ(S_{i},S_{r})] }where S

_{j}is some spin value other than S_{i}, soΔE

_{Potts,A}= -J_{Potts,A}.( |S_{2}| - |S_{1}| ) = J_{Potts,A}.( |S_{1}| - |S_{2}| )Since ΔE

_{Ising }= 2.ΔE_{Potts,A}the transition energies (and thus the transition probabilities) are different.Okano

et al.(1997, p.738) state: "It is known that the critical points [of the q-state Potts model] locate at J_{c}= log(1+√q)." Taking T_{c}= 1/J_{c}and q = 2 we obtainT

_{c}= 1/ln(1+√2) = 1.134593This is exactly one-half of the value usually given for the Ising model, namely, 2/ln(1+√2) = 2.269185 (see e.g. Stinchcombe 1983, p.177, taking J = 1 = k

_{B}), which value is obtained from simulations. Thus although the Ising model and the 2-state Potts model are commonly said to be equivalent, the values of the critical temperatures are not the same when Version A of the 2-state Potts model is used.Nevertheless, the critical exponents of the Ising model and the 2-state Potts model, Version A, are found to be the same, the only difference being the value of the critical temperature. Thus the difference between the two models is really only one of the scale used for the temperature.

Consider the Potts model, Version A, in which the interaction energy is J

_{Potts,A}= 2.J_{Ising}thenΔE

_{Potts,A}= 2.J_{Ising}.( |S_{2}| - |S_{1}| ) = ΔE_{Ising}so we can also say that the 2-state Potts model, Version A,

isstrictly equivalent to the Ising model provided the interaction energy of the 2-state Potts model is double that of the Ising model. Since the units of the critical temperature are related to units ofJ, and a J-unit in this modified Potts model, Version A, denotes twice the energy of a J-unit in the Ising model, the numerical value of the critical temperature in the former model is thus one-half of its value in the latter model, thus explaining the difference in the numerical values of the critical temperature in the Ising model and in the 2-state Potts model, Version A.

Computational Studies of Pure and Dilute Spin Models Computational Science Hermetic SystemsHome Page