---
title: 'Active XY Model: Nonequilibrium U(1) Theories'
url: https://www.emergentmind.com/topics/active-xy-model
type: topic
---

# Active XY Model: Nonequilibrium U(1) Theories

In the literature, the term **active XY model** refers not to a single canonical Hamiltonian but to a family of nonequilibrium extensions of the classical XY model in which the local \(U(1)\) angle field is driven by persistent torques, intrinsic frequencies, self-propelled particle motion, coupling to conserved diffusive species, or coupling to flow fields. Across these variants, the common structural element is an XY phase or spin variable \(\theta\), while the central physical question is how activity modifies Berezinskii–Kosterlitz–Thouless (BKT) physics, quasi-long-range order (QLRO), and the dynamics of topological defects [2401.11175][2103.12578][2507.17108][2104.10982][1911.00214][1107.0370].

## 1. Equilibrium reference model and defect taxonomy

The equilibrium reference point is the two-dimensional XY model with planar spins \(\mathbf{S}_i=(\cos\theta_i,\sin\theta_i)\) and nearest-neighbor Hamiltonian
\[
H_{\text{XY}} = -J\sum_{\langle ij\rangle}\cos(\theta_i-\theta_j).
\]
In two dimensions, the low-temperature phase has QLRO rather than true long-range order: the spin correlation decays algebraically, \(C(r)\sim r^{-\eta(T)}\), and the BKT transition occurs at \(\eta(T_{\mathrm{KT}})=1/4\). The topological defects are integer vortices, defined by a winding number \(m\in\mathbb{Z}\) around a loop [2307.13593].

A useful extension is the generalized XY model with a nematic interaction of period \(q=2\),
\[
H = - \left( \sum_{\langle i j \rangle} \Delta \cos(\theta_i - \theta_j)
+ \sum_{\langle i j \rangle} (1 - \Delta) \cos\big(2(\theta_i - \theta_j)\big) \right),
\]
which interpolates between the usual vector XY limit \(\Delta=1\) and a purely nematic model \(\Delta=0\). This model supports both integer and half-integer vortices and exhibits ferromagnetic, nematic, and disordered phases separated by BKT-like and Ising-like transition lines [2404.17504]. That defect taxonomy is directly relevant to active polar and active nematic generalizations, because activity changes defect energetics and kinetics before it changes the underlying topological classification.

## 2. Principal nonequilibrium formulations

One common construction is the **persistent-noise XY model**, in which passive XY spins are driven by an Ornstein–Uhlenbeck angular velocity. On a triangular lattice, the microscopic dynamics can be written as
\[
\dot{\theta}_i = -\partial_{\theta_i}U + \varpi_i,\qquad
\tau_0 \dot{\varpi}_i = - \varpi_i + \sqrt{2\tau_0 T}\,\xi_i,
\]
with exponentially correlated torque noise
\[
\langle \varpi_i(t)\varpi_j(t')\rangle = \delta_{ij}\, T\, e^{-|t-t'|/\tau_0}.
\]
For \(\tau_0\to 0\), the model reduces to the equilibrium XY Langevin dynamics, whereas any finite \(\tau_0\) breaks detailed balance [2602.22028].

A second formulation is the **self-driven rotor** or short-range noisy Kuramoto version,
\[
\dot{\theta}_i = \Omega_i + J \sum_{j\in\partial_i}\sin(\theta_j-\theta_i) + \sqrt{2D}\,\nu_i(t),
\]
with intrinsic frequencies \(\Omega_i\) drawn from a Gaussian of variance \(s^2\). After rescaling, the relevant activity parameter is \(\sigma=s/J\). Here the nonequilibrium drive is the sitewise energy injection by the intrinsic frequencies, which persists even at \(T\to 0\) [2103.12578].

A third class uses **mobile, self-propelled XY particles** on a lattice with exclusion and vacancies. The spin sector is governed by
\[
E = - J \sum_{\langle i,j\rangle} X_i\,X_j \,\cos(\theta_i - \theta_j),
\]
while particle motion is biased along the local spin direction by a self-propulsion parameter \(\epsilon\in[0,1]\). At \(\rho=1\) the model reduces to the standard XY model; for \(\rho<1\), it becomes a diluted XY system with spin-dependent motility [2507.17108].

A fourth formulation couples the XY phase to additional hydrodynamic fields. For active spinners on a substrate with a conserved diffusing species \(c(\mathbf{x},t)\), the hydrodynamic equations are
\[
\frac{\partial\theta}{\partial t} = \kappa\nabla^2\theta + \Omega_1 \delta c + \frac{\lambda}{2} (\boldsymbol\nabla\theta)^2 + \Omega_2 (\delta c)^2 + f_\theta,
\]
\[
\frac{\partial\delta c}{\partial t} = \lambda_0 \tilde\Omega_0 \nabla^2\theta + D\nabla^2\delta c + \lambda_0\tilde\Omega_1\,\boldsymbol\nabla\cdot(\delta c\,\boldsymbol\nabla\theta) + f_c,
\]
where the active terms violate detailed balance [2104.10982]. Closely related hydrodynamic theories couple \(\phi\) directly to a two-dimensional flow field \(\mathbf{v}\), with advection, vorticity, and compressibility couplings in the phase equation and an active Stokes equation for the flow [1911.00214].

A fifth formulation is the uniformly driven three-dimensional rotor model,
\[
d\varphi_x(t)= d\,dt - \sum_{y\sim x}\sin(\varphi_x(t)-\varphi_y(t))\,dt + \sqrt{2\beta^{-1}}\,dW_x(t),
\]
which represents a spatially extended XY medium under a uniform non-conservative torque [1107.0370].

## 3. Ordering phenomena under activity

Activity does not modify XY ordering in a universal way. In the persistent-noise model, the long-wavelength spin-wave sector remains Gaussian but acquires an effective large-scale temperature \(T_{\rm eff}^{\rm (spinwave)}=\tau_0T\). The equal-time spectrum becomes
\[
\langle \tilde{\Theta}(\mathbf{k}) \tilde{\Theta}(-\mathbf{k}) \rangle
= \frac{1}{1 + \kappa \tau_0 \ell_0^2 k^2}\,\frac{\tau_0 T/\kappa}{k^2},
\]
so that
\[
g_\theta(r)\sim \left(\frac{r}{\ell_\tau}\right)^{-\eta},\qquad
\eta=\frac{\tau_0 T}{2\pi\kappa}.
\]
This produces QLRO with exponents that grow linearly with persistence and can exceed the equilibrium BKT value \(1/4\) without immediate defect proliferation [2401.11175]. A later analysis concluded that the order–disorder transition remains of BKT type, with \(\nu=1/2\) and \(z\simeq2\), but with a persistence-dependent critical exponent \(\eta_c\) that increases with \(\tau_0\); for example, the reported values are \(\eta_c=0.232(2)\) at \(\tau_0=0\), \(0.342(4)\) at \(\tau_0=6\), \(0.43(1)\) at \(\tau_0=12\), and \(0.54(1)\) at \(\tau_0=24\) [2602.22028]. This resolves an earlier open question only partially: the essential singularities remain BKT-like, but the equilibrium relation \(\eta_c=1/4\) is lost.

The self-driven rotor model gives the opposite outcome. For any \(\sigma>0\), equilibrium QLRO is destroyed: at low temperature the correlation function becomes exponential,
\[
C(r)\sim e^{-r/\xi},
\]
and the system breaks into a mosaic of ordered domains with \(\xi\propto \sigma^{-1}\). A finite-size crossover line survives in the density of free vortices,
\[
T_c(\sigma)=T_{KT}(1-a\,\sigma\ln L),
\]
with \(a\simeq0.43\), but in the thermodynamic limit the true BKT transition survives only at \(\sigma=0\) [2103.12578]. The contrast with persistent-noise driving is structural: intrinsic-frequency disorder directly destroys the global phase field, whereas persistent Ornstein–Uhlenbeck forcing separates large-scale spin-wave fluctuations from local defect nucleation.

Hydrodynamic active XY theories on substrates and in flows produce additional ordered regimes. For active XY spins coupled to a conserved diffusing species, the stable phase-ordered states exhibit generalized QLRO in which
\[
\Delta_\theta^R \sim \big[\ln(L/a_0)\big]^{1-\eta},
\]
with nonuniversal \(\eta\). Depending on the nonlinear couplings, the ordered state can be stronger than equilibrium QLRO or weaker than equilibrium QLRO, with a corresponding suppression or enhancement of density fluctuations [2104.10982]. Flow-coupled live XY spins show an analogous classification: in the Malthusian case, sufficiently strong active flow–phase couplings can yield QLRO, stronger-than-QLRO, or weaker-than-QLRO, whereas other parameter choices lead to disorder of the same qualitative type as immobile active XY spins or the two-dimensional Kardar–Parisi–Zhang equation [1911.00214].

In three dimensions, a uniform torque produces yet another pattern. The driven XY model has essentially a unique translation-invariant stationary distribution for almost all temperatures, but at low temperature it also supports rotating states in which the ordered phase does not relax pointwise in time: the macroscopic magnetization simply rotates at the imposed drift rate [1107.0370].

## 4. Topological defects and defect-mediated dynamics

Topological defects remain central in all active XY variants, but their role becomes strongly model dependent. In equilibrium and generalized equilibrium XY models, vortices are classified by the topology of the angle field: integer vortices in the standard model, and both integer and half-integer vortices in the \(q=2\) generalized model [2307.13593].

In the self-driven rotor model, activity qualitatively changes defect motion. At equilibrium, an isolated vortex has a Gaussian displacement distribution and a mean-squared displacement
\[
\langle [\mathbf{r}(t)-\mathbf{r}(0)]^2\rangle = 4\mathcal{D}\,\frac{t}{\ln t}.
\]
For \(\sigma>0\), the displacement distribution remains Gaussian but the defects become superdiffusive,
\[
\langle [\mathbf{r}(t)-\mathbf{r}(0)]^2\rangle
= 4\mathcal{D}\,(\sigma t)^{3/2},
\]
and vortex–antivortex pairs generically unbind at any \(T>0\) once \(\sigma>0\). Their motion is guided by domain boundaries created by the active mosaic [2103.12578].

In the lattice gas of self-propelled XY particles, defects become direct organizers of phase separation. Vortex charge is measured by the plaquette winding
\[
m = \dfrac{1}{2\pi} \oint_{\square} \nabla \theta \cdot  \mathrm{d}\boldsymbol{\ell},
\]
and the total vorticity is \(N=N_{+1}-N_{-1}\). Positive-charge vortices can act as aggregation centers: particles accumulate into dense clusters around specific \(+1\) defects, whereas \(-1\) defects dissipate. Cluster growth proceeds by defect-mediated merging, and the late-stage relaxation follows a two-stage exponential form,
\[
\phi_{\mathrm{max}}(t)=\phi_{\mathrm{max}}(\infty)-A_1 e^{-t/\tau_1}-A_2 e^{-t/\tau_2},
\]
with a characteristic timescale
\[
\tau \sim L^3.
\]
The resulting motility-induced phase separation is therefore explicitly defect mediated [2507.17108].

Persistent-noise XY models modify defect physics less dramatically at low temperature. Their key mechanism is that large-scale spin waves are controlled by \(\tau T\) while local defect nucleation remains controlled by \(T\). That separation permits very large algebraic exponents without immediate vortex proliferation [2401.11175]. In the later finite-\(\tau_0\) analysis, vortices still drive the transition, but the critical fugacity–stiffness balance is shifted by persistence, producing a BKT-like transition with non-universal \(\eta_c(\tau_0)\) [2602.22028].

## 5. Nonequilibrium field theory and effective temperatures

A recurring theme in active XY models is the breakdown of any single-temperature equilibrium description. In the persistent-noise formulation, the exact spin-wave stationary functional is
\[
P[\Theta] \propto \exp\left[ -\frac{\kappa}{2\tau_0 T} \int d^2 r\,|\nabla\Theta|^2 - \tau_0 \kappa \ell_0^2 \int d^2 r\,|\nabla^2\Theta|^2 \right].
\]
The first term is equilibrium-like and identifies the large-scale effective temperature \(\tau_0T\); the second is a nonequilibrium ultraviolet regularization that generates the \(1/k^4\) tail at high wave number and has no equilibrium counterpart [2602.22028]. Earlier coarse graining reached the same qualitative conclusion in a slightly different notation: ultraviolet fluctuations remain controlled by \(T\), while long-wavelength spin waves are controlled by \(T_S=\tau T\), so local defect physics and large-scale orientational fluctuations decouple [2401.11175].

The substrate and flow-coupled theories make the nonequilibrium structure explicit at the hydrodynamic level. The phase equations contain KPZ-like nonlinearities, density- or flow-dependent frequencies, and active currents, none of which can be derived from a free-energy functional satisfying fluctuation–dissipation. In the substrate model, the sign of the renormalized logarithmic correction controls both orientational order and density statistics: stronger-than-equilibrium QLRO is paired with miniscule or hyperuniform number fluctuations, while weaker-than-equilibrium QLRO is paired with giant number fluctuations [2104.10982]. In the flow-coupled model, after eliminating the Stokes field, the phase sector reduces to a KPZ-like theory with additional marginal nonlinearities. Depending on the dimensionless active couplings \(\alpha_1\) or \(\alpha_2\), the effective nonlinearity can flow to zero, yielding QLRO, stronger-than-QLRO, or weaker-than-QLRO, or to strong coupling, yielding disorder [1911.00214].

Driven-dissipative XY-type systems can also be organized around fluxes rather than order parameters. A modified XY model with amplitude fluctuations and projected nearest-neighbor coupling supports a nonequilibrium steady state with inverse energy cascade, constant negative flux \(\Pi(K)\simeq -\epsilon\), and energy spectrum
\[
E(k)\propto k^{-3},
\]
showing that “active XY” behavior can include cascade transfer and non-Kolmogorov universality rather than only defect-mediated ordering transitions [2106.11670].

## 6. Methods, benchmarks, and scope of the term

The numerical and analytical toolkit is correspondingly diverse. Persistent-noise XY models have been simulated on triangular lattices of about \(230{,}000\) sites using explicit Euler integration with \(\Delta t=0.01\) for \(10^7\) steps, measuring the polar order parameter, \(g_\theta(r)\), and \(\langle|\tilde{\Theta}(\mathbf{k})|^2\rangle\) [2401.11175]. The later BKT analysis on the same class of models used system sizes up to \(L=3360\), very long runs up to \(10^7\)–\(5\times10^8\) time steps, susceptibility finite-size scaling, and quenches from perfectly ordered initial conditions to extract \(T_c\), \(\eta_c\), \(\lambda_c\), and \(z\) [2602.22028].

Self-driven rotor studies combine steady-state simulations with explicit vortex tracking on square lattices up to \(L=500\); free vortices are identified by winding number and paired using the Hungarian algorithm [2103.12578]. The active lattice-gas XY model uses asynchronous Monte Carlo updates, Hoshen–Kopelman cluster analysis, and time-averaged vorticity heatmaps over the \((\epsilon,\rho)\) plane [2507.17108]. Hydrodynamic active-spin theories rely on one-loop dynamic renormalization group, in which the marginal couplings remain visible at two dimensions and generate logarithmic rather than simple power-law corrections [2104.10982][1911.00214].

Equilibrium tensor-network calculations remain important because they define the passive baseline for defect content and transition structure. Higher-order tensor renormalization group and HOTRG calculations of generalized XY models resolve the nematic, ferromagnetic, and disordered phases, track integer and half-integer susceptibilities, and locate the BKT and Ising lines in the \(q=2\) model [2307.13593][2404.17504]. Programmable photonic simulators implement classical XY Hamiltonians directly through optical vector–matrix multiplication and provide a hardware benchmark for BKT physics and ground-state search in finite XY systems [2401.08055].

Taken together, these studies establish that the active XY model is best understood as a class of nonequilibrium \(U(1)\) theories rather than a single model. The decisive distinction from the equilibrium XY model is not merely the presence of noise, but the way nonequilibrium drive reorganizes the relation among spin waves, defects, transport, and density. Depending on the mechanism—persistent torque, intrinsic-frequency disorder, self-propulsion, coupling to conserved species, or coupling to flow—activity can extend QLRO beyond equilibrium bounds, destroy it immediately, generate rotating non-ergodic states, or reorganize defects into phase-separating structures [2401.11175][2103.12578][2507.17108][2104.10982][1911.00214][1107.0370].

Source: https://www.emergentmind.com/topics/active-xy-model