Cyclotron lines in subcritical X-ray pulsars: Monte Carlo simulations reveal the origin of the observed variability

Technische Samenvatting: Cyclotronlijnen in Subkritische Röntgenpulsars

Probleemstelling
Röntgenpulsars (XRPs) vertonen cyclotronresonante verstrooiingskenmerken (CRSF's) in hun spectra, die dienen als directe diagnose voor de sterkte van het magnetische veld van het neutronenster (NS). Waarnemingen tonen aanzienlijke variabiliteit in de centroid-energie ($E_{CRSF}$) en lijnbreedte ($\sigma_{CRSF}$) van de CRSF als functie van de röntgenluminositeit ($L$) en pulsfase. In het subkritische luminositeitregime ($L \lesssim 10^{37}$ erg s$^{-1}$) wordt vaak een positieve correlatie waargenomen tussen $E_{CRSF}$ en $L$, in contrast met de negatieve correlatie die bij superkritische bronnen wordt gezien. Hoewel kwalitatieve modellen suggereren dat dit voortkomt uit de Dopplerverschuiving van fotonen die verstrooien aan in vallend plasma in de accretietrechter, ontbreekt een kwantitatieve, rigoureuze stralingstransportbehandeling van deze variabiliteit. Specifiek vereist de wisselwerking tussen resonante verstrooiing, stralingsdruk-effecten op de stroomsnelheid en het kijkhoek-geometrie gedetailleerde simulaties om de waargenomen trends te verklaren.

Methodologie
De auteurs ontwikkelden een relativistische Monte Carlo (MC) code om gedetailleerde stralingstransportberekeningen uit te voeren binnen de accretietrechter boven het hete vlekje van de NS. Het model hanteert het fysieke raamwerk voorgesteld door Mushtukov et al. (2015b) voor subkritische bronnen:

• Geometrie: Een semi-oneindige, cilindrische accretiekanal gealigneerd met de magnetische as ($B = B\hat{z}$).
• Plasma-profielen: Analytische uitdrukkingen voor elektronendichtheid ($n_e$) en bulk-snelheid ($\beta$) worden afgeleid uit behoud van massa en impuls, waarbij de vertraging van het in vallende plasma door stralingsdruk wordt meegenomen. Het snelheidsprofiel hangt af van de luminositeit ten opzichte van een kritische luminositeit $L^*$.
• Stralingstransport: Zaad-fotonen worden isotroop ingebracht vanuit het hete vlekje met een zwartlichaamspectrum. Fotonenpakketten propageren door het plasma en ondergaan Comptonverstrooiing.
• Verstrooiingsfysica: De code maakt gebruik van benaderde, gepolarisatie-gemiddelde resonante magnetische verstrooiingsdoorsneden (Loudas et al. 2021) die geldig zijn voor $b \ll 1$ en $E \ll m_ec^2$. De verstrooiing wordt behandeld in het ruststelsel van het elektron, met Lorentz-transformaties toegepast om rekening te houden met de bulk-beweging van het plasma.
• Randvoorwaarden: Twee schema's werden getest: "Reflectie" (geometrische reflectie aan de basis, waarbij energie behouden blijft) en "Thermalisatie" (absorptie en heruitzending). De primaire resultaten maken gebruik van het Reflectie-schema.
• Parameters: Simulaties werden uitgevoerd voor een fiduciale NS ($M=1.4 M_\odot$, $R=10^6$ cm, $B \approx 10^{12}$ G) en toegepast op het specifieke geval van de XRP GX 304−1. De accretieluminositeit werd gevarieerd over het bereik $L/L^* \in (0, 1)$.

Belangrijkste Bijdragen en Resultaten
De studie biedt de eerste kwantitatieve MC-verificatie van het model van Mushtukov et al. (2015b), met de volgende resultaten:

1. Spectrale Morfologie: De emergente spectra vertonen een prominente, asymmetrische CRSF vergezeld van een brede blauwe vleugel. De CRSF is systematisch roodverschoven ten opzichte van de lokale klassieke cyclotron-energie ($E_c$). Deze roodverschuiving wordt veroorzaakt door het Dopplereffect tussen het laboratoriumstelsel en het in vallende elektronenstelsel.
2. Luminositeitsafhankelijkheid: Zowel $E_{CRSF}$ als $\sigma_{CRSF}$ tonen een positieve correlatie met luminositeit voor alle kijkhoeken. Naarmate $L$ toeneemt, vertraagt stralingsdruk de in vallende stroming effectiever, waardoor de bulk-snelheid ($\beta$) afneemt. Deze vermindering van de snelheid verlaagt de Doppler-roodverschuiving, waardoor $E_{CRSF}$ dichter bij $E_c$ komt te liggen. Tegelijkertijd verdiept de toegenomen optische diepte bij hogere luminositeiten de lijn en verbreedt de blauwe vleugel.
3. Kijkhoek-afhankelijkheid: De eigenschappen van de CRSF variëren aanzienlijk met de kijkhoek $\theta$ (ten opzichte van de magnetische as).
   • $E_{CRSF}$ neemt af naarmate $\cos \theta$ toeneemt (d.w.z. naarmate het zicht meer "face-on" wordt). Dit komt door de afhankelijkheid van de Dopplerverschuiving van de hoek tussen de fotonimpuls en de stroomsnelheid.
   • $\sigma_{CRSF}$ neemt toe met $\cos \theta$. Face-on-zichten bemonsteren een bredere reeks snelheidsgradiënten langs de gezichtlijn, wat leidt tot grotere Doppler-breedte.
   • Bijgevolg voorspelt het model een anticorrelatie tussen $E_{CRSF}$ en $\sigma_{CRSF}$ tijdens de puls cyclus (fase-opgeloste variabiliteit), terwijl fase-gemiddelde spectra een positieve correlatie met luminositeit vertonen.
4. Toepassing op GX 304−1: Wanneer toegepast op observationele data voor GX 304−1, reproduceert het model succesvol de waargenomen positieve correlaties van $E_{CRSF}$ en $\sigma_{CRSF}$ met luminositeit over bijna een orde van grootte. De beste fit wordt bereikt voor kijkhoeken waarbij de accretietrechter voornamelijk edge-on wordt bekeken ($\cos \theta \approx 0.0-0.1$), consistent met specifieke geometrische configuraties (bijvoorbeeld een bijna-gealigneerde rotator bij hoge inclinatie of een bijna-orthogonale rotator bij lage inclinatie).

Betekenis en Beweringen
Het artikel beweert dat resonante verstrooiing van straling die zich voortplant door het accreterende plasma in de magnetische trechter een natuurlijke, kwantitatieve verklaring biedt voor de waargenomen CRSF-variabiliteit in subkritische XRPs. De auteurs benadrukken dat hun resultaten:

• De hypothese valideren dat het accretiekanal het lijnvormende gebied is in subkritische bronnen.
• Aantonen dat bulk-beweging en de resulterende Dopplerverschuiving de primaire drijvers zijn van de luminositeitsafhankelijke variabiliteit, zonder dat de vorming van een collisionele schok (CS) als primaire lijnvormende locatie vereist is.
• De kritische rol van kijkgeometrie bij het interpreteren van fase-opgeloste data benadrukken, met voorspelling van specifieke anticorrelaties tussen lijnenergie en breedte die verschillen van fase-gemiddelde trends.
• Hun scenario onderscheiden van het CS-model (waarbij lijnvorming plaatsvindt tussen de schok en het oppervlak), met de opmerking dat het CS-scenario de sterke fase-afhankelijke variabiliteit die in hun simulaties wordt waargenomen, niet zou voorspellen.

De auteurs merken beperkingen op, waaronder het gebruik van vereenvoudigde cilindrische geometrie, 1D analytische profielen, en het verwaarlozen van algemene relativistische effecten en polarisatie. Zij betogen echter dat deze vereenvoudigingen de fundamentele conclusies met betrekking tot de oorsprong van de variabiliteit niet wijzigen. Het werk dient als een proof-of-concept voor het gebruik van resonante verstrooiing in de accretietrechter om spectra van subkritische XRPs te modelleren.