Classical XY model

The classical XY model (sometimes also called classical rotor (rotator) model or O(2) model) is a model of statistical mechanics. It is the special case of the n-vector model for n = 2.

Definition

Given a D-dimensional lattice Λ, per each lattice site $j\in\Lambda$ there is a two-dimensional, unit-length vector $\mathbf{s}_j=(\cos \theta_j, \sin\theta_j)$. A spin configuration, $\mathbf{s}=(\mathbf{s}_j)_{j\in \Lambda}$ is an assignment of the angle $\theta_j\in (-\pi,\pi]$ per each site j in the lattice.

Given a translation-invariant interaction J_ij = J(i − j) and a point dependent external field $\mathbf{h}_{i}=(h_j,0)$, the configuration energy is

$H(\mathbf{s}) = - \sum_{i\neq j} J_{ij}\; \mathbf{s}_i\cdot\mathbf{s}_j -\sum_{j} \mathbf{h}_j\cdot \mathbf{s}_j =- \sum_{i\neq j} J_{ij}\; \cos(\theta_i-\theta_j) -\sum_{j} h_j\cos\theta_j$

The case in which J_i,j = 0 except for ij nearest neighbor is called 'nearest neighbor case'.

The configuration probability is given by the Boltzmann distribution with inverse temperature $\beta\ge 0$:

$P(\mathbf{s}) ={e^{-\beta H(\mathbf{s})} \over Z} \qquad Z=\int_{[-\pi,\pi]^{\Lambda}}\!\prod_{j\in \Lambda}d\theta_j\;e^{-\beta H(\mathbf{s})} \,.$

where Z is the normalization, or partition function.^[1]

General properties

• The existence of the thermodynamic limit for the free energy and spin correlations were proved by Ginibre, extending to this case the Griffiths inequality.^[2]

• At high temperature, the spontaneous magnetization vanishes:

$M(\beta):=|\langle \mathbf{s}_i\rangle|=0$

Besides, cluster expansion shows that the spin correlations cluster exponentially fast: for instance

$|\langle \mathbf{s}_i\cdot \mathbf{s}_j\rangle| \le C(\beta)e^{-c(\beta)|i-j|}$

• Using the Griffiths inequality in the formulation of Ginibre, Aizenman and Simon^[3] proved that the two point spin correlation of the ferromagnetics XY model in dimension D, coupling J > 0 and inverse temperature β is dominated by (i.e. has upper bound given by) the two point correlation of the ferromagnetic Ising model in dimension D, coupling J > 0and inverse temperature β / 2

$\langle \mathbf{s}_i\cdot \mathbf{s}_j\rangle_{J,2\beta} \le \langle \sigma_i\sigma_j\rangle_{J,\beta}$

Hence the critical β of the XY model cannot be smaller than the double of the critical temperature of the Ising model

$\beta_c^{XY}\ge 2\beta_c^{\rm Is}$

One dimension

As in any 'nearest-neighbor' n-vector model with free boundary conditions, if the external field is zero, there exists a simple exact solution.

In the nearest-neighbor, free boundary conditions case, the Hamiltonian is

$H(\mathbf{s}) = - J [\cos(\theta_1-\theta_2)+\cdots+\cos(\theta_{L-1}-\theta_L)]$

therefore the partition function factorizes under the change of coordinates

$\theta_j=\theta_j'+\theta_{j-1} \qquad j\ge 2$

That gives

$\begin{align} Z&=\int_{-\pi}^{\pi}d\theta_1\cdots d\theta_L\; e^{\beta J \cos(\theta_1-\theta_2)}\cdots e^{\beta J \cos(\theta_{L-1}-\theta_L)} =2\pi \prod_{j=2}^L\int_{-\pi}^{\pi}d\theta'_j\;e^{\beta J \cos\theta'_j}= \\ &=2\pi\left[\int_{-\pi}^{\pi}d\theta'_j\;e^{\beta J \cos\theta'_j}\right]^{L-1} \end{align}$

Finally

$f(\beta, 0)=-\frac{1}{\beta}\ln \int_{-\pi}^{\pi}d\theta'_j\;e^{\beta J \cos\theta'_j}$

The same computation for periodic boundary condition (and still h = 0) requires the transfer matrix formalism.^[4]

Two Dimensions

• At low temperature, the spontaneous magnetization remains zero,

$M(\beta):=|\langle \mathbf{s}_i\rangle|=0$

but Fröhlich and Spencer proved that the decay of the correlations is only power law.^[5]

$|\langle \mathbf{s}_i\cdot \mathbf{s}_j\rangle| \ge\frac{C(\beta)}{1+|i-j|^{\eta(\beta)}}$

(a power law upper bound was found by McBryan and Spencer).

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 leads to a vortex-unbinding transition from the low-temperature phase to the high-temperature disordered phase. In two spatial dimensions the XY model exhibits a Kosterlitz-Thouless transition from the disordered high-temperature phase into the quasi-long range ordered low-temperature phase.

Three and Higher Dimensions

• At low temperature, infrared bound shows that the spontaneous magnetization is strictly positive:

$M(\beta):=|\langle \mathbf{s}_i\rangle|>0$

See also

• Goldstone boson
• Ising model
• Potts model
• Kosterlitz-Thouless transition
• Topological defect
• Superfluid film

References

1. ^ Lubensky, Chaikin (2000). Principles of Condensed Matter Physics. Cambridge University Press. pp. 699. ISBN 0521794501. http://books.google.com/books?id=P9YjNjzr9OIC&printsec=frontcover&cad=0#v=onepage&q&f=false. 
2. ^ Ginibre, J. (1970). "General formulation of Griffiths' inequalities". Comm. Math. Phys. 16 (4): 310–328. doi:10.1007/BF01646537. 
3. ^ Aizenman, M.; Simon, B. (1980). "A comparison of plane rotor and Ising models". Phys. Letts. A 76. doi:10.1016/0375-9601(80)90493-4. http://www.sciencedirect.com/science/article/pii/0375960180904934. 
4. ^ Mattis, D.C. (1984). "Transfer matrix in plane-rotator model". Phys.Letts 104 A. doi:10.1016/0375-9601(84)90816-8. http://www.sciencedirect.com/science/article/pii/0375960184908168. 
5. ^ Fröhlich, J.; Spencer, T. (1981). "The Kosterlitz–Thouless transition in two-dimensional abelian spin systems and the Coulomb gas". Comm. Math. Phys. 81 (4): 527–602. http://projecteuclid.org/euclid.cmp/1103920388. 

• Evgeny Demidov, Vortices in the XY model (2004)