{"id":"8a49376e-cf04-40c7-b820-b2a994f3b8e6","arxiv_id":"2411.19296","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Light quark fragmentation into S-wave fully-charmed tetraquarks is computed at leading order in NRQCD, giving production rates between the gluon and charm channels at the LHC and EIC.","lead":"This paper computes how often a light quark (up, down, or strange) turns into a fully charmed tetraquark, a hypothetical particle made of four charm quarks, inside high-energy collisions. It predicts that this new fragmentation channel is less common than the gluon channel but more common than the charm quark channel at the LHC and EIC.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Main concern: the central ordering claim relies on comparing Dq against gluon/charm curves whose LDME inputs and evolution setup are not stated; if those curves used different model inputs, the quoted ratios may be artifacts.","rationale":"The reader's conditional verdict is appropriate, and my analysis agrees that the LDME normalization is a serious limitation. However, the more load-bearing aspect for the paper's central claim is the uncontrolled comparison: the headline statement is a ratio-like ordering, not an absolute cross section. If the gluon and charm fragmentation inputs from [43,50] were obtained with different LDME values or a different evolution treatment, the ordering could be an artifact of inconsistent bookkeeping even if every individual calculation is correct. The perturbative short-distance part appears sound: the mq -> 0 limit reproduces the compact forms in Eq. (21), the log identities check out, and the fragmentation functions are positive in the physical region. The concern is therefore not a mathematical inconsistency in the SDC derivation but a missing control on the phenomenological comparison. I would keep the paper as conditional rather than reject or accept outright: the new calculation is a legitimate first step, but the quantitative claim needs a unified-input recomputation and a spread estimate across the available model LDMEs before the stated Dg > Dq > Dc ordering can be taken as a robust prediction.","tokens_in":13881,"tokens_out":25676,"duration_ms":231321,"concrete_test":"Regenerate the Dq, Dg, and Dc cross sections entering Table II using one common pipeline: CT14llo PDFs, mu0 = 4mc, the same alpha_s running, the same Runge-Kutta DGLAP solver, and the same Table I LDMEs. Then repeat with the other four LDME sets listed in [50]. If the ratios Dq/Dg and Dq/Dc change by more than the quoted 1-2 orders of magnitude, the ordering claim is not robust. If they persist, the central claim survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section V and Table II present the phenomenological payoff as a comparison: Dg > Dq > Dc, and this ordering is the paper's headline. The calculation of Dq itself (Eqs. (18)-(21)) has a clean mq -> 0 limit and appears internally consistent. The load-bearing step is the normalization and the comparison: the LDMEs in Table I come from a single potential model via vacuum saturation (Eq. (22)), with no uncertainty, even though the footnote acknowledges five model sets enumerated in [50]; and the Dg and Dc curves and cross sections are taken from [43,50] without stating which LDME values, initial scale, alpha_s scheme, or DGLAP solver they used. If the previously published curves used a different wavefunction set or omitted the interference term, the ratios in Table II (e.g., Dq/Dc = 37 for 0++ 1S at the LHC) could change substantially. The absolute scale of every channel is proportional to LDMEs, so unrepresentative wavefunction input shifts all absolute predictions; but the claimed ordering is only meaningful if all three channels are evaluated with the same inputs. The paper should at minimum state the inputs used for Dg and Dc and quantify the spread across the five model LDME sets.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper computes, within the NRQCD factorization framework, the fragmentation function of a light quark into S-wave fully-charmed tetraquarks at leading order in alpha_s and v. The authors derive closed-form short-distance coefficients for the 0++ and 2++ tetraquark channels, take the massless light-quark limit, and then use LDMEs extracted from a single potential model via the vacuum-saturation approximation. They evolve the fragmentation functions with the DGLAP equation and predict pT distributions and integrated cross sections for T4c production at the LHC (13 TeV) and EIC (140 GeV), comparing their new light-quark channel with the previously computed gluon and charm-quark fragmentation channels. The headline result is that the light-quark fragmentation contribution is smaller than the gluon one but larger than the charm-quark one at the LHC, while at the EIC it is comparable to the charm-quark contribution.","tokens_in":14129,"tokens_out":9040,"duration_ms":71851,"significance":"If the calculation is correct, the paper fills a missing channel in the NRQCD factorization treatment of fully-charmed tetraquark production and provides a useful step toward a complete high-pT production inventory. The perturbative matching is standard in structure, the mq-to-0 limit of the SDCs is well behaved and presented in closed form, and the evolution setup is clearly described at the level of the method. The paper honestly flags the model dependence of the LDMEs in a footnote, but the phenomenological comparison and the absolute predictions rest on inputs whose consistency and uncertainty are not currently documented. The value of the paper lies mainly in the new SDC calculation and the relative-size analysis; those elements are worth publishing once the comparison inputs and LDME uncertainties are made explicit.","major_comments":[{"comment":"The comparison curves and numbers for gluon and charm fragmentation are imported from Refs. [43,50], but the manuscript does not state which LDME set, initial scale, alpha_s scheme/value, DGLAP evolution treatment, and kinematic cuts were used for those channels. Since all channels are proportional to LDMEs, the claimed ordering D_g > D_q > D_c and the ratios in Table II are only meaningful if the three channels are evaluated with identical inputs. The footnote in Section V admits that five model LDME sets are enumerated in [50]; the paper should state the inputs used for D_g and D_c and show how the ratios in Table II change across those five sets.","section":"§V, Figs. 3–4 and Table II"},{"comment":"The LDMEs are taken from a single potential model [11] with no quoted uncertainty. Because the absolute cross sections in Table II and Figs. 3–4 scale linearly with these LDMEs, the choice of model is a load-bearing input for every numerical prediction. The authors should quote at least the range spanned by the five model sets mentioned in the footnote, and show how the absolute cross sections and the relative ordering of D_q, D_g, and D_c behave under that spread.","section":"§V, Table I and Eq. (22)"},{"comment":"It is not stated whether D_q in the phenomenological curves represents a single light quark flavor or the sum over u, d, s (and their antiquarks) with the corresponding PDFs. The integrated cross sections in Table II depend directly on this choice. The authors should specify the flavor-sum convention, list the active light flavors, and state how the subprocesses in Eq. (2) are combined in the numerical evaluation.","section":"§V, Eqs. (1)–(3) and Fig. 3 caption"}],"minor_comments":[{"comment":"The numerical value of alpha_s used in the SDCs and partonic cross sections is not given; since the SDCs scale as alpha_s^4, the scheme, order, and input value at the initial scale mu = 4m_c are needed for reproducibility.","section":"Section V"},{"comment":"The DGLAP evolution is described as a Runge-Kutta solution including all parton channels, but the number of active flavors, the starting scale, and the treatment of flavor thresholds are not specified.","section":"Section V"},{"comment":"For the 2++ state, only the LDME <O^(2)_{3,3}> is listed; a sentence explicitly stating that the other color-spin configurations are forbidden by Fermi statistics, consistent with the diquark-antidiquark spin assignment in Section III, would prevent confusion.","section":"Table I and Section III"},{"comment":"The factors 16 and 80 in the vacuum-saturation relation are stated without derivation; a brief reference to the origin of these factors, or an explicit demonstration, would help the reader verify the normalization of Table I.","section":"Eq. (22)"}],"recommendation":"major_revision","confidential_remarks":"The perturbative calculation appears sound and the paper is transparent about the model dependence of the LDMEs, but the headline ordering is presented without a documented common-input comparison across the three fragmentation channels. This is fixable within the manuscript's scope by adding the input specifications and an uncertainty scan over the five LDME sets; I would not reject on the current evidence, but the comparison needs to be made rigorous before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is Eqs. (21): the first explicit short-distance coefficients for light-quark fragmentation into S-wave fully-charmed tetraquarks at LO in alpha_s and v. The mq -> 0 limit is well behaved, the matching follows the same standard NRQCD procedure as the gluon and charm channels, and the formulas are compact enough to be usable. The paper also makes a clear quantitative claim: at high pT, light-quark fragmentation sits between gluon and charm fragmentation in size. That ordering, if correct, completes the leading-parton fragmentation picture for T4c.\n\nWhat the paper does well: the perturbative calculation is straightforwardly presented, the z-dependent expressions are explicit, and the comparison to prior channels is physically motivated. The central calculation looks internally consistent. I have no reason to doubt the SDCs.\n\nThe soft spots are real but not fatal. The absolute normalization comes from Table I, which uses a single potential model (Lu-Chen-Dong) via vacuum saturation, and no uncertainty is quoted even though the footnote acknowledges that five model LDME sets are enumerated in [50]. For a paper whose headline is a comparison between channels, that matters. More concretely, the gluon and charm curves in Figs. 3-4 and Table II are taken from [43,50] without stating which LDME values, initial scale, alpha_s scheme, or DGLAP solver were used for those curves. If those comparisons used different wavefunction inputs, the quoted ratios (e.g., Dq/Dc ~ 37 for 0++ 1S at the LHC) could shift substantially. The authors should state the inputs for all three channels and, ideally, show the spread across the five LDME sets. A second, minor issue: the DGLAP evolution is described only as 'Runge-Kutta'; a few lines about the grid and stopping criteria would help reproducibility.\n\nThe central perturbative result holds up. The phenomenological ordering is plausible but under-documented. This paper deserves a serious referee; it is a legitimate completion of a framework, not a desk reject. I would send it to review with the request that the authors clarify the comparison inputs and add a quantified LDME uncertainty.","headline":"First light-quark fragmentation functions for S-wave T4c at LO, with a clean mq->0 limit and solid perturbative matching; the phenomenological comparison needs documented inputs and LDME uncertainties.","tokens_in":14670,"tokens_out":2412,"would_cite":true,"duration_ms":21357,"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":"This paper establishes that light-quark fragmentation into the S-wave fully-charmed tetraquark $T_{4c}$ gives a larger high-$p_T$ production rate than charm-quark fragmentation at the LHC and a comparable rate at the EIC, so the…","keywords":["fragmentation function","fully-charmed tetraquark","NRQCD","short-distance coefficient","X(6900)","LHC","EIC","DGLAP evolution"],"falsifier":"An independent perturbative recomputation of the $O(\\alpha_s^4)$ short-distance coefficients for $q\\to T_{4c}$ that disagrees with the closed-form normalizations in Eqs. (21) would falsify the calculation; so would a lattice determination of the four-quark matrix elements that contradicts the vacuum-saturation/potential-model values by amounts large enough to break the predicted channel ordering.","tokens_in":13670,"feed_emoji":"⚛️","tokens_out":11391,"duration_ms":89553,"temperature":0.7,"pith_summary":"This paper establishes that a light quark can fragment into an S-wave fully-charmed tetraquark ($T_{4c}$), and that at high transverse momentum this process is more likely than charm-quark fragmentation while remaining below gluon fragmentation. Within the nonrelativistic QCD (NRQCD) factorization framework, the authors compute the leading-order (in $\\alpha_s$ and $v$) short-distance coefficients for light-quark fragmentation into the $0^{++}$ and $2^{++}$ tetraquark states, and combine them with long-distance matrix elements from a potential model. Using these fragmentation functions, they predict the $T_{4c}$ production cross sections at the LHC and the EIC. The light-quark channel fills the gap between the previously known gluon and charm channels, and its size matters for interpreting the $X(6900)$ resonance as a fully-charmed tetraquark.","feed_headline":"Light-quark channel beats charm in tetraquark production","feed_subtitle":"At LHC high pT, the new fragmentation channel sits between gluons and charm, refining X(6900) production rates.","key_machinery":"The carrying object is the NRQCD factorization formula $D_{q\\to H}(z)=\\sum_n d_n(z)\\langle O^H_n\\rangle$. For $T_{4c}$ the operators are the color-singlet four-quark operators $O^J_{3,3}$, $O^J_{6,6}$, and the interference term $O^J_{3,6}$, written in the diquark-antidiquark basis ($\\bar{\\mathbf{3}}\\otimes\\mathbf{3}$ and $\\mathbf{6}\\otimes\\bar{\\mathbf{6}}$). The new input is the set of perturbative short-distance coefficients $d_{3,3}$, $d_{6,6}$, $d_{3,6}$ for $J=0$ and $d_{3,3}$ for $J=2$, obtained by matching on-shell amplitudes for $q\\to c\\bar{c}c\\bar{c}$ onto the operator basis; in the $m_q\\to 0$ limit they reduce to the expressions in Eqs. (21). DGLAP evolution then mixes the quark and gluon fragmentation channels, which matters because the light-quark function gains a significant small-$z$ component from gluon splitting.","core_discovery":"The paper's central claim is that the fragmentation function $D_{q\\to T_{4c}}(z,\\mu)$ at leading order in $\\alpha_s$ and $v$ is fixed by the short-distance coefficients in Eqs. (21), with the nonperturbative normalization carried by four-quark NRQCD matrix elements. In the massless-light-quark limit these coefficients take compact closed forms, such as $(z-4)^2(z-1)/(z(z-2))$ times $\\pi^2\\alpha_s^4/864$ for the $(6,6)$ color channel of the $0^{++}$ state. Combining the evolved fragmentation functions with partonic cross sections and standard proton PDFs gives integrated cross sections at the 13 TeV LHC and 140 GeV EIC. The predicted ordering at the LHC is gluon $>$ light quark $>$ charm, while at the EIC the light-quark and charm-quark contributions are comparable and both sit below the gluon one.","pith_inferences":["Because the same long-distance matrix elements normalize all three fragmentation channels, the predicted ordering $D_g > D_q > D_c$ at the LHC is largely insensitive to the potential-model uncertainty; only the absolute rates, not the ratios, carry that normalization.","A lattice QCD evaluation of the four-quark matrix elements would convert the current single-model normalization into a first-principles prediction and would directly test the size of the light-quark contribution relative to the gluon one.","The $m_q\\to 0$ closed forms in Eqs. (21) have a distinctive $z$-dependence; measuring the $p_T$ shape of $T_{4c}$ at the LHC could in principle separate the light-quark contribution from the charm-quark one, since their evolved shapes differ.","The same short-distance coefficient matching technique transfers to fully-bottom tetraquarks, where the larger heavy-quark mass changes the coupling and the velocity expansion, potentially altering the ordering between quark and gluon fragmentation."],"forward_implications":["At the LHC with $p_T\\ge 20$ GeV, light-quark fragmentation yields $8.1\\times 10^3$ fb for the $1S$ $0^{++}$ state and $7.5\\times 10^3$ fb for the $1S$ $2^{++}$ state, roughly one to two orders of magnitude below gluon fragmentation but one to two orders above charm-quark fragmentation.","At the EIC (140 GeV), the light-quark and charm-quark fragmentation channels contribute comparably to $T_{4c}$ production, while gluon fragmentation remains the dominant source by at least an order of magnitude.","DGLAP evolution substantially increases the fragmentation function at small and intermediate $z$ through mixing with the gluon channel, so predictions at LHC and EIC scales must use the evolved functions rather than the initial-scale ones.","The same short-distance coefficients, with the light-quark mass replaced by the bottom mass, provide the initial-scale fragmentation function for $b\\to T_{4c}$; analogous substitutions give $c\\to T_{4b}$, so the formulas cover all quark-initiated channels.","A complete $T_{4c}$ high-$p_T$ production calculation must sum gluon, charm, and light-quark fragmentation together with their DGLAP mixing, since no single channel dominates the quark-initiated contribution."],"supporting_citations":[{"why":"gives the QCD factorization theorem for high-$p_T$ hadron production that turns fragmentation functions into cross sections in Eq. (1).","marker":"[51]"},{"why":"supplies the gluon fragmentation short-distance coefficients and the matching method for $T_{4c}$ that the light-quark calculation extends.","marker":"[43]"},{"why":"provides the charm and gluon fragmentation functions and the integrated cross sections that the light-quark results are compared against.","marker":"[50]"},{"why":"establishes the vacuum-saturation relations connecting the NRQCD LDMEs to the wavefunction at the origin, used in Eq. (22).","marker":"[48]"},{"why":"supplies the four-body wavefunction at the origin from which the LDME values in Table I are estimated.","marker":"[11]"},{"why":"defines the gauge-invariant quark fragmentation function used as the starting point in Eq. (7).","marker":"[55]"}],"fun_headline_variants":["Light quarks beat charm in tetraquark fragmentation","Tetraquark production: light quarks outpace charm at LHC","Light-quark fragmentation tops charm for X(6900) at LHC","Gluons lead, light quarks edge charm in tetraquark yields","Light quark fragmentation: a middle path for tetraquarks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the long-distance matrix elements, estimated from a single potential model via the vacuum saturation approximation, set the normalization of all three fragmentation channels; if that model is unrepresentative, the absolute cross sections in Table II shift by potentially orders of magnitude even though the perturbative calculation itself would remain valid.","fun_headline_variants_meta":{"raw":{"variants":["Light quarks beat charm in tetraquark fragmentation","Tetraquark production: light quarks outpace charm at LHC","Light-quark fragmentation tops charm for X(6900) at LHC","Gluons lead, light quarks edge charm in tetraquark yields","Light quark fragmentation: a middle path for tetraquarks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000171,"raw_usage":{"total_tokens":1234,"prompt_tokens":869,"completion_tokens":365,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":485,"completion_tokens_details":{"reasoning_tokens":273}},"tokens_in":485,"tokens_out":365,"duration_ms":3513,"temperature":1.0,"reasoning_tokens":273,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:20:06.773727+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent perturbative recomputation of the $O(\\alpha_s^4)$ short-distance coefficients for $q\\to T_{4c}$ that disagrees with the closed-form normalizations in Eqs. (21) would falsify the calculation; so would a lattice determination of the four-quark matrix elements that contradicts the vacuum-saturation/potential-model values by amounts large enough to break the predicted channel ordering.","supporting_citations":[],"review_version":1}