Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}_2$ from lattice QCD at physical quark masses

Cet article présente une première calcul en QCD sur réseau à masses de quarks physiques, dans le cadre de la théorie effective à grand moment (LaMET), de la fonction de distribution de parton (PDF) d'hélicité isovecteur du proton à l'ordre NNLO et de son moment twist-3 $\tilde{d}_2$, ce dernier étant proportionnel à la force de Lorentz de couleur moyenne subie par les quarks.

L'essentiel

🌌 Le Grand Mystère du Spin du Proton

Imaginez que le proton (la particule au cœur de chaque atome) est comme un tourbillon d'énergie. Il tourne sur lui-même, une propriété appelée "spin". Pendant longtemps, les physiciens pensaient que ce tourbillon venait uniquement de la rotation de ses trois pièces principales (les quarks), un peu comme si trois patineurs sur glace tenaient la main et tournaient ensemble.

Mais il y a 30 ans, une expérience a révélé un choc : les quarks ne contribuaient qu'à environ 30 % du spin total. Où est passée la moitié du tourbillon ? C'est ce qu'on appelle le "mystère du spin du proton".

🔍 La Nouvelle Enquête : Regarder à l'Intérieur

Cette nouvelle étude, menée par des chercheurs du Laboratoire national de Brookhaven et d'autres institutions, utilise une méthode très puissante appelée QCD sur réseau (une sorte de super-calculateur qui simule l'univers à l'échelle des particules).

Leur objectif ? Prendre une "photo" ultra-précise de la façon dont les quarks tournent à l'intérieur du proton.

1. L'Analogie du Train et des Passagers

Pour voir comment les quarks bougent, les chercheurs utilisent une technique appelée LaMET (Théorie effective à grand impulsion).

• L'idée : Imaginez que vous êtes dans un train qui roule très vite. Si vous regardez par la fenêtre, les passagers (les quarks) semblent figés ou se déplacent lentement par rapport à vous.
• Le problème : En physique quantique, on ne peut pas prendre de "photos" directes de l'intérieur d'un proton au repos car les règles du jeu changent (c'est trop flou).
• La solution : Les chercheurs accélèrent virtuellement le proton sur leur super-calculateur. En le regardant aller très vite, ils peuvent "décoder" la structure interne, un peu comme si on utilisait un stroboscope pour figer le mouvement d'une hélice de ventilateur.

2. La Carte des Quarks (La PDF d'hélicité)

Le résultat principal de l'article est une carte détaillée (appelée PDF d'hélicité) qui montre où se trouvent les quarks et comment ils tournent.

• L'analogie : C'est comme si on dessinait une carte de la densité de population d'une ville, mais au lieu de montrer où les gens vivent, on montre dans quelle direction ils regardent (vers le haut ou vers le bas) en fonction de leur vitesse.
• La découverte : Les chercheurs ont trouvé que dans la zone centrale (là où les quarks ont une vitesse moyenne), il y a plus de quarks qui tournent dans le sens du proton que ce que les anciennes estimations suggéraient. C'est comme si on découvrait qu'il y avait plus de cyclistes dans le sens de la circulation qu'on ne le pensait.

3. La Force Invisible (Le Moment ˜d2)

Au-delà de la simple rotation, l'étude s'intéresse à une interaction plus subtile : la force de Lorentz de couleur.

• L'analogie : Imaginez un quark comme un petit bateau dans une rivière turbulente. Le proton est la rivière. Les quarks ne tournent pas seulement sur eux-mêmes ; ils sont aussi poussés et tirés par les autres particules (les gluons) comme par un courant violent.
• Le résultat : Les chercheurs ont mesuré la "force moyenne" que ce courant exerce sur le quark. Ils ont découvert que cette force est presque nulle (très proche de zéro).
• Pourquoi c'est important ? Cela signifie que, dans cette configuration spécifique, les quarks ne sont pas trop secoués par les turbulences internes. C'est une surprise qui aide à comprendre pourquoi le proton est si stable.

🛠️ Comment ont-ils fait ? (Les Outils Magiques)

Pour obtenir ces résultats, ils ont dû surmonter deux obstacles majeurs :

1. Le Bruit de Fond (La Renormalisation) :
   Les calculs sur ordinateur génèrent beaucoup de "bruit" mathématique (des erreurs infinies). Les chercheurs ont utilisé une technique appelée renormalisation hybride.

   • L'analogie : C'est comme essayer d'écouter une conversation dans un stade bruyant. Ils ont d'abord éliminé le bruit constant (le fond sonore) pour isoler la voix claire de la conversation (le signal physique réel).

2. Le Flou aux Bords (Les Extrémités) :
   Leur méthode fonctionne très bien au centre de la carte (les vitesses moyennes), mais devient floue aux extrémités (très lent ou très rapide).

   • L'analogie : Imaginez une photo prise avec un objectif qui est très net au centre mais flou sur les bords. Pour combler les trous, ils ont utilisé des modèles mathématiques intelligents (basés sur des règles physiques connues) pour "dessiner" ce qui se passe sur les bords flous, en s'assurant que cela colle avec la partie nette du centre.

🏆 Pourquoi est-ce important ?

Cette étude est une victoire pour la précision.

• C'est la première fois qu'on calcule directement cette force interne (le moment ˜d2) avec une telle précision à partir des lois fondamentales de la physique, sans se fier uniquement à des modèles approximatifs.
• Cela nous rapproche de la résolution finale du mystère du spin du proton. On commence à comprendre comment la "partie manquante" du spin est partagée entre la rotation des quarks, celle des gluons (les "colles" qui tiennent le tout ensemble) et le mouvement orbital.

En résumé : Ces chercheurs ont utilisé un super-calculateur pour simuler un proton ultra-rapide, ont nettoyé le bruit mathématique pour voir clairement comment les quarks tournent, et ont découvert que les forces internes qui les secouent sont étonnamment faibles. C'est une pièce cruciale du puzzle pour comprendre de quoi est fait notre univers.

Résumé technique

1. Problématique et Contexte

La structure de spin du proton reste l'un des défis majeurs de la chromodynamique quantique (QCD). Bien que les mesures expérimentales suggèrent que le spin des quarks ne contribue qu'à environ 30 % au spin total du proton, la répartition précise de ce spin entre les quarks, les gluons et le moment angulaire orbital nécessite une détermination ab initio précise des fonctions de distribution de partons (PDF) d'hélicité, notées $g_1(x)$.

De plus, au-delà de l'approximation de twist-2 (dominante), la fonction de structure de spin transversal $g_T(x) = g_1(x) + g_2(x)$ contient des informations cruciales sur les corrélations quark-gluon. Un moment clé de cette structure, le moment réduit de twist-3 noté $\tilde{d}_2$, est proportionnel à la force de Lorentz de couleur moyenne subie par un quark dans le proton. Sa détermination expérimentale est difficile et imprécise.

Le défi théorique réside dans le fait que les PDF sont définies par des opérateurs non locaux sur le cône de lumière, ce qui les rend inaccessibles au calcul direct en QCD sur réseau (qui opère dans l'espace euclidien).

2. Méthodologie

Les auteurs utilisent la théorie effective à grand moment (LaMET - Large Momentum Effective Theory) combinée à des calculs de QCD sur réseau à masses de quarks physiques.

• Configuration de réseau :

  • Calcul effectué sur un réseau de taille $64^3 \times 64$ avec un espacement $a = 0.076$ fm.
  • Action des quarks : HISQ (Highly Improved Staggered Quark) pour les quarks de la mer ($N_f=2+1$) et action Clover-Wilson pour les quarks de valence.
  • Masses de quarks : Ajustées aux masses physiques (pion à $\approx 140$ MeV).
  • Impulsions du proton : États boostés à $P_z = \{0, 0.25, 1.02, 1.53\}$ GeV.

• Extraction des éléments de matrice :

  • Calcul des corrélateurs à trois points pour les opérateurs non locaux d'hélicité (canal twist-2) et transversaux (canal twist-3).
  • Utilisation de la méthode AMA (All Mode Averaging) et du lissage de moment en jauge de Coulomb pour améliorer le rapport signal/bruit.
  • Extraction des éléments de matrice nus $\tilde{h}^B$ via des ajustements à deux états des rapports de corrélateurs.

• Renormalisation et Matching :

  • Schéma hybride : Utilisation d'un schéma de renormalisation hybride pour traiter les divergences linéaires et logarithmiques des lignes de Wilson, incluant une resommation des renormalons (LRR) pour régulariser l'ambiguïté infrarouge.
  • Matching perturbatif : Conversion des quasi-PDF (définies sur le réseau) vers les PDF de cône de lumière dans le schéma $\overline{\text{MS}}$. Ce matching est effectué à l'ordre NNLO (Next-to-Next-to-Leading Order) avec des resommations de groupe de renormalisation (RG) et de seuil (threshold resummation) pour contrôler les logarithmes grands.
  • Extrapolation asymptotique : Les éléments de matrice à grande distance sont extrapolés avant la transformation de Fourier pour obtenir la dépendance en $x$.

• Approche complémentaire (SDF) :

  • Pour contraindre les régions de bout de spectre ($x \to 0$ et $x \to 1$) où l'expansion LaMET échoue, les auteurs utilisent l'analyse par factorisation à courte distance (SDF) en ajustant des modèles de comportement asymptotique aux éléments de matrice dans l'espace des coordonnées.

3. Contributions Clés

1. Calcul NNLO de la PDF d'hélicité isovectorielle : Première détermination complète de la PDF d'hélicité $g_1^{u-d}(x)$ à l'ordre NNLO avec des masses de quarks physiques et une extrapolation vers la limite du continu implicite via un maillage fin.
2. Extraction directe du moment $\tilde{d}_2$ : Première extraction du moment de twist-3 $\tilde{d}_2$ dans le schéma $\overline{\text{MS}}$ directement à partir des éléments de matrice non locaux, sans passer par des opérateurs locaux sujets à des mélanges complexes. Cela repose sur un développement OPE (Opérateur Product Expansion) amélioré qui isole la contribution véritable de twist-3.
3. Traitement des renormalons : Application systématique de la resommation des renormalons dans le schéma hybride pour éliminer les ambiguïtés de puissance dans le matching perturbatif.

4. Résultats Principaux

• Moment $\tilde{d}_2$ :
  Le moment réduit de twist-3 isovectoriel est extrait à l'échelle $\mu = 2$ GeV :
  $$\tilde{d}_2^{u-d}(2 \text{ GeV}) = 0.0024(46)$$
  Ce résultat est statistiquement compatible avec zéro. Cela indique une forte suppression de la corrélation quark-gluon véritable (au-delà de l'approximation de Wandzura-Wilczek) dans ce moment, en accord avec les contraintes phénoménologiques et les études de modèles récents.

• PDF d'hélicité $g_1(x)$ :

  • Les moments de la PDF (extraction via OPE) sont : $\langle x \rangle^{u-d} = 0.291(7)$ et $\langle x^2 \rangle^{u-d} = 0.101(10)$ à 2 GeV.
  • La forme de la PDF dans la région intermédiaire $x \in [0.25, 0.75]$ montre une distribution plus élevée que les ajustements globaux actuels (comme JAM22 ou BDSSV24), en particulier pour $x \in [0.4, 0.7]$.
  • Les incertitudes systématiques sont contrôlées par la variation d'échelle et la comparaison de deux formes d'extrapolation asymptotique.

• Comparaison avec d'autres méthodes :
  Une comparaison avec des calculs récents utilisant la jauge de Coulomb (sans ligne de Wilson) montre une bonne concordance pour $x < 0.6$, mais des écarts pour $x > 0.6$, potentiellement dus aux corrections de puissance à l'impulsion modérée utilisée ($P_z \approx 1.5$ GeV).

5. Signification et Perspectives

Ce travail représente une avancée majeure dans la compréhension ab initio de la structure de spin du proton :

• Il valide la faisabilité de l'extraction de moments de twist-3 (comme $\tilde{d}_2$) directement dans le schéma $\overline{\text{MS}}$ via des corrélateurs non locaux, offrant une alternative robuste aux méthodes basées sur des opérateurs locaux.
• La détermination de la PDF d'hélicité à l'ordre NNLO avec des masses physiques réduit les incertitudes théoriques et fournit une contrainte rigoureuse pour les analyses globales expérimentales.
• Le résultat $\tilde{d}_2 \approx 0$ suggère que les effets de corrélation quark-gluon dans ce moment spécifique sont négligeables par rapport à la structure de twist-2, ce qui simplifie l'interprétation physique de la force de Lorentz de couleur dans le proton.

Les auteurs soulignent que pour étendre la fenêtre de validité de LaMET vers les régions de $x$ extrêmes et réduire les corrections de puissance, de futurs calculs nécessiteront des impulsions de hadrons plus élevées et un contrôle accru des effets systématiques du réseau (plusieurs espacements de réseau, statistiques accrues).