Searching for missing direct photons in heavy-ion collisions with P and CP violation

Résumé technique : Recherche de photons directs manquants dans les collisions d'ions lourds avec violation de P et CP

Énoncé du problème
La production de photons directs dans les collisions relativistes d'ions lourds est traditionnellement modélisée en utilisant une combinaison de la QCD perturbative pour les processus à courte distance et de l'hydrodynamique pour le comportement collectif du plasma de quarks et de gluons (QGP). Bien que ce cadre ait été largement fructueux, il est confronté à un défi persistant connu sous le nom de « puzzle des photons directs manquants ». Plus précisément, les modèles conventionnels ne parviennent pas à décrire simultanément le spectre de photons dans la partie douce et son anisotropie azimutale (flux elliptique, $v_2$). Les approches standards sous-estiment soit le rendement des photons, soit prédisent une variation azimutale trop faible par rapport aux données expérimentales.

Des tentatives précédentes pour résoudre cela se sont concentrées sur le rayonnement synchrotron émis par le QGP dû aux champs magnétiques intenses. Bien que ce mécanisme fournisse des photons supplémentaires et augmente le $v_2$ en raison de l'émission anisotrope perpendiculaire au champ magnétique, le flux elliptique résultant dépasse souvent les observations expérimentales. De plus, l'inclusion de la rotation du plasma tend à amplifier cet effet de champ magnétique, exacerbant la divergence. Cet article postule que le flux elliptique $v_2$ dépend fortement des propriétés chirales du QGP, spécifiquement des coefficients de transport associés aux courants magnétiques chiraux et Hall anormaux, qui violent les symétries P et CP.

Méthodologie
Les auteurs développent un cadre théorique pour calculer le rayonnement de photons en présence de courants anormaux en résolvant l'équation de Dirac modifiée pour les fermions dans un champ magnétique constant $\mathbf{B}$ et en présence de paramètres violant P et CP. Ces paramètres sont introduits via un terme de théorie de champ effectif couplant le courant vectoriel chiral $j^\mu_5$ à un champ pseudoscalaire de fond $\theta$, défini par le vecteur $b_\mu = \partial_\mu \theta$. Le vecteur $b_\mu$ possède deux composantes non nulles dans cette étude : $b_0$ (le potentiel chimique chiral, $\mu_5$) et $b_3$ (le gradient chiral longitudinal le long du champ magnétique).

La contribution méthodologique centrale est la dérivation de fonctions d'onde exactes (spinors) pour les fermions satisfaisant l'équation :
$$(i \not{D} - \gamma_5 \not{b} - m)\psi = 0$$
où $D$ est la dérivée covariante. Les auteurs :

1. Découplent l'équation de Dirac : Ils réduisent l'équation de Dirac du premier ordre en une forme du second ordre pour résoudre les composantes des spinors dans un champ magnétique constant en utilisant la jauge de Landau symétrique.
2. Dériment les spinors exacts : Ils construisent des solutions exactes pour les cas $n=0$ (niveau de Landau le plus bas) et $n>0$, exprimant les spinors en termes de fonctions propres de l'oscillateur harmonique. Ces solutions tiennent compte du mélange des chiralités induit par $b_\mu$.
3. Calculent les taux de rayonnement : En utilisant ces fonctions d'onde exactes, ils calculent l'élément S-matrix pour le rayonnement synchrotron. Ils dérivent le taux d'émission différentiel pour un fermion unique et, par la suite, pour un plasma non interactif en appliquant les poids statistiques de Fermi-Dirac.
4. Simulation numérique : Les auteurs effectuent des calculs numériques pour l'intensité de rayonnement d'une particule unique et le rendement de photons intégré ($v_0$) ainsi que le flux elliptique ($v_2$) pour un plasma de quarks up et down. Ils utilisent des paramètres pertinents pour le QGP, notamment une masse de quark thermique de $m \approx 300$ MeV et des champs magnétiques typiques des collisions d'ions lourds ($|e|B \approx 18000$ MeV$^2$).

Principales contributions

• Fonctions d'onde exactes : L'article fournit les premières solutions exactes de l'équation de Dirac en présence d'un champ magnétique constant, d'un potentiel chimique chiral ($b_0$) et d'un gradient chiral longitudinal ($b_3$). Cela étend les travaux précédents qui étaient souvent restreints à la limite chirale ou ignoraient le gradient longitudinal.
• Inclusion du gradient chiral ($b_3$) : Contrairement à des calculs récents similaires (par exemple, Shovkovy et Wang) qui se concentraient principalement sur $b_0$, ce travail incorpore explicitement $b_3$, associé à la composante longitudinale du courant Hall anomal, et démontre son impact distinct sur les relations de dispersion et la structure des spinors.
• Dérivation de l'effet magnétique chiral : Les fonctions d'onde dérivées sont utilisées pour dériver indépendamment le courant magnétique chiral, vérifiant que le courant calculé est indépendant de $b_3$, comme attendu théoriquement.
• Résolution de la divergence de $v_2$ : L'étude démontre que si le rayonnement synchrotron seul prédit un flux elliptique trop important, l'inclusion des paramètres violant P et CP ($b_0$ et $b_3$) réduit significativement le $v_2$.

Résultats
Les résultats numériques pour un plasma QGP avec $T=300$ MeV et un champ magnétique $|e|B = m_\pi^2$ donnent les résultats suivants :

• Rendement des photons ($v_0$) : La présence de $b_0$ ou $b_3$ non nuls augmente le nombre total de photons émis d'un facteur d'environ 4 par rapport au cas standard. Cette augmentation sature lorsque les paramètres atteignent un ordre de grandeur de $10^{-2}m$ (environ 3 MeV).
• Flux elliptique ($v_2$) : La découverte la plus significative est la réduction du coefficient de flux elliptique. En l'absence de paramètres chiraux, le rayonnement synchrotron produit un $v_2$ qui excède les données expérimentales. Cependant, l'activation de $b_0$ ou $b_3$ supprime significativement le $v_2$.
• Comparaison avec les données : Les valeurs réduites de $v_2$, combinées au rendement de photons accru, amènent les prédictions théoriques dans une plage comparable aux données expérimentales de la collaboration PHENIX (centralité 20-40 %). Les auteurs notent que le résultat pour $v_2$ est robuste car il ne dépend pas des échelles de longueur et de temps arbitraires utilisées pour normaliser le rendement total.
• Scission d'hélicité : L'étude observe une scission de l'intensité de rayonnement basée sur les hélicités $s$ et $u$ des particules initiales, une caractéristique accentuée par la présence de $b_0$ et $b_3$.

Signification et revendications
L'article soutient que l'interaction entre le rayonnement synchrotron et les coefficients de transport chiral offre une solution viable au puzzle des photons directs manquants. Plus précisément, les auteurs affirment que le courant magnétique chiral, circulant parallèlement au champ magnétique, provoque l'émission de photons par les quarks relativistes principalement dans la direction vers l'avant, contrecarrant ainsi l'émission perpendiculaire caractéristique du rayonnement synchrotron standard. Ce mécanisme réduit naturellement le flux elliptique $v_2$ à des niveaux cohérents avec les observations tout en augmentant simultanément le rendement des photons.

Les auteurs présentent modestement leur travail comme un cadre théorique cohérent qui comble le fossé entre les modèles hydrodynamiques et les données expérimentales. Ils suggèrent qu'un plasma en rotation avec un déséquilibre chiral (un scénario réaliste dans les collisions d'ions lourds) pourrait pleinement résoudre le puzzle, notant que l'ajout de la rotation à leur hamiltonien est une extension directe. Ils reconnaissent que leur inclusion d'une masse de fermion finie (masse thermique) est un aspect novateur par rapport aux calculs de fermions sans masse, car elle empêche le découplage de l'équation de Dirac et permet une cinématique changeant la chiralité. L'article conclut en identifiant l'extension aux paramètres chiraux transverses ($b_1, b_2$) comme un objectif futur nécessitant une étude dédiée.