REVIEW 4 major objections 4 minor 4 cited by
On the quarkonium-in-jet collinear fragmentation at moderate-to-large transverse momentum
T0 review · 4 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read 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…
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (4)
- [Abstract; Sec. 1] 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.
- [Sec. 2; Fig. 1; Sec. 4] 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.
- [Sec. 4] 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.
- [Sec. 2] 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.
minor comments (4)
- [Sec. 3] 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'.
- [Sec. 3] The phrase 'de factotranslates' in the paragraph following Eq. (2) is missing a space and should read 'de facto translates'.
- [Fig. 1] 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.
- [Eq. (2)] 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.
Circularity Check
No circularity: the load-bearing equations are standard cited factorization results; NRFF1.0 is a proposed input set, not a fitted prediction.
full rationale
The paper's central equations are not derived from the quantities they are meant to predict. Eq. (1) is the standard SIFJF convolution quoted from Kang et al. (refs. [70,72]), and Eq. (2) is an explicitly labeled schematic for two-step DGLAP evolution. Neither equation is constructed so that its output reproduces NRFF1.0 by definition. The NRFF1.0 inputs are NLO NRQCD-based color-singlet initial-scale fragmentation functions; the paper performs no fit to data and makes no numerical prediction that could be forced by a fitted parameter. The assertion that leading-twist single-parton fragmentation prevails is a physics assumption, not a self-referential derivation, and the unsupported claim that NRFF1.0 'will supersede' ZCW19+/ZCFW22 is a statement of ambition rather than a circular reduction. The framework relies heavily on the author's own prior work (HF-NRevo refs. [39,40], symJETHAD refs. [45-54]), but self-citation alone is not circularity: the cited methodology is described in the paper, and no external 'uniqueness theorem' or fitted quantity is invoked to force the announced results. Because no central claim reduces by construction to its input, the circularity score is 0.
Assumptions & free parameters
free parameters (2)
- NRQCD color-singlet LDMEs
- Initial fragmentation scale mu0 per parton channel
assumptions (5)
- domain assumption Leading-twist collinear single-parton fragmentation prevails over higher-twist heavy-quark-pair fragmentation at moderate-to-large transverse momentum.
- domain assumption NRQCD factorization with a color-singlet Fock state truncation describes quarkonium fragmentation at the initial scale.
- standard math The SIFJF factorizes as a convolution of collinear fragmenting jet coefficients and standard FFs, Eq. (1), at leading power.
- domain assumption Time-like DGLAP evolution can be split in two steps, with matching scale muM approximately Q times RJ, so that energy resummation becomes jet-radius resummation.
- domain assumption The initial-scale FF inputs are obtainable from nonrelativistic NLO calculations.
Cite this review
Pith. "Pith review of On the quarkonium-in-jet collinear fragmentation at moderate-to-large transverse momentum." pith.science (2026). https://pith.science/paper/WIWRMU6A
@misc{pith2026241205661,
author = {Pith},
title = {Pith review of: On the quarkonium-in-jet collinear fragmentation at moderate-to-large transverse momentum},
year = {2026},
howpublished = {\url{https://pith.science/paper/WIWRMU6A}},
note = {Machine review of arXiv:2412.05661}
}
read the original abstract
We report progress on the Heavy-Flavor Non-Relativistic Evolution (HF-NRevo) setup, a novel methodology to address quarkonium formation within the fragmentation approximation. Our study sheds light on the moderate to large transverse-momentum sector, where 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. As for the initial energy-scale inputs, we rely on nonrelativistic next-to-leading calculations for all the parton-to-quarkonia fragmentation channels. Preliminary sets of variable-flavor number-scheme (VFNS) fragmentation functions, named NRFF1.0, are built via an evolution-threshold enhanced DGLAP scheme. Taking NRFF1.0 as a starting point, we use HF-NRevo to address the collinear fragmentation of quarkonia inside jets.
Figures
Forward citations
Cited by 4 Pith papers
-
Measurement of the fragmentation properties of jets containing $\Upsilon$(nS) mesons in proton-proton collisions at $\sqrt{s}$ = 13 TeV
First measurement of Upsilon(nS) jet-fragmentation profiles shows data have lower z and higher p_rel^T than PYTHIA 8.240/8.310 with CP1/CP5 tunes.
-
The Scaffold Effect: How Prompt Framing Drives Apparent Multimodal Gains in Clinical VLM Evaluation
Merely mentioning MRI availability in the prompt drives 70-80% of apparent multimodal F1 gains in clinical VLMs, even when no imaging is present.
-
High-energy dynamics of QCD: Theoretical and phenomenological results
Real NLO corrections to the BFKL Higgs impact factor are derived with a physical top-quark mass, and updated tetraquark fragmentation functions are applied to predict rates at 14 and 100 TeV.
-
All-charm tetraquarks at hadron colliders: A high-precision fragmentation perspective
TQ4Q2.0 supplies the first complete, uncertainty-quantified set of NRQCD-based fragmentation functions for all-heavy tetraquarks, including nonconstituent channels and public grids for jet-associated production.
Reference graph
Works this paper leans on
-
[8]
A. Arroyo Garcia et al., Phys. Lett. B 795, 569 (2019), 1904.04394
arXiv 2019
- [1]
-
[2]
N. Brambilla et al., Phys. Rev. D 75, 074014 (2007), hep-ph/0702079
arXiv 2007
- [3]
-
[4]
G. Altarelli et al., Nucl. Phys. B 534, 277 (1998), hep-ph/9806345
arXiv 1998
-
[5]
Candido et al., JHEP 11, 129 (2020), 2006.07377
A. Candido et al., JHEP 11, 129 (2020), 2006.07377
arXiv 2020
- [6]
- [7]
Show all 93 references
-
[9]
F. G. Celiberto et al., Phys. Lett. B786, 201 (2018), 1808.09511
2018 arXiv
-
[10]
A. D. Bolognino et al., Eur. Phys. J. C78, 1023 (2018), 1808.02395
2018 arXiv
-
[11]
A. D. Bolognino et al., Eur. Phys. J. C 81, 846 (2021), 2107.13415
2021 arXiv
-
[12]
F. G. Celiberto, Nuovo Cim. C42, 220 (2019), 1912.11313
2019 arXiv
- [13]
-
[14]
Z.-B. Kang, E. Li, and F. Salazar, JHEP 03, 027 (2024), 2310.12102
2024 arXiv
- [15]
- [16]
- [17]
- [18]
- [19]
-
[20]
F. G. Celiberto, Nuovo Cim. C44, 36 (2021), 2101.04630
2021 arXiv
- [21]
-
[22]
R. D. Ball et al. (NNPDF), Nature 608, 483 (2022), 2208.08372
2022 arXiv
- [23]
-
[24]
R. D. Ball et al., Phys. Rev. D 109, L091501 (2024), 2311.00743
2024 arXiv
-
[25]
W. E. Caswell et al., Phys. Lett. B 167, 437 (1986)
1986
-
[26]
G. T. Bodwin et al., Phys. Rev. D 51, 1125 (1995), hep-ph/9407339
1995 arXiv
-
[27]
B. Mele, P. Nason, Nucl. Phys. B 361, 626 (1991)
1991
- [28]
-
[29]
Cacciari et al., JHEP 10, 137 (2012), 1205.6344
M. Cacciari et al., JHEP 10, 137 (2012), 1205.6344
2012 arXiv
- [30]
- [31]
-
[32]
Czakon et al., JHEP 10, 216 (2021), 2102.08267
M. Czakon et al., JHEP 10, 216 (2021), 2102.08267
2021 arXiv
-
[33]
Czakon et al., JHEP 03, 251 (2023), 2210.06078
M. Czakon et al., JHEP 03, 251 (2023), 2210.06078
2023 arXiv
-
[34]
Generet, Ph.D
T. Generet, Ph.D. thesis, R WTH Aachen University (2023)
2023
-
[35]
Ghira et al., JHEP 11, 120 (2023), 2309.06139
A. Ghira et al., JHEP 11, 120 (2023), 2309.06139
2023 arXiv
-
[36]
Bonino et al., JHEP 06, 040 (2024), 2312.12519
L. Bonino et al., JHEP 06, 040 (2024), 2312.12519
2024 arXiv
- [37]
-
[38]
Cacciari et al., Phys
M. Cacciari et al., Phys. Rev. Lett. 73, 1586 (1994), hep-ph/9405241
1994 arXiv
-
[39]
F. G. Celiberto, Proceedings of Moriond QCD(2024), 2405.08221
2024 arXiv
-
[40]
F. G. Celiberto, PoS DIS2024, 168 (2024), 2406.10779
2024 arXiv
-
[41]
Forte et al., JHEP 05, 062 (2002), hep-ph/0204232
S. Forte et al., JHEP 05, 062 (2002), hep-ph/0204232
2002 arXiv
- [42]
-
[43]
M. G. Echevarria, JHEP 10, 144 (2019), 1907.06494
2019 arXiv
-
[44]
Boer et al., JHEP 08, 105 (2023), 2304.09473
D. Boer et al., JHEP 08, 105 (2023), 2304.09473
2023 arXiv
-
[45]
F. G. Celiberto, Ph.D. thesis, U. Calabria & INFN (2017), 1707.04315
2017 arXiv
-
[46]
F. G. Celiberto, Eur. Phys. J. C 81, 691 (2021), 2008.07378
2021 arXiv
-
[47]
F. G. Celiberto, Phys. Rev. D 105, 114008 (2022), 2204.06497
2022 arXiv
-
[48]
F. G. Celiberto, Eur. Phys. J. C 83, 332 (2023), 2208.14577. 6 REFERENCES
2023 arXiv
-
[49]
F. G. Celiberto et al., Eur. Phys. J. C 81, 293 (2021), 2008.00501
2021 arXiv
-
[50]
A. D. Bolognino et al., Phys. Rev. D 103, 094004 (2021), 2103.07396
2021 arXiv
-
[51]
F. G. Celiberto et al., Eur. Phys. J. C 81, 780 (2021), 2105.06432
2021 arXiv
-
[52]
F. G. Celiberto et al., Phys. Rev. D 104, 114007 (2021), 2109.11875
2021 arXiv
-
[53]
F. G. Celiberto et al., Phys. Rev. D 105, 114056 (2022), 2205.13429
2022 arXiv
-
[54]
F. G. Celiberto et al., Phys. Rev. D 106, 114004 (2022), 2207.05015
2022 arXiv
-
[55]
L. A. Harland-Lang et al., Eur. Phys. J. C 79, 225 (2019), 1811.08434
2019 arXiv
-
[56]
R. D. Ball et al. (NNPDF), Eur. Phys. J. C84, 517 (2024), 2401.10319
2024 arXiv
-
[57]
Kassabov et al., JHEP 03, 148 (2023), 2207.07616
Z. Kassabov et al., JHEP 03, 148 (2023), 2207.07616
2023 arXiv
- [58]
- [59]
- [60]
-
[61]
C. W. Bauer and E. Mereghetti, JHEP 04, 051 (2014), 1312.5605
2014 arXiv
-
[62]
Chien et al., JHEP 05, 125 (2016), 1512.06851
Y.-T. Chien et al., JHEP 05, 125 (2016), 1512.06851
2016 arXiv
- [63]
-
[64]
Kang et al., JHEP 11, 068 (2017), 1705.08443
Z.-B. Kang et al., JHEP 11, 068 (2017), 1705.08443
2017 arXiv
-
[65]
Metodiev et al., Phys
E. Metodiev et al., Phys. Rev. Lett. 120, 241602 (2018), 1802.00008
2018 arXiv
- [66]
-
[67]
Kasieczka et al., JHEP 09, 195 (2020), 2007.04319
G. Kasieczka et al., JHEP 09, 195 (2020), 2007.04319
2020 arXiv
- [68]
-
[69]
P. K. Dhani et al. (2024), 2410.05415
2024 arXiv
- [70]
-
[71]
Baumgart et al., JHEP 11, 003 (2014), 1406.2295
M. Baumgart et al., JHEP 11, 003 (2014), 1406.2295
2014 arXiv
-
[72]
Kang et al., JHEP 11, 155 (2016), 1606.07063
Z.-B. Kang et al., JHEP 11, 155 (2016), 1606.07063
2016 arXiv
- [73]
-
[74]
Caletti et al., JHEP 07, 076 (2021), 2104.06920
S. Caletti et al., JHEP 07, 076 (2021), 2104.06920
2021 arXiv
-
[75]
Reichelt et al., JHEP 03, 131 (2022), 2112.09545
D. Reichelt et al., JHEP 03, 131 (2022), 2112.09545
2022 arXiv
-
[76]
F. G. Celiberto et al., Eur. Phys. J. C 82, 929 (2022), 2202.12227
2022 arXiv
-
[77]
F. G. Celiberto, Universe 9, 324 (2023), 2305.14295
2023
-
[78]
F. G. Celiberto, Phys. Lett. B 835, 137554 (2022), 2206.09413
2022 arXiv
-
[79]
F. G. Celiberto, Eur. Phys. J. C 84, 384 (2024), 2401.01410
2024 arXiv
-
[80]
Chapon et al., Prog
E. Chapon et al., Prog. Part. Nucl. Phys. 122, 103906 (2022)
2022
-
[81]
Amoroso et al., Acta Phys
S. Amoroso et al., Acta Phys. Polon. B 53, A1 (2022), 2203.13923
2022 arXiv
-
[82]
R. A. Khalek et al., Nucl. Phys. A 1026, 122447 (2022), 2103.05419
2022 arXiv
-
[83]
R. A. Khalek et al. (2022), 2203.13199
2022 arXiv
-
[84]
Abir et al
R. Abir et al. (2023), 2305.14572
2023
-
[85]
Adachi et al
I. Adachi et al. (ILC International Community) (2022), 2203.07622
2022
-
[86]
Allaire et al., Comput
C. Allaire et al., Comput. Softw. Big Sci. 8, 5 (2024), 2307.08593
2024
-
[87]
Hekhorn, Proceedings of DIS(2024), 2406.06083
F. Hekhorn, Proceedings of DIS(2024), 2406.06083
2024 arXiv
-
[88]
Hammou et al., JHEP 11, 090 (2023), 2307.10370
E. Hammou et al., JHEP 11, 090 (2023), 2307.10370
2023 arXiv
-
[89]
M. N. Costantini et al. (2024), 2402.03308
2024 arXiv
-
[90]
Hammou, Proceedings of Moriond QCD(2024), 2405.09270
E. Hammou, Proceedings of Moriond QCD(2024), 2405.09270
2024 arXiv
-
[91]
F. G. Celiberto et al., Phys. Lett. B 848, 138406 (2024), 2308.00809
2024 arXiv
-
[92]
F. G. Celiberto, Symmetry 16, 550 (2024), 2403.15639
2024
-
[93]
F. G. Celiberto et al., Eur. Phys. J. C 84, 1071 (2024), 2405.14773
2024 arXiv
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.