Random energy model

In statistical physics of disordered systems, the random energy model is a toy model of a system with quenched disorder. It concerns the statistics of a system of $N$ particles, such that the number of possible states for the systems grow as $2^{N}$, while the energy of such states is a Gaussian stochastic variable. The model has an exact solution. Its simplicity makes this model suitable for pedagogical introduction of concepts like quenched disorder and replica symmetry.

Comparison with other disordered systems

The $r$-spin Infinite Range Model, in which all $r$-spin sets interact with a random, independent, identically distributed interaction constant, becomes the Random-Energy Model in a suitably defined $r\to \infty$ limit.

More precisely, if the Hamiltonian of the model is defined by

$H(\sigma )=\sum _{\{i_{1},\ldots ,i_{r}\}}J_{i_{1},\ldots i_{r}}\sigma _{i_{1}}\cdots \sigma _{i_{r}},$

where the sum runs over all ${N \choose r}$ distinct sets of $r$ indices, and, for each such set, $\{i_{1},\ldots ,i_{r}\}$, $J_{i_{1},\ldots ,i_{r}}$ is an independent Gaussian variable of mean 0 and variance $r!/(2N^{r-1})$, the Random-Energy model is recovered in the $r\to \infty$ limit.

Derivation of thermodynamical quantities

As its name suggests, in the REM each microscopic state has an independent distribution of energy. For a particular realization of the disorder, $P(E)=\delta (E-H[{\vec {\sigma }}])$ where ${\vec {\sigma }}$ refers to the set of individual spin configurations described by the state and $H[{\vec {\sigma }}]$ is the energy associated with it. Since there is disorder in the system, the final extensive variables like the free energy need to be averaged over all realizations of the disorder, just as in the case of the Edwards Anderson model. Averaging $P(E)$ over all possible realizations, we find that the probability of the disordered system having an energy $E$ is given by

$P(E)={\sqrt {\dfrac {1}{N\pi J^{2}}}}\exp \left(-{\dfrac {E^{2}}{J^{2}N}}\right).$

Moreover, the joint probability distribution of the energy values of two different microscopic configurations of the spins, $\sigma$ and $\sigma '$ factorizes:

$P(E,E')=P(E)\,P(E').$

It can be seen that the probability of a given spin configuration only depends on the energy of that state and not on the individual spin configuration.

The entropy of the REM is given by $S(E)=N\left[\log 2-\left({\dfrac {E}{NJ}}\right)^{2}\right]$ for $|E|<NJ{\sqrt {\log 2}}$.

Suppose a system is described by a total energy given by a sum of random energy

$E=\sum _{i}^{N}E_{i}$

suppose that these are independent and identical random variables with average ${\bar {E}}$ and standard deviation $\sigma$, then by the central limit theorem the energy E will be a random variable with gaussian distribution with mean $N{\bar {E}}$ and standard deviation ${\sqrt {N}}\sigma$.