Modelling Galactic neutrino emission: contributions from massive star clusters and interstellar cosmic rays

Technische Samenvatting: Modellering van Galactische Neutrino-emissie van Massieve Sterrenclusters en Interstellaire Kosmische Straling

Probleemstelling
De oorsprong van de meest energetische kosmische stralen (CR's) in het Universum, met name die met energieën die de PeV-schaal bereiken (PeVatrons), blijft onzeker. Hoewel supernova-restanten (SNRs) traditioneel worden beschouwd als de primaire versnellers, suggereren $\gamma$-straalobservaties dat zij moeite hebben om PeV-energieën te bereiken, wat de noodzaak voor alternatieve bronklassen creëert. Jonge Massieve Sterrenclusters (YMSCs) zijn voorgesteld als levensvatbare PeVatron-kandidaten vanwege de versterkte magnetische turbulentie gegenereerd door botsende sterrenwinden en meerdere SNR-gebeurtenissen. Gelijktijdig heeft het IceCube Observatory een hoogenergetische neutrinoflux gedetecteerd vanuit het Galactische Vlak. De huidige analyses vertrouwen echter op templates die uitsluitend rekening houden met de puur diffuse neutrino-emissie gegenereerd door de interactie van de CR-zee met interstellaire gas, waarbij deze modellen de bijdrage van onopgeloste hadronische versnellers, zoals de YMSC-populatie, niet expliciet verklaren, wat de specifieke signaturen van deze bronnen potentieel maskeert.

Methodologie
De auteurs presenteren een nieuw theoretisch kader om de verwachte neutrino-emissie van de Melkweg te berekenen, waarbij twee onderscheidende componenten worden geïntegreerd:

1. Synthetische YMSC-populatie: Om de beperking te overwinnen dat YMSCs alleen in de lokale zonnige nabijheid ($\lesssim 2$ kpc) geobserveerd kunnen worden, genereerden de auteurs 100 onafhankelijke realisaties van de Galactische YMSC-populatie. Deze synthetische populatie bemonster clustermassa's ($10^3 - 6,3 \times 10^4 M_\odot$), leeftijden ($<30$ Myr) en ruimtelijke verdelingen (volgens spiraalarmen, de Galactische staaf en de Lokale Spur) op basis van observationele beperkingen.
2. Deeltjesversnelling en Transport: Voor elke synthetische cluster simuleert het model de stellaire populatie (Hoofdreeks, Wolf-Rayet en Supernova-progenitoren) om de collectieve mechanische luminositeit van winden en supernovae te berekenen. De auteurs modelleren de resulterende wind-geblazen bellen, waarbij ze de gasdichtheid en magnetische veldprofielen berekenen.
   • Versnellingsmechanismen: Deeltjesversnelling wordt gemodelleerd op twee locaties: de collectieve windterminatieschok (WTS) en supernova-restant (SNR) schokken.
   • Diffusiescenario's: De propagatie van CR's binnen de clusterbellen wordt berekend onder drie verschillende turbulentie-aannames: Kolmogorov, Kraichnan en Bohm-diffusie.
   • Neutrino-productie: De neutrinofluxen worden berekend uit de hadronische botsingen van versnelde protonen met het doelgas (voornamelijk in de dichte, opgezette schil van de bellen).
3. Diffuse Galactische Component: Het model berekent ook de diffuse neutrino-emissie voortvloeiend uit de interactie van de Galactische CR-zee met interstellaire gas. Hiervoor wordt gebruikgemaakt van CR-propagatiemodellen (die spectrale verharding in het binnenste deel van de Melkweg inbegrepen) en gaskaarten (GALPROP en een alternatief 3D-model).
4. Template Generatie: De totale voorspelde flux is de som van de YMSC-bijdrage en de door CR geïnduceerde diffuse component. Deze worden gesynthetiseerd in ruimtelijke en spectrale templates voor vergelijking met IceCube-data.

Belangrijkste Bijdragen

• Nieuwe Neutrino-templates: Het artikel biedt de eerste uitgebreide neutrino-templates die expliciet de emissiviteit van de YMSC-populatie integreren naast de standaard CR-geïnduceerde diffuse component.
• Synthetische Populatiemodellering: Het biedt een statistisch robuuste schatting van de cumulatieve neutrino-output van de gehele Galactische YMSC-populatie, in plaats van enkel te vertrouwen op de enkele nabijgelegen, opgeloste clusters.
• Multi-component Analyse: Het raamwerk maakt onderscheid tussen de "in-situ" productie van neutrino's binnen clusters en de achtergrond diffuse emissie, waardoor een scheiding van bronbijdragen mogelijk is.
• Parametersensitiviteit: De studie verkent de afhankelijkheid van de voorspelde flux van deeltjesdiffusie-regimes (Kolmogorov, Kraichnan, Bohm) en CR-spectrale verharding in de binnenste Melkweg.

Resultaten

• Detecteerbaarheid van Individuele Bronnen: De voorspelde flux van individuele YMSCs blijkt marginaal detecteerbaar, zelfs door detectoren op kubieke kilometer-schaal (bijv. IceCube, KM3NeT). Bijgevolg wordt de cumulatieve emissie van de YMSC-populatie verwacht als een onopgeloste diffuse component te verschijnen.
• Flux Magnitude: De normalisatie van de modellen is consistent met de IceCube best-fit van bestaande Galactische templates. De auteurs stellen vast dat de onopgeloste bijdrage van cluster-emissie mogelijk niet verwaarloosbaar is.
• Diffusie-afhankelijkheid:
  • In het Bohm-diffusiescenario (meest efficiënte opsluiting), domineren YMSCs de neutrinoflux onder $\sim 50$ TeV, wat potentieel de flux overschat vergeleken met de Kraichnan- en Kolmogorov-scenario's als er geen spectrale verharding wordt aangenomen.
  • In de Kolmogorov- en Kraichnan-scenario's domineert de door CR geïnduceerde diffuse component over het algemeen of levert deze vergelijkbare bijdragen aan de YMSC-component, met name bij hogere energieën.
• Ruimtelijke Verdeling: De YMSC-bijdrage is geconcentreerd in het binnenste Galactische longitudinale bereik ($|l| \le 30^\circ$, de Galactische Ridge), terwijl de door CR geïnduceerde emissie domineert bij grotere longitudes.
• Spectrale Verharding: Het opnemen van spectrale verharding van de CR-zee in de binnenste Melkweg vermindert de relatieve bijdrage van YMSCs in de Kolmogorov- en Kraichnan-scenario's, waardoor de totale voorspelde spectrum beter in overeenstemming komt met de IceCube-observaties bij $\sim 100$ TeV.
• Consistentie: De voorspellingen van het model voor de $\gamma$-straal tegenhanger van de hadronische botsingen blijven consistent met de huidige diffuse $\gamma$-straalobservaties in het TeV-PeV bereik.

Betekenis en Claims
De auteurs beweren dat hun werk een noodzakelijke evolutie vormt in de theoretische modellering van de Galactische neutrino-emissie. Door verder te gaan dan louter diffuse templates en de specifieke emissiviteit van de YMSC-populatie op te nemen, bieden zij een instrument om de twee componenten van de Galactische neutrino-emissie te "ontwarren". Het artikel stelt dat het toepassen van deze nieuwe templates op toekomstige neutrino-datasets de gemeenschap in staat zal stellen de rol van sterrenclusters in extreme CR-versnelling en neutrino-productie te beperken. De auteurs benadrukken dat hun resultaten erop wijzen dat de onopgeloste bijdrage van clusters waarschijnlijk niet verwaarloosbaar is en dat een gezamenlijke passing van bronbijdragen en de diffuse component cruciaal is voor het identificeren van Galactische versnellers en het oplossen van de oorsprong van kosmische straling. Het werk claimt niet definitief de bron van het IceCube-signaal te hebben geïdentificeerd als YMSCs, maar biedt het theoretische kader dat vereist is om deze hypothese te toetsen aan observationele data.