{"id":"496fa26a-2108-4afe-b43e-2b60f9f106af","arxiv_id":"hep-ph/9902341","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Differential thrust and heavy-jet-mass distributions factorize into a convolution of a center-of-mass-energy-independent nonperturbative shape function with a resummed perturbative spectrum, organizing all 1/(eQ)^n power corrections.","lead":"In electron-positron collisions, jet-shape distributions near the two-jet limit receive nonperturbative corrections that this paper organizes into a single energy-independent shape function convolved with a perturbative prediction. The framework ties those corrections to universal energy-flow correlations, giving a template for extracting the strong coupling from jet-shape data.","discovery_kind":"unclear","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"Eqs. 26–28 and 37 appear sign-inconsistent: with λ2 = ∫₀^{µ²}dk² Γ > 0, the variance of f_T is −λ2 < 0, which no probability distribution can have and which contradicts the paper's own ansatz (Eqs. 44–45).","rationale":"The reader's weakest_assumption — the factorization premise and the eikonal replacement, with the Δf non-factorization of Eq. 13 — is a fair reading, but the paper acknowledges it explicitly: Eq. 13 states that f(εL,εR) does not factorize exactly, and Eq. 14 restricts the product form to NLL. The heavy-mass convolution Eq. 15/21 is therefore an NLL approximation, and the paper says so. That is an acknowledged limitation, not a new finding, so I do not re-present it as the load-bearing concern. The unacknowledged issue I find is the sign structure of the moment expansion: Eq. 26 defines λ2 as a positive integral, while Eq. 28 and Eq. 37 assign the variance the value −λ2, and the paper's own positive ansatz in Eqs. 44–45 requires the opposite sign. This is a genuine internal inconsistency in the 'all n' claim of abstract:s3 and in the connected-correlator relation Eq. 37, but it is a technical sign issue that is straightforward to repair and does not undermine the leading 1/Q shift, the convolution form of Eq. 19, or the Q-independence of f_T, which are the paper's central contributions. The reader's CONDITIONAL verdict already reflects a framework that is well-motivated but not a proof, with a non-reproducible two-parameter fit in Fig. 1; my concern adds a specific checkable correction rather than changing the verdict category. I therefore recommend UNCHANGED, with the sign check as a concrete verification step.","tokens_in":15091,"tokens_out":18960,"duration_ms":347878,"concrete_test":"Recompute the second moment directly from Eq. 27: expand exp(−S_NP) to O(ν²) using the literal Eq. 26, and compare the coefficient with Eq. 28's λ1² − λ2 and with the variance of the ansatz in Eqs. 44–45. If the literal Eq. 26 yields a negative variance, the sign in Eq. 26 (or in Eq. 28/37) needs correction; then check whether the corrected relation changes the identification of the connected correlator G(n1,n2) − G(n1)G(n2) in Eq. 37 and the heavy-mass variance in Eq. 38.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's 'all n' organization of power corrections (abstract:s3) is realized through the moment relations of Eqs. 26–28. Eq. 26 expands the nonperturbative exponent as S_NP = Σ_{n>0} (ν/Q)^n λ_n/n! with λ_n(µ) = ∫₀^{µ²} dk² k^{n−2} Γ(α_s(k²)). For n=2 this is λ2 = ∫₀^{µ²} dk² Γ, which is strictly positive because the cusp anomalous dimension Γ is positive. Eq. 27 defines f_T as the inverse Laplace transform of exp(−S_NP), and Eq. 28 then gives ∫dε ε² f_T = λ1² − λ2, so the variance is ⟨ε²⟩ − ⟨ε⟩² = −λ2 < 0. A probability distribution cannot have a negative variance. The paper's own two-parameter ansatz, Eq. 44, is a positive function whose moments satisfy Eq. 45, Λ²α = λ1² − λ2, and whose variance is Λ²[α − Γ(α+1/2)²/Γ(α)²] = −λ2 > 0; this requires λ2 < 0, the opposite of Eq. 26. The likely resolution is a missing (−1)^{n+1} in Eq. 26, since 1 − e^{−να} = Σ (−1)^{n+1}(να)^n/n!; with that factor Eq. 28 would read λ1² + λ2 and Eq. 37 would read +λ2, restoring positivity. As written, the n=2 moment and the connected-correlator identification in Eq. 37 carry the wrong sign. This is a concrete, checkable internal inconsistency in the 'all n' claim; it does not by itself challenge the leading 1/Q shift or the factorization/convolution structure, which are the paper's main contributions.","agreement_with_reader":"partial"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The paper proposes that the differential thrust (t = 1−T) and heavy-jet-mass (ρ) distributions of e+e− annihilation in the two-jet limit can be written as convolutions of Q-independent nonperturbative 'infrared shape functions' with resummed perturbative spectra. Section 2 motivates the convolution kinematically, isolating the soft light-cone momentum variables ε_R, ε_L of Eq. (8) and writing the radiation functions of Eqs. (9)–(10) as integrals over a common soft distribution f(ε_R, ε_L); the thrust reduces to a one-dimensional convolution with f_T(ε) (Eqs. (11)–(12), (17), (19)), while the heavy mass requires a two-hemisphere convolution (Eqs. (15), (21)). Section 3 rederives the same structure from the exponentiated Laplace-transformed cross section (Eqs. (24)–(25)), splitting the exponent at a factorization scale µ and identifying the nonperturbative remainder S_NP with the Laplace transform of f_T (Eq. (27)), which yields the moment sum rules of Eq. (28). Section 4 gives operator definitions of the shape functions as Wilson-line matrix elements of an energy-flow operator (Eqs. (31)–(32)), shows that their moments are universal correlators G(n̂_1…n̂_N) (Eqs. (34)–(35)), relates the leading 1/Q correction to the energy-energy correlation (Eq. (40)), and exhibits a two-parameter ansatz (Eqs. (44)–(45)) compared with thrust data from 14 to 161 GeV in Fig. 1.","tokens_in":15459,"tokens_out":15782,"duration_ms":286029,"significance":"If the factorization holds, the paper supplies a framework rather than a model: it organizes the entire series of 1/(eQ)^n corrections, not just the leading shift of Eq. (1), and it gives the shape functions an operator definition (Eqs. (31)–(32)) that makes the Q-independence claim falsifiable and connects event-shape power corrections to energy-flow correlators and to the EEC (Eq. (40)). The moment sum rules of Eqs. (28) and (37) and the explicit ansatz of Eqs. (44)–(45) turn the formalism into something that can be, and has been, confronted with data at several energies. I checked the normalization of Eq. (43) against Eq. (34) with the weight of Eq. (33) and it is correct, and I verified that the moments of the ansatz in Eq. (44) reproduce both relations of Eq. (45). These are genuine strengths: a parameter-free structural prediction (Q-independence of f_T), a worked phenomenological implementation, and a clean field-theoretic interpretation.","major_comments":[{"comment":"The heavy-mass convolution of Eq. (15) retains only the factorized part of the two-hemisphere shape function, f(ε_L,ε_R) ≈ f(ε_R)f(ε_L), while the paper itself states in Eq. (13) that the exact distribution contains a nonfactorizing term Δf, attributed to the splitting of an off-shell gluon into gluons in opposite hemispheres. The justification offered is that the perturbative expectation values factorize at NLL (Eq. (14), cited to Ref. [16]). That justifies the perturbative factor, but it does not establish that Δf is subleading at the same power in 1/(ρQ) as the terms kept in Eq. (21). Since the shape function dominates the distribution precisely for ρ ≲ µ/Q, where the neglected term is largest, the leading-power accuracy of Eq. (21) rests on an unstated bound on Δf. The paper should either estimate Δf or state explicitly that Eq. (21) is accurate only up to a nonfactorizing correction of the same parametric order, which would qualify the claim that the heavy-mass distribution is described by the same universal f.","section":"§2, Eqs. (13)–(15) and Eq. (21)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The abstract states that the infrared shape functions are independent of the center-of-mass energy Q, but the text states of Eq. (22) that 'in contrast with f_T, the shape function f_H depends both on the center-of-mass energy and the shape variable.' The two statements are reconcilable only because the Q-dependence of f_H enters through the upper limit ρQ and through R_J^{PT}, and is power-suppressed, as the expansion of Eq. (23) shows explicitly. As written, however, the abstract's blanket Q-independence claim is contradicted by the paper's own sentence about f_H. The abstract or the sentence in §2 should be amended to say that f_H is Q-independent up to corrections of order 1/(ρQ), so that the scope of the Q-independence claim is unambiguous.","section":"§2, Eqs. (22)–(23), versus abstract:s2","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The definition of the nonperturbative parameters λ_n(µ) should be written out explicitly, including the alternating sign and the 1/n factor carried by the intermediate expression of Eq. (26). The intermediate line contains (−ν/Q)^n/(n·n!), whereas the final line is written as (ν/Q)^n/n! λ_n(µ); the '≡' therefore defines λ_n = (−1)^{n+1} n^{−1} ∫₀^{µ²} dk⊥² k⊥^{n−2} Γ, not the integral alone. This matters because Eq. (28) then gives a variance ⟨ε²⟩−⟨ε⟩² = −λ2, which requires λ2 < 0, and Eq. (45) requires the same sign through Λ²α = λ1² − λ2. A reader who takes λ_n to be the bare integral (which is positive for n = 2, since Γ > 0) will conclude that f_T has a negative variance, which is impossible for the positive ansatz of Eq. (44). I do not believe this is an internal inconsistency: the sign is carried by the intermediate line, and Eqs. (28), (37) and (45) are mutually consistent once λ2 < 0. But the equation as printed invites the wrong reading, and since Eq. (26) is the vehicle for the 'all n' claim of abstract:s3, the definition should be displayed unambiguously.","section":"§3, Eq. (26)","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The overall factor in the definition of R_T^{PT}(t′,µ) is written θ(t) where it should be θ(t′); the argument of the step function should match the variable of the function being defined.","section":"§2, Eq. (18)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The statement that R_J^{PT}(0,µ) = R_J^{PT}(µ/Q,µ) would be easier to verify if the two arguments were shown explicitly in the integral, since the equality relies on the max(ρ,µ/Q) lower limit freezing the ρ-dependence below ρ = µ/Q.","section":"§2, Eq. (16)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The notation ⟨…⟩_PT is used both for an average over the full perturbative spectrum and, inside the ε_R,ε_L integrals, for an average over the double differential distribution in the two hemisphere masses. A sentence distinguishing the two usages, and stating that the perturbative factor in Eq. (15) is the double distribution before the NLL factorization of Eq. (14) is applied, would remove the ambiguity.","section":"§2, Eqs. (9)–(10)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The phrase 'for e^{1/2} ≫ Λ_QCD/Q ∼ e' is hard to parse; presumably the intended statement is that the hierarchy Qe^{1/2} ≫ Qe ≫ Λ_QCD defines the region in which only powers of 1/(Qe) need be retained. Rewording would help.","section":"§2, paragraph after Eq. (7)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The figure is reproduced from Ref. [10] and its caption defers the details to that reference. For a journal article the caption should be self-contained: state what 'scaled' means, what the vertical normalization is, and what the curves correspond to, so that the comparison across Q = 14–161 GeV can be read without consulting the proceedings.","section":"Fig. 1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"There are several typographical errors, among them 'extentions' in the last paragraph of §3 and 'truely' in §4; these should be corrected in the final version.","section":"§3, final paragraph, and §4","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is a letter-length exposition whose only quantitative confrontation with data (Fig. 1) is imported from the author's conference proceedings [10]. The physics is sound and the framework has been influential, but the editor may wish to confirm that reproduction of the figure is permissible and that the article is sufficiently self-contained for a journal audience, given how much of the phenomenological support lives in the cited proceedings."},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"This is the paper that put the infrared shape function on the map for event shapes, and the framework is the reason to read it. The convolution form for the differential thrust and heavy-jet-mass distributions (Eqs. 19, 21), the operator definition as a Wilson-line matrix element of an energy-flow operator (Eqs. 31, 32), and the tie to energy-energy correlations (Eq. 40) are the conceptual core. That core is right and it has aged well. The definitional element the reader flags is real — f_T is defined by Eq. 27 as the inverse Laplace transform of S_NP — but the content is the Q-independence of S_NP and the fact that the same f_T convolves with the full perturbative spectrum, not just the exponentiated form. That is more than a tautology, though it is asserted rather than proved in detail here.","headline":"The shape-function convolution framework is the real contribution and it holds up, but the moment relations in Eqs. 26–28 and 37 carry a sign error a referee should catch.","tokens_in":16147,"tokens_out":6630,"would_cite":true,"duration_ms":120508,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":["12.38.-t","13.87.-a"],"model":"selfhosted/deepseek-v4.1-flash","headline":"The paper argues that near the two-jet limit the thrust and heavy-jet-mass distributions factor into a convolution of one energy-independent soft shape function with a resummed perturbative spectrum, so that every power correction of the形式$","keywords":["event shapes","power corrections","factorization","infrared shape function","thrust","heavy jet mass","energy-energy correlation","Wilson lines"],"falsifier":"Fit the thrust shape function to data at $Q = M_Z$ and check whether it reproduces the thrust distribution at 14, 35, 91 and 161 GeV without refitting; separately test the parameter-free relations that the nonperturbative shift in $����$ is half that in $��t��$ and that the same $��_1$ fixes the $1/Q$ energy-energy correlation through $G(��) = ��_1/(4�� ��^3 ��)$.","tokens_in":16357,"feed_emoji":"⚛️","tokens_out":4180,"duration_ms":78469,"temperature":0.7,"pith_summary":"Near the two-jet limit of electron-positron annihilation, the thrust and heavy-jet-mass distributions are dominated by soft radiation carrying a small fraction of the beam energy. The paper argues that this soft radiation separates from the evolution of the narrow jets and can be packaged into a shape function: a probability distribution for how much light-cone energy the soft quanta deposit in each hemisphere. If that holds, the measured distribution is a convolution of this nonperturbative function with a perturbatively resummed spectrum, and every power correction of the form $1/(eQ)^n$ is generated by the same function rather than by a separate parameter for each $n$. The payoff is predictive: the shape function is fixed once, at a single energy, and must then describe the same physics at every energy, while its moments are vacuum expectation values of energy-flow operators that also control the energy-energy correlation.","feed_headline":"Jet-mass power corrections collapse into one shape function","feed_subtitle":"A single energy-flow function, fixed at one beam energy, is asked to predict thrust spectra from 14 to 161 GeV.","key_machinery":"The shape function itself. For thrust it is defined by the Laplace transform $��_0^�� d�� e^{-�ͧ�/Q} f_T(��;��) = ��(-S_{NP}(Q/��,��))$, which ties its moments to the ambiguous low-momentum integrals of the resummed exponent, and equivalently by the operator expression $f_T(��) = ��0|W^��(0) ��(�� - �� d�� n w_T(n) E(n)) W(0)|0��$, where $W$ is the product of eikonal Wilson lines of the two jets and $E(n)$ is an energy-flow operator measuring energy density in direction $n$, with weight $w_T(n) = 1-|��cos��|$. The weight vanishes along the jet directions, which is why the shape function has no collinear singularities. Taking moments of $f_T$ turns the weights into integrals of Green functions $G(n_1��n_N) = ��0|W^�� E(n_1)��E(n_N) W|0��$; these are the universal objects, independent of which event shape was chosen, and the single-energy-flow correlator $G(��)$ is what the leading $1/Q$ energy-energy correlation measures.","core_discovery":"The central claim is that the differential thrust and heavy-jet-mass distributions in the two-jet limit, $e = 1-T$ or $e = M^2/Q^2$, can be written as convolutions of a nonperturbative infrared shape function with a resummed perturbative cross section. The shape function depends on a factorization scale $��$ that separates soft from hard dynamics, and differs from one event shape to another, but is independent of the center-of-mass energy $Q$. The authors derive this twice: once from the kinematics of soft emission, writing the hemisphere masses as perturbative pieces plus soft light-cone sums $M^2_{R,L}/Q^2 + ��_{R,L}/Q$ and factorizing the soft distribution $f(��_R,��_L)$; and once from the resummed exponent, where the nonperturbative part of the Sudakov exponent $S_{NP} = ��_{n>0} (��/Q)^n ��_n(��)/n!$ is the Laplace transform of the thrust shape function $f_T(��)$, so that the moments of $f_T$ are exactly the $��_n$. The familiar shift $d��_{PT}(t)/dt �� d��_{PT}(t-��_1/Q)/dt$ reappears as the special case in which $f_T$ is a delta function, and the same limit gives half that shift for the heavy-jet mass, a parameter-free relation between the two observables. Beyond the leading shift, the shape function organizes all powers $1/(tQ)^n$ and $1/(��Q)^n$, and its variance controls the smearing that the simple shift misses.","pith_inferences":["If the shape function really carries no $Q$ dependence, a global fit across energies is over-determined, and the drift of fitted $��_1$ values with energy becomes a direct handle on the size of the neglected non-factorizing term $��f$ and of higher $1/Q$ corrections.","The Green functions are singular only where two directions coincide, and stereographic projection maps them onto a plane; that suggests the energy-flow correlators might be computable in a two-dimensional effective description, which the paper only gestures toward.","Because the construction needs only an additive weight on soft momenta that vanishes along the jets, the same convolution should exist for any infrared-safe observable of that class, making the shape-function idea as broad as the class of such weights.","The same structure should carry over to hadron-collider event shapes and jet masses, where soft radiation is emitted by more than two coloured legs; whether the shape function stays universal there is an open and testable question."],"forward_implications":["A shape function fitted to thrust data at one energy predicts the thrust distribution at all higher energies with no new parameters.","All corrections of the form $1/(tQ)^n$ follow from the moments of one function, with the simple shift of the perturbative spectrum as the delta-function limit.","The nonperturbative shift of the mean heavy-jet mass is exactly half that of the mean thrust, a relation that holds for any shape function.","Moments of the shape function are universal energy-flow matrix elements, so event-shape measurements constrain the same quantities that appear in energy-energy correlations.","The variance of the thrust distribution receives no $1/Q$ correction; its leading nonperturbative piece is $1/Q^2$ and is tied to connected two-point energy-flow correlations."],"supporting_citations":[{"why":"Supplies the resummed next-to-leading-logarithmic thrust and jet-mass spectra, and the exponentiation under Laplace transform, on which the convolution is built.","marker":"[16]"},{"why":"Supplies the eikonal Wilson-line description of soft radiation from the two jets and the leading $1/Q$ thrust shift that the shape function generalizes.","marker":"[4]"},{"why":"Supplies the argument that perturbative ambiguities fix the exponents of the leading power corrections to event-shape means.","marker":"[3]"},{"why":"Supplies the effective-coupling moment parameterization that the shape-function moments are meant to replace with a more general object.","marker":"[8]"},{"why":"Supplies the inclusive light-cone shape function for semileptonic $b$ decays, the template the event-shape shape function is modelled on.","marker":"[11]"},{"why":"Supplies the representation of the energy-flow operator at infinity in terms of the energy-momentum tensor.","marker":"[12]"},{"why":"Supplies the energy-flow correlation functions on a sphere of infinite radius used to define the universal Green functions.","marker":"[13]"},{"why":"Supplies the factorization review underlying the separation of jet and soft dynamics invoked in Sec. 2.","marker":"[17]"},{"why":"Supplies the phenomenological fit of the shape function to thrust data at several energies shown in Fig. 1.","marker":"[10]"}],"fun_headline_variants":["Thrust shape function absorbs all jet-mass power corrections","One shape function links thrust and heavy-jet-mass shifts","Soft shape function organizes every 1/(eQ)^n correction","Shape-function convolution resums thrust power corrections at any Q","Delta-function limit yields a parameter-free thrust–mass relation"],"cache_read_input_tokens":0,"weakest_assumption_plain":"Soft gluon emission from the two jets is incoherent and depends only on the jet directions and colour charges, so that the soft distribution factorizes from jet evolution and carries no dependence on the beam energy; the authors note this fails at the level of the $��f$ term from a gluon splitting into opposite hemispheres.","fun_headline_variants_meta":{"raw":{"variants":["Thrust shape function absorbs all jet-mass power corrections","One shape function links thrust and heavy-jet-mass shifts","Soft shape function organizes every 1/(eQ)^n correction","Shape-function convolution resums thrust power corrections at any Q","Delta-function limit yields a parameter-free thrust–mass relation"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":1536,"prompt_tokens":1033,"completion_tokens":503,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":434}},"tokens_in":1033,"tokens_out":503,"duration_ms":8829,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":434,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-29T18:34:41.323259+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"falsifier":"Fit the thrust shape function to data at $Q = M_Z$ and check whether it reproduces the thrust distribution at 14, 35, 91 and 161 GeV without refitting; separately test the parameter-free relations that the nonperturbative shift in $����$ is half that in $��t��$ and that the same $��_1$ fixes the $1/Q$ energy-energy correlation through $G(��) = ��_1/(4�� ��^3 ��)$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the resummed next-to-leading-logarithmic thrust and jet-mass spectra, and the exponentiation under Laplace transform, on which the convolution is built."},{"cited_title":"Korchemsky and G.Sterman, Nucl","cited_arxiv_id":null,"evidence_quote":"Supplies the eikonal Wilson-line description of soft radiation from the two jets and the leading $1/Q$ thrust shift that the shape function generalizes."},{"cited_title":"Webber, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the argument that perturbative ambiguities fix the exponents of the leading power corrections to event-shape means."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the effective-coupling moment parameterization that the shape-function moments are meant to replace with a more general object."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the inclusive light-cone shape function for semileptonic $b$ decays, the template the event-shape shape function is modelled on."},{"cited_title":"Jet observables and energy-momentum tensor","cited_arxiv_id":"hep-ph/9710349","evidence_quote":"Supplies the representation of the energy-flow operator at infinity in terms of the energy-momentum tensor."},{"cited_title":"Power Corrections and Nonlocal Operators","cited_arxiv_id":"hep-ph/9708346","evidence_quote":"Supplies the energy-flow correlation functions on a sphere of infinite radius used to define the universal Green functions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the factorization review underlying the separation of jet and soft dynamics invoked in Sec. 2."},{"cited_title":"Korchemsky, in proceedings of the 33th Rencontres de Mo riond, QCD and High Energy Hadronic Interactions, Les Arcs, Savoie, France, March 21-28, 1998; hep-ph/980653 7","cited_arxiv_id":null,"evidence_quote":"Supplies the phenomenological fit of the shape function to thrust data at several energies shown in Fig. 1."}],"review_version":1}