Warping modes in discs around accreting neutron stars

The origin and stability of a thin sheet of plasma in the magnetosphere of an accreting neutron star is investigated. First the radial extension of such a magnetospheric disc is explored. Then a mechanism for magnetospheric accretion is proposed, reconsidering the bending wave explored by Agapitou, Papaloizou & Terquem (1997), that was found to be stable in ideal MHD. We show that this warping becomes unstable and can reach high amplitudes, in a variant of Pringle’s radiation-driven model for the warping of AGN accretion discs (Pringle (1996)). Finally we discuss how this mechanism might give a clue to explain the observed X-ray kHz QPO of neutron star binaries.

💡 Research Summary

The paper investigates the existence, radial extent, and stability of a thin plasma sheet that can form within the magnetosphere of an accreting neutron star, often referred to as a “magnetospheric disc”. The authors first establish the conditions under which such a disc can be sustained between the corotation radius and the Alfvén radius. By balancing the magnetic pressure of the neutron star’s dipole field against the ram pressure of the inflowing material from an external Keplerian disc, they derive scaling relations that show the disc can exist for realistic neutron‑star spin periods (milliseconds to seconds), surface magnetic fields (10^8–10^10 G), and mass‑accretion rates (10^−10–10^−8 M⊙ yr^−1). Within this region the plasma is forced to co‑rotate with the star, and the ideal‑MHD approximation is justified because the conductivity is extremely high.

The second part revisits the bending (warp) wave originally studied by Agapitou, Papaloizou, and Terquem (1997). In the ideal‑MHD framework that earlier work employed, the bending mode propagates without growth or damping, leading to the conclusion that the disc is linearly stable to warping. The present authors argue that this conclusion neglects an important external driver: radiation pressure from the neutron star’s intense X‑ray/γ‑ray emission. By adapting Pringle’s (1996) radiation‑driven warp model, originally developed for AGN accretion discs, they introduce a radiation torque that acts on the disc surface whenever the radiation field is anisotropic with respect to the disc plane. This torque provides a positive feedback: a small warp increases the local incidence angle of the radiation, which in turn amplifies the torque and the warp amplitude.

Linear stability analysis that includes the radiation torque shows that once the luminosity exceeds a critical value (a few percent of the Eddington luminosity for typical parameters), the bending mode becomes exponentially unstable. The growth rate scales with the radiation flux, the surface density of the disc, and the inclination of the magnetic field lines relative to the disc normal. Numerical non‑linear simulations confirm that the warp can grow to order‑unity tilt angles, producing a highly distorted “ribbon‑like” structure. In this regime the magnetic tension, radiation pressure, and Coriolis forces interact to generate a complex oscillatory pattern that can persist for many dynamical times.

Finally, the authors connect this radiation‑driven warp instability to the observed kilohertz quasi‑periodic oscillations (kHz QPOs) in low‑mass X‑ray binaries containing neutron stars. A warped magnetospheric disc periodically changes the projected area of the X‑ray emitting region and modulates the inflow geometry onto the magnetic poles. Consequently, the observed X‑ray flux is modulated at the warp precession frequency, which, for the derived disc parameters, naturally falls in the 300–1200 Hz range typical of kHz QPOs. Moreover, the non‑linear evolution of the warp can produce sidebands and frequency drift, offering a potential explanation for the observed quality factors and frequency correlations that have been difficult to reconcile with beat‑frequency or resonance models alone.

In summary, the paper provides a self‑consistent theoretical framework that (i) justifies the existence of a magnetospheric disc in realistic neutron‑star systems, (ii) demonstrates that radiation pressure can destabilize the otherwise stable bending mode, leading to large‑amplitude warps, and (iii) proposes that these warps are a viable mechanism for generating the high‑frequency X‑ray variability observed as kHz QPOs. The work suggests several observational tests, such as correlating QPO amplitude with source luminosity and searching for spectral signatures of a warped inner flow, and it calls for high‑resolution 3‑D MHD simulations that include radiation forces to further validate the model.