Self-confinement of ultra-high-energy nuclei in cosmic filaments: implications for the UHECR spectrum and composition

Résumé technique : Confinement autonome des noyaux de très haute énergie dans les filaments cosmiques

Énoncé du problème
Les observations du spectre et de la composition en masse des rayons cosmiques d'ultra-haute énergie (UHECR) par l'Observatoire Pierre Auger indiquent que la population dominant au-dessus de la « cheville » ($E \sim 5 \times 10^{18}$ eV) nécessite un spectre de source exceptionnellement dur ($\gamma \lesssim 1$) et supprimé aux faibles rigidités (autour de l'échelle de l'EV). De tels spectres durs sont difficiles à concilier avec les mécanismes standards d'accélération par choc diffusif, qui prédisent généralement $\gamma \gtrsim 2$. Alors que des ajustements phénoménologiques antérieurs ont encodé cette caractéristique comme une propriété intrinsèque de l'accélérateur, des explications alternatives suggèrent que le spectre est remodelé par des effets de transport. Plus précisément, les scénarios de confinement autonome proposent que les UHECR s'échappant génèrent une turbulence magnétique via l'instabilité de streaming non résonante (NRSI), retardant leur propre libération de manière dépendante de la rigidité. Cependant, les explorations détaillées antérieures de ce mécanisme se sont concentrées principalement sur des compositions purement protoniques, laissant les implications pour la composition nucléaire mixte observée inexplorées.

Méthodologie
Les auteurs étendent le cadre de confinement autonome, précédemment développé pour les protons, à une composition nucléaire mixte. Le modèle suppose des sources d'UHECR intégrées dans le milieu intergalactique (IGM), spécifiquement au sein des amas de galaxies et des filaments cosmiques.

• Croissance de l'instabilité : Le courant électrique porté par les noyaux s'échappant excite la NRSI, amplifiant la turbulence magnétique. Le taux de croissance dépend de la luminosité de la source et de la luminosité effective dépendante de la composition ($L \equiv g_Z L$). La turbulence sature lorsque le champ amplifié atteint un niveau où l'échelle de diffusion résonante correspond au rayon de Larmor des particules.
• Modèle de transport : La région confinée est modélisée à l'aide d'une approche phénoménologique de « Boîte fuyante ». Le temps d'échappement est déterminé par la compétition entre la diffusion dans la turbulence auto-générée et l'advection par le plasma de fond. Une rigidité critique ($R_{cut}$) est définie où le temps de diffusion sur une longueur de cohérence magnétique ($\lambda_B$) dépasse l'âge de la source. En dessous de cette rigidité, le flux s'échappant est supprimé de manière exponentielle.
• Interactions : Le modèle inclut explicitement la photodésintégration des noyaux lourds et la production de paires Bethe-Heitler au sein de la région de confinement. Les noyaux lourds sont traités comme subissant une chaîne de désintégration en nombre de masse, produisant des protons secondaires.
• Propagation : Les noyaux s'échappant sont propagés jusqu'à la Terre en utilisant le code SimProp. Le spectre et la composition résultants (moyenne et variance de $X_{max}$) sont comparés aux données d'Auger.
• Signatures multi-messagers : Les auteurs calculent les flux de neutrinos et de rayons gamma cosmogéniques associés générés par les particules confinées dans les régions magnétisées, les comparant aux limites établies par IceCube, Auger, KM3NeT et Fermi-LAT.

Contributions clés

1. Extension à une composition mixte : L'étude généralise le scénario de confinement autonome pour inclure un mélange réaliste d'espèces nucléaires (protons, He, N, Ne, Si, Fe), démontrant comment le mécanisme affecte à la fois la forme spectrale et l'évolution de la composition en masse.
2. Production de protons secondaires : Les auteurs montrent que pendant le confinement prolongé des noyaux lourds, une photodésintégration efficace se produit. Cela produit un flux significatif de protons secondaires qui peuplent naturellement la gamme d'énergie en dessous de la cheville, aidant à expliquer l'évolution de la composition observée sans nécessiter de composante primaire protonique dominante aux énergies les plus élevées.
3. Contraintes sur l'espace des paramètres : L'article cartographie l'espace des paramètres viables de luminosité de source ($L$) et de longueur de cohérence du champ magnétique ($\lambda_B$) requis pour reproduire les données d'Auger, identifiant des plages spécifiques où des spectres d'accélération standards ($\gamma \simeq 2$) peuvent produire les spectres effectifs durs observés.
4. Contraintes multi-messagers : Pour la première fois dans ce contexte, l'étude évalue le fond diffus de rayons gamma et les flux de neutrinos associés au confinement autonome, fournissant de nouvelles vérifications de cohérence par rapport aux limites observationnelles actuelles.

Résultats

• Durcissement spectral : Le modèle reproduit avec succès les principales caractéristiques du spectre et de la composition observés des UHECR pour des longueurs de cohérence du champ magnétique de $\lambda_B \simeq 30\text{--}50$ Mpc. Dans ces scénarios, un spectre d'accélération standard ($\gamma \simeq 2$) est transformé en un spectre effectif dur déduit des observations en raison de la suppression dépendante de la rigidité du flux s'échappant en dessous de $R \sim 1$ EV.
• Composition : Pour $\lambda_B \simeq 30\text{--}50$ Mpc, le modèle prédit une contribution négligeable de protons primaires aux énergies les plus élevées. La composition observée est dominée par des noyaux lourds aux hautes énergies et des protons secondaires (issus de la photodésintégration de noyaux lourds confinés) en dessous de la cheville. Cette configuration reproduit l'évolution et la variance de $X_{max}$ observées par Auger sans nécessiter de spectre de source inversé.
• Convergence des contraintes :
  • Neutrinos : Le flux de neutrinos prédit provenant des particules confinées reste compatible avec les limites supérieures actuelles d'IceCube et d'Auger. Pour le cas de confinement le plus extrême ($\lambda_B = 50$ Mpc), le flux approche le niveau déduit d'un événement récent de KM3NeT, bien que cette comparaison soit notée comme phénoménologique plutôt que définitive.
  • Rayons gamma : Le fond diffus de rayons gamma prédit provenant des cascades électromagnétiques approche le fond isotrope mesuré aux hautes énergies. Cela suggère que le fond diffus de rayons gamma pourrait fournir une contrainte plus forte que les neutrinos sur les configurations de confinement les plus extrêmes.
• Viabilité astrophysique : La suppression requise à $R \sim 1$ EV est réalisable pour des luminosités de source dans la plage $10^{43}\text{--}10^{46}$ erg s$^{-1}$ et des longueurs de cohérence de $5\text{--}50$ Mpc. Ces paramètres sont cohérents avec des sources hébergées dans des amas de galaxies se propageant à travers les filaments cosmiques (WHIM), à condition que l'intensité initiale du champ magnétique soit $\gtrsim 1$ nG.

Signification
L'article conclut que le confinement autonome des UHECR dans les filaments cosmiques offre un mécanisme physiquement motivé pour transformer les spectres d'accélération standards en spectres effectifs durs requis par les observations. Dans cette perspective, la dureté apparente des sources d'UHECR n'est pas une propriété intrinsèque de l'accélérateur, mais l'empreinte observable de la turbulence magnétique auto-générée dans le réseau cosmique. Ce mécanisme explique naturellement la suppression aux faibles rigidités et l'évolution spécifique de la composition (noyaux lourds aux hautes énergies, protons secondaires aux énergies plus faibles) sans invoquer de physique d'accélération exotique. L'étude souligne que les futures mesures de la composition en masse (par exemple via AugerPrime) et une sensibilité améliorée aux neutrinos et rayons gamma diffus seront cruciales pour tester la viabilité de ce scénario et le rôle de la turbulence auto-générée dans la propagation des particules de plus haute énergie.