3XMM J185246.6+003317: another low magnetic field magnetar

We study the outburst of the newly discovered X-ray transient 3XMMJ185246.6+003317, re-analyzing all available XMM-Newton observations of the source to perform a phase-coherent timing analysis, and derive updated values of the period and period derivative. We find the source rotating at P = 11.55871346(6) s (90% confidence level; at epoch MJD 54728.7) but no evidence for a period derivative in the seven months of outburst decay spanned by the observations. This translates to a 3σ upper limit for the period derivative of ˙ P <1.4×10−13 s s−1, which, assuming the classical magneto-dipolar braking model, gives a limit on the dipolar magnetic field of Bdip < 4.1×1013 G. The X-ray outburst and spectral characteristics of 3XMM J185246.6+003317 confirm its identification as a magnetar, but the magnetic field upper limit we derive defines it as the third “low-B” magnetar discovered in the past 3 yr, after SGR 0418+5729 and Swift J1822.3−1606. We have also obtained an upper limit to the quiescent luminosity (<4×1033 erg s−1), in line with the expectations for an old magnetar. The discovery of this new low field magnetar reaffirms the prediction of about one outburst per year from the hidden population of aged magnetars.

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.

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

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.

Neutron Stars—Thermal Emitters

Confronting theoretical models with observations of thermal radiation emitted by neutron stars is one of the most important ways to understand the properties of both, superdense matter in the interiors of the neutron stars and dense magnetized plasmas in their outer layers. Here we review the theory of thermal emission from the surface layers of strongly magnetized neutron stars, and the main properties of the observational data. In particular, we focus on the nearby sources for which a clear thermal component has been detected, without being contaminated by other emission processes (magnetosphere, accretion, nebulae). We also discuss the applications of the modern theoretical models of the formation of spectra of strongly magnetized neutron stars to the observed thermally emitting objects.

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.

Pulsar science with the SKA

The unprecedented sensitivity and large field of view of SKA will be of paramount importance for pulsar science, and for many related research fields. In particular, beside the obvious discovery of many more pulsars (even those with very low luminosity), and the extremely accurate timing analysis of the current pulsar population, SKA will allow to use pulsars to measure or put strong constraints on gravitational waves, Galactic magnetism, planet masses, general relativity and nuclear physics.

Quiescent thermal emission from neutron stars in low-mass X-ray binaries

Context. We monitored the quiescent thermal emission from neutron stars in low-mass X-ray binaries after active periods of intense activity in X-rays (outbursts). Aims. The theoretical modeling of the thermal relaxation of the neutron star crust may be used to establish constraints on the crust composition and transport properties, depending on the astrophysical scenarios assumed. Methods. We numerically simulated the thermal evolution of the neutron star crust and compared them with inferred surface temperatures for five sources: MXB 1659−29, KS 1731−260, XTE J1701−462, EXO 0748−676  and IGR J17480−2446. Results. We find that the evolution of MXB 1659−29, KS 1731−260 and EXO 0748−676 can be well described within a deep crustal cooling scenario. Conversely, we find that the other two sources can only be explained with models beyond crustal cooling. For the peculiar emission of XTE J1701−462 we propose alternative scenarios such as residual accretion during quiescence, additional heat sources in the outer crust, and/or thermal isolation of the inner crust due to a buried magnetic field. We also explain the very recent reported temperature of IGR J17480−2446 with an additional heat deposition in the outer crust from shallow sources.

Magnetars: Properties, Origin and Evolution

Magnetars are neutron stars in which a strong magnetic field is the main energy source. About two dozens of magnetars, plus several candidates, are currently known in our Galaxy and in the Magellanic Clouds. They appear as highly variable X-ray sources and, in some cases, also as radio and/or optical pulsars. Their spin periods (2–12 s) and spin-down rates (∼10−13–10−10 s s−1) indicate external dipole fields of ∼1013−15 G, and there is evidence that even stronger magnetic fields are present inside the star and in non-dipolar magnetospheric components. Here we review the observed properties of the persistent emission from magnetars, discuss the main models proposed to explain the origin of their magnetic field and present recent developments in the study of their evolution and connection with other classes of neutron stars.

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.

Neutron Stars—Cooling and Transport

Observations of thermal radiation from neutron stars can potentially provide information about the states of supranuclear matter in the interiors of these stars with the aid of the theory of neutron-star thermal evolution. We review the basics of this theory for isolated neutron stars with strong magnetic fields, including most relevant thermodynamic and kinetic properties in the stellar core, crust, and blanketing envelopes.

Constraining the GRB-Magnetar Model by Means of the Galactic Pulsar Population

A large fraction of Gamma-ray bursts (GRBs) displays an X-ray plateau phase within <105 s from the prompt emission, proposed to be powered by the spin-down energy of a rapidly spinning newly born magnetar. In this work we use the properties of the Galactic neutron star population to constrain the GRB-magnetar scenario. We re-analyze the X-ray plateaus of all Swift GRBs with known redshift, between 2005 January and 2014 August. From the derived initial magnetic field distribution for the possible magnetars left behind by the GRBs, we study the evolution and properties of a simulated GRB-magnetar population using numerical simulations of magnetic field evolution, coupled with Monte Carlo simulations of Pulsar Population Synthesis in our Galaxy. We find that if the GRB X-ray plateaus are powered by the rotational energy of a newly formed magnetar, the current observational properties of the Galactic magnetar population are not compatible with being formed within the GRB scenario (regardless of the GRB type or rate at z = 0). Direct consequences would be that we should allow the existence of magnetars and "super-magnetars" having different progenitors, and that Type Ib/c SNe related to Long GRBs form systematically neutron stars with higher initial magnetic fields. We put an upper limit of ≤16 "super-magnetars" formed by a GRB in our Galaxy in the past Myr (at 99% c.l.). This limit is somewhat smaller than what is roughly expected from Long GRB rates, although the very large uncertainties do not allow us to draw strong conclusion in this respect.

Population synthesis of isolated neutron stars with magneto-rotational evolution – II. From radio-pulsars to magnetars

Population synthesis studies constitute a powerful method to reconstruct the birth distribution of periods and magnetic fields of the pulsar population. When this method is applied to populations in different wavelengths, it can break the degeneracy in the inferred properties of initial distributions that arises from single-band studies. In this context, we extend previous works to include X-ray thermal emitting pulsars within the same evolutionary model as radio-pulsars. We find that the cumulative distribution of the number of X-ray pulsars can be well reproduced by several models that, simultaneously, reproduce the characteristics of the radio-pulsar distribution. However, even considering the most favourable magneto-thermal evolution models with fast field decay, lognormal distributions of the initial magnetic field overpredict the number of visible sources with periods longer than 12 s. We then show that the problem can be solved with different distributions of magnetic field, such as a truncated lognormal distribution, or a binormal distribution with two distinct populations. We use the observational lack of isolated neutron stars (NSs) with spin periods P > 12 s to establish an upper limit to the fraction of magnetars born with B > 1015 G (less than 1 per cent). As future detections keep increasing the magnetar and high-B pulsar statistics, our approach can be used to establish a severe constraint on the maximum magnetic field at birth of NSs.