{"id":"6254ef6e-24ca-4f26-b69b-ab5f609ad0a2","arxiv_id":"2506.08589","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Thermal QCD sum rules predict that the X(6200) state, treated as an eta_c eta_c molecule, loses about 7% of its mass and 62% of its decay constant as the temperature rises to 0.14 GeV.","lead":"This paper calculates how the mass and coupling of the newly discovered X(6200) particle, assumed to be a pair of eta_c mesons, change with temperature. It predicts both decrease as the temperature approaches the quark-gluon plasma transition, which could guide searches in heavy-ion collisions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thermal prediction is imprinted by the ad hoc s(T) ansatz of Eq. (19); without independent constraints on the continuum threshold the 7% mass drop and 62% decay-constant drop are not established.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: Eq. (19) is ad hoc and the thermal drop is essentially encoded in it. My reading of the manuscript confirms this: the OPE input is standard (dimension-4 condensates, thermal coupling), the T=0 sum rule reproduces the X(6200) mass reasonably, but the temperature evolution is governed by the continuum threshold model. The paper even states that heavy-heavy systems 'deviate significantly' from the light-quark threshold treatments cited, and then provides no derivation for the (T/Tc)^8 form. This is a correctness risk rather than an internal inconsistency: under the stated assumptions the calculation is coherent, so the verdict should remain CONDITIONAL rather than REJECT. The concrete test I propose is the minimal decisive check: vary the exponent or freeze s(T), and quantify how much of the drop is ansatz-driven. That directly settles whether the quantitative claim (7% mass drop, 62% decay-constant drop) is a prediction or a parametrization. I do not treat the closing beyond-Standard-Model remark as load-bearing for the central claim, though it is unsupported and should be removed. I also note the paper provides full analytical expressions for the spectral densities, which is useful evidence that the T=0 calculation is reproducible in principle, but no code or machine-checked proof is provided, and the parameter count is modest (mc, condensate, s0, M^2 window), so novelty is in the application rather than the method.","tokens_in":13435,"tokens_out":3107,"duration_ms":27345,"concrete_test":"Recompute m(T) and f(T) for, e.g., s(T) = s0(1-(T/Tc)^4) + 16mc^2(T/Tc)^4 or a linear interpolation, keeping all other inputs and the Borel window fixed, and compare the values at T=0.14 GeV with Table 3. Also run the same calculation with the thermal threshold held constant (s(T)=s0) to isolate how much of the drop comes from the s(T) ansatz versus the thermal OPE and coupling. If the mass drop changes from ~7% to <2% or the decay-constant drop changes from ~62% to <20%, the headline thermal prediction is dominated by the unconstrained ansatz and the claim should be reframed as model-dependent.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central thermal claim is that m(T) and f(T) decrease substantially as T approaches Tc. Everything thermal in the sum rules flows through Eq. (19), s(T) = s0(1-(T/Tc)^8) + 16mc^2(T/Tc)^8. The exponent 8 and the endpoint 16mc^2 are asserted without derivation or citation to a first-principles result. The integrals in Eqs. (13)-(15) are cut at s(T), so lowering s(T) toward the partonic threshold 16mc^2 is exactly what pushes mass and decay constant down; this is visible in Figure 3, where smaller s0 at fixed T=0 already lowers m and f noticeably. Thus the predicted 7% mass drop and 62% decay-constant drop are largely a restatement of the assumed threshold evolution rather than an independent prediction of the OPE. The only cited justification (Refs. [46,47]) is for light-light and heavy-light systems; the paper itself acknowledges heavy-heavy thresholds deviate and need separate treatment, but no such treatment is given. Additionally, the quoted uncertainties in Table 3 are not propagated through the temperature dependence: Table 3 shows errors at T=0 and T=0.14, yet the temperature plots and the 7%/62% figures are single curves, so the statistical significance of the thermal drop is not quantified. A secondary but related concern is that the Borel window M^2=5-6.5 GeV^2 appears large relative to the lowest threshold 16mc^2≈25.8 GeV^2; with m=6.2 GeV and these M^2 values the exponential weight e^{-s/M^2} is steep but the pole contribution is stated to fall from 86% to 54% only by choosing s0≈44-45 GeV^2, which is far above the ηcηc threshold at ~35.6 GeV^2; the physical interpretation of the continuum in this window is not defended. The qualitative direction (thermal suppression) is plausible and consistent with similar TQCDSR studies, but the quantitative claim rests on an unvalidated ansatz.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the fully-charmed scalar state X(6200) observed by LHCb, assuming it is an eta_c eta_c molecular state with J^{PC}=0^{++}, and computes its mass and decay constant as functions of temperature using thermal QCD sum rules (TQCDSR) up to dimension four. The authors derive the correlation function with a di-pseudoscalar interpolating current, perform a Borel transformation, and apply a temperature-dependent continuum threshold. Their central result is that both the mass and the decay constant decrease as temperature approaches T_c = 155 MeV, with a benchmark at T = 0.14 GeV giving a mass of 5782.43 ± 37.42 MeV (about 7% below the T=0 value of 6201.74 ± 57.82 MeV) and a decay constant of 0.88 ± 0.05 × 10^{-2} GeV^4 (about 62% below the T=0 value of 2.31 ± 0.18 × 10^{-2} GeV^4). The T=0 results are claimed to be consistent with previous QCD sum rule literature.","tokens_in":13911,"tokens_out":4489,"duration_ms":54451,"significance":"If established, a robust prediction of the thermal melting of a fully-heavy molecular state would be of interest for heavy-ion phenomenology and for understanding the behavior of exotic hadrons in hot QCD matter. The paper has notable strengths: it provides explicit analytic expressions for the spectral densities and dimension-four contributions in Appendix A without truncation, which is valuable for reproducibility and independent checks; it also includes a T=0 benchmark that matches an earlier QCD sum rule result for the same state. However, the central thermal claim rests on an assumed parametrization of the continuum threshold, Eq. (19), whose functional form and endpoint are not derived or independently validated. The manuscript itself acknowledges that heavy-heavy systems require a separate treatment, yet applies the light-quark motivated ansatz anyway. Because the Borel integrals in Eqs. (13)–(15) are truncated at s(T), the predicted decreases are largely imprinted by this ansatz rather than by the OPE. The quoted uncertainties are also not propagated through the temperature dependence.","major_comments":[{"comment":"The temperature-dependent continuum threshold s(T) = s0(1-(T/Tc)^8) + 16mc^2(T/Tc)^8 is the key input that drives the thermal behavior, since the Borel integrals in Eqs. (13)–(15) are cut at s(T). The exponent 8 and the endpoint 16mc^2 are asserted without derivation or independent constraint. More importantly, the manuscript itself states, immediately after Eq. (19), that for heavy-heavy quark systems this behavior deviates significantly and necessitates a separate treatment as in Refs. [46,47], but no such separate treatment is given. As a consequence, the predicted 7% mass drop and 62% decay-constant drop are largely a restatement of the assumed threshold trajectory rather than a prediction of the OPE. The authors should either derive s(T) from first principles (e.g., from the thermal behavior of the relevant condensates), or treat the functional form as a model parameter and demonstrate the sensitivity of the results to its exponent and endpoint.","section":"Section 2, Eq. (19) and the paragraph following it"},{"comment":"The uncertainties quoted in Table 3 are not propagated through the temperature dependence. Figure 4 shows only single curves for m(T) and f(T) with no error bands, and the 7% and 62% changes are quoted without any uncertainty. Furthermore, the reported error at T = 0.14 GeV (37.42 MeV) is smaller than the error at T = 0 (57.82 MeV), which is counterintuitive unless the sources of uncertainty are not propagated through the thermal inputs. The authors should describe exactly how the errors in Table 3 were obtained, include all sources (mc, condensate value, s0 window, and the parameters of Eq. (19)), and present confidence bands for the thermal curves.","section":"Table 3 and Figure 4"},{"comment":"The Borel window M^2 = 5–6.5 GeV^2 is relatively low compared with the partonic threshold 16mc^2 ≈ 25.8 GeV^2, and the manuscript reports a pole contribution that falls from 86% to 54% over the analysis window. It is not stated whether this pole-contribution range refers to the T=0 M^2 window or to the entire temperature range used in Figure 4. If the 54% lower value applies at T = 0.14 GeV, then at the highest temperature the continuum contribution is comparable to the ground-state contribution, which substantially weakens the reliability of the extracted m(T) and f(T). The authors should specify the pole-contribution range as a function of temperature and either restrict the Borel window to maintain a robust pole dominance or discuss the implications of large continuum contamination.","section":"Section 2, pole contribution and Table 1"}],"minor_comments":[{"comment":"There are several typographical errors, including 'guarentee' (Section 2), 'condansates' (Section 2), and 'expicitly' (Section 2). These should be corrected.","section":"Throughout"},{"comment":"The manuscript cites Refs. [46,47] for the temperature dependence of s(T), but those references concern light-light and heavy-light systems. Given that the text explicitly notes heavy-heavy systems behave differently, the applicability of the same functional form should be justified explicitly or the discussion should be revised.","section":"Eq. (19) and Refs. [46,47]"},{"comment":"The concluding paragraph mentions implications for 'extra dimensions, supersymmetry, or dark matter.' This is unsupported by the analysis and should be removed or replaced with a more grounded statement about possible experimental consequences.","section":"Section 3, Discussion"},{"comment":"The notation Π(M^2, s(T), T) is used both for the Borel-transformed correlation function and, via the derivative Π', for the mass sum rule in Eq. (15). The definition of Π' as d/d(-1/M^2) is standard, but the authors should clarify which object is meant in each equation, especially since the derivative acts on a function that also depends on s(T) and M^2.","section":"Eqs. (14) and (15)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the T=0 benchmark is a useful consistency check. The central thermal claim, however, hinges entirely on the unvalidated parametrization of Eq. (19), which the authors themselves concede is not appropriate for heavy-heavy systems without a separate treatment. The lack of propagated uncertainties further weakens the quantitative claims. I believe these issues are fixable within the manuscript's scope, so major revision is appropriate rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look if you care about fully-heavy exotics in medium, but keep the caveats in hand. What is actually new: the temperature-dependent mass and decay constant of X(6200) modeled as an eta_c eta_c molecule with J^PC=0++, plus explicit thermal spectral densities in the appendix. The T=0 mass comes out at about 6.20 GeV, consistent with the LHCb peak and with the earlier Agaev et al. sum rule. The calculation follows standard TQCDSR steps, and the analytical expressions are given in full, which is more than many papers in this area do. Credit where due: this is a competent, transparent application of a known method to a newly observed state.\n\nThe soft spots are real and they are central. The thermal prediction is controlled by Eq. (19), the parametrization s(T) = s0(1-(T/Tc)^8) + 16mc^2(T/Tc)^8. The exponent 8 and the endpoint 16mc^2 are asserted, not derived. The paper itself notes that heavy-heavy thresholds deviate from the light-quark behavior cited for this form, but then uses it anyway without a separate treatment. Since the sum rule integrals are cut at s(T), the falling mass and decay constant are largely a restatement of the assumed threshold collapse. The stress-test's phrasing is fair: the 7% mass drop and 62% decay-constant drop are not independent predictions of the OPE. Second, the errors are not propagated into the temperature curves. Table 3 gives uncertainties at T=0 and T=0.14 only, and it is odd that the T=0.14 error is smaller than the T=0 error; the statistical significance of the drop is therefore not established. Third, the continuum window looks physically strained: s0 around 44-45 GeV^2 is far above the eta_c eta_c threshold near 35.6 GeV^2, and the paper does not defend the interpretation of the continuum in that region. The closing paragraph about extra dimensions, supersymmetry, and dark matter is unsupported and should go. The T=0 benchmark leans on a sum rule from an overlapping group, but that is not disqualifying given the explicit expressions here.\n\nThe qualitative direction—thermal suppression and melting—is plausible and consistent with other TQCDSR studies. The calculation is internally coherent under its stated assumptions. But the headline numbers are conditioned on an unvalidated ansatz. This paper deserves a serious referee: the question is meaningful, the method is standard, and the presentation is honest about the parametrization even if it does not justify it. A referee should ask for a sensitivity analysis of the thermal results to the s(T) functional form, proper error propagation, and removal of the speculative remarks. I would not desk-reject it, but I would send it back for revision.","headline":"A competent but assumption-driven thermal QCD sum rule study of X(6200) as an eta_c eta_c molecule: the T=0 mass matches, but the predicted thermal drop is largely written in by the ad hoc s(T) ansatz.","tokens_in":14460,"tokens_out":1751,"would_cite":false,"duration_ms":22573,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.-t","12.38.Lg"],"model":"deepseek-v4-flash","headline":"This paper predicts that X(6200), if an $\\eta_c\\eta_c$ molecule, gets lighter and more weakly coupled as temperature rises toward 155 MeV.","keywords":["X(6200)","fully-charmed exotic hadron","eta_c eta_c molecule","thermal QCD sum rules","temperature-dependent mass","gluon condensate","0++ scalar state","heavy-ion collisions"],"falsifier":"Recompute the sum rules with a continuum-threshold temperature dependence not of the assumed $s(T)=s_0(1-(T/T_c)^8)+16m_c^2(T/T_c)^8$ form, or compute the thermal spectral function directly on the lattice for the $J^{PC}=0^{++}$ $\\eta_c\\eta_c$ channel; if the pole mass moves by much less than the predicted ~7% between $T=0$ and $T=0.14$ GeV, the central claim fails.","tokens_in":13239,"feed_emoji":"🔥","tokens_out":11730,"duration_ms":125468,"temperature":0.7,"pith_summary":"The paper sets out to predict how the newly observed fully-charmed state X(6200) behaves in hot strongly interacting matter, under the hypothesis that it is a weakly bound molecule of two $\\eta_c$ mesons with $J^{PC}=0^{++}$. Using thermal QCD sum rules up to dimension four, it derives temperature-dependent sum rules for the state's mass and decay constant. The central result is that both quantities fall as the temperature rises: at $T=0$ the mass is $6201.74\\pm57.82$ MeV and the decay constant is $(2.31\\pm0.18)\\times10^{-2}$ GeV$^4$, while at $T=0.14$ GeV, close to $T_c=155$ MeV, the mass has dropped about 7% to $5782.43\\pm37.42$ MeV and the decay constant has dropped about 62% to $(0.88\\pm0.05)\\times10^{-2}$ GeV$^4$. If this is right, X(6200) is a sensitive probe of hot QCD matter, and its thermal softening should be visible in heavy-ion collisions.","feed_headline":"X(6200) mass drops 7% as temperature nears 155 MeV","feed_subtitle":"Sum-rule study predicts the charm-meson molecule softens sharply in hot QCD matter, with its decay constant falling 62%.","key_machinery":"The calculation rests on the thermal correlation function $\\Pi(p,T)=i\\int d^4x\\, e^{ip\\cdot x}\\langle\\Omega|\\mathcal{T}J(x)J^\\dagger(0)|\\Omega\\rangle$, built from the interpolating current $J(x)=(\\bar c_a i\\gamma_5 c_a)(\\bar c_b i\\gamma_5 c_b)$ for two pseudoscalar $\\eta_c$ mesons. Wick-contracting the charm fields produces a QCD-side spectral density, while the phenomenological side is a single pole plus continuum; Borel transformation and quark-hadron duality connect the two sides. The mass and decay constant are extracted as $m(T)=\\sqrt{\\Pi'/\\Pi}$ and $f^2(T)=e^{m^2/M^2}\\Pi/m^2$, where the Borel-transformed function $\\Pi$ is integrated up to a temperature-dependent continuum threshold $s(T)=s_0(1-(T/T_c)^8)+16m_c^2(T/T_c)^8$. This threshold ansatz, with its eighth-power interpolation between $s_0$ and $16m_c^2$, is the mechanism that carries the predicted thermal drop.","core_discovery":"The paper claims that the X(6200) resonance, interpreted as an $\\eta_c\\eta_c$ molecular state with $J^{PC}=0^{++}$, has a mass and a decay constant that both decrease monotonically as the temperature of the medium rises. At zero temperature the computed values reproduce earlier QCD sum-rule results for this molecule, and as the temperature approaches $T_c=155$ MeV the decrease accelerates, with the decay constant falling much faster than the mass. The thermal trend is read as a reduction of the state's binding strength in hot matter, making fully-charmed molecular states natural probes of the quark-gluon plasma and of hadronic matter under extreme conditions.","pith_inferences":["Because the threshold function $s(T)$ is chosen rather than derived, the quantitative sizes of the predicted drops are less secure than the qualitative falling trend; a physically motivated threshold could shift the numbers substantially.","A lattice-QCD calculation of the thermal spectral function in the $J^{PC}=0^{++}$ double-$\\eta_c$ channel could test the predicted pole-mass shift without relying on the threshold ansatz.","Repeating the same thermal sum-rule analysis for the alternative compact tetraquark assignment of X(6200) could show whether the melting pattern distinguishes molecular from compact fully-charmed structures."],"forward_implications":["At $T=0$ the sum rules reproduce the mass and decay constant obtained in earlier QCD sum-rule studies of the $\\eta_c\\eta_c$ molecule, which supports reading X(6200) as a molecular state.","Approaching $T_c=155$ MeV, the mass drops by about 7% and the decay constant by about 62%, with the decrease steepening near $T_c$.","In heavy-ion collisions, X(6200) should appear with a reduced mass and a much weaker coupling to the $\\eta_c\\eta_c$ channel than in vacuum.","Fully-charmed molecular states of this kind can serve as thermal probes of quark-gluon plasma and of matter in the early universe or in astrophysical environments."],"supporting_citations":[{"why":"Reports the experimental observation of the narrow ~6.2 GeV structure in double-J/psi events that motivates the whole analysis.","marker":"[7]"},{"why":"Previous QCD sum-rule study of the eta_c eta_c molecule giving the vacuum mass and width used as a comparison baseline.","marker":"[14]"},{"why":"Foundational paper introducing QCD sum rules on which the Borel sum-rule method is built.","marker":"[37]"},{"why":"Establishes the finite-temperature QCD sum-rule framework used throughout the paper.","marker":"[38]"},{"why":"Sets the adopted critical temperature Tc=155 MeV used in the temperature parametrization.","marker":"[39]"},{"why":"Provides the thermal gluon condensate decomposition into C1(T) and C2(T) used for the dimension-four contributions.","marker":"[43]"},{"why":"Justifies a separate temperature-dependent continuum threshold for heavy-heavy quark systems, the basis of the s(T) ansatz.","marker":"[46]"},{"why":"Supplies the charm quark and eta_c masses used as numerical inputs.","marker":"[48]"},{"why":"Supplies the vacuum gluon condensate value used at zero temperature.","marker":"[49]"}],"fun_headline_variants":["X(6200) mass and decay constant fall in hot QCD","Thermal sum rules predict softening of charmonium molecule","Hot QCD weakens fully-charmed molecular state","X(6200) loses binding as temperature rises","Sum-rule study: mass drops 7%, decay constant 62% at 155 MeV"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes a particular un-derived temperature dependence for the energy cutoff that separates the ground state from higher resonances, and the predicted 7% mass drop and 62% decay-constant drop are largely consequences of that choice.","fun_headline_variants_meta":{"raw":{"variants":["X(6200) mass and decay constant fall in hot QCD","Thermal sum rules predict softening of charmonium molecule","Hot QCD weakens fully-charmed molecular state","X(6200) loses binding as temperature rises","Sum-rule study: mass drops 7%, decay constant 62% at 155 MeV"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000184,"raw_usage":{"total_tokens":1262,"prompt_tokens":833,"completion_tokens":429,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":449,"completion_tokens_details":{"reasoning_tokens":339}},"tokens_in":449,"tokens_out":429,"duration_ms":4869,"temperature":1.0,"reasoning_tokens":339,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:07:50.412710+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the sum rules with a continuum-threshold temperature dependence not of the assumed $s(T)=s_0(1-(T/T_c)^8)+16m_c^2(T/T_c)^8$ form, or compute the thermal spectral function directly on the lattice for the $J^{PC}=0^{++}$ $\\eta_c\\eta_c$ channel; if the pole mass moves by much less than the predicted ~7% between $T=0$ and $T=0.14$ GeV, the central claim fails.","supporting_citations":[{"cited_title":"Aaij et al","cited_arxiv_id":null,"evidence_quote":"Reports the experimental observation of the narrow ~6.2 GeV structure in double-J/psi events that motivates the whole analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Previous QCD sum-rule study of the eta_c eta_c molecule giving the vacuum mass and width used as a comparison baseline."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the finite-temperature QCD sum-rule framework used throughout the paper."},{"cited_title":"Mallik, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the thermal gluon condensate decomposition into C1(T) and C2(T) used for the dimension-four contributions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Justifies a separate temperature-dependent continuum threshold for heavy-heavy quark systems, the basis of the s(T) ansatz."},{"cited_title":"Narison, Int","cited_arxiv_id":null,"evidence_quote":"Supplies the vacuum gluon condensate value used at zero temperature."}],"review_version":1}