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:2007.01526 · Physical review. B./Physical review. B

Interplay of spin waves and vortices in the two-dimensional XY model at small vortex-core energy

Ilaria Maccari, Nicolò Defenu, Lara Benfatto, Claudio Castellani, Tilman Enss

Published Sep 9, 2020• 23 citations
arXivHTMLPDF

Abstract

The Berezinskii-Kosterlitz-Thouless (BKT) mechanism describes universal vortex unbinding in many two-dimensional systems, including the paradigmatic XY model. However, most of these systems present a complex interplay between excitations at different length scales that complicates theoretical calculations of nonuniversal thermodynamic quantities. These difficulties may be overcome by suitably modifying the initial conditions of the BKT flow equations to account for noncritical fluctuations at small length scales. In this work, we perform a systematic study of the validity and limits of this two-step approach by constructing optimised initial conditions for the BKT flow. We find that the two-step approach can accurately reproduce the results of Monte Carlo simulations of the traditional XY model. To systematically study the interplay between vortices and spin-wave excitations, we introduce a modified XY model with increased vortex fugacity. We present large-scale Monte Carlo simulations of the spin stiffness and vortex density for this modified XY model and show that even at large vortex fugacity, vortex unbinding is accurately described by the nonperturbative functional renormalization group.

DOI
10.1103/physrevb.102.104505

Related papers

Topological Lattice Actions for the 2d XY Model

Wolfgang Bietenholz A, Bögli, Michael, Niedermayer, Ferenc +3 more

We consider the 2d XY Model with topological lattice actions, which are invariant against small deformations of the field configuration. These actions constrain the angle between neighbouring spins by an upper bound, or they explicitly suppress vortices (and anti-vortices). Although topological actions do not have a classical limit, they still lead to the universal behaviour of the Berezinskii-Kosterlitz-Thouless (BKT) phase transition — at least up to moderate vortex suppression. In the massive phase, the analytically known Step Scaling Function (SSF) is reproduced in numerical simulations. However, deviations from the expected universal behaviour of the lattice artifacts are observed. In the massless phase, the BKT value of the critical exponent ηc is confirmed. Hence, even though for some topological actions vortices cost zero energy, they still drive the standard BKT transition. In addition we identify a vortex-free transition point, which deviates from the BKT behaviour.

•Jan 1, 2013•6 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

Quench dynamics of the 2dXYmodel

Asja Jelić, Leticia F. Cugliandolo

We investigate the out of equilibrium dynamics of the two-dimensional XY model when cooled across the Berezinskii–Kosterlitz–Thouless (BKT) phase transition using different protocols. We focus on the evolution of the growing correlation length and the density of topological defects (vortices). By using Monte Carlo simulations we first determine the time and temperature dependence of the growing correlation length after an infinitely rapid quench from above the transition temperature to the quasi-long-range order region. The functional form is consistent with a logarithmic correction to the diffusive law and it serves to validate dynamic scaling in this problem. This analysis clarifies the different dynamic roles played by bound and free vortices. We then revisit the Kibble–Zurek mechanism in thermal phase transitions in which the disordered state is plagued with topological defects. We provide a theory of quenching rate dependence in systems with the BKT-type transition that goes beyond the equilibrium scaling arguments. Finally, we discuss the implications of our results for a host of physical systems with vortex excitations, including planar ferromagnets and liquid crystals.

arXiv:1012.0417•Journal of Statistical Mechanics Theory and Experiment•Feb 18, 2011•89 citations

Scaling of the magnetic permeability at the Berezinskii–Kosterlitz–Thouless transition from Coulomb gas simulations

Rogelio Díaz-Méndez, Jack Lidmar, Mats Wallin

A new approach to the Berezinskii–Kosterlitz–Thouless transition in the two-dimensional Coulomb gas model is explored by Monte Carlo simulation and finite size scaling. The usual mapping of a neutral two-dimensional superconductor in zero magnetic field to a Coulomb gas leads to an unscreened logarithmic interaction between the vortices, and with periodic boundary conditions vortex configurations are always vorticity neutral with an equal number of plus and minus vortices. We demonstrate that relaxing the neutrality condition has certain advantages. It leads to non-neutral vortex configurations that can appear in real systems with open boundary conditions and permits calculation of the compressibility, which for thin film superconductors corresponds to the magnetic permeability. The vortex-number fluctuation has remarkable scaling properties at and below the Berezinskii–Kosterlitz–Thouless transition. The fugacity variable becomes dangerously irrelevant in the low-temperature phase and leads to a multiplicative scaling correction to the mean-square vortex-number fluctuation and to the magnetic permeability. This multiplicative correction strongly affects the scaling properties of the vorticity fluctuation at and below the transition. Consequences of these findings are demonstrated using Monte Carlo simulations. Inclusion of the next-higher order correction to scaling is found to play an important role in the analysis of numerical data for the vortex number fluctuation and permits accurate determination of the critical properties.

arXiv:1808.07183•Journal of Statistical Mechanics Theory and Experiment•Dec 3, 2018•1 citation

Finite-size scaling method for the Berezinskii–Kosterlitz–Thouless transition

Yun-Da Hsieh, Ying-Jer Kao, Anders W. Sandvik

We test an improved finite-size scaling method for reliably extracting the critical temperature T BKT of a Berezinskii–Kosterlitz–Thouless (BKT) transition. Using known single-parameter logarithmic corrections to the spin stiffness ρ s at T BKT in combination with the Kosterlitz–Nelson relation between the transition temperature and the stiffness, ρ s ( T BKT ) = 2 T BKT /π, we define a size-dependent transition temperature T BKT ( L 1 , L 2 ) based on a pair of system sizes L 1 , L 2 , e.g., L 2 = 2 L 1 . We use Monte Carlo data for the standard two-dimensional classical XY model to demonstrate that this quantity is well behaved and can be reliably extrapolated to the thermodynamic limit using the next expected logarithmic correction beyond the ones included in defining T BKT ( L 1 , L 2 ). For the Monte Carlo calculations we use GPU (graphical processing unit) computing to obtain high-precision data for L up to 512. We find that the sub-leading logarithmic corrections have significant effects on the extrapolation. Our result T BKT = 0.8935(1) is several error bars above the previously best estimates of the transition temperature, T BKT ≈ 0.8929. If only the leading log-correction is used, the result is, however, consistent with the lower value, suggesting that previous works have underestimated T BKT because of the neglect of sub-leading logarithms. Our method is easy to implement in practice and should be applicable to generic BKT transitions.

arXiv:1302.2900•Journal of Statistical Mechanics Theory and Experiment•Sep 6, 2013•81 citations

Berezinskii-Kosterlitz-Thouless phase transitions in two-dimensional non-Abelian spin models

Oleg Borisenko, V. Chelnokov, Francesca Cuteri +1 more

It is argued that two-dimensional U(N) spin models for any N undergo a Berezinskii-Kosterlitz-Thouless (BKT)-like phase transition, similarly to the famous XY model. This conclusion follows from the Berezinskii-like calculation of the two-point correlation function in U(N) models, approximate renormalization group analysis, and numerical investigations of the U(2) model. It is shown, via Monte Carlo simulations, that the universality class of the U(2) model coincides with that of the XY model. Moreover, preliminary numerical results point out that two-dimensional SU(N) spin models with the fundamental and adjoint terms and N>4 exhibit two phase transitions of BKT type, similarly to Z(N) vector models.

arXiv:1512.05737•Physical review. E•Jul 7, 2016•2 citations
The Berezinskii-Kosterlitz-Thouless (BKT) mechanism describes universal vortex unbinding in many two-dimensional systems, including the paradigmatic XY model. However, most of these systems present a complex interplay between excitations at different length scales that complicates theoretical calculations of nonuniversal thermodynamic quantities. These difficulties may be overcome by suitably modifying the initial conditions of the BKT flow equations to account for noncritical fluctuations at small length scales. In this work, we perform a systematic study of the validity and limits of this two-step approach by constructing optimised initial conditions for the BKT flow. We find that the two-step approach can accurately reproduce the results of Monte Carlo simulations of the traditional XY model. To systematically study the interplay between vortices and spin-wave excitations, we introduce a modified XY model with increased vortex fugacity. We present large-scale Monte Carlo simulations of the spin stiffness and vortex density for this modified XY model and show that even at large vortex fugacity, vortex unbinding is accurately described by the nonperturbative functional renormalization group.
We consider the 2d XY Model with topological lattice actions, which are invariant against small deformations of the field configuration. These actions constrain the angle between neighbouring spins by an upper bound, or they explicitly suppress vortices (and anti-vortices). Although topological actions do not have a classical limit, they still lead to the universal behaviour of the Berezinskii-Kosterlitz-Thouless (BKT) phase transition — at least up to moderate vortex suppression. In the massive phase, the analytically known Step Scaling Function (SSF) is reproduced in numerical simulations. However, deviations from the expected universal behaviour of the lattice artifacts are observed. In the massless phase, the BKT value of the critical exponent ηc is confirmed. Hence, even though for some topological actions vortices cost zero energy, they still drive the standard BKT transition. In addition we identify a vortex-free transition point, which deviates from the BKT behaviour.
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.
We investigate the out of equilibrium dynamics of the two-dimensional XY model when cooled across the Berezinskii–Kosterlitz–Thouless (BKT) phase transition using different protocols. We focus on the evolution of the growing correlation length and the density of topological defects (vortices). By using Monte Carlo simulations we first determine the time and temperature dependence of the growing correlation length after an infinitely rapid quench from above the transition temperature to the quasi-long-range order region. The functional form is consistent with a logarithmic correction to the diffusive law and it serves to validate dynamic scaling in this problem. This analysis clarifies the different dynamic roles played by bound and free vortices. We then revisit the Kibble–Zurek mechanism in thermal phase transitions in which the disordered state is plagued with topological defects. We provide a theory of quenching rate dependence in systems with the BKT-type transition that goes beyond the equilibrium scaling arguments. Finally, we discuss the implications of our results for a host of physical systems with vortex excitations, including planar ferromagnets and liquid crystals.
A new approach to the Berezinskii–Kosterlitz–Thouless transition in the two-dimensional Coulomb gas model is explored by Monte Carlo simulation and finite size scaling. The usual mapping of a neutral two-dimensional superconductor in zero magnetic field to a Coulomb gas leads to an unscreened logarithmic interaction between the vortices, and with periodic boundary conditions vortex configurations are always vorticity neutral with an equal number of plus and minus vortices. We demonstrate that relaxing the neutrality condition has certain advantages. It leads to non-neutral vortex configurations that can appear in real systems with open boundary conditions and permits calculation of the compressibility, which for thin film superconductors corresponds to the magnetic permeability. The vortex-number fluctuation has remarkable scaling properties at and below the Berezinskii–Kosterlitz–Thouless transition. The fugacity variable becomes dangerously irrelevant in the low-temperature phase and leads to a multiplicative scaling correction to the mean-square vortex-number fluctuation and to the magnetic permeability. This multiplicative correction strongly affects the scaling properties of the vorticity fluctuation at and below the transition. Consequences of these findings are demonstrated using Monte Carlo simulations. Inclusion of the next-higher order correction to scaling is found to play an important role in the analysis of numerical data for the vortex number fluctuation and permits accurate determination of the critical properties.
We test an improved finite-size scaling method for reliably extracting the critical temperature T BKT of a Berezinskii–Kosterlitz–Thouless (BKT) transition. Using known single-parameter logarithmic corrections to the spin stiffness ρ s at T BKT in combination with the Kosterlitz–Nelson relation between the transition temperature and the stiffness, ρ s ( T BKT ) = 2 T BKT /π, we define a size-dependent transition temperature T BKT ( L 1 , L 2 ) based on a pair of system sizes L 1 , L 2 , e.g., L 2 = 2 L 1 . We use Monte Carlo data for the standard two-dimensional classical XY model to demonstrate that this quantity is well behaved and can be reliably extrapolated to the thermodynamic limit using the next expected logarithmic correction beyond the ones included in defining T BKT ( L 1 , L 2 ). For the Monte Carlo calculations we use GPU (graphical processing unit) computing to obtain high-precision data for L up to 512. We find that the sub-leading logarithmic corrections have significant effects on the extrapolation. Our result T BKT = 0.8935(1) is several error bars above the previously best estimates of the transition temperature, T BKT ≈ 0.8929. If only the leading log-correction is used, the result is, however, consistent with the lower value, suggesting that previous works have underestimated T BKT because of the neglect of sub-leading logarithms. Our method is easy to implement in practice and should be applicable to generic BKT transitions.
It is argued that two-dimensional U(N) spin models for any N undergo a Berezinskii-Kosterlitz-Thouless (BKT)-like phase transition, similarly to the famous XY model. This conclusion follows from the Berezinskii-like calculation of the two-point correlation function in U(N) models, approximate renormalization group analysis, and numerical investigations of the U(2) model. It is shown, via Monte Carlo simulations, that the universality class of the U(2) model coincides with that of the XY model. Moreover, preliminary numerical results point out that two-dimensional SU(N) spin models with the fundamental and adjoint terms and N>4 exhibit two phase transitions of BKT type, similarly to Z(N) vector models.