{"id":"78aed9b8-378e-4d1e-9c2f-f7712fddb109","arxiv_id":"2507.15923","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In the post-confinement regime the small-angle energy-energy correlator scales with the collision energy through the J=5 DGLAP anomalous dimension, and the light-ray OPE coefficients are moments of the dihadron fragmentation function.","lead":"This paper explains how a two-particle energy correlation behaves after quarks and gluons turn into hadrons, using light-ray operators, and shows the same description follows from dihadron fragmentation functions. It matters because it links two formalisms and suggests a new way to measure the strong coupling constant.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The J=5 plateau scaling rests on the unquantified free-hadron/odd-JL assumption; the MC validation does not isolate this exponent.","rationale":"Agree with the reader's weakest assumption. The J=5 exponent is obtained by assuming only JL=-6-2n operators contribute to the hadronic light-ray OPE; this is the crux because any odd-JL operator would either change the leading plateau exponent or contaminate the plateau at finite ζQ^2. The paper is admirably explicit about the assumption and even labels Eq. (12) as conjectural, which prevents a stronger verdict. The proposed fit is a direct check of whether the free-hadron spectrum is sufficient inside the same Pythia framework used for validation; it is more decisive than comparing to real data because it removes detector and systematic issues from the first test. This does not move the verdict: a conditional accept remains appropriate. One secondary point the reader did not stress is the mismatch between the formal ζQ^2 much less than Λ_QCD^2 regime and the O(50) GeV^2 fitting region; this strengthens, rather than replaces, the free-hadron concern. No code is shipped, so the conditional status is appropriate.","tokens_in":15832,"tokens_out":27250,"duration_ms":318834,"concrete_test":"Re-analyze the Pythia 8.2 track-based F_{1,1} data of Fig. 4 with the extended plateau ansatz F(ζ,Q)=c_0(Q)+c_{1/2}(Q)(ζQ^2)^{1/2}+c_1(Q)(ζQ^2)+..., using LL evolution c_{1/2}∝α_s(Q)^{γ_6/β0}, c_0∝α_s(Q)^{γ_5/β0}, and c_1∝α_s(Q)^{γ_7/β0}. If a 2σ nonzero c_{1/2} emerges and materially changes the extracted c_0 evolution, odd-JL hadron operators are present and the J=5-only description of the plateau is incomplete. If c_{1/2} is consistent with zero and c_0 matches the γ_5 exponent, the free-hadron spectrum is supported within the simulation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central scaling law Eq. (10) is derived by keeping only the free-hadron double-twist operators [EE]_{n,0} with JL=-6-2n (Eq. 6). The paper itself flags the caveat: the Fig. 2 caption and the supplementary section state that hadron interactions can induce odd-JL detectors, and the text says hadron interactions can modify the spectrum of JL values appearing. No bound or power-counting estimate is given for such operators. If an odd-JL operator at JL=-7 is generated, its celestial block is ζ^{1/2}, so it is absent from the ζ→0 limit but contaminates the plateau at the finite ζQ^2 used in the validation, and it would evolve with γ_6, not γ_5. More importantly, the formal derivation applies for ζQ^2 much less than Λ_QCD^2 (stated before Eq. 5), whereas the Monte Carlo fits in Figs. 4-5 extend to ζQ^2 ~ O(50) GeV^2, far outside that regime for any conventional Λ_QCD. The MC agreement therefore does not isolate the J=5 exponent from interaction-induced background; it validates an extrapolated Laurent ansatz. The J=5 claim may be correct, but its decisive assumption is currently unquantified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a light-ray operator product expansion (OPE) formalism for the small-angle energy-energy correlator (EEC) in the post-confinement regime. The authors argue that in this regime the hadronic EFT product E(n1)E(n2) decomposes into double-twist operators with boost weights JL = −6−2n, which match onto twist-2 DGLAP detectors in perturbative QCD. This yields the central prediction, Eq. (10): the energy-scaled correlator F_{α,β}(ζQ^2,Q^2) is a sum over (ζQ^2)^n times coefficients whose Q-dependence is controlled by the DGLAP anomalous dimension at J = −JL−1 = 5+2n, so the plateau (n=0) scales with the J=5 anomalous dimension. The authors also derive a correspondence between the non-perturbative OPE coefficients and moments of the dihadron fragmentation function, Eq. (16). They validate the LL-evolved prediction against Pythia 8.2 for e+e− and pp collisions and for track-based measurements, fitting at Q=500 GeV and evolving to 200 and 1000 GeV. The paper is clearly written and honest about its conjectural elements, including the contour deformation in Eq. (12) and the existence of QCD light-ray Regge trajectories.","tokens_in":16032,"tokens_out":11772,"duration_ms":128527,"significance":"If correct, the J=5 scaling is a crisp, testable prediction for track-based EEC measurements at LEP, RHIC, and the LHC, and the OPE/DFF dictionary provides a novel bridge between light-ray OPE and standard factorization approaches to hadronization. The paper includes a concrete validation pipeline (fit at one energy, evolve, compare), and the LL evolution is implemented with standard kernels. However, the strength of the validation is limited: the central J=5 conclusion rests on the free-hadron spectrum assumption, which the manuscript itself flags as potentially modified by hadron interactions, and the Monte Carlo test does not isolate the exponent because the same generator is used for the fit and the validation.","major_comments":[{"comment":"The derivation of the J=5 scaling assumes that the only operators in the hadronic light-ray OPE have boost weights JL = −6−2n. The manuscript itself acknowledges that hadron interactions may induce odd-JL operators (Fig. 2, orange points) and can modify the spectrum, but it provides no power-counting estimate of such contributions. Since an odd-JL operator with JL = −7 would produce a ζ^{1/2} celestial block that is subleading as ζ→0 but contaminates the finite ζQ^2 region used in the fits, the Monte Carlo agreement does not establish that the J=5 exponent is the one controlling the plateau. Please either quantify the interaction-induced corrections or clearly present the J=5 claim as conditional on the free-hadron spectrum.","section":"Light-ray OPE analysis, Fig. 2 caption and text after Eq. (5); Eq. (10)"},{"comment":"The validation is a soft consistency check: the OPE coefficients are obtained by fitting Pythia data at Q=500 GeV and the predictions are compared with Pythia data at Q=200 and 1000 GeV. Because Pythia's parton shower and the LL evolution in Eq. (11) are based on the same DGLAP kernels, good agreement is expected a priori. To make the test sensitive to the claimed non-perturbative scaling, the authors should isolate the leading exponent, for example by extracting the slope of the n=0 coefficient c^{(0)}_{α,β}(Q) as a function of Q and comparing it with bγ^{(0)}_5/β0, or by including an independent data set or experimental measurement.","section":"Phenomenological studies, Fig. 4 and Eqs. (23)-(25)"},{"comment":"The Laurent ansatz is used far outside the formal domain of validity of Eq. (10), which is derived for ζQ^2 ≪ Λ_QCD^2; the fits extend to ζQ^2 ~ 100 GeV^2, and the paper itself notes that the convergence breaks down around ~50 GeV^2. The agreement in the transition region is therefore testing the conjectural contour-deformation formula Eq. (12) rather than the OPE derivation. Please state how the finite truncation order N=8 controls the error and how the transition-regime prediction is justified.","section":"SM 'Details of Series Approximation and Scaling Prediction', Eq. (23); main text after Fig. 4"},{"comment":"The hadron-to-pQCD detector matching in Eq. (4), and the assertion that the twist-2 DGLAP trajectory has the largest −Δ_L for each JL, is imported from the companion paper [43]. Because the entire JL spectrum and the J=5 conclusion depend on this step, the manuscript should summarize the argument in the main text or state explicitly that Eq. (4) is a result of [43] that is not proven here, so that the present paper is self-contained on this load-bearing point.","section":"Light-ray OPE analysis, Eq. (4) and Fig. 2"}],"minor_comments":[{"comment":"The variable ζ is used throughout the paper (e.g., Eqs. (5)-(8), Figs. 4-8) but is never explicitly defined; please define ζ = (1 − n1·n2)/2 ≈ θ^2/4 at first use.","section":"General notation"},{"comment":"The vector notation in Eq. (13), especially the definition of ⃗H = (2H_q, H_g) and the factor of 2 for quark and anti-quark contributions, should be explained more explicitly for readers not familiar with the factorization of Ref. [16].","section":"Eq. (13)"},{"comment":"In Eq. (22), the same symbol J^{(n)} is used both for the expansion coefficient of the jet function and for the angular momentum labels elsewhere in the paper; please use a distinct notation to avoid confusion with the spin variable J.","section":"SM Eq. (22)"},{"comment":"The sentence 'the invariant mass of the hadron pair is related to their angular separation by m^2/4 = z1z2Q^2ζ' would benefit from a short derivation or a reference, since this relation fixes the normalization of ζ and the factor 4 in Eq. (14).","section":"Text before Eq. (13)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a well-written Letter with a clear central claim. It depends heavily on the companion paper [43] for the hadron-to-pQCD detector matching; editors may want to verify the status of [43] for the claims to be self-contained. The comparison with Pythia is honest, but the use of the same generator for fit and validation limits the strength of the evidence; in my view this is not grounds for rejection, but the authors should either temper the 'excellent agreement' language or provide the exponent-isolation test. The conjectural nature of Eq. (12) should be clearly marked in the published version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nShort version: this paper gives a light-ray OPE treatment of the post-confinement EEC, connects the OPE coefficients to moments of the dihadron fragmentation function, and argues the plateau height scales with Q via the J=5 DGLAP anomalous dimension. The formalism is elegant and the dictionary (Eq. 16) is genuinely new. But the numerical validation is softer than the prose suggests, and the load-bearing free-hadron assumption is unquantified. I'd send it out, and I'd want the authors to test the J=5 prediction on real track data before treating it as established.\n\nWhat's new: the hadronic light-ray OPE at fixed boost weight, matched onto perturbative DGLAP detectors, is a clean framework. The identification of R_{-6-2n} with 4^n times the n-th m^2-derivative moment of the DFF is a nice bridge between CFT-style light-ray OPE and standard factorization. The contour formula (12) is honestly labeled as a conjecture. The Q-evolution predictions for F_{1,1} and F_{2,2} across e+e- and jet/pp settings, with a fit at one energy and LL DGLAP evolution to others, is the right way to present the test.\n\nSoft spots: (1) The formal OPE expansion (6) is derived for ζQ^2 << Λ_QCD^2, but the fits in Fig. 4 cover ζQ^2 up to ~50 GeV^2, which is far outside that regime for any reasonable Λ_QCD. So the MC agreement validates an extrapolated truncated Laurent ansatz, not the strict asymptotic OPE. The paper notes the series breaks down around that scale, but the reader should not come away thinking the strict asymptotic regime was tested. (2) The J=5 conclusion assumes only double-twist hadron operators with JL=-6-2n appear. The paper itself flags that hadron interactions can induce odd-JL operators, but gives no power-counting estimate. If an odd-JL operator at JL=-7 is present, it evolves with γ_6, not γ_5, and contaminates the plateau at finite ζQ^2. The MC test does not isolate the exponent. (3) The matching input (4) comes from the companion paper [43]; it is plausible but not re-derived here. (4) The DFF dictionary RG check is partly circular because the hard-function moments obey the same DGLAP equation by construction.\n\nThe paper is honest about most of these limitations in the text and supplemental material. Nothing I saw breaks the central scaling argument internally. But the evidence is a self-consistent check against Pythia, not an external falsification.\n\nBottom line: this deserves a serious referee. The right referee will push on the odd-JL power counting and on confronting the prediction with CMS/ALICE track-based EEC data. I'd suggest engaging with it.","headline":"Plausible and elegantly framed, but the J=5 claim is not yet isolated by the Monte Carlo check; worth a serious referee.","tokens_in":16743,"tokens_out":2727,"would_cite":true,"duration_ms":24835,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The post-confinement EEC plateau scales with Q according to the J=5 DGLAP anomalous dimension.","keywords":["energy-energy correlator","light-ray OPE","hadronization","DGLAP anomalous dimension","dihadron fragmentation function","post-confinement plateau","quantum scaling","alpha_s determination"],"falsifier":"A track-based measurement of the EEC plateau height at two different collision energies $Q_1,Q_2$ with the same $\\zeta Q^2$ would settle the claim: if the ratio of plateau heights disagrees with $[\\alpha_s(\\kappa^2 Q_2^2)/\\alpha_s(\\kappa^2 Q_1^2)]^{\\hat\\gamma_5^{(0)}/\\beta_0}$ (up to scale uncertainties), the $J=5$ quantum-scaling prediction is falsified. A model in which hadron interactions produce sizeable odd-$J_L$ operators would make the effective exponent differ from the $J=5$ value, which is also testable in the same comparison.","tokens_in":15481,"feed_emoji":"⚛️","tokens_out":8495,"duration_ms":83669,"temperature":0.7,"pith_summary":"The paper aims to explain the $Q$-scaling of the small-angle Energy-Energy Correlator (EEC) after hadronization, where the correlator is dominated by nearly free hadrons rather than partons. Using the light-ray operator product expansion in the effective field theory of hadrons, it derives that the post-confinement plateau of the EEC grows as $Q^2$ times a logarithmic factor set by the time-like DGLAP anomalous dimension at Lorentz spin $J=5$. The same derivation maps the non-perturbative OPE coefficients onto moments of the dihadron fragmentation function, linking two previously separate formalisms for hadronization. The leading-log prediction is checked against Monte Carlo simulations of $e^+e^-$ and $pp$ collisions, including track-based measurements, and the authors argue that precise measurements of this quantum scaling could improve determinations of $\\alpha_s$.","feed_headline":"Post-confinement plateau of the energy correlator scales via J=5 DGLAP","feed_subtitle":"A calculable power law ties the hadronization plateau to the J=5 DGLAP anomalous dimension.","key_machinery":"The machinery is the light-ray OPE for a pair of energy detectors, $E(n_1)E(n_2)=\\sum_{J_L,k} A_{J_L,k}\\, \\mathcal C_{J_L}(n_1,n_2,\\partial_{n_2})\\, O^H_{J_L,k}(n_2)$, where each hadronic light-ray operator $O^H_{J_L,k}$ has a definite boost weight $J_L$ and the celestial block behaves as $|n_{12}|^{-6-J_L}$. In the free-hadron approximation the relevant operators are double-twist operators $[EE]_{n,0}$ with $J_L=-6-2n$. Each such operator is matched, by Lorentz symmetry alone, onto the perturbative twist-2 DGLAP light-ray detector with the same boost weight, whose renormalization-group evolution is the time-like DGLAP anomalous dimension; this is what converts the non-perturbative hadron spectrum into the $J=5$ quantum-scaling prediction for the leading plateau term.","core_discovery":"The central claim is that in the post-confinement regime $\\zeta Q^2 \\ll \\Lambda_{\\rm QCD}^2$, the reduced correlator $F_{\\alpha,\\beta}(\\zeta Q^2,Q^2)$ obeys the light-ray OPE expansion of Eq. (10): it is a sum over $n$ of $(\\zeta Q^2)^n$ times matching coefficients $\\vec R^{\\alpha,\\beta}_n$ and the leading-log DGLAP evolution operator evaluated at boost weight $J_L=-\\alpha-\\beta-4-2n$. For the standard small-angle EEC with $\\alpha=\\beta=1$, the leading plateau term ($n=0$) has $J_L=-6$, so its $Q$-dependence is controlled by the DGLAP anomalous dimension at $J=5$; the plateau height therefore scales as $Q^2$ times a mild logarithmic factor. The paper further claims that the non-perturbative matching coefficients are moments of the mass-dependent dihadron fragmentation function, $\\vec R_{-6-2n}=4^n\\,\\vec J^{(n)}_{1,1}$, and that this dictionary is consistent with the known evolution of the DFF. The authors also conjecture a contour-integral formula that interpolates between this post-confinement expansion and the perturbative pre-confinement expansion.","pith_inferences":["Beyond the paper, the dictionary with DFF moments suggests that existing DFF parametrizations could be used as non-perturbative input to predict the overall normalization of the EEC plateau, leaving only the $Q$-dependent exponent as the perturbative prediction.","Beyond the paper, measuring the generalized correlators $F_{\\alpha,\\beta}$ for several distinct weights would give a family of scaling exponents indexed by $\\alpha+\\beta$; agreement among them would be a sharper test of the single-$J_L$ assignment than the plateau height alone.","Going further, if the conjectured contour formula of Eq. (12) is correct, then the power corrections of the pre-confinement region and the post-confinement expansion are two sides of the same Regge spectrum, so hadron-level information could in principle predict the size of non-perturbative corrections in the perturbative regime."],"forward_implications":["The plateau height of the post-confinement EEC carries a calculable $Q$-dependence, so measuring it at two energies (e.g., 200 and 1000 GeV) directly tests the $J=5$ DGLAP anomalous dimension in a regime usually thought of as purely non-perturbative.","The identity $\\vec R_{-6-2n}=4^n\\,\\vec J^{(n)}_{1,1}$ makes the small-angle energy correlator a direct experimental window into the mass-dependent dihadron fragmentation function.","At leading-log accuracy the evolved prediction matches Monte Carlo simulations of both $e^+e^-$ and $pp$ collisions, and it continues to hold for track-based measurements, which are the most practical way to reach the post-confinement plateau.","Because the evolution kernel depends on $\\alpha_s$ through $\\beta_0$ and the one-loop anomalous dimension, a percent-level measurement of the plateau scaling would translate into a competitive $\\alpha_s$ determination at colliders."],"supporting_citations":[{"why":"Supplies the intrinsic pQCD light-ray operators and the matching of hadron operators onto twist-2 DGLAP detectors at equal boost weight, the bridge used in Eq. (4).","marker":"[43]"},{"why":"Provides the perturbative EEC factorization and the time-like DGLAP anomalous dimension that sets the J=5 exponent.","marker":"[16]"},{"why":"Establishes the light-ray OPE treatment of the EEC and its scaling violation, the basis for the quantum-scaling analysis.","marker":"[18]"},{"why":"Gives the celestial-block expansion of light-ray OPE used to derive Eq. (2).","marker":"[28]"},{"why":"Supplies the known evolution of the dihadron fragmentation function at non-zero angular separation, used to verify Eq. (18).","marker":"[32]"},{"why":"Is the modern mass-dependent dihadron fragmentation function framework that the OPE coefficients are matched to.","marker":"[41]"},{"why":"The Monte Carlo event generator used to validate the leading-log predictions in $e^+e^-$ and $pp$ collisions.","marker":"[50]"},{"why":"Establishes track-based EEC measurements, which the paper uses for its proposed plateau-scaling measurement.","marker":"[52]"}],"fun_headline_variants":["Quantum scaling of energy correlator set by J=5 DGLAP","Plateau scaling pinned to J=5 DGLAP anomalous dimension","Energy correlator plateau reveals fragmentation function moments","New OPE-DFF link sets energy correlator scaling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction rests on the approximation that hadrons in the post-confinement regime are nearly free massless particles, so only double-twist operators with even boost weights $J_L=-6-2n$ appear in the hadron light-ray OPE; as the authors note, hadron interactions could induce odd-$J_L$ contributions and shift the leading scaling. It also assumes the matching of hadronic operators onto perturbative twist-2 DGLAP detectors at equal boost weight, a result taken from the companion paper [43] rather than proven here.","fun_headline_variants_meta":{"raw":{"variants":["Quantum scaling of energy correlator set by J=5 DGLAP","Plateau scaling pinned to J=5 DGLAP anomalous dimension","Energy correlator plateau reveals fragmentation function moments","New OPE-DFF link sets energy correlator scaling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000477,"raw_usage":{"total_tokens":2398,"prompt_tokens":1009,"completion_tokens":1389,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":625,"completion_tokens_details":{"reasoning_tokens":1319}},"tokens_in":625,"tokens_out":1389,"duration_ms":11591,"temperature":1.0,"reasoning_tokens":1319,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:24:23.191919+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A track-based measurement of the EEC plateau height at two different collision energies $Q_1,Q_2$ with the same $\\zeta Q^2$ would settle the claim: if the ratio of plateau heights disagrees with $[\\alpha_s(\\kappa^2 Q_2^2)/\\alpha_s(\\kappa^2 Q_1^2)]^{\\hat\\gamma_5^{(0)}/\\beta_0}$ (up to scale uncertainties), the $J=5$ quantum-scaling prediction is falsified. A model in which hadron interactions produce sizeable odd-$J_L$ operators would make the effective exponent differ from the $J=5$ value, which is also testable in the same comparison.","supporting_citations":[],"review_version":1}