{"id":"8666789c-40cb-4b6b-951d-b3940d01e4f2","arxiv_id":"2505.06910","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A parameter-free formalism, based on one color-octet channel, reproduces LHC hadroproduction rates of χc1(3872) and supports its tetraquark interpretation.","lead":"This paper calculates how often the particle χc1(3872) is produced at the LHC under the assumption that it is a four-quark state. The computed rates match the measured ones, which the authors read as support for the tetraquark hypothesis.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Absolute cross-section predictions are not parameter-free: |phi(0)|^2 is fitted to LHCb branching data in Sec. IV.B, so only the nonprompt fraction tests the tetraquark hypothesis independently.","rationale":"The reader's weakest assumption correctly points to the potential-model dependence and the fitted normalization of |phi(0)|^2. My independent reading identifies essentially the same load-bearing issue, with a sharper formulation: the paper's own procedure in Sec. IV.B fixes the matrix element to an experimental branching ratio, so the absolute cross sections are not the advertised no-unknown-parameters predictions. This is a concrete, textually supported concern rather than a disagreement with the consensus. It does not invalidate the nonprompt-fraction test, which is genuinely parameter-free and does provide some support for a single-channel octet mechanism. But it means the abstract's central claim overstates the strength of the evidence from absolute rates. Since the reader already assigned a CONDITIONAL verdict with moderate confidence, this concern does not change the verdict; it reinforces the need for the authors to either present the potential-model prediction for |phi(0)|^2 as the primary absolute-rate prediction, including its spread, or revise the no-unknown-parameters wording. I therefore leave the verdict unchanged.","tokens_in":10842,"tokens_out":3828,"duration_ms":42271,"concrete_test":"Recompute the absolute prompt and nonprompt cross sections shown in Figs. 2 and 3 using |phi(0)|^2 = 0.69e-3 GeV^3 and |phi(0)|^2 = 2.0e-3 GeV^3, the two endpoints of the potential-model range quoted in Sec. IV.B, with all other inputs fixed. If either endpoint moves the predicted cross sections outside the experimental data bands (or changes them by more than the quoted uncertainties), the observed agreement in Figs. 2 and 3 is attributable to the LHCb-branching fit, not to the tetraquark potential model; the paper should then state explicitly that absolute-rate predictions contain one fitted parameter and that only the nonprompt fraction is parameter-free.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract and Sec. IV.A build the headline claim around a formalism with no unknown parameters, because Eq. (9) determines the NRQCD matrix element from the c-cbar wavefunction via potential-model input. But Sec. IV.B states that |phi(0)|^2 varies from 0.69e-3 to 2.0e-3 GeV^3 depending on the tetraquark potentials, and then says: 'for obtaining phenomenological results in this work we choose to fix |phi(0)|^2 against the experimental value' of Br(b -> chi_c1(3872)+X) x Br(chi_c1(3872) -> J/psi pi+ pi-), taking |phi(0)|^2 = 1.8e-3 GeV^3. Therefore the absolute prompt and nonprompt cross sections in Figs. 2 and 3 are not predictions of the tetraquark potential model; they are normalized to a measured decay-related rate. The agreement with hadroproduction data shown in those figures is then partly built in, and the only genuinely parameter-free result is the nonprompt fraction of Sec. III, where the matrix element cancels. Section V itself concedes that the nonperturbative matrix element was 'fixed against LHCb data' and that precision is limited by lack of knowledge of the potentials. This is an internal inconsistency between the stated no-unknown-parameters claim and the actual phenomenological procedure, and it weakens the evidence that hadroproduction data independently support the tetraquark hypothesis.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that the tetraquark (c\\bar c q\\bar q) interpretation of \\chi_{c1}(3872) proposed in Refs. [29-31] leads to an NRQCD factorization formula in which inclusive production is dominated by a single color-octet ^3S_1^{[8]} c\\bar c channel (Eq. (1)). The authors further claim that the corresponding nonperturbative matrix element \\langle O(^3S_1^{[8]}) \\rangle is fixed by the c\\bar c wavefunction at the origin through Eq. (9), leaving no unknown parameters. Using this setup, they compute the nonprompt fraction (Sec. III), prompt cross sections (Fig. 2), and nonprompt cross sections (Fig. 3) at LHC energies, with both fixed-order and threshold-resummed ingredients, and report agreement with CMS, ATLAS, and LHCb data, concluding that hadroproduction data support the tetraquark hypothesis.","tokens_in":11219,"tokens_out":5626,"duration_ms":56059,"significance":"The parameter-free test in Sec. III is genuinely valuable: because the matrix element cancels in Eq. (2), the nonprompt fraction is a direct prediction of the single-channel hypothesis, and its fair agreement with CMS and ATLAS data (Fig. 1) is a nontrivial success. The single-channel dominance itself is an interesting structural consequence of the tetraquark scenario. However, the broader claim that the formalism has no unknown parameters and that the absolute cross sections independently support the tetraquark hypothesis is not realized: the matrix element used in Figs. 2 and 3 is not taken from the potential-model wavefunction but is fixed to the measured b-decay branching-fraction product in Sec. IV.B. The paper's independent phenomenological content therefore rests mainly on the nonprompt fraction; the absolute cross-section agreement is a consistency check conditional on that input. This distinction substantially changes the strength of the conclusion.","major_comments":[{"comment":"The claim that the production formalism has no unknown parameters is not matched by the numerical procedure. Eq. (9) states \\langle O(^3S_1^{[8]}) \\rangle = 3|\\phi(0)|^2, and Sec. IV.A says the wavefunction is computable from a Schr\\\"odinger equation. But Sec. IV.B reports that the potential models of Refs. [29-31] give |\\phi(0)|^2 values between 0.69\\times10^{-3} and 2.0\\times10^{-3} GeV^3, and then \"for obtaining phenomenological results in this work we choose to fix |\\phi(0)|^2 against the experimental value\" Br(b\\to\\chi_{c1}(3872)+X)\\times Br(\\chi_{c1}(3872)\\to J/\\psi\\pi^+\\pi^-) = (4.3\\pm0.5)\\times10^{-5}, leading to |\\phi(0)|^2=1.8\\times10^{-3} GeV^3. Consequently the absolute prompt and nonprompt cross sections in Figs. 2 and 3 are normalized to a hadronic B-decay measurement rather than predicted from the tetraquark potentials. The agreement of those curves with hadroproduction data is therefore partly built in and does not independently test the tetraquark hypothesis in the way the abstract and Sec. I claim. The authors should either (a) quote the full potential-model spread as the prediction and compare it to data, or (b) explicitly demote the absolute cross sections to consistency checks and state that the no-unknown-parameters formalism has not yet been numerically realized.","section":"Sec. IV.B and Eq. (9)"},{"comment":"The concluding section itself concedes: \"the nonperturbative matrix element fixed against LHCb data [24] used in this work is within the range of potential model calculations based on Refs. [29-31]\". This admission confirms that the computation is not parameter-free and that the potential-model prediction is not precise enough to be used directly. The abstract and introduction should therefore not claim that the absolute cross sections provide a parameter-free test. What is currently demonstrated as parameter-free is the nonprompt fraction of Eq. (2), plus the single-channel dominance structure; the absolute rates are consistency checks conditioned on the measured B-decay branching-fraction product.","section":"Sec. V, Limitations"}],"minor_comments":[{"comment":"The caption of Fig. 1 does not define the meaning of the curves (FO versus resummed, and whether the nonprompt part uses FONLL with nonperturbative fragmentation). Please add a sentence describing the theoretical curves and the uncertainty band.","section":"Fig. 1 caption"},{"comment":"The text says \"We take the fixed order+next-to-leading log (FONLL) result for \\sigma_b [43,44], in which we include the correction from the nonperturbative fragmentation function\". Please clarify which nonperturbative fragmentation function is used and how it is matched to the NLO short-distance coefficient C(b\\to c\\bar c(^3S_1^{[8]})+X).","section":"Sec. III, FONLL description"},{"comment":"The statement that the mixing model of Refs. [27,28] leads to Br(b\\to\\chi_{c1}(3872)+X)Br(\\chi_{c1}(3872)\\to J/\\psi\\pi^+\\pi^-) less than about 1.1\\times10^{-5} is not accompanied by the derivation or the input values. Please include a brief explanation or a reference to where this bound is obtained.","section":"Sec. IV.C, Mixing model"},{"comment":"The dimensionless quantity S defined by the vacuum expectation value of adjoint Wilson lines should have its renormalization scale and scheme specified, since the text states S=1+O(\\alpha_s^2) and the NLO correction vanishes [46].","section":"Eq. (8)"}],"recommendation":"major_revision","confidential_remarks":"The paper's direct phenomenological evidence is narrower than the abstract claims: the nonprompt fraction is the only parameter-free test, while the absolute cross-section comparison is normalized to a measured B-decay branching-fraction product. The authors should reframe the title and conclusions accordingly. If they choose to present the absolute rates as consistency checks, the manuscript can become publishable after revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: the only genuinely parameter-free test in this paper works, but it is not the test the abstract advertises. The nonprompt fraction of χc1(3872) production, computed with the NRQCD matrix element cancelling out, matches the CMS and ATLAS data decently well once threshold and DGLAP resummation are included. The absolute prompt and nonprompt cross sections, however, are not parameter-free predictions: Sec. IV.B fixes |φ(0)|^2 to the LHCb branching-fraction product before comparing with hadroproduction data. The range of potential-model values (0.69–2.0 × 10^-3 GeV^3) is quoted, and then 1.8 × 10^-3 is chosen to match the B-decay measurement. So part of the agreement in Figs. 2 and 3 is built in. The paper's own Sec. V concedes this.\n\nWhat is actually new and good: the observation that the tetraquark hypothesis leads to single-color-octet-channel dominance in NRQCD, the derivation of ⟨O(3S1[8])⟩ = 3|φ(0)|^2, and the ratio prediction where that matrix element cancels. That is a real, falsifiable test, and it is the strongest part of the paper. The authors also do a careful job with NLO and resummed calculations, and they compare with multiple experiments at different collision energies.\n\nThe soft spots, in proportion: first, the abstract and Sec. IV.A overstate the case. There are no unknown parameters only if you accept that the potential model is predictive. By their own numbers, the potentials are tuned to the spectrum, and the value of |φ(0)|^2 varies over a factor of ~3 depending on the potential. Choosing the central value from B-decay data is a fit, not a prediction. Second, the conclusion that hadroproduction data \"well support\" the tetraquark hypothesis should be softened. The nonprompt fraction supports it; the absolute rates are consistent with it after one normalization is fixed from a different process. That is still evidence, but not independent confirmation. Third, the FO results sit below data and the resummed results overshoot at high pT; the authors mention this, and it is a minor issue given the scale uncertainties.\n\nWho this is for: hadron spectroscopists and anyone working on NRQCD factorization; the nonprompt fraction result should get attention. I would send this to referees, but I would tell the authors to fix the discrepancy between the abstract and the actual procedure. What they have is a worthwhile paper with an inflated headline, not a fatal flaw.","headline":"One clean parameter-free ratio test, but the absolute cross sections are normalized to B-decay data, so the 'no unknown parameters' claim is overstated.","tokens_in":11687,"tokens_out":2988,"would_cite":true,"duration_ms":27645,"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":"A single color-octet channel, fixed by the wavefunction, reproduces χc1(3872) production at the LHC.","keywords":["χc1(3872) tetraquark","NRQCD factorization","color-octet matrix element","hadroproduction","nonprompt fraction","potential-model wavefunction","threshold resummation","exotic hadron production"],"falsifier":"Compute the color-octet $c\\bar c$ wavefunction at the origin from first-principles QCD, for instance on the lattice or in the effective-field-theory tetraquark picture, and check whether $|\\varphi(0)|^2$ falls within the range $0.69\\times 10^{-3}$ to $2.0\\times 10^{-3}\\,\\text{GeV}^3$ spanned by the potential models used here; a value outside that range would eliminate the agreement with the measured absolute cross sections.","tokens_in":10655,"feed_emoji":"⚛️","tokens_out":15994,"duration_ms":137327,"temperature":0.7,"pith_summary":"The paper aims to show that the tetraquark interpretation of $\\chi_{c1}(3872)$ — a state in which the charm–anticharm pair sits in a color-octet $S$-wave at short distances — makes a parameter-free prediction for how often the particle is produced in high-energy proton collisions. In the nonrelativistic QCD (NRQCD) factorization formalism, only a single color-octet channel contributes at leading power, and its nonperturbative matrix element is determined by the $c\\bar c$ wavefunction at the origin. The authors compute prompt and nonprompt cross sections at the LHC and find general agreement with CMS, ATLAS, and LHCb data, including the nonprompt fraction, which is predicted without any fitted parameter. If this agreement holds, hadroproduction data would support the tetraquark hypothesis over molecular, charmonium, and mixing pictures of $\\chi_{c1}(3872)$.","feed_headline":"Parameter-free formula matches X(3872) LHC data","feed_subtitle":"One wavefunction-derived number, with no fitted parameters, reproduces both prompt and b-decay production rates.","key_machinery":"The central object is the NRQCD color-octet matrix element $\\langle O^{\\chi_{c1}(3872)}({}^3S_1^{[8]}) \\rangle$, which measures the probability that the charm–anticharm pair inside $\\chi_{c1}(3872)$ is in the color-octet spin-triplet $S$-wave state. The argument works because the tetraquark hypothesis places the $c\\bar c$ pair in exactly this Fock state at short distances, so only this one matrix element appears at leading power; the color-singlet $P$-wave and octet $D$-wave contributions are suppressed by powers of the ratio of the hadronic scale to the charm mass. The calculation then uses potential-model wavefunctions to convert this matrix element into $3|\\varphi(0)|^2$, with $\\varphi$ the $c\\bar c$ wavefunction at zero separation, so the production rate carries no unknown nonperturbative parameters once the tetraquark potential is specified.","core_discovery":"Under the tetraquark hypothesis examined here, the dominant Fock state of $\\chi_{c1}(3872)$ is a color-octet spin-triplet $S$-wave $c\\bar c$ pair, so at leading power in the nonrelativistic expansion the inclusive cross section factorizes as $\\sigma_{\\chi_{c1}(3872)} = \\sigma_{c\\bar c}({}^3S_1^{[8]}) \\times \\langle O({}^3S_1^{[8]}) \\rangle$, with no competing channels. The paper establishes that this single-channel structure makes the nonprompt fraction independent of the nonperturbative matrix element, and that the matrix element itself is $\\langle O({}^3S_1^{[8]}) \\rangle = 3|\\varphi(0)|^2$, computable from the wavefunction of the tetraquark potential model. Using a value of $|\\varphi(0)|^2$ within the potential-model range and fixed to LHCb $b$-decay data, the authors obtain prompt and nonprompt cross sections consistent with CMS, ATLAS, and LHCb measurements; the resummed predictions describe the nonprompt fraction well, whereas fixed-order results lie below data. The conclusion is that hadroproduction data support the tetraquark hypothesis.","pith_inferences":["A natural extension is to apply the same single-channel, wavefunction-determined formalism to other exotic candidates whose dominant Fock state is an $S$-wave color-octet heavy-quark pair, such as certain hybrid mesons.","The claim of no unknown parameters is conditional on the effective potentials tuned to the tetraquark spectrum; a first-principles calculation of the octet wavefunction at the origin would test it directly.","High-precision measurements of the nonprompt fraction at very large transverse momentum could discriminate the production mechanism even without absolute normalization, because the resummed and fixed-order predictions separate there.","Since the absolute normalization is tied to $|\\varphi(0)|^2$, more precise measurements of $\\chi_{c1}(3872)$ production could indirectly constrain the effective tetraquark potential."],"forward_implications":["The nonprompt fraction of $\\chi_{c1}(3872)$ production is predicted with no fitted parameter, because the matrix element cancels in the ratio; the resummed prediction agrees with CMS and ATLAS data.","Prompt hadroproduction cross sections at LHC energies are reproduced by a single color-octet channel, with resummation of threshold logarithms essential at large transverse momentum.","Nonprompt rates from $b$-hadron decays agree with ATLAS and LHCb measurements when the matrix element is fixed to the LHCb branching-fraction data.","The mixing model, which requires two nonperturbative unknowns, predicts $b$-decay production branching fractions more than three times smaller than LHCb observes, so the tetraquark picture fits the nonprompt data better.","If the tetraquark picture is right, $\\chi_{c1}(3872)$ production becomes a direct probe of perturbative color-octet $c\\bar c$ production in QCD."],"supporting_citations":[{"why":"Supplies the tetraquark hypothesis that χc1(3872) is a c-cbar q-qbar state whose c-cbar pair is color-octet S-wave at short distances.","marker":"[29–31]"},{"why":"Establishes the NRQCD factorization formalism in which the cross section becomes a short-distance coefficient times a nonperturbative matrix element.","marker":"[41]"},{"why":"Provides the NLO short-distance coefficient for b → c-cbar(3S1[8]) + X used for nonprompt production.","marker":"[42]"},{"why":"Provides the resummed b-quark production cross sections used for nonprompt rates and the nonprompt fraction.","marker":"[43, 44]"},{"why":"Provides the Fortran package used to compute the NLO prompt c-cbar(3S1[8]) cross section.","marker":"[45]"},{"why":"Supplies the threshold double-logarithm resummation that makes the prompt predictions agree with the measured rates.","marker":"[46]"},{"why":"Provides the LHCb production measurements and the branching-fraction input used to fix the wavefunction value |φ(0)|^2.","marker":"[24]"},{"why":"Supplies the potential-NRQCD techniques by which the matrix element is expressed as 3|φ(0)|^2.","marker":"[51–54]"}],"fun_headline_variants":["Tetraquark hypothesis passes X(3872) production test","No-parameter tetraquark model reproduces X(3872) rates","Hadroproduction data back tetraquark nature of X(3872)","Wavefunction-only formula nails X(3872) cross sections","X(3872) production rates support tetraquark interpretation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the $c\\bar c$ wavefunction obtained from a wave equation with potentials tuned to the tetraquark spectrum gives the true NRQCD matrix element through $\\langle O({}^3S_1^{[8]}) \\rangle = 3|\\varphi(0)|^2$; if those potentials are not reliable, the parameter-free cross-section predictions lose their meaning.","fun_headline_variants_meta":{"raw":{"variants":["Tetraquark hypothesis passes X(3872) production test","No-parameter tetraquark model reproduces X(3872) rates","Hadroproduction data back tetraquark nature of X(3872)","Wavefunction-only formula nails X(3872) cross sections","X(3872) production rates support tetraquark interpretation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000204,"raw_usage":{"total_tokens":1368,"prompt_tokens":902,"completion_tokens":466,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":518,"completion_tokens_details":{"reasoning_tokens":374}},"tokens_in":518,"tokens_out":466,"duration_ms":4548,"temperature":1.0,"reasoning_tokens":374,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:29:14.777325+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the color-octet $c\\bar c$ wavefunction at the origin from first-principles QCD, for instance on the lattice or in the effective-field-theory tetraquark picture, and check whether $|\\varphi(0)|^2$ falls within the range $0.69\\times 10^{-3}$ to $2.0\\times 10^{-3}\\,\\text{GeV}^3$ spanned by the potential models used here; a value outside that range would eliminate the agreement with the measured absolute cross sections.","supporting_citations":[{"cited_title":"FDCHQHP:A Fortran Package for Heavy Quarkonium HadroProduction","cited_arxiv_id":"1405.2143","evidence_quote":"Provides the Fortran package used to compute the NLO prompt c-cbar(3S1[8]) cross section."}],"review_version":1}