RAMSES-MCR: A consistent multi-group treatment of cosmic rays physics in momentum-space with the RAMSES code

Probleemstelling

Cosmische straling (CR) speelt een cruciale rol in de astrofysica, van het beïnvloeden van schokgolven en de thermodynamica van het interstellair medium (ISM) tot het reguleren van galactische winden. Traditionele numerieke simulaties behandelen CR's vaak als een "grijze" benadering: een enkelvoudig fluïdum met één effectieve energie en transportcoëfficiënt.

Deze benadering is echter ontoereikend omdat:

1. Energie-afhankelijkheid: Het transport (diffusie en streaming) en energieverliezen (Coulomb, hadronisch, synchrotron) van CR's zijn sterk afhankelijk van de deeltjessnelheid (momentum).
2. Spectrale evolutie: CR's hebben een breed vermogenswettig spectrum dat zich ruimtelijk en tijdelijk ontwikkelt door verschillende verliesmechanismen.
3. Koppeling met gas: De dynamische terugkoppeling van CR's op het gas (via druk en verwarming) hangt af van het specifieke spectrum, wat in grijze modellen verloren gaat.

Er is behoefte aan een methode die het CR-spectrum oplost in de impulsruimte (momentum space) binnen een adaptief mesh-verfijning (AMR) MHD-code, gekoppeld aan de gasdynamica.

Methodologie: RAMSES-MCR

De auteurs presenteren ramses-mcr, een implementatie van een consistente multi-groep spectrale methode voor CR-protonen in de AMR-code RAMSES. De methode is gebaseerd op een twee-moment formalisme (analoog aan stralingstransport) en werkt in het meebewegende fluïdumstelsel (fluid frame).

Kerncomponenten van de methode:

• Discretisatie in impulsruimte: Het CR-spectrum wordt opgedeeld in logaritmische bins rond een centrale impuls $p_i$. Voor elke bin worden zowel de energiedichtheid ($e_{c,i}$) als de aantaldichtheid ($n_{c,i}$) geëvolueerd.
• Twee-moment vergelijkingen: De methode lost vergelijkingen op voor de dichtheid en de flux (voor zowel aantal als energie). Dit omvat:
  • Advectie met het gas.
  • Anisotrope en isotrope diffusie (afhankelijk van het magnetisch veld).
  • Streaming (instabiliteit) langs magnetische veldlijnen.
  • Radiatieve verliezen (Coulomb, ionisatie, hadronisch) en adiabatische veranderingen.
• Spectrale koppeling: De term $Q_{i}^{\pm}$ beschrijft de overdracht van CR's tussen momentum-bins door verliezen en injectie. Dit zorgt ervoor dat het spectrum correct vervormt (bijv. afvlakken bij lage energieën door Coulomb-verliezen).
• Numerieke implementatie:
  • Gebruik van een operator-splitting benadering.
  • Hyperbolische termen (advectie/diffusie) worden expliciet opgelost, terwijl brontermen (verliezen) impliciet of semi-impliciet worden behandeld voor stabiliteit.
  • De methode is uitbreidbaar naar CR-elektronen (met synchrotronverliezen) en andere deeltjessoorten.

Belangrijkste Bijdragen

1. Consistente Multi-Groep Implementatie: De eerste implementatie van een volledig gekoppelde multi-groep CR-methode in RAMSES die zowel energiedichtheid als aantaldichtheid evolueert, in tegenstelling tot eerdere grijze benaderingen.
2. Fluid-Frame Formulier: De vergelijkingen zijn afgeleid in het fluïdumstelsel, wat nodig is voor consistentie onder adiabatische transformaties en nauwkeuriger is voor kleine waarden van de gereduceerde lichtsnelheid ($\tilde{c}$).
3. Geavanceerde Diffusiebehandeling: Implementatie van spectrale correcties voor diffusiecoëfficiënten. In plaats van één diffusiecoëfficiënt per bin, worden deze spectraal gemiddeld om de verandering in de interne helling van het spectrum door diffusie correct te vangen.
4. Validatie en Toepassing: Uitgebreide validatie tegen analytische oplossingen en toepassing op een complex 3D-supernova-remnant (SNR) scenario.

Resultaten

1. Validatietests (Ideale tests):

• Radiatieve verliezen: De methode reproduceert nauwkeurig de analytische evolutie van het spectrum onder Coulomb- en hadronische verliezen. Het toont de verwachte afvlakking bij lage energieën en behoud van de helling bij hoge energieën.
• Injectie en evenwicht: De simulatie bereikt correct het stationaire evenwichtsspectrum waar injectie en verliezen in balans zijn.
• Adiabatische veranderingen: Tests tonen aan dat een semi-impliciet tijdsintegratieschema de beste resultaten geeft met fouten < 1% bij het behoud van de spectrale vorm tijdens compressie/expansie.
• Diffusie: De methode toont aan dat spectrale gemiddelde diffusiecoëfficiënten nodig zijn om de vervorming van het spectrum correct te modelleren. Hogere energieën diffunderen sneller, wat leidt tot een ruimtelijk gescheiden spectrum.
• Shock Tube: De methode reproduceert correct de dynamiek van een CR-gemodificeerde schokbuis, inclusief de verdunning van de gasdichtheid door streaming-verliezen.

2. Toepassing: Supernova-Remnant (SNR):
De auteurs simuleerden de 3D-expansie van een SNR met anisotrope diffusie en energieverliezen.

• Spectrale redistributie: CR's met hoge energie (TeV) ontsnappen snel uit de restant, terwijl lagere energieën (GeV) langer worden ingesloten. Dit creëert een ruimtelijk variërend spectrum.
• Impulsversterking: De aanwezigheid van CR's zorgt voor een extra "snowplough"-fase waarin de impuls van het SNR blijft toenemen, in plaats van te stabiliseren zoals bij thermische gas alleen.
• Invloed van diffusie: Een lagere diffusiecoëfficiënt leidt tot een grotere impulsversterking (tot 135% ten opzichte van een run zonder CR's), omdat CR's langer in het systeem blijven en druk uitoefenen.
• Anisotropie: De CR's diffunderen langs magnetische veldlijnen, wat leidt tot smalle plumes buiten het SNR, in tegenstelling tot een bolvormige uitbreiding bij isotrope diffusie.

Significantie

RAMSES-MCR biedt een robuust en nauwkeurig instrument om de terugkoppeling van kosmische straling in astrofysische omgevingen te bestuderen. Door het spectrum expliciet op te lossen, kunnen onderzoekers:

• De impact van energie-afhankelijk transport op galactische winden en de vorming van sterren beter begrijpen.
• Realistische voorspellingen doen voor niet-thermische waarnemingen (zoals $\gamma$-straling en radio-synchrotron) die direct gekoppeld zijn aan het lokale CR-spectrum.
• De dynamische evolutie van schokgolven en SNR's modelleren met een fysiek correctere beschrijving van CR-druk en -verliezen dan grijze modellen mogelijk maken.

De methode opent de deur voor toekomstige studies naar de ruimtelijke variatie van diffusiecoëfficiënten en de interactie tussen CR's en koude gaswolken in galactische schijven.