The thermal X-ray spectra of several isolated neutron stars display deviations from a pure blackbody. The accurate physical interpretation of these spectral features bears profound implications for our understanding of the atmospheric composition, ma
gnetic field strength and topology, and equation of state of dense matter. With specific details varying from source to source, common explanations for the features have ranged from atomic transitions in the magnetized atmospheres or condensed surface, to cyclotron lines generated in a hot ionized layer near the surface. Here we quantitatively evaluate the X-ray spectral distortions induced by inhomogeneous temperature distributions of the neutron star surface. To this aim, we explore several surface temperature distributions, we simulate their corresponding general relativistic X-ray spectra (assuming an isotropic, blackbody emission), and fit the latter with a single blackbody model. We find that, in some cases, the presence of a spurious spectral line is required at a high significance level in order to obtain statistically acceptable fits, with central energy and equivalent width similar to the values typically observed. We also perform a fit to a specific object, RX J0806.4-4123, finding several surface temperature distributions able to model the observed spectrum. The explored effect is unlikely to work in all sources with detected lines, but in some cases it can indeed be responsible for the appearance of such lines. Our results enforce the idea that surface temperature anisotropy can be an important factor that should be considered and explored also in combination with more sophisticated emission models like atmospheres.
(Abridged) We present a systematic fit of a model of resonant cyclotron scattering (RCS) to the X-ray data of ten magnetars, including canonical and transient anomalous X-ray pulsars (AXPs), and soft gamma repeaters (SGRs). In this scenario, non-ther
mal magnetar spectra in the soft X-rays (i.e. below ~10 keV) result from resonant cyclotron scattering of the thermal surface emission by hot magnetospheric plasma. We find that this model can successfully account for the soft X-ray emission of magnetars, while using the same number of free parameters than the commonly used empirical blackbody plus power-law model. However, while the RCS model can alone reproduce the soft X-ray spectra of AXPs, the much harder spectra of SGRs below ~10 keV, requires the addition of a power-law component (the latter being the same component responsible for their hard X-ray emission). Although this model in its present form does not explain the hard X-ray emission of a few of these sources, we took this further component into account in our modeling not to overlook their contribution in the ~4-10 keV band. We find that the entire class of sources is characterized by magnetospheric plasma with a density which, at resonant radius, is about 3 orders of magnitudes higher than n_{GJ}, the Goldreich-Julian electron density. The inferred values of the intervening hydrogen column densities, are also in better agreement with more recent estimates inferred from the fit of single X-ray edges. For the entire sample of observations, we find indications for a correlation between the scattering depth and the electron thermal velocity, and the field strength. Moreover, in most transient anomalous X-ray pulsars the outburst state is characterized by a relatively high surface temperature which cools down during the decay.
