pith. sign in

arxiv: 1907.08056 · v1 · pith:4DLSSWEUnew · submitted 2019-07-18 · ⚛️ nucl-th

Searching for the QCD critical point with net-proton number fluctuations

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

classification ⚛️ nucl-th
keywords net-proton fluctuationsQCD critical pointcumulant ratiosbeam energy dependenceheavy ion collisionsSTAR experimentcritical fluctuationsphenomenological model
0
0 comments X

The pith

A phenomenological model coupling critical fluctuations to protons reproduces the beam-energy dependence of net-proton cumulant ratios seen in STAR data.

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

The paper examines net-proton number fluctuations measured in heavy-ion collisions as a potential probe of the QCD critical point. It develops a model in which critical mode fluctuations couple to protons and anti-protons through a phenomenological interaction. Ratios of the first four cumulants are computed within this framework as a function of beam energy. The model accounts for the monotonic trend in the lowest-order ratio and the non-monotonic trends in higher-order ratios that appear in the experimental measurements. Results vary with the strength of the coupling and the assumed location of the critical point, offering a direct link between theory parameters and observable trends.

Core claim

The central claim is that a phenomenologically motivated model in which critical mode fluctuations couple to protons and anti-protons qualitatively captures both the monotonic behavior of the lowest-order cumulant ratio and the non-monotonic behavior of higher-order ratios as a function of beam energy, matching the trends reported by the STAR Collaboration.

What carries the argument

Phenomenological coupling of critical mode fluctuations to protons and anti-protons, with adjustable coupling strength and critical-point location.

If this is right

  • Higher-order cumulant ratios serve as indicators of critical-point effects through their non-monotonic energy dependence.
  • Fitting the model parameters to data constrains possible locations of the QCD critical point.
  • The approach connects theoretical critical fluctuations to measurable net-proton observables in heavy-ion collisions.
  • Varying coupling strength shows how interaction details influence the visibility of critical signatures.

Where Pith is reading between the lines

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

  • The same coupling framework could be applied to other conserved charges such as net-charge or net-baryon fluctuations for consistency checks.
  • Targeted measurements at beam energies near the fitted critical point would provide a direct test of the predicted non-monotonic patterns.
  • If the model holds, it supplies a practical method to translate fluctuation data into constraints on critical-point properties.

Load-bearing premise

Critical mode fluctuations can be coupled to protons and anti-protons through a phenomenological model whose parameters are varied to match observed trends.

What would settle it

Absence of the predicted non-monotonic behavior in higher-order cumulant ratios at beam energies corresponding to the fitted critical-point location would falsify the model's ability to explain the data.

Figures

Figures reproduced from arXiv: 1907.08056 by Chihiro Sasaki, Krzysztof Redlich, Marcus Bluhm, Micha{\l} Szyma\'nski.

Figure 1
Figure 1. Figure 1: (Color online) The model setup used in this work. The filled band between the two dashed curves shows lattice QCD constraints for the chiral crossover transition. The green dot denotes the critical point with the spin model coordinate system attached to it and the first-order phase transition line for larger baryon chemical potential. The solid blue line corre￾sponds to the chemical freeze-out curve from [… view at source ↗
Figure 2
Figure 2. Figure 2: (Color online) QCD critical point locations from Tab. 1 plotted with the chemical freeze-out curve [15] used in this work. belongs to the same universality class as the three￾dimensional Ising model [16–18] we can identify the QCD order parameter, 𝜎, with the magnetization, 𝑀𝐼 , the order parameter of the spin model. Hence, the critical mode cumulants can be written as [10] ⟨(𝑉 𝛿𝜎) 𝑛 ⟩𝑐 = (︂ 𝑇 𝑉 𝐻0 )︂𝑛−1 ∂… view at source ↗
Figure 3
Figure 3. Figure 3: (Color online) The second to first net-proton number cumulant ratio for 𝑔 = 3 and 5 calculated following Ref. [10] (red solid and dashed lines, respectively) compared to refined model results [11] (blue solid and dashed lines, respectively). The preliminary STAR data on the net-proton number fluctua￾tions [5] (squares with the error bars containing both statistical and systematic errors) and HRG baseline r… view at source ↗
Figure 4
Figure 4. Figure 4: (Color online) Ratios of net-proton number cumulants calculated in the refined model [11] for fixed coupling 𝑔 = 5 and for different locations of the QCD critical point (listed in Tab. 1). gram on the refined model results, we consider three different locations of the CP listed in Tab. 1 and shown in [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (Color online) Ratios of net-proton number cumulants calculated in the refined model [11] with CP3 and for coupling strengths, 𝑔 = 3, 4 and 5 (orange solid, green long-dashed and red dash-dotted lines, respectively). The preliminary STAR data on the net-proton number fluctuations [5] (squares with the error bars containing both statistical and systematic errors) and HRG baseline results (black dotted lines… view at source ↗
Figure 3
Figure 3. Figure 3: Results obtained using the original model [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
read the original abstract

Net-proton number fluctuations can be measured experimentally and hence provide a source of important information about the matter created during relativistic heavy ion collisions. Particularly, they may give us clues about the conjectured QCD critical point. In this work the beam-energy dependence of ratios of the first four cumulants of the net-proton number is discussed. These quantities are calculated using a phenomenologically motivated model in which critical mode fluctuations couple to protons and anti-protons. Our model qualitatively captures both the monotonic behavior of the lowest-order ratio as well as the non-monotonic behavior of higher-order ratios, as seen in the experimental data from the STAR Collaboration. We also discuss the dependence of our results on the coupling strength and the location of the critical point.

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 presents a phenomenologically motivated model in which critical-mode fluctuations are coupled to protons and anti-protons. It computes the beam-energy dependence of ratios of the first four cumulants of the net-proton number and states that the model qualitatively reproduces both the monotonic behavior of the lowest-order ratio and the non-monotonic behavior of higher-order ratios observed in STAR data. The results are shown to depend on the coupling strength and the location of the critical point.

Significance. If the coupling were derived from first principles and the non-monotonicity emerged without parameter adjustment to the same data, the approach could provide a useful framework for interpreting cumulant ratios as potential signatures of the QCD critical point. As presented, the qualitative agreement is obtained by varying free parameters (coupling strength and critical-point coordinates) to match the observed trends, limiting the extent to which the calculation constitutes an independent test of the critical-point hypothesis.

major comments (1)
  1. [Abstract] Abstract: the central claim that the model 'qualitatively captures' the non-monotonic behavior of higher-order ratios rests on explicit variation of the coupling strength and the location of the critical point 'to discuss agreement with data.' Because these parameters are adjusted until the trends match the STAR measurements being explained, the non-monotonicity is at least partly engineered by choice of location rather than emerging as a robust consequence of critical physics; the lowest-order ratio is automatically monotonic once the coupling is non-zero.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the constructive feedback. We respond to the major comment below.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the central claim that the model 'qualitatively captures' the non-monotonic behavior of higher-order ratios rests on explicit variation of the coupling strength and the location of the critical point 'to discuss agreement with data.' Because these parameters are adjusted until the trends match the STAR measurements being explained, the non-monotonicity is at least partly engineered by choice of location rather than emerging as a robust consequence of critical physics; the lowest-order ratio is automatically monotonic once the coupling is non-zero.

    Authors: We agree that the model is phenomenological and that the coupling strength and critical-point location are varied to explore the conditions under which the observed trends in the STAR data can be reproduced. The non-monotonicity arises specifically from the coupling of protons to the critical mode when the critical point lies in a region traversed by the beam-energy scan; it is a dynamical consequence of the model rather than an arbitrary imposition. Nevertheless, we acknowledge that the agreement is obtained within a scanned parameter space and therefore does not constitute a parameter-free prediction or an independent test of the critical-point hypothesis. We will revise the abstract to state explicitly that the qualitative reproduction is demonstrated by exploring the model's parameter dependence rather than claimed as a robust, untuned outcome. The monotonicity of the lowest-order ratio is indeed a direct structural feature once the coupling is non-zero, while its detailed beam-energy dependence remains sensitive to the critical-point coordinates. revision: partial

Circularity Check

0 steps flagged

No circularity: phenomenological model with explicit parameter variation

full rationale

The paper presents calculations within a phenomenologically motivated model whose coupling strength and critical-point location are openly varied to explore agreement with STAR data. The central claim is qualitative capture of monotonic and non-monotonic trends, not a parameter-free prediction or derivation. No quoted step reduces by construction to its inputs, no self-citation is load-bearing for a uniqueness claim, and the model is self-contained as an exploratory tool rather than a closed loop. This is the expected honest non-finding for a tuned phenomenological study.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The model rests on a domain assumption that critical fluctuations can be coupled phenomenologically to protons; two free parameters (coupling strength and critical-point location) are introduced and varied to match data trends.

free parameters (2)
  • coupling strength
    Parameter controlling how strongly critical fluctuations affect proton production; varied to discuss dependence on results.
  • location of the critical point
    Coordinates in the phase diagram used to place the critical mode; varied to explore agreement with data.
axioms (1)
  • domain assumption Critical mode fluctuations couple to protons and anti-protons in a manner that can be modeled phenomenologically.
    Stated as the basis of the model in the abstract.

pith-pipeline@v0.9.0 · 5660 in / 1296 out tokens · 42261 ms · 2026-05-24T19:28:51.705664+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages · 1 internal anchor

  1. [1]

    Searching for the QCD critical point with net-proton number fluctuations

    Introduction The theoretical and experimental investigation of the phasediagramofstronglyinteractingmatterisanim- portant subject of modern high energy physics. One of the unresolved questions concerns the existence and location of the QCD critical point (CP) in the 𝑇 and 𝜇 plane. Strong fluctuations of the critical mode, 𝜎, in the vicinity of CP, althoug...

  2. [2]

    Model setup As a baseline model to calculate the net-proton num- ber cumulants we choose the hadron resonance gas (HRG) model in which the number density of each particle species is given by the ideal gas formula, 𝑛𝑖(𝑇, 𝜇𝑖) = 𝑑𝑖 ∫︁𝑑3𝑘 (2𝜋)3 𝑓0 𝑖 (𝑇, 𝜇𝑖) . (1) Here 𝑑𝑖 is the degeneracy factor and 𝑓0 𝑖 = 1 (−1)𝐵𝑖 + 𝑒(𝐸𝑖−𝜇𝑖)/𝑇 (2) is the equilibrium distribu...

  3. [3]

    Numerical results In this section we discuss numerical results on net- proton number cumulant ratios obtained within the current model. The set of model parameters includes thecouplingstrength 𝑔 between(anti)protonsandthe critical mode, the parameters of the magnetic equa- tion of state as well as the size of the critical region in the (𝑇, 𝜇) plane. Their...

  4. [4]

    Etoiles montantes en Pays de la Loire 2017

    Conclusions We presented ratios of net-proton number cumulants obtained within an effective model in which the cou- pling between (anti)protons and critical mode fluc- tuations is introduced by connecting particle masses to the order parameter. We modified the existing approach [10] to take into account the correct scal- ing properties of the baryon numbe...

  5. [5]

    M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. Lett. 81, 4816 (1998)

  6. [6]

    M. A. Stephanov, K. Rajagopal and E. V. Shuryak, Phys. Rev. D 60, 114028 (1999)

  7. [7]

    Luo [STAR Collaboration], PoS CPOD 2014, 019 (2015)

    X. Luo [STAR Collaboration], PoS CPOD 2014, 019 (2015)

  8. [8]

    Luo, Nucl

    X. Luo, Nucl. Phys. A956, 75 (2016)

  9. [9]

    Th¨ ader [STAR Collaboration], Nucl

    J. Th¨ ader [STAR Collaboration], Nucl. Phys. A956 320 (2016)

  10. [10]

    G. A. Almasi, B. Friman and K. Redlich, Phys. Rev. D96, no. 1, 014027 (2017)

  11. [11]

    Karsch, J

    F. Karsch, J. Phys. Conf. Ser.779, 012015 (2017)

  12. [12]

    Bzdak and V

    A. Bzdak and V. Koch, Phys. Rev. C99, no. 2, 024913 (2019)

  13. [13]

    Koch and A

    V. Koch and A. Bzdak, Acta Phys. Polon. B 47, 1867 (2016)

  14. [14]

    Bluhm, M

    M. Bluhm, M. Nahrgang, S. A. Bass and T. Sch¨ afer, Eur. Phys. J. C77, no. 4, 210 (2017)

  15. [15]

    Szyma´ nski, M

    M. Szyma´ nski, M. Bluhm, K. Redlich and C. Sasaki, arXiv:1905.00667 [nucl-th]

  16. [16]

    Sasaki, B

    C. Sasaki, B. Friman and K. Redlich, Phys. Rev. D75, 054026 (2007)

  17. [17]

    Sasaki, B

    C. Sasaki, B. Friman and K. Redlich, Phys. Rev. D77, 034024 (2008)

  18. [18]

    Zinn-Justin, Quantum Field Theory and Critical Phe- nomena, Clarendon Press, Oxford (2002)

    J. Zinn-Justin, Quantum Field Theory and Critical Phe- nomena, Clarendon Press, Oxford (2002)

  19. [19]

    Andronic, P

    A. Andronic, P. Braun-Munzinger, K. Redlich and J. Stachel, Nature561, no. 7723, 321 (2018)

  20. [20]

    Wilczek, Int

    F. Wilczek, Int. J. Mod. Phys. A07, 3911 (1992)

  21. [21]

    A. M. Halasz, A. D. Jackson, R. E. Shrock, M. A. Stephanov, and J. J. M. Verbaarschot, Phys. Rev. D58, 096007 (1998)

  22. [22]

    H.-T. Ding, F. Karsch, and S. Mukherjee, Int. J. Mod. Phys. E 24, 1530007 (2015)

  23. [23]

    Hatta and T

    Y. Hatta and T. Ikeda, Phys. Rev. D67, 014028 (2003)

  24. [24]

    Mukherjee, R

    S. Mukherjee, R. Venugopalan and Y. Yin, Phys. Rev. C 92, no. 3, 034912 (2015)

  25. [25]

    Nonaka and M

    C. Nonaka and M. Asakawa, Phys. Rev. C 71, 044904 (2005)

  26. [26]

    Guida and J

    R. Guida and J. Zinn-Justin, Nucl. Phys. B 489, 626 (1997)

  27. [27]

    Abelev et al

    B. Abelev et al. [ALICE Collaboration], Phys. Rev. C88, 044910 (2013)

  28. [28]

    B. B. Abelev et al. [ALICE Collaboration], Phys. Rev. Lett. 111, 222301 (2013)

  29. [29]

    B. B. Abelevet al. [ALICE Collaboration], Phys. Lett. B 728, 216 (2014), Erratum: [Phys. Lett. B734, 409 (2014)]

  30. [30]

    B. B. Abelev et al. [ALICE Collaboration], Phys. Rev. C 91, 024609 (2015)

  31. [31]

    Adamet al

    J. Adamet al. [ALICE Collaboration], Phys. Lett. B754, 360 (2016)

  32. [32]

    Adam et al

    J. Adam et al. [ALICE Collaboration], Phys. Rev. C93, no. 2, 024917 (2016). Received 18.07.19 6 ISSN 2071-0186. Ukr. J. Phys. ZZZZ. Vol. YY, No. XX