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

TL;DR

The paper tackles CRSF variability in subcritical X-ray pulsars by modeling resonant scattering in the magnetized accretion funnel above the hotspot with a relativistic Monte Carlo radiative-transfer code. It uses analytic, radiation-pressure–driven velocity and density profiles and approximate resonant cross-sections to compute angle-resolved spectra as the accretion luminosity $L$ varies, predicting a redshifted CRSF with a broad blue wing and positive $E_{\rm CRSF}$ and $\sigma_{\rm CRSF}$–$L$ correlations for all viewing angles. Phase-resolved analysis reveals an anticorrelation between the CRSF centroid energy and width within a given pulse phase due to viewing-angle dependence of Doppler shifts. When applied to GX 304-$p$1, the model reproduces observed CRSF variability over nearly an order of magnitude in $L$, favoring edge-on funnel geometries and supporting the accretion funnel as the line-forming region; the work also discusses boundary-condition effects, model limitations, and prospects for including RMHD, polarization, and general-relativistic effects in future studies.

Abstract

Observed cyclotron resonant scattering features (CRSFs) in X-ray pulsars (XRPs) exhibit strong variability. In the subcritical luminosity regime, the centroid energy ($E_{CRSF}$) and line width ($σ_{CRSF}$) often show positive correlations with the X-ray luminosity. We investigate the physical origin of the observed variability quantitatively, focusing on the effects of resonant scattering and Doppler shift induced by the plasma flow in the accretion funnel. We developed a relativistic Monte Carlo code to perform detailed radiative transfer calculations in the accretion funnel above the hotspot and derive angle-dependent spectra. Analytical plasma density and velocity profiles were adopted to account for the effects of radiation pressure on the flow. Approximate resonant scattering cross-sections were employed. We varied the accretion luminosity to explore the resulting variability of the CRSF properties. The emergent spectra exhibit a prominent, asymmetric CRSF accompanied by a broad blue wing. The CRSF is systematically redshifted relative to the classical cyclotron energy, with the magnitude of the redshift decreasing at higher luminosities and for larger viewing angles $θ$. Both $E_{CRSF}$ and $σ_{CRSF}$ correlate positively with luminosity for all viewing angles. Their absolute values, however, depend strongly on the viewing angle, indicating substantial variability over the pulse cycle and sensitivity to the system geometry. At fixed luminosity, $E_{CRSF}$ ($σ_{CRSF}$) decreases (increases) with increasing $\cosθ$. Consequently, phase-resolved observations are expected to reveal an anticorrelation between the CRSF centroid energy and width. When applied to the XRP GX 304$-$1, the model reproduces the observed CRSF variability over nearly an order of magnitude in luminosity for geometries in which the accretion funnel is predominantly viewed edge-on.

Figures (11)

• Figure 1: Illustration of the MC code's geometry, coordinates, and calculation of the line-of-sight optical depth. A photon packet located at $\vec{r}_0$ is moving in direction $\hat{n}$. The accumulated optical depth is computed by dividing the line of sight into a sequence of points $\vec{r}_i$ using an adaptive step-size algorithm. A target optical depth $\tau_{\rm rnd}$ is then sampled, and the position of the next scattering event is determined via linear interpolation between the discretized points. We assume a uniform magnetic field because the height of the line-forming region is much smaller than the NS radius.
• Figure 2: Representative emergent spectrum for an XRP with a fiducial luminosity $L=0.6L^*$. The solid black line represents the seed blackbody spectrum with temperature $k_{\mathrm{B}} T=5.9~\mathrm{keV}$. The solid red line shows the total normalized emergent spectrum. The dashed blue and green lines correspond to the angle-selected spectra observed at polar angles in the ranges $(0,\pi/2)$ and $(\pi/2,\pi)$, respectively. The vertical dashed gray line indicates the cyclotron energy at $E_{\mathrm{c}}=25.55~\mathrm{keV}$. The spectrum features a redshifted, but prominent CRSF around $E\sim 20\,\mathrm{keV}$ followed by a bump.
• Figure 3: Angle-resolved spectra for a typical XRP source with luminosity $L=0.6L^*$. The solid black line corresponds to the angle-averaged spectrum in the range $(0, \pi/2)$. The dashed colored lines show spectra extracted in five angular bins within this range. Each spectrum is normalized to its angular flux contribution. The vertical gray line marks the local cyclotron energy. With increasing $\cos\theta$, the CRSF centroid energy shifts to lower values and the blue wing becomes more pronounced.
• Figure 4: Angle-averaged emergent spectra for luminosities spanning the range $L/L^*\in (0,1)$. The color gradient from dark blue to yellow corresponds to increasing luminosity. The inset highlights the spectral region around the CRSF. The vertical gray dashed line denotes the local cyclotron energy. The positive correlation between $L$ and the CRSF energy $E_{\mathrm{CRSF}}$ is evident.
• Figure 5: Dependence of the CRSF centroid energy $E_{\mathrm{CRSF}}$ (panel a) and width $\sigma_{\mathrm{CRSF}}$ (panel b) on luminosity and viewing angle $\theta$ (measured with respect to the magnetic axis). Results are shown for five angular bins spanning the range $\theta \in (0,\pi/2)$, with $E_{\mathrm{c}} = 25.55\,\mathrm{keV}$. Solid black lines denote quantities extracted from angle-averaged spectra. Both $E_{\rm CRSF}$ and $\sigma_{\rm CRSF}$ exhibit a positive correlation with luminosity for all viewing angles. At fixed luminosity, $E_{\mathrm{CRSF}}$ ($\sigma_{\mathrm{CRSF}}$) decreases (increases) toward smaller viewing angles, indicating substantial variability over the pulse cycle.