Corrections QCD à l'ordre suivant-suivant-dominant pour la diffusion profonde semi-inclusive polarisée

Informations fondamentales

• ID de l'article: 2404.09959
• Titre: Next-to-Next-to-Leading Order QCD Corrections to Polarized Semi-Inclusive Deep-Inelastic Scattering
• Auteurs: Saurav Goyal, Roman N. Lee, Sven-Olaf Moch, Vaibhav Pathak, Narayan Rana, V. Ravindran
• Classification: hep-ph (physique des hautes énergies - phénoménologie), hep-ex (physique des hautes énergies - expérience)
• Date de publication: 15 avril 2024 (soumission arXiv), 31 décembre 2024 (dernière révision)
• Lien de l'article: https://arxiv.org/abs/2404.09959

Résumé

La diffusion profonde semi-inclusive polarisée (SIDIS) est un processus clé pour résoudre l'énigme du spin du proton. Cet article présente pour la première fois les résultats complets à l'ordre suivant-suivant-dominant (NNLO) de la chromodynamique quantique perturbative pour le processus SIDIS polarisé. L'analyse inclut tous les canaux de partons pour la diffusion lepton-hadron polarisée, ainsi que les hadrons de spin moyen identifiés dans l'état final. L'analyse numérique démontre l'importance des corrections NNLO et révèle une réduction de la dépendance d'échelle résiduelle dans la gamme cinématique explorée par le futur collisionneur électron-ion (EIC).

Contexte et motivation de la recherche

L'énigme du spin du proton

L'énigme du spin du proton est un problème de longue date non résolu en physique des particules. Traditionnellement, on pensait que le spin du proton provenait principalement du spin de ses quarks constitutifs, mais les expériences ont révélé que la contribution du spin des quarks au spin total du proton n'est que d'environ 30%, bien inférieure aux 100% attendus. Cela a incité les scientifiques à rechercher d'autres sources de contribution, notamment le spin des gluons et le moment angulaire orbital.

Importance de la SIDIS

1. Mesure des fonctions de structure: La diffusion profonde inélastique (DIS) fournit des informations essentielles sur la structure interne des hadrons, encodées par les fonctions de structure via les distributions de partons
2. Information de polarisation: La SIDIS polarisée permet de séparer et d'extraire les fonctions de distribution d'hélicité pour différentes saveurs de quarks
3. Factorisation QCD: Le processus SIDIS permet de séparer la dynamique perturbative à courte distance de la dynamique non-perturbative à longue distance

Limitations théoriques existantes

• Les prédictions théoriques précédentes pour la SIDIS polarisée n'atteignaient que la précision du prochain ordre dominant (NLO)
• Avec la construction de nouvelles installations expérimentales comme l'EIC, des prédictions théoriques de plus haute précision sont nécessaires
• Le développement des techniques de calcul NNLO rend possible le calcul complet à deux boucles

Contributions principales

1. Première réalisation du calcul QCD NNLO pour la SIDIS polarisée, comblant une lacune théorique importante dans ce domaine
2. Fourniture des coefficients de fonction complets pour tous les canaux de partons, incluant la diffusion quark-quark, quark-gluon, gluon-quark et gluon-gluon
3. Développement de nouvelles techniques pour traiter les intégrales complexes d'espace de phase, en particulier les intégrales multi-boucles avec contraintes supplémentaires
4. Vérification numérique montrant que les corrections NNLO réduisent significativement la dépendance d'échelle, améliorant la fiabilité des prédictions théoriques
5. Fourniture de prédictions théoriques de haute précision pour l'expérience EIC, soutenant les futures recherches sur la structure du spin du proton

Détails méthodologiques

Définition de la tâche

Calcul des corrections QCD NNLO pour le processus de diffusion profonde semi-inclusive polarisée: $l(k_l) + H(P) \rightarrow l(k'_l) + H'(P_H) + X$

où l'entrée est un lepton polarisé et un hadron polarisé, et la sortie est un lepton diffusé, un hadron marqué et d'autres particules.

Cadre théorique

Factorisation QCD

La fonction de structure polarisée $g_1$ peut être factorisée comme: $g_1 = \sum_{a,b} \int_x^1 \frac{dx_1}{x_1} \Delta f_a(x_1, \mu_F^2) \int_z^1 \frac{dz_1}{z_1} D_b(z_1, \mu_F^2) \times G_{1,ab}\left(\frac{x}{x_1}, \frac{z}{z_1}, \mu_F^2, Q^2\right)$

où:

• $\Delta f_a$: fonction de distribution de partons polarisés
• $D_b$: fonction de fragmentation de partons
• $G_{1,ab}$: fonction de coefficient (objet central du calcul de cet article)

Développement perturbatif

La fonction de coefficient se développe selon la constante de couplage fort: $G_{1,ab}(\mu_F^2) = \sum_{i=0}^{\infty} a_s^i(\mu_R^2) G_{1,ab}^{(i)}(\mu_F^2, \mu_R^2)$

Techniques de calcul

Schéma de régularisation

• Utilisation de la régularisation dimensionnelle pour traiter les divergences ultraviolettes et infrarouge
• Adoption du schéma de Larin pour traiter la définition de $\gamma_5$ dans l'espace-temps à $D$ dimensions
• Conversion vers le schéma $\overline{MS}$ via des constantes de renormalisation finies

Techniques d'intégration d'espace de phase

1. Méthode d'unitarité: conversion des intégrales d'espace de phase en intégrales de boucle
2. Identités IBP: réduction des intégrales complexes aux intégrales principales
3. Méthode des équations différentielles: résolution des équations différentielles systématiques des intégrales principales
4. Conditions aux limites: calcul des conditions aux limites dans la limite de seuil

Calcul des intégrales principales

• Processus RV (interférence réel-virtuel): 7 intégrales principales
• Processus RR (émission purement réelle): 20 intégrales principales
• Résultats exprimés en termes de polylogarithmes multiples et de polylogarithmes de Nielsen

Points d'innovation technique

Traitement des contraintes

Le processus SIDIS contient la contrainte supplémentaire $z/z_1 = p_a \cdot p_b/(p_a \cdot q)$, traitée comme suit: $\delta\left(z' - \frac{p_a \cdot p_b}{p_a \cdot q}\right) \rightarrow -\frac{1}{\pi} \text{Im}\left(\frac{1}{z' - \frac{p_a \cdot p_b}{p_a \cdot q} + i\epsilon}\right)$

Prolongement analytique des fonctions multivaluées

Les fonctions hypergéométriques rencontrées nécessitent un prolongement analytique dans différentes régions cinématiques: $_2F_1\left(1, \frac{\epsilon}{2}; 1+\frac{\epsilon}{2}; \frac{x'(1-x')}{z'(1-z')}\right)$

La connexion lisse entre différentes régions est réalisée via les transformations de Pfaff et la prescription $i\epsilon$ de Feynman.

Configuration expérimentale

Paramètres de calcul numérique

• Fonctions de distribution de partons:
  • LO/NLO: MAPPDF10NLO
  • NNLO: MAPPDF10NNLO
• Fonctions de fragmentation:
  • LO/NLO: MAPFF10NLO
  • NNLO: MAPFF10NNLO
• Nombre de saveurs de quarks actifs: $n_f = 3$
• Choix d'échelle: $\mu_R^2 = \mu_F^2 = Q^2$

Gamme cinématique

• Énergie COMPASS: $\sqrt{s} = 17,4$ GeV
• Énergie EIC: $\sqrt{s} = 45$ GeV
• Gamme d'intégration en z: $0,2 < z < 0,85$
• Contrainte en y: $0,5 \leq y \leq 0,9$

Étude de variation d'échelle

Adoption d'un schéma de variation d'échelle à 7 points: $(k_1, k_2) \in [1/2, 2]$ et $1/2 \leq k_1/k_2 \leq 2$

Résultats expérimentaux

Résultats principaux

Contributions des différents canaux de partons

Dans le processus de production de $\pi^+$, le canal de diffusion quark-quark domine, avec une contribution bien supérieure aux canaux quark-gluon et gluon-quark.

Impact numérique des corrections NNLO

1. Réduction significative de la dépendance d'échelle: la bande d'incertitude d'échelle des prédictions NNLO est nettement plus petite que celle du NLO
2. Corrections plus importantes dans la région de bas $Q^2$: à $Q^2 \sim 30$ GeV$^2$, il existe une différence notable entre les prédictions NNLO et NLO
3. Convergence dans la région de haut $Q^2$: avec l'augmentation de $Q^2$, les prédictions de différents ordres convergent progressivement

Comparaison avec les données expérimentales

La figure 3 montre la comparaison du rapport $g_1^{\pi^+}/F_1^{\pi^+}$ avec les données expérimentales de COMPASS:

• Les prédictions NNLO s'accordent bien avec les données expérimentales
• La réduction des incertitudes théoriques rend la comparaison avec l'expérience plus précise
• Fournit une référence théorique fiable pour les futures expériences EIC

Analyse de la dépendance d'échelle

L'étude de la variation d'échelle à différentes valeurs de $Q^2$ montre:

• $Q^2 = 30$ GeV$^2$: dépendance d'échelle la plus forte, corrections NNLO les plus significatives
• $Q^2 = 800$ GeV$^2$: dépendance d'échelle la plus faible, prédictions de différents ordres proches
• La tendance globale est conforme aux attentes de la QCD perturbative

Travaux connexes

Développement historique

1. Calculs NLO: de Florian, Stratmann, Vogelsang (1998) ont d'abord complété le calcul NLO pour la SIDIS polarisée
2. SIDIS non-polarisée: les calculs NNLO ont été récemment complétés, fournissant des références importantes pour ce travail
3. Resommation des logarithmes de seuil: développée jusqu'à la précision N³LL

Développement technique

• Techniques d'intégrales de boucle: maturation des méthodes IBP et des techniques d'équations différentielles
• Intégration d'espace de phase: application de la méthode d'unitarité
• Calculs multi-boucles: développement d'outils de calcul automatisés

Relation avec ce travail

Cet article applique pour la première fois ces techniques avancées au calcul NNLO de la SIDIS polarisée, réalisant une percée importante dans ce domaine.

Conclusions et discussions

Conclusions principales

1. Réalisation du premier calcul NNLO pour la SIDIS polarisée, fournissant des résultats complets pour tous les canaux de partons
2. Les corrections NNLO améliorent significativement la précision des prédictions théoriques, avec une réduction notable de la dépendance d'échelle
3. Fourniture d'un soutien théorique de haute précision pour l'expérience EIC, facilitant la mesure précise de la structure du spin du proton

Signification physique

• Amélioration de la précision des prédictions théoriques des fonctions de structure polarisées
• Fourniture d'un cadre théorique plus fiable pour l'extraction des fonctions de distribution de partons polarisés
• Avancement de la recherche théorique sur l'énigme du spin du proton

Directions futures

1. Corrections électrofaibles: inclusion des corrections d'interaction électrofaible
2. Calculs d'ordre supérieur: développement vers la précision N³LO
3. Analyse phénoménologique: ajustement global combinant les données expérimentales

Évaluation approfondie

Avantages

Innovation technique

1. Complexité de calcul: traitement réussi des intégrales d'espace de phase à deux boucles avec contraintes supplémentaires
2. Contribution méthodologique: développement d'une méthode systématique pour traiter le prolongement analytique des fonctions multivaluées
3. Stabilité numérique: réalisation de calculs numériques stables sur toute la gamme cinématique

Contributions théoriques

1. Complétude: fourniture de résultats NNLO complets pour tous les canaux de partons
2. Précision: vérification des résultats par plusieurs méthodes indépendantes
3. Praticabilité: les résultats ont été transformés en un programme numérique utilisable

Impact phénoménologique

1. Pertinence expérimentale: fourniture d'un soutien théorique important pour les expériences EIC imminentes
2. Amélioration de la précision: réduction significative des incertitudes théoriques
3. Perspectives d'application: établissement d'une base solide pour la détermination précise des PDF polarisés

Insuffisances

Limitations techniques

1. Complexité de calcul: le processus de calcul est extrêmement complexe, rendant la vérification des résultats difficile
2. Évaluation numérique: l'évaluation numérique dans certaines régions cinématiques nécessite encore une optimisation
3. Effets d'ordre supérieur: exclusion des corrections électrofaibles et autres effets

Limitations méthodologiques

1. Dépendance du schéma de régularisation: les résultats dépendent du choix spécifique du schéma de régularisation
2. Choix d'échelle: le choix optimal d'échelle nécessite une recherche supplémentaire
3. Effets non-perturbatifs: non-prise en compte des effets de torsion d'ordre supérieur

Évaluation de l'impact

Impact académique

1. Travail fondateur: comblage d'une lacune théorique importante dans ce domaine
2. Valeur méthodologique: fourniture de références techniques pour des calculs complexes similaires
3. Potentiel de citation: devrait devenir une référence importante dans ce domaine

Valeur pratique

1. Soutien expérimental: service direct aux grands projets expérimentaux comme l'EIC
2. Application phénoménologique: fourniture d'entrées importantes pour l'ajustement global des PDF
3. Établissement de normes: possibilité de devenir une norme de calcul dans ce domaine

Reproductibilité

1. Code ouvert: les auteurs s'engagent à fournir le code de calcul complet
2. Vérification des résultats: vérification par plusieurs méthodes indépendantes
3. Détails techniques: l'article fournit des détails techniques suffisants

Scénarios d'application

1. Expériences de physique des hautes énergies: expériences de diffusion polarisée EIC, COMPASS, etc.
2. Recherche théorique: analyse théorique de la structure du spin du proton
3. Recherche phénoménologique: analyse d'ajustement global des PDF polarisés
4. Recherche méthodologique: développement technique des calculs multi-boucles complexes

Évaluation générale: Cet article est un travail de haute qualité possédant une valeur théorique et une signification pratique importantes. Les auteurs ont réussi à compléter le premier calcul NNLO pour le processus SIDIS polarisé, réalisant non seulement une percée technique majeure, mais fournissant également un outil théorique important pour comprendre la structure du spin du proton. Bien que le calcul soit extrêmement complexe, les auteurs ont vérifié l'exactitude des résultats par plusieurs méthodes et ont démontré les importants effets physiques des corrections NNLO. Ce travail aura un impact profond sur la recherche théorique et expérimentale en physique des hautes énergies.