Correlating Particle Acceleration Rates with Plasma Conditions in Colliding Wind Binaries

Technische Samenvatting: Correlatie tussen deeltjesversnellingssnelheden en plasmacondities in botsende windbinairsystemen

Probleemstelling
Hoewel Supernova-restanten (SNe) traditioneel worden beschouwd als de primaire versnellers van Galactische kosmische straling tot $\sim$1 TeV, hebben recente waarnemingen door het LHAASO-experiment uitgebreide extreem hoogenergetische (VHE) emissies geïdentificeerd in stervormingsgebieden (SFRs) die niet geassocieerd zijn met SNe. Botsende windbinairsystemen (CWBs) zijn voorgesteld als potentiële bronnen van deze deeltjes, in staat om hadronen te versnellen tot TeV- en PeV-energieën. Echter, de specifieke fysische mechanismen die deze versnelling in de complexe, turbulente omgevingen van CWBs aansturen, blijven onduidelijk. Eerdere studies hebben vastgesteld dat CWBs niet-thermische emissie produceren, maar zij hebben de bijdragen van klassieke Diffusieve Schokversnelling (DSA) versus turbulente, magnetisch gedreven processen (zoals stochastische versnelling, magnetische reconnectie of het spiegel-effect) niet voldoende van elkaar kunnen onderscheiden. Dit werk beoogt te identificeren welke plasmacondities en fysische grootheden het sterkst correleren met deeltjesversnellingssnelheden om de dominante versnellingsmechanismen in CWBs te bepalen.

Methodologie
De auteurs voerden hoogresolutie, driedimensionale magnetohydrodynamische (MHD) simulaties uit van een CWB-systeem, gemodelleerd naar de $\eta$ Carinae-binair tijdens de periastronpassage. De simulatieopzet omvatte:

• MHD-modellering: Het systeem werd geëvolueerd met behulp van de ideale samendrukbare MHD-vergelijkingen, opgelost met de AMUN-code. De binair bestond uit een primaire en secundaire ster met afzonderlijke windparameters (snelheden van 500 km s$^{-1}$ en 2500 km s$^{-1}$, respectievelijk). De simulaties verkenden verschillende plasma-samendrukkingsregimes door de effectieve polytrope index ($\Gamma_{\text{eff}}$) te variëren van adiabatisch (5/3) tot sterk koelend (1.1).
• Resolutie: Adaptive Mesh Refinement (AMR) werd toegepast, waarbij een maximale effectieve resolutie van $960 \times 1920 \times 1920$ cellen werd bereikt om schokken, turbulentie en magnetische reconnectiestructuren nauwkeurig te vangen.
• Testdeeltjesintegratie: Na de MHD-evolutie werden $10^5$ testdeeltjes (die verschillende ionensoorten vertegenwoordigen, primair Fe$^+$ voor computationele efficiëntie, met resultaten gegeneraliseerd naar H$^+$) geïnjecteerd in de turbulente regio. Hun trajecten werden geïntegreerd met behulp van de relativistische bewegingsvergelijking, gekoppeld aan de gesimuleerde snelheid- en magnetische velden.
• Correlatieanalyse: De studie richtte zich op het correleren van de lokale deeltjesversnellingssnelheid ($dE/dt$) met vier specifieke plasmagrootheden die uit de simulatiesnapshots werden geëxtraheerd: magnetische veldsterkte ($|B|$), stroomdichtheid ($|J|$), vorticiteit ($|\nabla \times \mathbf{u}|$) en snelheidsdivergentie ($|\nabla \cdot \mathbf{u}|$). Statistische analyse werd uitgevoerd met behulp van 2D-waarschijnlijkheidsverdelingen (2D-PDF's) en Pearson-correlatiecoëfficiënten, waarbij alleen regio's werden beschouwd waar de versnellingssnelheid 100 eV/s overschreed.

Belangrijkste Resultaten

1. Versnellingsmechanismen: De analyse laat zien dat deeltjesversnelling wordt gedomineerd door turbulente, magnetisch gedreven processen in plaats van klassieke DSA.

   • Magnetisch veld en vorticiteit: Er werden sterke positieve correlaties gevonden tussen de versnellingssnelheid en zowel de magnetische veldsterkte ($|B|$) als de vorticiteit ($|\nabla \times \mathbf{u}|$). De Pearson-correlatiecoëfficiënten varieerden van 0,46 tot 0,55 voor deze variabelen. De regressiehellingen waren superlineair (tot 1,79), wat aangeeft dat sterkere magnetische velden en hogere turbulentie de energiewinst significant verhogen.
   • Snelheidsdivergentie: In contrast hiermee was de correlatie tussen de versnellingssnelheid en de snelheidsdivergentie (een proxy voor schokcompressie) zwak tot niet-bestaand (Pearson-coëfficiënten 0,18–0,33, met negatieve of bijna nul-hellingen). Dit suggereert dat goed gedefinieerde, planaire schokken een beperkte rol spelen bij de waargenomen versnelling.
   • Stroomdichtheid: De stroomdichtheid vertoonde een matige correlatie (0,44–0,51) maar lagere regressiehellingen. De auteurs interpreteren dit als een indirecte tracer van de complexiteit van het magnetische veld in plaats van een directe drijfveer, aangezien resistieve effecten verwaarloosbaar zijn in de ideale MHD-limiet.

2. Energieverdelingen:

   • Deeltjes evolueerden van thermische distributies naar relativistische machtswetten binnen uren tot dagen, waarbij energieën van honderden TeV werden bereikt.
   • De maximale energie en de positie van de "knee" in het energiespectrum zijn afhankelijk van de oppervlaktemagnetische veldsterkte ($B_*$), wat consistent is met het Hillas-confinement criterium.
     • De polytrope index ($\Gamma_{\text{eff}}$) beïnvloedde de spectrale vorm aanzienlijk. Modellen met sterker koeling (lagere $\Gamma_{\text{eff}}$) produceerden hardere spectra en hogere versnellingsefficiënties door verhoogde turbulentie en magnetische veldcomplexiteit. Adiabatische modellen ($\Gamma_{\text{eff}} = 5/3$) vertoonden steilere hellingen en minder efficiënte versnelling.

3. Ruimtelijke Verdeling: Regio's met de hoogste versnelling kwamen niet strikt overeen met de sterkste schokfronten (snelheidsdivergentie), maar eerder met de turbulente downstream-regio's waar magnetische velden worden versterkt en complexer worden.

Betekenis en Claims
Het artikel beweert een definitieve link te bieden tussen lokale plasmacondities en deeltjesversnellings-efficiëntie in CWBs, waarbij verder wordt gegaan dan stationaire aannames. De primaire betekenis ligt in het aantonen dat turbulentie en magnetische veldcomplexiteit de dominante drijvers zijn van deeltjesversnelling in deze systemen, waarbij zij de klassieke diffusieve schokversnelling overschaduwen.

De auteurs concluderen dat:

• CWBs hadronische deeltjes kunnen versnellen tot honderden TeV of zelfs PeV-energieën, waarbij meer dan 1% van de deeltjes het VHE-bereik bereikt, mits er sprake is van een matige magnetisatie ($>1$ G).
• Het versnellingsproces is inherent stochastisch en magnetisch gedreven, waarbij het vertrouwt op de interactie van deeltjes met turbulente structuren (Alfvén-golven, magnetische spiegels, reconnectie) in plaats van eenvoudige schokovergang.
• Deze bevindingen hebben directe implicaties voor het interpreteren van toekomstige hooggevoelige $\gamma$-ray waarnemingen (bijv. door LACT en CTAO), waarbij zij suggereren dat de detecteerbaarheid van CWBs als kosmische stralingsbronnen zwaar afhangt van het niveau van turbulentie en magnetische veldcomplexiteit, die worden gemoduleerd door koelingsefficiëntie en de orbitale fase.

Het werk stelt geen nieuwe experimentele opstellingen voor, maar biedt een theoretisch kader om bestaande en toekomstige observationele gegevens te interpreteren, waarbij wordt benadrukt dat de afwezigheid van sterke gamma-ray detectie in sommige CWBs kan komen door een lagere turbulentie of magnetisatie in plaats van een gebrek aan versnellingsvermogen.