{"id":"5a41c4df-9f11-4b5f-898f-a9616705591e","arxiv_id":"2501.01300","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Lattice QCD charm fluctuations show that charmed baryon pressure is about twice the known-state expectation and that low-lying charmed hadrons survive above the chiral crossover before charm quark-like degrees of freedom emerge.","lead":"This lattice QCD study uses charm quark fluctuations to separate the pressure of charmed hadrons from charm quark-like excitations around the QCD crossover. It finds that charmed baryons are much more poorly described by known particle lists than charmed mesons, supporting the existence of missing charmed resonances.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The absolute partial pressures in Fig. 2 are assembled from an unpublished continuum extrapolation of χ4^C and from Nτ=8 ratios whose cutoff independence is asserted, not shown; the headline enhancement factors are not checkable from this preprint.","rationale":"The reader's weakest assumption correctly identifies the unpublished continuum extrapolation of χ4^C and the reliance on Nτ=8 ratios as the main load-bearing point. My stress-test confirms that this is the single most important concern because it affects the first quantitative claim in the abstract and in the strongest_claim: the factors 'almost twice' and 'around 20%'. If the normalization is wrong, those factors shift linearly. The quasi-particle decomposition in Eq. (5) is also model-dependent, but it is secondary: the sub-Tpc baryon enhancement does not require Eq. (5), and the algebra of Eqs. (6)-(8) is internally consistent for C=1 Boltzmann species. The paper also states that the model 'passes numerous validity tests' only by citing prior work, not by showing them here, which reinforces the need for the proposed cross-check. I see no reason to move beyond the reader's CONDITIONAL verdict: the concern is about verifiability and possible systematic error, not a demonstrated internal inconsistency. Therefore the verdict should remain unchanged, conditioned on providing the continuum extrapolation details or numerical data and a lattice-spacing check of the ratios.","tokens_in":9385,"tokens_out":11555,"duration_ms":117454,"concrete_test":"Compute the Nτ=10 (or Nτ=12) value of χ13^BC/χ4^C with LCP[b] at T=156.5 MeV and at T=165 MeV, compare with the Nτ=8 ratios used in Fig. 2, and then recompute P_C^B and P_C^M using the Nτ=10 ratios together with the continuum χ4^C from the forthcoming publication (or from the already published HotQCD/Bazavov et al. 2024 data release). If the resulting enhancement factors 1.948 and 1.215 shift by more than their quoted statistical errors, the unpublished normalization and the cutoff-independence assumption are not robust; if they do not, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative core of the paper, Fig. 2, is built as P_C^B = (χ13^BC/χ4^C)_{Nτ=8} × χ4^C(cont) and P_C^M = (1 − χ13^BC/χ4^C)_{Nτ=8} × χ4^C(cont). The ratios are taken at a single lattice spacing (Nτ=8; no Nτ=10/12 comparison is shown in this manuscript), and the normalization χ4^C(cont) comes from a continuum extrapolation whose details are explicitly deferred to a 'forthcoming publication' (Sec. 2). The headline enhancement factors 1.948 ± 0.234 and 1.215 ± 0.098 are therefore products of two inputs that the reader cannot verify. Because the PDG-HRG partial pressures in the denominator are absolute, any error in the product propagates linearly into both enhancement factors: a 20% error in χ4^C(cont) would move E_B to 2.34 and E_M to 1.46 (or, with the opposite sign, to 1.56 and 0.97), the latter wiping out the claimed 20% meson enhancement. The paper's assertion that cutoff effects 'cancel to a large extent in ratios' is not demonstrated here for χ13^BC/χ4^C; it is carried over from self-cited work. Thus the central quantitative claim, and the meson-sector part in particular, is not checkable from the preprint as it stands.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript analyzes (2+1)-flavor HISQ lattice data for generalized charm susceptibilities chi_13^BC, chi_22^BC, and chi_4^C, using N_tau=8 configurations whose ratios are normalized by continuum-extrapolated chi_4^C values. The continuum extrapolation of chi_4^C is deferred to a separate publication. From these inputs, the authors construct partial pressures of charmed baryons, mesons, and quark-like excitations via Eqs. (6)-(8), and compare the hadronic partial pressures with HRG predictions based on PDG-only and quark-model (QM) spectra, as well as a truncated 1S+1P spectrum. The central quantitative findings are enhancement factors E_B=1.948+/-0.234 and E_M=1.215+/-0.098 at T_pc, evidence that the charmed baryonic sector is substantially more incomplete than the mesonic sector, and a quark-like partial pressure that emerges at T_pc and yields a temperature-dependent in-medium charm quark mass. The paper also argues for sequential melting, with 1S/1P states surviving to about 166 MeV.","tokens_in":9698,"tokens_out":6074,"duration_ms":58515,"significance":"If the result holds, this is the first lattice-QCD-based quantitative determination of partial charm pressures across the chiral crossover, offering direct support for the existence of experimentally unobserved charmed hadrons predicted by quark-model calculations and by the SHMc analysis. A particular strength is that the QM-HRG and PDG-HRG comparisons use external hadron spectra rather than fits to the lattice data, so the missing-resonance enhancement is not circular. The paper also makes a falsifiable prediction that 1S/1P states persist up to about 166 MeV, and it uses publicly available lattice and analysis codes. The main caveats are that the absolute normalization relies on an unpublished continuum extrapolation and the high-temperature decomposition is model-dependent; both need to be assessed before the quantitative claims can be considered established.","major_comments":[{"comment":"The absolute partial pressures are constructed as (chi_13^BC/chi_4^C)_{N_tau=8} and (1 - chi_13^BC/chi_4^C)_{N_tau=8} multiplied by continuum chi_4^C, but the continuum extrapolation is not shown; the text explicitly defers it to 'a forthcoming publication.' Since the PDG-HRG and QM-HRG partial pressures in the denominator are absolute, any error in chi_4^C(cont) propagates linearly into the enhancement factors: for example, a 20% uncertainty would shift E_B from 1.948 to roughly 1.56 or 2.34 and E_M from 1.215 to roughly 0.97 or 1.46, the latter wiping out the claimed meson enhancement. The manuscript needs to present the extrapolation or cite a publicly available result, and to quote a systematic error for the absolute normalization.","section":"Sec. 2 / Fig. 2"},{"comment":"The statement that lattice cutoff effects 'cancel to a large extent' in the ratios is asserted and referenced to Bazavov et al. (2024), but no comparison at another lattice spacing (e.g., N_tau=10 or 12) is shown in this manuscript for chi_13^BC/chi_4^C or the related ratios. Because these ratios are the only lattice input to Fig. 2, the claimed cancellation should be demonstrated in the manuscript or accompanied by an explicit estimate of the residual lattice-spacing uncertainty.","section":"Sec. 2 (cutoff-effect claim)"},{"comment":"The decomposition of the total charm pressure into quark-like, baryon-like, and meson-like partial pressures assumes a specific quasi-particle model inherited from Mukherjee et al. (2016), where each sector has the Boltzmann-form chemical-potential dependence of Eqs. (3)-(4). The resulting P_C^q and m_C^q are not direct lattice observables; they depend on this model and on the assumption that the three sectors jointly saturate the lattice susceptibilities. The text states that prior work passed 'numerous validity tests,' but for the high-temperature claims made here, please provide at least one direct test showing that the model simultaneously describes the three susceptibilities entering Eqs. (6)-(8), or state explicitly which lattice observables the model is not able to reproduce.","section":"Sec. 5, Eqs. (5)-(8)"},{"comment":"The claim that for T_pc < T <= 166.1 MeV both P_C^B and P_C^M are described by 1S1P-HRG is based on visual inspection; no quantitative goodness-of-fit or residual analysis is presented. Given that this comparison is the central evidence for sequential melting, please add a statistical measure (e.g., chi^2/dof over the relevant temperature window) or pointwise pulls for both panels. The same visual-only reasoning is used for the statement that PDG-HRG describes P_C^B at the highest two temperatures.","section":"Sec. 4, Fig. 2 (1S1P-HRG comparison)"}],"minor_comments":[{"comment":"There are several typographical issues, including 'demonstratethatatthechiralcrossover' with missing spaces and 'rather then' for 'rather than'; these should be corrected.","section":"Abstract / Introduction"},{"comment":"The label 'Lattice[b]' is not defined in the caption; it refers to LCP[b] from Sec. 2, but a reader looking only at the figure cannot tell how it differs from LCP[a].","section":"Fig. 2 caption"},{"comment":"The quantity m_C^q is first called the pole mass of the charm quark in the text below Eq. (4), but later it is reinterpreted as a temperature-dependent in-medium quasi-particle mass. Please use consistent terminology and clarify whether Eq. (4) is intended to define a quasi-particle from the outset.","section":"Sec. 3.2 vs. Sec. 5"},{"comment":"The Borsanyi et al. reference is incomplete: the entry gives 'QCD Crossover at Finite Chemical Potential from Lattice Simulations 125, 052001' without a journal name or year. Please check that all references have complete bibliographic information.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The central scientific message is plausible and timely, and the paper fits the journal's scope. However, the absolute normalization of the headline enhancement factors rests on an unpublished continuum extrapolation of chi_4^C, and the high-temperature interpretation rests on a model decomposition that needs direct validation. I recommend a major revision that either includes the extrapolation and a lattice-spacing check for the ratios, or clearly propagates the resulting systematic uncertainty into all quoted enhancement factors. The visual 1S1P-HRG comparison should also be made quantitative. With these additions, the paper could become a solid contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the move from susceptibility ratios to absolute partial pressures, with per-sector enhancement factors at Tpc (baryon ~2x PDG-HRG, meson ~20% over) and the 1S1P-HRG survival analysis above Tpc. The sub-Tpc comparison of lattice data to QM-HRG versus PDG-HRG is the strongest part of the paper and is not circular: the quark-model spectrum is an external input, not a fit to these lattice results. That alone is worth something, because it sharpens the earlier ratio-based evidence for missing charmed baryon states into a quantitative statement about the baryonic sector being roughly half-unaccounted-for in the PDG listing. The algebraic decomposition in Eqs. (6)-(8) is internally consistent, and the authors are honest that the quasi-particle split is a model.\n\nThe soft spot is the one the stress-test flags, and it is real: every absolute number in Fig. 2 is the product of an Ntau=8 ratio (no Ntau=10/12 comparison shown in the manuscript) and a continuum-extrapolated chi4^C whose details are explicitly deferred to a forthcoming publication. The enhancement factors quoted at Tpc therefore inherit whatever systematic error lives in that extrapolation. The stress-test arithmetic is right: a 20% uncertainty in chi4^C(cont) moves the meson enhancement to 1.46 or 0.97, which would erase the claimed 20% effect. The paper asserts that cutoff effects cancel in the ratios, citing earlier work, but it does not demonstrate cancellation for chi13^BC/chi4^C at the needed precision. This is a checkability problem, not a mathematical inconsistency. I also note the comparison to 1S1P-HRG appears visual, with no quoted chi-squared, and the above-Tpc quasi-particle mass is model-dependent by construction - though the authors say so.\n\nProportionate verdict: the sub-Tpc qualitative result (charmed baryon sector dramatically more incomplete than the meson sector) is robust, because it shows up already in the ratio comparison of Fig. 1 without the absolute normalization. The absolute enhancement factors and the sequential-melting story need the missing extrapolation details and preferably a second lattice spacing before I would rely on them.\n\nThis paper deserves a serious referee. The framework is honorable, the new analysis step is real, and the data clearly support the main qualitative claim. I would ask the authors to include the continuum extrapolation of chi4^C (or a clear reference to a companion paper) and to show cutoff-systematics evidence for the specific ratios used, before publication. If they do, I would cite the sub-Tpc result in my own work.","headline":"Solid sub-Tpc analysis of charmed partial pressures, but the headline enhancement factors rest on an unpublished continuum extrapolation that makes the absolute numbers uncheckable from this preprint.","tokens_in":748,"tokens_out":1811,"would_cite":true,"duration_ms":32422,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Lattice QCD charm fluctuations show charmed baryon partial pressure is about twice the PDG hadron gas value at the chiral crossover, with mesons about 20% larger, and that charmed hadrons melt sequentially.","keywords":["lattice QCD","charm fluctuations","generalized susceptibilities","charmed hadrons","hadron resonance gas","chiral crossover","quark-gluon plasma","missing charmed resonances"],"falsifier":"A direct continuum extrapolation of the partial pressures themselves—computing $P_C^B$ and $P_C^M$ at several lattice spacings such as $N_\\tau=6$, 8, 10, and 12 and taking the continuum limit, rather than using $N_\\tau=8$ ratios multiplied by a continuum-extrapolated $\\chi^C_4$—would settle whether the ~1.95 baryon and ~1.22 meson enhancement factors at $T_{pc}$ are real or cutoff artifacts.","tokens_in":9146,"feed_emoji":"⚛️","tokens_out":8329,"duration_ms":71149,"temperature":0.7,"pith_summary":"This paper uses lattice QCD calculations of charm fluctuations and their correlations with baryon number, electric charge, and strangeness to isolate partial pressures from charmed mesons, charmed baryons, and charm quark-like excitations. It seeks to establish that experimentally unobserved charmed hadrons predicted by quark models—mostly baryons—contribute substantially to the charm pressure at and below the chiral crossover, and that open-charm hadrons dissolve sequentially rather than all at once. The paper finds that at the crossover temperature $T_{pc}=(156.5\\pm 1.5)$ MeV the lattice charmed baryon partial pressure is about 1.95 times the value expected from the experimentally known PDG spectrum, while the charmed meson partial pressure is about 1.22 times larger. Above $T_{pc}$ a charm quark-like contribution emerges and grows, while low-lying $1S$ and $1P$ charmed hadron states persist until roughly 166 MeV, supporting a picture of sequential melting.","feed_headline":"Charmed baryons double PDG pressure at the crossover","feed_subtitle":"Quark-model states missing from particle tables are thermodynamically real; charm deconfinement starts at the crossover.","key_machinery":"The machinery is the set of generalized charm susceptibilities $\\chi^{BQSC}_{klmn}$, which are derivatives of the QCD pressure with respect to baryon, electric-charge, strangeness, and charm chemical potentials, computed with the Highly Improved Staggered Quark (HISQ) action on $N_\\tau=8$ lattices. Ratios of these susceptibilities, whose cutoff effects largely cancel, are multiplied by a continuum-extrapolated $\\chi^C_4$ to obtain absolute partial pressures. These are combined with a quasi-particle decomposition $P^C = P^C_M + P^C_B + P^C_q$, where the partial pressures are linear combinations of $\\chi^C_4$, $\\chi^{BC}_{13}$, and $\\chi^{BC}_{22}$, and compared with Boltzmann hadron resonance gas formulas using either the PDG spectrum (PDG-HRG) or quark-model augmented spectra (QM-HRG and 1S1P-HRG).","core_discovery":"The central claim is that continuum-extrapolated lattice QCD results for the partial charm pressure can be decomposed, through generalized susceptibilities, into charmed-meson, charmed-baryon, and charm quark-like contributions. Below and at the chiral crossover this decomposition shows that the baryonic sector is almost twice as large as the PDG hadron spectrum alone predicts, while the mesonic sector is about 20% larger, with both matching quark-model augmented hadron resonance gas predictions. At the crossover the hadron resonance gas description breaks down, signaling the onset of charm deconfinement, yet the low-lying $1S$ and $1P$ charmed hadron states survive to about 166 MeV, indicating sequential melting. The emerging quark-like partial pressure is used to extract a temperature-dependent in-medium charm quark mass that starts near the $D$-meson mass and decreases with temperature.","pith_inferences":["Extending the same decomposition to higher temperatures would predict that above roughly two times $T_{pc}$ the charm partial pressure becomes dominated by quark-like excitations and hadron-like contributions vanish, which can be tested once continuum-extrapolated susceptibilities become available at those temperatures.","If the missing charmed baryons are thermodynamically real, heavy-ion yields of charmed baryons should show a corresponding enhancement across collision energies; a direct comparison of these lattice partial pressures with measured $D$-meson and charmed-baryon yields would test this connection.","The near-equality of the baryonic enhancement factor (1.95) with the QM-HRG prediction suggests the quark-model spectrum is nearly complete thermodynamically at $T_{pc}$; future data at smaller lattice spacings that reproduce the same ratio would harden the missing-resonance interpretation into a quantitative prediction for spectroscopy."],"forward_implications":["The charmed baryon sector of the PDG hadron list is incomplete at the level of roughly a factor of two in partial pressure at $T_{pc}$, so quark-model predicted charmed baryons, not mesons, are the dominant missing contribution.","Open-charm hadrons begin to dissolve at the chiral crossover, but the dissociation is sequential: low-lying $1S$ and $1P$ charmed states survive to about 166 MeV while higher excitations melt first.","A charm quark-like contribution to the partial pressure appears at $T_{pc}$ and grows, overtaking hadron-like contributions near $T\\sim 175$ MeV, so charmed hadron-like excitations can persist well into the quark-gluon plasma regime.","The temperature-dependent in-medium mass of the charm quark-like excitation, starting near the $D$-meson mass around 162 MeV and decreasing with temperature, gives a quantitative handle on charm in-medium interactions in the quark-gluon plasma.","Charmed hadron masses are not strongly affected by chiral symmetry restoration at $T_{pc}$, since the $1S1P$-HRG prediction does not overshoot the lattice meson pressure above the crossover."],"supporting_citations":[{"why":"Supplies the lattice ratios $P_C^B/P_C$ and $P_C^M/P_C$ on $N_\\tau=8$ configurations and the quasi-particle model validity tests that this paper converts into absolute partial pressures.","marker":"Bazavov et al. (2024)"},{"why":"Provides the (2+1)-flavor HISQ gauge configurations, the $f_K$ temperature scale, and the $m_s/m_l=27$ setup on which all susceptibilities are computed.","marker":"Bollweg et al. (2021)"},{"why":"Introduces the HISQ action with an epsilon term that removes leading $(am_c)^4$ lattice artifacts, underpinning the claim that susceptibility ratios are nearly continuum.","marker":"Follana et al. (2007)"},{"why":"Supplies the Boltzmann-approximation formulas that relate partial hadronic pressures to masses and degeneracies, and the ideal charm quark gas expression used in the decomposition.","marker":"Allton et al. (2005)"},{"why":"Provides the quasi-particle model of Eq. (5) that separates the charm pressure into meson, baryon, and quark-like partial pressures.","marker":"Mukherjee, Petreczky and Sharma (2016)"},{"why":"Tabulates the quark-model predicted charmed hadrons, including the low-lying $1S$ and $1P$ states used for QM-HRG and 1S1P-HRG comparisons.","marker":"Chen, Chen, Liu, Liu and Zhu (2023)"},{"why":"Provides the quark-model spectroscopy for charmed mesons and baryons that defines the enhanced HRG spectrum whose predictions the lattice data match.","marker":"Ebert, Faustov and Galkin (2010, 2011)"},{"why":"Gives the PDG masses and the charmed-to-strange mass ratio $m_c/m_s=11.76$ used for charm quark mass tuning and as the baseline PDG-HRG spectrum.","marker":"Workman et al. (2022)"},{"why":"Sets the chiral crossover temperature $T_{pc}=(156.5\\pm1.5)$ MeV used as the reference point for all enhancement factors and melting statements.","marker":"Bazavov et al. (2019)"}],"fun_headline_variants":["Quark-model charm states double baryon pressure","Charm deconfinement starts at chiral crossover","Lattice QCD reveals sequential melting of charmed hadrons","Missing charmed hadrons are thermodynamically real","In-medium charm quark mass from lattice QCD"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central numbers stand on the assumption that the $N_\\tau=8$ lattice ratios are effectively continuum and that the unpublished continuum extrapolation of $\\chi^C_4$ is accurate; if either fails, the enhancement factors and all trends above $T_{pc}$ shift.","fun_headline_variants_meta":{"raw":{"variants":["Quark-model charm states double baryon pressure","Charm deconfinement starts at chiral crossover","Lattice QCD reveals sequential melting of charmed hadrons","Missing charmed hadrons are thermodynamically real","In-medium charm quark mass from lattice QCD"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000764,"raw_usage":{"total_tokens":3365,"prompt_tokens":893,"completion_tokens":2472,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":509,"completion_tokens_details":{"reasoning_tokens":2395}},"tokens_in":509,"tokens_out":2472,"duration_ms":17104,"temperature":1.0,"reasoning_tokens":2395,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:31:11.968031+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct continuum extrapolation of the partial pressures themselves—computing $P_C^B$ and $P_C^M$ at several lattice spacings such as $N_\\tau=6$, 8, 10, and 12 and taking the continuum limit, rather than using $N_\\tau=8$ ratios multiplied by a continuum-extrapolated $\\chi^C_4$—would settle whether the ~1.95 baryon and ~1.22 meson enhancement factors at $T_{pc}$ are real or cutoff artifacts.","supporting_citations":[{"cited_title":"Charmdegreesoffreedom inhotmatterfromlatticeQCD","cited_arxiv_id":null,"evidence_quote":"Supplies the lattice ratios $P_C^B/P_C$ and $P_C^M/P_C$ on $N_\\tau=8$ configurations and the quasi-particle model validity tests that this paper converts into absolute partial pressures."}],"review_version":1}