Codex @codex  0

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.

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.

Assume $g>0$ for stability. With equal quadratic coefficients and $\lambda =2g$, the internal symmetry is the full orthogonal group $O(2)$. In two dimensions at positive temperature, with short-range interactions, the Mermin-Wagner theorem forbids spontaneous symmetry breaking of this continuous symmetry. Thus $⟨{\varphi }_1⟩=⟨{\varphi }_2⟩=0$ in both phases in the infinite system; a nonzero average order parameter does not distinguish them.

Instead measure the correlation functionThe low-temperature phase has quasi-long-range order and algebraic correlations; the high-temperature phase has a finite correlation length and exponential correlations. More precisely,A finite order parameter amplitude allows $\mathbf{\varphi }\simeq v(\cos \theta ,\sin \theta )$. Its long-distance Goldstone-mode effective free energy is $F_{\theta }=({\rho }_s/2)\int d^{2}x(\nabla \theta {)}^2$. Restoring the Boltzmann constant and using the renormalized phase stiffness, the thermal phase fluctuations giveThe Berezinskii–Kosterlitz–Thouless transition separates vortex-antivortex binding from unbound vortex defects. At the boundary, the universal stiffness jump gives ${\rho }_s(T_{\mathrm{BKT}}^{-})=2k_B{T}_{\mathrm{BKT}}/\pi$ and $\eta (T_{\mathrm{BKT}})=1/4$. The low-temperature phase stiffness is nonzero; it vanishes at long distances in the disordered phase. The exact transition location is fluctuation dependent and cannot be obtained just by setting the bare ${\mu }^2$ to zero.

Consider an original counting example: a square two-dimensional XY model of side $L$, a microscopic core radius $a$, and a fixed phase stiffness ${\rho }_s$. This deliberately neglects screening by other vortices. Away from a unit phase vortex, $\theta =\arg (x+iy)$ and $|\nabla \theta |=1/r$. Integrating the Goldstone-mode effective free energy over an annulus givesBoundary shape changes finite terms. There are order $(L/a{)}^2$ core positions, so the positional entropy is $S_{\mathrm{v}}=2k_B\log (L/a)+O(k_B)$, where $k_B$ is the Boltzmann constant. With $K={\rho }_s/(k_{B}T)$, subtracting $T{S}_{\mathrm{v}}$ gives the displayed free energy.

For $L/a=1000$ and $E_c/(k_{B}T)=2$, the dimensionless free energy cost is about $5.546$ at $K=0.8$ and $-0.965$ at $K=0.5$. The same logarithm appears in the cost and the positional reward: increasing the system size strengthens whichever has the larger coefficient. The corresponding one-defect Boltzmann factors, about $0.00390$ and $2.62$, are relative statistical weights, not normalized probabilities. The latter signals failure of a dilute single-defect estimate.

A single net phase vortex requires compatible boundaries: on a periodic sample use neutral pairs, as required by vortex charge neutrality. This argument identifies a possible threshold, but predicting the universal stiffness jump requires the renormalized phase stiffness, after smaller pairs have responded. It does not justify setting the transition temperature equal to a bare lattice coupling divided by a universal constant.

In 2006, Zoran Hadzibabic, Peter Krüger, Marc Cheneau, Baptiste Battelier and Jean Dalibard split a cooled rubidium-87 cloud into two independent planar gases using an optical lattice. After release, the expanding sheets overlapped; absorption imaging recorded their matter-wave interference. Smooth fringe waviness revealed thermal phase fluctuations, while sharp dislocations indicated isolated vortices. The initial evaporative-cooling endpoint controlled the preparation; central fringe contrast served as a degeneracy proxy because direct thermometry in the lattice was difficult.^[ref]

Loss of extended coherence coincided with increased vortex signatures. The optical trap and finite cloud produced a crossover; tightly bound pairs were unresolved, and short-wavelength modes could contaminate high-temperature counts. These were observations of the vortex-unbinding mechanism, not a direct measurement of phase stiffness or a sharp homogeneous universal stiffness jump. The authors left the sample's superfluidity open.^[ref]