Off-centred force-free neutron star magnetospheres

Neutron star electromagnetic activity produces pairs that fill their magnetosphere represented to the zeroth order by the force-free approximation. Neither dissipation nor acceleration nor radiation from charged particles is expected from this simplified model. So far, only centred dipole magnetic fields have been studied in this limit. In this paper, we explore the consequences of a rotating off-centred dipole on the force-free magnetosphere, showing the new magnetic field geometry, its spin-down luminosity as well as the electromagnetic kick and torque felt by the neutron star. Solutions are obtained by time-dependent numerical simulations of the force-free regime using our pseudo-spectral code written in spherical coordinates. Our results are also compared to known analytical expressions found for the off-centred vacuum dipole by an expansion to lowest order in the parameter ε = d/R, where d is the displacement of the dipole from the stellar centre and R the neutron star radius. The presence of a force-free plasma enhances the spin-down luminosity as well as the electromagnetic kick and torque with respect to a centred force-free dipole. The impact on isolated and binary neutron stars is revised in light of these new results.

Relativistic rotating vector model for X-ray millisecond pulsars

The X-ray radiation produced on the surface of accreting magnetised neutron stars is expected to be strongly polarised. A swing of the polarisation vector with the pulsar phase gives a direct measure of the source inclination and magnetic obliquity. In the case of rapidly rotating millisecond pulsars, the relativistic motion of the emission region causes additional rotation of the polarisation plane. Here, we develop a relativistic rotating vector model, where we derive analytical expression for the polarisation angle as a function of the pulsar phase accounting for relativistic aberration and gravitational light bending in the Schwarzschild metric. We show that in the case of fast pulsars the rotation of the polarisation plane can reach tens of degrees, strongly influencing the observed shape of the polarisation angle’s phase dependence. The rotation angle grows nearly linearly with the spin rate but it is less sensitive to the neutron star radius. Overall, this angle is large even for large spots. Our results have implications with regard to the modelling of X-ray polarisation from accreting millisecond pulsars that are to be observed with the upcoming Imaging X-ray Polarimeter Explorer and the enhanced X-ray Timing and Polarimetry mission. The X-ray polarisation may improve constraints on the neutron star mass and radius coming from the pulse profile modelling.

Neutron star radius measurement from the ultraviolet and soft X-ray thermal emission of PSR J0437−4715

ABSTRACT We analysed the thermal emission from the entire surface of the millisecond pulsar PSR J0437−4715 observed in the ultraviolet and soft X-ray bands. For this, we calculated non-magnetized, partially ionized atmosphere models of hydrogen, helium, and iron compositions and included plasma frequency effects that may affect the emergent spectrum. This is particularly true for the coldest atmospheres composed of iron (up to a few per cent changes in the soft X-ray flux). Employing a Markov chain Monte Carlo method, we found that the spectral fits favour a hydrogen atmosphere, disfavour a helium composition, and rule out iron atmosphere and blackbody models. By using a Gaussian prior on the dust extinction, based on the latest 3D map of Galactic dust, and accounting for the presence of hot polar caps found in the previous work, we found that the hydrogen atmosphere model results in a well-constrained neutron star radius ${R_{\rm NS}}= 13.6^{+0.9}_{-0.8}{\, {\rm km}}$ and bulk surface temperature ${T_{\rm eff}^{\infty }}=\left(2.3\pm 0.1\right){\times 10^{5}}{\, {\rm K}}$. This relatively large radius favours a stiff equation of state and disfavours a strange quark composition inside neutron stars.

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

ABSTRACT We present a new method of matching observations of Type-I (thermonuclear) X-ray bursts with models, comparing the predictions of a semi-analytic ignition model with X-ray observations of the accretion-powered millisecond pulsar SAX J1808.4–3658 in outburst. We used a Bayesian analysis approach to marginalize over the parameters of interest and determine parameters such as fuel composition, distance/anisotropy factors, neutron star mass, and neutron star radius. Our study includes a treatment of the system inclination effects, inferring that the rotation axis of the system is inclined $\left(69^{+4}_{-2}\right)^\circ$ from the observers line of sight, assuming a flat disc model. This method can be applied to any accreting source that exhibits Type-I X-ray bursts. We find a hydrogen mass fraction of $0.57^{+0.13}_{-0.14}$ and CNO metallicity of $0.013^{+0.006}_{-0.004}$ for the accreted fuel is required by the model to match the observed burst energies, for a distance to the source of $3.3^{+0.3}_{-0.2}\, \mathrm{kpc}$. We infer a neutron star mass of $1.5^{+0.6}_{-0.3}\, \mathrm{M}_{\odot }$ and radius of $11.8^{+1.3}_{-0.9}\, \mathrm{km}$ for a surface gravity of $1.9^{+0.7}_{-0.4}\times 10^{14}\, \mathrm{cm}\, \mathrm{s}^{-2}$ for SAX J1808.4–3658.

Effects of Compton scattering on the neutron star radius constraints in rotation-powered millisecond pulsars

The aim of this work is to study the possible effects and biases on the radius constraints for rotation-powered millisecond pulsars when using Thomson approximation to describe electron scattering in the atmosphere models, instead of using exact formulation for Compton scattering. We compare the differences between the two models in the energy spectrum and angular distribution of the emitted radiation. We also analyse a self-generated, synthetic, phase-resolved energy spectrum, based on Compton atmosphere and the most X-ray luminous, rotation-powered millisecond pulsars observed by the Neutron star Interior Composition ExploreR (NICER). We derive constraints for the neutron star parameters using both the Compton and Thomson models. The results show that the method works by reproducing the correct parameters with the Compton model. However, biases are found in both the size and the temperature of the emitting hotspot, when using the Thomson model. The constraints on the radius are still not significantly changed, and therefore the Thomson model seems to be adequate if we are interested only in the radius measurements using NICER.