Energy-momentum tensor form factors and spin density distribution in the nucleon calculated in a quantized Skyrme model with vector mesons

Cette étude utilise un modèle de Skyrme quantifié avec des mésons vectoriels pour démontrer que le choix de la jauge pseudo, entre les tenseurs énergie-impulsion canonique et de Belinfante, entraîne des différences significatives dans la distribution spatiale locale du spin et de l'impulsion du nucléon, tout en préservant ses propriétés globales.

L'essentiel

🏗️ Le Nucléon : Une Maison avec deux plans d'architecte différents

Imaginez que le proton (ou le neutron, ensemble appelés "nucléons") est une petite maison très complexe, faite de briques invisibles (les quarks et les gluons). Les physiciens veulent comprendre comment cette maison tient debout : où est la pression ? Comment l'énergie est-elle répartie ? Et surtout, comment tourne-t-elle sur elle-même (le "spin") ?

Pour répondre à ces questions, les auteurs de cet article, Kenji Fukushima et Tomoya Uji, ont utilisé un modèle mathématique appelé le modèle Skyrmion (une sorte de "tourbillon" de champs qui imite le proton).

Mais voici le tour de magie : ils ont dessiné cette maison avec deux plans d'architecte différents pour voir si cela change la façon dont on la perçoit.

1. Les deux plans d'architecte (Les "Pseudogauges")

En physique, il existe deux façons principales de définir comment l'énergie et le mouvement sont stockés dans un objet. C'est un peu comme si vous deviez décrire le trafic dans une ville :

• Le Plan A (Canonique) : Il compte chaque voiture individuellement. Il distingue clairement qui conduit (le "spin", ou rotation intrinsèque) et qui se déplace sur la route (le "moment orbital").
• Le Plan B (Belinfante) : Il regroupe tout le trafic en un seul flux fluide. Il ne fait plus la différence entre le conducteur et la voiture ; tout est mélangé dans un seul grand mouvement.

En physique, ces deux plans sont mathématiquement équivalents pour le total (si vous additionnez tout, vous obtenez le même résultat). Mais ils diffèrent radicalement sur la répartition locale (ce qui se passe à un endroit précis de la maison).

2. La découverte : La maison semble différente selon le plan

Les chercheurs ont calculé la structure du proton avec ces deux plans et ont fait une découverte fascinante :

• La pression et la force (Le D-term) :
  Imaginez que le proton est un ballon de baudruche gonflé. Le plan A dit : "La pression est douce ici, forte là." Le plan B dit : "Non, la pression est très forte ici, presque nulle là."
  Résultat : Les deux plans donnent des valeurs très différentes pour la force qui maintient le proton ensemble. Cela signifie que si on essaie de cartographier la "pression" à l'intérieur du proton, la réponse dépend de quel "plan d'architecte" on utilise !

• La rotation (Le Spin) :
  C'est là que c'est le plus surprenant.

  • Avec le Plan A (Canonique), on voit clairement une partie du proton qui "tourne sur elle-même" (comme une toupie) et une autre partie qui "tourne autour" (comme la Lune autour de la Terre).
  • Avec le Plan B (Belinfante), la partie "toupie" disparaît complètement ! Tout le mouvement de rotation est décrit comme si tout le proton tournait autour de son centre, comme une planète.
  • L'analogie : C'est comme regarder une danseuse. Le Plan A vous dit : "Elle tourne sur ses talons ET elle tourne autour de la scène." Le Plan B vous dit : "Oubliez les talons, elle fait juste un grand tour autour de la scène." Les deux décrivent le même mouvement global, mais l'expérience visuelle locale est totalement différente.

3. Pourquoi est-ce important ?

Les physiciens vont bientôt utiliser un super microscope géant appelé le Collisionneur Électron-Ion (EIC) pour prendre des photos 3D des protons. Ils veulent savoir exactement où se trouve l'énergie et comment le proton tourne.

Ce papier nous met en garde : Il n'y a pas une seule "vraie" carte de la pression ou du spin à l'intérieur du proton. Il y a plusieurs cartes valables, selon comment on choisit de définir les règles du jeu mathématique.

• Si on regarde les données brutes de l'EIC, on ne peut pas savoir immédiatement quel "plan d'architecte" la nature utilise réellement.
• Les deux plans sont justes, mais ils racontent des histoires différentes sur la façon dont la matière est organisée à l'intérieur.

En résumé

C'est comme si vous essayiez de décrire la répartition de la chaleur dans une soupe.

• Une méthode dit : "La chaleur vient des légumes qui bougent."
• L'autre méthode dit : "La chaleur vient du bouillon qui tourne."
• Au final, la soupe est chaude partout de la même façon (le total est identique), mais si vous mettez une cuillère à un endroit précis, vous sentirez la chaleur différemment selon la méthode de calcul.

La leçon : La structure interne de la matière est plus subtile qu'il n'y paraît. Selon la "lunette" mathématique que vous portez, le proton peut sembler être un objet rigide qui tourne sur lui-même, ou un fluide qui tourne en bloc. Les deux sont vrais, mais ils nous donnent des intuitions différentes sur la mécanique quantique.

Résumé technique

Voici un résumé technique détaillé de l'article « Energy-momentum tensor form factors and spin density distribution in the nucleon calculated in a quantized Skyrme model with vector mesons » par Kenji Fukushima et Tomoya Uji.

1. Problématique et Contexte

L'objectif principal de cette étude est d'explorer l'ambiguïté dite de pseudogauge dans la définition du tenseur énergie-impulsion (TEP ou EMT en anglais) et son impact sur la description locale de la structure du nucléon, en particulier la densité de spin et de moment cinétique.

• Le problème de la non-unicité : En théorie des champs, le tenseur énergie-impulsion n'est pas unique. On peut ajouter des termes de divergence totale (super-potentiels) qui ne modifient pas les charges conservées globales (comme l'énergie totale ou le moment cinétique total), mais qui altèrent radicalement les densités locales de moment et de spin.
• Les deux constructions principales :
  • Le tenseur canonique ($T^{\mu\nu}_{can}$), dérivé du théorème de Noether, qui est généralement non symétrique et contient explicitement un courant de spin.
  • Le tenseur de Belinfante ($T^{\mu\nu}_{Bel}$), symétrisé, qui absorbe le courant de spin dans le tenseur lui-même, rendant le moment cinétique purement orbital dans sa forme locale.
• L'enjeu physique : Bien que les observables de premier twist (comme les GPDs non polarisées mesurées au futur collisionneur EIC) soient insensibles à cette ambiguïté (car elles ne sondent que la composante symétrique $T^{++}$), les distributions spatiales tridimensionnelles de pression, de force de cisaillement et de spin dépendent du choix de jauge. Il est crucial de comprendre comment ces choix affectent l'interprétation de la structure mécanique du nucléon.

2. Méthodologie

Les auteurs utilisent un modèle théorique contrôlé pour comparer explicitement les deux formulations dans un même cadre dynamique :

• Modèle utilisé : Un modèle de Skyrmion quantifié étendu aux mésons vectoriels ($\rho$ et $\omega$). L'ajout des mésons vectoriels est essentiel car ils portent une densité de spin intrinsèque, permettant de distinguer les contributions orbitales et de spin de manière non triviale.
• Approche semi-classique :
  1. Solution statique : Résolution numérique des équations du mouvement pour les profils de champ (pions, mésons $\rho$ et $\omega$) dans le cadre d'une ansatz "hérisson" (hedgehog) avec nombre baryonique $B=1$.
  2. Rotation collective : Introduction de coordonnées collectives pour la rotation adiabatique afin de générer des états de nucléon avec un spin et un isospin définis ($J=I=1/2$). Cela induit des champs vectoriels supplémentaires dépendant du temps.
  3. Quantification : Les modes de rotation sont quantifiés pour obtenir les éléments de matrice nucléaires.
• Calcul des TEP :
  • Construction du tenseur canonique via la formule de Noether.
  • Construction du tenseur de Belinfante via la variation métrique (ou symétrisation de Belinfante).
• Extraction des observables :
  • Calcul des éléments de matrice dans le cadre de Breit ($\Delta^0 = 0$).
  • Utilisation de formules d'inversion (transformées de Fourier-Bessel) pour passer des densités spatiales aux facteurs de forme $A(t)$, $D(t)$, $J(t)$ et, pour le cas canonique, au facteur de forme antisymétrique $G(t)$.

3. Contributions Clés et Résultats

L'étude fournit une réalisation concrète de l'ambiguïté de pseudogauge au niveau des densités locales et des facteurs de forme hors du front (off-forward).

A. Facteurs de forme du TEP

• Facteur de forme de l'impulsion $A(t)$ : Les deux formulations donnent $A(0)=1$, respectant la conservation de l'impulsion. Les courbes $A_{can}(t)$ et $A_{Bel}(t)$ sont très proches.
• Facteur de forme mécanique $D(t)$ (Terme D) : C'est ici que la différence est la plus marquée.
  • $D_{can}(0) \approx -1.30$
  • $D_{Bel}(0) \approx -3.88$
  • Conclusion : La valeur du terme D, qui encode la distribution de pression interne et la stabilité mécanique, dépend fortement du choix de la jauge. Les contraintes expérimentales actuelles (via les GPDs de twist-2) ne peuvent pas distinguer ces deux réalités, introduisant une incertitude systématique dans l'interprétation mécanique.
• Facteur de forme de moment cinétique $J(t)$ et $G(t)$ :
  • Pour Belinfante, $G(t) = 0$ par construction et $J_{Bel}(0) = 1/2$ (spin total).
  • Pour le cas canonique, $T^{0i} \neq T^{i0}$. Le facteur de forme $J_{can}(t)$ extrait de la partie symétrique ne satisfait pas $J_{can}(0) = 1/2$ seul.
  • Un nouveau facteur de forme antisymétrique $G(t)$ apparaît. La relation de conservation s'écrit : $L_{can}(0) + S_{can}(0) = J_{can}(0) + \frac{1}{2}G_{can}(0) = 1/2$.

B. Distributions de densité de spin et de moment

Les auteurs visualisent les densités dans l'espace 3D pour un nucléon polarisé selon l'axe $z$ :

• Cas Canonique : Le moment cinétique total est réparti entre une densité de moment cinétique orbital ($L_z$) et une densité de spin intrinsèque ($S_z$). La densité de spin $S_z$ reste non nulle au centre du nucléon.
• Cas Belinfante : Le courant de spin est absorbé. La densité de moment cinétique total $J_z$ prend une forme purement "orbitale" ($x^i T^{0j} - x^j T^{0i}$). Par conséquent, la densité s'annule au centre du nucléon ($r=0$), contrairement au cas canonique.
• Signification : Ces résultats montrent que la répartition spatiale du spin (où se trouve le spin ?) n'est pas une observable unique, mais dépend de la définition du tenseur énergie-impulsion utilisée.

C. Densités de boost

Les auteurs examinent également les densités liées au boost ($L^{00i}$ et $S^{00i}$). Ils confirment que dans le référentiel du repos, les charges totales associées s'annulent, mais les distributions locales diffèrent selon la jauge.

4. Importance et Implications

• Interprétation des données de l'EIC : L'article met en garde contre l'interprétation directe des reconstructions 3D de la pression et du spin comme des vérités absolues. Les données de diffusion Compton virtuel profond (DVCS) aux premier twist ne résolvent pas l'ambiguïté de pseudogauge. Une analyse complète nécessitera des observables de twist supérieur (GPDs d'ordre supérieur, GTMDs) qui sont sensibles aux termes de divergence.
• Modélisation théorique : Ce travail fournit un "laboratoire" théorique pour quantifier la taille de l'incertitude systématique liée au choix de jauge. Il démontre que la séparation entre spin et moment cinétique orbital n'est pas unique au niveau local.
• Rôle des mésons vectoriels : L'étude souligne l'importance des degrés de liberté vectoriels pour capturer la densité de spin intrinsèque, ce qui est manquant dans les modèles de Skyrmion à pions seuls.

En résumé, Fukushima et Uji démontrent que la structure mécanique et spinorielle du nucléon, telle qu'elle est décrite par les densités locales, est intrinsèquement ambiguë en raison de la liberté de pseudogauge. Leur modèle quantifie cette ambiguïté, offrant des repères cruciaux pour l'interprétation future des données de haute précision de l'Electron-Ion Collider (EIC).