Beyond the Local Void: A data-driven search for the origins of the Amaterasu particle

Technische Samenvatting: Voorbij de Local Void: Een datagestuurde zoektocht naar de oorsprong van het Amaterasu-deeltje

Probleemstelling
De oorsprong van ultra-hoogenergetische kosmische stralen (UHECR's), met name die boven de $10^{20}$ eV, blijft onbekend. Hoewel eerdere inspanningen een link hebben gelegd tussen UHECR's en catalogi van starburst-stelsels (SBG's) en actieve galactische kernen (AGN), hebben deze modellen moeite om gebeurtenissen bij de hoogste energieën te verklaren. Het "Amaterasu"-deeltje, gedetecteerd door de Telescope Array met een energie van $244 \pm 29$ (stat.) $^{+51}_{-76}$ (syst.) EeV, vormt een specifieke uitdaging: de aankomstrichting wijst naar de Local Void, een gebied met een lage dichtheid zonder bekende actieve bronnen. Eerdere analyses met behulp van 1D-backtracking of eenvoudige afwijkingsschattingen suggereren dat het deeltje niet afkomstig is van bekende actieve sterrenstelsels, wat leidt tot hypothesen over exotische fysica (bijv. schending van Lorentz-invariantie, superzware donkere materie) of transiënte bronnen. Een rigoureuze, datagestuurde methode is vereist om het volume van het universum in kaart te brengen dat consistent is met dergelijke extreme gebeurtenissen, terwijl complexe propagatie-effecten, magnetische afwijkingen en reconstructie-onzekerheden worden meegewogen.

Methodologie
De auteurs introduceren een op simulatie gebaseerd inferentie-framework dat gebruikmaakt van Approximate Bayesian Computation (ABC) om de oorsprong van individuele UHECR's te beperken. De methode integreert realistische 3D-propagatiemodellering met Bayesiaanse parameterschatting, waarbij de noodzaak om expliciet een likelihood-functie te definiëren wordt vermeden.

• Propagatiemodellering: Het framework maakt gebruik van CRPropa 3 om de propagatie van UHECR's door de 3D-ruimte te simuleren. Het bevat alle relevante interacties: foto-pionproductie, fotodisintegratie, elektron-paarproductie, nucleair verval en adiabatische verliezen.
• Magnetische velden:
  • Extragalactisch Magnetisch Veld (EGMF): Gemodelleerd als een Gaussisch willekeurig veld met een Kolmogorov-turbulentiespectrum op een rooster van 25 kpc. De parameters omvatten de wortel-kwadraat-gemiddelde veldsterkte ($B_{rms}$) en de coherentielengte ($L_c$), toegewezen met log-uniforme priors.
  • Galactisch Magnetisch Veld (GMF): Gemodelleerd met behulp van het best passende UF23-basismodel voor het coherente veld en JF12 Planck-getunede componenten voor de turbulente velden.
• Parameterruimte: De inferentie richt zich op zes vrije parameters: de positie van de bron (Galactische lengte $l$, breedtegraad $b$ en afstand $D_{src}$), de bronenergie ($E_{src}$) en de eigenschappen van het EGMF ($B_{rms}$, $L_c$).
  • Samenstelling van de bron: Er wordt vanuit gegaan dat het deeltje een ijzerkern is bij de bron om vergelijking met eerdere studies te vergemakkelijken en om de oorsprong buiten de Local Void te verkennen (wat grote afwijkingen vereist).
  • Priors: $E_{src}$ volgt een machtswet-spectrum ($E^{-1}$). De bronafstanden worden uniform gesampled binnen een sferisch volume tot $D_{max}$ (12 Mpc voor nominale energie, 15 Mpc voor lage energie). De prior voor de bronrichting is een von Mises-Fisher distributie gecentreerd op de teruggevolgde richting van een ijzerkern, rekening houdend met de verwachte GMF- en EGMF-afwijkingen.
• Inferentieproces (ABC):
  1. Parameters worden gesampled uit de priors.
  2. $10^6$ deeltjes worden isotroop uitgezonden vanuit de gesamplede bronlocatie en gepropageerd naar de aarde.
  3. Gesimuleerde gebeurtenissen worden vergeleken met de geobserveerde Amaterasu-data (aankomstrichting en gereconstrueerde energie) met behulp van een afstandmetriek.
  4. Parametersets worden geaccepteerd als de gesimuleerde gebeurtenis binnen $3\sigma$ van de geobserveerde energie en richting valt.
  5. Geaccepteerde sets worden gewogen op basis van de waarschijnlijkheid van de overeenkomst, wat resulteert in posterieure distributies.
  6. Twee energiescenario's worden geanalyseerd: een nominaal geval ($E_{det} = 244$ EeV) en een lagere systematische ondergrens ($E_{det} = 168$ EeV).

Belangrijkste Resultaten
Het toepassen van dit framework op het Amaterasu-deeltje levert de volgende bevindingen op:

• Bronkandidaten: De posterieure distributies onthullen een bredere set nabijgelegen bronkandidaten dan in eerdere analyses werd gevonden. Hoewel de Local Void een gebied van belang blijft, sluiten de resultaten bronnen binnen de Lokale Groep of nabijgelegen sterrenstelsels niet uit.
• Massa-afhankelijkheid: De geaccepteerde simulaties worden gedomineerd door zware kernen (massa-getal $A > 28$). Voor het nominale energiegeval komt 77,9% van de geaccepteerde gebeurtenissen overeen met een zware samenstelling, vergeleken met 16,2% voor gemiddelde en 5,9% voor lichte samenstellingen.
• Afstandsbeperkingen:
  • Voor de nominale energie ($244$ EeV) is de meest waarschijnlijke bronafstand $\sim 3$ Mpc. Een secundaire piek verschijnt tussen 10–12 Mpc voor zware kernen.
  • Voor het geval met lagere energie ($168$ EeV) verschuift de meest waarschijnlijke afstand naar $\sim 4$ Mpc voor lichte/gemiddelde massa's, met een langere staart voor zware kernen.
  • De aankomstrichting fungeert als het strengste selectiecriterium, waardoor de effectieve afstandshorizon wordt beperkt ondanks de grote energie-onzekerheden.
• Afwijking: De studie kwantificeert dat het Galactische Magnetische Veld (GMF) de afwijking domineert voor lichtere en gemiddelde massa's, terwijl het Extragalactische Magnetische Veld (EGMF) pas significant bijdraagt voor zware deeltjes bij grotere afstanden of lagere energieën.
• Rigiditeit: De voorkeursrigiditeit van de bron wordt gevonden als $R \approx 10^{19.3}$ V voor de nominale energie en $R \approx 10^{18.85}$ V voor het geval met lage energie.

Betekenis en Claims
Het artikel claimt een fundamenteel framework te bieden voor "simulatie-gebaseerde analyses van kosmische stralen bij de hoogste energieën." De primaire bijdrage is het vermogen om individuele UHECR-observaties direct te vertalen naar interpreteerbare bronbeperkingen door:

1. Gezamenlijk te beperken van energie en aankomstrichting, waarbij correlaties en reconstructie-onzekerheden en modellering natuurlijk worden geïntegreerd.
2. Te marginaliseren over onzekerheden in magnetische velden, broncompositie en reconstructiefouten om 3D-posterieure distributies van mogelijke bronlocaties te produceren.
3. De exclusiviteit van de "Local Void" uit te dagen: Door aan te tonen dat standaard propagatiemodellen met zware kernen consistente resultaten kunnen produceren voor bronnen buiten de Local Void, suggereren de auteurs dat verklaringen met exotische fysica niet strikt noodzakelijk zijn, hoewel ze deze ook niet uitsluiten.

De auteurs blijven bescheiden over de reikwijdte en merken op dat hun huidige implementatie beperkt wordt door computationele restricties (bijv. de keuze voor een specifiek GMF-model en een vaste broncompositie van ijzer) en dat toekomstig werk erop gericht zal zijn een breder scala aan GMF-modellen op te nemen en het deeltype als een vrije parameter te behandelen.