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.

Constraining the Properties of the Thermonuclear Burst Oscillation Source XTE J1814-338 Through Pulse Profile Modelling

Pulse Profile Modeling of Thermonuclear Burst Oscillations I: The Effect of Neglecting Variability

Pulse Profile Modelling of Thermonuclear Burst Oscillations II: Handling variability

Constraints on Neutron Star Parameters from Burst Oscillation Light Curves of the Accreting Millisecond Pulsar XTE J1814-338

Timing Features of the Accretion‐driven Millisecond X‐Ray Pulsar XTE J1807−294 in the 2003 March Outburst

Simulating X-ray bursts during a transient accretion event

Inferring system parameters from the bursts of the accretion-powered pulsar IGR J17498-2921

A Bayesian approach to matching thermonuclear X-ray burst observations with models

Diagnosing the burst influence upon accretion in the clocked burster GS 1826--238

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

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

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

A NICER View of PSR J1231-1411: A Complex Case

A NICER View of the Nearest and Brightest Millisecond Pulsar: PSR J0437$\unicode{x2013}$4715

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

Theoretical Analysis of the RX J0209.6−7427 X-Ray Spectrum during a Giant Outburst

A comprehensive study of type I (thermonuclear) bursts in the new transient SRGA J144459.2$-$604207

Models of Type I X-ray Bursts from GS 1826-24: A Probe of rp-Process Hydrogen Burning

Pulse Profiles of Accreting Neutron Stars from GRMHD Simulations

Numerical $\texttt{AXP4}$ Simulations of Pulse Profiles for Binary Accreting X-ray Pulsars $-$ II: A Case Study of Centaurus X-3

Pulse frequency variations and timing noise of MXB 0656–072 during the 2007–2008 type I outbursts and implications for its magnetic field