{"id":"3797f133-48a8-4d54-b753-f72ebd95f0e6","arxiv_id":"2412.05661","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A preliminary NLO set of quarkonium fragmentation functions, NRFF1.0, is introduced and proposed as input for quarkonium-in-jet fragmentation studies.","lead":"This conference proceedings describes HF-NRevo, a framework that builds quarkonium fragmentation functions from nonrelativistic QCD inputs and DGLAP evolution. It presents a preliminary set of such functions, NRFF1.0, and outlines how they would be used for quarkonium-in-jet measurements at future colliders.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Viability of NRFF1.0 in the SIFJF framework rests on unquantified leading-twist and color-octet assumptions; a jet-level comparison to full NRQCD would settle it.","rationale":"The reader's weakest_assumption already identifies the missing quantitative estimate of higher-twist and color-octet corrections. I agree and sharpen it: the CS-only nature of the NRFF1.0 inputs in Sec. 2 is directly visible in the figure caption, and the SIFJF convolution in Eq. (1) has no term for hard-scattering pair production. Together these omissions mean the central claim that NRFF1.0 is a viable input for the jet-fragmentation framework is not established by the paper. The proposed test is a concrete, existing-technique check: a full NLO NRQCD J/psi+jet calculation compared against the SIFJF convolution in the same kinematics would reveal whether the missing terms are numerically important. Since the paper is explicitly preliminary and proceedings-length, this does not warrant rejection; it supports the reader's conditional verdict. No internal mathematical error is visible in the quoted equations, and the claimed two-step evolution is standard, so the main issue is missing quantitative support rather than a demonstrated inconsistency.","tokens_in":6320,"tokens_out":8828,"duration_ms":89169,"concrete_test":"Compute the ratio R = sigma_full_NRQCD(J/psi+jet at sqrt(s)=13 TeV, anti-k_T, R_J=0.4, 30 < pT,jet < 120 GeV) / sigma_SIFJF[Eq. (1) with NRFF1.0_cs_nlo], using the same PDFs, scales, and CS LDME in both, and including CO LDMEs from a global fit in the full NRQCD calculation. If R differs from unity by more than the combined scale uncertainty across the z_J/psi bins, the leading-twist/CS-only framework misses numerically significant contributions.","verdict_should_be":"UNCHANGED","load_bearing_attack":"To support the claim that NRFF1.0 is a viable input for Eq. (1), the paper needs (i) leading-twist single-parton fragmentation to dominate over multi-parton/higher-twist contributions in the jet kinematics, and (ii) the initial-scale FFs to contain all channels that contribute at the quoted perturbative order. Condition (i) is asserted in the abstract and Sec. 1 ('leading-twist ... prevails') but never quantified. Because Eq. (1) is a convolution with standard single-parton FFs only, a hard-scattering Q-bar-Q pair produced outside the fragmentation chain is absent by construction; no estimate of that term is given. Condition (ii) is not met: Sec. 1 and the figure caption identify the inputs as 'NRFF1.0_cs_nlo' (color-singlet only), and Sec. 4 repeats 'CS initial-scale inputs'. Color-octet fragmentation channels, especially g -> J/psi via ^3S_1^[8], can be important in the moderate-pT regime; omitting them leaves the FF set incomplete even in the leading-twist sector. No comparison to data, to ZCW19+/ZCFW22, or to a full NRQCD computation is provided, so the supersession claim in Sec. 4 is unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This short proceedings article reports progress on the HF-NRevo framework and on NRFF1.0, a preliminary set of NLO NRQCD-based collinear fragmentation functions (FFs) for quarkonia, evolved in a variable-flavor-number scheme with DGLAP evolution. Two illustrative channels, b→Υ and g→Υ, are shown at three factorization scales. The manuscript then sketches how these FFs would enter the semi-inclusive fragmenting jet function (SIFJF) formula, Eq. (1), through the two-step evolution depicted in Eq. (2), and states that NRFF1.0 will supersede the ZCW19+ and ZCFW22 sets. The paper is explicitly a progress report, not a full phenomenological study.","tokens_in":6602,"tokens_out":5652,"duration_ms":55247,"significance":"If the NRFF1.0 set is delivered with the claimed NLO NRQCD inputs and validated against data or against existing FF sets, it could provide a useful new input for quarkonium-in-jet studies at moderate-to-large transverse momentum. The manuscript's strength is that it clearly connects an existing factorization formula, Eq. (1), to a concrete new FF-evolution pipeline, and it presents the qualitative behavior of two channels in a transparent way. The significance is conditional, however, because the paper contains no comparison to data, to earlier FF sets, or to a full NRQCD computation, and because two load-bearing assumptions—leading-twist dominance and the completeness of the partonic channel set—are stated but not quantified.","major_comments":[{"comment":"The manuscript repeatedly asserts that at moderate-to-large transverse momentum 'the leading-twist collinear fragmentation of a single parton prevails over the higher-twist fragmentation from a constituent heavy-quark pair produced in the hard scattering,' but no quantitative estimate or kinematic criterion is given. This premise is load-bearing because Eq. (1) convolves only standard single-parton FFs; any contribution from a hard-scattered QQbar pair not produced through a single-parton fragmentation chain is absent by construction. The manuscript should specify the transverse-momentum range in which this approximation is intended to hold and provide a concrete estimate of the neglected contribution, for example by comparing the leading-twist FF-based result with a full NRQCD calculation for a simple process such as J/psi+jet at LHC energies.","section":"Abstract; Sec. 1"},{"comment":"There is an inconsistency between the abstract, which claims NLO NRQCD inputs for 'all the parton-to-quarkonia fragmentation channels,' and the body of the paper, which states that the inputs are 'NRFF1.0_cs_nlo' and 'CS initial-scale inputs.' Color-octet channels, in particular the gluon fragmentation channel via the ^3S_1^[8] intermediate state, can be numerically important for vector quarkonia in the moderate-pT regime, and Eq. (1) sums over all partonic channels j. The manuscript gives no LDME values and no discussion of why color-octet contributions can be neglected. Either the color-octet channels should be included, or the claims should be explicitly restricted to color-singlet inputs with a quantitative justification for that restriction.","section":"Sec. 2; Fig. 1; Sec. 4"},{"comment":"The statement that NRFF1.0 'will supersede the ZCW19+ and ZCFW22 functions' is not supported by the evidence presented in this manuscript. Figure 1 shows only the new FFs at three scales, with no overlay of ZCW19+ or ZCFW22, no comparison to experimental data, and no comparison to a full NRQCD computation. The paper itself labels the FFs as 'preliminary,' so the supersession claim should be softened to an outlook item or, preferably, accompanied by a validation plot in the same figure. As written, the claim goes beyond what the data shown can establish.","section":"Sec. 4"},{"comment":"The manuscript states that MHOUs are quantified through a Monte-Carlo replica-like approach and through simultaneous variation of mu_F and mu_R by factors of 1/2 to 2, but Figure 1 contains no uncertainty bands, no replica spread, and no description of how the replica procedure is implemented or calibrated. Since the treatment of uncertainties is one of the three 'pillars' of the HF-NRevo methodology, the reader needs at least a sample uncertainty band or a precise reference to the companion papers where the replica construction is defined. Without this, the claim of quantified MHOUs is not checkable from the present manuscript.","section":"Sec. 2"}],"minor_comments":[{"comment":"The sentence 'Jet substructure measurements have recently gained prominence as a powerful sounds for core nature of the strong force' contains a typo: 'sounds' should be 'tools'.","section":"Sec. 3"},{"comment":"The phrase 'de factotranslates' in the paragraph following Eq. (2) is missing a space and should read 'de facto translates'.","section":"Sec. 3"},{"comment":"The figure caption would be clearer if it stated explicitly that the plotted quantity is z times the fragmentation function, and if the ratio panels were defined, for example as the ratio to a chosen reference scale.","section":"Fig. 1"},{"comment":"Equation (2) is only a schematic arrow diagram; a formal expression for the matching at mu_M, even in the simplest case, would make the two-step DGLAP statement in the surrounding text more checkable.","section":"Eq. (2)"}],"recommendation":"major_revision","confidential_remarks":"This is a conference-proceedings-style contribution that relies heavily on the author's own prior framework papers for HF-NRevo and symJETHAD. That self-citation pattern is understandable in a proceedings format, but the editor may wish to consider whether the journal expects a full-length research article; the current manuscript is essentially an extended abstract with two illustrative curves. The validation gaps identified in the major comments (missing comparison sets, unquantified leading-twist assumption, missing color-octet inputs) are substantive and should be addressed before publication. If the journal's proceedings format does not allow for a full validation study, at minimum the claims should be brought in line with the preliminary status of the results."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version. This is a conference proceedings, not a full paper, and it reads like one. The equations it quotes — the SIFJF factorization of Kang et al. and the two-step DGLAP evolution — are standard and quoted correctly. What is new is the NRFF1.0 set itself: two pilot curves for b→Υ and g→Υ, evolved from NLO NRQCD color-singlet inputs through the HF-NRevo machinery. That machinery is described in the author's earlier work; here we get a sketch plus two curves.\n\nThe good news: the paper is honest that these are preliminary. The HF-NRevo idea — interpret the low-scale input in a fixed-flavor scheme, match to a variable-flavor scheme with threshold-enhanced DGLAP, then evolve — is coherent. The two-step evolution for the SIFJF (first evolve the FF, then evolve the SIFJF with the jet radius as the matching scale) is a legitimate way to organize the resummation. If the method works out, it could give a more consistent NLO treatment of quarkonium-in-jet observables at HL-LHC and EIC. So the direction is sensible.\n\nThe soft spots are real. The central numerical object, NRFF1.0, is shown as two sets of curves without uncertainty bands, without comparison to the older ZCW19+/ZCFW22 sets, without data, and without the input LDMEs or the initial-scale values. The paper claims MHOUs were quantified, but no bands appear in Fig. 1. More substantively, the inputs are color-singlet only. The abstract and Sec. 1 assert that at moderate-to-large pT the leading-twist single-parton fragmentation dominates over higher-twist pair production, but no estimate of the higher-twist or color-octet contributions is given. Since the g→Υ channel is shown, omitting the color-octet g→Υ fragmentation is a noticeable gap — that channel is known to matter in this regime. And the Sec. 4 claim that NRFF1.0 \"will supersede\" the existing sets is an aspiration, not a demonstrated result.\n\nCitation-wise, the paper leans heavily on the author's own prior work, which is normal for proceedings and not itself a flaw. The reference list is broad and appropriate.\n\nWho gets value: people tracking quarkonium fragmentation or SIFJFs will want to know this program exists, and the two curves are a mild preview. As a proceedings write-up, it is acceptable. If it were aimed at a full journal paper, a referee should ask for validation: comparison to data, color-octet inputs, uncertainty bands, and a quantitative handle on the higher-twist assumption.\n\nRecommendation: I would send it to a referee if submitted as a regular paper, with the expectation of substantial revision. For a proceedings, it is fine as is.","headline":"A clean proceedings-style progress report that lays out a plausible HF-NRevo plan for quarkonium fragmentation and shows two unvalidated CS FF curves; the physics direction is sane but the quantitative support is not yet there.","tokens_in":7164,"tokens_out":3990,"would_cite":false,"duration_ms":34409,"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 claims that the NRFF1.0 fragmentation functions, built from NLO NRQCD inputs and evolved with the HF-NRevo threshold-aware DGLAP scheme, are a viable foundation for quarkonium-in-jet fragmentation and will supersede the…","keywords":["quarkonium fragmentation","NRQCD","DGLAP evolution","fragmenting jet functions","jet substructure","collinear factorization","heavy-flavor evolution","transverse momentum"],"falsifier":"Compute, in the same NRQCD framework, the ratio of the higher-twist heavy-quark-pair fragmentation contribution to the leading-twist single-parton one for quarkonium-in-jet kinematics with jet transverse momentum between about 30 and 120 GeV. If this ratio is not much smaller than one in that range, the central premise of the paper fails; alternatively, a precise LHC measurement of the $z_Q$ distribution of $J/\\psi$ or $\\Upsilon$ inside anti-$k_T$ jets that disagrees with NRFF1.0jet predictions by the size of the pair contribution would settle the question.","tokens_in":6115,"feed_emoji":"⚛️","tokens_out":6399,"duration_ms":54687,"temperature":0.7,"pith_summary":"This paper reports progress on HF-NRevo, a framework for quarkonium production in the fragmentation approximation. It builds preliminary next-to-leading-order fragmentation functions, named NRFF1.0, from non-relativistic QCD inputs for all parton-to-quarkonia channels, then uses them to set up collinear quarkonium-in-jet fragmentation via the semi-inclusive fragmenting jet function. The aim is to replace the currently used ZCW19+ and ZCFW22 sets with a unified, threshold-aware DGLAP-evolved set that is valid at moderate-to-large transverse momentum. If the framework is right, quarkonium-in-jet predictions at the HL-LHC and EIC will rest on a single consistent evolution scheme, and energy resummation will translate into a jet-radius resummation.","feed_headline":"New quarkonium fragmentation functions cover high-momentum jets","feed_subtitle":"NRFF1.0 merges NRQCD inputs with DGLAP evolution to sharpen quarkonium-in-jet predictions at the LHC and EIC.","key_machinery":"The central machinery is HF-NRevo, a three-pillar framework (interpretation, evolution, uncertainties) for quarkonium fragmentation functions. It treats short-distance formation as two-parton fragmentation in a fixed-flavor-number scheme, matches to a variable-flavor-number scheme, performs a symbolic expanded and decoupled DGLAP evolution before an all-order numerical evolution, and quantifies missing higher-order uncertainties with Monte-Carlo-like replicas. The second load-bearing object is the semi-inclusive fragmenting jet function (SIFJF) of Eq. (1), whose convolution of standard collinear FFs with NLO fragmenting-jet coefficients, followed by the two-step DGLAP evolution of Eq. (2), maps single-parton fragmentation into a quarkonium-inside-jet observable and converts energy resummation into jet-radius resummation.","core_discovery":"At the paper's center is a proposal to put quarkonium-inside-jet predictions on a single, self-consistent footing. The NRFF1.0 fragmentation-function sets combine NLO NRQCD color-singlet initial-scale inputs for all parton channels with a variable-flavor-number-scheme DGLAP evolution whose thresholds are handled first symbolically and then by all-order numerical evolution. Plugging NRFF1.0 into the semi-inclusive fragmenting jet function, Eq. (1), and running the two-step evolution of Eq. (2), from the initial scale $\\mu_0$ to a matching scale $\\mu_M \\approx Q R_J$ and then from $\\mu_M$ onward, turns energy resummation into a jet-radius resummation. The paper claims this construction yields the forthcoming NRFF1.0jet SIFJFs, which will supersede the ZCW19+ and ZCFW22 sets and give improved quarkonium-in-jet predictions at moderate-to-large transverse momentum.","pith_inferences":["One consequence the paper leaves implicit is that the same NRFF1.0 inputs could be reused for quarkonium production in processes other than jets, such as semi-inclusive deep inelastic scattering at the EIC, by swapping the fragmenting-jet coefficients for the appropriate collinear or transverse-momentum-dependent matching factors.","If the higher-twist heavy-quark-pair contribution is indeed negligible only above some transverse-momentum threshold, NRFF1.0jet predictions will carry a $p_T$-dependent systematic floor; a dedicated calculation of that contribution would turn the paper's regime assumption into a tested boundary.","The Monte-Carlo replica approach for missing higher-order uncertainties points toward a data-driven extraction of quarkonium fragmentation functions from future LHC and EIC measurements, with NRFF1.0 as the prior.","A testable extension is to measure the ratio of $\\Upsilon$ to $J/\\psi$ yields inside jets as a function of $z_Q$ and jet $p_T$; since NRFF1.0 treats the $b$ and gluon channels separately, this ratio would discriminate the framework's flavor thresholds from the ZCW19+ and ZCFW22 baselines."],"forward_implications":["NRFF1.0 will supersede ZCW19+ and ZCFW22 as the default collinear fragmentation input for vector quarkonia and charmed B mesons.","The two-step evolution yields NRFF1.0jet SIFJFs, enabling NLO quarkonium-in-jet predictions for anti-$k_T$ and cone jets with jet-radius resummation.","The same framework extends to quarkonium-in-jet angularities and other resummation-sensitive jet-substructure observables.","The threshold-aware VFNS DGLAP scheme removes ambiguities in quarkonium FF evolution by treating thresholds of all parton species symbolically before all-order evolution.","Missing higher-order uncertainties are quantified by simultaneous renormalization- and factorization-scale scans, aligned with modern PDF uncertainty methodology."],"supporting_citations":[{"why":"Defines NRQCD factorization, giving the Fock-state expansion and long-distance matrix elements that underpin the initial-scale fragmentation inputs.","marker":"[25, 26]"},{"why":"Supplies the NLO parton-to-quarkonium fragmentation calculations used as the short-distance initial-scale inputs.","marker":"[27, 28]"},{"why":"Establishes that single-parton fragmentation dominates at moderate-to-high transverse masses, justifying the fragmentation approximation.","marker":"[38]"},{"why":"Introduces the HF-NRevo framework and its interpretation, evolution, and uncertainty pillars.","marker":"[39, 40]"},{"why":"Provides the Monte-Carlo replica method adapted to quantify missing higher-order uncertainties in the fragmentation functions.","marker":"[41]"},{"why":"Gives the fixed-flavor-number-scheme to variable-flavor-number-scheme matching procedure used to cross heavy-flavor thresholds.","marker":"[42]"},{"why":"Contains the semi-inclusive fragmenting jet function master formula of Eq. (1) for a hadron inside a jet.","marker":"[70]"},{"why":"Provides the NLO collinear fragmenting jet coefficients for anti-$k_T$ and cone jet algorithms used in the SIFJF convolution.","marker":"[72]"},{"why":"Defines the ZCW19+ vector-quarkonium fragmentation functions that NRFF1.0 is designed to supersede.","marker":"[76, 77]"}],"fun_headline_variants":["NRFF1.0: unified quarkonium fragmentation for jets","NRFF1.0 set sharpens quarkonium-in-jet predictions","Quarkonium jets get sharper fragmentation functions","New NRQCD-DGLAP fragmentation functions for quarkonium jets"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that at moderate-to-large transverse momentum the quarkonium yield is dominated by leading-twist fragmentation of a single parton, with the contribution from a heavy-quark pair produced in the hard scattering being negligible; the paper gives no numerical estimate of how small that higher-twist term actually is.","fun_headline_variants_meta":{"raw":{"variants":["NRFF1.0: unified quarkonium fragmentation for jets","NRFF1.0 set sharpens quarkonium-in-jet predictions","Quarkonium jets get sharper fragmentation functions","New NRQCD-DGLAP fragmentation functions for quarkonium jets"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000507,"raw_usage":{"total_tokens":2455,"prompt_tokens":912,"completion_tokens":1543,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":528,"completion_tokens_details":{"reasoning_tokens":1482}},"tokens_in":528,"tokens_out":1543,"duration_ms":10822,"temperature":1.0,"reasoning_tokens":1482,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T20:29:42.694449+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute, in the same NRQCD framework, the ratio of the higher-twist heavy-quark-pair fragmentation contribution to the leading-twist single-parton one for quarkonium-in-jet kinematics with jet transverse momentum between about 30 and 120 GeV. If this ratio is not much smaller than one in that range, the central premise of the paper fails; alternatively, a precise LHC measurement of the $z_Q$ distribution of $J/\\psi$ or $\\Upsilon$ inside anti-$k_T$ jets that disagrees with NRFF1.0jet predictions by the size of the pair contribution would settle the question.","supporting_citations":[],"review_version":1}