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

Technische Samenvatting: Gecombineerde Multi-Periode Fermi-LAT en LHAASO Beperkingen op Axion-achtige Deeltjes van Mrk 421

Probleemstelling
Axion-achtige deeltjes (ALPs) zijn goed gemotiveerde uitbreidingen van het Standaardmodel die foton-ALP-oscillaties in externe magnetische velden kunnen induceren, wat potentieel de spectra van hoogenergetische gammastraling van verre bronnen kan veranderen. Hoewel eerdere studies de ALP-parameters hebben beperkt, blijven significante delen van de parameterruimte onbeperkt. Een grote uitdaging in deze zoektochten is de behandeling van systematische onzekerheden, met name met betrekking tot het Extragalactisch Achtergrondlicht (EBL), dat gammastraling dempt. Traditionele benaderingen behandelen EBL-absorptie en foton-ALP-oscillaties vaak als onafhankelijke, multiplicatieve factoren, waarbij ze de niet-triviale interactie tussen beide verwaarlozen. Bovendien gaan analyses over meerdere tijdperken (epochs) er vaak van uit dat observatieperioden statistisch onafhankelijk zijn, waarbij zij zwakke residuele correlaties negeren die voortkomen uit gedeelde bronvariabiliteit, instrumentele systematiek en modelveronderstellingen. Deze simplificaties kunnen leiden tot bevoordeelde beperkingen of kunstmatig verscherpte betrouwbaarheidsintervallen. Bovendien vertoont de blazar Mrk 421 een zwakke correlatie tussen zijn GeV (Fermi-LAT) en TeV (LHAASO) emissies, wat duidt op verschillende fysieke oorsprongen die de interpretatie van breedbandige data voor ALP-zoektochten bemoeilijken.

Methodologie
De auteurs presenteren een uitgebreid kader om ALP-parameters ($m_a$ en $g_{a\gamma}$) te beperken met behulp van vijf gelijktijdige observationele tijdperken van Mrk 421 van maart 2021 tot maart 2024, waarbij data van Fermi-LAT en LHAASO worden gecombineerd.

1. Propagatieformalisme: In tegen tegenstelling tot gefactoriseerde benaderingen, gebruikt deze studie een zelfconsistent foton-propagatiekader waarbij EBL-absorptie direct in de foton-ALP-mengingshamiltoniaan wordt opgenomen via een absorptieterm ($i/2\lambda_\gamma$). De EBL-intensiteit wordt behandeld met een normalisatiefactor ($K_{EBL}$) als een nuisance-parameter met een Gaussische prior ($\mu=1.0, \sigma=0.2$) om onzekerheden te accounten.
2. Bronmodellering: Twee verschillende jet-emissiemodellen worden getest om de zwakke GeV-TeV correlatie in Mrk 421 aan te pakken:
   • Een twee-zone hybride model met een compacte binnenste zone (flaring VHE-emissie) en een uitgebreide buitenste zone (quiescente emissie).
   • Een single-zone hadronisch model, fysiek plausibel gezien de zwakke correlatie en de verschillende variabiliteitsamplitudes.
     Beide modellen maken gebruik van helicale en verwarde (tangled) magnetische veldcomponenten en propageren fotonen door de bron, het intergalactische medium (IGMF) en het Galactische magnetische veld (GMF).
3. Statistisch Kader: De auteurs ontwikkelen een gezamenlijke profile-likelihood analyse. Om de beperking van het behandelen van tijdperken als onafhankelijk aan te pakken, introduceren zij een Gaussische copula formalisme met een conservatieve correlatiecoëfficiënt ($\rho = 0.03$). Dit modelleert zwakke inter-epoch correlaties via getransformeerde $\Delta\chi^2$ excess distributies. De analyse optimaliseert nuisance-parameters (intrinsieke spectrale parameters en $K_{EBL}$) globaal onder de no-ALP hypothese voordat de ALP-massa en koppeling worden geprofileerd.
4. Dataverwerking: Fotonoverlevingskansen worden vooraf berekend op een vierdimensionale grid (energie, ALP-parameters, EBL-normalisatie) met behulp van de gammaALPs code en geïnterpoleerd tijdens de likelihood-evaluatie. Intrinsieke spectrale modellen worden per epoch geselecteerd met behulp van het Akaike Informatie Criterium gecorrigeerd voor eindige steekproefomvang (AICc).

Belangrijkste Bijdragen

• Verenigd Propagatiemodel: Het artikel implementeert een zelfconsistente behandeling van foton-ALP-menging en EBL-absorptie binnen de propagatievergelijking, waardoor de biases van multiplicatieve factorisatie worden vermeden.
• Gecorreleerde Multi-Epoch Analyse: De introductie van een Gaussische copula om zwakke correlaties tussen observationele perioden te behandelen, biedt een statistisch consistente methode voor het combineren van multi-epoch datasets, wat kunstmatige verscherping van betrouwbaarheidsintervallen voorkomt.
• Breedbandige Beperkingen: De studie maakt gebruik van de unieke breedbandige dekking van Fermi-LAT (GeV) en LHAASO (multi-TeV tot PeV) om ALP-parameters over een breed massabereik (0.1–500 neV) te verkennen.
• Robuuste Modellering: Door zowel twee-zone als single-zone hadronische scenario's te testen, houdt het werk rekening met de complexe emissiefysica van Mrk 421, waardoor de beperkingen niet overmatig afhankelijk zijn van een enkele emissie-aanname.

Resultaten
Met behulp van het gezamenlijke profile-likelihood kader leiden de auteurs 95% betrouwbaarheidsniveau (CL) bovengrenzen af voor de foton-ALP-koppeling ($g_{a\gamma}$):

• Lage-massa regime ($m_a \lesssim 1$ neV): De beperkingen zijn zeer strikt, met een ondergrens van $5.0 \times 10^{-14}$ GeV$^{-1}$ en een bovengrens van $8.33 \times 10^{-14}$ GeV$^{-1}$. Dit gebied is gevoelig voor CAST-beperkingen.
• Intermediaire tot Hoge Massa ($0.1 - 500$ neV):
  • Voor het twee-zone hybride model is de conservatieve bovengrens $g_{a\gamma} < 6.50 \times 10^{-12}$ GeV$^{-1}$.
  • Voor het single-zone hadronische model is de conservatieve bovengrens $g_{a\gamma} < 7.30 \times 10^{-12}$ GeV$^{-1}$.
• Gedrag: De uitsluitingscontouren vertonen een karakteristiek "jump-plateau" kenmerk. Een scherpe transitie vindt plaats rond $m_a \approx 1$ neV, waar de bovengrens met ongeveer een orde van grootte toeneemt, gevolgd door een geleidelijke relaxatie bij hogere massa's. De ondergrens blijft stabiel over het gehele massabereik.
• Robuustheid: Tests met alternatieve correlatiecoëfficiënten ($\rho = 0, 0.15, 0.3$) tonen aan dat variaties in de beperkingen binnen minder dan één orde van grootte blijven, wat de stabiliteit van de resultaten bevestigt.

Betekenis
Het artikel beweert dat deze resultaten aanzienlijk profiteren van de breedbandige zeer-hoge-energie dekking van LHAASO en de langdurige monitoring van Mrk 421. Het ontwikkelde kader biedt een robuuste en schaalbare pipeline die geschikt is voor toekomstige multi-messenger ALP-zoektochten met volgende generatie gamma-waarnemingsinstrumenten. Door de zwakke correlaties expliciet te modelleren en de EBL-absorptie zelfconsistent te behandelen, biedt de studie betrouwbaardere en conservatievere beperkingen op de ALP-parameterruimte, waarmee de systematische biases aanwezig in eerdere factoriserende benaderingen worden geadresseerd. Het werk demonstreert dat de zwakke GeV-TeV correlatie in Mrk 421 complexe emissiemodellen (twee-zone of hadronisch) vereist voor nauwkeurige ALP-beperkingen, en dat de voorgestelde statistische methodologie effectief de daarmee gepaard gaande onzekerheden beheert.