Unifying the observational diversity of isolated neutron stars via magneto-thermal evolution models

Observations of magnetars and some of the high magnetic field pulsars have shown that their thermal luminosity is systematically higher than that of classical radio-pulsars, thus confirming the idea that magnetic fields are involved in their X-ray emission. Here we present the results of 2D simulations of the fully coupled evolution of temperature and magnetic field in neutron stars, including the state-of-the-art kinetic coefficients and, for the first time, the important effect of the Hall term. After gathering and thoroughly re-analysing in a consistent way all the best available data on isolated, thermally emitting neutron stars, we compare our theoretical models to a data sample of 40 sources. We find that our evolutionary models can explain the phenomenological diversity of magnetars, high-B radio-pulsars, and isolated nearby neutron stars by only varying their initial magnetic field, mass and envelope composition. Nearly all sources appear to follow the expectations of the standard theoretical models. Finally, we discuss the expected outburst rates and the evolutionary links between different classes. Our results constitute a major step towards the grand unification of the isolated neutron star zoo.

The imprint of the crustal magnetic field on the thermal spectra and pulse profiles of isolated neutron stars

Isolated neutron stars (NSs) show a bewildering variety of astrophysical manifestations, presumably shaped by the magnetic field strength and topology at birth. Here, using state-of-the-art calculations of the coupled magnetic and thermal evolution of NSs, we compute the thermal spectra and pulse profiles expected for a variety of initial magnetic field configurations. In particular, we contrast models with purely poloidal magnetic fields to models dominated by a strong internal toroidal component. We find that, while the former displays double-peaked profiles and very low pulsed fractions, in the latter, the anisotropy in the surface temperature produced by the toroidal field often results in a single pulse profile, with pulsed fractions that can exceed the 50–60 per cent level even for perfectly isotropic local emission. We further use our theoretical results to generate simulated ‘observed’ spectra, and show that blackbody (BB) fits result in inferred radii that can be significantly smaller than the actual NS radius, even as low as ∼1–2 km for old NSs with strong internal toroidal fields and a high absorption column density along their line of sight. We compute the size of the inferred BB radius for a few representative magnetic field configurations, NS ages and magnitudes of the column density. Our theoretical results are of direct relevance to the interpretation of X-ray observations of isolated NSs, as well as to the constraints on the equation of state of dense matter through radius measurements.

Accreting magnetars: a new type of high-mass X-ray binaries?

The discovery of very slow pulsations (Pspin =5560 s) has solved the long-standing question of the nature of the compact object in the high-mass X-ray binary 4U 2206+54 but has posed new ones. According to spin evolutionary models in close binary systems, such slow pulsations require a neutron star magnetic field strength larger than the quantum critical value of 4.4 × 1013 G, suggesting the presence of a magnetar. We present the first XMM–Newton observations of 4U 2206+54 and investigate its spin evolution. We find that the observed spin-down rate agrees with the magnetar scenario. We analyse Integral Spacecraft Gamma-Ray Imager (ISGRI)/INTErnational Gamma-RAy Laboratory (INTEGRAL) observations of 4U 2206+54 to search for the previously suggested cyclotron resonance scattering feature at ∼30 keV. We do not find a clear indication of the presence of the line, although certain spectra display shallow dips, not always at 30 keV. The association of these dips with a cyclotron line is very dubious because of its apparent transient nature. We also investigate the energy spectrum of 4U 2206+54 in the energy range 0.3–10 keV with unprecedented detail and report for the first time the detection of very weak 6.5 keV fluorescence iron lines. The photoelectric absorption is consistent with the interstellar value, indicating very small amount of local matter, which would explain the weakness of the florescence lines. The lack of matter locally to the source may be the consequence of the relatively large orbital separation of the two components of the binary. The wind would be too tenuous in the vicinity of the neutron star.

Central compact objects and the hidden magnetic field scenario

Central compact objects (CCOs) are X-ray sources lying close to the centre of supernova remnants, with inferred values of the surface magnetic fields significantly lower (≲1011 G) than those of standard pulsars. In this paper, we revise the hidden magnetic field scenario, presenting the first 2D simulations of the submergence and re-emergence of the magnetic field in the crust of a neutron star. A post-supernova accretion stage of about 10−4–10−3 M⊙ over a vast region of the surface is required to bury the magnetic field into the inner crust. When accretion stops, the field re-emerges on a typical time-scale of 1–100 kyr, depending on the submergence conditions. After this stage, the surface magnetic field is restored close to its birth values. A possible observable consequence of the hidden magnetic field is the anisotropy of the surface temperature distribution, in agreement with observations of several of these sources. We conclude that the hidden magnetic field model is viable as an alternative to the antimagnetar scenario, and it could provide the missing link between CCOs and the other classes of isolated neutron stars.

The evolution and masses of the neutron star and donor star in the high mass X-ray binary OAO 1657−415

We report near-infrared radial velocity (RV) measurements of the recently identified donor star in the high mass X-ray binary (HMXB) system OAO 1657−415 obtained in the H band using ISAAC on the Very Large Telescope. Cross-correlation methods were employed to construct a RV curve with a semi-amplitude of 22.1 ± 3.5 km s−1. Combined with other measured parameters of this system it provides a dynamically determined neutron star (NS) mass of 1.42 ± 0.26 M⊙ and a mass of 14.3 ± 0.8 M⊙ for the Ofpe/WN9 highly evolved donor star. OAO 1657−415 is an eclipsing HMXB pulsar with the largest eccentricity and orbital period of any within its class. Of the 10 known eclipsing X-ray binary pulsars OAO 1657−415 becomes the ninth with a dynamically determined NS mass solution and only the second in an eccentric system. Furthermore, the donor star in OAO 1657−415 is much more highly evolved than the majority of the supergiant donors in other HMXBs, joining a small but growing list of HMXBs donors with extensive hydrogen depleted atmospheres. Considering the evolutionary development of OAO 1657−415, we have estimated the binding energy of the envelope of the mass donor and find that there is insufficient energy for the removal of the donor’s envelope via spiral-in, ruling out a common envelope evolutionary scenario. With its non-zero eccentricity and relatively large orbital period the identification of a definitive evolutionary pathway for OAO 1657−415 remains problematic, we conclude by proposing two scenarios which may account for OAO 1657−415 current orbital configuration.

Fast X-ray transients towards the Galactic bulge with the Rossi X-ray Timing Explorer

In X-ray binaries, rapid variability in X-ray flux of greater than an order of magnitude on time-scales of a day or less appears to be a signature of wind accretion from a supergiant companion. When the variability takes the form of rare, brief, bright outbursts with only faint emission between them, the systems are called supergiant fast X-ray transients (SFXTs). We present data from twice-weekly scans of the Galactic bulge by the Rossi X-ray Timing Explorer that allow us to compare the behaviour of known SFXTs and possible SFXT candidates with the persistently bright supergiant X-ray binary 4U 1700−377. We independently confirm the orbital periods reported by other groups for SFXTs SAX J1818.6−1703 and IGR J17544−2619. The new data do not independently reproduce the orbital period reported for XTE J1739−302, but slightly improve the significance of the original result when the data are combined. The bulge source XTE J1743−363 shows a combination of fast variability and a long-term decline in activity, the latter behaviour not being characteristic of supergiant X-ray binaries. A far-red spectrum of the companion suggests that it is a symbiotic neutron star binary rather than a high-mass binary, and the reddest known of this class: the spectral type is approximately M8 III.

Quiescent state and outburst evolution of SGR 0501+4516

We report on the quiescent state of the soft gamma repeater SGR 0501+4516 observed by XMM–Newton on 2009 August 30. The source exhibits an absorbed flux ∼75 times lower than that measured at the peak of the 2008 outburst, and a rather soft spectrum, with the same value of the blackbody temperature observed with ROSAT back in 1992. This new observation is put into the context of all existing X-ray data since its discovery in 2008 August, allowing us to complete the study of the timing and spectral evolution of the source from outburst until its quiescent state. The set of deep XMM–Newton observations performed during the few years time-scale of its outburst allows us to monitor the spectral characteristics of this magnetar as a function of its rotational period, and their evolution along these years. After the first ∼10 d, the initially hot and bright surface spot progressively cooled down during the decay. We discuss the behaviour of this magnetar in the context of its simulated secular evolution, inferring a plausible dipolar field at birth of 3 × 1014 G, and a current (magnetothermal) age of ∼10 kyr.

Spatiotemporal characteristics of the Huangtupo landslide in the Three Gorges region (China) constrained by radar interferometry

The Huangtupo landslide is one of the largest in the Three Gorges region, China. The county-seat town of Badong, located on the south shore between the Xiling and Wu gorges of the Yangtze River, was moved to this unstable slope prior to the construction of the Three Gorges Project, since the new Three Gorges reservoir completely submerged the location of the old city. The instability of the slope is affecting the new town by causing residential safety problems. The Huangtupo landslide provides scientists an opportunity to understand landslide response to fluctuating river water level and heavy rainfall episodes, which is essential to decide upon appropriate remediation measures. Interferometric Synthetic Aperture Radar (InSAR) techniques provide a very useful tool for the study of superficial and spatially variable displacement phenomena. In this paper, three sets of radar data have been processed to investigate the Huangtupo landslide. Results show that maximum displacements are affecting the northwest zone of the slope corresponding to Riverside slumping mass I#. The other main landslide bodies (i.e. Riverside slumping mass II#, Substation landslide and Garden Spot landslide) exhibit a stable behaviour in agreement with in situ data, although some active areas have been recognized in the foot of the Substation landslide and Garden Spot landslide. InSAR has allowed us to study the kinematic behaviour of the landslide and to identify its active boundaries. Furthermore, the analysis of the InSAR displacement time-series has helped recognize the different displacement patterns on the slope and their relationships with various triggering factors. For those persistent scatterers, which exhibit long-term displacements, they can be decomposed into a creep model (controlled by geological conditions) and a superimposed recoverable term (dependent on external factors), which appears closely correlated with reservoir water level changes close to the river's edge. These results, combined with in situ data, provide a comprehensive analysis of the Huangtupo landslide, which is essential for its management.

Pulse phase-coherent timing and spectroscopy of CXOU J164710.2−45521 outbursts

We present a long-term phase-coherent timing analysis and pulse-phase resolved spectroscopy for the two outbursts observed from the transient anomalous X-ray pulsar CXOU J164710.2−455216. For the first outburst we used 11 Chandra and XMM–Newton observations between 2006 September and 2009 August, the longest baseline yet for this source. We obtain a coherent timing solution with P = 10.61065583(4) s, Ṗ = 9.72(1) × 10−13 s s−1 and P̈ = –1.05(5) × 10−20 s s−2. Under the standard assumptions this implies a surface dipolar magnetic field of ∼1014 G, confirming this source as a standard B magnetar. We also study the evolution of the pulse profile (shape, intensity and pulsed fraction) as a function of time and energy. Using the phase-coherent timing solution we perform a phase-resolved spectroscopy analysis, following the spectral evolution of pulse-phase features, which hints at the physical processes taking place on the star. The results are discussed from the perspective of magnetothermal evolution models and the untwisting magnetosphere model. Finally, we present similar analysis for the second, less intense, 2011 outburst. For the timing analysis we used Swift data together with 2 XMM–Newton and Chandra pointings. The results inferred for both outbursts are compared and briefly discussed in a more general framework.

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).

On the decorrelation filtering of RL05 GRACE data for global applications

In autumn 2012, the new release 05 (RL05) of monthly geopotencial spherical harmonics Stokes coefficients (SC) from GRACE (Gravity Recovery and Climate Experiment) mission was published. This release reduces the noise in high degree and order SC, but they still need to be filtered. One of the most common filtering processing is the combination of decorrelation and Gaussian filters. Both of them are parameters dependent and must be tuned by the users. Previous studies have analyzed the parameters choice for the RL05 GRACE data for oceanic applications, and for RL04 data for global application. This study updates the latter for RL05 data extending the statistics analysis. The choice of the parameters of the decorrelation filter has been optimized to: (1) balance the noise reduction and the geophysical signal attenuation produced by the filtering process; (2) minimize the differences between GRACE and model-based data; (3) maximize the ratio of variability between continents and oceans. The Gaussian filter has been optimized following the latter criteria. Besides, an anisotropic filter, the fan filter, has been analyzed as an alternative to the Gauss filter, producing better statistics.

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.

First EURONEAR NEA discoveries from La Palma using the INT

Since 2006, the European Near Earth Asteroids Research (EURONEAR) project has been contributing to the research of near-Earth asteroids (NEAs) within a European network. One of the main aims is the amelioration of the orbits of NEAs, and starting in 2014 February we focus on the recovery of one-opposition NEAs using the Isaac Newton Telescope (INT) in La Palma in override mode. Part of this NEA recovery project, since 2014 June EURONEAR serendipitously started to discover and secure the first NEAs from La Palma and using the INT, thanks to the teamwork including amateurs and students who promptly reduce the data, report discoveries and secure new objects recovered with the INT and few other telescopes from the EURONEAR network. Five NEAs were discovered with the INT, including 2014 LU14, 2014 NL52 (one very fast rotator), 2014 OL339 (the fourth known Earth quasi-satellite), 2014 SG143 (a quite large NEA), and 2014 VP. Another very fast moving NEA was discovered but was unfortunately lost due to lack of follow-up time. Additionally, another 14 NEA candidates were identified based on two models, all being rapidly followed-up using the INT and another 11 telescopes within the EURONEAR network. They include one object discovered by Pan-STARRS, two Mars crossers, two Hungarias, one Jupiter trojan, and other few inner main belt asteroids (MBAs). Using the INT and Sierra Nevada 1.5 m for photometry, then the Gran Telescopio de Canarias for spectroscopy, we derived the very rapid rotation of 2014 NL52, then its albedo, magnitude, size, and its spectral class. Based on the total sky coverage in dark conditions, we evaluate the actual survey discovery rate using 2-m class telescopes. One NEA is possible to be discovered randomly within minimum 2.8 deg2 and maximum 5.5 deg2. These findings update our past statistics, being based on double sky coverage and taking into account the recent increase in discovery.

The X-ray outburst of the Galactic Centre magnetar SGR J1745−2900 during the first 1.5 year

In 2013 April a new magnetar, SGR 1745−2900, was discovered as it entered an outburst, at only 2.4 arcsec angular distance from the supermassive black hole at the centre of the Milky Way, Sagittarius A*. SGR 1745−2900 has a surface dipolar magnetic field of ∼2 × 1014 G, and it is the neutron star closest to a black hole ever observed. The new source was detected both in the radio and X-ray bands, with a peak X-ray luminosity LX ∼ 5 × 1035 erg s−1. Here we report on the long-term Chandra (25 observations) and XMM–Newton (eight observations) X-ray monitoring campaign of SGR 1745−2900 from the onset of the outburst in 2013 April until 2014 September. This unprecedented data set allows us to refine the timing properties of the source, as well as to study the outburst spectral evolution as a function of time and rotational phase. Our timing analysis confirms the increase in the spin period derivative by a factor of ∼2 around 2013 June, and reveals that a further increase occurred between 2013 October 30 and 2014 February 21. We find that the period derivative changed from 6.6 × 10−12 to 3.3 × 10−11 s s−1 in 1.5 yr. On the other hand, this magnetar shows a slow flux decay compared to other magnetars and a rather inefficient surface cooling. In particular, starquake-induced crustal cooling models alone have difficulty in explaining the high luminosity of the source for the first ∼200 d of its outburst, and additional heating of the star surface from currents flowing in a twisted magnetic bundle is probably playing an important role in the outburst evolution.