pith. machine review for the scientific record. sign in

arxiv: 2605.08543 · v1 · submitted 2026-05-08 · ✦ hep-ph

Recognition: no theorem link

Scaling properties of net-baryon number fluctuations at the deconfinement critical point

Authors on Pith no claims yet

Pith reviewed 2026-05-12 01:23 UTC · model grok-4.3

classification ✦ hep-ph
keywords net-baryon number fluctuationsdeconfinement critical pointcritical exponentsPolyakov loop modelmean-field approximationGinzburg-Landau criterion3D Ising modelheavy quark limit
0
0 comments X

The pith

Net-baryon fluctuations follow mean-field scaling near the critical point

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

The paper investigates the scaling of the first four cumulants of net-baryon number near the deconfinement critical point in the heavy quark limit of QCD. It links these fluctuations to Polyakov loop susceptibilities and derives mean-field critical exponents from Landau theory, confirming them with numerical results in an effective model. A Ginzburg-Landau criterion is used to show that the critical region becomes smaller as baryon density increases. Beyond mean-field effects are estimated by applying the scaling function from the 3D Ising model.

Core claim

In the mean-field approximation, critical exponents for net-baryon cumulants are obtained from Landau theory and verified through numerical computations in the Polyakov loop model. The Ginzburg-Landau criterion indicates that the critical region shrinks with increasing baryon chemical potential. Critical exponents outside the mean-field regime are estimated using the exact scaling function of the 3D Ising model.

What carries the argument

The connection between baryon-number fluctuations and Polyakov loop susceptibilities, together with Landau theory for deriving exponents and the Ginzburg-Landau criterion for assessing the critical region.

If this is right

  • The cumulants of net-baryon number follow specific scaling laws determined by mean-field theory near the critical point.
  • The size of the critical region in the phase diagram decreases as the baryon density rises.
  • Beyond-mean-field critical exponents can be approximated by those of the 3D Ising universality class via its scaling function.
  • These properties hold at both zero and nonzero baryon chemical potentials in the heavy-quark limit.

Where Pith is reading between the lines

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

  • This implies that searches for the critical point in heavy-ion collisions may need to focus on narrower ranges of parameters at higher densities.
  • The method could be applied to other models or observables to test universality in QCD phase transitions.
  • If validated, it provides a way to estimate fluctuation behaviors without full non-perturbative calculations.

Load-bearing premise

The effective Polyakov loop model in the heavy-quark limit faithfully reproduces the critical behavior of full QCD, and the Ginzburg-Landau criterion reliably estimates the critical region without higher-order corrections.

What would settle it

Numerical simulations of full QCD at finite baryon density that measure the net-baryon cumulants and compare the observed scaling exponents or the extent of the critical region to the predictions would falsify or confirm the results.

Figures

Figures reproduced from arXiv: 2605.08543 by Chihiro Sasaki, Krzysztof Redlich, Micha{\l} Szyma\'nski, Pok Man Lo.

Figure 1
Figure 1. Figure 1: Green solid lines correspond to the crossover with [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The same as in Fig. 1 but for [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Fourth order baryon number cumulant in the proximity of the critical point in function of quark mass and temperature [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Third order baryon number susceptibility in the prox [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. 4th order susceptibility of net-baryon number in the [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: Red squares correspond to the approach from the [PITH_FULL_IMAGE:figures/full_fig_p007_9.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. 2nd (red squares), 3rd (green dots) and 4th (blue triangles) order susceptibilities of net-baryon number in the proximity [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. The 4th order susceptibility of the net-baryon number [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. 2nd (red squares), 3rd (green dots) and 4th (blue triangles) order susceptibilities of net-baryon number in the [PITH_FULL_IMAGE:figures/full_fig_p009_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Estimate of critical region based on GL ratio for [PITH_FULL_IMAGE:figures/full_fig_p009_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Contours of the GL ratio for [PITH_FULL_IMAGE:figures/full_fig_p010_12.png] view at source ↗
read the original abstract

We investigate the critical behavior of the first four cumulants of the net-baryon number near the deconfinement critical point in QCD in the limit of heavy quarks. By connecting baryon-number fluctuations to Polyakov loop susceptibilities, we analyze their mean-field scaling properties at zero and non-zero baryon chemical potentials. In the mean-field approximation we derive critical exponents via the Landau theory and validate them through explicit numerical calculations in an effective Polyakov loop model. Using a Ginzburg-Landau criterion as a diagnostic of beyond-mean-field effects, we estimate the size of the critical region and find that it shrinks with increasing baryon density. By utilizing the exact scaling function in the 3D Ising model, we estimate critical exponents beyond the mean-field approximation.

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

2 major / 2 minor

Summary. The paper investigates the critical behavior of the first four cumulants of net-baryon number near the deconfinement critical point in the heavy-quark limit of QCD. It connects baryon-number fluctuations to Polyakov-loop susceptibilities, derives mean-field critical exponents from Landau theory and validates them numerically in an effective Polyakov-loop model, applies a Ginzburg-Landau criterion to estimate the size of the critical region (finding shrinkage with increasing baryon density), and uses the exact 3D Ising scaling function to estimate exponents beyond mean-field.

Significance. If the central results hold, the work provides a useful mapping of net-baryon cumulants onto Polyakov-loop observables and a concrete estimate of how the critical region contracts with baryon chemical potential. The explicit numerical validation of Landau-derived mean-field exponents and the direct import of the 3D Ising scaling function are methodological strengths that could inform searches for critical behavior in heavy-ion collisions. The shrinkage claim, however, rests on an approximation whose mu dependence has not been cross-checked.

major comments (2)
  1. Ginzburg-Landau criterion (the section applying the criterion to bound the critical region): The criterion equates the quartic coefficient to the one-loop fluctuation integral but omits higher-order cumulants and the mu-dependent shift in the fluctuation spectrum induced by the Polyakov-loop background. This approximation is load-bearing for the claim that the critical region shrinks with increasing baryon density; the mu dependence has not been cross-validated, so the reported shrinkage may be overestimated.
  2. Numerical validation in the effective Polyakov-loop model (the section presenting the explicit calculations): The mean-field exponents are confirmed within the model, but the model is restricted to the heavy-quark limit. The central claim that the results inform full QCD critical behavior therefore depends on an unquantified assumption that this effective model faithfully reproduces the relevant critical properties at finite density.
minor comments (2)
  1. Abstract: The statement that the critical region 'shrinks with increasing baryon density' should be accompanied by a brief qualifier that this follows from the Ginzburg-Landau estimate, to avoid overstatement.
  2. Notation: The mapping between net-baryon cumulants and Polyakov-loop susceptibilities is introduced without an explicit equation relating the two sets of observables; adding this relation early would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and have made revisions to clarify limitations and strengthen the presentation where appropriate.

read point-by-point responses
  1. Referee: Ginzburg-Landau criterion (the section applying the criterion to bound the critical region): The criterion equates the quartic coefficient to the one-loop fluctuation integral but omits higher-order cumulants and the mu-dependent shift in the fluctuation spectrum induced by the Polyakov-loop background. This approximation is load-bearing for the claim that the critical region shrinks with increasing baryon density; the mu dependence has not been cross-validated, so the reported shrinkage may be overestimated.

    Authors: We agree that the Ginzburg-Landau criterion is an approximation that neglects higher-order cumulants and the full mu-dependent shift in the fluctuation spectrum due to the Polyakov-loop background. This is the standard diagnostic used in effective models to estimate the critical region size, with the leading mu dependence arising from the explicit factors in the Landau potential and the one-loop integral. We have added a dedicated paragraph in the revised manuscript discussing these limitations and arguing that the qualitative shrinkage trend with increasing baryon density remains robust under the leading-order treatment. A complete cross-validation would require a more elaborate model incorporating additional degrees of freedom, which lies outside the present scope. revision: partial

  2. Referee: Numerical validation in the effective Polyakov-loop model (the section presenting the explicit calculations): The mean-field exponents are confirmed within the model, but the model is restricted to the heavy-quark limit. The central claim that the results inform full QCD critical behavior therefore depends on an unquantified assumption that this effective model faithfully reproduces the relevant critical properties at finite density.

    Authors: The effective Polyakov-loop model is employed in the heavy-quark limit precisely because it realizes a first-order deconfinement transition at finite baryon density with the Polyakov loop serving as the order parameter, permitting an explicit mapping between net-baryon cumulants and Polyakov-loop susceptibilities. The mean-field scaling relations themselves follow from the general Landau theory and are independent of the microscopic details. We have revised the introduction and conclusions to state explicitly that the results apply to the heavy-quark limit and to discuss the possible implications for full QCD, noting that while the universality class may change with light quarks, the qualitative structure of the cumulant scaling near the critical point is expected to persist. revision: partial

Circularity Check

0 steps flagged

No circularity: standard derivations and imported scaling functions

full rationale

The paper derives mean-field critical exponents directly from Landau theory (a standard textbook approach independent of the model), validates them via explicit numerical solution of the effective Polyakov-loop model, applies the conventional Ginzburg-Landau criterion as a diagnostic, and imports the known 3D Ising scaling function rather than fitting it to the same data. The mapping of net-baryon cumulants to Polyakov-loop susceptibilities is a model definition, not a self-referential prediction. No load-bearing step reduces by construction to a fitted parameter, self-citation chain, or ansatz smuggled from prior work by the same authors. The derivation chain remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Only the abstract is available, so the ledger is limited to assumptions explicitly invoked there; no free parameters or invented entities are named in the provided text.

axioms (2)
  • domain assumption Mean-field (Landau) theory applies near the deconfinement critical point
    Used to derive critical exponents for the cumulants
  • domain assumption The effective Polyakov-loop model captures essential critical physics in the heavy-quark limit
    Basis for numerical validation of the mean-field exponents

pith-pipeline@v0.9.0 · 5436 in / 1539 out tokens · 51395 ms · 2026-05-12T01:23:16.713767+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

61 extracted references · 61 canonical work pages

  1. [1]

    The last equation can be obtained from the requirement that at the critical point, the inverse of the curvature matrix of the original poten- tial becomes singular, i.e

    Two additional constraints arise from the gap equa- tions (7) (which setsh= 0). The last equation can be obtained from the requirement that at the critical point, the inverse of the curvature matrix of the original poten- tial becomes singular, i.e. det(C) = ∂2(U/T 4) ∂ℓ2 ∂2(U/T 4) ∂¯ℓ2 − ∂2(U/T 4) ∂ℓ∂ ¯ℓ 2 = 0. (B7) Using these equations,T CP ,M CP ,ℓ CP...

  2. [2]

    Properties of hot and dense matter from relativistic heavy ion collisions,

    P. Braun-Munzinger, V. Koch, T. Sch¨ afer and J. Stachel, “Properties of hot and dense matter from relativistic heavy ion collisions,” Phys. Rept. 621, 76-126 (2016) doi:10.1016/j.physrep.2015.12.003 [arXiv:1510.00442 [nucl-th]]

  3. [3]

    Search for the QCD Critical Point with Fluctuations of Conserved Quantities in Relativistic Heavy-Ion Collisions at RHIC : An Overview,

    X. Luo and N. Xu, “Search for the QCD Critical Point with Fluctuations of Conserved Quantities in Relativistic Heavy-Ion Collisions at RHIC : An Overview,” Nucl. Sci. Tech.28, no.8, 112 (2017) doi:10.1007/s41365-017-0257- 0 [arXiv:1701.02105 [nucl-ex]]

  4. [4]

    Decoding the phase structure of QCD via particle production at high energy,

    A. Andronic, P. Braun-Munzinger, K. Redlich and J. Stachel, “Decoding the phase structure of QCD via particle production at high energy,” Nature561, no.7723, 321-330 (2018) doi:10.1038/s41586-018-0491-6 [arXiv:1710.09425 [nucl-th]]

  5. [5]

    Mapping the Phases of Quantum Chro- modynamics with Beam Energy Scan,

    A. Bzdak, S. Esumi, V. Koch, J. Liao, M. Stephanov and N. Xu, “Mapping the Phases of Quantum Chro- modynamics with Beam Energy Scan,” Phys. Rept. 853, 1-87 (2020) doi:10.1016/j.physrep.2020.01.005 [arXiv:1906.00936 [nucl-th]]

  6. [6]

    The phase diagram of dense QCD,

    K. Fukushima and T. Hatsuda, “The phase diagram of dense QCD,” Rept. Prog. Phys.74, 014001 (2011) 13 doi:10.1088/0034-4885/74/1/014001 [arXiv:1005.4814 [hep-ph]]

  7. [7]

    The phase diagram of nuclear and quark matter at high baryon den- sity,

    K. Fukushima and C. Sasaki, “The phase diagram of nuclear and quark matter at high baryon den- sity,” Prog. Part. Nucl. Phys.72, 99-154 (2013) doi:10.1016/j.ppnp.2013.05.003 [arXiv:1301.6377 [hep- ph]]

  8. [8]

    Overview of the QCD phase dia- gram: Recent progress from the lattice,

    J. N. Guenther, “Overview of the QCD phase dia- gram: Recent progress from the lattice,” Eur. Phys. J. A 57, no.4, 136 (2021) doi:10.1140/epja/s10050-021-00354- 6 [arXiv:2010.15503 [hep-lat]]

  9. [9]

    Equation of state for QCD from lattice sim- ulations,

    C. Ratti, “Equation of state for QCD from lattice sim- ulations,” Prog. Part. Nucl. Phys.129, 104007 (2023) doi:10.1016/j.ppnp.2022.104007

  10. [10]

    Phase Transitions in Particle Physics: Results and Perspectives from Lattice Quantum Chromo-Dynamics,

    G. Aarts, J. Aichelin, C. Allton, A. Athenodorou, D. Bachtis, C. Bonanno, N. Brambilla, E. Bratkovskaya, M. Bruno and M. Caselle,et al.“Phase Transitions in Particle Physics: Results and Perspectives from Lattice Quantum Chromo-Dynamics,” Prog. Part. Nucl. Phys. 133, 104070 (2023) doi:10.1016/j.ppnp.2023.104070 [arXiv:2301.04382 [hep-lat]]

  11. [11]

    Phase diagram of QCD in a magnetic field: A review,

    J. O. Andersen, W. R. Naylor and A. Tranberg, “Phase diagram of QCD in a magnetic field: A review,” Rev. Mod. Phys.88, 025001 (2016) doi:10.1103/RevModPhys.88.025001 [arXiv:1411.7176 [hep-ph]]

  12. [12]

    Quantum field theory in a magnetic field: From quantum chromody- namics to graphene and Dirac semimetals,

    V. A. Miransky and I. A. Shovkovy, “Quantum field theory in a magnetic field: From quantum chromody- namics to graphene and Dirac semimetals,” Phys. Rept. 576, 1-209 (2015) doi:10.1016/j.physrep.2015.02.003 [arXiv:1503.00732 [hep-ph]]

  13. [13]

    QCD with background electromagnetic fields on the lattice: A review,

    G. Endrodi, “QCD with background electromagnetic fields on the lattice: A review,” Prog. Part. Nucl. Phys. 141, 104153 (2025) doi:10.1016/j.ppnp.2024.104153 [arXiv:2406.19780 [hep-lat]]

  14. [14]

    Adhikari, M

    P. Adhikari, M. Ammon, S. S. Avancini, A. Ay- ala, A. Bandyopadhyay, D. Blaschke, F. L. Braghin, P. Buividovich, R. P. Cardoso and C. Cartwright,et al. “Strongly interacting matter in extreme magnetic fields,” [arXiv:2412.18632 [nucl-th]]

  15. [15]

    Thermody- namics of SU(3) lattice gauge theory,

    G. Boyd, J. Engels, F. Karsch, E. Laermann, C. Leg- eland, M. Lutgemeier and B. Petersson, “Thermody- namics of SU(3) lattice gauge theory,” Nucl. Phys. B 469, 419-444 (1996) doi:10.1016/0550-3213(96)00170-8 [arXiv:hep-lat/9602007 [hep-lat]]

  16. [16]

    Critical Behavior at Finite Temperature Confinement Transitions,

    B. Svetitsky and L. G. Yaffe, “Critical Behavior at Finite Temperature Confinement Transitions,” Nucl. Phys. B 210, 423-447 (1982) doi:10.1016/0550-3213(82)90172-9

  17. [17]

    The Confinement problem in lat- tice gauge theory,

    J. Greensite, “The Confinement problem in lat- tice gauge theory,” Prog. Part. Nucl. Phys.51, 1 (2003) doi:10.1016/S0146-6410(03)90012-3 [arXiv:hep- lat/0301023 [hep-lat]]

  18. [18]

    Thermal Properties of Gauge Fields and Quark Liberation,

    A. M. Polyakov, “Thermal Properties of Gauge Fields and Quark Liberation,” Phys. Lett. B72, 477-480 (1978) doi:10.1016/0370-2693(78)90737-2

  19. [19]

    On the Phase Transition Towards Per- manent Quark Confinement,

    G. ’t Hooft, “On the Phase Transition Towards Per- manent Quark Confinement,” Nucl. Phys. B138, 1-25 (1978) doi:10.1016/0550-3213(78)90153-0

  20. [20]

    Mean Field Analysis of SU(N) Deconfining Transitions in the Presence of Dy- namical Quarks,

    F. Green and F. Karsch, “Mean Field Analysis of SU(N) Deconfining Transitions in the Presence of Dy- namical Quarks,” Nucl. Phys. B238, 297-306 (1984) doi:10.1016/0550-3213(84)90452-8

  21. [21]

    End point of the first-order phase transition of QCD in the heavy quark region by reweighting from quenched QCD,

    S. Ejiriet al.[WHOT-QCD], “End point of the first-order phase transition of QCD in the heavy quark region by reweighting from quenched QCD,” Phys. Rev. D101, no.5, 054505 (2020) doi:10.1103/PhysRevD.101.054505 [arXiv:1912.10500 [hep-lat]]

  22. [22]

    Deconfinement critical point of lattice QCD withN f=2 Wilson fermions,

    F. Cuteri, O. Philipsen, A. Sch¨ on and A. Scia- rra, “Deconfinement critical point of lattice QCD withN f=2 Wilson fermions,” Phys. Rev. D103, no.1, 014513 (2021) doi:10.1103/PhysRevD.103.014513 [arXiv:2009.14033 [hep-lat]]

  23. [23]

    High-precision analysis of the critical point in heavy-quark QCD at Nt=6,

    R. Ashikawa, M. Kitazawa, S. Ejiri and K. Kanaya, “High-precision analysis of the critical point in heavy-quark QCD at Nt=6,” Phys. Rev. D110, no.7, 074508 (2024) doi:10.1103/PhysRevD.110.074508 [arXiv:2407.09156 [hep-lat]]

  24. [24]

    Phase structure of QCD for heavy quarks,

    C. S. Fischer, J. Luecker and J. M. Pawlowski, “Phase structure of QCD for heavy quarks,” Phys. Rev. D91, no.1, 014024 (2015) doi:10.1103/PhysRevD.91.014024 [arXiv:1409.8462 [hep-ph]]

  25. [25]

    Crit- ical endpoint for deconfinement in matrix and other effective models,

    K. Kashiwa, R. D. Pisarski and V. V. Skokov, “Crit- ical endpoint for deconfinement in matrix and other effective models,” Phys. Rev. D85, 114029 (2012) doi:10.1103/PhysRevD.85.114029 [arXiv:1205.0545 [hep- ph]]

  26. [26]

    Polyakov loop fluctuations and deconfinement in the limit of heavy quarks,

    P. M. Lo, B. Friman and K. Redlich, “Polyakov loop fluctuations and deconfinement in the limit of heavy quarks,” Phys. Rev. D90, no.7, 074035 (2014) doi:10.1103/PhysRevD.90.074035 [arXiv:1406.4050 [hep- ph]]

  27. [27]

    Deconfinement in the presence of a strong mag- netic field,

    P. M. Lo, M. Szyma´ nski, C. Sasaki and K. Redlich, “Deconfinement in the presence of a strong mag- netic field,” Phys. Rev. D102, no.3, 034024 (2020) doi:10.1103/PhysRevD.102.034024 [arXiv:2004.04138 [hep-ph]]

  28. [28]

    De- confinement critical point of heavy quark effective lattice theories,

    A. Q. Pham, J. Kim and O. Philipsen, “De- confinement critical point of heavy quark effective lattice theories,” PoSLA TTICE2021, 113 (2022) doi:10.22323/1.396.0113 [arXiv:2112.06505 [hep-lat]]

  29. [29]

    Heavy-quark corner of the Columbia plot from the center-symmetric Curci- Ferrari model,

    V. T. Mari Surkau and U. Reinosa, “Heavy-quark corner of the Columbia plot from the center-symmetric Curci- Ferrari model,” Phys. Rev. D112, no.1, 014002 (2025) doi:10.1103/7dx2-vmj6 [arXiv:2503.22373 [hep-ph]]

  30. [30]

    Konrad and O

    C. Konrad and O. Philipsen, “From deconfinement to nuclear matter: mean-field approaches for ef- fective Polyakov loop theories of lattice QCD,” [arXiv:2509.03172 [hep-lat]]

  31. [31]

    Polyakov loop in 2+1 flavor QCD from low to high temperatures,

    A. Bazavov, N. Brambilla, H. T. Ding, P. Pe- treczky, H. P. Schadler, A. Vairo and J. H. We- ber, “Polyakov loop in 2+1 flavor QCD from low to high temperatures,” Phys. Rev. D93, no.11, 114502 (2016) doi:10.1103/PhysRevD.93.114502 [arXiv:1603.06637 [hep-lat]]

  32. [32]

    Lattice calculation of the Polyakov loop and Polyakov loop cor- relators,

    J. H. Weber [TUMQCD], “Lattice calculation of the Polyakov loop and Polyakov loop cor- relators,” EPJ Web Conf.137, 07028 (2017) doi:10.1051/epjconf/201713707028

  33. [33]

    Freeze-out parame- ters: lattice meets experiment,

    S. Borsanyi, Z. Fodor, S. D. Katz, S. Krieg, C. Ratti and K. K. Szabo, “Freeze-out parame- ters: lattice meets experiment,” Phys. Rev. Lett. 111, 062005 (2013) doi:10.1103/PhysRevLett.111.062005 [arXiv:1305.5161 [hep-lat]]

  34. [34]

    Strangeness at high temperatures: from hadrons to quarks,

    A. Bazavov, H. T. Ding, P. Hegde, O. Kaczmarek, F. Karsch, E. Laermann, Y. Maezawa, S. Mukherjee, H. Ohno and P. Petreczky,et al.“Strangeness at high temperatures: from hadrons to quarks,” Phys. Rev. Lett. 111, 082301 (2013) doi:10.1103/PhysRevLett.111.082301 14 [arXiv:1304.7220 [hep-lat]]

  35. [35]

    Skewness, kur- tosis, and the fifth and sixth order cumulants of net baryon-number distributions from lattice QCD con- front high-statistics STAR data,

    A. Bazavov, D. Bollweg, H. T. Ding, P. Enns, J. Goswami, P. Hegde, O. Kaczmarek, F. Karsch, R. Larsen and S. Mukherjee,et al.“Skewness, kur- tosis, and the fifth and sixth order cumulants of net baryon-number distributions from lattice QCD con- front high-statistics STAR data,” Phys. Rev. D101, no.7, 074502 (2020) doi:10.1103/PhysRevD.101.074502 [arXiv:20...

  36. [36]

    Continuum-extrapolated high- order baryon fluctuations,

    S. Borsanyi, Z. Fodor, J. N. Guenther, S. D. Katz, P. Parotto, A. Pasztor, D. Pesznyak, K. K. Szabo and C. H. Wong, “Continuum-extrapolated high- order baryon fluctuations,” Phys. Rev. D110, no.1, L011501 (2024) doi:10.1103/PhysRevD.110.L011501 [arXiv:2312.07528 [hep-lat]]

  37. [37]

    Fluctuations of con- served charges in relativistic heavy ion collisions: An in- troduction,

    M. Asakawa and M. Kitazawa, “Fluctuations of con- served charges in relativistic heavy ion collisions: An in- troduction,” Prog. Part. Nucl. Phys.90, 299-342 (2016) doi:10.1016/j.ppnp.2016.04.002 [arXiv:1512.05038 [nucl- th]]

  38. [38]

    Thermodynamics of fermionic excitations in heavy-quark QCD,

    K. Tohme, T. M. Doi, M. Kitazawa, K. Redlich and C. Sasaki, “Thermodynamics of fermionic excitations in heavy-quark QCD,” [arXiv:2508.09927 [hep-lat]]

  39. [39]

    Non-Gaussian fluctuations near the QCD critical point,

    M. A. Stephanov, “Non-Gaussian fluctuations near the QCD critical point,” Phys. Rev. Lett.102, 032301 (2009) doi:10.1103/PhysRevLett.102.032301 [arXiv:0809.3450 [hep-ph]]

  40. [40]

    Susceptibilities near the QCD (tri)critical point,

    B. J. Schaefer and J. Wambach, “Susceptibilities near the QCD (tri)critical point,” Phys. Rev. D 75, 085015 (2007) doi:10.1103/PhysRevD.75.085015 [arXiv:hep-ph/0603256 [hep-ph]]

  41. [41]

    QCD critical region and higher moments for three flavor models,

    B. J. Schaefer and M. Wagner, “QCD critical region and higher moments for three flavor models,” Phys. Rev. D85, 034027 (2012) doi:10.1103/PhysRevD.85.034027 [arXiv:1111.6871 [hep-ph]]

  42. [42]

    Lectures on phase transitions and the renormalization group,

    N. Goldenfeld, “Lectures on phase transitions and the renormalization group,”

  43. [43]

    Introduction to the functional renormalization group,

    P. Kopietz, L. Bartosch and F. Sch¨ utz, “Introduction to the functional renormalization group,” Lect. Notes Phys. 798, 1-380 (2010) doi:10.1007/978-3-642-05094-7

  44. [44]

    Universality, the QCD critical / tricritical point and the quark num- ber susceptibility,

    Y. Hatta and T. Ikeda, “Universality, the QCD critical / tricritical point and the quark num- ber susceptibility,” Phys. Rev. D67, 014028 (2003) doi:10.1103/PhysRevD.67.014028 [arXiv:hep- ph/0210284 [hep-ph]]

  45. [45]

    Beyond mean field: Fluctuation diagnos- tics and fixed-point behavior,

    P. M. Lo, “Beyond mean field: Fluctuation diagnos- tics and fixed-point behavior,” J. Subatomic Part. Cos- mol.5, 100379 (2026) doi:10.1016/j.jspc.2026.100379 [arXiv:2602.21095 [nucl-th]]

  46. [46]

    Quark Libera- tion at High Temperature: A Monte Carlo Study of SU(2) Gauge Theory,

    L. D. McLerran and B. Svetitsky, “Quark Libera- tion at High Temperature: A Monte Carlo Study of SU(2) Gauge Theory,” Phys. Rev. D24, 450 (1981) doi:10.1103/PhysRevD.24.450

  47. [47]

    Polyakov loop fluctuations in SU(3) lattice gauge theory and an effective gluon potential,

    P. M. Lo, B. Friman, O. Kaczmarek, K. Redlich and C. Sasaki, “Polyakov loop fluctuations in SU(3) lattice gauge theory and an effective gluon potential,” Phys. Rev. D88, 074502 (2013) doi:10.1103/PhysRevD.88.074502 [arXiv:1307.5958 [hep-lat]]

  48. [48]

    Polyakov loop model- ing for hot QCD,

    K. Fukushima and V. Skokov, “Polyakov loop model- ing for hot QCD,” Prog. Part. Nucl. Phys.96, 154-199 (2017) doi:10.1016/j.ppnp.2017.05.002 [arXiv:1705.00718 [hep-ph]]

  49. [49]

    Dynamical Model of Elementary Particles Based on an Analogy with Su- perconductivity. 1.,

    Y. Nambu and G. Jona-Lasinio, “Dynamical Model of Elementary Particles Based on an Analogy with Su- perconductivity. 1.,” Phys. Rev.122, 345-358 (1961) doi:10.1103/PhysRev.122.345

  50. [50]

    Dynamical model of elementary particles based on an analogy with su- perconductivity. II.,

    Y. Nambu and G. Jona-Lasinio, “Dynamical model of elementary particles based on an analogy with su- perconductivity. II.,” Phys. Rev.124, 246-254 (1961) doi:10.1103/PhysRev.124.246

  51. [51]

    The Nambu-Jona-Lasinio model of quantum chromodynamics,

    S. P. Klevansky, “The Nambu-Jona-Lasinio model of quantum chromodynamics,” Rev. Mod. Phys.64, 649- 708 (1992) doi:10.1103/RevModPhys.64.649

  52. [52]

    NJL model analysis of quark matter at large density,

    M. Buballa, “NJL model analysis of quark matter at large density,” Phys. Rept.407, 205-376 (2005) doi:10.1016/j.physrep.2004.11.004

  53. [53]

    Interplay between chiral and axial symmetries in a SU(2) Nambu-Jona-Lasinio Model with the Polyakov loop,

    M. C. Ruivo, M. Santos, P. Costa and C. A. de Sousa, “Interplay between chiral and axial symmetries in a SU(2) Nambu-Jona-Lasinio Model with the Polyakov loop,” Phys. Rev. D85, 036001 (2012)

  54. [54]

    Causality violation and the speed of sound of hot and dense quark matter in the Nambu–Jona-Lasinio model,

    A. E. B. Pasqualotto, R. L. S. Farias, W. R. Tavares, S. S. Avancini and G. Krein, “Causality violation and the speed of sound of hot and dense quark matter in the Nambu–Jona-Lasinio model,” Phys. Rev. D107, no.9, 9 (2023)

  55. [55]

    Susceptibilities and the Phase Structure of a Chiral Model with Polyakov Loops,

    C. Sasaki, B. Friman and K. Redlich, “Susceptibilities and the Phase Structure of a Chiral Model with Polyakov Loops,” Phys. Rev. D75, 074013 (2007)

  56. [56]

    Deconfinement of heavy quarks at finite density and strong magnetic field,

    M. Szyma´ nski, “Deconfinement of heavy quarks at finite density and strong magnetic field,” Eur. Phys. J. ST229, no.22-23, 3387-3394 (2020) doi:10.1140/epjst/e2020- 000115-3

  57. [57]

    Hadronic fluctu- ations at the QCD phase transition,

    S. Ejiri, F. Karsch and K. Redlich, “Hadronic fluctu- ations at the QCD phase transition,” Phys. Lett. B 633, 275-282 (2006) doi:10.1016/j.physletb.2005.11.083 [arXiv:hep-ph/0509051 [hep-ph]]

  58. [58]

    Scaling functions of the three-dimensional Z(2), O(2), and O(4) models and their finite-size dependence in an external field,

    F. Karsch, M. Neumann and M. Sarkar, “Scaling functions of the three-dimensional Z(2), O(2), and O(4) models and their finite-size dependence in an external field,” Phys. Rev. D108, no.1, 014505 (2023) doi:10.1103/PhysRevD.108.014505 [arXiv:2304.01710 [hep-lat]]

  59. [59]

    Bootstrapping the 3d Ising stress tensor,

    C. H. Chang, V. Dommes, R. S. Erramilli, A. Homrich, P. Kravchuk, A. Liu, M. S. Mitchell, D. Poland and D. Simmons-Duffin, “Bootstrapping the 3d Ising stress tensor,” JHEP03, 136 (2025) doi:10.1007/JHEP03(2025)136 [arXiv:2411.15300 [hep- th]]

  60. [60]

    Buryak, P.D

    A. Dumitru and R. D. Pisarski, “Event-by-event fluc- tuations from decay of a Polyakov loop condensate,” Phys. Lett. B504, 282-290 (2001) doi:10.1016/S0370- 2693(01)00286-6 [arXiv:hep-ph/0010083 [hep-ph]]

  61. [61]

    Confinement and Chiral Phase Tran- sitions: The Role of Polyakov Loop Kinetics Terms,

    B. Hua and J. Zhu, “Confinement and Chiral Phase Tran- sitions: The Role of Polyakov Loop Kinetics Terms,” [arXiv:2601.19270 [hep-ph]]