pith. sign in

arxiv: 2312.12130 · v2 · submitted 2023-12-19 · ✦ hep-lat · nucl-ex· nucl-th

Finite-size behavior of higher-order cumulant ratios near criticality in two-dimensional Potts models

Rajiv V. Gavai , Bedangadas Mohanty , Jaydev Singh Rao , Swati Saha This is my paper

Pith reviewed 2026-05-24 05:40 UTC · model grok-4.3

classification ✦ hep-lat nucl-exnucl-th
keywords cumulant ratiosPotts modelssecond-order phase transitionsfinite-size effectscritical fluctuationsQCD phase diagramMonte Carlo simulations
0
0 comments X

The pith

Cumulant ratios in two-dimensional Potts models do not realize the QCD-predicted hierarchy except in a narrow window above the critical temperature.

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

The paper tests whether a hierarchy of cumulant ratios predicted for the QCD phase transition appears generically in other finite systems near second-order transitions. It does this by performing Monte Carlo simulations of the two-state and three-state Potts models in two dimensions, calculating magnetization cumulants up to sixth order on lattices of different sizes. The simulations show that the expected ordering among the ratios holds only in a narrow temperature range just above the critical temperature. Outside this range, the hierarchy is not observed, nor is its reverse, across the temperatures studied.

Core claim

Theoretical considerations predict a hierarchy χ6/χ2 < χ5/χ1 < χ4/χ2 < χ3/χ1 among net-baryon number cumulant ratios near the QCD transition at small baryon chemical potential. Monte Carlo simulations on two-dimensional Potts models find that neither this complete hierarchy nor its exact reverse is realized over the studied temperature range in either model, except within a narrow temperature window above the critical temperature.

What carries the argument

Ratios of higher-order cumulants of the total magnetization in finite-size two-dimensional Potts models simulated near their critical temperatures.

If this is right

  • Higher-order cumulants exhibit extrema whose magnitudes increase with cumulant order and lattice size.
  • The hierarchy appears only in a narrow window above the critical temperature.
  • This behavior is consistent across both the two-state and three-state Potts models.
  • Finite-size effects enhance critical fluctuations reflected in the cumulants.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The hierarchy may be specific to the QCD transition rather than a general feature of second-order phase transitions in finite systems.
  • Studies in three dimensions or other spin models could test if dimensionality affects the emergence of the ordering.
  • If the finding holds, it suggests that cumulant hierarchies are not reliable generic indicators of criticality in all statistical systems.

Load-bearing premise

The cumulant ratios from magnetization in two-dimensional Potts models act as a faithful proxy for whether the QCD hierarchy emerges in any finite system at a second-order phase transition.

What would settle it

Finding the full hierarchy persisting across a broad temperature range above criticality in these or similar models would falsify the main claim.

Figures

Figures reproduced from arXiv: 2312.12130 by Bedangadas Mohanty, Jaydev Singh Rao, Rajiv V. Gavai, Swati Saha.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) The second order susceptibility of magnetization [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) The second order susceptibility of magnetization [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The third and fourth order susceptibilities of magnetization, [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The fifth and sixth order susceptibilities of magnetization, [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Ratios of the susceptibilities, [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Ratios of the susceptibilities, [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
read the original abstract

Theoretical considerations predict a specific hierarchy among ratios of net-baryon number cumulants ($\chi_n$, where $n$ is the order of cumulant) in the vicinity of the transition from the low-temperature hadronic phase to the high temperature quark-gluon plasma phase at small baryon chemical potential, $\mu_\mathrm{B}$, in the QCD phase diagram. This hierarchy, $\frac{\chi_6}{\chi_2} < \frac{\chi_5}{\chi_1} < \frac{\chi_4}{\chi_2} < \frac{\chi_3}{\chi_1}$, has been observed by the STAR experiment in net-proton number (a proxy of net-baryon number) cumulant ratios over a broad range of collision energies. Motivated by these findings, we investigate whether similar ordering emerges generically in finite statistical systems undergoing second-order phase transitions. We employ two different spin models: the two-state and three-state Potts models in two dimensions, both exhibiting a transition from an ordered phase to a disordered phase at their respective critical temperatures. Monte Carlo simulations are performed on square lattices of varying sizes using the Wolff cluster algorithm. Cumulants of the total magnetization are calculated up to sixth order in both of these models in a temperature range near their corresponding critical temperatures. Higher-order cumulants exhibit extrema (peaks/troughs) whose magnitudes grow with both cumulant order and lattice size, reflecting enhanced critical fluctuations. Except within a narrow temperature window above the critical temperature, neither the complete hierarchy nor its exact reverse is realized over the studied temperature range in either 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

1 major / 0 minor

Summary. The manuscript reports Monte Carlo simulations of the two-dimensional two-state and three-state Potts models using the Wolff cluster algorithm to compute cumulants of the magnetization up to sixth order near their critical points on finite lattices. The central finding is that the hierarchy of cumulant ratios χ6/χ2 < χ5/χ1 < χ4/χ2 < χ3/χ1 predicted for QCD near the critical point is not observed in these models except in a narrow temperature window above Tc, and neither is the reverse hierarchy realized over the broader temperature range studied.

Significance. Should the numerical results hold, the work indicates that the observed ordering in net-proton cumulants from STAR does not arise as a generic consequence of finite-size effects near a second-order phase transition in simple spin models. This could point to the importance of conserved charges or other QCD-specific features in producing the hierarchy. The study employs standard, well-established methods for 2D Potts models, providing a clear numerical test of the generic hypothesis.

major comments (1)
  1. [Abstract] Abstract (paragraph 2): The claim that magnetization cumulants in these Potts models test whether the QCD hierarchy 'emerges generically in finite statistical systems undergoing second-order phase transitions' is undermined by the mismatch in fluctuation type. QCD χ_n are obtained from derivatives of the pressure w.r.t. the conserved baryon chemical potential, whereas the Potts cumulants are of the non-conserved order parameter (magnetization) under an external field. This distinction is load-bearing for the generic interpretation, as the two classes of fluctuations have different scaling and conservation properties.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment and the constructive comment on the abstract. We agree that the distinction between fluctuation types merits clarification and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Abstract] Abstract (paragraph 2): The claim that magnetization cumulants in these Potts models test whether the QCD hierarchy 'emerges generically in finite statistical systems undergoing second-order phase transitions' is undermined by the mismatch in fluctuation type. QCD χ_n are obtained from derivatives of the pressure w.r.t. the conserved baryon chemical potential, whereas the Potts cumulants are of the non-conserved order parameter (magnetization) under an external field. This distinction is load-bearing for the generic interpretation, as the two classes of fluctuations have different scaling and conservation properties.

    Authors: We acknowledge the referee's point that the cumulants studied here are those of a non-conserved order parameter (magnetization) in the presence of an external field, whereas the QCD cumulants χ_n derive from derivatives of the pressure with respect to a conserved chemical potential. Our study employs the 2D Potts models as standard, well-studied examples of second-order transitions to test whether the reported hierarchy can arise from finite-size effects near criticality in a simple statistical system. Nevertheless, to prevent any ambiguity in the generic interpretation, we will revise the abstract (and, if appropriate, the introduction) to explicitly note that the fluctuations considered are those of a non-conserved order parameter. This change clarifies the scope without altering the numerical findings or the conclusion that the hierarchy is not realized over most of the temperature range studied. revision: yes

Circularity Check

0 steps flagged

No circularity: direct Monte Carlo computation of cumulant ratios

full rationale

The paper performs explicit Monte Carlo simulations with the Wolff algorithm on finite 2D Potts lattices, extracts cumulants of total magnetization up to order 6, and reports their ratios as a function of temperature and volume. The central claim (neither hierarchy nor reverse realized except in a narrow window above Tc) is an empirical observation from these computed values, not a derivation that reduces to a fitted parameter, self-definition, or self-citation chain. The interpretive step linking Potts magnetization cumulants to QCD net-baryon cumulants is an external-motivation assumption, not a load-bearing circular step inside the reported results.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper relies on standard statistical mechanics without introducing new free parameters, axioms beyond domain conventions, or invented entities.

axioms (1)
  • domain assumption Cumulants of magnetization in Potts models can serve as a proxy for testing generic features of net-baryon cumulant ratios near second-order transitions.
    Invoked in the motivation paragraph to justify the choice of models.

pith-pipeline@v0.9.0 · 5840 in / 1389 out tokens · 35947 ms · 2026-05-24T05:40:27.089538+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

34 extracted references · 34 canonical work pages

  1. [1]

    The three-state Potts model that belongs to the same universality class as the two-dimensional Ising model, undergoes a second order phase transition at TC = 1 ln(1+ √

    J /kB ≈ 2.2692 J /kB [23]. The three-state Potts model that belongs to the same universality class as the two-dimensional Ising model, undergoes a second order phase transition at TC = 1 ln(1+ √

    work page

  2. [2]

    J /kB ≈ 0.9950 J /kB [22]. A. Cumulants In this paper, the main objects of interest are the higher-order cumulants (also referred as susceptibilities in the paper), which reflect the higher-order correlations amongst the spins. In spin models, the pth order suscep- tibility normalized by the total number of spin sites in the lattice is defined [24] by χp(...

    work page

  3. [3]

    lattice volume

    By fitting χ2 as a Gaussian function of the temper- ature, T , the peak location ( TC,L, χ2(TC,L)) is deter- mined for each lattice size. The obtained values of TC,L and χ2(TC,L) as a function of L are tabulated in Tables I and II for q = 2 and q = 3 state Potts model respec- tively. Also, we define the critical region around TC,L for each lattice size by...

    work page

  4. [4]

    Edward V. Shuryak. Quantum Chromodynamics and the Theory of Superdense Matter. Phys. Rept. , 61:71–158, 1980

    work page 1980

  5. [5]

    Thermodynamics of strong-interaction matter from Lattice QCD

    Heng-Tong Ding, Frithjof Karsch, and Swagato Mukher- jee. Thermodynamics of strong-interaction matter from Lattice QCD. Int. J. Mod. Phys. E , 24(10):1530007, 2015

    work page 2015

  6. [6]

    Katz, Stefan Krieg, and Kalman K

    Szabocls Borsanyi, Zoltan Fodor, Christian Hoelbling, Sandor D. Katz, Stefan Krieg, and Kalman K. Szabo. Full result for the QCD equation of state with 2+1 fla- vors. Phys. Lett. B , 730:99–104, 2014

    work page 2014

  7. [7]

    Decoding the phase struc- ture of QCD via particle production at high energy

    Anton Andronic, Peter Braun-Munzinger, Krzysztof Redlich, and Johanna Stachel. Decoding the phase struc- ture of QCD via particle production at high energy. Na- ture, 561(7723):321–330, 2018

    work page 2018

  8. [8]

    Y. Aoki, G. Endrodi, Z. Fodor, S. D. Katz, and K. K. Sz- abo. The Order of the quantum chromodynamics transi- tion predicted by the standard model of particle physics. Nature, 443:675–678, 2006

    work page 2006

  9. [9]

    Color super- conductivity and chiral symmetry restoration at nonzero baryon density and temperature

    Juergen Berges and Krishna Rajagopal. Color super- conductivity and chiral symmetry restoration at nonzero baryon density and temperature. Nucl. Phys. B , 538:215– 232, 1999

    work page 1999

  10. [10]

    M. A. Halasz, A. D. Jackson, R. E. Shrock, M. A. Stephanov, and J. J. M. Verbaarschot. Phase diagram of qcd. Phys. Rev. D , 58:096007, Sep 1998

    work page 1998

  11. [11]

    Asakawa and K

    M. Asakawa and K. Yazaki. Chiral Restoration at Finite Density and Temperature. Nucl. Phys. A , 504:668–684, 1989

    work page 1989

  12. [12]

    Pandav, D

    A. Pandav, D. Mallick, and B. Mohanty. Search for the QCD critical point in high energy nuclear collisions.Prog. Part. Nucl. Phys. , 125:103960, 2022

    work page 2022

  13. [13]

    Mapping the Phases of Quantum Chromodynamics with Beam Energy Scan

    Adam Bzdak, Shinichi Esumi, Volker Koch, Jinfeng Liao, Mikhail Stephanov, and Nu Xu. Mapping the Phases of Quantum Chromodynamics with Beam Energy Scan. Phys. Rept., 853:1–87, 2020

    work page 2020

  14. [14]

    Beam Energy Dependence of Fifth and Sixth-Order Net-proton Number Fluctua- tions in Au+Au Collisions at RHIC

    Bassam Aboona et al. Beam Energy Dependence of Fifth and Sixth-Order Net-proton Number Fluctua- tions in Au+Au Collisions at RHIC. Phys. Rev. Lett. , 130(8):082301, 2023

    work page 2023

  15. [15]

    Bazavov, D

    A. Bazavov, D. Bollweg, H.-T. Ding, P. Enns, J. Goswami, P. Hegde, O. Kaczmarek, F. Karsch, R. Larsen, Swagato Mukherjee, H. Ohno, P. Petreczky, C. Schmidt, S. Sharma, and P. Steinbrecher. Skewness, kurtosis, and the fifth and sixth order cumulants of net baryon-number distributions from lattice qcd confront high-statistics star data. Phys. Rev. D , 101:0...

    work page 2020

  16. [16]

    Bazavov, H.-T

    A. Bazavov, H.-T. Ding, P. Hegde, O. Kaczmarek, F. Karsch, E. Laermann, Swagato Mukherjee, H. Ohno, P. Petreczky, E. Rinaldi, H. Sandmeyer, C. Schmidt, Chris Schroeder, S. Sharma, W. Soeldner, R. A. Soltz, P. Steinbrecher, and P. M. Vranas. Skewness and kurto- sis of net baryon-number distributions at small values of the baryon chemical potential. Phys. R...

    work page 2017

  17. [17]

    Pawlowski, Fabian Rennecke, Rui Wen, and Shi Yin

    Wei-jie Fu, Xiaofeng Luo, Jan M. Pawlowski, Fabian Rennecke, Rui Wen, and Shi Yin. Hyper-order baryon number fluctuations at finite temperature and density. Phys. Rev. D , 104:094047, Nov 2021

    work page 2021

  18. [18]

    Guenther, San- dor K

    Szabolcs Borsanyi, Zoltan Fodor, Jana N. Guenther, San- dor K. Katz, Kalman K. Szabo, Attila Pasztor, Israel Portillo, and Claudia Ratti. Higher order fluctuations and correlations of conserved charges from lattice QCD. JHEP, 10:205, 2018

    work page 2018

  19. [19]

    Benjamin Svetitsky and Laurence G. Yaffe. Critical Be- havior at Finite Temperature Confinement Transitions. Nucl. Phys. B , 210:423–447, 1982

    work page 1982

  20. [20]

    R. V. Gavai, F. Karsch, and B. Petersson. A Study of the Correlation Length Near a First Order Phase Tran- sition: The 3- D Three State Potts Model. Nucl. Phys. B, 322:738–758, 1989

    work page 1989

  21. [21]

    L. G. Yaffe and B. Svetitsky. First-order phase transition in the su(3) gauge theory at finite temperature. Phys. Rev. D, 26:963–965, Aug 1982

    work page 1982

  22. [22]

    Mean field analysis of su (n) deconfining transitions in the presence of dynamical quarks

    F Green and Frithjof Karsch. Mean field analysis of su (n) deconfining transitions in the presence of dynamical quarks. Nuclear Physics B , 238(2):297–306, 1984

    work page 1984

  23. [23]

    H. W. J. Bl¨ ote and R. H. Swendsen. First-order phase transitions and the three-state potts model. Phys. Rev. Lett., 43:799–802, Sep 1979

    work page 1979

  24. [24]

    The Three- dimensional, three state Potts model in an external field

    Frithjof Karsch and Sven Stickan. The Three- dimensional, three state Potts model in an external field. Phys. Lett. B , 488:319–325, 2000

    work page 2000

  25. [25]

    F. Y. Wu. The Potts model. Reviews of Modern Physics , 54(1):235–268, January 1982

    work page 1982

  26. [26]

    Crystal statistics

    Lars Onsager. Crystal statistics. i. a two-dimensional model with an order-disorder transition. Phys. Rev. , 65:117–149, Feb 1944

    work page 1944

  27. [27]

    R. K. Pathria and Paul D. Beale. Statistical Mechanics. Elsevier, 3rd edition, 2011

    work page 2011

  28. [28]

    R. G. M. Rodrigues, B. V. Costa, and L. A. S. M´ ol. Moment-generating function zeros in the study of phase transitions. Phys. Rev. E , 104:064103, Dec 2021

    work page 2021

  29. [29]

    A recursive formulation of the old problem of obtaining moments from cumulants and vice versa

    Peter J Smith. A recursive formulation of the old problem of obtaining moments from cumulants and vice versa. The American Statistician , 49(2):217–218, 1995

    work page 1995

  30. [30]

    Privman and M

    V. Privman and M. E. Fisher. Universal critical ampli- tudes in finite-size scaling. Phys. Rev. B , 30:322–327, Jul 1984

    work page 1984

  31. [31]

    Pathria and Paul D

    R.K. Pathria and Paul D. Beale. 14 - phase transitions: The renormalization group approach. In R.K. Pathria and Paul D. Beale, editors, Statistical Mechanics (Third Edition), pages 539–581. Academic Press, Boston, third edition edition, 2011

    work page 2011

  32. [32]

    Critical phenomena and renormalization group theory

    Andrea Pelissetto and Ettore Vicari. Critical phenomena and renormalization group theory. Phys. Rept., 368:549– 727, 2002

    work page 2002

  33. [33]

    Collective monte carlo updating for spin sys- tems

    Ulli Wolff. Collective monte carlo updating for spin sys- tems. Phys. Rev. Lett. , 62:361–364, Jan 1989

    work page 1989

  34. [34]

    M. E. J. Newman and G. T. Barkema. Monte Carlo methods in statistical physics . Clarendon Press, Oxford, 1999

    work page 1999