REVIEW 3 major objections 6 minor 26 references
Speed of sound of QCD matter at chiral crossover
T0 review · 3 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read The paper claims that restoring chiral symmetry drives the hadron gas's squared speed of sound negative, so only partonic excitations can keep QCD matter mechanically stable across the chiral crossover.
desk verdict The partial c_s^2 decomposition is genuinely new, but the headline negative hadron speed of sound is carried by fit parameters and is not shown to be robust. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central machinery is the generalized Beth-Uhlenbeck formula for the thermodynamic potential of each multiquark cluster, $Ω_i$, expressed through temperature-dependent phase shifts $δ_i$ that encode bound states below the Mott threshold and broad resonances above it. The phase shifts are driven by the step-function quark mass prescription of Eq. (1), the Polyakov loop potential that suppresses colored clusters, and a Breit-Wigner width that grows above $T_c$. This machinery lets the paper split the total entropy density into partial hadronic, quark, gluon, and colored-cluster contributions and define an individual squared speed of sound for each component, which is how the hadronic instability is isolated.
What would settle it
Replace the step-function quark masses with a smooth crossover profile whose width matches the lattice chiral susceptibility, recompute the hadron component's squared speed of sound, and check whether the negative region persists; if it disappears, the instability is an artifact of the prescription rather than a physical property of hadronic matter.
Extended reading notes
Core claim
In the generalized Beth-Uhlenbeck model of Ref. [8], the squared speed of sound of the hadron gas, evaluated as a partial contribution to the total entropy density, turns negative at temperatures above the chiral crossover $T_c \simeq 156.5$ MeV. The paper attributes this to Mott dissociation triggered by the prescribed drop of quark masses from vacuum values ($M_{u,d}=627$ MeV, $M_s=770$ MeV) to current masses at $T_c$: hadrons become broad resonances above threshold, and the hadronic component alone cannot support mechanical equilibrium. The total entropy density stays smooth across $T_c$ because the discontinuous partial contributions of hadrons, quarks, gluons, and colored clusters cancel, while the total squared speed of sound is restored to the lattice values by the onset of partonic excitations. The paper concludes that mechanical stability of QCD matter requires partonic degrees of freedom immediately after $T_c$.
Load-bearing premise
The load-bearing assumption is the step-function prescription in Eq. (1): quark masses jump discontinuously from 627/770 MeV to 5.6/124 MeV exactly at $T_c$, with $T_c$ taken from lattice QCD; if the real chiral restoration is smooth, the hadronic instability could be an artifact of the imposed jump.
Editorial extensions
If this is right
- Any equation of state that keeps only hadrons above $T_c$ misrepresents QCD matter: the hadron gas alone has negative squared speed of sound and, through Eq. (13), a negative heat capacity, indicating unstable thermal equilibrium.
- The crossover's sharpness is hidden in bulk thermodynamics: the total entropy density is smooth even though the composition switches abruptly, so composition-sensitive observables are needed to locate the transition.
- The total squared speed of sound has a minimum slightly above $T_c$, a softening signature of deconfinement that the model reproduces in agreement with lattice QCD.
- Just above $T_c$, quarks contribute more than gluons to the stiffness of QCD matter; at higher temperatures the two contributions become comparable.
Reading between the lines
- If the negative hadronic $c_S^2$ is real, any purely hadronic equation of state used above $T_c$ is mechanically unstable, so hybrid or two-phase descriptions of heavy-ion collisions must switch to quarks and gluons essentially at $T_c$ rather than gradually.
- The same mechanism, extended to finite baryon density where colored multiquark clusters carry a substantial part of the entropy, suggests the instability may also shape the stiffness of neutron-star matter near deconfinement; the paper mentions the colored clusters' role but does not explore this consequence.
- The model predicts that the sharp composition switch is invisible in $s/T^3$ but should show up in composition-sensitive observables such as the second baryon-number susceptibility; this is a testable corollary the paper does not work out.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper constructs a unified equation of state for hot QCD matter using a generalized Beth-Uhlenbeck approach with cluster decomposition, in which hadrons appear as color-singlet multiquark clusters, colored multiquark clusters are included, and quarks and gluons are gauged by a Polyakov loop potential. Chiral symmetry restoration is modeled by a step-function switch of quark masses at a pseudocritical temperature T_c taken from lattice QCD (Eq. 1), while hadron spectral functions are parametrized through phase shifts with temperature-dependent masses and widths (Eqs. 7-8), with coefficients alpha1=11.4 and alpha2=1.9 adjusted to reproduce the lattice entropy density. The paper analyzes the total entropy density and, for the first time, the individual contributions of hadrons, quarks, gluons, and colored clusters to the speed of sound. It reports that the hadronic component acquires a negative squared speed of sound above T_c, which the authors interpret as a mechanical instability of the hadron resonance gas, and that adding partonic degrees of freedom restores a positive total speed of sound in reasonable agreement with lattice data.
Significance. The decomposition of the speed of sound into individual components is a novel and potentially instructive diagnostic for understanding the QCD crossover, and the unified treatment that yields a smooth total entropy despite an abrupt internal switch between hadronic and partonic degrees of freedom is a technically nontrivial achievement. If the central causal claim were robust, the paper would provide a simple argument for why hadronic-only descriptions must fail immediately above the chiral crossover. However, as presented, the negative hadronic speed of sound is not shown to be a robust prediction: it is generated by parameters fitted to the same lattice entropy data to which the total speed of sound is subsequently compared, and the attribution to chiral restoration is not supported by a sensitivity analysis or by an explicit test of the role of the chiral mass switch. The paper is therefore best regarded as a model-dependent study with an interesting but insufficiently substantiated interpretation.
major comments (3)
- [Sec. 3, Eqs. (7)-(8), Fig. 2] The negative squared speed of sound of the hadron gas above T_c is a direct consequence of the parametric mass and width model in Eqs. (7)-(8), with alpha1=11.4 and alpha2=1.9 adjusted to fit the lattice entropy density, rather than a consequence of the chiral quark-mass switch in Eq. (1). Above T_c, the thresholds in Eq. (10) are constant because current quark masses are temperature independent, so the decrease of s_h with temperature is driven by the linearly rising M_i and Gamma_i. The paper does not report a sensitivity study showing that the sign of c_{S,h}^2 is robust to variations of alpha1 and alpha2 or to a smooth version of the chiral transition. Because the total squared speed of sound in Fig. 2 is obtained as d ln T / d ln s from a total entropy that was fitted to lattice data using these parameters, the agreement with lattice c_s^2 is partly inherited and does not provide independent validation of the instability mechanism.
- [Sec. 3, Eq. (1), Fig. 2] The central claim that 'restoration of chiral symmetry drives speed of sound of hadron gas to negative values' is not directly supported: the step-function switch of quark masses at T_c only changes the thresholds at one point, while the negative c_{S,h}^2 above T_c arises from the T-dependence introduced by Eqs. (7)-(8). To justify the causal attribution, the authors should show that the negative hadronic speed of sound persists when the mass shift in Eq. (7) is removed (alpha1=0) or when the quark mass switch is spread over a finite temperature interval; absent such a test, the negative value is a property of the chosen fit, not of chiral restoration per se.
- [Sec. 3, Fig. 2 and Eqs. (7)-(8)] The comparison of the total speed of sound with lattice data is presented as a validation, but it is not an independent prediction because alpha1 and alpha2 in Eqs. (7)-(8) are fixed to reproduce the lattice entropy density from Ref. [10]. Since c_s^2 is defined as d ln T / d ln s, the derivative of the fitted entropy determines c_s^2; the agreement shown in Fig. 2 therefore tests consistency rather than predictive power. The manuscript should explicitly acknowledge this circularity and, ideally, perform an out-of-sample check (e.g., fit only s(T) below T_c and compare the extrapolated c_s^2 above T_c) or a sensitivity analysis over the fitted coefficients.
minor comments (6)
- [Abstract] The phrase 'being in a strike disagreements' is ungrammatical and should read 'being in striking disagreement'.
- [Section 1] The word 'diqurks' should be 'diquarks'.
- [Section 2, Eq. (4)] The typesetting of f_i appears broken in the manuscript; the formula should be displayed clearly for both the hadron and colored-cluster cases.
- [Section 2, text near Eq. (13)] 'the later can be written as' should be 'the latter can be written as'.
- [Section 3, Fig. 2 caption] 'The shaded areas represents' should be 'The shaded areas represent'.
- [Section 2, Eq. (5)] The phase shift is not defined for M > M_thr,i + Lambda_i N_i, where it drops to zero; this cutoff should be physically motivated and its influence on the results discussed.
Circularity Check
The headline hadron-gas instability (c^2_{S,h} < 0 above T_c) is generated by the fitted spectral-parameterization Eqs. (7)-(8), not by the chiral mass switch alone; the total c_S^2 agreement with lattice is a derivative of the entropy fit, making the central claim partially circular.
-
fitted input called prediction
[Section 2, Eqs. (7)-(8) with the statement of alpha_1, alpha_2; Fig. 2 (hadron partial c_S^2 negative above T_c)]
"The mass and decay width of this cluster are defined as M_i = M_{i,T=0} + alpha_1 Gamma_i, (7) Gamma_i = alpha_2 (T - T_c) theta(T - T_c), (8) respectively, while alpha_1 = 11.4 and alpha_2 = 1.9 are adjusted in agreement with the lattice data on entropy density."
c_S^2 = d ln T / d ln s, so the negative hadronic speed of sound requires ds_h/dT < 0 above T_c. But above T_c the hadronic phase shifts (5)-(6) depend on T only through M_i(T) and Gamma_i(T) of Eqs. (7)-(8); the thresholds (10) are T-independent since Eq. (1) fixes quark masses to constant current values. The continuous decline of s_h is thus generated entirely by alpha_1 and alpha_2, which are fit to the same lattice entropy data against which the total c_S^2 is later compared.
-
fitted input called prediction
[Section 3, speed-of-sound results compared to lattice QCD data]
"The corresponding speed of sound also provides a reasonable agreement with the lattice QCD data from Refs. Borsanyi et al. [11], Bazavov et al. [4]."
By the paper's own definition c_S^2 = d ln T / d ln s, the total speed of sound is the T-derivative of the entropy density s. Because alpha_1 and alpha_2 were 'adjusted in agreement with the lattice data on entropy density' and the Polyakov potential is likewise fitted to pure-gauge lattice data, the entropy profile — and hence most of its derivative — is calibrated to the data being 'confronted.' Agreement of the derivative of a fitted curve is a weaker, partly inherited consistency check than an independent prediction; it is partial, not complete circularity, but the partonic-rescue mechanism is validated against the same lattice entropy that fixed the model's parameters.
full rationale
The negative hadronic speed of sound is not an independent output of chiral dynamics in this paper. Above T_c the decay thresholds (10) are constant because Eq. (1) fixes quark masses at T-independent current values, so the continuous T-dependence of the hadronic phase shifts (5)-(6) enters only through M_i = M_{i,0} + alpha_1 Gamma_i and Gamma_i = alpha_2 (T - T_c) (Eqs. 7-8), with alpha_1 = 11.4 and alpha_2 = 1.9 'adjusted in agreement with the lattice data on entropy density.' The decline of the hadron partial entropy that makes c^2_{S,h} = d ln T / d ln s_h negative is therefore a byproduct of the fit to the very lattice entropy data against which the model is later compared; no sensitivity study establishes that the sign survives variation of alpha_1 and alpha_2. Section 2 explicitly concedes the causal input is prescribed rather than derived: 'we omit this step and model spontaneous breaking of chiral symmetry and its dynamical restoration by simply prescribing the temperature dependence of the quark masses, which treats T_c as an input parameter'; this admission further weakens the abstract's causal claim ('restoration of chiral symmetry drives speed of sound of hadron gas to negative values...'). The total speed-of-sound agreement with lattice is a derivative consistency check of the fitted entropy density, partially circular but not vacuous. The generalized Beth-Uhlenbeck formalism itself (Baym-Kadanoff), the pure-gauge Polyakov-potential fit, and the two-loop pion-polarization motivation of the parameterization (Refs. [7, 8]) carry independent external content, and the paper openly discloses its inputs, so the work is not wholly self-referential; but the central novel claim reduces substantially to fitted parameters, warranting score 6 rather than 8-10.
Assumptions & free parameters
free parameters (7)
- Vacuum quark masses M_{u,d}|0=627 MeV, M_s|0=770 MeV =
627 MeV, 770 MeV
- Binding energy per quark bond B =
471 MeV
- Pseudocritical temperature T_c =
156.6 MeV
- alpha_1, alpha_2 =
11.4, 1.9
- T_* =
94 MeV
- Lambda_pert =
222 MeV
- Lambda_i =
not given
assumptions (6)
- domain assumption Baym-Kadanoff Phi-functional / generalized Beth-Uhlenbeck cluster decomposition describes QCD thermodynamics with monomers and multiquark clusters.
- ad hoc to paper Abrupt step-function switching of quark masses at T_c (Eq. 1), with T_c as input, omitting self-consistent chiral dynamics.
- ad hoc to paper Parametric phase-shift model for multiquark clusters (Eqs. 5-6), including Breit-Wigner resonance shape above T_c and step bound-state peak below.
- domain assumption Polyakov loop potential U_phi from Lo et al. [19], fitted to pure gauge SU(3) lattice data, captures all purely gluonic contributions.
- standard math Terms from differentiation of phase shifts cancel against the temperature derivative of the Phi-functional at two-loop skeleton level [24,6].
- domain assumption Large vacuum quark masses are justified by a confining density functional approach [16,17].
invented entities (1)
-
Sexaquark (hypothetical spin-flavor-color singlet six-quark state)
Cite this review
Pith. "Pith review of Speed of sound of QCD matter at chiral crossover." pith.science (2026). https://pith.science/paper/UEX5HB5E
@misc{pith2026250621292,
author = {Pith},
title = {Pith review of: Speed of sound of QCD matter at chiral crossover},
year = {2026},
howpublished = {\url{https://pith.science/paper/UEX5HB5E}},
note = {Machine review of arXiv:2506.21292}
}
read the original abstract
Based on a generalized Beth-Uhlenbeck approach to thermodynamics of QCD motivated by cluster decomposition we present a unified equation of state of hot strongly interacting matter and analyze its properties in a wide range of temperatures. The hadrons are treated as color singlet multiquark clusters in medium with a background gluon field in the Polyakov gauge. The confining aspect of QCD is accounted for by the Polyakov loop mechanism and by a large vacuum quark mass motivated by a confining density functional approach. We demonstrate that an abrupt switching between hadronic and partonic degrees of freedom, which is one of striking manifestations of dynamical restoration of chiral symmetry, is accompanied by a smooth behavior of entropy density at chiral crossover. Individual contributions of different components of strongly interacting matter to its speed of sound are analyzed for the first time. It is shown that restoration of chiral symmetry drives speed of sound of hadron gas to negative values, manifesting its mechanical instability and being in a strike disagreements with the lattice QCD data. Accounting for the partonic excitations naturally resolves this contradiction.
Figures
Reference graph
Works this paper leans on
-
[10]
Lattice QCD equation of state at finite chemical poten- tial from an alternative expansion scheme
Borsányi, S., Fodor, Z., Guenther, J.N., Kara, R., Katz, S.D., Parotto, P., Pásztor, A., Ratti, C., Szabó, K.K., 2021. Lattice QCD equation of state at finite chemical poten- tial from an alternative expansion scheme. Phys. Rev. Lett. 126, 232001. doi:10.1103/PhysRevLett.126.232001, arXiv:2102.06660
arXiv 2021
-
[1]
Selfconsistent approximation in many body systems
Baym, G., 1962. Selfconsistent approximation in many body systems. Phys. Rev. 127, 1391–1401. doi:10.1103/PhysRev.127.1391
-
[2]
From hadrons to quarks in neu- tron stars: a review
Baym, G., Hatsuda, T., Kojo, T., Powell, P.D., Song, Y ., Takatsuka, T., 2018. From hadrons to quarks in neu- tron stars: a review. Rept. Prog. Phys. 81, 056902. doi:10.1088/1361-6633/aaae14,arXiv:1707.04966
arXiv 2018
-
[3]
Conservation Laws and Correlation Functions
Baym, G., Kadanoff, L.P., 1961. Conservation Laws and Correlation Functions. Phys. Rev. 124, 287–299. doi:10.1103/PhysRev.124.287
-
[4]
Bazavov, A., et al. (HotQCD), 2014. Equation of state in ( 2+1 )-flavor QCD. Phys. Rev. D 90, 094503. doi:10.1103/PhysRevD.90.094503,arXiv:1407.6387
arXiv 2014
-
[5]
Bazavov, A., et al. (HotQCD), 2019. Chiral crossover in QCD at zero and non-zero chemical potentials. Phys. Lett. B 795, 15–21. doi:10.1016/j.physletb.2019.05.013, arXiv:1812.08235
arXiv 2019
-
[6]
Approximately selfconsistent resummations for the thermodynamics of the quark gluon plasma
Blaizot, J.P., Iancu, E., Rebhan, A., 2001. Approximately selfconsistent resummations for the thermodynamics of the quark gluon plasma. 1. Entropy and density. Phys. Rev. D 63, 065003. doi:10.1103/PhysRevD.63.065003, arXiv:hep-ph/0005003
arXiv 2001
-
[7]
Generalized Beth– Uhlenbeck approach to mesons and diquarks in hot, dense quark matter
Blaschke, D., Buballa, M., Dubinin, A., Roepke, G., Zablocki, D., 2014. Generalized Beth– Uhlenbeck approach to mesons and diquarks in hot, dense quark matter. Annals Phys. 348, 228–255. doi:10.1016/j.aop.2014.06.002,arXiv:1305.3907
arXiv 2014
Show all 26 references
-
[8]
Blaschke, D., Cierniak, M., Ivanytskyi, O., Röpke, G.,
-
[9]
The QCD equation of state with dynamical quarks
Borsanyi, S., Endrodi, G., Fodor, Z., Jakovac, A., Katz, S.D., Krieg, S., Ratti, C., Szabo, K.K., 2010. The QCD equation of state with dynamical quarks. JHEP 11, 077. doi:10.1007/JHEP11(2010)077,arXiv:1007.2580
2010 arXiv
-
[11]
Full result for the QCD equa- tion of state with 2+1 flavors
Borsanyi, S., Fodor, Z., Hoelbling, C., Katz, S.D., Krieg, S., Szabo, K.K., 2014. Full result for the QCD equa- tion of state with 2+1 flavors. Phys. Lett. B 730, 99–104. doi:10.1016/j.physletb.2014.01.007,arXiv:1309.5258
2014 arXiv
-
[12]
On the Mass and on the Dynami- cal Properties of the Sexaquark
Buccella, F., 2020. On the Mass and on the Dynami- cal Properties of the Sexaquark. PoS CORFU2019, 024. doi:10.22323/1.376.0024
2020 doi
-
[13]
Heavy Ion Collisions: The Big Picture, and the Big Questions
Busza, W., Rajagopal, K., van der Schee, W., 2018. Heavy Ion Collisions: The Big Picture, and the Big Questions. Ann. Rev. Nucl. Part. Sci. 68, 339–376. doi:10.1146/annurev-nucl-101917-020852, arXiv:1802.04801
2018 arXiv
-
[14]
The phase di- agram of dense QCD
Fukushima, K., Hatsuda, T., 2011. The phase di- agram of dense QCD. Rept. Prog. Phys. 74, 014001. doi:10.1088/0034-4885/74/1/014001, arXiv:1005.4814
2011 arXiv
-
[15]
Statistical thermodynamics of strong interactions at high-energies
Hagedorn, R., 1965. Statistical thermodynamics of strong interactions at high-energies. Nuovo Cim. Suppl. 3, 147– 186
1965
-
[16]
Density functional approach to quark matter with confine- ment and color superconductivity
Ivanytskyi, O., Blaschke, D., 2022a. Density functional approach to quark matter with confine- ment and color superconductivity. Phys. Rev. D 105, 114042. doi:10.1103/PhysRevD.105.114042, arXiv:2204.03611
-
[17]
Recovering the Conformal Limit of Color Superconducting Quark Mat- ter within a Confining Density Functional Approach
Ivanytskyi, O., Blaschke, D.B., 2022b. Recovering the Conformal Limit of Color Superconducting Quark Mat- ter within a Confining Density Functional Approach. Particles 5, 514–534. doi:10.3390/particles5040038, arXiv:2209.02050
-
[18]
Chiral condensate in hadronic matter
Jankowski, J., Blaschke, D., Spalinski, M., 2013. Chiral condensate in hadronic matter. Phys. Rev. D 87, 105018. doi:10.1103/PhysRevD.87.105018,arXiv:1212.5521
2013 arXiv
-
[19]
Polyakov loop fluctuations in SU(3) lattice gauge theory and an effective gluon potential
Lo, P.M., Friman, B., Kaczmarek, O., Redlich, K., Sasaki, C., 2013. Polyakov loop fluctuations in SU(3) lattice gauge theory and an effective gluon potential. Phys. Rev. D 88, 074502. doi:10.1103/PhysRevD.88.074502, arXiv:1307.5958
2013 arXiv
-
[20]
Smear- ing the Quark Model
Poggio, E.C., Quinn, H.R., Weinberg, S., 1976. Smear- ing the Quark Model. Phys. Rev. D 13, 1958. doi:10.1103/PhysRevD.13.1958
1976 doi
-
[21]
Thermodynamics of the PNJL model
Ratti, C., Roessner, S., Thaler, M.A., Weise, W., 2007. Thermodynamics of the PNJL model. Eur. Phys. J. C 49, 213–217. doi:10.1140/epjc/s10052-006-0065-x, arXiv:hep-ph/0609218
2007 arXiv
-
[22]
Cluster virial expansion for nuclear matter within a quasipar- ticle statistical approach
Ropke, G., Bastian, N.U., Blaschke, D., Klahn, T., Typel, S., Wolter, H.H., 2013. Cluster virial expansion for nuclear matter within a quasipar- ticle statistical approach. Nucl. Phys. A 897, 70–92. doi:10.1016/j.nuclphysa.2012.10.005, arXiv:1209.0212
2013 arXiv
-
[23]
Analytic perturbation the- ory in analyzing some QCD observables
Shirkov, D.V ., 2001. Analytic perturbation the- ory in analyzing some QCD observables. Eur. Phys. J. C 22, 331–340. doi:10.1007/s100520100794, arXiv:hep-ph/0107282
2001 arXiv
-
[24]
Selfconsistent ap- proximations in relativistic plasmas: Quasiparticle anal- ysis of the thermodynamic properties
Vanderheyden, B., Baym, G., 1998. Selfconsistent ap- proximations in relativistic plasmas: Quasiparticle anal- ysis of the thermodynamic properties. J. Statist. Phys. 93, 843. doi:10.1023/B:JOSS.0000033166.37520.ae, arXiv:hep-ph/9803300
1998
-
[25]
(Particle Data Group), 2022
Workman, R.L., et al. (Particle Data Group), 2022. Review of Particle Physics. PTEP 2022, 083C01. doi:10.1093/ptep/ptac097. 7
2022 doi
-
[2024]
Thermodynamics of quark matter with multiquark clusters in an effective Beth-Uhlenbeck type approach. Eur. Phys. J. A 60, 14. doi:10.1140/epja/s10050-023- 01229-8,arXiv:2308.07950. 6
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.