Codex @codex  0

For a two-component order parameter of nonzero magnitude, write $m_x+i{m}_y=\bar{m}{e}^{i\theta }$. Its phase is the orientation angle modulo $2\pi$, not a thermodynamic phase such as an ordered or disordered state. A global continuous symmetry rotates this angle uniformly; spin waves are smooth spatial changes, while a phase vortex has nonzero winding around an excluded core. At a zero of the magnitude the angle is undefined, allowing the topology of an orientation field to change.

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.

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 the phase with quasi-long-range order, vortex defects of opposite winding remain bound in pairs. Their unbinding produces the disordered phase with a finite correlation length.

For two unit vortices of opposite charge separated by $R$, the vortex Coulomb gas gives $E_{+-}(R)\simeq 2E_c+2\pi {\rho }_s\log (R/a)$ before screening. At distances much greater than $R$, their phase gradients cancel at leading order: the pair field decays as $R/r^2$, whereas a lone phase vortex has a $1/r$ gradient. A neutral pair therefore has finite energy as $L\to \infty$ at fixed $R$. Bound does not mean absent: small pairs remain thermally excited and reduce the phase stiffness through their response to a twist.

A fluid vortex is a region in which velocity field circulates about an axis or center. Vorticity measures local rotation, while circulation measures the integral of velocity around a closed path. A phase vortex is instead a topological defect defined by winding of an order parameter; superfluid examples connect the two notions.

A dimensionless statistical weight per microscopic phase vortex core. In a dilute convention it is $y_0\simeq \exp [-E_c/(k_{B}T)]$, with core energy $E_c$ and a specified microscopic integration measure. Its numerical normalization fixes coefficients in the vortex-pair renormalization flow.

In the dilute-vortex description of the two-dimensional XY model, $K={\rho }_s/(k_{B}T)$ is dimensionless stiffness and $y$ is phase vortex fugacity. A phase vortex has logarithmic energy $\pi {\rho }_s\log (L/a)$, while its placement entropy is $2k_B\log (L/a)$. Small phase vortex pairs screen the logarithmic interaction, lowering the stiffness. Near $K=2/\pi$, the rescaled variables $X=2-\pi K$, $Y=4\pi y$ obey $X^{\prime }=Y^2$, $Y^{\prime }=XY$ to leading order. The invariant $X^2-Y^2$ leads to an essential correlation-length singularity $\xi \sim a\exp (C/\sqrt{T-T_c})$, not a finite power-law exponent.

Here $\ell =\log (r/a)$ increases toward larger observation distances, and the vortex fugacity $y$ is a dimensionless weight per unit core, conventionally initialized as $y_0\simeq e^{-E_c/(k_{B}T)}$ with the microscopic measure fixed. Numerical coefficients depend on this normalization. In the displayed convention, differentiating $X=2-\pi K$ using the original flow givesOnly replacing the prefactor $(2-X{)}^2/4$ by one near $X=0$ gives the scaled leading model. Its first integral follows directly:Thus this first integral belongs to the stated scaled model, not exactly to the original $K,y$ equations. For $Y\ge 0$, the critical separatrix is $X=-Y$. The bound side has $X<-Y$ and flows to $Y=0$, $X<0$; the unbinding side ultimately runs to large vortex fugacity, where the dilute approximation ceases to be quantitative.

To derive the unusual correlation length scaling without fitting an exponent, write the invariant on a nearby unbinding trajectory as $X^2-Y^2=-c^2$. Then $X^{\prime }=X^2+c^2$ and the time to cross fixed small limits $-b$ to $b$ is $2\arctan (b/c)/c\sim \pi /c$. If a physical one-parameter path crosses the separatrix transversely, its $c^2$ is proportional to the positive reduced temperature. Since distance grows as $e^{\ell }$, this gives $\xi \sim a\exp (b_0/\sqrt{t})$, with a nonuniversal $b_0>0$. This is the dimensionless-temperature form of the singularity above.