Glauber dynamics

The Ising model is an abstract model for the magnetic interaction of neighboring atoms. It is conventionally considered on a two-dimensional square lattice, with magnetic interactions occurring only between nearest neighbors. In this model, each lattice site is given a spin $\backslash sigma\_\{x,y\}$ that is either up (+1) or down (-1); the x and y are the grid coordinates. Glauber's algorithm becomes:{{cite journal|title= An introduction to Monte Carlo methods |year= 2015 |doi= 10.1016/j.physa.2014.06.014 |arxiv= 1404.0209 |language=en-US|last1= Walter |first1= J.-C. |last2= Barkema |first2= G.T. |journal= Physica A: Statistical Mechanics and Its Applications |volume= 418 |pages= 78–87 |s2cid= 118589022 }}

1. Choose a location $x,y$ at random.
2. Sum the spins of the nearest-neighbors. For a two-D square lattice, there are four: $S\; =\; \backslash sigma\_\{x+1,y\}\; +\; \backslash sigma\_\{x-1,y\}\; +\; \backslash sigma\_\{x,y+1\}\; +\; \backslash sigma\_\{x,y-1\}$.
3. Compute the change in energy if the spin at x, y were to flip. This is given by the Hamiltonian for the Ising model; it is $\backslash Delta\; E\; =\; 2\backslash sigma\_\{x,\; y\}\; S.$
4. Flip the spin with probability given by the Fermi function $p(\backslash Delta\; E)\; =\; 1/(1\; +\; e^\{\backslash Delta\; E/T\}),$ where T is the temperature.
5. Go To Step 1.

In Glauber algorithm, if the energy change in flipping a spin is zero, $\backslash Delta\; E\; =\; 0$, then the spin would flip with probability $p(\backslash Delta\; E)\; =\; 0.5$. Similarly, if the temperature is very high, then $e^\{\backslash Delta\; E/T\}\backslash approx\; 1$ and again, the probability of a flip is one-half. If the temperature is very low, then a flip to a higher-energy state almost never happens, while a flip to a lower energy state almost always happens.

File:Glauber Acceptance Probability.png

The Glauber algorithm can be compared to the Metropolis–Hastings algorithm. These two differ in how a spin site is selected (step 1), and in the probability of a spin-flip (step 4). In the Glauber dynamic, every spin has an equal chance of being chosen at each time step, and the decision to flip that spin, or not, is given by the Fermi function, as given above.

By contrast, the Metropolis algorithm considers a spin site with a probability given by the Boltzmann weight $e^\{-\backslash Delta\; E/T\}$, but if it is accepted, then it always flips a spin in favor of lowering the energy. Thus, the total spin-flip probability is:

:$p(\backslash Delta\; E)\; =\; \backslash left\backslash \{\backslash begin\{matrix\}$

1, \ \ \ \ \ \ \ \ \ \Delta E \leqslant 0\\

e^{-\Delta E/T}, \ \Delta E > 0\\

\end{matrix}\right.

File:Metropolis-Hastings Acceptance Probability.pngAlthough both of the acceptance probabilities approximate a step curve and they are almost indistinguishable at very low temperatures, they differ when temperature gets high. For an Ising model on a 2d lattice, the critical temperature is $T\; =\; 2.27$.

At thermal equilibrium, the Glauber and Metropolis algorithms should give identical results. In general, at equilibrium, any MCMC algorithm should produce the same distribution, as long as the algorithm satisfies ergodicity and detailed balance. In both algorithms, for any change in energy, $p(\backslash Delta\; E)\; \backslash neq\; 0$, meaning that transition between the states of the system is always possible despite being very unlikely at some temperatures. So, the condition for ergodicity is satisfied for both of the algorithms. Detailed balance, which is a requirement of reversibility, states that if you observe the system for a long enough time, the system goes from state $A$ to $B$ with the same frequency as going from $B$ to $A$. In equilibrium, the probability of observing the system at state A is given by the Boltzmann weight, $e^\{-E\_A/T\}$. So, the amount of time the system spends in low energy states is larger than in high energy states and there is more chance that the system is observed in states where it spends more time. Meaning that when the transition from $A$ to $B$ is energetically unfavorable, the system happens to be at $A$ more frequently, counterbalancing the lower intrinsic probability of transition. Therefore, both, Glauber and Metropolis–Hastings algorithms exhibit detailed balance.

• Simulation package IsingLenzMC provides simulation of Glauber Dynamics on 1D lattices with external field. [https://cran.r-project.org/web/packages/isingLenzMC/index.html CRAN].

• Metropolis algorithm
• Ising model
• Monte Carlo algorithm
• Simulated annealing

Category:Monte Carlo methods

Category:Spin models