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

Technische Samenvatting: Zelfopsluiting van ultra-hoge-energie kernen in kosmische filamenten

Probleemstelling
Waarnemingen van het spectrum en de massasamenstelling van kosmische straling met ultra-hoge energie (UHECR) door het Pierre Auger-observatorium geven aan dat de populatie die dominant is boven de "enkel" ($E \sim 5 \times 10^{18}$ eV) een bron-spectrum vereist dat ongewoon hard is ($\gamma \lesssim 1$) en onderdrukt bij lage rigiditeiten (rond de EV-schaal). Dergelijke harde spectra zijn moeilijk te verenigen met standaard diffussieve schokversnellingmechanismen, die doorgaans $\gamma \gtrsim 2$ voorspellen. Hoewel eerdere fenomenologische fits dit kenmerk hebben gecodeerd als een intrinsieke eigenschap van de versneller, suggereren alternatieve verklaringen dat het spectrum wordt hervormd door transporteffecten. Specifiek stellen scenario's van zelfopsluiting voor dat ontsnappende UHECR's magnetische turbulentie genereren via de niet-resonante stromingsinstabiliteit (NRSI), waardoor hun eigen vrijlating op een manier die afhankelijk is van de rigiditeit wordt vertraagd. Eerdere gedetailleerde verkenningen van dit mechanisme richtten zich echter voornamelijk op composities bestaande uit alleen protonen, waardoor de implicaties voor de waargenomen gemengde nucleaire samenstelling ongetest bleven.

Methodologie
De auteurs breiden het zelfopsluitingskader, dat eerder voor protonen is ontwikkeld, uit naar een gemengde nucleaire samenstelling. Het model neemt UHECR-bronnen aan die zijn ingebed in het intergalactische medium (IGM), specifiek binnen sterrenstelselclusters en kosmische filamenten.

• Groei van instabiliteit: De elektrische stroom die wordt gedragen door ontsnappende kernen exciteert de NRSI, waardoor magnetische turbulentie wordt versterkt. De groeisnelheid hangt af van de bronluminositeit en de samenstellingsafhankelijke effectieve luminositeit ($L \equiv g_Z L$). De turbulentie verzadigt wanneer het versterkte veld een niveau bereikt waarbij de resonante verstrooiingsschaal overeenkomt met de Larmorstraal van het deeltje.
• Transportmodel: Het opgesloten gebied wordt gemodelleerd met een fenomenologische "Leaky Box"-benadering. De ontsnappingstijd wordt bepaald door de competitie tussen diffusie in de zelfgegenereerde turbulentie en advectie door het achtergrondplasma. Een kritieke rigiditeit ($R_{cut}$) wordt gedefinieerd waarbij de diffusietijd over een magnetische coherentielengte ($\lambda_B$) de leeftijd van de bron overschrijdt. Onder deze rigiditeit wordt de ontsnappende flux exponentieel onderdrukt.
• Interacties: Het model bevat expliciet fotodesintegratie van zware kernen en Bethe-Heitler-paarproductie binnen het opsluitingsgebied. Zware kernen worden behandeld als ondergaand een vervalrij in massagetal, waarbij secundaire protonen worden geproduceerd.
• Propagatie: Ontsnappende kernen worden met de code SimProp naar de Aarde gepropageerd. Het resulterende spectrum en de samenstelling (gemiddelde en variantie van $X_{max}$) worden vergeleken met Auger-gegevens.
• Multi-messenger handtekeningen: De auteurs berekenen de bijbehorende kosmogene neutrino- en gammastralingfluxen die worden gegenereerd door deeltjes die zijn opgesloten in de gemagnetiseerde gebieden, en vergelijken deze met limieten van IceCube, Auger, KM3NeT en Fermi-LAT.

Belangrijkste bijdragen

1. Uitbreiding naar gemengde samenstelling: De studie generaliseert het zelfopsluitingsscenario om een realistische mix van nucleaire soorten te omvatten (protonen, He, N, Ne, Si, Fe), en demonstreert hoe het mechanisme zowel de spectrale vorm als de evolutie van de massasamenstelling beïnvloedt.
2. Productie van secundaire protonen: De auteurs tonen aan dat tijdens de langdurige opsluiting van zware kernen efficiënte fotodesintegratie optreedt. Dit produceert een aanzienlijke flux van secundaire protonen die op natuurlijke wijze de energiebereik onder de enkel bevolken, wat helpt de waargenomen evolutie van de samenstelling te verklaren zonder dat een dominante primaire protoncomponent bij de hoogste energieën vereist is.
3. Beperkingen van de parameterruimte: Het artikel in kaart de levensvatbare parameterruimte van bronluminositeit ($L$) en magnetische-veld coherentielengte ($\lambda_B$) die nodig is om Auger-gegevens te reproduceren, en identificeert specifieke bereiken waar standaard versnellingsspectra ($\gamma \simeq 2$) de waargenomen harde effectieve spectra kunnen opleveren.
4. Multi-messenger beperkingen: Voor het eerst in deze context evalueert de studie het diffuse gammastralingachtergrond en neutrinofluxen die geassocieerd zijn met zelfopsluiting, en biedt nieuwe consistentiecontroles tegenover huidige waarnemingslimieten.

Resultaten

• Spectrale harding: Het model reproduceert succesvol de belangrijkste kenmerken van het waargenomen UHECR-spectrum en de samenstelling voor magnetische-veld coherentielengten van $\lambda_B \simeq 30\text{--}50$ Mpc. In deze scenario's wordt een standaard versnellingsspectrum ($\gamma \simeq 2$) omgezet in het harde effectieve spectrum dat uit waarnemingen wordt afgeleid, als gevolg van de rigiditeitsafhankelijke onderdrukking van de ontsnappende flux onder $R \sim 1$ EV.
• Samenstelling: Voor $\lambda_B \simeq 30\text{--}50$ Mpc voorspelt het model een verwaarloosbare bijdrage van primaire protonen bij de hoogste energieën. De waargenomen samenstelling wordt gedomineerd door zware kernen bij hoge energieën en secundaire protonen (van fotodesintegratie van opgesloten zware kernen) onder de enkel. Deze configuratie reproduceert de door Auger waargenomen evolutie en variantie van $X_{max}$ zonder dat een omgekeerd bronspectrum vereist is.
• Convergentie van beperkingen:
  • Neutrino's: De voorspelde neutrinoflux van opgesloten deeltjes blijft compatibel met huidige bovengrenzen van IceCube en Auger. Voor het meest extreme opsluitingsgeval ($\lambda_B = 50$ Mpc) nadert de flux het niveau dat wordt afgeleid uit een recente KM3NeT-gebeurtenis, hoewel deze vergelijking wordt genoteerd als fenomenologisch en niet als definitief.
  • Gammastraling: De voorspelde diffuse gammastralingachtergrond van elektromagnetische cascades nadert de gemeten isotrope achtergrond bij hoge energieën. Dit suggereert dat de diffuse gammastralingachtergrond een sterkere beperking kan opleggen dan neutrino's aan de meest extreme opsluitingsconfiguraties.
• Astrofysische levensvatbaarheid: De vereiste onderdrukking bij $R \sim 1$ EV is haalbaar voor bronluminositeiten in het bereik $10^{43}\text{--}10^{46}$ erg s$^{-1}$ en coherentielengten van $5\text{--}50$ Mpc. Deze parameters zijn consistent met bronnen die zijn gehost in sterrenstelselclusters die zich voortbewegen door kosmische filamenten (WHIM), mits de initiële magnetische veldsterkte $\gtrsim 1$ nG is.

Betekenis
Het artikel concludeert dat UHECR-zelfopsluiting in kosmische filamenten een fysiek gemotiveerd mechanisme biedt om standaard versnellingsspectra om te zetten in de harde effectieve spectra die door waarnemingen worden vereist. In dit beeld is de schijnbare hardheid van UHECR-bronnen geen intrinsieke eigenschap van de versneller, maar de waarneembare afdruk van zelfgegenereerde magnetische turbulentie in het kosmische web. Dit mechanisme verklaart op natuurlijke wijze de onderdrukking bij lage rigiditeit en de specifieke evolutie van de samenstelling (zware kernen bij hoge energie, secundaire protonen bij lagere energie) zonder exotische versnellingsfysiek aan te roepen. De studie benadrukt dat toekomstige metingen van de massasamenstelling (bijvoorbeeld via AugerPrime) en een verbeterde gevoeligheid voor diffuse neutrino's en gammastraling cruciaal zullen zijn voor het testen van de levensvatbaarheid van dit scenario en de rol van zelfgegenereerde turbulentie in de propagatie van deeltjes met de hoogste energie.