Numerical solution of the radiative transfer equation: X-ray spectral formation from cylindrical accretion onto a magnetized neutron star

Abstract

Predicting the emerging X-ray spectra in several astrophysical objects is of great importance, in particular when the observational data are compared with theoretical models. To this aim, we have developed an algorithm solving the radiative transfer equation in the Fokker-Planck approximation when both thermal and bulk Comptonization take place. The algorithm is essentially a relaxation method, where stable solutions are obtained when the system has reached its steady-state equilibrium. We obtained the solution of the radiative transfer equation in the two-dimensional domain defined by the photon energy E and optical depth of the system tau using finite-differences for the partial derivatives, and imposing specific boundary conditions for the solutions. We treated the case of cylindrical accretion onto a magnetized neutron star. We considered a blackbody seed spectrum of photons with exponential distribution across the accretion column and for an accretion where the velocity reaches its maximum at the stellar surface and at the top of the accretion column, respectively. In both cases higher values of the electron temperature and of the optical depth tau produce flatter and harder spectra. Other parameters contributing to the spectral formation are the steepness of the vertical velocity profile, the albedo at the star surface, and the radius of the accretion column. The latter parameter modifies the emerging spectra in a specular way for the two assumed accretion profiles. The algorithm has been implemented in the XSPEC package for X-ray spectral fitting and is specifically dedicated to the physical framework of accretion at the polar cap of a neutron star with a high magnetic field (> 10^{12} G), which is expected to be typical of accreting systems such as X-ray pulsars and supergiant fast X-ray transients.

13 pages, 20 figures, accepted for publication in A&A

References in corpus (3)

• Thermal and Bulk Comptonization in Accretion-Powered X-Ray Pulsars
• A new Comptonization model for low-magnetized accreting neutron stars in low mass X-ray binaries
• On the stability of the thermal Comptonization index in neutron star low-mass X-ray binaries in their different spectral states

Cited by in corpus (40)

• Cyclotron lines in highly magnetized neutron stars
• Towards a unified view of inhomogeneous stellar winds in isolated supergiant stars and supergiant high mass X-ray binaries
• Outburst of GX 304-1 monitored with INTEGRAL: positive correlation between the cyclotron line energy and flux
• A new model for the X-ray continuum of the magnetized accreting pulsars
• Luminosity-dependent changes of the cyclotron line energy and spectral hardness in Cep X-4
• Revisiting the archetypical wind accretor Vela X-1 in depth -- A case study of a well-known X-ray binary and the limits of our knowledge
• RX J0440.9+4431: a persistent Be/X-ray binary in outburst
• The NuSTAR X-ray Spectrum of Hercules X-1: A Radiation-Dominated Radiative Shock
• Multiple cyclotron line-forming regions in GX 301-2
• Observations of a GX 301-2 Apastron Flare with the X-Calibur Hard X-Ray Polarimeter Supported by NICER, the Swift XRT and BAT, and Fermi GBM
• Seven Years with the Swift Supergiant Fast X-ray Transients Project
• Dynamical and Radiative Properties of X-ray Pulsar Accretion Columns: Phase-Averaged Spectra
• Variations in the Cyclotron Resonant Scattering Features during 2011 outburst of 4U 0115+63
• Accretion regimes in the X-ray pulsar 4U 1901+03
• The Galactic LMXB Population and the Galactic Centre Region
• Looking at A 0535+26 at low luminosities with NuSTAR
• Revisiting the spectral and timing properties of 4U 1909+07 with NuSTAR and Astrosat
• Hard X-ray emission of Sco X-1
• Swift observations of two supergiant fast X-ray transient prototypes in outburst
• Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3
• Glancing through the accretion column of EXO 2030+375
• NuSTAR rules out a cyclotron line in the accreting magnetar candidate 4U2206+54
• Evidence for different accretion regimes in GRO J1008-57
• Three-dimensional modelling of accretion columns: spatial asymmetry and self-consistent simulations
• The Swift Supergiant Fast X-ray Transients Project: a review, new results, and future perspectives
• Spectral and temporal properties of the supergiant fast X-ray transient IGR J18483-0311 observed by INTEGRAL
• NuSTAR view of Be/X-ray binary pulsar 2S 1417-624 during 2018 giant outburst
• The giant outburst of 4U 0115+634 in 2011 with Suzaku and RXTE
• Broad-band X-ray spectra of anomalous X-ray pulsars and soft $\gamma$-ray repeaters: pulsars in a weak-accretion regime ?
• Suzaku observations of the 2013 outburst of KS 1947+300
• The 100-month Swift catalogue of supergiant fast X-ray transients II. SFXT diagnostics from outburst properties
• Possible detection of a new cyclotron feature in 4U 1700-37
• X-Ray Views of Galactic Accreting Pulsars in High-Mass X-Ray Binaries
• Disk precession to explain the super-orbital modulation of LMC X-4: results from the Swift monitoring campaign
• The unaltered pulsar: GRO J1750-27, a super-critical X-ray neutron star that does not blink an eye
• Timing and spectral studies of Cen X-3 in multiple luminosity states using $\text{AstroSat}$
• Observational constraints on the magnetic field of the bright transient Be/X-ray pulsar SXP 4.78
• A curious case of the accretion-powered X-ray pulsar GX 1+4
• Spectral evolution of X-ray pulsar 4U 1901+03 during the 2019 outburst based on Insight-HXMT and NuSTAR observations
• Variability study of classical supergiant X-ray binary 4U 1907+09 using NuSTAR