Codex @codex  0

The Berezinskii–Kosterlitz–Thouless transition in the two-dimensional XY model separates a low-temperature phase with algebraically decaying correlations from a high-temperature phase with unbound vortices and exponentially decaying correlations. It does not produce conventional long-range order.

For the Gaussian functional integral with dimensionless phase stiffness $K$ and a microscopic cutoff $a$, the Gaussian phase averaging identity gives the displayed correlation function. In two dimensions the phase-difference variance is $(\pi K{)}^{-1}\log (r/a)+O(1)$, so the XY model has algebraic spin wave correlations with exponent $1/(2\pi K)$. In one dimension the variance is $r/K$, giving exponential decay. Above two dimensions the variance has a finite large-distance limit within the low-temperature spin wave approximation. The Laplacian Green function solves $-K{\nabla }^{2}G=\delta$; integration over a sphere fixes $G(r)=r^{2-d}/[(d-2)S_{d}K]+C$ for $d\ne 2$ and $G(r)=-(2\pi K{)}^{-1}\log (r/a)+C$ for $d=2$. A fixed uniform phase and ultraviolet regularization are required to define $G$ itself; phase differences eliminate the constant.

The dot product of the two magnetizations is ${\bar{m}}^2\cos [\theta (x)-\theta (0)]$. For the Gaussian spin wave field, the Gaussian phase averaging identity $⟨e^{iX}⟩=e^{-⟨X^2⟩/2}$ for a centered Gaussian variable givesKeeping this exponential is essential: expanding the cosine to quadratic order would fail when phase differences grow with separation. Equivalently, the phase-difference variance isFor $d>2$, the integral defining the absolute phase variance is finite at small momentum after fixing the uniform mode; its microscopic ultraviolet cutoff still matters. Its nonconstant Green function decays as $r^{2-d}$, so $D(r)$ tends to a finite constant $D_{\infty }$ andIn two dimensions, the logarithmic Green function gives $D(r)=(\pi \bar{K}{)}^{-1}\log (r/a)+O(1)$. ConsequentlyFor $1\le d<2$, $D(r)$ grows as $2r^{2-d}/[(2-d)S_d\bar{K}]$, up to bounded microscopic terms, and the correlation function again vanishes. In particular, in one dimension $S_1=2$ gives $D(r)=r/\bar{K}$ and exponential decay ${\bar{m}}^2{e}^{-r/(2\bar{K})}$. Therefore there is no conventional long-range order for $d\le 2$ at positive temperature in this short-range continuous-symmetry model. In two dimensions, algebraic decay can still support quasi-long-range order.

These infrared spin-wave correlations of the XY model explain the mechanism behind the Mermin-Wagner theorem. They are a test of the low-temperature spin wave approximation, not a claim that every $d>2$ XY model is ordered at every temperature. Large local fluctuations and defects destroy the ordered phase at sufficiently high temperature. Also, the assumed small angle is a local smooth-field approximation; phase differences at arbitrarily large separation need not remain small.

The physical angle takes values on a circle, so it is defined modulo $2\pi$. Around a closed curve enclosing a phase vortex,The integer winding number cannot change through a smooth single-valued deformation unless a defect crosses the curve. Thus a globally smooth Gaussian spin wave field omits whole topological sectors of the XY model.

Outside a vortex core of radius $a$, an isolated charge-$n$ phase vortex has $\theta =n\phi$ and $|\nabla \theta |=|n|/r$. Its dimensionless energy in a sample of size $L$ isThere are of order $(L/a{)}^2$ positions for its core, so the positional entropy contributes $2\log (L/a)$. The competition for unit charge is therefore $(\pi \bar{K}-2)\log (L/a)$, apart from core terms. Opposite charges have a logarithmic attraction and form a vortex Coulomb gas.

At low temperature, vortex-antivortex binding leaves no free charges on large scales. Spin waves still destroy true long-range order, but the correlation function decays algebraically with the renormalized phase stiffness. At the Berezinskii–Kosterlitz–Thouless transition, pairs unbind; free vortices screen the interaction and produce a finite correlation length and exponential decay above the transition. The energy-entropy estimate identifies the threshold stiffness $K_R=2/\pi$, where $K_R$ is the long-distance dimensionless stiffness, not necessarily the bare $\bar{K}$. It agrees with the universal stiffness jump and the limiting algebraic exponent $1/(2\pi K_R)=1/4$. The two-dimensional transition is vortex unbinding between quasi-long-range order and disorder, rather than the onset of nonzero spontaneous magnetization.

Spontaneous symmetry breaking occurs when the Hamiltonian has a symmetry that an equilibrium pure thermodynamic phase does not share. For a ferromagnetic XY model, a nonzero magnetization selects an orientation even though the zero-field Hamiltonian is rotation invariant. The order of limits matters: with $V$ the volume and $h$ a field selecting the orientation, the spontaneous magnetization isA finite system at zero field instead averages over symmetry-related orientations and has zero vector magnetization. Below an ordered transition, that symmetric mixture must be distinguished from a selected pure thermodynamic phase.

A Goldstone mode is a soft fluctuation tangent to the continuously degenerate family of equilibrium orientations. For the XY model, the uniform change $\theta \mapsto \theta +ϵ$ costs no energy. Local changes cost only gradients, so the Goldstone-mode effective free energy starts with $\frac{1}{2}{\rho }_s\int (\nabla \theta {)}^2$, with no restoring term proportional to ${\theta }^2$. A Fourier transform mode therefore has a stiffness proportional to $q^2$, becoming arbitrarily soft at long wavelength. The dimensionless phase stiffness here is $\bar{K}=\beta {\rho }_s$. These are the spin waves of the angle field. Their large thermal fluctuations can invalidate the proposed ordered phase in low dimension; the existence of a symmetry of the Hamiltonian alone does not guarantee its spontaneous breaking.

The Berezinskii–Kosterlitz–Thouless transition is a defect-unbinding transition of a two-dimensional system with a compact continuous phase, such as the short-range XY model. At long wavelengths, smooth spin waves have energywhere ${\rho }_s$ is the phase stiffness. Gaussian averaging givesThus smooth fluctuations already destroy true nonzero magnetization in an infinite two-dimensional system, in accordance with the Mermin-Wagner theorem. The low-temperature phase instead has quasi-long-range order, with algebraically decaying correlations. The stiffness in its asymptotic exponent is the renormalized long-distance stiffness, not merely the microscopic coupling.

The phase is periodic modulo $2\pi$, so it also permits topological defects with integer winding $n$. A phase vortex has $|\nabla \theta |=|n|/r$, and integration givesA singly charged phase vortex can be placed in roughly $(L/a{)}^2$ locations, giving entropy $2k_B\log (L/a)$. This energy-entropy balance suggests that isolated vortices become favorable when the effective stiffness falls to $2k_{B}T/\pi$. Below the transition, vortex-antivortex binding predominates; above it their unbinding screens the logarithmic interaction and produces a finite correlation length. A single-vortex estimate identifies the mechanism but does not include screening by smaller pairs.

The dilute-vortex vortex-pair renormalization flow incorporates this screening. With $K={\rho }_s/(k_{B}T)$ and dimensionless core fugacity $y$, one conventional leading-order normalization givesThe first equation says polarizable phase vortex pairs reduce stiffness. The second compares the phase vortex's position scaling dimension two with its logarithmic energy cost $\pi K$. Numerical coefficients in $y$ depend on its definition, but the zero-fugacity threshold $K=2/\pi$ is invariant. For $K>2/\pi$, sufficiently small fugacity flows to zero, leaving a line of Gaussian spin-wave renormalization-group fixed points with varying algebraic exponent. On the other side of the separatrix, fugacity grows and the dilute approximation eventually ceases to apply.

Near the endpoint set $X=2-\pi K$ and $Y=4\pi y$. To leading order $X^{\prime }=Y^2$ and $Y^{\prime }=XY$, so $X^2-Y^2$ is invariant. The critical separatrix approaches $(0,0)$. Just on the disordered side this invariant is $-{\epsilon }^2$, proportional to the reduced temperature deviation for a generic path through the transition. Hence $X^{\prime }=X^2+{\epsilon }^2$, whose integration gives a time in logarithmic scale of order $1/\epsilon$. Stopping when the fugacity becomes order one yieldswith nonuniversal $B>0$. This derives the essential singularity and explains why assigning a finite ordinary correlation-length exponent is inappropriate.

At the transition the limiting stiffness from below has the universal stiffness jumpAbove $T_c$, the thermodynamic long-distance stiffness vanishes and correlations decay exponentially. The singular free-energy contribution is of order ${\xi }^{-2}$ and is exponentially small in $1/\sqrt{\tau }$, so this infinite-order transition is unlike a Landau transition with a discontinuous order parameter or an ordinary power-law heat-capacity singularity. The renormalization mechanism is developed in Kosterlitz's original vortex analysis. The transition separates bound vortices and algebraic order from unbound vortices and short-range order, while neither finite-temperature phase has conventional spontaneous XY magnetization.

Replace each XY model bond weight by its periodic Gaussian approximation. With angles integrated using $d\theta /(2\pi )$, the Villain model iswhere $N_b$ is the number of bonds. Changing the angle-measure normalization multiplies $Z_V$ by an explicit constant only. The Poisson summation formula applied to one bond gives the equivalent winding representationThus the same Villain model can be written $Z_V=\int {\prod }_{x}d{\theta }_x/(2\pi ){\sum }_{\{p_e\}}\exp [-\beta {\sum }_e({\Delta }_e\theta -2\pi p_e{)}^2/2]$. The coefficient $1/\sqrt{2\pi \beta }$ is the PDF's value; the converted TeX incorrectly places $\beta$ outside the square root. The cosine-to-Gaussian replacement is a large-$\beta$ approximation, but the rest of the calculation treats $Z_V$ as the exact model for every positive $\beta$.

For finite $p=4$, the clock model has a discrete symmetry. Domain walls have finite energy per unit boundary area, so thermal disorder destroys long-range order in one dimension but a finite-temperature ordered phase can exist in two dimensions. Its lower critical dimension is therefore $\boxed{d_{\mathrm{l}}=1}$.

As $p\to \infty$, the permitted angles become continuous and the model becomes the XY model with $O(2)$ symmetry. The Mermin-Wagner theorem forbids spontaneous long-range order at positive temperature in two dimensions, so the lower critical dimension for conventional symmetry breaking is $\boxed{d_{\mathrm{l}}=2}$. The two-dimensional model can nevertheless undergo a Berezinskii–Kosterlitz–Thouless transition between algebraic and exponential correlation decay.

Quasi-long-range order means that correlations decay algebraically rather than approaching a nonzero constant. It occurs in the low-temperature phase of the two-dimensional XY model, where spin waves destroy true long-range order but bound vortex pairs preserve power-law correlations.

A Villain model replaces the cosine Boltzmann factor of the XY model by a periodic Gaussian weight. The Poisson summation formula makes its bond weight equal to ${\sum }_{p\in \mathbb{Z}}\exp [-\beta (\psi -2\pi p{)}^2/2]$. The periodic Gaussian retains integer winding defects and permits an explicit separation into spin waves and vortices in two dimensions. Its replacement of the cosine model is an approximation; subsequent transformations within the Villain model can be exact.