REVIEW 3 major objections 5 minor 1 cited by
Dense QC$_2$D. What's up with that?!?
T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read In dense two-color QCD, the squared speed of sound exceeds the conformal limit 1/3.
desk verdict A plausible result that overreaches on error control: the ∂κ/∂a uncertainty is a load-bearing gap, so 'clearly breaches' is not yet backed. 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 object is the squared speed of sound $C_s^2 = dP/d\varepsilon$, built from the pressure $P$ and the energy density $\varepsilon = T^\mu_\mu + 3P$, where the trace anomaly $T^\mu_\mu$ is split into gluonic and fermionic parts renormalised by the $\beta$ functions $\partial\beta/\partial a$ and $\partial\kappa/\partial a$. The paper's new ingredient is a scale-setting chain that fixes the lattice spacing from the Cornell form of the static quark potential, using Coulomb-gauge Wilson lines rather than Wilson loops, and determines the $\beta$ functions from fits to $\beta(a)$ and $\kappa(a)$ along the line of constant physics. The thermodynamic observables are extrapolated to zero diquark source $j$ after simulations with $aj = 0.01$--$0.03$ that lift the low-lying modes of the superfluid phase. The conformal limit $C_s^2 = 1/3$ serves as the benchmark against which the dense-matter equation of state is compared.
What would settle it
A calculation of the same speed-of-sound curve on a second, finer lattice spacing on the same line of constant physics, extrapolated to the continuum at fixed $\mu_q/m_\pi$, would settle the claim: if the continuum value does not exceed $1/3$ for $\mu_q/m_\pi$ around 1.0--1.4, the reported breach is a lattice artefact.
Extended reading notes
Core claim
Using unimproved Wilson fermions and gauge action on a lattice with $\beta=2.1$, $\kappa=0.1577$, and spacing $a=0.130$ fm set by the Cornell static-quark potential, the authors compute the equation of state of dense two-color QCD at quark masses on the line of constant physics $m_\pi/m_\rho \simeq 0.81$. The pressure is obtained by integrating the quark number density in $\mu_q$, with lattice artefacts mitigated by a Stefan--Boltzmann quotient (Scheme II), and the energy density follows from the trace anomaly using $\beta$ functions determined from the new scale setting. The central result, stated in Section 3.5, is that $C_s^2 = dP/d\varepsilon$ rises sharply above the onset chemical potential $\mu_0 = m_\pi/2$ and clearly breaches the conformal limit $1/3$, before later falling. The authors report that this behaviour is consistent with other recent QC2D calculations and with isospin QCD simulations, while remaining below the upper bound derived from relativistic hydrodynamics.
Load-bearing premise
The physical scale is set by assuming the Cornell static-quark potential with string tension $\sigma=(440\text{ MeV})^2$ holds for two-color QCD, and that Coulomb-gauge Wilson lines match Wilson-loop results within $2\sigma$; if either assumption fails, the $\mu$ axis of the speed-of-sound curve shifts and the comparison with the conformal bound could change.
Editorial extensions
If this is right
- If the lattice result survives the continuum limit, the speed of sound in dense two-color QCD exceeds $1/3$ for a range of chemical potentials above onset, meaning the conformal bound is not a universal property of QCD-like matter.
- The consistency of the finer-lattice result with earlier coarser-lattice results [8-10] supports the conclusion that the excess over $1/3$ is not a lattice-spacing artefact.
- The new scale-setting procedure based on the Cornell potential gives beta functions consistent with Karsch coefficients on the coarse lattice, validating the fine-lattice thermodynamics.
- The computed $C_s^2$ remains below the upper bound derived from relativistic hydrodynamics [22], so the two constraints on the speed of sound are compatible.
- The diquark condensate extrapolated to $j=0$ is non-zero above onset and the chiral perturbation theory exponent $C_2 = 1/3$ fails at high densities, indicating that the superfluid phase is not $\chi$PT-dominated there.
Reading between the lines
- If the same breach of the conformal bound persists in the continuum limit, the equation of state of dense strongly interacting matter would be stiffer than conformal in a window above nuclear saturation, with potential consequences for neutron-star radii inferred from gravitational-wave and X-ray observations.
- The Wilson-line method for extracting the static potential, once cross-checked against Wilson loops, could reduce the cost of scale setting on finer lattices and make continuum extrapolations of dense-matter thermodynamics more tractable.
- The observed failure of $\chi$PT at high density suggests the superfluid may cross over to a BEC-BCS-type regime; this could be tested on the same ensembles by measuring the diquark pair radius or the fermion dispersion relation.
- Applying the same thermodynamic analysis to isospin-asymmetric QCD could reveal whether the $C_s^2 > 1/3$ window near onset is a universal feature of dense QCD-like theories.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings paper reports new lattice simulations of two-color QCD (QC2D) at non-zero baryon chemical potential, with the stated goal of computing the speed of sound C_s^2 = dP/dε on a finer lattice (β = 2.1, a ≃ 0.130 fm) than earlier studies. The authors tune to a line of constant physics m_π/m_ρ ≃ 0.81, extract beta functions from fits to scale-setting data, compute the pressure by integrating the quark number density with scheme-II lattice-artefact corrections, and obtain the energy density from the trace anomaly via Eqs. (2)-(3). The central result, stated in §3.5, is that C_s^2 rises sharply above the onset chemical potential and 'clearly breaches' the conformal limit 1/3. The paper also describes a new Wilson-line-based static potential method for scale setting and releases the data and analysis code on Zenodo.
Significance. If confirmed, the finding that C_s^2 exceeds the conformal bound in dense two-color QCD supports and extends similar observations in independent lattice studies and in isospin QCD, with potential implications for neutron-star equations of state. The paper's use of a finer lattice and a different scale-setting procedure than earlier works is a valuable cross-check. The explicit release of ensembles, analysis code, and speed-of-sound data (Zenodo records [23], [25], [26]) is a significant strength for reproducibility. However, the headline claim is not yet backed by a quantitative error budget: the uncertainties in the beta functions, the diquark-source extrapolation, and the derivative of the interpolated pressure-energy curve are not propagated into C_s^2. The authors themselves acknowledge in §4 that the error analysis is 'crude', so the current manuscript should be treated as a preliminary report rather than a definitive measurement.
major comments (3)
- [§3.4, Eq. (13); Table 2] The fermionic trace anomaly in Eq. (13) is proportional to a (∂κ/∂a). For the fine lattice (β = 2.1) Table 2 quotes ∂κ/∂a = 0.152 ± 0.32 fm^{-1}, i.e. a relative uncertainty above 200% and a value consistent with zero. With κ^{-1} ≃ 6.34 and the subtracted chiral condensate of Fig. 5b (magnitude ~0.008 in lattice units), the 1σ uncertainty in T_q^μμ is comparable to or larger than the total trace-anomaly signal shown in Fig. 4. The manuscript does not propagate this uncertainty through Eq. (2) into ε or into C_s^2 = dP/dε. If ∂κ/∂a lies near the lower end of its error bar, the fermionic contribution changes sign and magnitude, and the breach of the conformal bound in §3.5 may not survive. This is the central load-bearing issue: the claim 'C_s^2 clearly breaches' is not supported by the current error accounting. The authors should either (a) propagate all fit and statistical uncertainties into ε and C_s^2, or (b) constrain ∂κ/∂a with additional scale-setting data, or (c) explicitly weaken the claim to a qualitative observation pending a full error budget.
- [§4; §3.2] The Discussion acknowledges that measurements are taken every trajectory with possible autocorrelations, that errors are bootstrapped rather than jackknifed, and that fit errors are 'currently read off of the fitting function'. Separately, §3.2 states that the lowest diquark-source runs were performed on the smaller L_s = 16^3 spatial volume, while the highest source used N_s = 24. These are exactly the conditions that control the diquark extrapolation (Eq. (8)) used to obtain the zero-source observables that feed into the pressure and trace anomaly. Without an estimate of the systematic error from the mixed-volume diquark extrapolation, and without a confidence band on the derivative of the cubic-spline interpolation in Fig. 6b, the reported C_s^2 values in Fig. 7 have no demonstrated validity. The authors should provide the central values and uncertainties of C_s^2 as a function of μ_q, even if preliminary, so the reader can judge whether the conformal-limit breach is statistically significant.
- [§2.2, Eq. (5)] The lattice spacing is set by assuming the Cornell static-quark potential with the three-color QCD string tension σ = (440 MeV)^2, and the Wilson-line method in Coulomb gauge is claimed to agree with Wilson-loop results 'within 2σ' but only as 'results in preparation'. While the dimensionless ratio C_s^2 is invariant under a global change of scale, the physical units of the chemical potential and the comparison with other groups' scale-setting choices (e.g. T = 200 MeV in refs. [8-10]) depend on this assumption. The manuscript should either provide the validation of the Wilson-line/Wilson-loop equivalence as an appendix or reference, or explicitly state the sensitivity of the speed-of-sound curve to the assumed string tension.
minor comments (5)
- [§3.5] The sentence 'There is a sharp increase above at the onset chemical potential' appears to be missing a word; it should read 'above the onset chemical potential'.
- [Figure 4] The y-axis notation ('−10×') is unclear; please specify the power of ten and the units of T^μμ (e.g., lattice units a^{-4} or normalized by μ_c^4).
- [§3.4, Eq. (13)] The text says the μ = 0 values are subtracted to obtain Fig. 5; please clarify whether the subtraction is applied to the entire operator (4N_f N_c − ⟨ψ̄ψ⟩) or only to ⟨ψ̄ψ⟩, since the constant term may cancel in the difference.
- [Figure 6b] Please specify whether the 'cubic spline interpolation' is a smoothing spline or an interpolating spline through all points, and how many knots or the smoothing parameter were used, since the speed of sound is obtained by differentiating this curve.
- [Abstract] The abstract describes the paper as 'recent updates and results'; since this is a proceedings contribution, it may be helpful to state that these are preliminary results and that a full error analysis is in progress, consistent with the Discussion.
Circularity Check
No significant circularity: the speed-of-sound result is derived from measured thermodynamic quantities, not from a fitted parameter renamed as a prediction.
full rationale
The central claim, C_s^2 > 1/3, is not equivalent to any input of the calculation. Equation (1) defines C_s^2 = dP/dε. The pressure P is obtained in Eq. (10) by integrating the measured quark-number density n_q with a parameter-free Stefan-Boltzmann correction, and the energy density ε is obtained from Eqs. (2), (3), (12), and (13) using measured plaquette and chiral-condensate values multiplied by the beta functions ∂β/∂a and ∂κ/∂a. The beta functions are fitted to lattice spacings along the line of constant physics, but the fitted quantities are the beta functions themselves, not the speed of sound; C_s^2 only emerges after differentiating an interpolated P(ε) curve. The conformal bound 1/3 is not a fit parameter. The paper also checks its results against independent recent QC2D works [8-10] and isospin QCD simulations, and against Stefan-Boltzmann benchmarks. The self-citations [1-7,16,23-26] provide code, prior methodology, and a standard proof of the absence of a sign problem; none is load-bearing for the headline claim, and the independent works [8-10] are by different authors. Two caveats are genuine concerns but are not circularity: the lattice spacing relies on an assumed (440 MeV)^2 string tension and on a Wilson-line/Wilson-loop equivalence that the paper states is still 'results in preparation' (Sec. 2.2), and the fermionic beta function ∂κ/∂a carries a large uncertainty (Table 2). These are calibration and error-accounting risks; they do not make the derivation reduce to its own inputs. No circular step is therefore established.
Assumptions & free parameters
free parameters (5)
- Cornell potential parameters alpha_s and V0 =
not reported
- beta(a) fit coefficients =
quadratic coefficients from Fig. 1a
- kappa(a) fit coefficients =
fit coefficients from Fig. 1b
- diquark condensate fit parameters C0, C1, C2 =
per-chemical-potential fits; C2 ~ 1/3 near onset
- linear diquark extrapolation slopes for n_q, plaquette, chiral condensate =
not reported
assumptions (6)
- standard math SU(2) gauge theory with an even number of quark flavours has no complex action problem.
- domain assumption The string tension sigma=(440 MeV)^2 provides the physical scale for QC2D.
- domain assumption The static quark potential is described by the Cornell form in Eq. (5).
- domain assumption Wilson lines fixed to Coulomb gauge give the same static quark potential as Wilson loops.
- domain assumption The Stefan-Boltzmann correction Scheme II removes lattice IR and UV artefacts from the pressure.
- domain assumption The quark number density correction can be evaluated on a larger spatial volume (4 N_s) than the simulated volume.
Cite this review
Pith. "Pith review of Dense QC$_2$D. What's up with that?!?." pith.science (2026). https://pith.science/paper/ZJTH2STO
@misc{pith2026241215872,
author = {Pith},
title = {Pith review of: Dense QC$_2$D. What's up with that?!?},
year = {2026},
howpublished = {\url{https://pith.science/paper/ZJTH2STO}},
note = {Machine review of arXiv:2412.15872}
}
abstract
We present recent updates and results from QC$_2$D (Two Colour QCD) simulations at non-zero baryon density, including progress toward determining the speed of sound.
Figures
Figures from the paper (4 more)
Forward citations
Cited by 1 Pith paper
-
Phase and equation of state of finite density QC$_2$D at lower temperature
Proceedings summarizing lattice QC2D results: the conformal bound c_s^2/c^2 = 1/3 is exceeded in the BCS phase at T = 40 and 80 MeV, with a rich hadronic/BEC/BCS phase structure.
Reference graph
Works this paper leans on
-
[23]
D. Lawlor, S. Hands, K. Seyong and J.-I. Skullerud,su2hmc, July, 2024. 10.5281/zenodo.12910604
-
[25]
J.-I. Skullerud and D. Lawlor,QC2D Analysis Suite, Dec., 2024. 10.5281/zenodo.14503846
-
[26]
D. Lawlor and J.-I. Skullerud,Dense QC2D. What’s up with that?!? The dataset., Nov., 2024. 10.5281/zenodo.14201453. 10
- [1]
- [2]
-
[3]
T. Boz, S. Cotter, L. Fister, D. Mehta and J.-I. Skullerud,Phase transitions and gluodynamics in 2-colour matter at high density, Eur. Phys. J. A49 (2013) 87 [1303.3223]
arXiv 2013
- [4]
-
[5]
S. Cotter,Non-perturbative determination of Karsch Coefficients in 2 Colour QCD, philosophiæ doctor, National University of Ireland Maynooth, 2015
work page 2015
Show all 26 references
-
[6]
Boz,Quark and Gluon Propagation in Two-Colour Quantum Chromodynamics at Finite Density, philosophiæ doctor, National University of Ireland Maynooth, 2018
T. Boz,Quark and Gluon Propagation in Two-Colour Quantum Chromodynamics at Finite Density, philosophiæ doctor, National University of Ireland Maynooth, 2018
2018
-
[7]
T. Boz, P. Giudice, S. Hands and J.-I. Skullerud,Dense two-color QCD towards continuum and chiral limits,Phys. Rev. D101(2020) 074506 [1912.10975]
2020 arXiv
-
[8]
K. Iida, E. Itou, K. Murakami and D. Suenaga,Lattice study on finite density QC2D towards zero temperature, 2405.20566
-
[9]
Itou and K
E. Itou and K. Iida,Speed of sound exceeding the conformal bound in dense 2-color QCD, in40th International Symposium on Lattice Field Theory, 11, 2023 [2311.15259]
2023 arXiv
-
[10]
Iida and E
K. Iida and E. Itou,Velocity of sound beyond the high-density relativistic limit from lattice simulation of dense two-color QCD,PTEP 2022 (2022) 111B01 [2207.01253]
2022 arXiv
-
[11]
Begun, V.G
A. Begun, V.G. Bornyakov, V.A. Goy, A. Nakamura and R.N. Rogalyov,Study of two color QCD on large lattices,Phys. Rev. D105 (2022) 114505 [2203.04909]
2022 arXiv
-
[12]
and Cuteri, Francesca and Endrődi, Gergely,Equation of state and speed of sound of isospin-asymmetric QCD on the lattice,JHEP 07(2023) 055 [2212.14016]
Brandt, Bastian B. and Cuteri, Francesca and Endrődi, Gergely,Equation of state and speed of sound of isospin-asymmetric QCD on the lattice,JHEP 07(2023) 055 [2212.14016]
2023 arXiv
-
[13]
Brandt, B. B. and Endrődi, G. and Schmalzbauer, S.,QCD phase diagram for nonzero isospin-asymmetry,Phys. Rev. D97(2018) 054514 [1712.08190]
2018 arXiv
-
[14]
and Endrődi, Gergely and Fraga, Eduardo S
Brandt, Bastian B. and Endrődi, Gergely and Fraga, Eduardo S. and Hippert, Mauricio and Schaffner-Bielich, Jurgen and Schmalzbauer, Sebastian,New class of compact stars: Pion stars,Phys. Rev. D98(2018) 094510 [1802.06685]
2018 arXiv
-
[15]
Nagata,Finite-density lattice QCD and sign problem: Current status and open problems,Prog
K. Nagata,Finite-density lattice QCD and sign problem: Current status and open problems,Prog. Part. Nucl. Phys.127 (2022) 103991 [2108.12423]
2022 arXiv
-
[16]
Hands, I
S. Hands, I. Montvay, S. Morrison, M. Oevers, L. Scorzato and J. Skullerud,Numerical study of dense adjoint matter in two color QCD, Eur. Phys. J. C17(2000) 285 [hep-lat/0006018]
2000 arXiv
-
[17]
Koehn et al.,An overview of existing and new nuclear and astrophysical constraints on the equation of state of neutron-rich dense matter, 2402.04172
H. Koehn et al.,An overview of existing and new nuclear and astrophysical constraints on the equation of state of neutron-rich dense matter, 2402.04172. 9 Dense QC2D. What’s up with that?!? Dale Lawlor
-
[18]
Riley et al.,A𝑁𝐼𝐶𝐸𝑅 View of PSR J0030+0451: Millisecond Pulsar Parameter Estimation, Astrophys
T.E. Riley et al.,A𝑁𝐼𝐶𝐸𝑅 View of PSR J0030+0451: Millisecond Pulsar Parameter Estimation, Astrophys. J. Lett.887 (2019) L21 [1912.05702]
2019 arXiv
-
[19]
Miller et al.,PSR J0030+0451 Mass and Radius from𝑁𝐼𝐶𝐸𝑅 Data and Implications for the Properties of Neutron Star Matter,Astrophys
M.C. Miller et al.,PSR J0030+0451 Mass and Radius from𝑁𝐼𝐶𝐸𝑅 Data and Implications for the Properties of Neutron Star Matter,Astrophys. J. Lett.887 (2019) L24 [1912.05705]
2019 arXiv
-
[20]
LIGO Scientific, Virgo collaboration,GW170817: Measurementsofneutronstar radiiandequation of state,Phys. Rev. Lett.121(2018) 161101 [1805.11581]
2018 arXiv
-
[21]
LIGO Scientific, Virgo collaboration, GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral, Phys. Rev. Lett.119 (2017) 161101 [1710.05832]
2017 arXiv
-
[22]
Hippert, J
M. Hippert, J. Noronha and P. Romatschke,Upper Bound on the Speed of Sound in Nuclear Matter from Transport, 2402.14085
-
[24]
Lüscher, Martin,A Portable high quality random number generator for lattice field theory simulations, Comput. Phys. Commun.79 (1994) 100 [hep-lat/9309020]
1994 arXiv
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.