Pith. sign in

REVIEW 2 major objections 3 minor 85 references

High-precision Monte Carlo study of several models in the three-dimensional U(1) universality class

T0 review · 2 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Using worm-type Monte Carlo on a torus, this paper reports that wrapping probabilities of directed flows and particle world-lines locate the critical points of the 3D XY, Villain, and Bose-Hubbard models at…

desk verdict High-precision benchmark critical points for the 3D U(1) class, with a clever new use of wrapping probabilities; the error bars rest partly on assumed correction exponents, but the central results look solid. read the letter →

arxiv 1908.10990 v1 pith:7KIDJNXE submitted 2019-08-29 cond-mat.stat-mech cond-mat.quant-gascond-mat.str-elhep-lat

Wanwan Xu , Yanan Sun , Jian-Ping Lv , Youjin Deng This is my paper · ORCID
classification cond-mat.stat-mechcond-mat.quant-gascond-mat.str-elhep-lat MSC 82B2082B2782B80 PACS 05.10.Ln64.60.Fr
keywords 3DXYmodelVillainBose-Hubbardwrappingprobabilitywormalgorithmfinite-sizescalingcriticalexponentsU(1)universalityclass
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Using worm-type Monte Carlo on a torus, this paper argues that wrapping probabilities—whether the directed flow of the classical XY or Villain models, or the particle world-lines of the 2D Bose-Hubbard model, winds around the lattice—are the most reliable high-precision observables for the three-dimensional U(1) universality class. From finite-size scaling of these dimensionless quantities, it determines the critical points $T_c(\mathrm{XY}) = 2.201\,844\,1(5)$, $T_c(\mathrm{Villain}) = 0.333\,067\,04(7)$, and $(t/U)_c(\mathrm{BH}) = 0.059\,729\,1(8)$, improving on earlier estimates significantly. At the Villain critical point, the derivative of the one-direction wrapping probability with respect to temperature shows negligible leading finite-size corrections, which yields the correlation-length exponent $\nu = 0.671\,83(18)$; a susceptibility-like worm-return-time quantity gives $\eta = 0.038\,53(48)$. A sympathetic reader would care because these numbers are benchmarks that can be tested against conformal bootstrap calculations, other Monte Carlo methods, and experiments such as the superfluid transition of helium.

What carries the argument

The central object is the dimensionless wrapping probability $R_\kappa$, equal to one when a directed flow (XY and Villain models) or a particle world-line (Bose-Hubbard model) winds around the torus in direction $\kappa$. Its temperature derivative $G_{R_\kappa E} = \mathrm{d}R_\kappa/\mathrm{d}T$ is estimated from the covariance of $R_\kappa$ with the energy and obeys the finite-size scaling form $G_{R_\kappa E} = L^{1/\nu}(Q_0 + \sum_m b_m L^{-\omega_m})$. The simulations use worm-type updates, which sample the directed-flow or world-line configurations efficiently and let the wrapping numbers be read from the movement of the two defect points; the Bose-Hubbard model is simulated in the imaginary-time world-line representation with $\beta = 2L$. The load-bearing feature is that the leading correction amplitude $b_1$ for $R_x$ is extremely small, so $\nu$ can be determined without fitting that correction amplitude.

What would settle it

At the Villain critical point, refit the $G_{R_x E}$ data to $L^{1/\nu}(Q_0 + b_1 L^{-\omega_1})$ with $b_1$ and $\omega_1$ left free and see whether the resulting $\nu$ stays within $0.671\,83(18)$ when the smallest included lattice size $L_{\min}$ is raised; a stable shift beyond one quoted error would falsify the claim of negligible leading corrections. A second, model-level check is to repeat the Bose-Hubbard critical-point analysis with a different inverse-temperature contour, say $\beta = 4L$ instead of $2L$; if $(t/U)_c$ moves outside $0.059\,729\,1(8)$, the quantum-to-classical scaling assumption would need revision.

Watch

Extended reading notes

Core claim

The central discovery is that the topology of Monte Carlo configurations—measured by which directions the directed flows or world-lines wrap around the periodic boundaries—carries the critical information of the U(1) transition more cleanly than conventional observables. In the Villain model at criticality, $G_{R_x E} = dR_x/dT$, computed as a covariance of the wrapping indicator $R_x$ with the energy, scales as $L^{1/\nu}$ with a leading correction amplitude consistent with zero, while $G_{R_a E}$, $G_{R_2 E}$, and the stiffness derivative carry visible corrections. This lets $\nu = 0.671\,83(18)$ be extracted from fits that omit the leading correction term, avoiding the corresponding parameter uncertainty. The same wrapping observables, together with the superfluid stiffness measured through winding-number fluctuations, locate the three critical points with uncertainty at or below $10^{-7}$, and the critical wrapping probabilities $R_x^c = 0.3787(2)$, $R_a^c = 0.6889(4)$, and $R_2^c = 0.2640(3)$ are found to be universal between the XY and Villain models. In addition, the mean worm return time $T_w$, scaling as $L^{2-\eta}$, gives $\eta = 0.038\,53(48)$.

Load-bearing premise

The fits fix the two correction exponents $\omega_1 = 0.789$ and $\omega_2 = 1.77$ from earlier renormalization-group and Monte Carlo work and assume they also apply to the new wrapping observables and to the Bose-Hubbard world-line data; the second of these is not independently verified here, and the analysis also assumes that no smooth, non-scaling background term enters at the fitted orders.

Editorial extensions

If this is right

  • If correct, $T_c(\mathrm{XY})=2.201\,844\,1(5)$, $T_c(\mathrm{Villain})=0.333\,067\,04(7)$, and $(t/U)_c=0.059\,729\,1(8)$ supersede previous critical-point estimates for all three models.
  • $\nu = 0.671\,83(18)$ is consistent with the earlier Monte Carlo values $0.6717(1)$ and $0.6717(3)$, and the paper's numbers make the space-shuttle helium value $\nu = 0.6709(1)$ unlikely.
  • The critical wrapping probabilities $R_x^c$, $R_a^c$, and $R_2^c$ reported for the XY and Villain models provide new universal dimensionless numbers for this universality class.
  • For the Bose-Hubbard model, the quantum critical point estimate improves on prior results by more than a factor of 40 in precision.
  • The negligible leading corrections in $G_{R_x E}$ provide a practical route to $\nu$ with fewer fitting parameters, useful for other models in the same universality class.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural test is whether $G_{R_x E}$ keeps its near-zero $b_1$ in other realizations of the (2+1)-dimensional U(1) class, such as the quantum rotor model or compact lattice gauge theory; if it does, the route to $\nu$ becomes model-independent and could exceed the current precision.
  • The small corrections may be a feature of the directed-flow representation rather than of the universality class; simulating the XY model in its spin representation with the same wrapping criterion would separate representation effects from universal ones.
  • The reported universal values $R_x^c$, $R_a^c$, and $R_2^c$ could be sharpened by adapting high-precision cluster algorithms to the directed-flow picture, providing independent cross-checks.
  • The same covariance technique could be applied to other geometric observables, such as the probability of double winding or the distribution of winding numbers, to see whether similarly small corrections persist and to yield independent estimates of $\omega_1$.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 3 minor

Summary. This paper presents a worm-type Monte Carlo study of three models in the three-dimensional U(1) universality class: the classical XY model, the Villain model, and the two-dimensional Bose-Hubbard model with unitary filling. From finite-size scaling analyses of wrapping probabilities and the superfluid stiffness, the authors determine the critical points Tc(XY)=2.2018441(5), Tc(Villain)=0.33306704(7), and (t/U)c=0.0597291(8). They further measure the correlation-length exponent nu=0.67183(18) from the temperature derivative of a wrapping probability, which is found to have negligible leading corrections, and the exponent eta=0.03853(48) from a susceptibility-like quantity Tw. The paper introduces wrapping-probability derivatives as a high-precision observable for U(1) criticality and reports universal critical wrapping probabilities.

Significance. If the quoted precision is substantiated, these results provide the most accurate numerical benchmarks for the 3D U(1) universality class to date, improving on the best existing Monte Carlo estimates. The method of using GRxE, the temperature derivative of a wrapping probability, for determining nu is a useful new technique that avoids reliance on correction-to-scaling amplitudes. The paper is exemplary in its transparency: the finite-size scaling fits are documented in full detail, with chi2/dof, Lmin stability checks, and error bars derived from the spread of many fit choices. The consistency of the universal amplitudes across the XY and Villain models (Table V) is a persuasive check. The main weakness is that the quoted systematic errors depend on assumptions about correction-to-scaling exponents and the absence of analytic backgrounds that are only partially tested; this is the focus of the major comments.

major comments (2)
  1. [Secs. III C and IV B (Eq. (32) and Table VI)] The assertion that the analytic background in Tw is 'effectively higher-order' is not demonstrated, and the same issue applies to GRxE. For Tw, the analytic-background exponent relative to the leading behavior is 2-eta about 1.961, which is close to the assumed omega2=1.77; for GRxE, an analytic background would contribute with relative exponent 1/nu about 1.488. If the corresponding amplitude is not tiny, it can bias eta or nu by more than the quoted uncertainties (0.00048 and 0.00018, respectively). I request a quantitative check, for example by adding a term proportional to L^{-(2-eta)} or L^{-1/nu} in Eqs. (32) and (33) and reporting the fitted amplitude, or an equivalent argument showing that the background is negligible.
  2. [Secs. III C and IV A (Tables III, IV, VII, IX, XI)] All finite-size scaling analyses fix the correction exponents omega1=0.789 and omega2=1.77. While omega1 is partially tested (Tables VII and VIII yield omega1 about 0.77(13) and about 0.7), omega2 is never independently verified for the wrapping-probability observables or for the Bose-Hubbard world-line data. Since all fits share the same omega2, a systematic error in this exponent would shift the quoted Tc, nu, and eta coherently, and the spread-based error estimates would not cover the shift. I ask for fits with omega2 treated as free, or with a different subleading correction structure, for at least the Villain-model observables that determine the headline results, to demonstrate that the quoted precision is robust.
minor comments (3)
  1. [Eq. (24) and Eq. (26)] It would be helpful to state explicitly that the correction amplitudes b_m are assumed independent of temperature, so that the derivative with respect to T in Eq. (26) does not act on them; this assumption underlies the claimed negligibility of corrections for GRxE and should be stated in the text.
  2. [Fig. 4] The y-axis label 'Rx-b1L^{-omega1}-b2L^{-omega2}' is somewhat ambiguous; please define in the caption that the plotted quantity is the wrapping probability with the fitted correction terms subtracted.
  3. [Table VII] The header of Table VII lists columns b1, omega1, and b2, but the text says 'omega2=1.77 is adopted'; a note clarifying that omega2 is fixed while omega1 is sometimes free (and sometimes fixed at 0.789) would improve readability.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: critical points and exponents are measured by fitting Monte Carlo data to standard finite-size scaling forms, with correction exponents taken from external RG/MC work and partially re-checked internally.

full rationale

The derivation chain is self-contained. The central results (Tc values, nu, eta) come from least-squares fits of raw worm-Monte-Carlo measurements to standard finite-size scaling forms (Eqs. 24, 28, 32, 33); no reported prediction is also a fit input in a way that forces the answer by construction. The correction exponents omega1=0.789 and omega2=1.77 are adopted from external RG/MC sources ([56] and [11]), not from the authors' own prior work, and omega1 is independently re-estimated in Tables VII and VIII (0.77(13) and approximately 0.7), so the leading correction is not merely imported via self-citation. Using the paper's own high-precision Tc in the exponent fits is a standard two-stage procedure: the Tc errors are roughly two orders of magnitude smaller than the exponent uncertainties, and the fits are stable under variation of L_min; no equation reduces nu or eta to the Tc input. The statement that an analytic background in Tw is 'effectively higher-order' (Sec. III C, after Eq. 32) is an explicitly stated robustness assumption rather than a circular step; if incorrect it would bias eta, but that is a systematic-accuracy concern, not a self-referential reduction. Self-citations ([23], [38], [43], [50]-[52]) support algorithmic methodology and historical comparisons and are not load-bearing. The results are also cross-checked against external benchmarks: Tc values against prior independent Monte Carlo estimates, and nu, eta against [11], [12], and the conformal bootstrap [17]. Hence the paper is self-contained against external data.

Assumptions & free parameters 7 free parameters · 4 assumptions · 0 invented entities

The central claims rest on standard lattice-model Monte Carlo and finite-size scaling. The only external inputs are the correction exponents omega1=0.789 and omega2=1.77, the universality mapping among the three models, and the standard worm-algorithm detailed balance. No new physical entities are introduced. The free parameters are the physical quantities the paper aims to determine; they are fitted from the raw Monte Carlo data rather than imposed by the analysis.

free parameters (7)
  • Tc(XY) = 2.2018441(5)
    Fitted critical temperature of the 3D XY model from finite-size scaling of wrapping probabilities and superfluid stiffness (Tables IX and X).
  • Tc(Villain) = 0.33306704(7)
    Fitted critical temperature of the 3D Villain model from scaling of Rx, Ra, R2, rho_s L, and Tw (Tables III and IV).
  • (t/U)c(BH) = 0.0597291(8)
    Fitted quantum critical point of the 2D unitary-filling Bose-Hubbard model (Tables XI and XII).
  • nu = 0.67183(18)
    Correlation length exponent fitted from the scaling of the wrapping-probability derivative G_RxE at the estimated Tc (Table VI).
  • eta = 0.03853(48)
    Critical exponent fitted from the worm return time Tw at the estimated Tc (Table VII).
  • Universal amplitudes Q0 = Rc_x=0.3787(2), Rc_a=0.6889(4), Rc_2=0.2640(3), rho_s^c L=0.5156(3)
    Critical dimensionless wrapping probabilities and scaled superfluid stiffness fitted for the Villain and XY models (Table V).
  • Per-fit correction amplitudes a1, b1, b2 = Listed individually in Tables III-XII
    Nuisance parameters in the finite-size scaling fits; their values are used for stability checks and do not enter the headline physical conclusions directly.
assumptions (4)
  • domain assumption The 3D XY model, 3D Villain model, and 2D unitary-filling Bose-Hubbard model share the same 3D U(1) universality class with dynamical critical exponent z=1.
    Invoked in Secs. II and III C to combine classical and quantum finite-size scaling; standard in the literature but not derived within this paper.
  • domain assumption Correction-to-scaling exponents omega1=0.789 and omega2=1.77 from RG and prior Monte Carlo apply to all observables, including the new wrapping probabilities and their derivatives.
    Used in fits to Eqs. (24), (26), (28), (30), (32), and (33); partially verified for omega1 in Tables VII and VIII, but omega2 is assumed without independent verification.
  • domain assumption The finite-size scaling ansatze (23)-(32) hold with no dangerous irrelevant variables, and the analytic background in the Tw scaling is negligible.
    Stated in Sec. III C; the analytic-background claim is argued to be higher-order but not proven.
  • domain assumption Wrapping probabilities defined via net flow in the directed-flow and world-line representations satisfy the same universal finite-size scaling as cluster-based wrapping probabilities.
    Assumed in Sec. III B; empirical evidence is presented, but no derivation is given for the extension beyond percolation-type cluster representations.

how reviews work

0 comments
Cite this review

Pith. "Pith review of High-precision Monte Carlo study of several models in the three-dimensional U(1) universality class." pith.science (2026). https://pith.science/paper/7KIDJNXE

@misc{pith2026190810990,
  author       = {Pith},
  title        = {Pith review of: High-precision Monte Carlo study of several models in the three-dimensional U(1) universality class},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7KIDJNXE}},
  note         = {Machine review of arXiv:1908.10990}
}
abstract

We present a worm-type Monte Carlo study of several typical models in the three-dimensional (3D) U(1) universality class, which include the classical 3D XY model in the directed flow representation and its Villain version, as well as the 2D quantum Bose-Hubbard (BH) model with unitary filling in the imaginary-time world-line representation. From the topology of the configurations on a torus, we sample the superfluid stiffness $\rho_s$ and the dimensionless wrapping probability $R$. From the finite-size scaling analyses of $\rho_s$ and of $R$, we determine the critical points as $T_c ({\rm XY}) =2.201\, 844 \,1(5)$ and $T_c ({\rm Villain})=0.333\, 067\, 04(7)$ and $(t/U)_c ({\rm BH})=0.059 \, 729 \,1(8)$, where $T$ is the temperature for the classical models, and $t$ and $U$ are respectively the hopping and on-site interaction strength for the BH model. The precision of our estimates improves significantly over that of the existing results. Moreover, it is observed that at criticality, the derivative of a wrapping probability with respect to $T$ suffers from negligible leading corrections and enables a precise determination of the correlation length critical exponent as $\nu=0.671 \, 83(18)$. In addition, the critical exponent $\eta$ is estimated as $\eta=0.038 \, 53(48)$ by analyzing a susceptibility-like quantity. We believe that these numerical results would provide a solid reference in the study of classical and quantum phase transitions in the 3D U(1) universality, including the recent development of the conformal bootstrap method.

Figures

Figures reproduced from arXiv: 1908.10990 by the authors.

Figure 1
Figure 1. FIG. 1. A directed flow configuration of the XY and Villain mode [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Illustration of wrappings for directed flow configura [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 4
Figure 4. FIG. 4. Scaled [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figures from the paper (3 more)
Figure 5
Figure 5. Figure 5: FIG. 5. Corrections to leading scaling revealed by [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Wrapping probabilities [PITH_FULL_IMAGE:figures/full_fig_p009_6.png]
Figure 8
Figure 8. Figure 8: FIG. 8. (a) Scaled quantity [PITH_FULL_IMAGE:figures/full_fig_p011_8.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

85 extracted references · 80 canonical work pages

  1. [1]

    for the directed flow representation. We begin with the partitio n function ZXY = ( 1 2π )N ∫ e− HXY T ∏ r dθr, (2) 3 where N is the number of lattice sites on the simple-cubic lattice and the exponential e− HXY T can be expanded as e− HXY T = ∏ ⟨rr′⟩ e cos(θ r −θ r′ ) T . (3) Next, we combine ( 2) and ( 3) with the Fourier transform e cos(θ r −θ r′ ) T = ...

    work page

  2. [2]

    III A 1 applies to the Villain model once a substitute of Step 5 is taken as fol- lows

    W orm Algorithm for the Villain Model The worm algorithm formulated in Sec. III A 1 applies to the Villain model once a substitute of Step 5 is taken as fol- lows. Step 5. Accept the proposal with probability Pacc = min(1, e −(C′2 IIN −C2 IIN ) 2T ) according to the Metropolis scheme. The definition of Tw (12) applies to the Villain model as well. B. Sampl...

    work page

  3. [3]

    We simulate the model in canonical ensemble with the worm QMC method within the imaginary-time path-integral representation

    2D Unitary-filling BH Model We extend the applicability of wrapping probability ap- proach to the 2D unitary-filling BH model ( 8), aiming to precisely locate the U(1) QCP . We simulate the model in canonical ensemble with the worm QMC method within the imaginary-time path-integral representation. The simula tions are performed on periodic L × L lattices wi...

    work page

  4. [4]

    MC+HTE 2001 0.671 55(27) 0.038 0(4)

    2001

  5. [5]

    Deconfined quantum critical points,

    T. Senthil, A. Vishwanath, L. Balents, S. Sachdev, and M. P . A. Fisher, “Deconfined quantum critical points,” Science 303, 1490–1494 (2004)

    work page 2004

  6. [6]

    Quantum criticality beyond the landau- ginzburg-wilson paradigm,

    T. Senthil, L. Balents, S. Sachdev, A. Vishwanath, and M. P . A. Fisher, “Quantum criticality beyond the landau- ginzburg-wilson paradigm,” Phys. Rev. B 70, 144407 (2004)

    work page 2004

  7. [7]

    Quantum phase transition from a superfluid to a mot t insulator in a gas of ultracold atoms,

    M. Greiner, O. Mandel, T. Esslinger, T. W. H¨ ansch, and I. Bloch, “Quantum phase transition from a superfluid to a mot t insulator in a gas of ultracold atoms,” Nature 415, 39 (2002)

    work page 2002

  8. [8]

    MC (GPU) 2012 0.671 38(11)

    work page 2012

Show all 85 references
  1. [9]

    MC 2001 0.671 6(5) 0.038 0(5)

    2001

  2. [10]

    Directed geometrical worm algorithm applied to the quantum rotor model,

    F. Alet and E. S. Sørensen, “Directed geometrical worm algorithm applied to the quantum rotor model,” Phys. Rev. E 68, 026702 (2003)

    2003

  3. [11]

    MC+HTE 2006 0.671 7(1) 0.038 1(2)

    2006

  4. [12]

    High-precision measurement of the thermal ex- ponent for the three-dimensional xy universality class,

    E. Burovski, J. Machta, N. Prokofev, and B. Svis- tunov, “High-precision measurement of the thermal ex- ponent for the three-dimensional xy universality class,” Phys. Rev. B 74, 132502 (2006)

    2006

  5. [13]

    MC (GPU) 2012 0.670 98(16)

    2012

  6. [14]

    MC (GPU) 2014 0.672(4)

    2014

  7. [15]

    MC 1993 2.201 67(10)

    1993

  8. [16]

    HTE 2000 0.671 66(55) 0.038 1(3)

    2000

  9. [17]

    Estimates of the critical temperatures for the 3D X Y and Villain models and the critical hopping amplitude for th e two- dimensional (2D) unitary-filling BH model

    CB 2016 0.671 9(11) 0.038 52(64) this work MC 2019 0.671 83(18) 0.038 53(48) TABLE II. Estimates of the critical temperatures for the 3D X Y and Villain models and the critical hopping amplitude for th e two- dimensional (2D) unitary-filling BH model. ‘SCE’ is the abbr evia- ti...

    2016

  10. [18]

    2003 0.670 9(1)

    Exp. 2003 0.670 9(1)

    2003

  11. [19]

    MC (GPU) 2012 2.201 831 2(6)

    2012

  12. [20]

    MC (GPU) 2012 2.201 852(1)

    2012

  13. [21]

    MC (GPU) 2014 2.201 836(6) this work MC 2019 2.201 844 1(5) Villain

    2014

  14. [22]

    W orm Algorithm for the XY model Extending state space.— A character of worm algorithm is enlarging state space. It extends the original directed fl ow space ∆ C = 0 by including two additional degrees of free- dom, namely, two defects I and M individually on a site, by defining...

  15. [23]

    MC 2014 0.333 067 0(2) this work MC 2019 0.333 067 04(7) BH

    2014

  16. [24]

    Sampled Quantities for the 3D XY and Villain Models Some wrapping-related quantities for 3D XY and Villain models are defined as follows. The wrapping probabilities in the directed flow representation are given by Rx = ⟨Rx⟩ = ⟨Ry⟩ = ⟨Rz⟩, (13) Ra = 1 − ⟨(1 − Rx)(1 − Ry)(1 − Rz)⟩...

  17. [25]

    For a given spatial direction (say κ), the event wrapping (namely, Rκ = 1 ) relates to a non-zero winding number Wκ ⁄= 0 of particle lines

    Sampled Quantities for the 2D BH Model For the 2D BH model, the wrapping probabilities of particle lines in the world-line representation read Rx = ⟨Rx⟩ = ⟨Ry⟩, (19) Ra = 1 − ⟨(1 − Rx)(1 − Ry)⟩, (20) R2 = ⟨RxRy⟩, (21) where Rx, Ra and R2 define the probabilities that the wrap- ...

  18. [26]

    conventional

    and (32). For the present case, this analytic background is effectively higher-order corrections compared with the correc- tion terms taken into account explicitly throughout this wo rk. IV . NUMERICAL RESULTS AND FINITE-SIZE SCALING ANALYSES In this section, we present Monte ...

  19. [27]

    conventional

    3D Villain Model We simulate the 3D Villain model on periodic L × L × L simple-cubic lattices with linear sizes L = 16 , 24, 32, 64, 128, 256, 384 and 512 for different T around T = 0 .333 067. The simulations use the worm Monte Carlo algorithm de- scribed in Sec. III A 2. The...

  20. [28]

    SCE 1993 0.059 74(4)

    1993

  21. [29]

    Albeit these estimates are all based on Monte Carlo simula- tions, they are not completely consistent with each other

    MC 2008 0.059 74(3) this work MC 2019 0.059 729 1(8) 2.201 831 2(6) [13], 2.201 852(1) [13] and 2.201 836(6) [14]. Albeit these estimates are all based on Monte Carlo simula- tions, they are not completely consistent with each other. We perform a high-temperature expansion on model (

    2008

  22. [30]

    1996 0.670 19(13)

    Exp. 1996 0.670 19(13)

    1996

  23. [31]

    MC 2005 2.201 840 5(48)

    2005

  24. [32]

    MC 1996 2.201 843(19)

    1996

  25. [33]

    MC 2002 2.201 833(19)

    2002

  26. [34]

    3D XY model By means of the worm algorithm formulated in Sec. III A 1, we simulate the 3D XY model ( 1) on periodic simple-cubic lattices with linear sizes L = 8 , 16, 32, 64, 128, 256, 384 0.375 0.378 0.381 Rx-b1L-ω 1-b2L-ω 2 L T=0.3330667 T=0.333067039 T=0.3330674 0.685 0.69...

  27. [35]

    Ordering, metasta- bility and phase transitions in two-dimensional systems,

    J. M. Kosterlitz and D. J. Thouless, “Ordering, metasta- bility and phase transitions in two-dimensional systems,” J. Phys. C 6, 1181 (1973)

    1973

  28. [36]

    B. V . Svistunov, E. S. Babaev, and N. V . Prokof’ev, Superfluid states of matter (CRC Press, 2015)

    2015

  29. [37]

    Pyrochlore p ho- tons: The u(1) spin liquid in a s=1/2 three-dimensional frus - trated magnet,

    M. Hermele, M. P . A. Fisher, and L. Balents, “Pyrochlore p ho- tons: The u(1) spin liquid in a s=1/2 three-dimensional frus - trated magnet,” Phys. Rev. B 69, 064404 (2004)

    2004

  30. [38]

    Stability of u(1) spin liquids in two dimen- sions,

    M. Hermele, T. Senthil, M. P . A. Fisher, P . A. Lee, N. Nagaosa, and X.-G. Wen, “Stability of u(1) spin liquids in two dimen- sions,” Phys. Rev. B 70, 214437 (2004)

    2004

  31. [39]

    L. D. Faddeev, Gauge fields: an introduction to quantum theory (CRC Press, 2018)

    2018

  32. [40]

    Critical behavior of the three-dimensional xy u ni- versality class,

    M. Campostrini, M. Hasenbusch, A. Pelissetto, P . Rossi, and E. Vicari, “Critical behavior of the three-dimensional xy u ni- versality class,” Phys. Rev. B 63, 214503 (2001)

    2001

  33. [41]

    Theoretical estimates of the critical exponent s of the superfluid transition in he4 by lattice methods,

    M. Campostrini, M. Hasenbusch, A. Pelissetto, and E. Vicari, “Theoretical estimates of the critical exponent s of the superfluid transition in he4 by lattice methods,” Phys. Rev. B 74, 144506 (2006)

    2006

  34. [42]

    High-precision monte carlo study of the three-dimensional xy model on gpu,

    T.-Y . Lan, Y .-D. Hsieh, and Y .-J. Kao, “High-precision monte carlo study of the three-dimensional xy model on gpu,” arXiv preprint arXiv:1211.0780 (2012)

    2012 arXiv

  35. [43]

    Cuda programs for the gpu computing of the swendsen–wang multi-cluster spin flip algorithm: 2d and 3d ising, potts, and xy models,

    Y . Komura and Y . Okabe, “Cuda programs for the gpu computing of the swendsen–wang multi-cluster spin flip algorithm: 2d and 3d ising, potts, and xy models,” Comp. Phys. Commu. 185, 1038–1043 (2014)

    2014

  36. [44]

    Criti- cal behaviourof the 3d xy-model: A monte carlo study,

    A. P . Gottlob, M. Hasenbusch, and S. Meyer, “Criti- cal behaviourof the 3d xy-model: A monte carlo study,” Nucl. Phys. B (PS) 30, 838–841 (1993)

    1993

  37. [45]

    Deter- mination of the critical exponents for the λ transition of 4he by high-temperature expansion,

    M. Campostrini, A. Pelissetto, P . Rossi, and E. Vicari, “Deter- mination of the critical exponents for the λ transition of 4he by high-temperature expansion,” Phys. Rev. B 61, 5905 (2000)

    2000

  38. [46]

    Precision islands in the ising and o(n) models,

    F. Kos, D. Poland, D. Simmons-Duffin, and A. Vichi, “Precision islands in the ising and o(n) models,” JHEP 08, 036 (2016)

    2016

  39. [47]

    Specific heat of liquid helium in zero gravity very near the lambda point,

    J. A. Lipa, J. A. Nissen, D. A. Stricker, D. R. Swanson, an d T. C. P . Chui, “Specific heat of liquid helium in zero gravity very near the lambda point,” Phys. Rev. B 68, 174518 (2003)

    2003

  40. [48]

    Goldman, Percolation, localization, and superconductivity , V ol

    A. Goldman, Percolation, localization, and superconductivity , V ol. 109 (Springer Science & Business Media, 2013)

    2013

  41. [49]

    Nonuniversal critical d ynam- ics in monte carlo simulations,

    R. H. Swendsen and J.-S. Wang, “Nonuniversal critical d ynam- ics in monte carlo simulations,” Phys. Rev. Lett. 58, 86 (1987)

    1987

  42. [50]

    Collective monte carlo updating for spin sys tems,

    U. Wolff, “Collective monte carlo updating for spin sys tems,” Phys. Rev. Lett. 62, 361 (1989)

    1989

  43. [51]

    Worm algorithms for classical statistical models,

    N. V . Prokof’ev and B. V . Svistunov, “Worm algorithms for classical statistical models,” Phys. Rev. Lett. 87, 160601 (2001)

    2001

  44. [52]

    Un i- versal conductivity in a two-dimensional superfluid-to-in sulator quantum critical system,

    K. Chen, L. Liu, Y . Deng, L. Pollet, and N. Prokof’ev, “Un i- versal conductivity in a two-dimensional superfluid-to-in sulator quantum critical system,” Phys. Rev. Lett. 112, 030402 (2014)

    2014

  45. [53]

    The dy- namics of quantum criticality revealed by quantum monte car lo and holography,

    W. Witczak-Krempa, E. S. Sørensen, and S. Sachdev, “The dy- namics of quantum criticality revealed by quantum monte car lo and holography,” Nat. Phys. 10, 361 (2014)

    2014

  46. [54]

    Non- standard hubbard models in optical lattices: a review,

    O. Dutta, M. Gajda, P . Hauke, M. Lewenstein, D. L¨ uhmann , B. A. Malomed, T. Sowi´ nski, and J. Zakrzewski, “Non- standard hubbard models in optical lattices: a review,” Rep. Prog. Phys. 78, 066001 (2015)

    2015

  47. [55]

    Suppression o f the critical temperature for superfluidity near the mott tra nsi- tion,

    S. Trotzky, L. Pollet, F. Gerbier, U. Schnorrberger, I. Bloch, N. V . Prokof’ev, B. Svistunov, and M. Troyer, “Suppression o f the critical temperature for superfluidity near the mott tra nsi- tion,” Nat. Phys. 6, 998 (2010)

    2010

  48. [56]

    Observation of correlated particle-hole pairs and string order in low-dimensional mott insulators,

    M. Endres, M. Cheneau, T. Fukuhara, C. Weitenberg, P . Schauß, C. Gross, L. Mazza, M. C. Banuls, L. Pol- let, I. Bloch, et al. , “Observation of correlated particle-hole pairs and string order in low-dimensional mott insulators, ” Science 334, 200–203 (2011)

    2011

  49. [57]

    Dynamics and thermodynamics of the bose-hubbard model,

    N. Elstner and H. Monien, “Dynamics and thermodynamics of the bose-hubbard model,” Phys. Rev. B 59, 12184 (1999)

    1999

  50. [58]

    Monte carlo study of the two-dimensional bose- hubbard model,

    B. Capogrosso-Sansone, S ¸ . S¨ oyler, N. V . Prokof’ev, and B. V . 13 TABLE IX. Fits of the wrapping probabilities Rx, Ra, R2 to (24) and the scaled SF stiffness ρsL to (28) for the 3D XY model. The correction exponents ω 1 = 0. 789 and ω 2 = 1. 77 are adopted. Qua. Lmin χ 2/...

    2008

  51. [59]

    Heat capacity and thermal relaxation of bulk helium very near the lambda point,

    J. A. Lipa, D. R. Swanson, J. A. Nissen, T. C. P . Chui, and U. E. Israelsson, “Heat capacity and thermal relaxation of bulk helium very near the lambda point,” Phys. Rev. Lett. 76, 944 (1996)

    1996

  52. [60]

    Sur- face and bulk transitions in three-dimensional o(n) models ,

    Y . Deng, H. W. J. Bl¨ ote, and M. P . Nightingale, “Sur- face and bulk transitions in three-dimensional o(n) models ,” Phys. Rev. E 72, 016128 (2005)

    2005

  53. [61]

    Finite size effects on measures of critical ex po- nents in d=3 o (n) models,

    H. G. Ballesteros, L. A. Fernandez, V . Martin-Mayor, an d A. M. Sudupe, “Finite size effects on measures of critical ex po- nents in d=3 o (n) models,” Phys. Lett. B 387, 125–131 (1996)

    1996

  54. [62]

    Universal amplitude ratios from nu- merical studies of the three-dimensional o(2) model,

    A. Cucchieri, J. Engels, S. Holtmann, T. Mendes, and T. Schulze, “Universal amplitude ratios from nu- merical studies of the three-dimensional o(2) model,” J. Phys. A: Math. and Gen. 35, 6517 (2002)

    2002

  55. [63]

    Cluster monte carlo algorit hm for the quantum rotor model,

    F. Alet and E. S. Sørensen, “Cluster monte carlo algorit hm for the quantum rotor model,” Phys. Rev. E 67, 015701 (2003)

    2003

  56. [64]

    Universal conductivity of two- dimensional films at the superconductor-insulator transit ion,

    M. C. Cha, M. P . A. Fisher, S. M. Girvin, M. Wallin, and A. P . Y oung, “Universal conductivity of two- dimensional films at the superconductor-insulator transit ion,” Phys. Rev. B 44, 6883 (1991)

    1991

  57. [65]

    Universal scaling of the conductivity at the superfluid-insulator phase transition ,

    J. ˇSmakov and E. Sørensen, “Universal scaling of the conductivity at the superfluid-insulator phase transition ,” Phys. Rev. Lett. 95, 180603 (2005)

    2005

  58. [66]

    Nishimori and G

    H. Nishimori and G. Ortiz, Elements of phase transitions and critical phenomena (OUP Oxford, 2010)

    2010

  59. [67]

    Dynamic crit- ical behavior of the worm algorithm for the ising model,

    Y . Deng, T. M. Garoni, and A. D. Sokal, “Dynamic crit- ical behavior of the worm algorithm for the ising model,” Phys. Rev. Lett. 99, 110601 (2007)

    2007

  60. [68]

    Ex - act, complete, and universal continuous-time worldline mo nte carlo approach to the statistics of discrete quantum system s,

    N. V . Prokof’ev, B. V . Svistunov, and I. S. Tupitsyn, “Ex - act, complete, and universal continuous-time worldline mo nte carlo approach to the statistics of discrete quantum system s,” JETP 87, 310–321 (1998)

    1998

  61. [69]

    worm algorithm in quantum monte carlo simulations,

    N. V . Prokof’ev, B. V . Svistunov, and I. S. Tupit- syn, “worm algorithm in quantum monte carlo simulations,” Phys. Lett. A 238, 253–257 (1998)

    1998

  62. [70]

    Worm algo- rithm for problems of quantum and classical statistics,

    N. V . Prokof’ev and B. V . Svistunov, “Worm algo- rithm for problems of quantum and classical statistics,” arXiv preprint arXiv:0910.1393 (2009)

    2009 arXiv

  63. [71]

    Recent developments in quantum monte carlo simulations with applications for cold gases,

    L. Pollet, “Recent developments in quantum monte carlo simulations with applications for cold gases,” 14 TABLE X. Fits of the wrapping probabilities Rx, Ra, R2 to (24) and the scaled SF stiffness ρsL to (28) for the 3D XY model. The critical exponent ν is fixed at our final esti...

    2012

  64. [72]

    Worm-type monte carl o simulation of the ashkin-teller model on the triangular lat tice,

    J.-P . Lv, Y . Deng, and Q.-H. Chen, “Worm-type monte carl o simulation of the ashkin-teller model on the triangular lat tice,” Phys. Rev. E 84, 021125 (2011)

    2011

  65. [73]

    On the universality of crossing probabilities in two-dimensi onal percolation,

    R. P . Langlands, C. Pichet, P . Pouliot, and Y . Saint-Aubin, “On the universality of crossing probabilities in two-dimensi onal percolation,” J. Stat. Phys. 67, 553–574 (1992)

    1992

  66. [74]

    Critical percolation on the torus,

    H. T. Pinson, “Critical percolation on the torus,” J. Stat. Phys. 75, 1167–1177 (1994)

    1994

  67. [75]

    Shape-dependen t universality in percolation,

    R. M. Ziff, C. D. Lorenz, and P . Kleban, “Shape-dependen t universality in percolation,” Physica A 266, 17–26 (1999)

    1999

  68. [76]

    Fast monte carlo algorith m for site or bond percolation,

    M. E. J. Newman and R. M. Ziff, “Fast monte carlo algorith m for site or bond percolation,” Phys. Rev. E 64, 016706 (2001)

    2001

  69. [77]

    Homology of fortuin–kasteleyn clusters of potts models on the torus,

    L. P . Arguin, “Homology of fortuin–kasteleyn clusters of potts models on the torus,” J. Stat. Phys. 109, 301–310 (2002)

    2002

  70. [78]

    Percolation on two- a nd three-dimensional lattices,

    P . H. L. Martins and J. A. Plascak, “Percolation on two- a nd three-dimensional lattices,” Phys. Rev. E 67, 046119 (2003)

    2003

  71. [79]

    Bond and site percolation in three dimensions,

    J. Wang, Z. Zhou, W. Zhang, T. M. Garoni, and Y . Deng, “Bond and site percolation in three dimensions,” Phys. Rev. E 87, 052107 (2013)

    2013

  72. [80]

    Simulta- neous analysis of three-dimensional percolation models,

    X. Xu, J. Wang, J.-P . Lv, and Y . Deng, “Simulta- neous analysis of three-dimensional percolation models,” Front. Phys. 9, 113–119 (2014)

    2014

  73. [81]

    Geometric properties of the fortuin-kasteleyn representation of the ising model,

    P . Hou, S. Fang, J. Wang, H. Hu, and Y . Deng, “Geometric properties of the fortuin-kasteleyn representation of the ising model,” Phys. Rev. E 99, 042150 (2019)

    2019

  74. [82]

    Path-integral comput ation of superfluid densities,

    E. L. Pollock and D. M. Ceperley, “Path-integral comput ation of superfluid densities,” Phys. Rev. B 36, 8343 (1987)

    1987

  75. [83]

    Boson localization and the superfluid-insulator transiti on,

    Weichman P . B. Grinstein G. Fisher, M. P . A. and D. S. Fish er, “Boson localization and the superfluid-insulator transiti on,” Phys. Rev. B 40, 546 (1989)

    1989

  76. [84]

    Strongly disor- dered superfluids: Quantum fluctuations and critical behavi or,

    M. Ma, B. I. Halperin, and P . A. Lee, “Strongly disor- dered superfluids: Quantum fluctuations and critical behavi or,” Phys. Rev. B 34, 3136 (1986)

    1986

  77. [85]

    Critical exponents of the n-vector model,

    R. Guida and J. Zinn-Justin, “Critical exponents of the n-vector model,” J Phys. A: Math. and Gen. 31, 8103 (1998) . 15 TABLE XI. Fits of the wrapping probabilities Rx, Ra, R2 to (24) and the scaled SF stiffness ρsL to (28) for the 2D unitary-filling BH model. The correction ex...

    1998

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.