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

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.

💡 Research Summary

This paper addresses the long‑standing problem of why cyclotron resonant scattering features (CRSFs) in subcritical X‑ray pulsars (those with luminosities below the critical value L* ≈ 10³⁷ erg s⁻¹) show a positive correlation between the line centroid energy (E_CRSF) and the X‑ray luminosity, as well as a similar correlation for the line width (σ_CRSF). The authors build upon the analytical framework proposed by Mushtukov et al. (2015b), which posits that radiation pressure from the hot‑spot emission decelerates the infalling plasma above the neutron‑star surface, thereby reducing the bulk Doppler red‑shift of photons that undergo resonant scattering in the magnetic funnel.

To test this idea quantitatively, the authors develop a relativistic Monte Carlo radiative‑transfer code that follows photons emitted isotropically from a black‑body hot spot at the base of a cylindrical accretion funnel. The funnel is assumed to have a uniform magnetic field aligned with the cylinder axis. The plasma density and velocity profiles are taken from analytical expressions: the velocity β(h) follows a non‑linear function of height h that depends on the free‑fall speed, the terminal surface speed β₀, and the funnel radius d; β₀ itself scales as β_ff √(1 − L/L*), linking it directly to the luminosity. Mass conservation yields a density profile n_e(h) ∝ L^{3/5} B^{1/2} β(h)^{-1}.

The code includes both non‑resonant Thomson scattering and resonant magnetic scattering, the latter using an approximate cross‑section up to ~10⁶ σ_T near the cyclotron resonance. Photons are tracked through multiple scatterings, allowing the authors to capture the combined effects of Doppler shifts, resonant trapping, and angular redistribution. Simulations are performed for a range of luminosities (10³⁶–10³⁸ erg s⁻¹) and viewing angles θ (0–π/2).

Key results are:

1. The emergent spectra display an asymmetric CRSF with a pronounced blue wing. The blue wing originates from photons that are Doppler‑blue‑shifted in the upper part of the funnel where the plasma flow is directed away from the observer.

2. For every viewing angle, the CRSF centroid energy increases with luminosity. The increase is modest (≈5–10 % when L rises by a factor of ten) but systematic, reflecting the reduced Doppler red‑shift as radiation pressure slows the flow.

3. The line width also grows with luminosity because higher luminosities increase the resonant optical depth, leading to more scatterings before photons escape.

4. Viewing‑angle dependence is strong. At fixed luminosity, larger cos θ (i.e., more face‑on viewing) yields lower E_CRSF and larger σ_CRSF, while edge‑on views (θ ≈ 90°) produce higher centroid energies and narrower lines. This is a direct consequence of the projection of the bulk flow velocity onto the line of sight.

5. When the model is applied to the well‑studied pulsar GX 304‑1, which spans nearly an order of magnitude in luminosity, the simulated E_CRSF–L and σ_CRSF–L trends match the observations if the accretion funnel is viewed predominantly edge‑on. The model also predicts an anticorrelation between centroid energy and width over the pulse phase, consistent with phase‑resolved data.

The authors compare their findings with alternative explanations. Collisionless shock models predict a positive E_CRSF–L correlation because the shock height decreases with increasing accretion rate, but they require a detailed kinetic treatment of shock formation that is not yet available. Radiation‑mediated shock models, appropriate for super‑critical sources, predict a negative correlation, opposite to what is observed in subcritical systems. The present work demonstrates that a simple radiation‑pressure‑induced deceleration of the flow, combined with resonant scattering, suffices to reproduce the observed phenomenology without invoking shocks.

Limitations are acknowledged: the plasma temperature and electron distribution are treated as isothermal and non‑relativistic; the magnetic field is assumed uniform, ignoring possible multipolar components; and the funnel geometry is simplified to a cylinder. Future work should incorporate full RMHD simulations to obtain self‑consistent temperature and velocity fields, explore more realistic magnetic topologies, and include the effects of electron‑positron pair production at higher luminosities.

In conclusion, the paper provides the first quantitative Monte Carlo verification that Doppler effects in a radiation‑decelerated accretion funnel naturally generate the observed positive correlations of CRSF centroid energy and width with luminosity in subcritical X‑ray pulsars. The model offers a unified framework for interpreting both luminosity‑dependent and pulse‑phase‑dependent CRSF variability, and it sets the stage for more sophisticated simulations that can further constrain neutron‑star magnetic field geometry and accretion physics.