Expansion à court temps d'écoulement des opérateurs de twist-deux non-singlets à l'ordre suivant-suivant dominant en QCD

Informations fondamentales

• ID de l'article : 2511.17145
• Titre : Short-flow-time expansion of non-singlet twist-two operators at next-to-next-to-leading order QCD
• Auteurs : Robert V. Harlander, Jonas T. Kohnen, Andrea Shindler (RWTH Aachen University, Lawrence Berkeley National Laboratory, UC Berkeley)
• Classification : hep-ph (physique des hautes énergies - phénoménologie), hep-lat (théorie des champs sur réseau)
• Date de publication : Novembre 2025
• Lien de l'article : https://arxiv.org/abs/2511.17145

Résumé

Le formalisme du flux de gradient (gradient-flow formalism) fournit un cadre théorique permettant de déterminer directement à partir des calculs QCD sur réseau les moments des fonctions de distribution des partons (PDFs). La conversion du schéma de flux de gradient vers le schéma $\overline{\text{MS}}$ nécessite les coefficients d'appariement de l'expansion à court temps d'écoulement (short-flow-time expansion), qui peuvent être calculés en théorie des perturbations. Cet article détermine les coefficients d'appariement pour les six premiers moments non-singlets des PDFs avec une précision atteignant l'ordre suivant-suivant dominant (NNLO) de la constante de couplage fort.

Contexte et motivation de la recherche

Problèmes fondamentaux

1. Difficulté du calcul ab initio des PDFs : Les fonctions de distribution des partons (PDFs) décrivent la probabilité de trouver un parton portant une fraction d'impulsion x donnée dans un hadron. Elles constituent un input non-perturbatif essentiel pour les prédictions théoriques des processus de diffusion hadronique. Traditionnellement, les PDFs sont obtenues par ajustement aux données expérimentales, ce qui manque d'un calcul direct à partir des premiers principes de la QCD.
2. Obstacles techniques rencontrés par la QCD sur réseau :
   • Les PDFs sont définies sur le cône de lumière, tandis que la QCD sur réseau est formulée dans l'espace-temps euclidien, ce qui pose une difficulté fondamentale
   • Les moments de Mellin des PDFs peuvent être exprimés comme des éléments de matrice d'opérateurs de twist deux localisés, qui peuvent en principe être calculés sur réseau
   • Cependant, la discrétisation sur réseau brise la symétrie de rotation de l'espace-temps continu, entraînant un mélange d'opérateurs de haute dimension avec des opérateurs de basse dimension, produisant des divergences en puissance inversement proportionnelles à l'espacement du réseau a
   • Cela rend le calcul pratique des moments d'ordre supérieur (n>4) impossible, et même les moments n=3,4 présentent d'énormes incertitudes statistiques
3. Percée du schéma de flux de gradient :
   • La Ref. 33 propose d'utiliser le formalisme du flux de gradient (GFF) pour éviter complètement les divergences en puissance
   • Le champ d'écoulement au temps d'écoulement t>0 supprime exponentiellement les modes de haute impulsion, rendant les éléments de matrice des opérateurs composites ultraviolet-finis
   • Le GFF préserve la symétrie de rotation complète de l'espace euclidien, éliminant le mélange avec les opérateurs de basse dimension

Motivation de la recherche

Bien que le GFF fournisse des éléments de matrice finis à t>0, physiquement on a besoin des résultats à t=0 dans le schéma $\overline{\text{MS}}$. Cela nécessite d'étendre les opérateurs d'écoulement en une série d'opérateurs localisés via l'expansion à court temps d'écoulement (SFTX), dont les coefficients sont des coefficients d'appariement calculables en perturbation. Les coefficients d'appariement NLO sont connus, mais d'autres applications montrent que les corrections d'ordre supérieur sont cruciales pour réduire les erreurs systématiques et stabiliser l'extrapolation t→0. Par conséquent, cet article pousse le calcul des coefficients d'appariement à la précision NNLO.

Contributions principales

1. Premier calcul des coefficients d'appariement NNLO : Détermination des coefficients d'expansion à court temps d'écoulement ζₙ(t,μ) pour les opérateurs de twist-deux non-singlets avec n=1 à 6, avec une précision de $\mathcal{O}(\alpha_s^2)$
2. Cadre de calcul perturbatif complet :
   • Traitement de jusqu'à 805 diagrammes de Feynman NNLO (cas n≥5)
   • Application systématique de la méthode des opérateurs de projection pour extraire les coefficients d'appariement
   • Réduction par IBP ramenant les résultats à 6 intégrales maîtres (1 seule au NLO)
3. Vérification d'application pratique : Ces coefficients ont été utilisés par la collaboration OpenLat pour la première détermination à partir de la QCD sur réseau des moments de la PDF du quark de valence du pion 49,50
4. Vérifications croisées multiples :
   • Cohérence avec les résultats NNLO connus pour n=1 (courant vectoriel) et n=2 (opérateur d'énergie cinétique du quark)
   • Vérification de l'indépendance du paramètre de jauge en jauge générale Rξ
   • Vérification de la structure des termes logarithmiques via les équations du groupe de renormalisation

Détails méthodologiques

Définition de la tâche

Objectif : Calculer les coefficients d'appariement nus $\zeta_n^B(t)$ dans la relation d'expansion à court temps d'écoulement $\hat{O}^{rs}_{\{\mu_1\cdots\mu_n\}}(t,x) = \zeta_n^B(t) O_{\{\mu_1\cdots\mu_n\}}(x) + \mathcal{O}(t)$ et les coefficients finis renormalisés $\zeta_n(t,\mu) = (Z_n \zeta_n^B(t))(\mu)$

Entrées :

• Définition des opérateurs de twist-deux : $\hat{O}^{rs}_{\{\mu_1\cdots\mu_n\}}(t,x) = \mathring{Z}_\chi \bar{\chi}^r(t,x) \gamma_{\{\mu_1} \overleftrightarrow{D}_{\mu_2}(t,x) \cdots \overleftrightarrow{D}_{\mu_n\}}(t,x)\chi^s(t,x)$
• Équations du flux de gradient définissant les champs d'écoulement $B_\mu^a(t)$ et $\chi(t)$
• Précision du calcul en boucles : $\mathcal{O}(\alpha_s^2)$

Sorties : Expressions analytiques $\zeta_n(t,\mu)$ incluant les dépendances en facteurs de couleur $C_F, C_A, T_R$ et nombre de saveurs $n_f$

Architecture de calcul

1. Méthode des opérateurs de projection

Définition de l'opérateur de projection (Eq. 2.9) : $P^{\mu_1\cdots\mu_n}[O] = N_n(\partial_p - \partial_k)^{n-1}_{\mu_2\cdots\mu_n} \text{Tr}\left[\gamma_{\mu_1} \langle\psi_q(p)|O(p+k)|\psi_{q'}(k)\rangle\right]\bigg|_{\substack{k=p=0\\m_q=m_{q'}=0}}$

Le facteur de normalisation $N_n$ est choisi de sorte qu'au niveau de l'arbre : $P^{\mu_1\cdots\mu_n}[O_{\{\mu_1\cdots\mu_m\}}] = \delta_{nm}$

Propriété clé : Toutes les intégrales scalaires dans les diagrammes en boucles d'ordre supérieur sont des intégrales sans échelle, qui s'annulent en régularisation dimensionnelle, donc cette relation est valide à tous les ordres.

2. Extraction des coefficients d'appariement

Application de l'opérateur de projection à la relation d'expansion à court temps d'écoulement : $\zeta_n^B(t) = P^{\mu_1\cdots\mu_n}[\hat{O}_{\{\mu_1\cdots\mu_n\}}(t)]$

Cela transforme le problème en calcul d'éléments de matrice hors-couche de l'opérateur d'écoulement.

3. Procédure de calcul des diagrammes de Feynman

Génération des diagrammes : Utilisation de qgraf pour générer les diagrammes de Feynman

• Nombre de diagrammes contributifs (Tableau 1) :
  • LO : 1 (pour tous les n)
  • NLO : 10-15
  • NNLO : 375-805

Identification topologique : tapir/exp identifie les topologies intégrales et insère les règles de Feynman

Simplification algébrique : FORM effectue la simplification algébrique de Lorentz et Dirac jusqu'à réduction aux seules intégrales scalaires

Réduction IBP : Kira et FireFly effectuent la réduction intégrale

• Au niveau NNLO, réduction à 6 intégrales maîtres (dont les solutions analytiques sont connues)
• Au NLO, une seule intégrale maître est nécessaire

Points d'innovation technique

1. Cadre automatisé : Génération automatique des règles de Feynman jusqu'aux opérateurs de projection, traitant les structures tensorielles de Lorentz complexes
2. Traitement perturbatif des champs d'écoulement :
   • Les champs d'écoulement $B_\mu(t), \chi(t)$ satisfont des équations différentielles du premier ordre (Appendice A)
   • Conditions aux limites : $B_\mu(t=0)=A_\mu, \chi(t=0)=\psi$
   • Méthode systématique de construction des solutions perturbatives établie dans les Refs. 36,37,53
3. Choix du schéma de renormalisation :
   • Constante de couplage fort en schéma $\overline{\text{MS}}$
   • Champ de quark d'écoulement en "schéma annelé" (ringed scheme), défini par la condition (Eq. 2.5) : $\mathring{Z}_\chi \sum_{q=1}^{n_f} \langle\bar{\chi}_q(t)\overleftrightarrow{\not{D}}(t)\chi_q(t)\rangle = -\frac{n_c n_f}{(4\pi t)^2}$
   • Ceci est une renormalisation non-minimale incluant une partie finie $\zeta_\chi(t,\mu)$
4. Vérification de l'invariance de jauge : Calcul en jauge générale $R_\xi$, vérifiant que le paramètre de jauge ξ s'annule complètement dans le résultat final

Configuration expérimentale

Stratégie de vérification du calcul

1. Comparaison avec les résultats connus

• n=1 (courant vectoriel) : Cohérence avec le résultat NNLO de la Ref. 73
• n=2 (opérateur d'énergie cinétique du quark) : Cohérence avec le résultat NNLO de la Ref. 56
• Ces deux cas particuliers fournissent des vérifications d'exactitude indépendantes

2. Vérification par les équations du groupe de renormalisation

Les coefficients d'appariement satisfont l'équation RG : $\mu^2 \frac{d}{d\mu^2}\zeta_n(t,\mu) = \gamma_n \zeta_n(t,\mu)$

où la dimension anormale $\gamma_n = -\alpha_s \sum_{i\geq0} \alpha_s^i \gamma_{n,i}$ est connue jusqu'à trois boucles (Tableau 2).

Méthode de vérification : Vérifier que le coefficient de $L_{\mu t} = \ln(2\mu^2 t) + \gamma_E$ dans les résultats calculés est cohérent avec $\gamma_n$.

3. Invariance de jauge

Calcul complet en jauge générale $R_\xi$, vérifiant que la dépendance du paramètre de jauge ξ s'annule complètement dans le résultat final.

Détails d'implémentation

Régularisation dimensionnelle : $D = 4 - 2\epsilon$

Chaîne d'outils de calcul :

• qgraf 58,59 : Génération des diagrammes de Feynman
• tapir/exp 60-63 : Identification topologique
• FORM 64,65 : Algèbre symbolique
• Kira 67,68 : Réduction IBP
• FireFly 69,70 : Méthode du corps fini

Intégrales maîtres : 6 intégrales maîtres au niveau NNLO, dont les expressions analytiques sont connues 56

Résultats expérimentaux

Résultats principaux

Coefficients d'appariement NNLO complets

L'article fournit les expressions analytiques complètes pour n=1 à 6 (Eqs. 3.4-3.9). Pour n=2 par exemple :

$\zeta_2(t,\mu) = 1 + \alpha_s C_F\left(\frac{1}{72} - \ln 2 - \frac{3}{4}\ln 3 + \frac{2}{3}L_{\mu t}\right)$$+ \alpha_s^2\left\{\frac{c_\chi^{(2)}}{16} + C_F^2[\cdots] + C_A C_F[\cdots] + C_F n_f T_R[\cdots]\right.$$\left. + L_{\mu t}[\cdots] + L_{\mu t}^2\left[\frac{2}{9}C_F^2 + \frac{11}{36}C_A C_F - \frac{1}{9}C_F n_f T_R\right]\right\} + \mathcal{O}(\alpha_s^3)$

Incluant :

• Fonctions transcendantes : $\ln 2, \ln 3, \text{Li}_2(1/4), \zeta(2)=\pi^2/6$
• Structure des facteurs de couleur : $C_F^2, C_A C_F, C_F T_R n_f$
• Dépendance logarithmique : $L_{\mu t}, L_{\mu t}^2$

Analyse des caractéristiques structurelles

1. Confirmation NLO : Les coefficients NLO pour tous les n sont en accord complet avec l'expression générale en n de la Ref. 33
2. Structure logarithmique :
   • Le coefficient de $L_{\mu t}$ est contrôlé par la dimension anormale à une boucle $\gamma_{n,0}$
   • Le coefficient de $L_{\mu t}^2$ est déterminé conjointement par $\gamma_{n,0}$ et le coefficient de fonction bêta $\beta_0$
   • Vérification via l'équation RG (Eq. 3.10)
3. Partie finie : Contient des coefficients numériques complexes, par exemple pour n=1 : $C_F^2: -\frac{41}{128} - \frac{5}{32}\zeta(2) + \frac{3}{8}\ln 2 + \frac{1}{4}\ln^2 2 - \frac{3}{32}\ln 3 + \frac{3}{2}\text{Li}_2(1/4)$

Caractéristiques numériques

Pour la QCD ($C_F=4/3, C_A=3, T_R=1/2$), avec des paramètres typiques :

• Amplitude des corrections NLO : $\mathcal{O}(10\%)$ (dépendant de $\alpha_s$ et $L_{\mu t}$)
• Corrections NNLO : Fournissant des corrections supplémentaires de $\mathcal{O}(\alpha_s^2) \sim 1-5\%$
• Les termes logarithmiques améliorés $L_{\mu t}^2$ deviennent importants aux temps d'écoulement courts

Applications pratiques

Étude de la PDF du pion de la collaboration OpenLat 49,50 :

• Première utilisation de ces coefficients NNLO pour déterminer à partir de la QCD sur réseau les moments d'ordre supérieur de la PDF du quark de valence du pion
• Utilisation de configurations de champs de jauge générées avec des fermions de Wilson stabilisés 44-47
• Les résultats préliminaires ont été rapportés dans la Ref. 48

Cas particulier : Courant vectoriel (n=1)

La conservation du courant vectoriel exige $\gamma_1=0$ (à tous les ordres), donc : $\zeta_1(t,\mu) = 1 + \alpha_s C_F\left(\frac{1}{8} - \ln 2 - \frac{3}{4}\ln 3\right) + \mathcal{O}(\alpha_s^2)$

Le terme NNLO ne contient pas $L_{\mu t}$, uniquement des termes constants, ce qui est une propriété particulière des courants conservés.

Travaux connexes

Méthodes de calcul des PDFs en QCD sur réseau

1. Méthode traditionnelle des moments 27-32 :
   • Calcul direct des éléments de matrice d'opérateurs localisés
   • Limitée par les divergences en puissance, pratiquement n≤4
   • La collaboration ETM a calculé $\langle x \rangle, \langle x^2 \rangle, \langle x^3 \rangle$ pour les pions et kaons 31,32
2. Méthode des quasi-PDFs 20 :
   • Utilisation de la théorie effective à grande impulsion
   • Nécessite un contrôle soigneux des effets de courte distance et de l'impulsion du hadron
3. Méthode des pseudo-PDFs 21,22 :
   • Connexion systématique via les fonctions de corrélation euclidiennes à la physique du cône de lumière
   • Nécessite une extrapolation à grande impulsion
4. Méthode OPE 19,25 :
   • Cadre du développement en produit d'opérateurs
   • La Ref. 25 étudie la méthode du champ amélioré en jauge de Coulomb

Applications du flux de gradient en QCD

1. Théorie fondamentale 34-37 :
   • Narayanan-Neuberger et Lüscher établissent les fondations du GFF
   • Preuve de la finitude ultraviolette des opérateurs composites à t>0
2. Tenseur énergie-impulsion 38,41,56 :
   • Suzuki et al. établissent le cadre général de l'expansion à court temps d'écoulement
   • Coefficients d'appariement NNLO du tenseur énergie-impulsion 56
3. Autres opérateurs 42,43,73 :
   • Appariement à deux boucles des opérateurs dipolaires de couleur 42
   • Renormalisation du flux de gradient pour les mélanges de mésons et les durées de vie 43
   • Appariement NNLO des opérateurs bilinéaires de quark 73

Contribution unique de cet article

• Première extension des coefficients d'appariement des opérateurs de twist-deux à NNLO
• Couverture de n=1 à 6, fournissant des inputs théoriques pour le calcul précis des moments de PDFs d'ordre supérieur
• Formation d'un cadre NNLO complet avec le travail NLO de la Ref. 33

Conclusions et discussion

Conclusions principales

1. Calcul NNLO achevé : Calcul réussi des coefficients d'expansion à court temps d'écoulement pour les six premiers opérateurs de twist-deux non-singlets, avec une précision de $\mathcal{O}(\alpha_s^2)$
2. Vérification de l'utilité pratique : Ces coefficients ont déjà été appliqués dans des calculs QCD sur réseau réels pour déterminer les moments de la PDF du pion 49,50
3. Universalité : Les résultats s'appliquent à toutes les combinaisons de saveurs non-singlets (comme $u-d, u+d-2s$, etc.) et à tous les systèmes hadroniques (proton, pion, etc.)
4. Extensions potentielles : Puisque les éléments de matrice non-avant impliquent également des opérateurs de twist-deux, les résultats pourraient être pertinents pour les études des distributions de partons généralisées (GPDs)

Limitations

1. Limité à la partie non-singlet : La partie singlet (impliquant les distributions de gluons) n'a pas encore été calculée en raison de l'expansion de la base d'opérateurs et des modes de mélange plus complexes
2. Plage limitée de n : Actuellement jusqu'à n=6, les moments supérieurs nécessiteraient de traiter plus de diagrammes de Feynman (bien que cela soit en principe possible)
3. Convergence du développement perturbatif :
   • Dans la région de couplage fort ou pour grand $L_{\mu t}$, la série perturbative peut converger lentement
   • Nécessite une vérification numérique réelle dans les calculs sur réseau
4. Sélection du temps d'écoulement : L'expansion à court temps d'écoulement exige que $t$ soit suffisamment petit, mais les calculs sur réseau nécessitent $\sqrt{8t} \gg a$ (espacement du réseau), créant une fenêtre pratique limitée

Directions futures

1. Extension à la partie singlet :
   • Les auteurs affirment explicitement : "ne pas anticiper d'obstacles conceptuels"
   • Nécessite de traiter le mélange quark-gluon et une base d'opérateurs plus grande
2. Moments d'ordre supérieur : Extension à n>6, bien que le nombre de diagrammes de Feynman augmente, la méthode est établie
3. Calcul à trois boucles (N³LO) :
   • Extrêmement difficile techniquement
   • Nécessite l'attente de plus de résultats analytiques pour les intégrales maîtres
4. Applications numériques :
   • Application à plus de systèmes hadroniques (proton, nucléon, etc.)
   • Comparaison avec d'autres méthodes de détermination des PDFs (quasi-PDFs, etc.)
5. Effets non-perturbatifs : Étude des contributions non-perturbatives des corrections en puissance $\mathcal{O}(t^n)$

Évaluation approfondie

Points forts

1. Rigueur technique :
   • Vérifications multiples indépendantes (résultats connus pour n=1,2, équations RG, invariance de jauge)
   • Calcul complet à deux boucles NNLO, traitant des centaines de diagrammes de Feynman
   • Fourniture de fichiers supplémentaires lisibles par ordinateur (formats Mathematica et Python)
2. Impact pratique :
   • Déjà adopté par les calculs QCD sur réseau réels 49,50
   • Comble une lacune théorique cruciale dans l'application du GFF au calcul des PDFs
   • La précision NNLO est essentielle pour réduire les erreurs systématiques
3. Contributions méthodologiques :
   • Application systématique de la méthode des opérateurs de projection
   • Établit une feuille de route claire pour les calculs de la partie singlet et d'ordres supérieurs
   • La chaîne d'outils automatisée peut s'appliquer à d'autres opérateurs
4. Présentation claire :
   • Structure bien organisée, cadre théorique, méthodes de calcul et résultats clairement séparés
   • L'appendice fournit les détails complets de la définition du flux de gradient et des constantes de renormalisation
   • Les expressions analytiques complètes sont fournies, facilitant l'utilisation dans les recherches ultérieures

Insuffisances

1. Analyse numérique limitée :
   • Pas de comparaison de la taille des coefficients d'appariement sous différents paramètres
   • Manque d'évaluation quantitative des corrections NNLO par rapport au NLO
   • Pas de discussion sur les indicateurs de convergence de la série perturbative
2. Intuition physique insuffisante :
   • Pas de discussion approfondie sur la raison pour laquelle certaines fonctions transcendantes (comme $\text{Li}_2(1/4)$) apparaissent
   • Manque d'explication physique de la structure des facteurs de couleur
   • Pas d'exploration des motifs systématiques entre différents n
3. Guidance pour l'application pratique :
   • Pas de recommandations quantitatives pour la sélection du temps d'écoulement
   • Manque de discussion sur l'estimation des erreurs (erreurs de troncature, dépendance du schéma, etc.)
   • Pas de guidance sur le choix de l'échelle de renormalisation optimale
4. Absence de la partie singlet : Bien que les auteurs reconnaissent cette limitation, pour une détermination complète des PDFs, la partie singlet (contribution des gluons) est tout aussi importante

Évaluation de l'impact

1. Contribution au domaine :
   • Élevée : Ceci est un input théorique clé pour la détermination précise des moments de PDFs d'ordre supérieur en QCD sur réseau
   • Rend possible pour la première fois le calcul des moments n>4
   • Fournit une base pour la phénoménologie QCD précise future
2. Valeur pratique :
   • Application directe : Déjà utilisée dans l'étude du pion d'OpenLat
   • Applicabilité large : Peut s'appliquer à tous les hadrons (proton, neutron, nucléon, etc.)
   • Valeur à long terme : Avec l'amélioration de la précision des calculs sur réseau, les coefficients NNLO deviendront de plus en plus importants
3. Reproductibilité :
   • Excellente : Fourniture de fichiers supplémentaires complets (Zeta.m et Zeta.py)
   • Description détaillée des méthodes de calcul et de la chaîne d'outils
   • Les vérifications croisées multiples avec les résultats connus renforcent la crédibilité
4. Potentiel de citation :
   • Devrait devenir une référence standard pour les calculs de PDFs sur réseau utilisant le flux de gradient
   • La méthodologie servira de modèle pour les calculs d'ordre supérieur d'autres opérateurs

Scénarios d'application

1. Calculs QCD sur réseau :
   • Toute recherche calculant les moments des PDFs utilisant la méthode du flux de gradient
   • Études de structure hadronique nécessitant une précision NNLO
   • Détermination des PDFs pour les hadrons légers (proton, pion, kaon, etc.)
2. Phénoménologie QCD :
   • Fournir des contraintes sur réseau pour les ajustements globaux des PDFs
   • Vérification précise du théorème de factorisation QCD
   • Prédictions théoriques pour les futures installations expérimentales (EIC, etc.)
3. Développement méthodologique :
   • Comme référence et modèle pour les calculs de la partie singlet
   • Calcul des coefficients d'appariement pour d'autres opérateurs composites (opérateurs liés aux GPDs)
   • Référence technique pour les calculs à plus haute boucle
4. Scénarios non applicables :
   • Partie singlet (gluon) des moments de PDFs (nécessite un travail ultérieur)
   • Situations nécessitant des corrections en puissance non-perturbatives
   • Temps d'écoulement très grand (dépassant la plage de validité de l'expansion à court temps)

Références clés

33 A. Shindler, "Moments of parton distribution functions of any order from lattice QCD," Phys. Rev. D 110 (2024) L051503 - Travail original proposant la méthode du flux de gradient pour les moments des PDFs

36 M. Lüscher and P. Weisz, "Perturbative analysis of the gradient flow in non-abelian gauge theories," JHEP 02 (2011) 051 - Fondations de la théorie des perturbations du flux de gradient

49,50 A. Francis et al., "Moments of parton distributions functions of the pion from lattice QCD using gradient flow" - Première application des résultats de cet article

56 R. V. Harlander, Y. Kluth, and F. Lange, "The two-loop energy–momentum tensor within the gradient-flow formalism," Eur. Phys. J. C 78 (2018) 944 - Calcul NNLO antérieur pour le cas n=2

73 J. Borgulat et al., "Short-flow-time expansion of quark bilinears through NNLO QCD," JHEP 05 (2024) 179 - Calcul NNLO antérieur pour le cas n=1

Évaluation globale : Ceci est un article de calcul théorique techniquement solide avec une haute valeur pratique. Il fournit les coefficients d'appariement NNLO critiques pour la détermination précise des moments de PDFs d'ordre supérieur en QCD sur réseau, comblant une lacune importante dans l'application de la méthode du flux de gradient. Bien que l'analyse numérique et l'intuition physique pourraient être améliorées, le calcul rigoureux, les vérifications multiples et la preuve d'application pratique démontrent son importance. Avec l'amélioration de la précision des calculs sur réseau et la construction de nouvelles installations expérimentales comme l'EIC, l'impact à long terme de ce travail continuera à se manifester.