{"id":"49818eb8-e5c1-4ad1-8449-8fa0667d6aa0","arxiv_id":"2412.15872","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"On a finer lattice than previous studies, the speed of sound in dense two-color QCD again rises above the conformal bound near the superfluid transition.","lead":"This paper reports new computer simulations of dense two-color QCD, a toy version of the strong nuclear force, and uses them to estimate how the speed of sound changes as matter gets denser. The results agree with earlier studies suggesting the speed of sound can exceed the 'conformal' value thought to apply to dense matter.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim that C_s^2 clearly exceeds 1/3 is not yet supported because the fermionic trace anomaly is multiplied by ∂κ/∂a, which is measured only to 0.152 ± 0.32 (Table 2); propagating this uncertainty could move ε enough to erase the breach.","rationale":"I read the paper as an early-results proceedings whose central claim is that C_s^2 rises sharply above onset and clearly exceeds the conformal bound on a finer lattice with a different scale-setting scheme. For that claim to hold, the derivative dP/dε must be robust against the stated uncertainties. The most load-bearing weakness I find is not the overall lattice scale, which mostly cancels in the dimensionless ratio C_s^2, but the poorly determined beta-function coefficient ∂κ/∂a. Equation (13) multiplies the measured chiral-condensate shift by this coefficient; Table 2 quotes ∂κ/∂a = 0.152 ± 0.32 for the fine lattice, so the coefficient is consistent with zero and even changes sign within 1σ. Since the fermionic trace anomaly enters ε linearly, and C_s^2 is the inverse slope of ε versus P, a sign or magnitude change in T_q^μμ can directly alter whether C_s^2 exceeds 1/3. The paper's own 'Discussion' section concedes that autocorrelations are untreated, low-diquark-source runs are on small volumes, and the error analysis is crude; these are honest limitations, but they reinforce the concern that the headline claim is not yet backed by a quantitative uncertainty band. The reader's conditional verdict is appropriate; I would keep it, but with the specific requirement that the ∂κ/∂a uncertainty be propagated into ε and C_s^2 before 'clearly breaches' is stated. Independent support from Refs. [8–10] and the released data/code are genuine strengths, and I do not see an internal inconsistency or a fatal flaw; the issue is the gap between the strength of the claim and the propagated errors.","tokens_in":7896,"tokens_out":8226,"duration_ms":78608,"concrete_test":"Recompute the fermionic trace anomaly and ε with ∂κ/∂a set to 0.152, −0.168, and 0.472 fm^{-1} in Eq. (13), keep all other inputs fixed, and re-derive C_s^2 from the spline derivative dP/dε. If the lower edge of the resulting C_s^2 band crosses 1/3 or if the sign of the deviation changes, the headline claim requires either a tighter determination of ∂κ/∂a or a softened statement. The same check should be repeated with jackknife blocks that account for autocorrelations, not just bootstrap resampling of individual configurations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result, C_s^2 > 1/3, is computed as dP/dε, with ε = Tμμ + 3P. The fermionic part of the trace anomaly, Eq. (13), is T_q^μμ = -a (∂κ/∂a) κ^{-1} (4 N_f N_c - ⟨ψ̄ψ⟩), after μ = 0 subtraction. For the fine lattice (β = 2.1, a = 0.130 fm), Table 2 gives ∂κ/∂a = 0.152 ± 0.32 fm^{-1}; this is a relative uncertainty above 200% and is consistent with zero. In the combination a(∂κ/∂a), the fractional error is the same. Multiplying by κ^{-1} ≈ 6.34 and by the measured condensate shift (Fig. 5b, magnitude ∼ 0.008 in lattice units), the 1σ uncertainty in T_q^μμ is comparable to the entire trace-anomaly signal shown in Fig. 4. The paper does not state that this uncertainty is propagated into ε or into C_s^2; the Discussion admits that errors from fits are 'currently read off of the fitting function' and that autocorrelations are untreated. If ∂κ/∂a is at its lower end, or negative, T_q^μμ changes sign and magnitude, ε shifts, and dP/dε changes; the breach above the conformal bound may not survive. The reader's scale-setting concern is related, but a global scale error largely cancels in the dimensionless ratio C_s^2; the ∂κ/∂a uncertainty is a more direct and numerically inescapable path from a poorly determined input to the headline claim. This is not a claim that the result is wrong: independent groups report similar behavior, and the released data and code are valuable. But 'clearly breaches' is not supported by the current error accounting until this coefficient is constrained or its uncertainty is propagated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":8266,"tokens_out":5787,"duration_ms":51902,"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":[{"comment":"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.","section":"§3.4, Eq. (13); Table 2"},{"comment":"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.","section":"§4; §3.2"},{"comment":"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.","section":"§2.2, Eq. (5)"}],"minor_comments":[{"comment":"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'.","section":"§3.5"},{"comment":"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).","section":"Figure 4"},{"comment":"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.","section":"§3.4, Eq. (13)"},{"comment":"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.","section":"Figure 6b"},{"comment":"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.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The paper is a proceedings contribution and the authors are transparent about the preliminary state of the analysis. However, the headline claim that C_s^2 'clearly breaches' the conformal limit is not backed by a quantitative error budget, and the specific ∂κ/∂a uncertainty in Table 2 is large enough that the result could be an artefact of the input. This is fixable within the scope of the manuscript: propagating the existing uncertainties, reporting C_s^2 with error bars, and softening the claim would make the paper acceptable. I would not recommend rejection because the data release and the comparison with independent work are valuable, and the central physical finding is plausible and consistent with other groups. The editor may wish to ask the authors to include the validation of the Wilson-line method or to clearly mark it as preliminary."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Simon Hands and colleagues have put out a short, honest proceedings paper. The headline—C_s^2 exceeds the conformal bound in dense two-color QCD—is not new; Itou and Iida already reported it, and the authors say their finer-lattice data are consistent with refs [8–10]. What is new: a fine lattice at β=2.1 with a=0.130 fm, beta functions determined on a line of constant physics rather than Karsch coefficients, a Wilson-line scale-setting method, and a released analysis suite with data. That last item is real value; the code and datasets are on Zenodo.\n\nThe central calculation follows the standard thermodynamic-integration route, with a check against the Stefan-Boltzmann pressure and good consistency for the coarse-lattice beta functions. The paper is transparent about its weaknesses in Section 4: autocorrelations untreated, errors bootstrapped rather than jackknifed, fit errors read off the fitting function, low diquark-source runs at small volume. Those are all acknowledged.\n\nThe soft spot that matters more is one the authors do not flag. The fermionic trace anomaly, Eq. (13), is multiplied by a ∂κ/∂a measured as 0.152 ± 0.32 fm^{-1} on the fine lattice—relative uncertainty above 200%, consistent with zero. That uncertainty is not propagated into ε or into C_s^2. The stress-test concern holds up: if ∂κ/∂a sits at its lower end, the fermionic contribution to the trace anomaly changes sign and magnitude, and the breach above 1/3 may not survive. This is a load-bearing gap, not a cosmetic one. The same issue affects the comparison with other groups, since the same coefficient enters the scale setting.\n\nThere are two smaller reservations. The lattice spacing relies on a Cornell string tension tuned to real QCD and on a Wilson-line/Wilson-loop equivalence that is still 'in preparation' (Section 2.2). And the choice of Scheme II over Scheme I is justified mostly by matching SB at high density, which is reasonable but post hoc.\n\nSo the paper is worth reading by anyone tracking the dense-QCD EoS, and the released data will be reused. But 'clearly breaches' is not supported by the error accounting as it stands. I would not press for the claim to be withdrawn; I would press for the uncertainty in ∂κ/∂a to be propagated and for the matter to be checked on larger volume before anyone builds on the number. For a proceedings paper this level of candor is good; for a central claim it needs another pass. I'd accept this for peer review because the question is important and the group is credible, but with the expectation of revision.","headline":"A plausible result that overreaches on error control: the ∂κ/∂a uncertainty is a load-bearing gap, so 'clearly breaches' is not yet backed.","tokens_in":8896,"tokens_out":2802,"would_cite":false,"duration_ms":25597,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["11.15.Ha","12.38.Gc","12.38.Mh"],"model":"deepseek-v4-flash","headline":"In dense two-color QCD, the squared speed of sound exceeds the conformal limit 1/3.","keywords":["two-color QCD","lattice QCD","speed of sound","equation of state","conformal limit","diquark superfluid","finite-density lattice simulation","static quark potential"],"falsifier":"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.","tokens_in":7622,"feed_emoji":"⚛️","tokens_out":16529,"duration_ms":116229,"temperature":0.7,"pith_summary":"The paper reports lattice simulations of two-color QCD (QC2D) at non-zero baryon density, a sign-problem-free testbed for dense QCD matter. The authors compute the equation of state on a finer lattice than earlier work, with a new scale-setting procedure, and extract the speed of sound. They find that the squared speed of sound $C_s^2$ rises sharply after the onset chemical potential and clearly exceeds the conformal limit $1/3$ over a range of densities. Because the result is consistent with earlier coarser-lattice determinations and with isospin QCD simulations, it supports the idea that the breach of the conformal bound is a robust property of dense strongly interacting matter rather than a lattice artefact.","feed_headline":"Dense two-color QCD's sound speed tops the 1/3 bound","feed_subtitle":"A finer lattice with new scale setting confirms the squared speed of sound rises past 1/3 near onset.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Earlier coarser-lattice QC2D speed-of-sound result that this work extends and confirms at finer spacing.","marker":"[8]"},{"why":"Previous lattice calculation reporting the same conformal-bound breach, the direct predecessor this paper updates.","marker":"[9]"},{"why":"Earlier QC2D result with a different scale-setting procedure, used for consistency comparison.","marker":"[10]"},{"why":"Relativistic-hydrodynamics upper bound on the speed of sound that the computed $C_s^2$ is required to respect.","marker":"[22]"},{"why":"Isospin QCD simulation that also sees $C_s^2$ exceed the conformal bound, cited as corroborating behaviour.","marker":"[12]"},{"why":"Provides the Stefan-Boltzmann correction formulae and quark number density prescription used for the pressure integration.","marker":"[7]"},{"why":"Karsch coefficients used in earlier beta-function determinations, against which the new beta functions are cross-checked.","marker":"[5]"}],"fun_headline_variants":["Two-color QCD sound speed breaches 1/3 bound","Dense QC2D sound speed surpasses conformal limit","Sound speed in two-color QCD tops 1/3","QC2D equation of state: sound speed >1/3 near onset","New lattice calculation: two-color QCD sound speed exceeds 1/3"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Two-color QCD sound speed breaches 1/3 bound","Dense QC2D sound speed surpasses conformal limit","Sound speed in two-color QCD tops 1/3","QC2D equation of state: sound speed >1/3 near onset","New lattice calculation: two-color QCD sound speed exceeds 1/3"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000587,"raw_usage":{"total_tokens":2674,"prompt_tokens":779,"completion_tokens":1895,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":395,"completion_tokens_details":{"reasoning_tokens":1803}},"tokens_in":395,"tokens_out":1895,"duration_ms":11925,"temperature":1.0,"reasoning_tokens":1803,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:00:27.392792+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Speed of sound exceeding the conformal bound in dense 2-color QCD","cited_arxiv_id":"2311.15259","evidence_quote":"Previous lattice calculation reporting the same conformal-bound breach, the direct predecessor this paper updates."},{"cited_title":"Cotter,Non-perturbative determination of Karsch Coefficients in 2 Colour QCD, philosophiæ doctor, National University of Ireland Maynooth, 2015","cited_arxiv_id":null,"evidence_quote":"Karsch coefficients used in earlier beta-function determinations, against which the new beta functions are cross-checked."}],"review_version":1}