Numerical investigation of particle acceleration at interplanetary shocks: diffusive and superdiffusive scenarios

Probleemstelling

Energetische deeltjes zijn overal aanwezig in ruimtelijke en astrofysische plasma's, en interplanetaire schokgolven (IP-schokken) worden beschouwd als een van de belangrijkste versnellingsmechanismen in de heliosfeer. Hoewel waarnemingen vaak een piek in de deeltjesflux bij de schok tonen, wat wijst op lokale versnelling, is het transport van deze deeltjes stroomopwaarts en stroomafwaarts complex.
Traditionele modellen, zoals Diffusive Shock Acceleration (DSA), gaan uit van normale diffusie (Brownse beweging). Echter, waarnemingen tonen vaak stroomopwaartse fluxen die niet exponentieel afnemen (zoals bij DSA voorspeld), maar eerder een machtsverval vertonen. Dit suggereert "anomal transport" of superdiffusie, waarbij deeltjes sneller bewegen dan bij normale diffusie ($\langle \Delta r^2 \rangle \propto t^\alpha$ met $\alpha > 1$). Een cruciale vraag is hoe deze superdiffusie de versnellingsprocessen beïnvloedt en of dit leidt tot energie-spectra die beter overeenkomen met waarnemingen dan het standaard DSA-model.

Methodologie

De auteurs hebben een nieuw numeriek test-deeltjesmodel ontwikkeld, voortbouwend op eerdere werken (Prete et al., 2019, 2021), om de versnelling van deeltjes bij een oneindig vlakke schok te simuleren.

• Simulatieopzet: Het model simuleert een 1D doos van $-10^{10}$ m tot $10^{10}$ m met een schok bij $x=0$. Deeltjes worden geïnjecteerd met een initiële energie van 70 keV (zaaddeeltjes).
• Transportmechanismen: Er worden twee scenario's vergeleken:
  1. Normale Diffusie (DSA): Deeltjes ondergaan een random walk met een vaste verstrooiingstijd ($\tau_s$).
  2. Superdiffusie (SSA): Deeltjes volgen een Lévy-walk. Hierbij zijn de vrije padlengtes en de verstrooiingstijden gekoppeld via een machtsverdeling met exponent $\mu$. Dit staat lange "sprongen" toe, wat leidt tot superdiffusief gedrag ($\alpha = 4 - \mu$ voor $2 < \mu < 3$).
• Versnellingsmechanisme: In tegenstelling tot eerdere versies van dit model, krijgen deeltjes nu impulswinst bij elke overgang door de schok (eerste-orde Fermi-versnelling). De impulsverandering wordt berekend via Galilei-transformaties tussen de stroomopwaartse en stroomafwaartse referentiekaders.
• Data-gedreven validatie: De simulaties worden vergeleken met waarnemingen van de ACE-ruimtesonde tijdens een schokovergang op 14 december 2006. De parameters voor de simulaties ($\tau_s$ voor diffusie; $\mu$ en $\tau_0$ voor superdiffusie) worden afgeleid uit de waarnemingsdata van deze specifieke gebeurtenis.

Belangrijkste Bijdragen

1. Integratie van versnelling in Lévy-walk modellen: Voor het eerst wordt een test-deeltjescode voor superdiffusie uitgebreid met een Fermi-achtig versnellingsmechanisme bij schokovergangen, wat een directe vergelijking mogelijk maakt met waargenomen energie-spectra.
2. Data-gedreven simulatie: De auteurs hebben een methode ontwikkeld om de modelparameters ($\mu$ en $\tau_0$) direct af te leiden uit de vorm van de deeltjesflux in de waarnemingen, waardoor de simulaties specifiek zijn afgestemd op een echte interplanetaire schok.
3. Vergelijking DSA vs. SSA: Het biedt een robuuste numerieke vergelijking tussen standaard diffusie en superdiffusie in termen van versnellingstijden en resulterende energie-spectra.

Resultaten

• Deeltjesdichtheid en Trajecten:
  • Bij normale diffusie blijven deeltjes in de vroege fase dicht bij de schok, maar worden ze later gedomineerd door advektie stroomafwaarts, met een lage kans om terug te keren.
  • Bij superdiffusie blijven deeltjes langer in de buurt van de schok en hebben ze door de lange sprongen (Lévy-statistiek) een veel hogere kans om terug te keren stroomopwaarts. Dit versnelt het versnellingsproces aanzienlijk.
• Energie-spectra:
  • De numerieke spectra komen uitstekend overeen met theoretische voorspellingen.
  • Superdiffusie leidt tot hardere spectra (minder steil) dan DSA, omdat de versnellingsindex $\gamma$ afhankelijk is van zowel de compressieratio $r$ als de superdiffusie-exponent $\mu$.
  • De simulaties tonen aan dat superdiffusie de tijd verkort die nodig is om hoge energieën te bereiken.
• Validatie met ACE-data (14/12/2006):
  • De fluxen stroomopwaarts van de schok vertonen een machtsverval (power-law decay) in plaats van een exponentieel verval. Een fit met een machtsverval (superdiffusie) levert een veel betere overeenkomst op ($\tilde{\chi}^2$ is veel lager) dan een exponentiële fit (DSA).
  • De afgeleide parameters zijn $\mu \approx 2.66$ en $\tau_0 \approx 4$ s.
  • De simulaties met deze superdiffusieve parameters reproduceren de waargenomen fluxprofielen zowel stroomopwaarts als stroomafwaarts zeer nauwkeurig.
  • Het berekende differentieel energie-spectrum stroomafwaarts komt overeen met de ACE-metingen en bevestigt de voorspelde steilte voor superdiffusie.

Betekenis en Conclusie

De studie concludeert dat anomal transport (superdiffusie) een cruciale rol speelt bij de versnelling van energetische deeltjes bij interplanetaire schokken.

• Superdiffusie verklaart niet alleen de waargenomen machtsverval-profielen stroomopwaarts, maar versnelt ook het versnellingsproces, waardoor deeltjes sneller hoge energieën bereiken dan voorspeld door standaard DSA.
• Het model, hoewel vereenvoudigd (bijv. geen zelfgegenereerde turbulentie of 3D-structuur), slaagt erin de belangrijkste observatiekenmerken van IP-schokken te reproduceren wanneer superdiffusie als hoofdtransportmechanisme wordt aangenomen.
• Dit onderstreept het belang van het gebruik van Lévy-walk modellen in plaats van alleen Gaussische diffusie voor het begrijpen van deeltjesversnelling in de heliosfeer en andere astrofysische omgevingen.

De auteurs pleiten voor verdere analyse van een groot aantal schokgebeurtenissen (bijv. met Solar Orbiter-data) om te bepalen onder welke specifieke omstandigheden (Mach-getallen, geometrie, turbulentie) superdiffusie optreedt.