Spectral features in isolated neutron stars induced by inhomogeneous surface temperatures

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, magnetic 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.

Population synthesis of isolated neutron stars with magneto-rotational evolution

We revisit the population synthesis of isolated radio-pulsars incorporating recent advances on the evolution of the magnetic field and the angle between the magnetic and rotational axes from new simulations of the magneto-thermal evolution and magnetosphere models, respectively. An interesting novelty in our approach is that we do not assume the existence of a death line. We discuss regions in parameter space that are more consistent with the observational data. In particular, we find that any broad distribution of birth spin periods with P0 ≲ 0.5 s can fit the data, and that if the alignment angle is allowed to vary consistently with the torque model, realistic magnetospheric models are favoured compared to models with classical magneto-dipolar radiation losses. Assuming that the initial magnetic field is given by a lognormal distribution, our optimal model has mean strength 〈log B0[G]〉 ≈ 13.0–13.2 with width σ(log B0) = 0.6–0.7. However, there are strong correlations between parameters. This degeneracy in the parameter space can be broken by an independent estimate of the pulsar birth rate or by future studies correlating this information with the population in other observational bands (X-rays and γ-rays).

The many lives of magnetized neutron stars

The magnetic field strength at birth is arguably one of the most important properties to determine the evolutionary path of a neutron star. Objects with very high fields, collectively known as magnetars, are characterized by high X-ray quiescent luminosities, occurrence of outbursts, and, for some of them, sporadic giant flares. While the magnetic field strength is believed to drive their collective behaviour, however, the diversity of their properties, and, especially, the observation of magnetar-like bursts from “low-field” pulsars, has been a theoretical puzzle. In this review, we discuss results of long-term simulations following the coupled evolution of the X-ray luminosity and the timing properties for a large, homogeneous sample of X-ray emitting isolated neutron stars, accounting for a range of initial magnetic field strengths, envelope compositions, and neutron star masses. In addition, by following the evolution of magnetic stresses within the neutron star crust, we can also relate the observed magnetar phenomenology to the physical properties of neutron stars, and in particular to their age and magnetic field strength and topology. The dichotomy of “high-B” field pulsars versus magnetars is naturally explained, and occasional outbursts from old, low B-field neutron stars are predicted. We conclude by speculating on the fate of old magnetars, and by presenting observational diagnostics of the neutron star crustal field topology.

Discussing the physical meaning of the absorption feature at 2.1 keV in 4U 1538–52

High resolution X-ray spectroscopy is a powerful tool for studying the nature of the matter surrounding the neutron star in X-ray binaries and its interaction between the stellar wind and the compact object. In particular, absorption features in their spectra could reveal the presence of atmospheres of the neutron star or their magnetic field strength. Here we present an investigation of the absorption feature at 2.1 keV in the X-ray spectrum of the high mass X-ray binary 4U 1538–52 based on our previous analysis of the XMM-Newton data. We study various possible origins and discuss the different physical scenarios in order to explain this feature. A likely interpretation is that the feature is associated with atomic transitions in an O/Ne neutron star atmosphere or of hydrogen and helium like Fe or Si ions formed in the stellar wind of the donor.

Probing large-scale wind structures in Vela X–1 using off-states with INTEGRAL

Vela X–1 is the prototype of the class of wind-fed accreting pulsars in high-mass X-ray binaries hosting a supergiant donor. We have analysed in a systematic way 10 years of INTEGRAL data of Vela X–1 (22–50 keV) and we found that when outside the X-ray eclipse, the source undergoes several luminosity drops where the hard X-rays luminosity goes below ∼3 × 1035 erg s−1, becoming undetected by INTEGRAL. These drops in the X-ray flux are usually referred to as ‘off-states’ in the literature. We have investigated the distribution of these off-states along the Vela X–1 ∼ 8.9 d orbit, finding that their orbital occurrence displays an asymmetric distribution, with a higher probability to observe an off-state near the pre-eclipse than during the post-eclipse. This asymmetry can be explained by scattering of hard X-rays in a region of ionized wind, able to reduce the source hard X-ray brightness preferentially near eclipse ingress. We associate this ionized large-scale wind structure with the photoionization wake produced by the interaction of the supergiant wind with the X-ray emission from the neutron star. We emphasize that this observational result could be obtained thanks to the accumulation of a decade of INTEGRAL data, with observations covering the whole orbit several times, allowing us to detect an asymmetric pattern in the orbital distribution of off-states in Vela X–1.

Orbital phase-resolved spectroscopy of 4U 1538−52 with MAXI

Context. 4U 1538−52, an absorbed high mass X-ray binary with an orbital period of ~3.73 days, shows moderate orbital intensity modulations with a low level of counts during the eclipse. Several models have been proposed to explain the accretion at different orbital phases by a spherically symmetric stellar wind from the companion. Aims. The aim of this work is to study both the light curve and orbital phase spectroscopy of this source in the long term. In particular, we study the folded light curve and the changes in the spectral parameters with orbital phase to analyse the stellar wind of QV Nor, the mass donor of this binary system. Methods. We used all the observations made from the Gas Slit Camera on board MAXI of 4U 1538−52 covering many orbits continuously. We obtained the good interval times for all orbital phase ranges, which were the input for extracting our data. We estimated the orbital period of the system and then folded the light curves, and we fitted the X-ray spectra with the same model for every orbital phase spectrum. We also extracted the averaged spectrum of all the MAXI data available. Results. The MAXI spectra in the 2–20 keV energy range were fitted with an absorbed Comptonisation of cool photons on hot electrons. We found a strong orbital dependence of the absorption column density but neither the fluorescence iron emission line nor low energy excess were needed to fit the MAXI spectra. The variation in the spectral parameters over the binary orbit were used to examine the mode of accretion onto the neutron star in 4U 1538−52. We deduce a best value of Ṁ/v∞ = 0.65 × 10-9M⊙ yr-1/ (km s-1) for QV Nor.

On the formation and evolution of the first Be star in a black hole binary MWC 656

We find that the formation of MWC 656 (the first Be binary containing a black hole) involves a common envelope phase and a supernova explosion. This result supports the idea that a rapidly rotating Be star can emerge out of a common envelope phase, which is very intriguing because this evolutionary stage is thought to be too fast to lead to significant accretion and spin up of the B star. We predict ∼10–100 of B-BH binaries to currently reside in the Galactic disc, among which around 1/3 contain a Be star, but there is only a small chance to observe a system with parameters resembling MWC 656. If MWC 656 is representative of intrinsic Galactic Be-BH binary population, it may indicate that standard evolutionary theory needs to be revised. This would pose another evolutionary problem in understanding black hole (BH) binaries, with BH X-ray novae formation issue being the prime example. Future evolution of MWC 656 with an ∼5 M⊙ BH and with an ∼13 M⊙ main-sequence companion on an ∼60 d orbit may lead to the formation of a coalescing BH–NS (neutron star) system. The estimated Advanced LIGO/Virgo detection rate of such systems is up to ∼0.2 yr−1. This empirical estimate is a lower limit as it is obtained with only one particular evolutionary scenario, the MWC 656 binary. This is only a third such estimate available (after Cyg X-1 and Cyg X-3), and it lends additional support to the existence of so far undetected BH–NS binaries.

On the Radial Onset of Clumping in the Wind of the B0I Massive Star QV Nor

We present an analysis of a 78 ks Chandra high-energy transmission gratings observation of the B0I star QV Nor, the massive donor of the wind-accreting pulsar 4U1538−52. The neutron star (NS) orbits its companion in a very close orbit (r < 1.4R*, in units of the stellar radii), thereby allowing probing of the innermost wind regions. The flux of the Fe Kα line during eclipse reduces to only ∼30% of the flux measured out of eclipse. This indicates that the majority of Fe fluorescence must be produced in regions close to the NS, at distances smaller than 1R* from its surface. The fact that the flux of the continuum decreases to only ∼3% during eclipse allows for a high contrast of the Fe Kα line fluorescence during eclipse. The line is not resolved and centered at 1.9368 0.0018 l = 0.0032 - + Å. From the inferred plasma speed limit of v < c l < 800 l D km s−1 and range of ionization parameters of log 1, 2 x = [- ], together with the stellar density profile, we constrain the location of the cold, dense material in the stellar wind of QV Nor using simple geometrical considerations. We then use the Fe Kα line fluorescence as a tracer of wind clumps and determine that these clumps in the stellar wind of QV Nor (B0I) must already be present at radii r < 1.25R*, close to the photosphere of the star.

Modelling of the surface emission of the low magnetic field magnetar SGR 0418+5729

We perform a detailed modelling of the post-outburst surface emission of the low magnetic field magnetar SGR 0418+5729. The dipolar magnetic field of this source, B=6×1012G estimated from its spin-down rate, is in the observed range of magnetic fields for normal pulsars. The source is further characterized by a high pulse fraction and a single-peak profile. Using synthetic temperature distribution profiles, and fully accounting for the general-relativistic effects of light deflection and gravitational redshift, we generate synthetic X-ray spectra and pulse profiles that we fit to the observations. We find that asymmetric and symmetric surface temperature distributions can reproduce equally well the observed pulse profiles and spectra of SGR 0418. None the less, the modelling allows us to place constraints on the system geometry (i.e. the angles ψ and ξ that the rotation axis makes with the line of sight and the dipolar axis, respectively), as well as on the spot size and temperature contrast on the neutron star surface. After performing an analysis iterating between the pulse profile and spectra, as done in similar previous works, we further employed, for the first time in this context, a Markov-Chain Monte Carlo approach to extract constraints on the model parameters from the pulse profiles and spectra, simultaneously. We find that, to reproduce the observed spectrum and flux modulation: (a) the angles must be restricted to 65° ≲ ψ + ξ ≲ 125° or 235° ≲ ψ + ξ ≲ 295°; (b) the temperature contrast between the poles and the equator must be at least a factor of ∼6, and (c) the size of the hottest region ranges between 0.2 and 0.7 km (including uncertainties on the source distance). Lastly, we interpret our findings within the context of internal and external heating models.

The outburst decay of the low magnetic field magnetar SWIFT J1822.3−1606: phase-resolved analysis and evidence for a variable cyclotron feature

We study the timing and spectral properties of the low-magnetic field, transient magnetar SWIFT J1822.3−1606 as it approached quiescence. We coherently phase-connect the observations over a time-span of ∼500 d since the discovery of SWIFT J1822.3−1606 following the Swift-Burst Alert Telescope (BAT) trigger on 2011 July 14, and carried out a detailed pulse phase spectroscopy along the outburst decay. We follow the spectral evolution of different pulse phase intervals and find a phase and energy-variable spectral feature, which we interpret as proton cyclotron resonant scattering of soft photon from currents circulating in a strong (≳1014 G) small-scale component of the magnetic field near the neutron star surface, superimposed to the much weaker (∼3 × 1013 G) magnetic field. We discuss also the implications of the pulse-resolved spectral analysis for the emission regions on the surface of the cooling magnetar.

The Goodness of Simultaneous Fits in ISIS

In a previous work, we introduced a tool for analyzing multiple datasets simultaneously, which has been implemented into ISIS. This tool was used to fit many spectra of X-ray binaries. However, the large number of degrees of freedom and individual datasets raise an issue about a good measure for a simultaneous fit quality. We present three ways to check the goodness of these fits: we investigate the goodness of each fit in all datasets, we define a combined goodness exploiting the logical structure of a simultaneous fit, and we stack the fit residuals of all datasets to detect weak features. These tools are applied to all RXTE-spectra from GRO 1008−57, revealing calibration features that are not detected significantly in any single spectrum. Stacking the residuals from the best-fit model for the Vela X-1 and XTE J1859+083 data evidences fluorescent emission lines that would have gone undetected otherwise.

The discovery, monitoring and environment of SGR J1935+2154

We report on the discovery of a new member of the magnetar class, SGR J1935+2154, and on its timing and spectral properties measured by an extensive observational campaign carried out between 2014 July and 2015 March with Chandra and XMM–Newton (11 pointings). We discovered the spin period of SGR J1935+2154 through the detection of coherent pulsations at a period of about 3.24 s. The magnetar is slowing down at a rate of P˙=1.43(1)×10−11 s s−1 and with a decreasing trend due to a negative P¨ of −3.5(7) × 10−19 s s−2. This implies a surface dipolar magnetic field strength of ∼2.2 × 1014 G, a characteristic age of about 3.6 kyr and a spin-down luminosity Lsd ∼1.7 × 1034 erg s−1. The source spectrum is well modelled by a blackbody with temperature of about 500 eV plus a power-law component with photon index of about 2. The source showed a moderate long-term variability, with a flux decay of about 25 per cent during the first four months since its discovery, and a re-brightening of the same amount during the second four months. The X-ray data were also used to study the source environment. In particular, we discovered a diffuse emission extending on spatial scales from about 1 arcsec up to at least 1 arcmin around SGR J1935+2154 both in Chandra and XMM–Newton data. This component is constant in flux (at least within uncertainties) and its spectrum is well modelled by a power-law spectrum steeper than that of the pulsar. Though a scattering halo origin seems to be more probable we cannot exclude that part, or all, of the diffuse emission is due to a pulsar wind nebula.