Classical XY model on arvau

Wikipedia / Classical XY model

The classical XY model (sometimes also called classical rotor (rotator) model or O(2) model) is a lattice model of statistical mechanics. In general, the XY model can be seen as a specialization of Stanley's n-vector model for .

Definition

Given a -dimensional lattice , at each lattice site there is a two-dimensional, unit-length vector . A spin configuration, is an assignment of the angle for each .

Given an interaction (often assumed to be translation-invariant, so , the configuration energy of a given spin configuration is

Often one considers the nearest neighbor case, in which except when and are nearest neighbors.

The configuration probability is given by the Boltzmann distribution with inverse temperature :

where is the normalization, or partition function. indicates the expectation of the random variable in the infinite volume limit, after periodic boundary conditions have been imposed.

The classical XY model is the semiclassical limit of the spin- quantum XY model in the limit as becomes large.

Rigorous results

The existence of the thermodynamic limit for the free energy and spin correlations were proved by Ginibre, extending to this case the Griffiths inequality.
Using the Griffiths inequality in the formulation of Ginibre, Aizenman and Simon proved that the two point spin correlation of the ferromagnetics XY model in dimension , coupling and inverse temperature is dominated by (i.e. has an upper bound given by) the two point correlation of the ferromagnetic Ising model in dimension , coupling and inverse temperature Hence the critical of the XY model cannot be smaller than the double of the critical

One dimension

As in any 'nearest-neighbor' n-vector model is the modified Bessel function), the free energy Using the properties of the modified Bessel functions, the specific heat (per spin) can be expressed as is the short-range correlation function,

Even in the thermodynamic limit, there is no divergence in the specific heat. Indeed, like the one-dimensional Ising model, the one-dimensional XY model has no phase transitions at finite temperature.

The same computation for periodic boundary condition (and still ) requires the transfer matrix formalism, though the result is the same.

is small enough that it can be treated as a perturbation to the system in zero-field, then the magnetic susceptibility can be estimated. This is done by using the eigenstates computed by the transfer matrix approach and computing the energy shift with second-order perturbation theory. One finds is the Curie constant.

Two dimensions

The two-dimensional XY model with nearest-neighbor interactions is an example of a two-dimensional system with continuous symmetry that does not have long-range order as required by the Mermin–Wagner theorem. Likewise, there is not a conventional phase transition present that would be associated with symmetry breaking. However, as will be discussed later, the system does show signs of a transition from a disordered high-temperature state to a quasi-ordered state below some critical temperature, called the Kosterlitz-Thouless transition which occurs at low temperatures. For example, Mattis (1984 The 2D XY model has also been studied in great detail using Monte Carlo simulations, for example with the Metropolis algorithm (often, it can be discretized into finitely-many angles, like in the related Potts model, the Boltzmann factor for the energy change. The Monte Carlo method has been used to verify, with various methods, the critical temperature of the system, and is estimated to be. The Monte Carlo method can also compute average values that are used to compute thermodynamic quantities like magnetization, spin-spin correlation, correlation lengths, and specific heat. These are important ways to characterize the behavior of the system near the critical temperature. The magnetization and squared magnetization, for example, can be computed as are the number of spins. The mean magnetization characterizes the magnitude of the net magnetic moment of the system; in many magnetic systems this is zero above a critical temperature and becomes non-zero spontaneously at low temperatures. Similarly the mean-squared magnetization characterizes the average of the square of net components of the spins across the lattice. Either of these are commonly used to characterize the order parameter of a system. Rigorous analysis of the XY model shows the magnetization in the thermodynamic limit is zero, and that the square magnetization approximately follows, which vanishes in the thermodynamic limit. Indeed, at high temperatures this quantity approaches zero since the components of the spins will tend to be randomized and thus sum to zero. However at low temperatures for a finite system, the mean-square magnetization increases, suggesting there are regions of the spin space that are aligned to contribute to a non-zero contribution. The magnetization shown (for a 25x25 lattice) is one example of this, that appears to suggest a phase transition, while no such transition exists in the thermodynamic limit.

. There is no feature in the specific heat consistent with critical behavior (like a divergence) at this predicted temperature. Indeed, estimating the critical temperature comes from other methods, like from the helicity modulus, or the temperature dependence of the divergence of susceptibility.. This peak position and height have been shown not to depend on system size, for lattices of linear size greater than 256; indeed, the specific heat anomaly remains rounded and finite for increasing lattice size, with no divergent peak.

must vary smoothly over changes in position. Expanding the original cosine as a Taylor series The continuous version of the XY model is often used to model systems that possess order parameters with the same kinds of symmetry, e.g. superfluid helium, hexatic liquid crystals. This is what makes them peculiar from other phase transitions which are always accompanied with a symmetry breaking. Topological defects in the XY model lead to a vortex-unbinding transition from the low-temperature phase to the high-temperature disordered phase. Indeed, the fact that at high temperature correlations decay exponentially fast, while at low temperatures decay with power law, even though in both regimes , is called Kosterlitz–Thouless transition. Kosterlitz and Thouless provided a simple argument of why this would be the case: this considers the ground state is an effective length scale. This indicates that at low temperatures, any vortices that arise will want to annihilate with antivortices to lower the system energy. Indeed, this will be the case qualitatively if one watches 'snapshots' of the spin system at low temperatures, where vortices and antivortices gradually come together to annihilate. Thus, the low-temperature state will consist of bound vortex-antivortex pairs. Meanwhile, at high temperatures, there will be a collection of unbound vortices and antivortices that are free to move about the plane.

To visualize the Ising model, one can use an arrow pointing up or down, or represented as a point colored black/white to indicate its state. To visualize the XY spin system, the spins can be represented as an arrow pointing in some direction, or as being represented as a point with some color. Here it is necessary to represent the spin with a spectrum of colors due to each of the possible continuous variables. This can be done using, for example, a continuous and periodic red-green-blue spectrum. As shown on the figure, cyan corresponds to a zero angle (pointing to the right), whereas red corresponds to a 180 degree angle (pointing to the left). One can then study snapshots of the spin configurations at different temperatures to elucidate what happens above and below the critical temperature of the XY model. At high temperatures, the spins will not have a preferred orientation and there will be unpredictable variation of angles between neighboring spins, as there will be no preferred energetically favorable configuration. In this case, the color map will look highly pixellated. Meanwhile, at low temperatures, one possible ground-state configuration has all spins pointed in the same orientation (same angle); these would correspond to regions (domains) of the color map where all spins have roughly the same color.

corresponds to a vortex (or antivortex). These vortexes are topologically non-trivial objects that come in vortex-antivortex pairs, which can separate or pair-annihilate. In the colormap, these defects can be identified in regions where there is a large color gradient where all colors of the spectrum meet around a point. Qualitatively, these defects can look like inward- or outward-pointing sources of flow, or whirlpools of spins that collectively clockwise or counterclockwise, or hyperbolic-looking features with some spins pointing toward and some spins pointing away from the defect. As the configuration is studied at long time scales and at low temperatures, it is observed that many of these vortex-antivortex pairs get closer together and eventually pair-annihilate. It is only at high temperatures that these vortices and antivortices are liberated and unbind from one another.

In the continuous XY model, the high-temperature spontaneous magnetization Besides, cluster expansion At low temperatures, i.e. , the spontaneous magnetization remains zero (see the Mermin–Wagner theorem but the decay of the correlations is only power law: Fröhlich and Spencer found the lower bound

Three and higher dimensions

Independently of the range of the interaction, at low enough temperature the magnetization is positive.

. Besides, cluster expansion.
At low temperature, infrared bound but, conjecturally, in each of these extremal states the truncated correlations decay algebraically.

Phase transition

As mentioned above in one dimension the XY model does not have a phase transition, while in two dimensions it has the Berezinski-Kosterlitz-Thouless transition between the phases with exponentially and powerlaw decaying correlation functions.

In three and higher dimensions the XY model has a ferromagnet-paramagnet phase transition. At low temperatures the spontaneous magnetization is nonzero: this is the ferromagnetic phase. As the temperature is increased, spontaneous magnetization gradually decreases and vanishes at a critical temperature. It remains zero at all higher temperatures: this is the paramagnetic phase.

In four and higher dimensions the phase transition has mean field theory critical exponents (with logarithmic corrections in four dimensions).

Three dimensional case: the critical exponents

The three dimensional case is interesting because the critical exponents at the phase transition are nontrivial. Many three-dimensional physical systems belong to the same universality class as the three dimensional XY model and share the same critical exponents, most notably easy-plane magnets and liquid Helium-4. The values of these critical exponents are measured by experiments, Monte Carlo simulations, and can also be computed by theoretical methods of quantum field theory, such as the renormalization group and the conformal bootstrap. Renormalization group methods are applicable because the critical point of the XY model is believed to be described by a renormalization group fixed point. Conformal bootstrap methods are applicable because it is also believed to be a unitary three dimensional conformal field theory.

Most important critical exponents in the Ginzburg–Landau. According to a conformal bootstrap computation, these three dimensions are given by:This gives the following values of the critical exponents:Monte Carlo methods give compatible determinations:.

Notes

References

Evgeny Demidov, Vortices in the XY model (2004)