{"id":"51b6d1f9-ba20-44e3-a260-6e7da56637ee","arxiv_id":"2505.01734","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Charm fluctuations from lattice QCD match a hadron resonance gas only when quark-model-predicted charmed hadrons are added, and above the crossover a charm-quark-like partial pressure emerges.","lead":"Lattice QCD calculations of charm fluctuations show that below the QCD crossover temperature charm behaves like a gas of hadrons, including many charmed particles not yet discovered, while just above it a charm-quark-like component appears alongside hadron-like excitations. The continuum estimates give quantitative reference points for interpreting charmed-hadron measurements from heavy-ion collisions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim rests on the unvalidated quasi-particle projection in Eq. (15): it vanishes for all ordinary charmed hadrons by construction, so only the self-cited three-operator test of Ref. [22] supports interpreting its non-zero value above Tpc as a charm quark.","rationale":"The reader's weakest_assumption is exactly the load-bearing point. The lattice work itself is careful: multiple lattice spacings, two LCPs, bootstrap errors, Nτ=8/12 cross-checks of ratios, and a public data repository. The QM-HRG-vs-PDG-HRG observations below Tpc, and the breakdown of QM-HRG in selected ratios above Tpc, are well supported and would survive even if the quasi-particle interpretation were abandoned. What does not survive is the abstract's specific phrase 'charm quark pressure becomes non-zero': Eq. (15) is a linear combination of generalized susceptibilities, and its vanishing below Tpc is automatic for any mixture of B=0 mesons and B=1 baryons, so it cannot validate the operator. The paper explicitly says the only support for the projection is the consistency test in the authors' previous Ref. [22], which is self-cited and not included here. There is also a genuine quantitative gap: B=1/3 and B=2/3 carriers contribute identically to the combination χ_13^{BC} − χ_22^{BC}, so the operator does not uniquely select a charm quark. These considerations do not justify rejection because the paper is transparent about the model dependence and the observational core (missing charmed hadrons, QM-HRG breakdown) is solid; they do justify keeping the verdict CONDITIONAL rather than upgrading it. I therefore agree with the reader and leave the verdict unchanged.","tokens_in":22472,"tokens_out":9183,"duration_ms":97948,"concrete_test":"Using the published dataset [33], recompute P_C^q on the Nτ=12 and Nτ=16 configurations from Eq. (15) and from the two alternative projections defined in Ref. [22] (operators with |Q|=2/3,|C|=1 and |B|=1/3,|Q|=2/3,|C|=1) over the temperature range 157-177 MeV, with jackknife errors and identical LCP[b] masses. If the three operators do not agree within errors above Tpc, the quasi-particle interpretation is not unique and the central claim fails. If they do agree on the finer lattices, the self-cited consistency check is independently confirmed and the concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline claim that the charm quark pressure becomes non-zero above Tpc is not a direct lattice measurement; it is the output of the quasi-particle decomposition in Eq. (14). The projection Eq. (15), P_C^q = 9(χ_13^{BC} − χ_22^{BC})/2, is constructed to give zero for any C=1 carrier with B=0 or B=1, so its vanishing below Tpc is guaranteed in any HRG-like phase and carries no independent evidence. Its non-zero value above Tpc identifies a B=1/3, C=1 excitation only if every other contribution has B=0 or B=1 and the gas is non-interacting. The operator cannot even distinguish B=1/3 from B=2/3 carriers, since B−B^2 = 2/9 for both. The only direct validation cited is the three-operator consistency test in the authors' previous work [22], a self-cited result not reproduced in this paper; the lattice data in Sec. VI.A are the same susceptibilities whose HRG deviation is the evidence. The continuum estimate additionally multiplies Nτ=8 ratios from [22] by the Nτ=16[D] continuum χC4 rather than continuum-extrapolating Eq. (15) itself. If hadronic or B=2/3 contributions overlap with this projection, the abstract's 'charm quark pressure' is not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a lattice QCD study of generalized charm susceptibilities in 2+1 flavor QCD with quenched charm, using HISQ configurations at Nτ = 8, 12, and 16. A new line of constant physics (LCP[D]) is constructed by requiring the physical D-meson mass, which reduces the dominant cutoff effects in the charm sector. The authors provide a continuum estimate of the quartic charm fluctuation χ_4^C and study ratios of baryon-charm, charge-charm, and strangeness-charm susceptibilities. They show that below T_pc the results are described by a quark-model based hadron resonance gas (QM-HRG) that includes missing charmed states, while above T_pc the QM-HRG description breaks down for several observables. Using a quasi-particle decomposition proposed in Refs. [11] and [22], the paper converts the susceptibility ratios into partial pressures of charmed mesons, charmed baryons, and charm quarks, and reports that the charm quark partial pressure becomes nonzero above the chiral crossover, with a temperature-dependent in-medium mass m_C^q(T). The paper also decomposes charmed pressures by strangeness and discusses the sensitivity of χ_22^{SC} to the spectrum of doubly-strange charmed baryons.","tokens_in":22754,"tokens_out":5842,"duration_ms":59660,"significance":"If the central interpretation is correct, the paper provides important evidence that charm thermodynamics across the QCD crossover is described by coexisting charmed hadron-like and quark-like excitations, with a charmed baryon spectrum substantially richer than the PDG tables. The work has several concrete strengths: the LCP[D] construction demonstrably reduces a large source of cutoff effect; the Nτ = 8[b] and Nτ = 12[b] ratios agree where both are available, supporting the ratio method; errors are propagated with bootstrap procedures; and all data are publicly released. The paper is also appropriately cautious in several places, explicitly labeling the HRG description of χ_4^C above T_pc as accidental and the 1S1P-HRG comparison as rough guidance. However, the headline result is not a direct lattice measurement: the nonzero charm quark pressure is an output of a model decomposition whose central operator, Eq. (15), has not been validated in the present paper beyond a self-cited consistency test from Ref. [22]. The significance of the paper therefore depends on whether this model dependence is made fully explicit and whether the projection is supported by independent tests.","major_comments":[{"comment":"The operator P_C^q = 9(χ_13^{BC} − χ_22^{BC})/2 vanishes by construction for any non-interacting species with C = 1 and B = 0 or B = 1, so the observed vanishing below T_pc is a built-in property of the HRG-like assignment and does not by itself provide independent evidence for the absence of quark-like excitations. Moreover, the same operator yields identical values for B = 1/3 and B = 2/3 carriers because B − B^2 = 2/9 in both cases. Interpreting P_C^q as a charm quark partial pressure therefore requires the additional assumptions that the system is a non-interacting mixture of mesons, baryons, and B = 1/3 quarks, and that no charmed hadron contributes to the combination. The only direct validation cited is the three-operator consistency test in the authors' previous work [22], which is not reproduced or described in sufficient detail here. Unless this test is shown or an independent validation is provided, the abstract's statement that 'the charm quark pressure becomes non-zero above the chiral crossover' is not established by the present lattice data.","section":"Sec. VI.A, Eq. (15)"},{"comment":"The continuum estimates of P_C^M, P_C^B, and P_C^q shown in Figs. 12 and 13 are formed by multiplying the Nτ = 8[b] normalized ratios from Ref. [22] by the Nτ = 16[D] continuum estimate of χ_4^C, rather than by continuum-extrapolating the combinations in Eqs. (15)–(17) themselves. Cutoff cancellation is demonstrated for the ratios, but the product with a separately continuum-extrapolated overall normalization has no documented systematic uncertainty. A direct continuum extrapolation of P_C^q/χ_4^C (or of P_C^q itself) using the Nτ = 8 and Nτ = 12 data presented in this paper would be needed to support the quoted absolute values and the associated enhancement factors.","section":"Sec. VI.A, Figs. 12 and 13"},{"comment":"The decomposition P_C = P_M + P_B + P_q assumes that the three channels are non-interacting and independent, with all in-medium effects encoded in the temperature-dependent mass m_C^q in Eq. (10). If interactions mix the channels, the linear relations (15)–(17) do not define physical partial pressures. The manuscript should state this limitation explicitly in the abstract or conclusions, or soften the headline claim from 'the charm quark pressure becomes non-zero' to 'the charm-quark-like projection of the susceptibilities becomes non-zero in the quasi-particle model'.","section":"Sec. VI.A, Eq. (14)"}],"minor_comments":[{"comment":"In the sentence following Eq. (13), the symbol P_{C,S=2}^{M} should be P_{C,S=2}^{B}, since the text is describing the partial pressure of strange charmed baryons with strangeness two, not mesons.","section":"Eq. (13)"},{"comment":"The caption of Fig. 14 states 'Solid: QM-HRG, Dotted: PDG-HRG, Dashed: 1S1P-HRG', while the main text and the figure keys elsewhere describe 'Dashed: QM-HRG, Dotted: PDG-HRG, Solid: 1S1P-HRG'. The line styles should be made consistent.","section":"Fig. 14 caption"},{"comment":"The phrase 'quarks-antiquarks' in the note about m_C^q differing from Ref. [29] should read 'quarks and antiquarks'.","section":"Sec. VI.A, Fig. 13 caption"},{"comment":"The sentence 'the QM-HRG description breaks down just above T_pc signaling the appearance of new degrees of freedom' uses stronger language than the earlier 'signals the possible appearance'; since the quasi-particle interpretation is model-dependent, 'possibly signaling' would be more accurate.","section":"Sec. V, concluding paragraph"}],"recommendation":"major_revision","confidential_remarks":"The central result is model-dependent, and the key validation is a self-cited consistency test from Ref. [22] that is not reproduced here. The editors may wish to ensure the paper is read as a lattice study plus a quasi-particle interpretation rather than as a direct measurement of a charm quark pressure. No concerns about data availability or attribution; the data release and the careful LCP[D] construction are commendable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take on Kaczmarek et al., arXiv:2505.01734. The lattice work is careful and the new numbers are useful; the headline \"charm quark pressure becomes non-zero above Tpc\" is a model output, not a direct measurement, and the abstract overstates it.\n\nWhat's actually new: the LCP[D] construction (fixing the D-meson mass rather than charmonium) is a genuine step forward — it visibly reduces the dominant cutoff effect in χC4, and the Nτ=8 versus Nτ=12 agreement of the ratios justifies using ratios to cancel cutoffs. The continuum χC4 from Nτ=16[D] is a quantity the community will use. Data are in a public repository. The paper is also candid in the right places: it flags that χC4's agreement with HRG above Tpc is accidental, that the 1S1P-HRG comparison is only rough guidance, and that the QM-HRG underprediction of χSC22/χC4 may just mean the |S|=2 spectrum is wrong.\n\nSoft spots, in order of size. First, Eq. (15): the operator P_C^q = 9(χ_13^BC − χ_22^BC)/2 vanishes for any B=0 or B=1 carrier by construction, so its vanishing below Tpc is guaranteed and carries no evidence. What the lattice shows is that a B−B²-weighted combination rises above Tpc. The identification of that with a charm quark rests on the three-operator consistency test in the authors' own Ref [22], not reproduced here; plausible, but the headlines are built on it. The stress-test note is also right that the operator cannot distinguish B=1/3 from B=2/3 carriers, since both give B−B²=2/9. I would not call this circular — the model interprets the deviation rather than generating it — but it is model-dependent, and the abstract's \"charm quark pressure\" goes beyond what the operator alone establishes.\n\nSecond, the \"continuum estimates\" of partial pressures are Nτ=8 ratios multiplied by the Nτ=16[D] continuum χC4. That is reasonable if the ratios are cutoff-free — and the Nτ=12 checks support that — but it is not the same as continuum-extrapolating the partial pressures, and the residual systematic is unquoted. Third and minor: the conclusions contain a sentence saying results above Tpc are \"well described by\" QM-HRG a paragraph before saying QM-HRG breaks down above Tpc; that contradicts the abstract and needs fixing. The 100 MeV downward shift of |S|=2 baryons is hand-adjusted, but the authors present it as an estimate and Fig. 15 shows it does the job.\n\nBottom line: the qualitative story — QM-HRG works below Tpc only with missing charmed states, and specific ratios deviate from it above — is well supported by the lattice data. The quark-like interpretation is plausible but should be framed as inference. The paper is for lattice practitioners, heavy-ion phenomenologists, and anyone building HRG or quasi-particle models of charm thermodynamics. It deserves a serious referee, and would be citable after a revision that reframes the headline, reproduces or forwards the three-operator test, and fixes the conclusions paragraph.","headline":"The lattice data are careful and the continuum χC4 is usable, but the 'charm quark pressure' headline is a quasi-particle model output, not a measurement — worth publishing after revision.","tokens_in":23441,"tokens_out":8571,"would_cite":true,"duration_ms":75317,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["11.10.Wx","11.15.Ha","12.38.Aw","12.38.Gc","12.38.Mh","24.60.Ky","25.75.Gz","25.75.Nq"],"model":"deepseek-v4-flash","headline":"Charm thermodynamics across the QCD crossover splits into hadron-like and quark-like pieces.","keywords":["charm susceptibilities","chiral crossover","hadron resonance gas","lattice QCD","quasi-particle model","open charm hadrons","generalized susceptibilities","charmed baryons"],"falsifier":"Compute the spectral content of the operator $\\chi^{BC}_{13} - \\chi^{BC}_{22}$ above $T_{pc}$: if lattice correlation functions built from these current combinations show no $|B| = 1/3$ excitation peak that survives the continuum limit — or if charmed-hadron states with $|B| = 1$ contribute to this combination — then the claimed charm quark partial pressure is an artifact of the projection. A simpler complementary test: measure the isolated meson and baryon partial pressures independently from hadronic correlation functions and check that they sum with $P_q^C$ to the total charm pressure.","tokens_in":22110,"feed_emoji":"⚛️","tokens_out":5891,"duration_ms":53665,"temperature":0.7,"pith_summary":"This paper uses lattice QCD to track how charm degrees of freedom behave as matter heats through the chiral crossover temperature $T_{pc} \\approx 156.5$ MeV. It establishes that a hadron resonance gas description that includes quark-model-predicted charmed hadrons works below $T_{pc}$ but breaks down just above it. It then decomposes the total charm pressure into partial pressures of charmed mesons, charmed baryons, and charm quarks, and finds that the charm quark partial pressure, zero below the crossover, becomes nonzero above it. On the way, it shows that charmed hadron tables are incomplete: at $T_{pc}$ only 51(6)% of the charmed baryon pressure and 88(7)% of the charmed meson pressure come from experimentally known states. If correct, charm matter across the crossover is a coexistence of hadron-like and quark-like excitations.","feed_headline":"Charm quark pressure switches on above the QCD crossover","feed_subtitle":"Lattice QCD: at Tpc, known charmed baryons supply just 51(6)% of the pressure; missing states and a quark-like part carry the rest.","key_machinery":"The basic identity is the decomposition of total charm pressure into three partial pressures, $P^C = P_M^C + P_B^C + P_q^C$, with the quark piece projected by $P_q^C = 9(\\chi^{BC}_{13} - \\chi^{BC}_{22})/2$, an operator carrying $|B| = 1/3$ and $|C| = 1$. These are expressed as combinations of generalized charm susceptibilities $\\chi^{BC}_{mn}$, and the continuum limit is reached by computing the quartic charm fluctuation $\\chi^C_4$ on $N_\\tau = 8, 12, 16$ lattices using a line of constant physics fixed by the $D$-meson mass. The hadronic baseline is QM-HRG, a hadron resonance gas that supplements PDG states with quark-model-predicted charmed hadrons; it is the benchmark against which the breakdown above $T_{pc}$ is measured.","core_discovery":"The central claim is that the QCD crossover in the charm sector is not a single transition from hadrons to free quarks but a gradual coexistence: below about $T_{pc}$ everything is charmed hadrons, while just above $T_{pc}$ a charm-quark-like excitation with baryon number $1/3$ and charm $1$ carries a nonzero share of the pressure, coexisting with charmed meson and baryon excitations that still dominate up to about 176 MeV. Quantitatively, continuum-limit charmed meson and baryon partial pressures are enhanced by factors $1.13(9)$ and $1.95(23)$ relative to PDG-based HRG, meaning half the charmed baryon pressure comes from states not in the tables. The quasi-particle model of Ref. [11], with pressures built from generalized susceptibilities, is the interpretive frame: charmed hadron pressures fall below quark-model HRG predictions above $T_{pc}$, while the extracted charm quark pressure rises from zero.","pith_inferences":["A sharp way to test the quark-like assignment is to compute, on the lattice, the overlap of the operator in Eq. (15) with charmed hadron states; if that overlap is nonzero, some of the 'quark pressure' is actually hadron pressure relabeled.","If the coexistence picture is right, charm-flow observables in heavy-ion collisions should show a gradual onset of quark transport rather than an abrupt appearance.","The same susceptibility-ratio technology could be applied to bottom quarks, whose larger mass should push the analogous crossover to higher temperatures.","Because $\\chi^C_4$ agrees with QM-HRG even above $T_{pc}$ while some ratios do not, total charm yield alone is a poor probe of deconfinement; only charge- and strangeness-resolved correlations separate the degrees of freedom."],"forward_implications":["At $T_{pc}$ the charmed baryon spectrum known to experiment supplies only about half the charmed baryon pressure; the other half must be composed of charmed baryon resonances predicted by quark models but not yet observed.","Above the crossover the QM-HRG description fails, so heavy-ion phenomenology must not assume a purely hadronic charm yield even a few MeV above $T_{pc}$.","The nonzero charm quark partial pressure above $T_{pc}$ gives a concrete quasi-particle mass for charm that drops with temperature.","The ratio $\\chi^{SC}_{22}/\\chi^C_4$ is sensitive to doubly-strange charmed baryons; matching it may require shifting their masses down by about 100 MeV compared to current quark-model predictions."],"supporting_citations":[{"why":"Introduced the quasi-particle decomposition of charm pressure into meson, baryon, and quark partial pressures and supplied the operator constructions used here.","marker":"[11]"},{"why":"High-statistics $N_\\tau = 8$ calculation that established the three-operator consistency test and gave the fractional partial pressures converted to absolute values in this paper.","marker":"[22]"},{"why":"Earlier lattice study of generalized charm susceptibilities that first showed hadron-gas behavior below $T_{pc}$ and new degrees of freedom above; provides high-temperature $N_\\tau = 8$ data used here.","marker":"[10]"},{"why":"Supplies the (2+1)-flavor HISQ gauge configurations, the temperature scale, and the $m_s/m_l = 27$ line of constant physics underlying all lattice calculations.","marker":"[17]"},{"why":"Relativistic quark model predictions for the masses of not-yet-observed charmed hadrons used to build the QM-HRG particle list.","marker":"[37]"},{"why":"Particle Data Group compilation of experimentally known charmed hadrons defines the PDG-HRG baseline and the experimentally known fraction of charm pressure.","marker":"[35]"},{"why":"Supports treating the charm quark in the quenched approximation by showing dynamical charm effects become significant only above about 300 MeV.","marker":"[45]"},{"why":"Lattice studies of charmed meson and baryon correlation functions that support the existence of hadron-like excitations above $T_{pc}$ and in-medium mass modification.","marker":"[23, 24]"}],"fun_headline_variants":["Charm quark pressure emerges at the QCD crossover","Half of charmed baryon pressure from unlisted states","Charm quarks share pressure with hadrons above QCD crossover","Lattice QCD reveals a hadron-quark charm coexistence"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the operator $P_q^C = 9(\\chi^{BC}_{13} - \\chi^{BC}_{22})/2$ isolates a charm-quark-like excitation with $|B| = 1/3$ and $|C| = 1$ and receives no contribution from charmed hadron states; if hadronic states overlap the operator, the nonzero 'quark pressure' above $T_{pc}$ would not establish a new quark-like degree of freedom.","fun_headline_variants_meta":{"raw":{"variants":["Charm quark pressure emerges at the QCD crossover","Half of charmed baryon pressure from unlisted states","Charm quarks share pressure with hadrons above QCD crossover","Lattice QCD reveals a hadron-quark charm coexistence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000339,"raw_usage":{"total_tokens":1851,"prompt_tokens":905,"completion_tokens":946,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":521,"completion_tokens_details":{"reasoning_tokens":877}},"tokens_in":521,"tokens_out":946,"duration_ms":9261,"temperature":1.0,"reasoning_tokens":877,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:14:31.252317+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the spectral content of the operator $\\chi^{BC}_{13} - \\chi^{BC}_{22}$ above $T_{pc}$: if lattice correlation functions built from these current combinations show no $|B| = 1/3$ excitation peak that survives the continuum limit — or if charmed-hadron states with $|B| = 1$ contribute to this combination — then the claimed charm quark partial pressure is an artifact of the projection. A simpler complementary test: measure the isolated meson and baryon partial pressures independently from hadronic correlation functions and check that they sum with $P_q^C$ to the total charm pressure.","supporting_citations":[{"cited_title":"The melting and abundance of open charm hadrons","cited_arxiv_id":"1404.4043","evidence_quote":"Earlier lattice study of generalized charm susceptibilities that first showed hadron-gas behavior below $T_{pc}$ and new degrees of freedom above; provides high-temperature $N_\\tau = 8$ data used here."},{"cited_title":"Kaczmarek, F","cited_arxiv_id":null,"evidence_quote":"Relativistic quark model predictions for the masses of not-yet-observed charmed hadrons used to build the QM-HRG particle list."},{"cited_title":"For instance, we added the unobserved D∗ + in addition to the observed D∗","cited_arxiv_id":null,"evidence_quote":"Particle Data Group compilation of experimentally known charmed hadrons defines the PDG-HRG baseline and the experimentally known fraction of charm pressure."}],"review_version":1}