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:1010.3075 · Physical Review E

Berezinskii-Kosterlitz-Thouless-like percolation transitions in the two-dimensional XY model

Hao Hu, Youjin Deng, Henk W. J. Blöte

Published Jan 25, 2011• 8 citations
arXivHTMLPDF

Abstract

We study a percolation problem on a substrate formed by two-dimensional XY spin configurations using Monte Carlo methods. For a given spin configuration, we construct percolation clusters by randomly choosing a direction x in the spin vector space, and then placing a percolation bond between nearest-neighbor sites i and j with probability p(ij)=max(0,1-e(-2Ks(i)(x)s(j)(x))), where K>0 governs the percolation process. A line of percolation thresholds K(c)(J) is found in the low-temperature range J≥J(c), where J>0 is the XY coupling strength. Analysis of the correlation function g(p)(r), defined as the probability that two sites separated by a distance r belong to the same percolation cluster, yields algebraic decay for K≥K(c)(J), and the associated critical exponent depends on J and K. Along the threshold line K(c)(J), the scaling dimension for g(p) is, within numerical uncertainties, equal to 1/8. On this basis, we conjecture that the percolation transition along the K(c)(J) line is of the Berezinskii-Kosterlitz-Thouless type.

DOI
10.1103/physreve.83.011124

Related papers

Conducting-angle-based percolation in the XY model

Yan-Cheng Wang, Wenan Guo, Bernard Nienhuis +1 more

We define a percolation problem on the basis of spin configurations of the two-dimensional XY model. Neighboring spins belong to the same percolation cluster if their orientations differ less than a certain threshold called the conducting angle. The percolation properties of this model are studied by means of Monte Carlo simulations and a finite-size scaling analysis. Our simulations show the existence of percolation transitions when the conducting angle is varied, and we determine the transition point for several values of the XY coupling. It appears that the critical behavior of this percolation model can be well described by the standard percolation theory. The critical exponents of the percolation transitions, as determined by finite-size scaling, agree with the universality class of the two-dimensional percolation model on a uniform substrate. This holds over the whole temperature range, even in the low-temperature phase where the XY substrate is critical in the sense that it displays algebraic decay of correlations.

arXiv:0912.3343•Physical Review E•Mar 17, 2010•3 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

Critical behavior of the fully frustrated two dimensional XY model

Anderson Barbosa Lima, B. V. Costa

Using Monte Carlo simulations we have investigated the critical behavior of the classical fully frustrated XY model in two dimensions in a square lattice. There are two phase transitions in the model, one of the Berezinskii-Kosterlitz-Thouless type and a Z2 transition at higher temperature. We show that the vortex anti-vortex density has a clear signature of the Z2 phase transition at exactly the percolation threshold.

Brazilian Journal of Physics•Jun 1, 2004•3 citations

Percolation between vacancies in the two-dimensional Blume-Capel model

Youjin Deng, Wenan Guo, Henk W. J. Blöte

Using suitable Monte Carlo methods and finite-size scaling, we investigate the Blume-Capel model on the square lattice. We construct percolation clusters by placing nearest-neighbor bonds between vacancies with a variable bond probability p(b) . At the tricritical point, we locate the percolation threshold of these vacancy clusters at p(bc) =0.706 33 (6) . At this point, we determine the fractal dimension of the vacancy clusters as Xf =0.1308 (5) approximately equal to 21/160, and the exponent governing the renormalization flow in the p(b) direction as y(p) =0.426 (2) approximately equal to 17/40 . For bond probability p(b) > p(bc) , the vacancy clusters maintain strong critical correlations; the fractal dimension is Xf =0.0750 (2) approximately equal to 3/40 and the leading correction exponent is y(p) =-0.45 (2) approximately equal to -19/40 . The above values fit well in the Kac table for the tricritical Ising model. These vacancy clusters have much analogy with those consisting of Ising spins of the same sign, although the associated quantities rho and magnetization m are energylike and magnetic quantities, respectively. However, along the critical line of the Blume-Capel model, the vacancies are more or less uniformly distributed over the whole lattice. In this case, no critical percolation correlations are observed in the vacancy clusters, at least in the physical region p(b) < or = 1 .

Physical Review E•Jul 1, 2005•17 citations

Spin-spin correlation function of the two-dimensional XY model with weak site or bond dilution

O. Kapikranian

The spin-spin correlation function of the two-dimensional (2D) $XY$ model decays as a power law at all temperatures below the Berezinskii-Kosterlitz-Thouless transition point with a temperature-dependent exponent $\ensuremath{\eta}=\ensuremath{\eta}(T/J)$ ($J$ is the ferromagnetic coupling strength). It is known from computer experiments that in the 2D $XY$ model with site or bond dilution this exponent depends on the concentration $p$ of removed sites or bonds as well. Knowing the slope $\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p$ at point $p=0$, one can predict the value of the exponent for small dilution concentrations: ${\ensuremath{\eta}(p)\ensuremath{\simeq}\ensuremath{\eta}(0)+p(\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p)|}_{p=0}$. As shown in this paper, the spin-wave Hamiltonian allows one to obtain exact results for this slope: ${(\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p)|}_{p=0}=T/(2J)+O({(T/J)}^{2})$ and $T/(\ensuremath{\pi}J)+O({(T/J)}^{2})$ for site and for bond dilution, respectively.

arXiv:1111.3939•Physical Review B•Mar 5, 2012•3 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
We study a percolation problem on a substrate formed by two-dimensional XY spin configurations using Monte Carlo methods. For a given spin configuration, we construct percolation clusters by randomly choosing a direction x in the spin vector space, and then placing a percolation bond between nearest-neighbor sites i and j with probability p(ij)=max(0,1-e(-2Ks(i)(x)s(j)(x))), where K>0 governs the percolation process. A line of percolation thresholds K(c)(J) is found in the low-temperature range J≥J(c), where J>0 is the XY coupling strength. Analysis of the correlation function g(p)(r), defined as the probability that two sites separated by a distance r belong to the same percolation cluster, yields algebraic decay for K≥K(c)(J), and the associated critical exponent depends on J and K. Along the threshold line K(c)(J), the scaling dimension for g(p) is, within numerical uncertainties, equal to 1/8. On this basis, we conjecture that the percolation transition along the K(c)(J) line is of the Berezinskii-Kosterlitz-Thouless type.
We define a percolation problem on the basis of spin configurations of the two-dimensional XY model. Neighboring spins belong to the same percolation cluster if their orientations differ less than a certain threshold called the conducting angle. The percolation properties of this model are studied by means of Monte Carlo simulations and a finite-size scaling analysis. Our simulations show the existence of percolation transitions when the conducting angle is varied, and we determine the transition point for several values of the XY coupling. It appears that the critical behavior of this percolation model can be well described by the standard percolation theory. The critical exponents of the percolation transitions, as determined by finite-size scaling, agree with the universality class of the two-dimensional percolation model on a uniform substrate. This holds over the whole temperature range, even in the low-temperature phase where the XY substrate is critical in the sense that it displays algebraic decay of correlations.
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.
Using suitable Monte Carlo methods and finite-size scaling, we investigate the Blume-Capel model on the square lattice. We construct percolation clusters by placing nearest-neighbor bonds between vacancies with a variable bond probability p(b) . At the tricritical point, we locate the percolation threshold of these vacancy clusters at p(bc) =0.706 33 (6) . At this point, we determine the fractal dimension of the vacancy clusters as Xf =0.1308 (5) approximately equal to 21/160, and the exponent governing the renormalization flow in the p(b) direction as y(p) =0.426 (2) approximately equal to 17/40 . For bond probability p(b) > p(bc) , the vacancy clusters maintain strong critical correlations; the fractal dimension is Xf =0.0750 (2) approximately equal to 3/40 and the leading correction exponent is y(p) =-0.45 (2) approximately equal to -19/40 . The above values fit well in the Kac table for the tricritical Ising model. These vacancy clusters have much analogy with those consisting of Ising spins of the same sign, although the associated quantities rho and magnetization m are energylike and magnetic quantities, respectively. However, along the critical line of the Blume-Capel model, the vacancies are more or less uniformly distributed over the whole lattice. In this case, no critical percolation correlations are observed in the vacancy clusters, at least in the physical region p(b) < or = 1 .
The spin-spin correlation function of the two-dimensional (2D) XYXYXY model decays as a power law at all temperatures below the Berezinskii-Kosterlitz-Thouless transition point with a temperature-dependent exponent $\ensuremath{\eta}=\ensuremath{\eta}(T/J)$ (JJJ is the ferromagnetic coupling strength). It is known from computer experiments that in the 2D XYXYXY model with site or bond dilution this exponent depends on the concentration ppp of removed sites or bonds as well. Knowing the slope $\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p$ at point p=0p=0p=0, one can predict the value of the exponent for small dilution concentrations: ${\ensuremath{\eta}(p)\ensuremath{\simeq}\ensuremath{\eta}(0)+p(\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p)|}_{p=0}$. As shown in this paper, the spin-wave Hamiltonian allows one to obtain exact results for this slope: ${(\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p)|}_{p=0}=T/(2J)+O({(T/J)}^{2})$ and $T/(\ensuremath{\pi}J)+O({(T/J)}^{2})$ for site and for bond dilution, respectively.
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.