2025-11-25T13:22:18.679604

Characterising hadronisation across the phase space

Banfi, El-Menoufi
Leading hadronisation corrections to two-jet global event shapes amount to a shift in the corresponding perturbative distributions. It has been recently established that this shift depends significantly on the value of the considered event shape. These analyses consider perturbative configurations with only three partons emitting an ultra-soft non-perturbative gluon. These are dominant in the three-jet region. However, multiple soft and/or collinear emissions need to be considered to accurately describe event shape distributions near their Sudakov peak. In this region, non-perturbative shifts are usually computed by considering only two hard emitters. In this paper, we find that this approximation misses an important correction to the shift due to the emission of an additional soft wide-angle gluon. We then compute its contribution, and embed it in a general treatment for the shift that is valid in both the two-jet and three-jet regions.
academic

Caractérisation de l'hadronisation dans l'espace des phases

Informations de base

  • ID de l'article: 2511.17708
  • Titre: Characterising hadronisation across the phase space
  • Auteurs: Andrea Banfi (University of Sussex), Basem Kamal El-Menoufi (Monash University)
  • Classification: hep-ph (physique des hautes énergies - phénoménologie), hep-th (physique des hautes énergies - théorie)
  • Date de soumission: 21 novembre 2025
  • Lien de l'article: https://arxiv.org/abs/2511.17708

Résumé

Cet article étudie les corrections d'hadronisation au premier ordre pour les variables de forme d'événement globales (event shapes) dans l'annihilation e+e−. La théorie existante suggère que les corrections d'hadronisation se manifestent comme un décalage (shift) de la distribution perturbative, et ce décalage dépend fortement de la valeur de la variable de forme d'événement. Les recherches antérieures considéraient principalement les configurations d'états finaux à trois particules dans la région à trois jets, ou seulement deux émetteurs durs dans la région à deux jets. Cet article découvre que dans la région à deux jets proche du pic de Sudakov, la méthode traditionnelle omet des corrections importantes dues à l'émission de gluons mous à grand angle. Les auteurs calculent cette contribution et établissent un cadre théorique unifié valable dans les régions à deux et trois jets.

Contexte de recherche et motivation

Problèmes fondamentaux à résoudre

  1. Dépendance en espace des phases des corrections non-perturbatives: Les corrections d'hadronisation ne sont pas constantes, mais dépendent fortement de la valeur de la variable de forme d'événement v
  2. Description unifiée des régions à deux et trois jets: La théorie existante s'applique soit strictement à la limite à deux jets (τ→0), soit uniquement à la région à trois jets, sans cadre unifié
  3. Contribution manquante des gluons mous à grand angle: Le calcul traditionnel de la limite à deux jets ne considère que les émissions molles colinéaires, négligeant les effets importants des gluons mous à grand angle

Importance du problème

  • Détermination précise de la constante de couplage fort αs: L'analyse des données de LEP montre que les valeurs de αs extraites de différentes variables de forme d'événement s'écartent de la moyenne mondiale de jusqu'à 3 écarts-types, rendant crucial une compréhension précise des corrections d'hadronisation pour éliminer ces divergences
  • Cohérence théorique: Nécessité de fournir des prédictions théoriques cohérentes dans toute l'étendue de l'espace des phases, en particulier près du pic de la distribution de forme d'événement (τ≃0,025 à l'énergie LEP1)
  • Compréhension de la dynamique non-perturbative de la QCD: L'hadronisation est l'un des phénomènes QCD les plus fondamentaux mais les plus difficiles à comprendre

Limitations des méthodes existantes

  1. Limite strictement à deux jets: La méthode traditionnelle basée sur le modèle dispersif ne considère que le dipôle qq̄ plus un nombre arbitraire de gluons mous colinéaires, prédisant un décalage non-perturbatif ζV(0) constant
  2. Calculs à trois jets à ordre fixe: Les travaux récents 22,23 ont calculé à la précision O(αs) le décalage ζV(v) pour les états finaux à trois particules, découvrant sa forte dépendance en v, mais cela s'applique uniquement à la région à trois jets
  3. Terme constant O(αs) manquant: Le calcul de la limite à deux jets ne fournit que le terme αs ln(1/τ), impossible de déterminer le terme constant O(αs), rendant impossible de contrôler les corrections où ζV(τ)/ζV(0) s'écarte de 1

Contributions principales

  1. Découverte et calcul de la contribution des gluons mous à grand angle: Démonstration qu'une émission de gluon mou à grand angle supplémentaire produit une nouvelle contribution proportionnelle à CA, dont le terme constant O(αs) est 16(ln2-1)CA
  2. Établissement de la dimension anormale γ: Découverte que le décalage non-perturbatif possède une évolution de dimension anormale, γ=4(ln2-1), avec son expression de resommation
  3. Construction d'un cadre unifié: Proposition d'une formule générale valable dans les régions à deux et trois jets (équations 5.8 et 5.23), réalisant une interpolation lisse
  4. Clarification de la structure de couleur des corrections O(αs):
    • Contribution colinéaire dure: ∝CF, coefficient -3
    • Contribution mole à grand angle: ∝CA, coefficient 16(ln2-1)
  5. Fourniture d'un schéma d'appariement: Schéma spécifique pour l'appariement des résultats de resommation avec les calculs à ordre fixe

Détails de la méthode

Définition de la tâche

Étude des corrections non-perturbatives de la variable de forme d'événement V dans le processus e+e−→hadrons. Définition:

  • Variable de forme d'événement: V({pi}), fonction de tous les moments des particules finales
  • Variable thrust: T = (1/Q)max_{n̂T}Σi|p⃗i·n̂T|, τ=1-T
  • Décalage non-perturbatif: ⟨δVNP⟩, tel que la distribution cumulative satisfasse δΣNP(v) = -⟨δVNP⟩∫dZ{ki}δ(v-V({p̃},{ki}))

Cadre théorique

1. Fondements du modèle Gluer

Les effets non-perturbatifs sont modélisés par l'émission d'un gluon "gluer" avec impulsion transverse bien inférieure à l'échelle perturbative:

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

où les facteurs de couleur: Cqq̄ = 1-CA/(2CF), Cqg = Cgq̄ = CA/(2CF)

2. Fonction centrale hV

Définition de la quantité clé du changement de forme d'événement:

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

Pour la variable thrust, la forme explicite est:

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. Décomposition de la limite à deux jets

L'innovation clé est la décomposition de la contribution O(αs):

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

où:

  • c1: provient de la région mole colinéaire
  • c0: contribution colinéaire dure (∝CF) + contribution mole à grand angle (∝CA)
  • c(τ): terme supprimé en puissance, disparaît quand τ→0

4. Paramétrisation de la région mole à grand angle

Utilisation de la paramétrisation impulsion transverse-rapidité:

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

Définition de la fonction de rapidité:

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

son intégrale donne la dimension anormale:

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

Points techniques innovants

1. Identification de la contribution mole à grand angle

Découverte que dans la limite mole kt→0, les fonctions individuelles A, B, C peuvent diverger, mais la combinaison ⟨hij/n̂T_τ⟩ reste finie:

A^(qg)_q = (1/2)√(Q/kt)e^{η/2}  (diverge)
mais |AB| + C² est fini

2. Découverte de la dimension anormale

Par analyse des contributions à différentes échelles de moment kt, découverte de deux régions:

  • kt ~ τQ: donne une constante O(αs) finie
  • kt ≪ τQ: produit un renforcement logarithmique ln(τQ/κus), conduisant à la dimension anormale

Forme finale:

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

3. Combinaison avec la resommation mole colinéaire

Intégration du gluon mole à grand angle dans l'ensemble mole colinéaire:

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. Fonction d'appariement χT(R')

L'effet du gluon mole à grand angle dans l'ensemble mole colinéaire est codé par la fonction:

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

où R' = -τdR/dτ ≃ (2αs/π)CF ln(1/τ) est le paramètre de resommation.

Configuration expérimentale

Configuration du calcul théorique

Cet article est un travail purement théorique ne nécessitant pas d'ajustement de données expérimentales, mais la configuration comprend:

  1. Échelle d'énergie: Q = MZ = 91,2 GeV (énergie LEP1)
  2. Constante de couplage fort: αs(MZ) = 0,118
  3. Précision de resommation: NLL (logarithme suivant l'ordre principal)
  4. Schéma d'appariement: appariement additif, incluant la constante colinéaire dure O(αs)

Méthode de calcul numérique

  1. Intégration de l'espace des phases: paramétrisation par variables de Dalitz (x1,x2), xi = 2Ei/Q
  2. Intégration de rapidité: calcul numérique des intégrales elliptiques complètes K et E
  3. Commutateur de resommation: utilisation du logarithme modifié pour éviter les divergences aux points cinématiques
    ln(1/τ) → (1/p)ln(1 + 1/τ^p - 1/τ^p_max)
    
    où p∈1/2, 4/3 est utilisé pour estimer l'incertitude théorique

Schémas de comparaison

  1. Résultats à trois jets à ordre fixe: calcul O(αs) de ref.23
  2. Limite naïve à deux jets: ζT(0) = 2 (constante)
  3. Résultats de resommation appariés: formule complète de cet article

Résultats expérimentaux

Résultats numériques principaux

1. Contribution du dipôle qq̄

La figure 1 montre ζqq̄(τ)/ζT(0):

  • environ 0,4 à τ=0,3
  • approche rapidement 1 en coordonnées logarithmiques
  • très proche de la limite à deux jets pour τ<10^{-2}

Décomposition:

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

où cqq̄(τ) disparaît rapidement quand τ→0 (figure 3).

2. Contribution du dipôle qg

La figure 2 montre ζqg(τ)/ζT(0):

  • environ 0,33 à τ=0,3
  • converge asymptotiquement vers 1/2 mais très lentement
  • encore significativement éloigné de la valeur limite à τ=10^{-3}

Décomposition:

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

Découverte clé: le terme constant 16(ln2-1)≈2,09 provient de la région mole à grand angle

3. Vérification numérique de la dimension anormale

Intégration de la fonction de rapidité fT(η) (appendice B):

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

Ceci conduit à l'exposant de dimension anormale:

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

Résultats d'appariement complets

La figure 7 présente trois calculs de ζT(τ):

  1. Ordre fixe (bleu): décroît rapidement de 2 à environ 1,6 pour τ<0,1
  2. Resommation appariée (bande rouge):
    • inférieure à l'ordre fixe d'environ 10-15% dans la plage τ∈0,05,0,3
    • augmente rapidement pour τ<0,03 en raison de χT(R')∼-lnR' et de la croissance de αs(τQ)
    • largeur de la bande d'incertitude environ ±5%
  3. Limite naïve à deux jets (ligne pointillée horizontale): ζT(0)=2

Observations clés:

  • Dans la plage d'ajustement typique 0,05<τ<0,3, la resommation pousse le résultat d'ordre fixe vers la limite à deux jets
  • La théorie échoue près du pic τ≃0,025 (αs(τQ) proche du pôle de Landau)

Analyse d'ablation

Importance relative des différentes contributions

À la limite τ→0:

ζT(τ) ∝ (αs(Q)/αs(τQ))^{0.134} × [terme de base + correction mole CA + correction colinéaire dure CF]
  1. Facteur de dimension anormale: (αs(Q)/αs(τQ))^{0.134} environ 1,02 à τ=0,1, effet modéré
  2. Constante mole CA à grand angle: 16(ln2-1)≈2,09, environ 17% du terme logarithmique à τ=0,1
  3. Constante colinéaire dure CF: -3, environ -13% du terme logarithmique

Travaux connexes

Développement historique

  1. Travaux fondateurs des années 1990:
    • Dokshitzer-Webber 6,10: établissement du modèle gluer et de la méthode dispersive
    • Akhoury-Zakharov 7: méthode des renormalons
    • Conclusion principale: les corrections de puissance linéaire ∝ΛQCD/Q sont les plus importantes
  2. Détermination précise des années 2010:
    • Gehrmann et al. 17: extraction de αs à partir de la distribution de thrust
    • Hoang et al. 18: extraction de αs à partir du paramètre C
    • Problème: les différentes formes d'événement donnent des αs différant de 3σ
  3. Percées récentes:
    • Caola et al. 22,27: premier calcul du décalage dépendant de v dans la région à trois jets
    • Nason-Zanderighi 23: utilisation pour l'ajustement de αs, découverte de problèmes graves dans la région du pic
    • Dasgupta-Hounat 24: étude de la dimension anormale de la moyenne
    • Farren-Colloty et al. 25: vérification numérique de la dimension anormale avec le générateur PanScales

Avantages relatifs de cet article

  1. Cadre unifié: première théorie cohérente dans tout l'espace des phases
  2. Calcul analytique: expression analytique explicite de la contribution mole à grand angle
  3. Dimension anormale: dérivée à partir des premiers principes, non ajustement numérique
  4. Généralité: la méthode s'applique à toute variable de forme d'événement

Conclusions et discussion

Conclusions principales

  1. Les gluons moles à grand angle ne sont pas négligeables: à la précision O(αs), les gluons moles à grand angle contribuent un terme constant ∝CA, d'importance égale au terme colinéaire dur ∝CF
  2. Existence de la dimension anormale: le décalage non-perturbatif possède une équation d'évolution, dimension anormale γ=4(ln2-1)
  3. Formule unifiée: les équations (5.23) et (6.1) fournissent une expression valable dans les régions à deux et trois jets
  4. Impact phénoménologique: dans la plage d'ajustement standard 0,05<τ<0,3, l'effet de resommation réduit le décalage de 10-15%, pouvant améliorer la cohérence de l'extraction de αs

Limitations

  1. Échec dans la région du pic: pour τ<0,03, αs(τQ) trop grand rend l'expansion perturbative invalide
  2. Approximation de couleur dominante: l'hypothèse d'exponentiation de la dimension anormale (équations 5.21-5.22) repose sur l'observation numérique dans la limite de grand Nc
  3. Limitation de précision NLL: seule la resommation NLL est incluse, les corrections NNLL peuvent être importantes
  4. Effets de recul non considérés: pour certaines formes d'événement (comme broadening), le choix du schéma de recul peut affecter les résultats
  5. Incertitude d'appariement: le choix du paramètre p (1/2 à 4/3) introduit une incertitude d'environ 5%

Directions futures

  1. Généralisation à d'autres observables: C-parameter, jet broadenings, etc.
  2. Précision NNLL: inclusion des logarithmes d'ordre supérieur et des fonctions dures
  3. Application phénoménologique: réajustement des données LEP pour extraire αs
  4. Effets de recul: étude de l'impact des différents schémas de conservation de l'impulsion
  5. Structure de couleur complète: calcul complet au-delà de la couleur dominante

Évaluation approfondie

Points forts

1. Profondeur théorique

  • Découverte de nouvelle physique: la contribution CA des gluons moles à grand angle est un effet complètement négligé auparavant
  • Rigueur mathématique: calcul complet et fiable de la décomposition de Sudakov aux intégrales elliptiques
  • Clarté conceptuelle: distinction claire par stratification en kt entre les contributions colinéaire dure, mole colinéaire et mole à grand angle

2. Innovation méthodologique

  • Dérivation de la dimension anormale: émerge naturellement de l'intégration de l'espace des phases, plutôt que d'être imposée
  • Cadre unifié: l'équation (5.8) interpole élégamment les deux limites
  • Schéma d'appariement: la technique du logarithme modifié pour traiter les points cinématiques est pratique

3. Technique de calcul

  • Combinaison analytique et numérique: évaluation numérique fiable des intégrales elliptiques complexes
  • Fonction de rapidité fT(η): construction ingénieuse permettant l'intégration sur l'axe réel entier
  • Séparation de structure de couleur: séparation explicite des termes CF et CA facilitant la compréhension physique

4. Qualité de rédaction

  • La section introduction expose clairement la motivation, les équations (1.1-1.3) se déploient progressivement
  • Figures abondantes (7 figures principales), comparaison entre coordonnées linéaires et logarithmiques efficace
  • Appendices complets fournissant les détails techniques

Insuffisances

1. Hypothèses théoriques

  • Exponentiation de la dimension anormale (équation 5.21): basée uniquement sur l'observation numérique de ref.25, manque de preuve analytique
  • Approximation de couleur dominante: la séparation des termes CA/CF peut avoir des corrections en théorie de couleur complète
  • Validité du modèle gluer lui-même: l'approximation de gluon unique n'est pas rigoureusement justifiée quand ΛQCD/Q≪1

2. Précision numérique

  • Paramètre d'appariement p: le choix de la plage 1/2,4/3 manque de justification systématique
  • Traitement du couplage courant: pour τ<0,03, αs(τQ) diverge, la prédiction théorique échoue
  • Coupure infrarouge µI: le choix de µI dans l'équation (5.22) est arbitraire, pas suffisamment discuté

3. Application phénoménologique

  • Pas d'ajustement de données réelles: seules les prédictions théoriques sont fournies, pas de vérification si cela améliore l'extraction de αs
  • Comparaison avec Monte Carlo manquante: pas de comparaison avec les générateurs d'événements PYTHIA/HERWIG
  • Estimation d'incertitude incomplète: seule l'incertitude d'appariement est considérée, pas les variations d'échelle, etc.

4. Discussion insuffisante de la généralisation

  • Seul thrust est calculé en détail, les valeurs de γ pour d'autres observables (C-parameter, etc.) ne sont pas données
  • Le traitement des observables sensibles au recul n'est que brièvement mentionné en conclusion
  • La généralisation aux collisionneurs hadroniques (pp→jets) n'est pas discutée

Évaluation de l'impact

Impact académique

  • Contribution théorique fondamentale: la découverte de la dimension anormale peut ouvrir de nouvelles directions de recherche en QCD non-perturbative
  • Valeur méthodologique: le cadre unifié fournit un exemple pour l'étude d'autres corrections de puissance (comme 1/Q²)
  • Potentiel de citation: devrait devenir une référence standard pour les corrections d'hadronisation des formes d'événement

Valeur pratique

  • Amélioration de l'extraction de αs: si l'application phénoménologique réussit, peut résoudre la divergence de 3σ de longue date
  • Ajustement des générateurs d'événements: fournit de nouvelles contraintes théoriques pour PYTHIA, etc.
  • Test QCD précis: base théorique pour les mesures de haute précision des futurs collisionneurs (FCC-ee)

Reproductibilité

  • Formules analytiques complètes: l'appendice A fournit les expressions explicites de toutes les fonctions
  • Méthode numérique explicite: les méthodes de calcul des intégrales elliptiques et des intégrales de rapidité sont claires
  • Implémentabilité élevée: facile de programmer basé sur les équations (5.23) et (6.1)
  • Manque: pas de code public ou de tableaux numériques fournis

Scénarios d'application

Application directe

  1. Réanalyse des données LEP: réextraction de αs(MZ), plage d'ajustement 0,05<τ<0,3
  2. Prédictions théoriques: fournir des prédictions NLL+NP pour les futurs collisionneurs e+e− comme FCC-ee
  3. Étude des formes d'événement: observables similaires comme C-parameter, jet masses, etc.

Nécessitant une extension

  1. Précision NNLL: combinaison avec les codes de resommation NNLL existants (comme ARES)
  2. Collisionneur hadronique: pp→dijets nécessite de considérer le rayonnement initial et les PDF
  3. Production de saveurs lourdes: les b-jets peuvent avoir des effets de masse supplémentaires

Non applicable

  1. Région du pic τ<0,03: échec perturbatif, nécessite une méthode complètement non-perturbative
  2. Observables multi-jets (N≥4 jets): topologie de couleur plus complexe
  3. Observables non-globales: comme jet substructure, les effets non-linéaires sont importants

Références (sélection)

Travaux fondateurs

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

Travaux connexes récents

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

Références méthodologiques

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


Évaluation globale

Ceci est un article de physique théorique de haut niveau, réalisant des progrès substantiels sur le problème fondamental des corrections non-perturbatives de la QCD. La découverte de la dimension anormale et l'établissement du cadre unifié possèdent une valeur académique importante. La technique de calcul est solide et l'image physique est claire. Les principales insuffisances résident dans l'absence de vérification par ajustement de données réelles et la preuve rigoureuse de certaines hypothèses théoriques (comme l'exponentiation de la dimension anormale). Ce travail fournit des outils théoriques importants pour la détermination précise future de αs et la compréhension de l'hadronisation, et devrait produire un impact académique significatif. Travaux recommandés pour la suite: (1) comparaison avec les données LEP pour vérification; (2) généralisation à la précision NNLL; (3) calcul de la dimension anormale pour d'autres formes d'événement.