Classical, large scale 3D MHD simulations of interacting pulsar wind nebulae

Probleemstelling en Context

Pulsar-windnevels (PWN's) zijn extreme omgevingen die ontstaan wanneer de relativistische wind van een snel roterende, gemagnetiseerde neutronenster (pulsar) interacteert met zijn omgeving. De morfologie en emissie-eigenschappen van deze nevels worden sterk beïnvloed door het omringende medium, dat bestaat uit supernova-uitwerpselen en het circumstellaire medium (CSM) van de voorouderster.

Eerdere studies (zoals Meyer et al. 2024a) hebben aangetoond dat bij statische massieve sterren het CSM, gevormd door sterwinden, een langwerpige holte creëert die uitgelijnd is met het lokale interstellair magnetisch veld (ISM). Als de rotatieas van de pulsar samenvalt met dit magnetische veld, ontstaat er een PWN met een karakteristieke vorm (vaak rechthoekig of diamantvormig). Echter, de 2D-modellen uit eerdere werken dwongen de auteurs om aan te nemen dat de richting van het ISM-magnetische veld en de rotatieas van de pulsar colineair (in één lijn) waren.

Het centrale probleem dat dit artikel adresseert, is het ontbreken van kennis over hoe de hoek tussen de rotatieas van de pulsar en de richting van het lokale ISM-magnetische veld de vorming, structuur en niet-thermische radio-emissie van de PWN beïnvloedt in een volledig 3D-context.

Methodologie

De auteurs voeren uitgebreide 3D-magnetohydrodynamische (MHD) simulaties uit met de PLUTO-code om de evolutie van een PWN binnen een supernovarest (SNR) na te bootsen. De aanpak bestaat uit een multi-stap workflow:

1. Modellering van het Circumstellaire Medium (CSM):

   • Een statische $35 M_\odot$-ster wordt gemodelleerd van de nul-leeftijd op de hoofdreeks tot de pre-supernova-fase.
   • Dit gebeurt in 2,5D (cylindrische coördinaten) om de interactie van de sterwind met het ISM te simuleren, waarbij rekening wordt gehouden met sterrenwindbellen van verschillende fasen (Rode Reus, Wolf-Rayet).
   • Het ISM wordt gemodelleerd als een homogeen medium (warm fase, $T=8000$ K, $n=0.79$ cm$^{-3}$) met een uniform magnetisch veld van $7 \mu$G.

2. Supernova-explosie en 1D-fase:

   • Een sferisch symmetrische kerninstortings-supernova-explosie ($10^{51}$ erg) wordt in 1D gesimuleerd totdat de schokgolf een straal van 2 pc bereikt.

3. 3D Globale Simulatie:

   • De 2,5D CSM-structuur wordt omgezet naar een 3D-structuur (via rotatie om de symmetrie-as).
   • De 1D-oplossing wordt hierop gemapt als initiële condities voor een 3D Cartesische simulatie (grid van $2048^3$ zones).
   • De pulsar wordt ingebracht na ~20 jaar met een initiële spinperiode van 0,3 s en een spin-down vermogen van $10^{38}$ erg/s.
   • Variabele Parameter: De hoek $\theta_{mag}$ tussen de rotatieas van de pulsar (vastgesteld langs de Oz-as) en de richting van het ISM-magnetisch veld (dat de vorm van de CSM-holte bepaalt) wordt gevarieerd: $0^\circ, 30^\circ, 45^\circ$ en $90^\circ$.

4. Stralingstransport:

   • Er worden niet-thermische radiokaartjes (synchrotronstraling bij 1,4 GHz) gegenereerd met RADMC-3D om de simulaties direct vergelijkbaar te maken met astronomische waarnemingen.

Belangrijkste Bijdragen

• Eerste 3D MHD-studie: Dit is een van de eerste studies die de volledige evolutie van een PWN in een statische SNR in 3D simuleert, waarbij expliciet rekening wordt gehouden met de misalignement tussen de pulsar-rotatieas en het ISM-magnetisch veld.
• Mechanisme voor vervorming: Het artikel identificeert een nieuw mechanisme voor de vervorming van polaire uitlopers van PWN's. In plaats van alleen instabiliteiten (zoals kink-instabiliteiten) als oorzaak te zien, tonen de auteurs aan dat de geometrie van de door de sterwind uitgeholde holte de polaire uitlopers fysiek kan buigen en kan leiden tot complexe structuren.
• Synthetische Kaarten: Het leveren van synthetische radiokaartjes voor verschillende projectiehoeken en hoeken van misalignement, wat cruciaal is voor de interpretatie van waarnemingen van complexe PWN's.

Resultaten

1. Invloed van de Hoek ($\theta_{mag}$) op de Morfologie:

   • $\theta_{mag} = 0^\circ$ (Gealigneerd): De PWN ontwikkelt een equatoriale schijf en polaire uitlopers die uitgelijnd zijn met het magnetisch veld. De vorm is grotendeels symmetrisch en kan rechthoekig lijken in projectie.
   • $\theta_{mag} = 30^\circ - 45^\circ$ (Gedeeltelijk misaligneerd): De polaire uitlopers buigen in tegenovergestelde richtingen door de interactie met de wanden van de CSM-holte. Dit leidt tot complexe, onregelmatige vormen (oblong, "wobbly" schijven) en een gebroken bipolariteit. De schijf groeit in een vlak dat een hoek maakt met het equatoriale vlak van de pulsar.
   • $\theta_{mag} = 90^\circ$ (Perpendiculair): De polaire uitlopers botsen tegen de wanden van de holte en worden dikker. De schijf wordt zeer uitgebreid en neemt een rechthoekige vorm aan met uitsteeksels langs de y-as. De diamantvormige structuur die in eerdere 2D-studies werd gezien, verschijnt opnieuw maar met dikkere componenten.

2. Tijdsontwikkeling:

   • Naarmate de SNR veroudert (40 tot 80 kyr), worden de structuren complexer. De interactie tussen de gereflecteerde schokgolven van de supernova en de pulsarwind leidt tot het loslaten van klonten materiaal en een toenemende verhouding tussen de grootte van de schijf en de polaire uitlopers.

3. Magnetische Veldstructuur:

   • Bij misalignement verliest het magnetisch veld zijn symmetrie. De veldlijnen reconecteren aan het einde van de holte, wat leidt tot een complexe, verweven structuur die sterk afwijkt van de geordende toroidale velden bij gealigneerde gevallen.

4. Menging van Materialen:

   • Een verrassend resultaat is dat de hoek $\theta_{mag}$ geen significant effect heeft op de menging van materialen (sterwind, supernova-uitwerpselen, pulsarwind) op oude tijdschalen (50-80 kyr).
   • De mengingsefficiëntie nadert 100% voor hoofdreeks- en Wolf-Rayet-materiaal, terwijl rode reus-wind en supernova-uitwerpselen een menging van ~20% behouden. Dit suggereert dat de interne chemische samenstelling op lange termijn niet sterk afhankelijk is van de initiële hoek tussen de magnetische velden.

5. Radio-emissie:

   • De synthetische kaarten tonen een enorme variatie aan vormen: van rechthoekig tot circulair en onregelmatig, afhankelijk van de projectiehoek en de waarde van $\theta_{mag}$. Dit verklaart de diversiteit aan waargenomen PWN-vormen in het heelal.

Significantie en Conclusie

Dit onderzoek bevestigt dat de hoek tussen de rotatieas van de pulsar en het lokale ISM-magnetisch veld een bepalende parameter is voor de morfologie en de waarnemingskenmerken van PWN's in statische systemen. Het biedt een alternatieve verklaring voor de vaak onregelmatige en gebogen polaire uitlopers die in waarnemingen worden gezien, zonder noodzaak te hebben tot extreme instabiliteiten of hoge snelheden van de pulsar zelf.

Hoewel de vorming van de nevel sterk wordt beïnvloed door deze hoek, blijkt de chemische menging op lange termijn robuust te zijn tegenover deze variatie. De studie onderstreept het belang van 3D-simulaties om de complexe interacties tussen sterwinden, supernova-uitwerpselen en magnetische velden correct te begrijpen, en biedt een raamwerk voor het interpreteren van toekomstige radio-observaties van pulsar-windnevels.