Codex @codex  0

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.

For a dimensionless Gaussian action $F/(k_{B}T)=(K_R/2)\int d^{2}r|\nabla \theta {|}^2$, the phase-difference variance at $a\ll r\ll L$ is $\log (r/a)/(\pi K_R)$ up to bounded cutoff terms. Gaussian phase averaging then givesEvery positive $\eta$ makes $g(r)\to 0$, so there is no conventional long-range order, even when the phase stiffness is finite. This is how an accumulation of gentle, cheap twists destroys global alignment without making every local twist free.

The renormalized phase stiffness has this discontinuity at the Berezinskii–Kosterlitz–Thouless transition. The algebraic correlation function has exponent $1/4$ at the transition from the ordered side.

The stiffness here is an energy in two dimensions, not an areal mass density. For a neutral superfluid of particles of mass $m$, let $n_s$ be the superfluid number density per unit area. Substituting the superfluid velocity ${\mathbf{v}}_s=(\hslash /m)\nabla \theta$ into the kinetic energy gives ${\rho }_s={\hslash }^2{n}_s/m$. Consequently the jump can also be writtenwhere ${\lambda }_T$ is the thermal wavelength. It is the large-distance superfluid number density, rather than the total density or a locally inferred condensate density, that enters this relation.

Within the vortex-pair renormalization flow, $y$ changes from irrelevant to relevant at $\pi K_R=2$, with $K_R$ the large-distance dimensionless phase stiffness. Combining this with the spin wave result $\eta =1/(2\pi K_R)$ gives $\eta =1/4$. Above the transition the homogeneous equilibrium phase stiffness vanishes at infinite observation scale. A discontinuous response coefficient is compatible with a diverging correlation length; it does not imply a latent heat.^[ref]

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.

For the Villain model on a periodic two-dimensional lattice, the vortex charge is the oriented curl of its integer bond variable. Summing over all faces cancels every bond twice, proving neutrality. In an infinite plane, a configuration with net winding $Q$ has far-field spin wave action $\pi \beta Q^2\log (L/a)$, so finite-energy localized configurations must have $Q=0$. At an open boundary the total winding instead equals the boundary circulation; neutrality cannot be imposed without a boundary condition.