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:2310.20314 · Chinese Physics B

Emergent topological ordered phase for the Ising-XY model revealed by cluster-updating Monte Carlo method

Heyang Ma, Wanzhou 万舟 Zhang 张, Yanting 彦婷 Tian 田, Chengxiang 成祥 Ding 丁, Youjin Deng

Published Jan 11, 2024• 4 citations
arXivHTMLPDF

Abstract

The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy. At low temperatures, theoretical predictions [Phys. Rev. A 72 053604 (2005)] and [arXiv: 0706.1609] indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering. However, due to ergodic difficulties faced by Monte Carlo methods at low temperatures, this topological phase has not been numerically explored. We propose a linear cluster updating Monte Carlo method, which flips spins without rejection in the anisotropy limit but does not change the energy. Using this scheme and conventional Monte Carlo methods, we succeed in revealing the nature of topological phases with half-vortices and domain walls. In the constructed global phase diagram, Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier. We also propose and explore a wide range of quantities, including magnetism, superfluidity, specific heat, susceptibility, and even percolation susceptibility, and obtain consistent and reliable results. Furthermore, we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes, as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections. The critical exponents of different types of phase transitions are reasonably fitted. The results are useful to help cold atom experiments explore the half-vortex topological phase.

DOI
10.1088/1674-1056/ad1d4d

Related papers

Phase transitions in coupled two dimensional XY systems with spatial anisotropy

Cenke Xu

We study phase transitions of coupled two dimensional XY systems with spatial anisotropy and $U(1) \times \mathbb{Z}_2$ symmetry, motivated by spinless bosonic atoms trapped in square optical lattice on the metastable first excited $p-$level orbitals with anisotropic hopping amplitudes. The phase transitions of the system are generally split into an Ising transition and an XY transition, but the sequence and the nature of the transitions depend on the ratio between the anisotropic couplings. In the isotropic limit the XY variables are expected to be disordered before the Ising variables when thermal or quantum fluctuations are turned on gradually. In the anisotropic limit with zero perpendicular hoppings, the finite temperature transition is a Kosterlitz-Thouless transition driven by proliferation of hybrid half vortices, and the zero temperature quantum phase transition is split into a bond order transition and a 3D XY transition, which can be driven by the condensation of either single vortices or half vortices. After the condensation of half vortices the resultant state is a Mott Insulator of paired bosons. A small perpendicular hopping $J_b$ leads to a 2D Ising transition at low temperature and a 2+1d quantum Ising transition with a small charging energy at zero temperature. Global phase diagrams for both classical and quantum phase transitions are drawn. The analytical results obtained in this work are expected to be checked both numerically and experimentally.

arXiv:0706.1609•Jun 12, 2007•1 citation

Critical Behavior Simulation of the Two-Dimensional Ising Model—A Computational Physics Study Based on the Monte Carlo Method

佳钢 李

To demonstrate the application of computational physics in modern statistical physics, this paper investigates the critical behavior of a finite-size ferromagnetic system based on the two-dimensional square-lattice Ising model.A nearest-neighbor Ising Hamiltonian is constructed, and periodic boundary conditions are adopted to reduce finite-boundary effects.The Metropolis Monte Carlo algorithm is used to sample spin configurations at different temperatures, while energy, magnetization, specific heat, magnetic susceptibility and Binder cumulant are calculated to describe the transition from an ordered state to a disordered state.The simulation results show that the system exhibits clear long-range order at low temperatures, where the order parameter remains high.As temperature increases, thermal fluctuations become stronger, magnetization decreases rapidly, and the system gradually evolves into a disordered state.The specific heat curve shows a pronounced peak at approximately (T = 2.269), while the susceptibility peak shifts slightly toward the high-temperature side due to finite-size effects and the discrete temperature grid.Overall, the simulated critical region is consistent with the theoretical critical temperature of the two-dimensional square-lattice Ising model.Spin-configuration analysis further reveals enhanced fluctuations and domain reconstruction near the critical region.In addition, a thermalization criterion and an extended sampling test for L = 48 near the critical region are added to examine the influence of finite sampling and equilibration.These results indicate that Monte Carlo simulation of the two-dimensional Ising model can effectively reproduce continuous phase transition, critical fluctuation and finite-size effects, providing a basic computational case for understanding critical phenomena in complex systems and data-driven physical analysis.

Modern Physics•Jan 1, 2026•0 citations

Monte Carlo simulation with Wolff algorithm for scaling behavior of two dimensional XY model with KT phase transition

Yuntai Song

Abstract The property of XY model with its KT phase transition has been discussed in many aspects such as specific heat and magnetic susceptibility along with correlation function. Also, this is one of the original models for topological phase transition which gives a new type of phase transition without symmetry breaking. Since the calculation of these properties above can be achieved by Monte Carlo Method with Wolff algorithm, the visualization of vortex and antivortex is also an interesting task which is accomplished in this article. By successfully modeling the system and getting the step-by-step configuration, the relation between the generation of vortex and the variation of temperature can be concluded in a qualitative way. The calculation of spin correlation function shows the phase transition between disorder and quasi long-range order which represents the topological phase transition. As a reliable method, the Monte Carlo method with Wolff algorithm can be a choice of generating the training set for neural network which can calculated much more sophisticated case in a relatively short time and the transition between classical simulation and quantum simulation can be approached as mentioned in the end of this article.

Journal of Physics Conference Series•Nov 1, 2023•0 citations

Percolation of the two-dimensional XY model in the flow representation

Bao‐Zong Wang, Pengcheng Hou, Chun-Jiong Huang +1 more

We simulate the two-dimensional XY model in the flow representation by a worm-type algorithm, up to linear system size L=4096, and study the geometric properties of the flow configurations. As the coupling strength K increases, we observe that the system undergoes a percolation transition K_{perc} from a disordered phase consisting of small clusters into an ordered phase containing a giant percolating cluster. Namely, in the low-temperature phase, there exhibits a long-ranged order regarding the flow connectivity, in contrast to the quasi-long-range order associated with spin properties. Near K_{perc}, the scaling behavior of geometric observables is well described by the standard finite-size scaling ansatz for a second-order phase transition. The estimated percolation threshold K_{perc}=1.1053(4) is close to but obviously smaller than the Berezinskii-Kosterlitz-Thouless (BKT) transition point K_{BKT}=1.1193(10), which is determined from the magnetic susceptibility and the superfluid density. Various interesting questions arise from these unconventional observations, and their solutions would shed light on a variety of classical and quantum systems of BKT phase transitions.

arXiv:2010.14427•Physical review. E•Jun 21, 2021•8 citations

Extinction of BKT transition by nonmagnetic disorder in planar-symmetry spin models

Gary Mathew Wysin, Alejandro Pereira, Ivo A. Marques +2 more

The Berezinskii-Kosterlitz-Thouless (BKT) transition in two-dimensional planar rotator and XY models on a square lattice, diluted by randomly placed vacancies, is studied here using hybrid Monte Carlo simulations that combine single spin flip, cluster and over-relaxation techniques. The transition temperature $T_c$ is determined as a function of vacancy density $ρ_{vac}$ by calculations of the helicity modulus and the by finite-size scaling of the in-plane magnetic susceptibility. The results for $T_c$ are consistent with those from the much less precise fourth-order cumulant of Binder. $T_c$ is found to decrease monotonically with increasing $ρ_{vac}$, and falls to zero close to the square lattice percolation limit, $ρ_{vac}\approx 0.41$ . The result is physically reasonable: the long-range orientational order of the low-temperature phase cannot be maintained in the absence of sufficient spin interactions across the lattice.

arXiv:cond-mat/0504145•Apr 6, 2005•1 citation

Monte Carlo Investigation of Magnetic Phase Transitions in the Two-Dimensional Ising Model

Satyendra Lohar, Neha Srivas, Akhilesh Dwivedi

Many problems in statistical physics resist exact solution once the number of interacting particles becomes large, and the two-dimensional Ising model is a classic case where numerical simulation fills that gap. This paper reports a Monte Carlo study of the two-dimensional Ising model using the Metropolis algorithm, undertaken to examine how magnetization, internal energy, specific heat, and magnetic susceptibility change with temperature. Spin configurations were generated on a square lattice under periodic boundary conditions, with each proposed spin flip accepted or rejected according to the Boltzmann weight of the resulting energy change. The simulation reproduces the expected transition from an ordered ferromagnetic state at low temperature to a disordered paramagnetic state at high temperature, with magnetization falling sharply and both specific heat and susceptibility peaking near a critical temperature of Tc ≈ 2.3–2.4 (in units of J/kB), consistent with the exact value obtained analytically for this system. These results confirm that the Metropolis algorithm, despite its simplicity, remains a dependable tool for locating second-order phase transitions and extracting thermodynamic quantities in lattice spin systems, and the same computational framework can be adapted to more elaborate models used in condensed matter and materials research.

•Aug 14, 2026•0 citations
The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy. At low temperatures, theoretical predictions [Phys. Rev. A 72 053604 (2005)] and [arXiv: 0706.1609] indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering. However, due to ergodic difficulties faced by Monte Carlo methods at low temperatures, this topological phase has not been numerically explored. We propose a linear cluster updating Monte Carlo method, which flips spins without rejection in the anisotropy limit but does not change the energy. Using this scheme and conventional Monte Carlo methods, we succeed in revealing the nature of topological phases with half-vortices and domain walls. In the constructed global phase diagram, Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier. We also propose and explore a wide range of quantities, including magnetism, superfluidity, specific heat, susceptibility, and even percolation susceptibility, and obtain consistent and reliable results. Furthermore, we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes, as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections. The critical exponents of different types of phase transitions are reasonably fitted. The results are useful to help cold atom experiments explore the half-vortex topological phase.
We study phase transitions of coupled two dimensional XY systems with spatial anisotropy and U(1)×Z2U(1) \times \mathbb{Z}_2U(1)×Z2​ symmetry, motivated by spinless bosonic atoms trapped in square optical lattice on the metastable first excited p−p-p−level orbitals with anisotropic hopping amplitudes. The phase transitions of the system are generally split into an Ising transition and an XY transition, but the sequence and the nature of the transitions depend on the ratio between the anisotropic couplings. In the isotropic limit the XY variables are expected to be disordered before the Ising variables when thermal or quantum fluctuations are turned on gradually. In the anisotropic limit with zero perpendicular hoppings, the finite temperature transition is a Kosterlitz-Thouless transition driven by proliferation of hybrid half vortices, and the zero temperature quantum phase transition is split into a bond order transition and a 3D XY transition, which can be driven by the condensation of either single vortices or half vortices. After the condensation of half vortices the resultant state is a Mott Insulator of paired bosons. A small perpendicular hopping JbJ_bJb​ leads to a 2D Ising transition at low temperature and a 2+1d quantum Ising transition with a small charging energy at zero temperature. Global phase diagrams for both classical and quantum phase transitions are drawn. The analytical results obtained in this work are expected to be checked both numerically and experimentally.
To demonstrate the application of computational physics in modern statistical physics, this paper investigates the critical behavior of a finite-size ferromagnetic system based on the two-dimensional square-lattice Ising model.A nearest-neighbor Ising Hamiltonian is constructed, and periodic boundary conditions are adopted to reduce finite-boundary effects.The Metropolis Monte Carlo algorithm is used to sample spin configurations at different temperatures, while energy, magnetization, specific heat, magnetic susceptibility and Binder cumulant are calculated to describe the transition from an ordered state to a disordered state.The simulation results show that the system exhibits clear long-range order at low temperatures, where the order parameter remains high.As temperature increases, thermal fluctuations become stronger, magnetization decreases rapidly, and the system gradually evolves into a disordered state.The specific heat curve shows a pronounced peak at approximately (T = 2.269), while the susceptibility peak shifts slightly toward the high-temperature side due to finite-size effects and the discrete temperature grid.Overall, the simulated critical region is consistent with the theoretical critical temperature of the two-dimensional square-lattice Ising model.Spin-configuration analysis further reveals enhanced fluctuations and domain reconstruction near the critical region.In addition, a thermalization criterion and an extended sampling test for L = 48 near the critical region are added to examine the influence of finite sampling and equilibration.These results indicate that Monte Carlo simulation of the two-dimensional Ising model can effectively reproduce continuous phase transition, critical fluctuation and finite-size effects, providing a basic computational case for understanding critical phenomena in complex systems and data-driven physical analysis.
Abstract The property of XY model with its KT phase transition has been discussed in many aspects such as specific heat and magnetic susceptibility along with correlation function. Also, this is one of the original models for topological phase transition which gives a new type of phase transition without symmetry breaking. Since the calculation of these properties above can be achieved by Monte Carlo Method with Wolff algorithm, the visualization of vortex and antivortex is also an interesting task which is accomplished in this article. By successfully modeling the system and getting the step-by-step configuration, the relation between the generation of vortex and the variation of temperature can be concluded in a qualitative way. The calculation of spin correlation function shows the phase transition between disorder and quasi long-range order which represents the topological phase transition. As a reliable method, the Monte Carlo method with Wolff algorithm can be a choice of generating the training set for neural network which can calculated much more sophisticated case in a relatively short time and the transition between classical simulation and quantum simulation can be approached as mentioned in the end of this article.
We simulate the two-dimensional XY model in the flow representation by a worm-type algorithm, up to linear system size L=4096, and study the geometric properties of the flow configurations. As the coupling strength K increases, we observe that the system undergoes a percolation transition K_{perc} from a disordered phase consisting of small clusters into an ordered phase containing a giant percolating cluster. Namely, in the low-temperature phase, there exhibits a long-ranged order regarding the flow connectivity, in contrast to the quasi-long-range order associated with spin properties. Near K_{perc}, the scaling behavior of geometric observables is well described by the standard finite-size scaling ansatz for a second-order phase transition. The estimated percolation threshold K_{perc}=1.1053(4) is close to but obviously smaller than the Berezinskii-Kosterlitz-Thouless (BKT) transition point K_{BKT}=1.1193(10), which is determined from the magnetic susceptibility and the superfluid density. Various interesting questions arise from these unconventional observations, and their solutions would shed light on a variety of classical and quantum systems of BKT phase transitions.
The Berezinskii-Kosterlitz-Thouless (BKT) transition in two-dimensional planar rotator and XY models on a square lattice, diluted by randomly placed vacancies, is studied here using hybrid Monte Carlo simulations that combine single spin flip, cluster and over-relaxation techniques. The transition temperature TcT_cTc​ is determined as a function of vacancy density ρvacρ_{vac}ρvac​ by calculations of the helicity modulus and the by finite-size scaling of the in-plane magnetic susceptibility. The results for TcT_cTc​ are consistent with those from the much less precise fourth-order cumulant of Binder. TcT_cTc​ is found to decrease monotonically with increasing ρvacρ_{vac}ρvac​, and falls to zero close to the square lattice percolation limit, ρvac≈0.41ρ_{vac}\approx 0.41ρvac​≈0.41 . The result is physically reasonable: the long-range orientational order of the low-temperature phase cannot be maintained in the absence of sufficient spin interactions across the lattice.
Many problems in statistical physics resist exact solution once the number of interacting particles becomes large, and the two-dimensional Ising model is a classic case where numerical simulation fills that gap. This paper reports a Monte Carlo study of the two-dimensional Ising model using the Metropolis algorithm, undertaken to examine how magnetization, internal energy, specific heat, and magnetic susceptibility change with temperature. Spin configurations were generated on a square lattice under periodic boundary conditions, with each proposed spin flip accepted or rejected according to the Boltzmann weight of the resulting energy change. The simulation reproduces the expected transition from an ordered ferromagnetic state at low temperature to a disordered paramagnetic state at high temperature, with magnetization falling sharply and both specific heat and susceptibility peaking near a critical temperature of Tc ≈ 2.3–2.4 (in units of J/kB), consistent with the exact value obtained analytically for this system. These results confirm that the Metropolis algorithm, despite its simplicity, remains a dependable tool for locating second-order phase transitions and extracting thermodynamic quantities in lattice spin systems, and the same computational framework can be adapted to more elaborate models used in condensed matter and materials research.