Deconfinement and freezeout boundaries in equilibrium thermal models

Abdel Nasser Tawfik (Nile U.; A. H. El-Kateb; ECTP; Johann Wolfgang Goethe-Universitat); Muhammad Maher; Sara Abdelaziz (Helwan U)

Deconfinement and freezeout boundaries in equilibrium thermal models

T0 review · 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 Equilibrium thermal model curves for chemical freezeout and deconfinement reproduce heavy-ion data, coincide with lattice QCD at low baryon density, and separate at high baryon density, suggesting a mixed hadron-QGP window.

desk verdict A workmanlike phenomenological map of freezeout and deconfinement boundaries whose low-density claims hold up, but whose new high-density mixed-phase window is a gap between two unquantified fitted curves. read the letter →

arxiv 1908.00426 v2 pith:FRAZFQFX submitted 2019-08-01 hep-ph

Abdel Nasser Tawfik (Nile U. , ECTP , Johann Wolfgang Goethe-Universitat) , Muhammad Maher , A. H. El-Kateb , Sara Abdelaziz (Helwan U) This is my paper · ORCID

classification hep-ph

keywords boundariesfreezeoutmodelthermaldeconfinementdensityagreementbaryonic

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

Heavy-ion experiments smash nuclei together to create a short-lived fireball of quarks and gluons, which expands and cools into ordinary hadrons. Two boundaries matter in such collisions. The chemical freezeout boundary is the moment when the types and numbers of particles stop changing. The deconfinement boundary is where quarks and gluons stop being bound inside hadrons and move freely. This paper uses an equilibrium thermal model, an ideal gas of known hadrons, to draw both boundaries on a chart with temperature on one axis and baryon density on the other.

The freezeout boundary from the model matches the experimental points across the whole chart. At low baryon density, the freezeout and deconfinement boundaries sit almost on top of each other, and both agree with lattice QCD simulations. At high baryon density, the two curves separate. The paper interprets the space between them as a region where hadrons and quark-gluon plasma may coexist, and it gives a rough window for this region, baryon chemical potential between about 320 and 560 MeV.

The limitations are that the constants that define the boundaries come from earlier fits, the curves have no error bars, and the comparison with the Polyakov linear-sigma model is partly circular because that model was evaluated using the same freezeout condition. This is a useful comparison study, not a derivation from first principles.

Extended reading notes

Core claim

The central claim is that at low baryon density the deconfinement and chemical freezeout boundaries coincide and agree with lattice QCD, while at large baryon density they separate, creating a window (mu_b roughly 320 to 560 MeV) where hadrons and quark-gluon plasma likely coexist. The paper states: 'Along the entire freezeout boundary, there is an excellent agreement between the thermal model calculations and the experiments.'

Load-bearing premise

The comparison rests on the assumption that chemical freezeout occurs at a universal constant value of entropy density over T^3 and that deconfinement occurs at a universal constant energy density, with the numeric constants inherited from earlier fits (refs [22,23,28]); the values are not restated here, and the mixed-phase window is read from the gap between these two fitted-condition curves.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

Quick take: this is a workmanlike phenomenological paper that does a useful job of assembling the T-μ_b plane comparison, but the one genuinely new claim—a mixed-phase window at high baryon density—is read off the gap between two fitted curves whose defining constants are never stated and whose uncertainties are not shown. The low-density claim, that deconfinement and freezeout coincide and agree with lattice, is much better supported.

What's actually new: the paper brings together experimental freezeout points, lattice bands, PLSM results, and HRG curves under two constant-condition hypotheses, and it explicitly highlights the separation at large μ_b. That is a reasonable map for the FAIR/NICA energy range. The HRG machinery is standard, the data compilation is useful, and the authors are honest about the PLSM freezeout condition and about the missing-states sensitivity. Credit where due: the comparison with experimental freezeout points across the whole μ_b range is genuine external validation, and the low-density coincidence with lattice is a real check.

The soft spots are real but not fatal. The constant s/T^3 and constant energy density values are inherited from earlier fits and not restated here, so the reader cannot reproduce Fig. 1 without going to the earlier papers. The curves have no uncertainties, so the separation between them at μ_b ~ 320–560 MeV is visually appealing but quantitatively unestablished. And the agreement with PLSM is partly built in: both sides use the same constant s/T^3 freezeout condition, so the agreement is not an independent confirmation. The paper's own caveat about missing states is the right worry, but it is left unquantified; since entropy density and energy density respond differently to additional resonances, the two boundaries could shift by different amounts and the window could close. I think the stress-test note is right that this is the load-bearing weakness.

That said, the central low-density argument holds up, and the paper does not contain a mathematical error. It deserves a serious referee, but I would want the constants stated, error bands on the curves, and a sensitivity study for missing states before the mixed-phase claim is taken as established. The paper is for phenomenologists working on the QCD phase diagram and heavy-ion freezeout; it is a reasonable reference map, not a breakthrough.

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

The central curves are not derived parameter-free: the freezeout and deconfinement conditions are predefined constants, and the mu_b mapping in Eq (4) uses earlier fits. No new entities are introduced. The PLSM agreement is partly circular because the same freezeout condition defines both curves.

free parameters (4)

constant s/T^3 for freezeout = not stated in this paper (from refs [22,23])
Defines the chemical freezeout curve in the HRG; previously fitted to reproduce particle yields.
constant rho/T^4 for deconfinement = not stated in this paper (from ref [28])
Defines the deconfinement curve; previously chosen as a line-of-constant-physics condition.
mu_b vs sqrt(s_NN) coefficients a, b = a = 1.245 +/- 0.049 GeV, b = 0.244 +/- 0.028 GeV^-1
Eq (4) used to map experimental energies to baryon chemical potential; fitted to data in ref [20].
hadron mass cutoff = 2.5 GeV
PDG hadrons up to 2.5 GeV included; a truncation choice that affects the thermodynamics.

assumptions (6)

standard math Grand canonical ideal gas partition function describes the hadron resonance gas
Eqs (1)-(3), standard statistical mechanics.
ad hoc to paper Freezeout occurs at constant s/T^3
Assumed universal freezeout condition, not derived from QCD; the constant is fit.
ad hoc to paper Deconfinement occurs at constant energy density
Line-of-constant-physics condition, not derived; the constant is fit.
domain assumption PDG hadron list up to 2.5 GeV captures relevant thermodynamics
Truncation of the resonance spectrum; missing states separately tested (refs [57-59]).
domain assumption Strangeness chemical potential fixed by net strangeness neutrality
Used to close the system of chemical potentials.
domain assumption Polyakov linear-sigma model is a valid effective theory for QCD thermodynamics
Used for comparison; from refs [63,65-70].

how reviews work

0 comments

Cite this review

Pith. "Pith review of Deconfinement and freezeout boundaries in equilibrium thermal models." pith.science (2026). https://pith.science/paper/FRAZFQFX

@misc{pith2026190800426,
  author       = {Pith},
  title        = {Pith review of: Deconfinement and freezeout boundaries in equilibrium thermal models},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FRAZFQFX}},
  note         = {Machine review of arXiv:1908.00426}
}

read the original abstract

In different approaches, the temperature-baryon density plane of QCD matter is studied for deconfinement and chemical freezeout boundaries. Results from various heavy-ion experiments are compared with the recent lattice simulations, the effective QCD-like Polyakov linear-sigma model, and the equilibrium thermal models. Along the entire freezeout boundary, there is an excellent agreement between the thermal model calculations and the experiments. Also, the thermal model calculations agree well with the estimations deduced from the Polyakov linear-sigma model (PLSM). At low baryonic density or high energies, both deconfinement and chemical freezeout boundaries are likely coincident and therefore the agreement with the lattice simulations becomes excellent as well, while at large baryonic density, the two boundaries become distinguishable forming a phase where hadrons and quark-gluon plasma likely coexist.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

71 extracted references · 62 canonical work pages

[1]

Banks and A

T. Banks and A. Ukawa, Nucl. Phys. B225, 145 (1983)

work page 1983
[2]

A. N. Tawﬁk, Int. J. Mod. Phys. A29, 1430021 (2014)

work page 2014
[3]

B. P. Abbott et al., Phys. Rev. Lett. 119, 161101 (2017)

2017
[4]

D. D. Ivanenko and D. F. Kurdgelaidze, Astrophysics 1, 251 (1965)

work page 1965
[5]

Gyulassy and L

M. Gyulassy and L. McLerran, Nucl. Phys. A750, 30 (2005)

2005
[6]

Busza, K

W. Busza, K. Rajagopal, and W. van der Schee, Ann. Rev. Nuc l. Part. Sci. 68, 339 (2018)

work page 2018
[7]

Greiner, P

C. Greiner, P. Koch-Steinheimer, F. M. Liu, I. A. Shovkov y, and H. Stoecker, J. Phys. G31, S725 (2005)

work page 2005
[8]

Why does the thermal model for hadron production in heavy ion collisions work?

B. Mueller and A. Schaefer, arXiv:1712.03567 [nucl-th] (2017)

work page Pith review arXiv 2017

Show all 71 references

[9]

Vovchenko, M

V. Vovchenko, M. I. Gorenstein, C. Greiner, and H. Stoeck er, Phys. Rev. C99, 045204 (2019). 7

2019
[10]

D’Alessio, Y

L. D’Alessio, Y. Kafri, A. Polkovnikov, and M. Rigol, Ad v. Phys. 65, 239 (2016)

2016
[11]

G. Fast, R. Hagedorn, and L. W. Jones, Nuovo Cim. 27, 856 (1963)

1963
[12]

Fast and R

G. Fast and R. Hagedorn, Nuovo Cim. 27, 208 (1963)

1963
[13]

Karsch, K

F. Karsch, K. Redlich, and A. Tawﬁk, Eur. Phys. J. C29, 549 (2003)

2003
[14]

Karsch, K

F. Karsch, K. Redlich, and A. Tawﬁk, Phys. Lett. B571, 67 (2003)

2003
[15]

Redlich, F

K. Redlich, F. Karsch, and A. Tawﬁk, J. Phys. G30, S1271 (2004)

2004
[16]

Tawﬁk, Phys

A. Tawﬁk, Phys. Rev. D71, 054502 (2005)

2005
[17]

Tawﬁk and D

A. Tawﬁk and D. Toublan, Phys. Lett. B623, 48 (2005)

2005
[18]

Tanabashi et al., Phys

M. Tanabashi et al., Phys. Rev. D98, 030001 (2018)

2018
[19]

The strangeness degrees of freedo m play an important role, especially at low µb or high √sNN

are taken into account. The strangeness degrees of freedo m play an important role, especially at low µb or high √sNN. For the sake of completeness, we have also checked the same calcu lations but with the inclusion of a large number of possible missing states [57, 58]. We fou...
[20]

Beringer et al., Phys

J. Beringer et al., Phys. Rev. D86, 010001 (2012)

2012
[21]

A. N. Tawﬁk and E. Abbas, Phys. Part. Nucl. Lett. 12, 521 (2015)

2015
[22]

Tawﬁk, M

A. Tawﬁk, M. Y. El-Bakry, D. M. Habashy, M. T. Mohamed, an d E. Abbas, Int. J. Mod. Phys. E25, 1650018 (2016)

2016
[23]

Tawﬁk, Nucl

A. Tawﬁk, Nucl. Phys. A764, 387 (2006)

2006
[24]

Tawﬁk, Europhys

A. Tawﬁk, Europhys. Lett. 75, 420 (2006)

2006
[25]

Tawﬁk, Nucl

A. Tawﬁk, Nucl. Phys. A922, 225 (2014)

2014
[26]

Tawﬁk, Adv

A. Tawﬁk, Adv. High Energy Phys. 2013, 574871 (2013)

2013
[27]

Tawﬁk, Phys

A. Tawﬁk, Phys. Rev. C88, 035203 (2013)

2013
[28]

A. N. Tawﬁk, H. Yassin, and E. R. A. Elyazeed, Phys. Rev. D92, 085002 (2015)

2015
[29]

Tawﬁk, J

A. Tawﬁk, J. Phys. G31, S1105 (2005)

2005
[30]

Agakishiev et al., Physical Review C 84, 014902 (2011)

G. Agakishiev et al., Physical Review C 84, 014902 (2011)

2011
[31]

V. V. Begun, V. Vovchenko, and M. I. Gorenstein, J. Phys. Conf. Ser. 779, 012080 (2017)

2017
[32]

Lorenz, Nucl

M. Lorenz, Nucl. Phys. A931, 785 (2014)

2014
[33]

Abgrall et al., Eur

N. Abgrall et al., Eur. Phys. J. C74, 2794 (2014)

2014
[34]

Agakishiev et al., Phys

G. Agakishiev et al., Phys. Rev. C90, 015202 (2014)

2014
[35]

Agakishiev et al., Phys

G. Agakishiev et al., Phys. Rev. C92, 024903 (2015)

2015
[36]

Anticic et al., Phys

T. Anticic et al., Phys. Rev. C80, 034906 (2009)

2009
[37]

Anticic et al., Physical Review C 79, 044904 (2009)

T. Anticic et al., Physical Review C 79, 044904 (2009)

2009
[38]

Vovchenko, M

V. Vovchenko, M. Gorenstein, L. Satarov, and H. St¨ ocke r, Chemical freeze-out conditions in hadron resonance gas, in New Horizons in Fundamental Physics , pages 127–137, Springer, 2017

2017
[39]

Alt et al., Phys

C. Alt et al., Phys. Rev. Lett. 94, 192301 (2005)

2005
[40]

Kraus, J

I. Kraus, J. Phys. G31, S147 (2005)

2005
[41]

Floris, Nucl

M. Floris, Nucl. Phys. A931, 103 (2014)

2014
[42]

Andronic, P

A. Andronic, P. Braun-Munzinger, K. Redlich, and J. Sta chel, Nature 561, 321 (2018)

2018
[43]

Abelev et al., Physical Review Letters 111, 222301 (2013)

B. Abelev et al., Physical Review Letters 111, 222301 (2013)

2013
[44]

Alexandre et al., Nuclear Physics A 931, 1093 (2014)

D. Alexandre et al., Nuclear Physics A 931, 1093 (2014)

2014
[45]

Abelev et al., Physical Review C 91, 024609 (2015)

B. Abelev et al., Physical Review C 91, 024609 (2015)

2015
[46]

Collaboration et al., arXiv preprint arXiv:1506.08 453 (2015)

A. Collaboration et al., arXiv preprint arXiv:1506.08 453 (2015)

2015
[47]

Collaboration et al., arXiv preprint arXiv:1710.07 531 (2017)

A. Collaboration et al., arXiv preprint arXiv:1710.07 531 (2017)

2017
[48]

Abelev et al., Physical Review Letters 109, 252301 (2012)

B. Abelev et al., Physical Review Letters 109, 252301 (2012)

2012
[49]

Cleymans, H

J. Cleymans, H. Oeschler, K. Redlich, and S. Wheaton, Ph ysical Review C 73, 034905 (2006)

2006
[50]

Adler et al., Physical review letters 89, 202301 (2002)

C. Adler et al., Physical review letters 89, 202301 (2002)

2002
[51]

Adler et al., Physical Review Letters 87, 112303 (2001)

C. Adler et al., Physical Review Letters 87, 112303 (2001)

2001
[52]

Adams et al., Physics Letters B 637, 161 (2006)

J. Adams et al., Physics Letters B 637, 161 (2006)

2006
[53]

Velkovska, Nucl

J. Velkovska, Nucl. Phys. A698, 507 (2002)

2002
[54]

Adams et al., Nucl

J. Adams et al., Nucl. Phys. A757, 102 (2005)

2005
[55]

Oeschler, Lect

H. Oeschler, Lect. Notes Phys 516, 1 (1999)

1999
[56]

V. I. Yukalova and E. P. Yukalov, PoS Baldin-ISHEPP-XXI, 046 (2012)

2012
[57]

K. A. Bugaev et al., Eur. Phys. J. A52, 227 (2016)

2016
[58]

Man Lo, M

P. Man Lo, M. Marczenko, K. Redlich, and C. Sasaki, Eur. P hys. J. A52, 235 (2016)

2016
[59]

J. Noronha-Hostler, Implications of Missing Resonanc es in Heavy Ions Collisions, in Workshop on Excited Hyperons in QCD Thermodynamics at Freeze-Out (YSTAR2016) Mini-Proc eedings, pages 118–127, 2016

2016
[60]

Capstick and N

S. Capstick and N. Isgur, Phys. Rev. D34, 2809 (1986), [AIP Conf. Proc.132,267(1985)]

1986
[61]

Bazavov et al., Phys

A. Bazavov et al., Phys. Rev. Lett. 113, 072001 (2014)

2014
[62]

Bellwied et al., Phys

R. Bellwied et al., Phys. Lett. B751, 559 (2015)

2015
[63]

Bazavov et al., Phys

A. Bazavov et al., Phys. Lett. B795, 15 (2019)

2019
[64]

Tawﬁk, N

A. Tawﬁk, N. Magdy, and A. Diab, Phys. Rev. C89, 055210 (2014)

2014
[65]

Philipsen, Lattice QCD at non-zero temperature and b aryon density, in Modern perspectives in lattice QCD: Quantum ﬁeld theory and high performance computing

O. Philipsen, Lattice QCD at non-zero temperature and b aryon density, in Modern perspectives in lattice QCD: Quantum ﬁeld theory and high performance computing. Procee dings, International School, 93rd Session, Les Houches, France, August 3-28, 2009 , pages 273–330, 2010

2009
[66]

A. N. Tawﬁk and A. M. Diab, Phys. Rev. C91, 015204 (2015)

2015
[67]

A. M. Abdel Aal Diab and A. N. Tawﬁk, EPJ Web Conf. 177, 09005 (2018). 8

2018
[68]

Nasser Tawﬁk and A

A. Nasser Tawﬁk and A. Magied Diab, Int. J. Mod. Phys. A30, 1550059 (2015)

2015
[69]

A. N. Tawﬁk, A. M. Diab, and M. T. Hussein, J. Phys. G45, 055008 (2018)

2018
[70]

A. N. Tawﬁk, A. M. Diab, and M. T. Hussein, Int. J. Mod. Phy s. A31, 1650175 (2016)

2016
[71]

A. N. Tawﬁk, A. M. Diab, and M. T. Hussein, Chin. Phys. C43, 034103 (2019)

2019

Pith tools

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

[1] [1]

Banks and A

T. Banks and A. Ukawa, Nucl. Phys. B225, 145 (1983)

work page 1983

[2] [2]

A. N. Tawﬁk, Int. J. Mod. Phys. A29, 1430021 (2014)

work page 2014

[3] [3]

B. P. Abbott et al., Phys. Rev. Lett. 119, 161101 (2017)

2017

[4] [4]

D. D. Ivanenko and D. F. Kurdgelaidze, Astrophysics 1, 251 (1965)

work page 1965

[5] [5]

Gyulassy and L

M. Gyulassy and L. McLerran, Nucl. Phys. A750, 30 (2005)

2005

[6] [6]

Busza, K

W. Busza, K. Rajagopal, and W. van der Schee, Ann. Rev. Nuc l. Part. Sci. 68, 339 (2018)

work page 2018

[7] [7]

Greiner, P

C. Greiner, P. Koch-Steinheimer, F. M. Liu, I. A. Shovkov y, and H. Stoecker, J. Phys. G31, S725 (2005)

work page 2005

[8] [8]

Why does the thermal model for hadron production in heavy ion collisions work?

B. Mueller and A. Schaefer, arXiv:1712.03567 [nucl-th] (2017)

work page Pith review arXiv 2017