Joint constraint on the propagation origin of the cosmic-ray spectral knee from energy spectrum and anisotropy observations

1. Problématique et Contexte

L'origine du « genou » (knee) du spectre des rayons cosmiques, une rupture de pente observée autour de quelques Pétaélectronvolts (PeV) où l'indice spectral passe d'environ -2,7 à -3,1, reste l'un des mystères non résolus de l'astrophysique des particules. Trois mécanismes principaux sont généralement proposés pour expliquer cette structure :

1. Limites d'accélération : Les sources (comme les restes de supernova) ne peuvent pas accélérer les particules au-delà d'une certaine énergie, dépendante de la charge $Z$.
2. Effets de propagation : Les propriétés de diffusion des rayons cosmiques dans le champ magnétique galactique changent à cette énergie (variation du coefficient de diffusion).
3. Interactions : Des interactions avec le fond de particules ou une nouvelle physique dans l'atmosphère.

La difficulté réside dans le fait que les modèles basés sur l'accélération et ceux basés sur la propagation prédisent tous deux une dépendance à la rigidité ($R = p/Z$) pour le genou, rendant la distinction difficile via le seul spectre énergétique. L'article vise à trancher cette question en testant spécifiquement l'hypothèse d'une origine par propagation, en utilisant conjointement les données du spectre énergétique et de l'anisotropie.

2. Méthodologie

Les auteurs ont développé une approche combinant des modèles de propagation et des données observationnelles récentes (notamment de l'expérience LHAASO et de l'observatoire Pierre Auger).

• Modélisation de la Propagation :

  • Utilisation du code DRAGON pour résoudre l'équation de propagation des rayons cosmiques.
  • Adoption d'un modèle de propagation dépendant de l'espace (SDP - Spatially Dependent Propagation) pour reproduire les structures spectrales à basse énergie (durcissement à ~200 GV, adoucissement à ~14 TeV).
  • Introduction d'une source locale proche pour expliquer les anomalies spectrales à basse énergie.
  • Hypothèse clé du genou : Le modèle postule que le genou provient d'une variation brutale du coefficient de diffusion $D(R)$. Une fonction de diffusion à double loi de puissance est construite avec deux paramètres libres : la rigidité du genou ($R_{knee}$) et la variation de l'indice de diffusion ($\delta_{knee}$).
  • Contrairement aux modèles d'accélération, le spectre d'injection des sources ne contient pas de coupure exponentielle ($R_c \to \infty$).

• Contraintes Observationnelles :

  • Les paramètres $R_{knee}$ et $\delta_{knee}$ sont contraints par les mesures précises du spectre des protons par LHAASO ($R_{knee} \approx 3,3$ PV, $\delta_{knee} \approx 0,99$).
  • Le modèle génère ensuite les prédictions théoriques pour l'anisotropie des rayons cosmiques galactiques dans la région du genou et au-delà.

• Traitement des Données d'Anisotropie :

  • Les données d'anisotropie à ultra-haute énergie (UHE) de l'observatoire Pierre Auger contiennent un mélange de composantes galactiques (GCR) et extragalactiques (EGCR).
  • Les auteurs développent une méthode de soustraction (Appendice A et B) pour isoler la composante galactique pure des données d'Auger, en utilisant les phases opposées des anisotropies galactiques et extragalactiques et les flux mesurés.
  • Trois scénarios de correction sont testés : Modified I (correction sur 3 points hauts), Modified II (correction sur 4 points hauts) et un scénario MCS (Most Conservative Scenario) utilisant les limites supérieures.

• Test Statistique :

  • Un test de Kolmogorov-Smirnov est appliqué pour comparer la distribution des amplitudes d'anisotropie prédites par le modèle de propagation avec les amplitudes observées (corrigées) dans la gamme d'énergie de 2,7 PeV à 2,48 EeV.

3. Résultats Clés

• Reproduction du Spectre : Le modèle de propagation modifié réussit à reproduire avec précision la structure du genou observée dans le spectre des protons par LHAASO et le spectre global (all-particle) des expériences KASCADE, IceTop et LHAASO. Cela confirme que la variation du coefficient de diffusion est une hypothèse viable pour expliquer la forme du spectre.
• Discordance de l'Anisotropie :
  • Bien que le modèle corresponde bien aux données d'anisotropie en dessous de 10 PeV, il existe une déviation significative entre les prédictions du modèle et les observations au-dessus du genou (gamme PeV-EeV).
  • Le modèle de propagation prédit une amplitude d'anisotropie beaucoup plus élevée que celle qui est effectivement observée par Auger après correction de la contamination extragalactique.
• Validation Statistique :
  • Le test de Kolmogorov rejette l'hypothèse nulle (qui stipule que le genou est d'origine propagationnelle) avec un niveau de confiance de 95 % ($CL$).
  • Ce rejet est robuste : il persiste même dans le scénario le plus conservateur (MCS), où l'on maximise l'amplitude possible de l'anisotropie galactique en minimisant la soustraction de la composante extragalactique.

4. Contributions et Significativité

• Preuve contre l'origine par propagation : C'est la première étude à utiliser conjointement le spectre de protons de haute précision de LHAASO et les données d'anisotropie d'Auger pour contraindre rigoureusement l'origine du genou. Le résultat indique que, sous les hypothèses actuelles, un mécanisme purement lié à la propagation (changement de diffusion) ne peut pas expliquer simultanément le spectre et l'anisotropie observés.
• Méthodologie de séparation des composantes : L'article propose une méthode robuste pour extraire la composante galactique de l'anisotropie des rayons cosmiques ultra-énergétiques à partir des données d'Auger, en tenant compte des phases opposées et des flux relatifs.
• Implications pour la physique des rayons cosmiques :
  • Si l'origine par propagation est exclue à 95 % de confiance, cela renforce la nécessité d'envisager des scénarios d'accélération (limites des sources) ou d'interactions pour expliquer le genou.
  • Cela suggère que la transition entre les rayons cosmiques galactiques et extragalactiques pourrait être plus complexe ou se produire à des énergies différentes de ce que les modèles de diffusion simples ne le prédisent.

5. Limites et Perspectives

Les auteurs notent une incertitude majeure : l'extrapolation des anisotropies extragalactiques de l'hémisphère sud (observées par Auger) vers l'hémisphère nord (où se trouve LHAASO) pourrait ne pas être parfaite. Si cette extrapolation est erronée, la soustraction de la composante extragalactique pourrait être incorrecte, affectant la conclusion.

Cependant, même avec des corrections conservatrices, le rejet du modèle de propagation reste significatif. Les auteurs concluent que des mesures futures de plus haute précision, notamment par LHAASO sur les spectres et l'anisotropie de différents noyaux dans la région du genou, seront nécessaires pour confirmer définitivement ce résultat et identifier le mécanisme réel à l'œuvre.