Abstract:Recent observations of X-ray pulsars at low luminosities allow, for the first time, to compare theoretical models for the emission from highly magnetized neutron star atmospheres at low mass accretion rates ($\dot{M} \lesssim 10^{15}$ g s$^{-1}$) with the broadband X-ray data. The purpose of this paper is to investigate the spectral formation in the neutron star atmosphere at low $\dot{M}$ and to conduct a parameter study of physical properties of the emitting region. We obtain the structure of the static atmosphere, assuming that Coulomb collisions are the dominant deceleration process. The upper part of the atmosphere is strongly heated by the braking plasma, reaching temperatures of 30-40 keV, while its denser isothermal interior is much cooler (~2 keV). We numerically solve the polarized radiative transfer in the atmosphere with magnetic Compton scattering, free-free processes, and non-thermal cyclotron emission due to possible collisional excitations of electrons. The strongly polarized emitted spectrum has a double-hump shape that is observed in low-luminosity X-ray pulsars. A low-energy "thermal" component is dominated by extraordinary photons that can leave the atmosphere from deeper layers due to their long mean free path at soft energies. We find that a high-energy component is formed due to resonant Comptonization in the heated non-isothermal part of the atmosphere even in the absence of collisional excitations. The latter, however, affect the ratio of the two components. A strong cyclotron line originates from the optically thin, uppermost zone. A fit of the model to NuSTAR and Swift/XRT observations of GX 304-1 provides an accurate description of the data with reasonable parameters. The model can thus reproduce the characteristic double-hump spectrum observed in low-luminosity X-ray pulsars and provides insights into spectral formation.

A broadband self-consistent modelling of the X-ray spectrum of 4U 1626--67

Broadband X-Ray Spectroscopy of the Pulsar Wind Nebula in HESS J1640-465

A BeppoSAX observation of the massive X-ray binary 4U1700-37

Spectral Formation in X-Ray Pulsars: Bulk Comptonization in the Accretion Shock

Looking into the Theory of Pulsar Accretion: Cen X-3 and XTE J1946+274

Spectrum formation in X-ray pulsars at very low mass accretion rate: Monte-Carlo approach

Probing the emission mechanism and nature of the pulsating compact object in the X-ray binary SAX J1324.4-6200

AstroSat and Insight-HXMT observations of the long-period X-ray pulsar 4U 2206+54

X-ray emission from magnetized neutron star atmospheres at low mass accretion rates. I. Phase-averaged spectrum

A Monte Carlo Simulation of the Broad Band X-ray Emission of the Accreting Millisecond X-ray pulsar MAXI J1816-195

X-Ray spectral and temporal properties of LMXB 4U 1608-52- observed with AstroSat and NICER

Spectral Formation in X-Ray Pulsar Accretion Columns

Study of the X-ray Pulsar XTE J1946+274 with NuSTAR

Theoretical Analysis of the RX J0209.6-7427 X-ray Spectrum During Giant Outburst

Dynamical and Radiative Properties of X-Ray Pulsar Accretion Columns: Phase-averaged Spectra

Broadband X-ray timing and spectral characteristics of the accretion-powered millisecond X-ray pulsar MAXI J1816$-$195

Peculiarities of super-Eddington flares from the X-ray pulsar LMC X-4 with NuSTAR

Spectral evolution of the X-ray nova XTE J1859+226 during its outburst observed by BeppoSAX and RXTE

Quasi-simultaneous Integral, Swift, And Nustar Observations Of The New X-ray Clocked Burster 1rxsj180408.9-342058

$\texttt{AXP4}$, a Numerical Simulator of Pulse Profiles for Binary Accreting X-ray Pulsars $-$ I: Beam Geometries