Replica trick

From Infogalactic: the planetary knowledge core
Jump to: navigation, search

Lua error in package.lua at line 80: module 'strict' not found. In statistical physics of spin glasses and other systems with quenched disorder, the replica trick is a mathematical technique based on the application of the formula

\lim_{n\to 0} {Z^n-1\over n}=\ln Z.

Mathematical Trick

This mathematical trick is used in computation involving functions of a variable that can be expressed as a power series in that variable. The crux of this technique is to reduce the function of a variable, say f(z), into powers of z or, in other words, replicas of z, and perform the same computation which is to be done on f(z), using the powers of z.

A particular case which is of great use in physics is in averaging the free energy, or -\ln Z[J_{ij}], over values of J_{ij} with a certain probability distribution, typically Gaussian,[books on spin glasses 1] and the function Z[J_{ij}] \sim e^{-\beta J_{ij}}. Notice that if it were Z[J_{ij}](or more generally, any power of J_{ij}) and not its logarithm which we wanted to average, the resulting integral (assuming a Gaussian distribution) would be of the form \int dJ_{ij} e^{-\beta J - \alpha J^{2}}, which can be performed by completing squares and carrying out the standard gaussian integration. But we have the special property or the limit form expression for the logarithm function, given by:

\ln Z_{ij} = \lim_{n\to 0}\dfrac{Z^{n}-1}{n}

which clearly reduces the task of averaging to solving a relatively simpler Gaussian integral.[1] The replica trick involves extending this argument to the case where n is no longer constrained to be an integer, by positing that if Z^n can be calculated for all positive integers n then this may be sufficient to allow the limiting behaviour as n\to0 to be calculated.

Clearly, such an argument poses many mathematical questions, and the resulting formalism for performing the limit n\to0 typically introduces many subtleties (see Mezard et al.). When using mean field theory to perform one's calculations, taking this limit often requires introducing extra order parameters, in consequence of 'replica symmetry breaking' which is closely related to ergodicity breaking and slow dynamics within disorder systems.

Physical applications

The replica trick is used in determining ground states of statistical mechanical systems, in the mean field approximation. Typically, for systems in which the determination of ground state is easy, one can analyze fluctuations near the ground state. But in cases where, for some reason the determination of ground state is hard, one uses the replica method.[papers on spin glasses 1] An example is the case of a quenched disorder in a spin system like spin glass with different types of magnetic links between spin sites, thereby causing many configurations to have the same energy. Hence finding a particular ground state is hard.[papers on spin glasses 1]

In statistical physics of quenched disorder systems, any two states (set of configurations) with the same realization of the disorder, or in case of Spin glasses, with the same distribution of ferromagnetic and antiferromagnetic bonds, are called replicas of each other.[papers on spin glasses 2] For systems with quenched disorder, one typically expects that macroscopic quantities will be self-averaging, whereby any macroscopic quantity for a specific realization of the disorder will be indistinguishable from the same quantity calculated by averaging over all possible realizations of the disorder. Hence replicas are introduced for integrating out the disorder[books on spin glasses 1] in a system.

In the case of a spin glass, we expect the free energy per spin (or any self averaging quantity) in the thermodynamic limit to be independent of the particular values of ferromagnetic and antiferromagnetic couplings between individual sites, across the lattice. So, we explicitly find the free energy as a function of the disorder parameter (in this case, parameters of the distribution of ferromagnetic and antiferromagnetic bonds) and average the free energy over all realizations of the disorder (all values of the coupling between sites, each with its corresponding probability, given by the distribution function). As free energy takes the form: 
F = \overline{F[J_{ij}]} = -k_{B}T\overline{\ln Z[J]}
where J_{ij} describes the disorder (for spin glasses, it describes the nature of magnetic interaction between each of the individual sites i and j) and [\cdots ] denotes the average over all values of the couplings described in J, weighted with a given distribution. To perform the averaging over the logarithm function, the replica trick come in handy, in replacing the logarithm with its limit form mentioned above. In this case, the quantity Z^n represents the joint partition function of n identical systems.

REM: The easiest Replica problem

The Random Energy Model (REM) is one of the simplest models of statistical mechanics of disordered systems, and probably the simplest model to show the meaning and power of the Replica Trick to the level 1 of Replica Symmetry Breaking. The model is especially suitable for this introduction because an exact result by a different procedure is known, and the Replica Trick can be proved to work by crosschecking of results.

Proof of initial formula

We prove the formula

\lim_{n\to 0} {x^n-1\over n}=\ln x.

Start from the observation

\exp y = 1 + y + {1 \over 2!} y^2 + {1 \over 3!} y^3 + \dots = \lim_{N\to \infty} \sum_{r=0}^N {N! \over r! (N-r)!} ({y \over N})^r = \lim_{N\to \infty} (1 + {y \over N})^N.

We write y=\ln x, then

 x = \lim_{N\to \infty} (1+{\ln x\over N})^N

or rearranged

\ln x = \lim_{N\to \infty} {(x^{1/N} - 1) \over 1/N} = \lim_{n\to 0} {x^n-1\over n}.

A much easier and straighforward proof can be obtained by using the L'Hôpital's rule.

In fact, perhaps the easiest is to note that

 \lim_{n\to 0} {x^n-1\over n} = \left. \frac{\mathrm{d} x^n}{\mathrm{d} n} \right|_{n=0} = \left. x^n \ln x \right|_{n=0} = \ln x.

See also

The cavity method is an alternative method, often of simpler use than the replica method, for studying disordered mean field problems. It has been devised to deal with models on locally tree-like graphs

References

  • M Mezard, G Parisi & M Virasoro, "Spin Glass Theory and Beyond", World Scientific, 1987

Papers on Spin Glasses

  1. 1.0 1.1 Lua error in package.lua at line 80: module 'strict' not found.
  2. Lua error in package.lua at line 80: module 'strict' not found.

Books on Spin Glasses

  1. 1.0 1.1 Lua error in package.lua at line 80: module 'strict' not found.
  1. Lua error in package.lua at line 80: module 'strict' not found.