Joint Multi-Period Fermi-LAT and LHAASO Constraints on Axion-Like Particles from Mrk 421 Using Profile Likelihood with Gaussian Copula Correlation

Résumé technique : Contraintes conjointes multi-périodes Fermi-LAT et LHAASO sur les axions-comme des particules à partir de Mrk 421

Énoncé du problème
Les axions-comme des particules (ALPs) sont des extensions bien motivées du Modèle Standard qui peuvent induire des oscillations photon-ALP dans des champs magnétiques externes, modifiant potentiellement les spectres des rayons gamma de haute énergie provenant de sources lointaines. Bien que des études antérieures aient contraint les paramètres des ALPs, des régions significatives de l'espace des paramètres restent non contraintes. Un défi majeur dans ces recherches est le traitement des incertitudes systématiques, particulièrement celles liées au Fond Diffus Extragalactique (EBL), qui atténue les rayons gamma. Les approches traditionnelles traitent souvent l'absorption par l'EBL et les oscillations photon-ALP comme des facteurs indépendants et multiplicatifs, négligeant leur interaction non triviale. De plus, les analyses multi-époques supposent fréquemment que les périodes d'observation sont statistiquement indépendantes, ignorant les faibles corrélations résiduelles découlant de la variabilité commune de la source, des systématiques instrumentales et des hypothèses de modèles. Ces simplifications peuvent conduire à des contraintes biaisées ou à des intervalles de confiance artificiellement resserrés. Par ailleurs, le blazar Mrk 421 présente une faible corrélation entre ses émissions GeV (Fermi-LAT) et TeV (LHAASO), suggérant des origines physiques distinctes qui compliquent l'interprétation des données à large bande pour les recherches d'ALPs.

Méthodologie
Les auteurs présentent un cadre complet pour contraindre les paramètres des ALPs ($m_a$ et $g_{a\gamma}$) en utilisant cinq époques d'observation simultanées de Mrk 421, de mars 2021 à mars 2024, combinant les données de Fermi-LAT et LHAASO.

1. Formalisme de propagation : Contrairement aux approximations factorisées, l'étude emploie un cadre de propagation de photons auto-cohérent où l'absorption par l'EBL est incorporée directement dans l'Hamiltonien de mélange photon-ALP via un terme absorbant ($i/2\lambda_\gamma$). L'intensité de l'EBL est traitée avec un facteur de normalisation ($K_{EBL}$) agissant comme un paramètre de nuisance avec un a priori gaussien ($\mu=1.0, \sigma=0.2$) pour tenir compte des incertitudes.
2. Modélisation de la source : Deux modèles d'émission de jet distincts sont testés pour répondre à la faible corrélation GeV-TeV de Mrk 421 :
   • Un modèle hybride à deux zones comprenant une zone interne compacte (émission VHE en éruption) et une zone externe étendue (émission quiescente).
   • Un modèle hadronique à zone unique, physiquement plausible compte tenu de la faible corrélation et des amplitudes de variabilité distinctes.
     Les deux modèles utilisent des composantes de champs magnétiques hélicoïdaux et emmêlés et propagent les photons à travers la source, le milieu intergalactique (IGMF) et le champ magnétique galactique (GMF).
3. Cadre statistique : Les auteurs développent une analyse de profil de vraisemblance conjointe. Pour remédier à la limitation consistant à traiter les époques comme indépendantes, ils introduisent un formalisme de copule gaussienne avec un coefficient de corrélation conservateur ($\rho = 0.03$). Cela modélise les faibles corrélations inter-époques à travers des distributions d'excès $\Delta\chi^2$ transformées. L'analyse optimise globalement les paramètres de nuisance (paramètres spectraux intrinsèques et $K_{EBL}$) sous l'hypothèse de l'absence d'ALPs avant de profiler la masse et le couplage des ALPs.
4. Traitement des données : Les probabilités de survie des photons sont précalculées sur une grille quadridimensionnelle (énergie, paramètres ALP, normalisation EBL) à l'aide du code gammaALPs et interpolées lors de l'évaluation de la vraisemblance. Les modèles spectraux intrinsèques sont sélectionnés par époque à l'aide du critère d'information d'Akaike corrigé pour la taille finie de l'échantillon (AICc).

Contributions clés

• Modèle de propagation unifié : L'article implémente un traitement auto-cohérent du mélange photon-ALP et de l'absorption par l'EBL au sein de l'équation de propagation, évitant les biais de la factorisation multiplicative.
• Analyse multi-époques corrélée : L'introduction d'une copule gaussienne pour gérer les faibles corrélations entre les périodes d'observation fournit une méthode statistiquement cohérente pour combiner les ensembles de données multi-époques, empêchant le resserrement artificiel des intervalles de confiance.
• Contraintes à large bande : L'étude tire parti de la couverture unique en large bande de Fermi-LAT (GeV) et LHAASO (multi-TeV à PeV) pour sonder les paramètres des ALPs sur une large gamme de masses (0.1–500 neV).
• Modélisation robuste : En testant les scénarios de modèles tant hybrides à deux zones que hadroniques à zone unique, ce travail tient compte de la physique d'émission complexe de Mrk 421, garantissant que les contraintes ne dépendent pas excessivement d'une seule hypothèse d'émission.

Résultats
En utilisant le cadre de profil de vraisemblance conjoint, les auteurs dérivent des limites supérieures de confiance (CL) à 95 % sur le couplage photon-ALP ($g_{a\gamma}$) :

• Régime de faible masse ($m_a \lesssim 1$ neV) : Les contraintes sont très strictes, avec une limite inférieure de $5.0 \times 10^{-14}$ GeV$^{-1}$ et une limite supérieure de $8.33 \times 10^{-14}$ GeV$^{-1}$. Cette région est sensible aux contraintes de CAST.
• Masse intermédiaire à élevée ($0.1 - 500$ neV) :
  • Pour le modèle hybride à deux zones, la limite supérieure conservatrice est $g_{a\gamma} < 6.50 \times 10^{-12}$ GeV$^{-1}$.
  • Pour le modèle hadronique à zone unique, la limite supérieure conservatrice est $g_{a\gamma} < 7.30 \times 10^{-12}$ GeV$^{-1}$.
• Comportement : Les contours d'exclusion présentent une caractéristique de type « saut-plateau ». Une transition abrupte se produit autour de $m_a \approx 1$ neV, où la limite supérieure augmente d'environ un ordre de grandeur, suivie d'un relâchement graduel à des masses plus élevées. La limite inférieure reste stable sur toute la gamme de masses.
• Robustesse : Des tests avec des coefficients de corrélation alternatifs ($\rho = 0, 0.15, 0.3$) montrent que les variations des contraintes restent inférieures à un ordre de grandeur, confirmant la stabilité des résultats.

Signification
L'article affirme que ces résultats bénéficient considérablement de la couverture à large bande très haute énergie de LHAASO et de la surveillance à long terme de Mrk 421. Le cadre développé offre un pipeline robuste et évolutif adapté aux futures recherches d'ALPs avec les observatoires de gamma de prochaine génération. En modélisant explicitement les faibles corrélations et en traitant l'absorption par l'EBL de manière auto-cohérente, l'étude fournit des contraintes plus fiables et conservatrices sur l'espace des paramètres des ALPs, répondant aux biais systématiques présents dans les approches factorisées précédentes. Le travail démontre que la faible corrélation GeV-TeV de Mrk 421 nécessite des modèles d'émission complexes (deux zones ou hadroniques) pour des contraintes d'ALPs précises, et que la méthodologie statistique proposée gère efficacement les incertitudes associées.