Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2606.25340

Kosterlitz--Thouless Criticality in a Dipole-Conserving XY Model

Han-Xie Wang, Shuai A. Chen, Zheng Yan, Peng Ye

Published Jun 24, 2026• 0 citations
arXivHTMLPDF

Abstract

In this Letter, we study finite-temperature phase transitions in a two-dimensional classical statistical model called ``dipole-conserving XY model''. Unlike the conventional XY model, the original phase field $θ$ has no quasi-long-range order, conventional phase vortices have finite self-energy, and the standard helicity modulus vanishes identically. Analytic vortex energetics and Gaussian continuum theory show that Kosterlitz--Thouless (KT) criticality is instead controlled by phase-gradient (PG) vortices defined in the two compact dipole fields $χ_α=a\partial_αθ$, whose self-energies grow logarithmically with system size. We determine the phase diagram using parallel-tempering Metropolis Monte Carlo simulations and three diagnostics tailored to dipole conservation. KT finite-size scaling of dipole-field correlation-ratio crossings locates the critical temperatures; generalized helicity moduli defined through quadratic phase twists measure the stiffness of individual dipole channels; and PG-vortex densities obtained from plaquette winding numbers identify the proliferating vortex species. At isotropic couplings, the two PG-vortex species unbind simultaneously at a single KT transition. Spatial anisotropy separates their unbinding temperatures, yielding two KT transitions and an intermediate phase with quasi-long-range order in only one dipole channel. The transitions merge again when the mixed-derivative channel is removed. Our results establish the finite-temperature phase structure of the dipole-conserving XY model and identify thermal PG-vortex unbinding as a novel route to KT criticality in systems with higher-moment conservation.

DOI
10.48550/arxiv.2606.25340

Related papers

Phase transitions in the two-dimensional XY model with random phases: A Monte Carlo study

J. Maucourt, Daniel R. Grempel

We study the two-dimensional $\mathrm{XY}$ model with quenched random phases by Monte Carlo simulation and finite-size scaling analysis. We determine the phase diagram of the model and study its critical behavior as a function of disorder and temperature. If the strength of the randomness is less than a critical value, ${\ensuremath{\sigma}}_{c},$ the system has a Kosterlitz-Thouless (KT) phase transition from the paramagnetic phase to a state with quasi-long-range order. Our data suggest that the latter exists down to $T=0$ in contradiction with theories that predict the appearance of a low-temperature reentrant phase. At the critical disorder ${T}_{\mathrm{KT}}\ensuremath{\rightarrow}0$ and for $\ensuremath{\sigma}>{\ensuremath{\sigma}}_{c}$ there is no quasiordered phase. At zero temperature there is a phase transition between two different glassy states at ${\ensuremath{\sigma}}_{c}.$ The functional dependence of the correlation length on $\ensuremath{\sigma}$ suggests that this transition corresponds to the disorder-driven unbinding of vortex pairs.

arXiv:cond-mat/9703109•Physical review. B, Condensed matter•Aug 1, 1997•30 citations

Berezinskii-Kosterlitz-Thouless transition in the XY model on the honeycomb lattice: a comprehensive Monte Carlo analysis

Felipe E. F. de Andrade, L. N. Jorge, Cláudio José DaSilva

Abstract In this paper, we thoroughly examined the Berezinskii-Kosterlitz-Thouless (BKT) phase transition in the two-dimensional XY model on the honeycomb lattice. To address its thermodynamical behavior, we combined standard numerical Monte Carlo simulations with the simulated annealing (SA) protocol and entropic simulations based on the Wang-Landau (WL) algorithm. The transition temperature was determined using the second (ϒ) and fourth-order (ϒ4) helicity modulus as the order parameter. Our best finite-size scaling analysis suggests T BKT = 0.575(8) from SA and T BKT = 0.576(3) from WL. These values deviate significantly from the expected theoretical value of 1 / 2 . We believe that this discrepancy suggests that the theoretical assumptions regarding the analytical calculation may need to be revisited. Additionally, we calculated the vortex density and the formation energy of the vortex-antivortex pairs, where the obtained vortex formation energy is 2μ = 5.80(12). Upon comparison with the square lattice, our results support the notion of instability of the honeycomb lattice to support the spin long-range order and give additional backing to the critical behavior we found.

arXiv:2406.12076•Physica Scripta•May 8, 2025•3 citations

Numerical Studies of Vortices and Helicity Modulus in the Two-Dimensional Generalized XY Model

Yunzhou Sun, Qin Wu, Yang Xiao-li +4 more

Two-dimensional generalized XY spin model on a triangular lattice is studied by means of Monte-Carlo simulations. The critical temperatures of Berezinskii-Kosterlitz-Thouless (BKT) phase transition are obtained by the method of helicity modulus. It is found that the results are consistent with those obtained by other methods. The vortex density and the vortex-antivortex pair formation energy are also obtained. The result shows that the critical temperature decreases with the increase of the generalization parameter q . While the vortex-antivortex pair formation energy increases with the increase of q when q>1 .

Frontiers in Physics•Mar 4, 2022•3 citations

Critical behavior of the two-dimensionalXYmodel

Rajan Gupta, CLIVE F. BAILLIE

We present detailed Monte Carlo results for the susceptibility \ensuremath{\chi}, correlation length \ensuremath{\xi}, and specific heat ${\mathit{C}}_{\mathit{v}}$ for the XY model. The simulations are done on ${64}^{2}$, ${128}^{2}$, ${256}^{2}$, and ${512}^{2}$ lattices over the temperature range 0.98\ensuremath{\le}T\ensuremath{\le}1.43 corresponding to 270. Fits to \ensuremath{\chi} and \ensuremath{\xi} data favor a Kosterlitz-Thouless (KT) singularity over a second-order transition; however, unconstrained four-parameter KT fits do not confirm the predicted values \ensuremath{\nu}=0.5 and \ensuremath{\eta}=0.25. Our best estimate, ${\mathit{T}}_{\mathit{c}}$=0.894(5), is obtained using KT fits with \ensuremath{\nu} fixed at 0.5. The exponent \ensuremath{\eta} is calculated as a function of temperature in the spin-wave phase using Monte Carlo renormalization-group and finite-size-scaling methods. Both methods give consistent results and we find \ensuremath{\eta}\ensuremath{\approxeq}0.235 at T=0.894. We also present results for the behavior of vortex density across the transition and exhibit how the dilute-gas approximation breaks down.

Physical review. B, Condensed matter•Feb 1, 1992•196 citations

The two-dimensional XY model at the transition temperature: a high-precision Monte Carlo study

Martin Hasenbusch

We study the classical XY (plane rotator) model at the Kosterlitz–Thouless phase transition. We simulate the model using the single-cluster algorithm on square lattices of a linear size up to L = 2048. We derive the finite-size behaviour of the second moment correlation length over the lattice size ξ 2nd / L at the transition temperature. This new prediction and the analogous one for the helicity modulus ϒ are confronted with our Monte Carlo data. This way β KT = 1.1199 is confirmed as inverse transition temperature. Finally, we address the puzzle of logarithmic corrections of the magnetic susceptibility χ at the transition temperature.

arXiv:cond-mat/0502556•Journal of Physics A Mathematical and General•Jun 15, 2005•165 citations

Universality of the Berezinskii–Kosterlitz–Thouless type of phase transition in the dipolar XY-model

A Yu Vasiliev, Andrei E. Tarkhov, L I Menshikov +2 more

We investigate the nature of the phase transition occurring in a planar XY-model spin system with dipole–dipole interactions. It is demonstrated that a Berezinskii–Kosterlitz–Thouless (BKT) type of phase transition always takes place at a finite temperature separating the ordered (ferro) and the disordered (para) phases. The low-temperature phase corresponds to an ordered state with thermal fluctuations, composed of a 'gas' of bound vortex–antivortex pairs, which would, when considered isolated, be characterized by a constant vortex–antivortex attraction force which is due to the dipolar interaction term in the Hamiltonian. Using a topological charge model, we show that small bound pairs are easily polarized, and screen the vortex–antivortex interaction in sufficiently large pairs. Screening changes the linear attraction potential of vortices to a logarithmic one, and leads to the familiar pair dissociation mechanism of the BKT type phase transition. The topological charge model is confirmed by numerical simulations, in which we demonstrate that the transition temperature slightly increases when compared with the BKT result for short-range interactions.

arXiv:1303.4915•New Journal of Physics•May 7, 2014•11 citations
In this Letter, we study finite-temperature phase transitions in a two-dimensional classical statistical model called ``dipole-conserving XY model''. Unlike the conventional XY model, the original phase field θθθ has no quasi-long-range order, conventional phase vortices have finite self-energy, and the standard helicity modulus vanishes identically. Analytic vortex energetics and Gaussian continuum theory show that Kosterlitz--Thouless (KT) criticality is instead controlled by phase-gradient (PG) vortices defined in the two compact dipole fields χα=a∂αθχ_α=a\partial_αθχα​=a∂α​θ, whose self-energies grow logarithmically with system size. We determine the phase diagram using parallel-tempering Metropolis Monte Carlo simulations and three diagnostics tailored to dipole conservation. KT finite-size scaling of dipole-field correlation-ratio crossings locates the critical temperatures; generalized helicity moduli defined through quadratic phase twists measure the stiffness of individual dipole channels; and PG-vortex densities obtained from plaquette winding numbers identify the proliferating vortex species. At isotropic couplings, the two PG-vortex species unbind simultaneously at a single KT transition. Spatial anisotropy separates their unbinding temperatures, yielding two KT transitions and an intermediate phase with quasi-long-range order in only one dipole channel. The transitions merge again when the mixed-derivative channel is removed. Our results establish the finite-temperature phase structure of the dipole-conserving XY model and identify thermal PG-vortex unbinding as a novel route to KT criticality in systems with higher-moment conservation.
We study the two-dimensional XY\mathrm{XY}XY model with quenched random phases by Monte Carlo simulation and finite-size scaling analysis. We determine the phase diagram of the model and study its critical behavior as a function of disorder and temperature. If the strength of the randomness is less than a critical value, ${\ensuremath{\sigma}}_{c},$ the system has a Kosterlitz-Thouless (KT) phase transition from the paramagnetic phase to a state with quasi-long-range order. Our data suggest that the latter exists down to T=0T=0T=0 in contradiction with theories that predict the appearance of a low-temperature reentrant phase. At the critical disorder ${T}_{\mathrm{KT}}\ensuremath{\rightarrow}0$ and for $\ensuremath{\sigma}>{\ensuremath{\sigma}}_{c}$ there is no quasiordered phase. At zero temperature there is a phase transition between two different glassy states at ${\ensuremath{\sigma}}_{c}.$ The functional dependence of the correlation length on $\ensuremath{\sigma}$ suggests that this transition corresponds to the disorder-driven unbinding of vortex pairs.
Abstract In this paper, we thoroughly examined the Berezinskii-Kosterlitz-Thouless (BKT) phase transition in the two-dimensional XY model on the honeycomb lattice. To address its thermodynamical behavior, we combined standard numerical Monte Carlo simulations with the simulated annealing (SA) protocol and entropic simulations based on the Wang-Landau (WL) algorithm. The transition temperature was determined using the second (ϒ) and fourth-order (ϒ4) helicity modulus as the order parameter. Our best finite-size scaling analysis suggests T BKT = 0.575(8) from SA and T BKT = 0.576(3) from WL. These values deviate significantly from the expected theoretical value of 1 / 2 . We believe that this discrepancy suggests that the theoretical assumptions regarding the analytical calculation may need to be revisited. Additionally, we calculated the vortex density and the formation energy of the vortex-antivortex pairs, where the obtained vortex formation energy is 2μ = 5.80(12). Upon comparison with the square lattice, our results support the notion of instability of the honeycomb lattice to support the spin long-range order and give additional backing to the critical behavior we found.
Two-dimensional generalized XY spin model on a triangular lattice is studied by means of Monte-Carlo simulations. The critical temperatures of Berezinskii-Kosterlitz-Thouless (BKT) phase transition are obtained by the method of helicity modulus. It is found that the results are consistent with those obtained by other methods. The vortex density and the vortex-antivortex pair formation energy are also obtained. The result shows that the critical temperature decreases with the increase of the generalization parameter q . While the vortex-antivortex pair formation energy increases with the increase of q when q>1 .
We present detailed Monte Carlo results for the susceptibility \ensuremath{\chi}, correlation length \ensuremath{\xi}, and specific heat Cv{\mathit{C}}_{\mathit{v}}Cv​ for the XY model. The simulations are done on 642{64}^{2}642, 1282{128}^{2}1282, 2562{256}^{2}2562, and 5122{512}^{2}5122 lattices over the temperature range 0.98\ensuremath{\le}T\ensuremath{\le}1.43 corresponding to 270. Fits to \ensuremath{\chi} and \ensuremath{\xi} data favor a Kosterlitz-Thouless (KT) singularity over a second-order transition; however, unconstrained four-parameter KT fits do not confirm the predicted values \ensuremath{\nu}=0.5 and \ensuremath{\eta}=0.25. Our best estimate, Tc{\mathit{T}}_{\mathit{c}}Tc​=0.894(5), is obtained using KT fits with \ensuremath{\nu} fixed at 0.5. The exponent \ensuremath{\eta} is calculated as a function of temperature in the spin-wave phase using Monte Carlo renormalization-group and finite-size-scaling methods. Both methods give consistent results and we find \ensuremath{\eta}\ensuremath{\approxeq}0.235 at T=0.894. We also present results for the behavior of vortex density across the transition and exhibit how the dilute-gas approximation breaks down.
We study the classical XY (plane rotator) model at the Kosterlitz–Thouless phase transition. We simulate the model using the single-cluster algorithm on square lattices of a linear size up to L = 2048. We derive the finite-size behaviour of the second moment correlation length over the lattice size ξ 2nd / L at the transition temperature. This new prediction and the analogous one for the helicity modulus ϒ are confronted with our Monte Carlo data. This way β KT = 1.1199 is confirmed as inverse transition temperature. Finally, we address the puzzle of logarithmic corrections of the magnetic susceptibility χ at the transition temperature.
We investigate the nature of the phase transition occurring in a planar XY-model spin system with dipole–dipole interactions. It is demonstrated that a Berezinskii–Kosterlitz–Thouless (BKT) type of phase transition always takes place at a finite temperature separating the ordered (ferro) and the disordered (para) phases. The low-temperature phase corresponds to an ordered state with thermal fluctuations, composed of a 'gas' of bound vortex–antivortex pairs, which would, when considered isolated, be characterized by a constant vortex–antivortex attraction force which is due to the dipolar interaction term in the Hamiltonian. Using a topological charge model, we show that small bound pairs are easily polarized, and screen the vortex–antivortex interaction in sufficiently large pairs. Screening changes the linear attraction potential of vortices to a logarithmic one, and leads to the familiar pair dissociation mechanism of the BKT type phase transition. The topological charge model is confirmed by numerical simulations, in which we demonstrate that the transition temperature slightly increases when compared with the BKT result for short-range interactions.