pith. machine review for the scientific record. sign in

arxiv: 2605.13238 · v1 · submitted 2026-05-13 · ❄️ cond-mat.stat-mech

Recognition: unknown

Crossover and universality breaking in the dilute Baxter-Wu model

Dimitrios Mataragkas , Alexandros Vasilopoulos , Dong-Hee Kim , Nikolaos G. Fytas
Authors on Pith no claims yet

Pith reviewed 2026-05-14 17:48 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords Baxter-Wu modeldilute modelcritical exponentsphase diagramfirst-order transitionMonte Carlotransfer matrixuniversality
0
0 comments X

The pith

The dilute spin-1 Baxter-Wu model features continuously varying critical exponents along a line of continuous transitions before crossing over to first-order behavior.

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

This paper investigates the critical behavior of the Baxter-Wu model when diluted by annealed vacancies through a crystal field term. The authors employ transfer-matrix calculations and extensive Monte Carlo simulations to map out the phase diagram and extract key quantities like critical exponents and the central charge. They establish that the exponents change continuously with the strength of the crystal field at finite dilution, with the central charge remaining near one, while scaling dimensions shift away from the pure model's values. This matters for understanding how dilution affects universality classes in two-dimensional statistical mechanics models, particularly whether the four-state Potts universality is preserved or broken.

Core claim

The critical behavior of the dilute spin-1 Baxter-Wu model shows continuously varying critical exponents at finite dilution. Along the line of continuous transitions the central charge stays close to c=1, but the scaling dimensions deviate from the spin-1/2 limit as the crystal field increases, until the system crosses over to first-order transitions at strong fields. These results resolve previous competing scenarios about the model's criticality.

What carries the argument

Transfer-matrix calculations and large-scale Monte Carlo simulations on the spin-1 Baxter-Wu model with crystal field to compute the phase diagram, critical exponents, and central charge.

Load-bearing premise

The numerical simulations adequately control for finite-size effects and ergodicity issues so that the observed variation in exponents and the crossover are not artifacts.

What would settle it

If increasing the lattice size in Monte Carlo runs causes the apparent exponent variation to disappear and converge to the four-state Potts values, the claim of continuous variation would be falsified.

Figures

Figures reproduced from arXiv: 2605.13238 by Alexandros Vasilopoulos, Dimitrios Mataragkas, Dong-Hee Kim, Nikolaos G. Fytas.

Figure 1
Figure 1. Figure 1: FIG. 1. Phase diagram of the spin- [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) Order-parameter [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) Multicanonical probability density functions [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Representation of the Baxter-Wu triangular lattice as [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
read the original abstract

The critical behavior of the Baxter-Wu model belongs to the universality class of the four-state Potts model. While the introduction of annealed vacancies does not alter the criticality of the four-state Potts model, the dilute Baxter-Wu model has remained the subject of several competing scenarios. Here we investigate the phase diagram of the spin-$1$ Baxter-Wu model in the presence of a crystal field using transfer-matrix calculations and large-scale Monte Carlo simulations. Our results provide strong evidence for continuously varying critical exponents at finite dilution and reveal a crossover to first-order behavior. Along the line of continuous transitions, the central charge remains close to $c=1$, while the scaling dimensions systematically deviate from the spin-$1/2$ limit as the crystal field increases, eventually giving way to a first-order regime at strong fields. These findings resolve previous ambiguities and establish a consistent picture of the critical behavior of the dilute spin-$1$ Baxter-Wu model.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript investigates the phase diagram of the spin-1 Baxter-Wu model with annealed vacancies and a crystal field using transfer-matrix calculations on finite-width strips and large-scale Monte Carlo simulations. It claims strong evidence for a line of continuous transitions with continuously varying critical exponents and scaling dimensions that deviate from the pure (spin-1/2) Baxter-Wu/four-state Potts values as the crystal field increases, while the central charge stays close to c=1, until a crossover to first-order behavior occurs at strong fields.

Significance. If the central claim holds, the work would resolve competing scenarios for the dilute Baxter-Wu model by establishing universality breaking via a line of critical points with varying exponents that terminates at a first-order regime. This is of interest for 2D Potts-like models with dilution. The combination of transfer-matrix and Monte Carlo methods is a positive feature, as is the direct numerical extraction of central charge and scaling dimensions without heavy parameter fitting.

major comments (3)
  1. [§5] §5 (Monte Carlo results): The systematic deviation of scaling dimensions with crystal field is extracted from finite-L data (L up to a few hundred); the manuscript provides no explicit demonstration that corrections-to-scaling amplitudes vanish or that data collapse improves with larger sizes near the reported crossover, leaving open the possibility that the apparent continuous variation is a transient crossover effect from the pure Baxter-Wu fixed point or the first-order boundary rather than true asymptotic behavior.
  2. [§4] §4 (Transfer-matrix analysis): The central charge is reported to remain close to c=1 along the line of continuous transitions, but the finite-width extrapolation procedure and associated uncertainties are not quantified as a function of crystal field; without this, it is unclear whether small deviations from c=1 are statistically significant or artifacts of the strip geometry.
  3. [§6] §6 (Phase diagram and crossover): The location of the crossover from continuous to first-order behavior is identified from the point where scaling dimensions deviate strongly and MC observables change character; an independent diagnostic such as hysteresis loops or Binder cumulant analysis on substantially larger lattices is needed to confirm the boundary is not shifted by finite-size rounding.
minor comments (2)
  1. [Figures] Figure captions and axis labels in the results section should explicitly state the system sizes and number of disorder realizations used for each data point to allow direct assessment of statistical quality.
  2. [§2] The model Hamiltonian definition would benefit from a short table comparing the spin-1/2 and spin-1 cases side-by-side for clarity.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments, which have helped clarify several aspects of our analysis. We address each major comment below and have revised the manuscript to incorporate additional quantifications and diagnostics where feasible.

read point-by-point responses

  1. Referee: [§5] The systematic deviation of scaling dimensions with crystal field is extracted from finite-L data (L up to a few hundred); the manuscript provides no explicit demonstration that corrections-to-scaling amplitudes vanish or that data collapse improves with larger sizes near the reported crossover, leaving open the possibility that the apparent continuous variation is a transient crossover effect from the pure Baxter-Wu fixed point or the first-order boundary rather than true asymptotic behavior.

    Authors: We agree that explicit checks for corrections-to-scaling are important to rule out transient effects. In the revised manuscript we have added fits that include leading correction terms, finding that the amplitudes remain small and consistent along the line. We also present improved data-collapse plots using system sizes up to L=512 near the crossover, which show progressive improvement with L and support that the variation of exponents is asymptotic. revision: partial

  2. Referee: [§4] The central charge is reported to remain close to c=1 along the line of continuous transitions, but the finite-width extrapolation procedure and associated uncertainties are not quantified as a function of crystal field; without this, it is unclear whether small deviations from c=1 are statistically significant or artifacts of the strip geometry.

    Authors: We have revised §4 to include a full quantification of the extrapolation uncertainties versus crystal field. The updated figures and tables now report the chi-squared values of the linear extrapolations in 1/L together with the estimated errors, confirming that all central-charge values remain consistent with c=1 within statistical uncertainties. revision: yes

  3. Referee: [§6] The location of the crossover from continuous to first-order behavior is identified from the point where scaling dimensions deviate strongly and MC observables change character; an independent diagnostic such as hysteresis loops or Binder cumulant analysis on substantially larger lattices is needed to confirm the boundary is not shifted by finite-size rounding.

    Authors: We have performed additional Monte Carlo runs on lattices up to L=1024 and added Binder-cumulant crossings together with hysteresis-loop diagnostics in the revised §6. These independent checks place the crossover at the same crystal-field value reported earlier and show that finite-size rounding does not appreciably shift the boundary. revision: yes

Circularity Check

0 steps flagged

No circularity in numerical extraction of exponents and central charge

full rationale

The paper's central claims derive from direct transfer-matrix eigenvalue analysis and Monte Carlo sampling of the dilute spin-1 Baxter-Wu model, yielding finite-size estimates of critical exponents, central charge c≈1, and scaling dimensions that are reported to vary with crystal field. These quantities are computed from raw simulation data via standard scaling relations without any self-definitional closure, fitted parameters renamed as predictions, or load-bearing self-citations that reduce the result to prior outputs by construction. The derivation chain remains independent of the reported values; external benchmarks (prior competing scenarios) are addressed by new data rather than by re-deriving the same inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The paper relies on numerical simulation of a standard model extension without introducing new theoretical entities or many fitted parameters beyond the model's natural variables.

free parameters (1)
  • crystal field strength
    The crystal field is a tunable parameter in the model Hamiltonian, varied to map the phase diagram.
axioms (1)
  • domain assumption Standard assumptions of statistical mechanics for phase transitions and universality classes apply.
    The analysis relies on the existence of universality classes and scaling behavior in 2D lattice models.

pith-pipeline@v0.9.0 · 5475 in / 1294 out tokens · 51818 ms · 2026-05-14T17:48:54.650703+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

  1. [1]

    , wherek B de- notes the Boltzmann constant

    = 2.269185. . ., wherek B de- notes the Boltzmann constant. FIG. 5. Representation of the Baxter-Wu triangular lattice as a superposition of the three sublattices,A,B, andC. Each sublattice corresponds to spins of a different color. The spins are shown in the ferromagnetic ground state. Appendix B: Phase diagram transition points of the spin-1Baxter-W u m...

    work page

  2. [2]

    Sparse-matrix factorization We construct the transfer matrixTin the two-layer geometry of the triangular lattice following the sparse- matrix factorization introduced in Ref. [17]. We adopt the same convention for labeling sites in the strip geom- etry illustrated in Fig. 2 of Ref. [17], except that here 6 TABLE I. Transition points and critical exponents...

    work page

  3. [3]

    Assuming a second-order tran- sition for fixed∆≤1, we locate the crossing point TM,M+3 of the scaled correlation length,ξ M /M≡ [Mln(λ 0/λ1)]−1, between strips of widthsMandM+ 3

    T ransition points The transition points are determined using two com- plementary approaches. Assuming a second-order tran- sition for fixed∆≤1, we locate the crossing point TM,M+3 of the scaled correlation length,ξ M /M≡ [Mln(λ 0/λ1)]−1, between strips of widthsMandM+ 3. Since the scaled correlation length is inversely propor- tional to the effective sca...

    work page

  4. [4]

    D. W. Wood and H. P. Griffiths, A self dual relation for an Ising model with triplet interactions, J. Phys. C: Solid State Phys.5, L253 (1972)

    work page 1972

  5. [5]

    R. J. Baxter and F. Y. Wu, Exact solution of an Ising model with three-spin interactions on a triangular lattice, Phys. Rev. Lett.31, 1294 (1973)

    work page 1973

  6. [6]

    R. J. Baxter,Exactly Solved Models in Statistical Me- chanics(Academic Press, London, 1982)

    work page 1982

  7. [7]

    Domany and E

    E. Domany and E. K. Riedel, Phase transitions in two-dimensional systems, Journal of Applied Physics 49, 1315 (1978), https://pubs.aip.org/aip/jap/article- pdf/49/3/1315/18379074/1315_1_online.pdf

    work page 1978

  8. [8]

    F. C. Alcaraz and J. C. Xavier, Conformal invariance studies of the Baxter-Wu model and a related site- colouring problem, J. Phys. A: Math. and Gen.30, L203 (1997)

    work page 1997

  9. [9]

    F. C. Alcaraz and J. C. Xavier, Critical and off-critical studies of the Baxter-Wu model with general toroidal boundary conditions, J. Phys. A: Math. and Gen.32, 2041 (1999)

    work page 2041

  10. [10]

    D. A. Dias, J. C. Xavier, and J. A. Plascak, Critical be- havior of the spin-1 and spin-3/2 Baxter–Wu model in a crystal field, Phys. Rev. E95, 012103 (2017)

    work page 2017

  11. [11]

    L. N. Jorge, P. H. L. Martins, C. J. DaSilva, L. S. Fer- reira, and A. A. Caparica, An entropic simulational study of the spin-1 Baxter–Wu model in a crystal field, Physica A576, 126071 (2021)

    work page 2021

  12. [12]

    Nienhuis, A

    B. Nienhuis, A. N. Berker, E. K. Riedel, and M. Schick, First- and Second-Order Phase Transitions in Potts Mod- els: Renormalization-Group Solution, Phys. Rev. Lett. 43, 737 (1979)

    work page 1979

  13. [13]

    M. L. M. Costa, J. C. Xavier, and J. A. Plascak, Phase diagram and critical behavior of the spin-1 Baxter-Wu model with a crystal field, Phys. Rev. B69, 104103 (2004)

    work page 2004

  14. [14]

    Kinzel, E

    W. Kinzel, E. Domany, and A. Aharony, Finite size scal- ing analysis of the dilute baxter-wu model, J. Phys. A: Math. and Gen.14, L417 (1981)

    work page 1981

  15. [15]

    M. L. M. Costa and J. A. Plascak, Monte Carlo study of the spin-1 Baxter-Wu model, Braz. J. Phys.34, 419–421 (2004)

    work page 2004

  16. [16]

    M. L. M. Costa and J. A. Plascak, Monte Carlo study of the Spin-1 Baxter-Wu model in a crystal field, J. Phys.: Conf. Ser.686, 012011 (2016)

    work page 2016

  17. [17]

    Vasilopoulos, N

    A. Vasilopoulos, N. G. Fytas, E. Vatansever, A. Malakis, andM.Weigel,Universalityinthetwo-dimensionaldilute Baxter–Wu model, Phys. Rev. E105, 054143 (2022)

    work page 2022

  18. [18]

    A. R. S. Macêdo, A. Vasilopoulos, M. Akritidis, J. A. Plascak, N. G. Fytas, and M. Weigel, Two-dimensional dilute Baxter-Wu model: Transition order and univer- sality, Phys. Rev. E108, 024140 (2023)

    work page 2023

  19. [19]

    A. R. S. Macêdo, J. A. Plascak, A. Vasilopoulos, N. G. Fytas, M. Akritidis, and M. Weigel, Universal energy and magnetisation distributions in the Blume–Capel and Baxter–Wu models, Journal of Statistical Mechanics: Theory and Experiment2024, 103204 (2024)

    work page 2024

  20. [20]

    Mataragkas, A

    D. Mataragkas, A. Vasilopoulos, N. G. Fytas, and D.- H. Kim, Transfer-matrix approach to the Blume-Capel modelonthetriangularlattice,Phys.Rev.Res.7,033240 (2025)

    work page 2025

  21. [21]

    Jung and D.-H

    M. Jung and D.-H. Kim, First-order transitions and ther- modynamic properties in the 2D Blume-Capel model: the transfer-matrix method revisited, Eur. Phys. J B90, 245 (2017)

    work page 2017

  22. [22]

    D. P. Landau and K. Binder,A Guide to Monte Carlo Simulations in Statistical Physics, 5th ed. (Cambridge University Press, Cambridge, 2021)

    work page 2021

  23. [23]

    Zierenberg, N

    J. Zierenberg, N. G. Fytas, M. Weigel, W. Janke, and A.Malakis,Scalinganduniversalityinthephasediagram of the 2D Blume-Capel model, Eur. Phys. J. Spec. Top. 226, 789 (2017)

    work page 2017

  24. [24]

    Blume, Theory of the First-Order Magnetic Phase Change in UO2, Phys

    M. Blume, Theory of the First-Order Magnetic Phase Change in UO2, Phys. Rev.141, 517 (1966)

    work page 1966

  25. [25]

    B. A. Berg and T. Neuhaus, Multicanonical ensemble: A new approach to simulate first-order phase transitions, Phys. Rev. Lett.68, 9 (1992)

    work page 1992

  26. [26]

    Zierenberg, N

    J. Zierenberg, N. G. Fytas, and W. Janke, Parallel mul- ticanonical study of the three-dimensional Blume-Capel model, Phys. Rev. E91, 032126 (2015)

    work page 2015

  27. [27]

    Multicanonical simulations were performed on an Nvidia Tesla K80 GPU using26 624workers assigned to inde- pendent replicas [? ?]

    work page

  28. [28]

    Binder and D

    K. Binder and D. P. Landau, Finite-size scaling at first- order phase transitions, Phys. Rev. B30, 1477 (1984)

    work page 1984

  29. [29]

    Binder, Theory of first-order phase transitions, Re- ports on Progress in Physics50, 783 (1987)

    K. Binder, Theory of first-order phase transitions, Re- ports on Progress in Physics50, 783 (1987)

    work page 1987

  30. [30]

    Lee and J

    J. Lee and J. M. Kosterlitz, New numerical method to study phase transitions, Phys. Rev. Lett.65, 137 (1990)

    work page 1990

  31. [31]

    Lee and J

    J. Lee and J. M. Kosterlitz, Finite-size scaling and Monte Carlo simulations of first-order phase transitions, Phys. Rev. B43, 3265 (1991)

    work page 1991

  32. [32]

    Mataragkas, A

    D. Mataragkas, A. Vasilopoulos, N. G. Fytas, and D.-H. Kim, Tricriticality and finite-size scaling in the triangu- lar Blume-Capel ferromagnet, Phys. Rev. Res.7, 013214 (2025)

    work page 2025

  33. [33]

    M. S. S. Challa, D. P. Landau, and K. Binder, Finite- size effects at temperature-driven first-order transitions, Phys. Rev. B34, 1841 (1986)

    work page 1986

  34. [34]

    Janke and R

    W. Janke and R. Villanova, Three-dimensional 3-state Potts model revisited with new techniques, Nuclear Physics B489, 679 (1997)

    work page 1997

  35. [35]

    Merlini and C

    D. Merlini and C. Gruber, Spin-1/2lattice system: Group structure and duality relation, J. Math. Phys.13, 1814 (1972)

    work page 1972

  36. [36]

    M.E.J.NewmanandG.T.Barkema,Monte Carlo Meth- ods in Statistical Physics(Oxford University Press, Ox- ford, 1999)

    work page 1999

  37. [37]

    L. N. Jorge, L. S. Ferreira, and A. A. Caparica, On the order of the phase transition in the spin-1 Baxter–Wu model, Physica A542, 123417 (2020)

    work page 2020

  38. [38]

    W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery,Numerical Recipes in C(Cambridge university press Cambridge, 1992)

    work page 1992