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

Technische Zusammenfassung: Zyklotronlinien in subkritischen Röntgenpulsaren

Problemstellung
Röntgenpulsare (XRPs) zeigen in ihren Spektren zyklotrone resonante Streueigenschaften (CRSFs), die als direkte Diagnose für die Stärke des Magnetfelds des Neutronensterns (NS) dienen. Beobachtungen zeigen eine signifikante Variabilität der CRSF-Zentrenenergie ($E_{CRSF}$) und der Linienbreite ($\sigma_{CRSF}$) als Funktion der Röntgenleuchtkraft ($L$) und der Pulsphase. Im subkritischen Leuchtkraftbereich ($L \lesssim 10^{37}$ erg s$^{-1}$) wird häufig eine positive Korrelation zwischen $E_{CRSF}$ und $L$ beobachtet, was im Gegensatz zur negativen Korrelation in superkritischen Quellen steht. Während qualitative Modelle nahelegen, dass dies aus dem Dopplereffekt von Photonen resultiert, die an abfallendem Plasma im Akkretionstrichter streuen, hat es an einer quantitativen, rigorosen Strahlungstransportbehandlung dieser Variabilität gefehlt. Insbesondere erfordert das Zusammenspiel zwischen resonanter Streuung, Strahlungsdruckeffekten auf die Strömungsgeschwindigkeit und der Beobachtungsgeometrie detaillierte Simulationen, um die beobachteten Trends zu erklären.

Methodik
Die Autoren entwickelten einen relativistischen Monte-Carlo-(MC-)Code, um detaillierte Strahlungstransportberechnungen innerhalb des Akkretionstrichters oberhalb des NS-Hotspots durchzuführen. Das Modell übernimmt den von Mushtukov et al. (2015b) für subkritische Quellen vorgeschlagenen physikalischen Rahmen:

• Geometrie: Ein halbunendlicher, zylindrischer Akkretionskanal, der mit der magnetischen Achse ausgerichtet ist ($B = B\hat{z}$).
• Plasmaprofile: Analytische Ausdrücke für die Elektronendichte ($n_e$) und die Bulk-Geschwindigkeit ($\beta$) werden aus der Erhaltung von Masse und Impuls abgeleitet und beinhalten die Verzögerung des abfallenden Plasmas durch Strahlungsdruck. Das Geschwindigkeitsprofil hängt von der Leuchtkraft relativ zu einer kritischen Leuchtkraft $L^*$ ab.
• Strahlungstransport: Seed-Photonen werden isotrop vom Hotspot mit einem Schwarzkörperspektrum injiziert. Photonenpakete propagieren durch das Plasma und unterliegen der Compton-Streuung.
• Streuphysik: Der Code verwendet angenäherte, polarisationsgemittelte resonante magnetische Streuquerschnitte (Loudas et al. 2021), die für $b \ll 1$ und $E \ll m_ec^2$ gültig sind. Die Streuung wird im Elektronenruhesystem behandelt, wobei Lorentz-Transformationen angewendet werden, um die Bulk-Bewegung des Plasmas zu berücksichtigen.
• Randbedingungen: Zwei Schemata wurden getestet: „Reflexion" (geometrische Reflexion an der Basis unter Energieerhaltung) und „Thermalisierung" (Absorption und Wiederemission). Die primären Ergebnisse nutzen das Reflexionsschema.
• Parameter: Simulationen wurden für einen fiduziellen NS ($M=1.4 M_\odot$, $R=10^6$ cm, $B \approx 10^{12}$ G) durchgeführt und auf den spezifischen Fall des XRP GX 304−1 angewendet. Die Akkretionsleuchtkraft wurde im Bereich $L/L^* \in (0, 1)$ variiert.

Hauptbeiträge und Ergebnisse
Die Studie liefert die erste quantitative MC-Verifikation des Modells von Mushtukov et al. (2015b) und ergibt folgende Ergebnisse:

1. Spektrale Morphologie: Die emittierten Spektren zeigen eine ausgeprägte, asymmetrische CRSF, begleitet von einem breiten blauen Flügel. Die CRSF ist systematisch gegenüber der lokalen klassischen Zyklotronenergie ($E_c$) rotverschoben. Diese Rotverschiebung wird durch den Dopplereffekt zwischen dem Laborsystem und dem abfallenden Elektronensystem verursacht.
2. Leuchtkraftabhängigkeit: Sowohl $E_{CRSF}$ als auch $\sigma_{CRSF}$ zeigen für alle Betrachtungswinkel eine positive Korrelation mit der Leuchtkraft. Mit steigendem $L$ verzögert der Strahlungsdruck den abfallenden Fluss effektiver und reduziert die Bulk-Geschwindigkeit ($\beta$). Diese Verringerung der Geschwindigkeit mindert die Doppler-Rotverschiebung und verschiebt $E_{CRSF}$ näher an $E_c$. Gleichzeitig vertieft die erhöhte optische Tiefe bei höheren Leuchtkräften die Linie und verbreitert den blauen Flügel.
3. Betrachtungswinkel-Abhängigkeit: Die CRSF-Eigenschaften variieren signifikant mit dem Betrachtungswinkel $\theta$ (relativ zur magnetischen Achse).
   • $E_{CRSF}$ nimmt ab, wenn $\cos \theta$ zunimmt (d. h., wenn die Sicht „frontaler" wird). Dies liegt an der Winkelabhängigkeit des Dopplereffekts zwischen dem Photonenimpuls und der Strömungsgeschwindigkeit.
   • $\sigma_{CRSF}$ nimmt mit $\cos \theta$ zu. Frontale Ansichten erfassen einen breiteren Bereich von Geschwindigkeitsgradienten entlang der Sichtlinie, was zu einer stärkeren Dopplerverbreiterung führt.
   • Folglich sagt das Modell eine Antikorrelation zwischen $E_{CRSF}$ und $\sigma_{CRSF}$ während des Pulszyklus (phasenaufgelöste Variabilität) voraus, wohingegen phasengemittelte Spektren eine positive Korrelation mit der Leuchtkraft zeigen.
4. Anwendung auf GX 304−1: Bei Anwendung auf Beobachtungsdaten für GX 304−1 reproduziert das Modell erfolgreich die beobachteten positiven Korrelationen von $E_{CRSF}$ und $\sigma_{CRSF}$ mit der Leuchtkraft über fast eine Größenordnung. Die beste Anpassung wird für Betrachtungswinkel erreicht, bei denen der Akkretionstrichter überwiegend von der Kante gesehen wird ($\cos \theta \approx 0.0-0.1$), was mit spezifischen geometrischen Konfigurationen übereinstimmt (z. B. ein nahezu ausgerichteter Rotator bei hoher Inklination oder ein nahezu orthogonaler Rotator bei niedriger Inklination).

Bedeutung und Behauptungen
Die Arbeit behauptet, dass die resonante Streuung von Strahlung, die durch das akkretierende Plasma im magnetischen Trichter propagiert, eine natürliche, quantitative Erklärung für die beobachtete CRSF-Variabilität in subkritischen XRPs liefert. Die Autoren betonen, dass ihre Ergebnisse:

• die Hypothese validieren, dass der Akkretionskanal der linienbildende Bereich in subkritischen Quellen ist;
• zeigen, dass Bulk-Bewegung und der daraus resultierende Dopplereffekt die Haupttreiber der leuchtkraftabhängigen Variabilität sind, ohne dass die Bildung eines kollisionslosen Schocks (CS) als primärer linienbildender Ort erforderlich ist;
• die kritische Rolle der Betrachtungsgeometrie bei der Interpretation phasenaufgelöster Daten hervorheben und spezifische Antikorrelationen zwischen Linienenergie und -breite vorhersagen, die sich von phasengemittelten Trends unterscheiden;
• ihr Szenario vom CS-Modell (bei dem die Linienbildung zwischen Schock und Oberfläche erfolgt) unterscheiden, wobei das CS-Szenario keine starke phasenabhängige Variabilität vorhersagen würde, wie sie in ihren Simulationen beobachtet wird.

Die Autoren weisen auf Einschränkungen hin, darunter die Verwendung einer vereinfachten zylindrischen Geometrie, 1D-analytischer Profile sowie die Vernachlässigung allgemeinrelativistischer Effekte und Polarisation. Sie argumentieren jedoch, dass diese Vereinfachungen die grundlegenden Schlussfolgerungen bezüglich des Ursprungs der Variabilität nicht verändern. Die Arbeit dient als Proof-of-Concept für die Verwendung resonanter Streuung im Akkretionstrichter zur Modellierung subkritischer XRP-Spektren.