Characterising hadronisation across the phase space

Basic Information

• Paper ID: 2511.17708
• Title: Characterising hadronisation across the phase space
• Authors: Andrea Banfi (University of Sussex), Basem Kamal El-Menoufi (Monash University)
• Categories: hep-ph (High Energy Physics - Phenomenology), hep-th (High Energy Physics - Theory)
• Submission Date: November 21, 2025
• Paper Link: https://arxiv.org/abs/2511.17708

Abstract

This paper investigates leading-order hadronisation corrections to global event shape variables in e+e− annihilation. Existing theory indicates that hadronisation corrections manifest as shifts in the perturbative distribution, with the shift strongly dependent on the event shape variable value. Previous studies primarily considered three-particle final states in the three-jet region or only two hard emitters in the two-jet region. The authors discover that in the two-jet region near the Sudakov peak, conventional methods miss important corrections from additional soft wide-angle gluon emission. They calculate this contribution and establish a unified theoretical framework valid in both two-jet and three-jet regions.

Research Background and Motivation

Core Problems to Address

1. Phase-space dependence of non-perturbative corrections: Hadronisation corrections are not constant but strongly depend on the event shape variable v
2. Unified description of two-jet and three-jet regions: Existing theory applies either strictly to the two-jet limit (τ→0) or only to the three-jet region, lacking a unified framework
3. Missing soft wide-angle gluon contributions: Traditional two-jet limit calculations only consider soft collinear emission, neglecting important soft wide-angle gluon effects

Importance of the Problem

• Precise determination of strong coupling αs: LEP data analysis shows deviations up to 3 standard deviations between αs values extracted from different event shapes and the world average. Accurate understanding of hadronisation corrections is crucial for eliminating these discrepancies
• Theoretical consistency: Need to provide consistent theoretical predictions across the entire phase space, particularly near the peak of event shape distributions (τ≃0.025 at LEP1 energy)
• Understanding QCD non-perturbative dynamics: Hadronisation is one of the most fundamental yet least understood phenomena in QCD

Limitations of Existing Methods

1. Strict two-jet limit: Traditional methods based on the dispersive model only consider qq̄ dipoles plus arbitrary soft collinear gluons, predicting constant non-perturbative shift ζV(0)
2. Fixed-order three-jet calculations: Recent work 22,23 computed the shift ζV(v) for three-particle final states at O(αs) precision, finding strong v-dependence, but only applicable to the three-jet region
3. Missing O(αs) constant terms: Two-jet limit calculations only provide αs ln(1/τ) terms, unable to determine O(αs) constants, making corrections to ζV(τ)/ζV(0) uncontrollable

Core Contributions

1. Discovery and calculation of soft wide-angle gluon contributions: Proves that additional soft wide-angle gluon emission produces new contributions proportional to CA, with O(αs) constant term 16(ln2-1)CA
2. Establishment of anomalous dimension γ: Discovers anomalous dimension evolution of non-perturbative shift, γ=4(ln2-1), and provides its resummation expression
3. Construction of unified framework: Proposes general formulas (equations 5.8 and 5.23) valid in both two-jet and three-jet regions, achieving smooth interpolation
4. Clarification of O(αs) correction color structure:
   • Hard collinear contribution: ∝CF, coefficient -3
   • Soft wide-angle contribution: ∝CA, coefficient 16(ln2-1)
5. Provision of matching scheme: Provides specific scheme for matching resummation results with fixed-order calculations

Methodology Details

Task Definition

Study non-perturbative corrections to event shape variable V in e+e−→hadrons process. Define:

• Event shape variable: V({pi}), function of all final-state particle momenta
• Thrust variable: T = (1/Q)max_{n̂T}Σi|p⃗i·n̂T|, τ=1-T
• Non-perturbative shift: ⟨δVNP⟩, such that cumulative distribution satisfies δΣNP(v) = -⟨δVNP⟩∫dZ{ki}δ(v-V({p̃},{ki}))

Theoretical Framework

1. Gluer Model Foundation

Non-perturbative effects modeled through emission of a "gluer" gluon with transverse momentum far below the perturbative scale:

[dkus]M²NP(kus) = Σ_dipoles C_dip (dκus/κus)M²NP(κus)dηus(dϕus/2π)

where color factors: Cqq̄ = 1-CA/(2CF), Cqg = Cgq̄ = CA/(2CF)

2. Core Function hV

Define key quantity for event shape change:

hV(ηus,ϕus,{p̃},k) = lim_{κus→0} (1/κus)[V({p̃},kus,{ki}) - V({p̃},{ki})]

For thrust variable, explicit form:

hT = lim_{κus→0} (1/κus)(-|k⃗us·n⃗T| + |k⃗+us·n⃗T| + |k⃗-us·n⃗T|)
    = -|Ae^ηus + Be^{-ηus} + Ccosϕus| + |A|e^ηus + |B|e^{-ηus}

3. Decomposition of Two-Jet Limit

Key innovation: decompose O(αs) contribution as:

ζ̄ij(τ) = (CFαs/2π)[c1 ln(1/τ) + c0 + c(τ)]

where:

• c1: from soft collinear region
• c0: hard collinear contribution (∝CF) + soft wide-angle contribution (∝CA)
• c(τ): power-suppressed terms, vanishing as τ→0

4. Parametrization of Soft Wide-Angle Region

Using transverse momentum-rapidity parametrization:

k = (kt/√sqq̄)e^η pq + (kt/√sqq̄)e^{-η} pq̄ + k⊥
1-x1 = (kt/Q)e^η,  1-x2 = (kt/Q)e^{-η}

Define rapidity function:

fT(η) = gT(η) - 2Θ(η<0)

whose integral gives anomalous dimension:

γ = ∫_{-∞}^{∞} dη fT(η) = 4(ln2-1)

Technical Innovations

1. Identification of Soft Wide-Angle Contribution

Discovery that in soft limit kt→0, individual functions A,B,C may diverge, but combination ⟨hij/n̂T_τ⟩ remains finite:

A^(qg)_q = (1/2)√(Q/kt)e^{η/2}  (divergent)
but |AB| + C² combination is finite

2. Discovery of Anomalous Dimension

Through analysis of contributions at different momentum scales kt, identify two regions:

• kt ~ τQ: gives finite O(αs) constant
• kt ≪ τQ: produces logarithmic enhancement ln(τQ/κus), leading to anomalous dimension

Final form:

ζ̄T(τ,κus) = ζT(τ)exp[2CAγ∫_{κus}^Q (dkt/kt)(αs(kt)/π)]

3. Combination with Soft Collinear Resummation

Embed soft wide-angle gluons into soft collinear ensemble:

v(dσ/dv)ζ̄V = ∫dZsc[{ki}](lim_{k→0}⟨h^{qq̄}_V⟩)vδ(v-Vsc)
             + (CA/2CF)∫_{κus}[dk]M²(k)∫dZsc[⟨h^{qg}_V+h^{gq̄}_V-h^{qq̄}_V⟩]×...

4. Matching Function χT(R')

Effect of soft wide-angle gluons in soft collinear ensemble encoded by function:

χT(R') = ∫_0^1 (dζ/ζ)[(1-ζ)^{-1+R'} - 1] = 1/R' + d/dR' lnF(R')

where R' = -τdR/dτ ≃ (2αs/π)CF ln(1/τ) is resummation parameter.

Experimental Setup

Theoretical Calculation Setup

This is pure theoretical work without experimental data fitting, but setup includes:

1. Energy scale: Q = MZ = 91.2 GeV (LEP1 energy)
2. Strong coupling constant: αs(MZ) = 0.118
3. Resummation precision: NLL (next-to-leading logarithm)
4. Matching scheme: additive matching, including O(αs) hard collinear constant

Numerical Calculation Methods

1. Phase space integration: parametrization using Dalitz variables (x1,x2), xi = 2Ei/Q
2. Rapidity integration: numerical calculation of complete elliptic integrals K and E
3. Resummation switch: use modified logarithm to avoid kinematic endpoint divergence

   ln(1/τ) → (1/p)ln(1 + 1/τ^p - 1/τ^p_max)
   

   where p∈1/2, 4/3 used for theoretical uncertainty estimation

Comparison Schemes

1. Fixed-order three-jet results: O(αs) calculation from ref.23
2. Naive two-jet limit: ζT(0) = 2 (constant)
3. Matched resummation results: complete formula from this work

Experimental Results

Main Numerical Results

1. qq̄ Dipole Contribution

Figure 1 shows ζqq̄(τ)/ζT(0):

• approximately 0.4 at τ=0.3
• rapidly approaches 1 on logarithmic scale
• very close to two-jet limit for τ<10^{-2}

Decomposed as:

ζ̄qq̄(τ) = (αsCF/2π)[2(4ln(1/τ) - 3) + cqq̄(τ)]

where cqq̄(τ) rapidly vanishes as τ→0 (Figure 3).

2. qg Dipole Contribution

Figure 2 shows ζqg(τ)/ζT(0):

• approximately 0.33 at τ=0.3
• asymptotically approaches 1/2 but converges very slowly
• still significantly deviates from limit value at τ=10^{-3}

Decomposed as:

ζ̄qg(τ) = (αsCF/2π)(4ln(1/τ) - 3 + 16(ln2-1) + cqg(τ))

Key finding: constant term 16(ln2-1)≈2.09 originates from soft wide-angle region

3. Numerical Verification of Anomalous Dimension

Integration of rapidity function fT(η) (Appendix B):

γ = ∫_{-∞}^{∞} dη fT(η) = 4(ln2-1) ≈ 1.227

This leads to anomalous dimension exponent:

CAγ/(πβ0ζT(0)) = CA/(πβ0) × 2(ln2-1) ≈ 0.134 (Nc=3, nf=5)

Complete Matching Results

Figure 7 displays three calculations of ζT(τ):

1. Fixed-order (blue): rapidly decreases from 2 to approximately 1.6 for τ<0.1
2. Matched resummation (red band):
   • 10-15% below fixed-order in range τ∈0.05,0.3
   • rapidly increases for τ<0.03 due to χT(R')∼-lnR' and αs(τQ) growth
   • uncertainty band width approximately ±5%
3. Naive two-jet limit (horizontal dashed line): ζT(0)=2

Key observations:

• In typical fitting range 0.05<τ<0.3, resummation pushes fixed-order results toward two-jet limit
• Theory fails near peak region τ≃0.025 (αs(τQ) approaches Landau pole)

Ablation Analysis

Relative Importance of Different Contributions

In τ→0 limit:

ζT(τ) ∝ (αs(Q)/αs(τQ))^{0.134} × [base term + CA soft wide-angle correction + CF hard collinear correction]

1. Anomalous dimension factor: (αs(Q)/αs(τQ))^{0.134} approximately 1.02 at τ=0.1, mild effect
2. CA soft wide-angle constant: 16(ln2-1)≈2.09, approximately 17% relative to logarithmic term 4ln(1/τ) at τ=0.1
3. CF hard collinear constant: -3, approximately -13% of logarithmic term

Related Work

Historical Development

1. Foundational work in 1990s:
   • Dokshitzer-Webber 6,10: established gluer model and dispersive method
   • Akhoury-Zakharov 7: renormalon method
   • Main conclusion: linear power corrections ∝ΛQCD/Q most important
2. Precision measurements in 2010s:
   • Gehrmann et al. 17: αs extraction from thrust distribution
   • Hoang et al. 18: αs extraction from C-parameter
   • Problem: different event shapes yield αs differing by 3σ
3. Recent breakthroughs:
   • Caola et al. 22,27: first calculation of v-dependent shift in three-jet region
   • Nason-Zanderighi 23: applied to αs fitting, discovered severe peak region problems
   • Dasgupta-Hounat 24: studied anomalous dimensions of averages
   • Farren-Colloty et al. 25: numerical verification of anomalous dimension with PanScales shower

Advantages Relative to This Work

1. Unified framework: first consistent theory across entire phase space
2. Analytical calculations: explicit analytical expressions for soft wide-angle contributions
3. Anomalous dimension: derived from first principles rather than numerical fitting
4. Generalizability: method applicable to arbitrary event shape variables

Conclusions and Discussion

Main Conclusions

1. Soft wide-angle gluons cannot be neglected: at O(αs) precision, soft wide-angle gluons contribute constant term ∝CA, equally important as hard collinear term ∝CF
2. Anomalous dimension exists: non-perturbative shift possesses evolution equation with anomalous dimension γ=4(ln2-1)
3. Unified formula: equations (5.23) and (6.1) provide expressions valid in both two-jet and three-jet regions
4. Phenomenological impact: in standard fitting range 0.05<τ<0.3, resummation effects reduce shift by 10-15%, potentially improving αs extraction consistency

Limitations

1. Peak region failure: for τ<0.03, αs(τQ) too large causes perturbative expansion breakdown
2. Leading color approximation: anomalous dimension exponentiation assumption (equations 5.21-5.22) based on numerical observation in large Nc limit
3. NLL precision limit: only includes NLL resummation, NNLL corrections may be important
4. Recoil effects not considered: for certain event shapes (e.g., broadening), momentum conservation scheme choice may affect results
5. Matching uncertainty: parameter p choice (1/2 to 4/3) introduces approximately 5% uncertainty

Future Directions

1. Generalization to other observables: C-parameter, jet broadenings, etc.
2. NNLL precision: include higher-order logarithms and hard functions
3. Phenomenological applications: refit LEP data to extract αs
4. Recoil effects: study impact of different momentum conservation schemes
5. Full color structure: complete calculations beyond leading color

In-Depth Evaluation

Strengths

1. Theoretical Depth

• Discovery of new physics: CA contribution from soft wide-angle gluons is completely novel effect
• Mathematical rigor: calculations from Sudakov decomposition to elliptic integrals complete and reliable
• Clear concepts: kt stratification clearly distinguishes hard collinear, soft collinear, and soft wide-angle contributions

2. Methodological Innovation

• Anomalous dimension derivation: emerges naturally from phase space integration rather than external assumption
• Unified framework: equation (5.8) elegantly interpolates between two limits
• Matching scheme: modified logarithm technique for handling kinematic endpoints practical

3. Computational Technique

• Analytical-numerical combination: reliable numerical evaluation of complex elliptic integrals
• Rapidity function fT(η): clever construction enabling integration over entire real axis
• Color structure separation: explicit CF and CA separation aids physical understanding

4. Writing Quality

• Introduction clearly motivates problem, equations (1.1-1.3) progressively developed for accessibility
• Abundant figures (7 main plots), linear and logarithmic coordinate comparisons effective
• Complete appendices provide technical details

Weaknesses

1. Theoretical Assumptions

• Anomalous dimension exponentiation (equation 5.21): based only on numerical observation from ref.25, lacks analytical proof
• Leading color approximation: CA/CF term separation may have corrections in full color theory
• Gluer model itself: single-gluon approximation validity for ΛQCD/Q≪1 not rigorously proven

2. Numerical Precision

• Matching parameter p: 1/2,4/3 range selection lacks systematic justification
• Running coupling treatment: theory fails for τ<0.03 when αs(τQ) diverges
• Infrared cutoff µI: choice in equation (5.22) somewhat arbitrary, insufficient discussion

3. Phenomenological Application

• No actual data fitting: only provides theoretical predictions, does not verify whether αs extraction consistency improves
• Monte Carlo comparison missing: no comparison with PYTHIA/HERWIG event generators
• Incomplete uncertainty estimation: only considers matching uncertainty, excludes scale variation, etc.

4. Insufficient Generalization Discussion

• Only thrust calculated in detail, γ values for other observables (C-parameter, etc.) not provided
• Recoil-sensitive observable treatment only briefly mentioned in conclusions
• Hadron collider (pp→jets) generalization not discussed

Impact Assessment

Academic Impact

• Fundamental theory contribution: anomalous dimension discovery may open new research directions in non-perturbative QCD
• Methodological value: unified framework serves as paradigm for other power correction studies (e.g., 1/Q²)
• Citation potential: expected to become standard reference for event shape hadronisation corrections

Practical Value

• αs extraction improvement: if phenomenological application succeeds, may resolve long-standing 3σ discrepancy
• Event generator tuning: provides new theoretical constraints for PYTHIA, etc.
• Precision QCD tests: theoretical foundation for high-precision measurements at future colliders (FCC-ee)

Reproducibility

• Complete analytical formulas: Appendix A provides explicit expressions for all functions
• Clear numerical methods: elliptic integral and rapidity integration calculation methods transparent
• High code implementability: formulas (5.23) and (6.1) easily programmed
• Missing: no public code or numerical tables provided

Applicable Scenarios

Direct Application

1. LEP data reanalysis: reextract αs(MZ), fitting range 0.05<τ<0.3
2. Theoretical predictions: provide NLL+NP predictions for future e+e− colliders like FCC-ee
3. Event shape studies: C-parameter, jet masses and similar observables

Requiring Extension

1. NNLL precision: combine with existing NNLL resummation codes (e.g., ARES)
2. Hadron colliders: pp→dijets requires initial-state radiation and PDF considerations
3. Heavy flavor production: b-jets may have additional mass effects

Not Applicable

1. τ<0.03 peak region: perturbative failure, requires fully non-perturbative methods
2. Multi-jet observables (N≥4 jets): more complex color topology
3. Non-global observables: e.g., jet substructure where non-linear effects important

Selected References

Foundational Work

6 Y. L. Dokshitzer, B. R. Webber, Calculation of power corrections to hadronic event shapes, Phys. Lett. B 352 (1995) 451 9 Y. L. Dokshitzer et al., Dispersive approach to power behaved contributions, Nucl. Phys. B 469 (1996) 93

Most Recent Related Work

22 F. Caola et al., Linear power corrections to e+e− shape variables in the three-jet region, arXiv:2204.02247 23 P. Nason, G. Zanderighi, Fits of αs using power corrections in the three-jet region, arXiv:2301.03607 24 M. Dasgupta, F. Hounat, Exploring soft anomalous dimensions, JHEP 09 (2025) 060 25 C. Farren-Colloty et al., Anomalous scaling of linear power corrections, arXiv:2507.18696

Methodological References

14 A. Banfi et al., Interplay between perturbative and non-perturbative effects with the ARES method, JHEP 08 (2023) 221 29 A. Banfi et al., A general method for the resummation of event-shape distributions, JHEP 05 (2015) 102

Overall Assessment

This is a high-quality theoretical physics paper achieving substantial progress on the fundamental QCD non-perturbative corrections problem. The discovery of anomalous dimension and establishment of unified framework possess significant academic value. Computational techniques are solid and physical picture clear. Main weaknesses are lack of actual data fitting verification and rigorous proof of certain theoretical assumptions (e.g., anomalous dimension exponentiation). This work provides important theoretical tools for future precise αs determination and hadronisation understanding, with expected significant academic impact. Recommended follow-up work: (1) comparison with LEP data to verify; (2) extension to NNLL precision; (3) calculation of anomalous dimensions for other event shapes.