{"id":"0430b2a4-5e6b-4f5b-a775-70bdab4fe067","arxiv_id":"2412.12282","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The standard operator definition of fragmentation functions is shown to follow from collinear factorization for small-mass n-hadron final states, contrary to a recently proposed modified definition.","lead":"This paper re-derives QCD factorization for e+e- annihilation into a small cluster of hadrons and concludes that the standard single-hadron fragmentation function definition, with its usual prefactor, is the correct one for n-hadron final states. It thereby defends the bulk of past dihadron analyses against a recent proposal that would modify the definition with nonuniversal prefactors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The standard-prefactor claim is shown at zeroth order in a non-gauge theory; the assertion that full QCD with Wilson lines preserves the same hard parts and DGLAP kernels is not demonstrated.","rationale":"The paper is a careful retracing of the SIA factorization argument and gives a real zeroth-order derivation plus a scalar Yukawa cross-check of the renormalization argument, so there is no basis for rejection. However, the abstract's claim about \"the same hard parts and evolution kernels\" in QCD goes beyond what is explicitly shown: the detailed derivation is restricted to a non-gauge, single-flavor theory at zeroth order in the hard part, and the gauge-theory extension with Wilson lines is asserted rather than computed. The reader's weakest assumption identified exactly this gap, and I agree with it. The missing one-loop QCD check would settle whether the standard prefactor remains truly universal for n-hadron final states or whether hidden n-dependent (or relative-momentum-dependent) counterterms appear. The paper's own admission that it cannot retrace the derivation of Ref. [29], and its caveat that reduced Mh-integrated dihadron FFs can receive extra evolution terms, reinforces that the unintegrated operator statement, while plausible and well motivated, is not yet demonstrated at the order needed for the full claim. The appropriate verdict therefore remains CONDITIONAL: accept the derivation for what it proves and require the explicit NLO QCD check before treating the hard-part and evolution-kernel universality as fully established.","tokens_in":61539,"tokens_out":5779,"duration_ms":62570,"concrete_test":"Compute the one-loop (O(α_s)) ultraviolet renormalization of the collinear dihadron fragmentation function defined by Eq. (79) with the QCD Wilson line, at fixed z and fixed relative momentum (e.g., fixed Mh and ζ), using MS. Extract the operator counterterm Z_jj' and verify that it is independent of Mh, ζ, and n and equals the standard single-hadron DGLAP kernel; if a counterterm depending on the dihadron internal variables is required, the central claim fails. This directly tests the \"straightforward\" extension asserted in Sec. II and the \"trend continues\" statement in Sec. IV.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that the standard single-hadron operator definition, with prefactor 1/(4ξ), remains the universal collinear fragmentation function for small-mass n-hadron final states, with hard parts and evolution kernels unchanged. What is actually derived (Secs. IV–V) is the factorization formula in a non-gauge, single-flavor, zeroth-order-in-the-hard-part theory. The extension to full QCD is asserted in Sec. II as \"straightforward, based on existing derivations\" and in Sec. IV as \"the trend continues to higher orders,\" not demonstrated. This matters most for the evolution-kernel part of the claim: the only explicit renormalization check (Appendix C) is a scalar Yukawa model, and the h-independence of Z in Eq. (134) is an argument, not a QCD one-loop calculation with Wilson lines and final-state cuts. If the MS counterterm for Eq. (79) with a Wilson line acquired any dependence on the n-hadron relative variables (Mh, ζ) or on n, the evolution kernel would no longer be universal and hard parts would need compensating factors. The paper itself scopes \"same evolution kernels\" to the unintegrated operator definition: Sec. IX notes that reduced FFs integrated over Mh can acquire extra evolution terms [45], so past phenomenology using reduced dihadron FFs is not automatically covered by the unintegrated statement. The concern is a missing check, not a demonstrated error.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits collinear factorization for e+e- annihilation into a small-mass n-hadron system plus unobserved particles. Working in a deliberately simplified non-gauge, single-flavor theory with a zeroth-order hard part, the authors retrace the derivation of the factorized hadronic tensor and show that the resulting fragmentation function is the standard single-hadron operator definition, with the usual 1/(4\\xi) prefactor, extended to an n-hadron final state. They verify this for two different choices of dihadron kinematical variables, compare their result in detail with the alternative definition of Ref. [29] (Appendix B), and argue in Appendix C, using a scalar Yukawa model, that introducing n-hadron-dependent prefactors would modify the renormalization factor and hence the evolution kernel. The paper concludes that the standard operator definition remains universal for multihadron states and that previous dihadron phenomenology based on that definition is not compromised.","tokens_in":61798,"tokens_out":5813,"duration_ms":56133,"significance":"If the full-QCD generalization holds, the paper settles an important formal controversy: the standard Bacchetta-Radici-style dihadron fragmentation function, rather than the modified definition advocated by the JAM collaboration in Ref. [29], is the object that follows from collinear factorization with unchanged hard parts and evolution kernels. The paper has genuine strengths: the leading-power derivation in Sec. IV is careful about which momentum components can be approximated; the identification of the factorized coefficient with the operator definition in Sec. V is explicit; Sec. VII provides a thorough treatment of Jacobian factors for two common phase-space variable choices; and Appendix B gives a concrete translation between the competing operator definitions and shows how the alternative prefactors arise from imposing a sum rule rather than from factorization. The cross-section manipulations in the simplified theory are internally coherent.","major_comments":[{"comment":"The claim that the standard n-hadron operator definition applies with the same hard parts and evolution kernels is derived only in a non-gauge, single-flavor theory at zeroth order in the hard part. The extension to full QCD is justified by the statements in Sec. II that it is \"straightforward, based on existing derivations\" and in Sec. IV that \"the trend continues to higher orders,\" but no explicit QCD all-order or one-loop calculation with Wilson lines and final-state cuts is provided. Since this extension is exactly the load-bearing part of the paper's principal conclusion, the manuscript should either supply that check or explicitly narrow the theorem to the simplified theory and present the QCD statement as an expectation rather than a verified result.","section":"Secs. II and IV, Eq. (70)"},{"comment":"The h-independence of the renormalization factor Zj'j, which is the basis for the \"same evolution kernels\" claim, is demonstrated only in a scalar Yukawa toy model with an overall constant prefactor ratio. No QCD one-loop computation of the counterterm for the operator in Eq. (79) with the Wilson line is shown. The argument that a nontrivial prefactor N(\\xi) would change Z is convincing in the toy model, but the converse statement for real QCD is not established; if the MS counterterm for Eq. (79) acquired dependence on Mh, \\zeta, or n, the evolution kernel would no longer be universal and the hard parts would require compensating factors. This is a missing check rather than a demonstrated error, but it is central to the paper's claim of unchanged DGLAP kernels.","section":"Sec. VIII, Eq. (134), and Appendix C"},{"comment":"The abstract claims that past dihadron phenomenology is reaffirmed with \"the same hard parts and evolution kernels,\" but Sec. IX explicitly acknowledges that reduced fragmentation functions obtained by integrating over Mh can acquire extra evolution terms (citing Ref. [45]). Since most phenomenological applications, including the cited works, use reduced dihadron fragmentation functions rather than the unintegrated operator definition, the paper's conclusion is broader than what the derivation establishes. The manuscript should clearly separate the universal statement for the unintegrated operator definition from the reduced quantities used in fits, for which the evolution equations are not shown to be unchanged.","section":"Sec. IX and Sec. VII A, Eqs. (108)-(112)"}],"minor_comments":[{"comment":"There is a typo in the opening paragraph: \"semi-inclusive annhilation\" should read \"semi-inclusive annihilation.\"","section":"Sec. II"},{"comment":"The heading \"Variable choice B\" is typeset as \"V aryable choice B\" in the manuscript; this should be corrected.","section":"Sec. VII"},{"comment":"The distinction between \"hatted\" approximations (\\hat{k}, \\hat{k}_2), used for hard partons, and \"tilde\" approximations (\\tilde{k}, \\tilde{k}_2), used inside collinear subgraphs, is central to the derivation. A short summary table of frames and approximations would make this section noticeably easier to follow.","section":"Sec. IV, Eqs. (46)-(52)"},{"comment":"The prefactor 1/(64\\pi^3 \\xi_1 \\xi_2) is clear after the equality \\xi_1\\xi_2=\\xi^2 x_1 x_2 is written, but for readability it would help to define x_1 and x_2 immediately before Eq. (132) rather than relying on Eq. (97).","section":"Sec. VIII, Eq. (132)"},{"comment":"The notation Z^{[O(g^2)]}_{q/s}(\\xi) and the use of the dummy index j' in Eq. (C4) is slightly confusing; the authors may wish to clarify that the sum over j' collapses to the quark-in-scalar channel at this order and that no sum over hadron types is implied.","section":"Appendix C, Eqs. (C4)-(C6)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a serious and mostly well-executed contribution to an active controversy, and the Appendix B comparison is valuable. My reservation is not about the simplified derivation, which is sound as far as it goes, but about the gap between the demonstrated result (zeroth order, non-gauge, unintegrated definition) and the advertised conclusion (full QCD, all orders, unchanged evolution kernels, reliability of reduced-FF phenomenology). I would be comfortable with acceptance after the authors either provide an explicit QCD one-loop check with Wilson lines or carefully rescope the claims and make the reduced-FF caveat prominent in the abstract."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, careful paper that makes the case that the standard single-hadron fragmentation function definition, 1/(4ξ) prefactor and all, survives extension to small-mass n-hadron final states, and that the recently proposed nonuniversal-prefactor definitions are not compatible with collinear factorization. I largely agree with the reader's conditional verdict, and I think the paper earns a serious referee.\n\nWhat's new: the explicit leading-power derivation for the multidifferential n-hadron case, and the concrete argument showing where the alternative definition goes wrong—conflating parton momentum fractions with external kinematical variables, and imposing an invalid multiplicity sum rule. The step-by-step factorization in Secs. IV–V is internally coherent. The matching of the operator definition to the factorized cross section, rather than assuming it, is done properly. The Yukawa model check in Sec. VIII and Appendix C is a nice touch: it shows that if you put an n-dependent ξ-prefactor into the operator definition, the MS renormalization factor changes and the evolution kernel would not stay universal. That is real evidence, even if it is a toy model.\n\nSoft spots: the derivation is zeroth order in a non-gauge, single-flavor theory. The claim that the extension to full QCD with Wilson lines preserves the same hard parts and DGLAP kernels is asserted, not demonstrated. The only explicit renormalization computation is in scalar Yukawa, not QCD with final-state cuts and Wilson lines. But this is a missing check, not a demonstrated error, and the paper is transparent about the simplifications. The paper also honestly says it cannot retrace the alternative derivation and that practical phenomenological impact is unclear. Those admissions are to its credit. The h-independence of the evolution kernel is argued by standard reasoning; a referee might ask for a QCD one-loop computation to close the gap, but I would not treat the absence as fatal.\n\nWho should read it: anyone working on dihadron fragmentation for transversity, and the JAM collaboration specifically. It deserves a proper peer review; the topic is contested and the response is substantive. I'd accept it for review and let the referees decide how much additional evidence, if any, is needed.","headline":"A careful, honest derivation that the standard 1/(4ξ) fragmentation function definition holds for small-mass n-hadron final states at leading power, with the all-orders QCD extension asserted rather than proven; deserves serious peer review.","tokens_in":62328,"tokens_out":1834,"would_cite":true,"duration_ms":20075,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.-t"],"model":"deepseek-v4-flash","headline":"Standard fragmentation functions hold for multi-hadron final states.","keywords":["QCD factorization","fragmentation functions","dihadron fragmentation","n-hadron final states","semi-inclusive e+e- annihilation","collinear factorization","operator definition","evolution kernels"],"falsifier":"A concrete check would be to compute the two-loop virtual correction to the n-hadron fragmentation function in a Yukawa model and see whether the same ultraviolet counterterm that renormalizes the single-hadron function also cancels all collinear poles for n=2; if a hadron-dependent prefactor is required, the standard definition is not universal.","tokens_in":1719,"feed_emoji":"⚛️","tokens_out":6122,"duration_ms":114752,"temperature":0.7,"pith_summary":"The paper retraces the collinear-factorization derivation for unpolarized semi-inclusive e+e- annihilation into a small-mass cluster of n hadrons and concludes that the standard operator definition of a fragmentation function, with its usual 1/(4\\xi) prefactor, applies unchanged when the single hadron is replaced by an n-hadron state. If correct, this preserves the universality of hard parts and evolution kernels across single- and multi-hadron fragmentation, and it removes a formal objection raised against most past dihadron phenomenology. The paper also argues that alternative definitions with extra momentum-fraction-dependent prefactors are not consistent with factorization, because they would force the hadron-dependent factors into the hard part and change the evolution kernels. The weight of the argument falls on the factorization chain being independent of the observed final state, with the extension to full QCD and higher orders asserted as straightforward.","feed_headline":"Standard fragmentation functions hold for multi-hadron final states","feed_subtitle":"Retracing e+e- factorization shows no new prefactors or evolution kernels are needed for small-mass hadron clusters.","key_machinery":"The load-bearing object is the light-cone fragmentation function operator, whose collinear form is the standard single-hadron definition with the observed state changed to an n-hadron system, namely d(\\xi,\\{ph\\}) ~ (\\xi/4) \\int [dk^- $d^{2}$kT/(2\\pi)^4] Tr[(\\gamma^+/2) J1(\\tilde{k},\\{ph\\})]. The derivation relies on separating the partonic momentum into 'hatted' collinear momenta \\hat{k},\\hat{k}_2 for the hard part and 'tilde' momenta \\tilde{k},\\tilde{k}_2 for the collinear subgraphs, and on keeping the distinction between the partonic momentum fractions \\xi_i and the external kinematical variables z_i; conflating the two is the identified error in the alternative definition. This machinery shows the 1/(4\\xi) prefactor comes from the parton momentum, while the phase-space Jacobians for the observed hadrons belong to the cross section rather than to the universal fragmentation function.","core_discovery":"The central claim is that for the process l^+l^- -> {h1,...,hn}+X with $M_h^{2}$ << $Q^{2}$, no new physics enters in the fragmentation function: the same light-cone operator that defines single-hadron fragmentation, with prefactor 1/(4\\xi), defines the n-hadron fragmentation function after inserting the n-hadron state in the matrix element. The \\xi prefactor arises from the fragmenting parton momentum, not from the hadronic phase space, so the factorization formula takes the standard form with the same partonic hard tensor and the same DGLAP kernels. Retracing the steps with a recently proposed nonuniversal prefactor 1/(\\xi1...\\xi n) shows that it must be accompanied by extra powers of \\xi in the hard part, and the paper demonstrates in a Yukawa model that such a modified prefactor changes the renormalization factor and therefore the evolution kernel. Consequences for a multiplicity sum rule are diagnosed as invalid outside the factorization region and in QCD.","pith_inferences":["A direct consequence of the paper's logic is that any phenomenological implementation that sets the parton momentum fraction equal to the observed hadron momentum fraction before factorization must be re-examined at higher orders, because that identification is only exact at zeroth order.","One testable extension is to compute the two-loop collinear fragmentation function for a two-hadron final state in a simple Yukawa theory under the modified prefactor definition and check whether its evolution kernel departs from the single-hadron DGLAP kernel; the paper only sketches this at one loop.","If the standard definition is correct, then extractions that use the same cross-section formula but cite different operator definitions should agree within uncertainties; the paper leaves that comparison as future work.","The argument suggests that exact multiplicity sum rules of the type used to motivate the modified prefactors should not be imposed as constraints on fragmentation functions, since such sum rules fail in QCD and outside the factorization region."],"forward_implications":["Existing extractions of dihadron fragmentation functions built on the standard definition remain consistent with collinear QCD factorization, with no correction factor applied to their hard parts.","The same DGLAP evolution kernels apply to n-hadron fragmentation functions as to single-hadron ones, because the renormalization factor does not depend on the observed hadron state.","Definitions that multiply the fragmentation operator by extra factors of \\xi would force those factors into the hard part and would alter the evolution kernels, as shown by the paper's model calculation.","When the cross section is written in variables such as M_h, \\zeta, or z1,z2,l_T, the accompanying Jacobian factors belong to the phase-space convention, not to the universal fragmentation function.","The same reasoning covers polarization-sensitive quantities such as interference fragmentation functions, since the argument does not depend on spin."],"supporting_citations":[{"why":"Supplies the canonical collinear-factorization derivation and the standard fragmentation function definition that the paper retraces for n hadrons.","marker":"[22]"},{"why":"Originates the operator definition with the usual prefactor that the paper argues is unchanged.","marker":"[23]"},{"why":"Proposes the alternative nonuniversal prefactor definition that the paper argues is inconsistent with factorization.","marker":"[29]"},{"why":"Defines the dihadron fragmentation function that the paper identifies with the standard prefactor and compares to the proposed modification.","marker":"[13]"},{"why":"Provides a standard dihadron fragmentation function definition used in past phenomenology whose consistency with factorization is reaffirmed.","marker":"[10]"},{"why":"Exemplifies the phenomenological SIA dihadron cross-section formula whose validity the paper confirms.","marker":"[4]"},{"why":"Supports the paper's identification of the origin of the discrepancy between standard and modified definitions.","marker":"[42]"},{"why":"Backs the argument that exact multiplicity sum rules of the type used to motivate modified prefactors fail in QCD.","marker":"[43]"}],"fun_headline_variants":["Same fragmentation operator works for hadron clusters","No new prefactors or kernels for multihadron fragmentation","Nonuniversal prefactor claims refuted for multihadron fragmentation","Multihadron final states keep standard fragmentation functions","QCD factorization retraced for multihadron clusters"],"cache_read_input_tokens":64512,"weakest_assumption_plain":"The load-bearing premise is that the simplified, non-gauge, one-flavour derivation extends without change to full QCD with gauge links and to all orders; the paper says this is straightforward and that the trend continues, but does not prove it explicitly.","fun_headline_variants_meta":{"raw":{"variants":["Same fragmentation operator works for hadron clusters","No new prefactors or kernels for multihadron fragmentation","Nonuniversal prefactor claims refuted for multihadron fragmentation","Multihadron final states keep standard fragmentation functions","QCD factorization retraced for multihadron clusters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000927,"raw_usage":{"total_tokens":4031,"prompt_tokens":1062,"completion_tokens":2969,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":678,"completion_tokens_details":{"reasoning_tokens":2891}},"tokens_in":678,"tokens_out":2969,"duration_ms":19550,"temperature":1.0,"reasoning_tokens":2891,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:14:27.220824+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete check would be to compute the two-loop virtual correction to the n-hadron fragmentation function in a Yukawa model and see whether the same ultraviolet counterterm that renormalizes the single-hadron function also cancels all collinear poles for n=2; if a hadron-dependent prefactor is required, the standard definition is not universal.","supporting_citations":[{"cited_title":"Flexible Parametrization of Generalized Parton Distributions: The Chiral-Odd Sector","cited_arxiv_id":"1311.0483","evidence_quote":"Exemplifies the phenomenological SIA dihadron cross-section formula whose validity the paper confirms."},{"cited_title":"Metz, Nucleon transversity and tensor charge from di-hadron production (2024), URL https://indico.cfnssbu","cited_arxiv_id":null,"evidence_quote":"Supports the paper's identification of the origin of the discrepancy between standard and modified definitions."},{"cited_title":"Pitonyak, Transversity from single-hadron tssas and dihadron fragmentation theory developments (2024), URL https: //agenda.infn.it/event/38132/contributions/234378/","cited_arxiv_id":null,"evidence_quote":"Backs the argument that exact multiplicity sum rules of the type used to motivate modified prefactors fail in QCD."}],"review_version":1}