# Monte Carlo Simulation of the Ising Model

## Introduction

The simplest, nontrivial model that has a phase transition is the Ising model, one of most important models in statistical mechanics. The model consists of spins located on a lattice such that each spin can take on one of two values designated as up and down or ±1. The interaction energy between two neighboring spins is -J if the two spins are in the same state and +J if they are in opposite states. The Ising model undergoes a phase transition between an ordered and a disordered phase in two dimensions or more. We will consider the Ising model in the context of magnetism, but many other systems can be understood in terms of the Ising model.

The model was proposed by Wilhelm Lenz (1888–1957) in 1920 and was solved exactly for the one-dimensional case by his student Ernst Ising in 1925. The one-dimensional case does not have a phase transition. Lars Onsager (1903–1976) solved the Ising model exactly in 1944 for two dimensions in the absence of an external magnetic field and showed that there was a phase transition in two dimensions. No exact solution is available in two dimensions in an arbitrary magnetic field and in three dimensions. Hence, simulations of the Ising model are essential.

The total energy can be expressed in the form

,

where si =±1 and J is the exchange constant. We will refer to s as the spin.

### References

Each year many hundreds of papers are published that apply the Ising model to problems in such diverse fields as neural networks, protein folding, biological membranes, and social behavior.

1. A biographical note about Ising's life is at <www.bradley.edu/las/phy/personnel/ising.html>. Also see S. Kobe, "Ernst Ising 1900–1998," Braz. J. Phys. 30 (4), 649–654 (2000). Available from various Web sites.
2. See <en.wikipedia.org/wiki/Lars_Onsager> for a summary of Onsager's life and contributions.
3. Stephen G. Brush, "History of the Lenz-Ising model," Rev. Mod. Phys. 39, 883–893 (1967).
4. Many textbooks on statistical mechanics discuss the Ising model. See for example, M. E. J. Newman and G. T. Barkema, Monte Carlo Methods in Statistical Physics, Oxford University Press (1999).
5. H. Gould, J. Tobochnik, and W. Christian, An Introduction to Computer Simulation Methods: Applications to Physical Systems, Addison-Wesley (2006), Chapter 15.

### Note:

This simulation was created by Wolfgang Christian using the Easy Java Simulations (Ejs) modeling tool. It is based on a Java program from An Introduction to Computer Simulation Methods. You can modify this simulation if you have Ejs installed by right-clicking within a plot and selecting "Open Ejs Model" from the pop-up menu. Information about Ejs is available at: <http://www.um.es/fem/Ejs/>.