A new code for the Hall-driven magnetic evolution of neutron stars

Over the past decade, the numerical modeling of the magnetic field evolution in astrophysical scenarios has become an increasingly important field. In the crystallized crust of neutron stars the evolution of the magnetic field is governed by the Hall induction equation. In this equation the relative contribution of the two terms (Hall term and Ohmic dissipation) varies depending on the local conditions of temperature and magnetic field strength. This results in the transition from the purely parabolic character of the equations to the hyperbolic regime as the magnetic Reynolds number increases, which presents severe numerical problems. Up to now, most attempts to study this problem were based on spectral methods, but they failed in representing the transition to large magnetic Reynolds numbers. We present a new code based on upwind finite differences techniques that can handle situations with arbitrary low magnetic diffusivity and it is suitable for studying the formation of sharp current sheets during the evolution. The code is thoroughly tested in different limits and used to illustrate the evolution of the crustal magnetic field in a neutron star in some representative cases. Our code, coupled to cooling codes, can be used to perform long-term simulations of the magneto-thermal evolution of neutron stars.

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 outburst decay of the low magnetic field magnetar SGR 0418+5729

We report on the long-term X-ray monitoring of the outburst decay of the low magnetic field magnetar SGR 0418+5729 using all the available X-ray data obtained with RXTE, Swift, Chandra, and XMM-Newton observations from the discovery of the source in 2009 June up to 2012 August. The timing analysis allowed us to obtain the first measurement of the period derivative of SGR 0418+5729: ˙ P = 4(1) × 10−15 s s−1, significant at a ∼3.5σ confidence level. This leads to a surface dipolar magnetic field of Bdip 6 × 1012 G. This measurement confirms SGR 0418+5729 as the lowest magnetic field magnetar. Following the flux and spectral evolution from the beginning of the outburst up to ∼1200 days, we observe a gradual cooling of the tiny hot spot responsible for the X-ray emission, from a temperature of ∼0.9 to 0.3 keV. Simultaneously, the X-ray flux decreased by about three orders of magnitude: from about 1.4 × 10−11 to 1.2 × 10−14 erg s−1 cm−2. Deep radio, millimeter, optical, and gamma-ray observations did not detect the source counterpart, implying stringent limits on its multi-band emission, as well as constraints on the presence of a fossil disk. By modeling the magneto-thermal secular evolution of SGR 0418+5729, we infer a realistic age of ∼550 kyr, and a dipolar magnetic field at birth of ∼1014 G. The outburst characteristics suggest the presence of a thin twisted bundle with a small heated spot at its base. The bundle untwisted in the first few months following the outburst, while the hot spot decreases in temperature and size. We estimate the outburst rate of low magnetic field magnetars to be about one per year per galaxy, and we briefly discuss the consequences of such a result in several other astrophysical contexts.

A strongly magnetized pulsar within the grasp of the milky way’s supermassive black hole

The center of our Galaxy hosts a supermassive black hole, Sagittarius (Sgr) A∗. Young, massive stars within 0.5 pc of Sgr A∗ are evidence of an episode of intense star formation near the black hole a few million years ago, which might have left behind a young neutron star traveling deep into Sgr A∗’s gravitational potential. On 2013 April 25, a short X-ray burst was observed from the direction of the Galactic center. With a series of observations with the Chandra and the Swift satellites, we pinpoint the associated magnetar at an angular distance of 2.4±0.3 arcsec from Sgr A∗, and refine the source spin period and its derivative (P = 3.7635537(2) s and ˙ P = 6.61(4) × 10−12 s s−1), confirmed by quasi simultaneous radio observations performed with the Green Bank Telescope and Parkes Radio Telescope, which also constrain a dispersion measure of DM = 1750 ± 50 pc cm−3, the highest ever observed for a radio pulsar. We have found that this X-ray source is a young magnetar at ≈0.07–2 pc from Sgr A∗. Simulations of its possible motion around Sgr A∗ show that it is likely (∼90% probability) in a bound orbit around the black hole. The radiation front produced by the past activity from the magnetar passing through the molecular clouds surrounding the Galactic center region might be responsible for a large fraction of the light echoes observed in the Fe fluorescence features.

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.

The influence of magnetic field geometry on magnetars X-ray spectra

Nowadays, the analysis of the X-ray spectra of magnetically powered neutron stars or magnetars is one of the most valuable tools to gain insight into the physical processes occurring in their interiors and magnetospheres. In particular, the magnetospheric plasma leaves a strong imprint on the observed X-ray spectrum by means of Compton up-scattering of the thermal radiation coming from the star surface. Motivated by the increased quality of the observational data, much theoretical work has been devoted to develop Monte Carlo (MC) codes that incorporate the effects of resonant Compton scattering (RCS) in the modeling of radiative transfer of photons through the magnetosphere. The two key ingredients in this simulations are the kinetic plasma properties and the magnetic field (MF) configuration. The MF geometry is expected to be complex, but up to now only mathematically simple solutions (self-similar solutions) have been employed. In this work, we discuss the effects of new, more realistic, MF geometries on synthetic spectra. We use new force-free solutions [14] in a previously developed MC code [9] to assess the influence of MF geometry on the emerging spectra. Our main result is that the shape of the final spectrum is mostly sensitive to uncertain parameters of the magnetospheric plasma, but the MF geometry plays an important role on the angle-dependence of the spectra.

GRO J1008−57: an (almost) predictable transient X-ray binary

A study of archival RXTE, Swift, and Suzaku pointed observations of the transient high-mass X-ray binary GRO J1008−57 is presented. A new orbital ephemeris based on pulse arrival-timing shows the times of maximum luminosities during outbursts of GRO J1008−57 to be close to periastron at orbital phase − 0.03. This makes the source one of a few for which outburst dates can be predicted with very high precision. Spectra of the source in 2005, 2007, and 2011 can be well described by a simple power law with high-energy cutoff and an additional black body at lower energies. The photon index of the power law and the black-body flux only depend on the 15–50 keV source flux. No apparent hysteresis effects are seen. These correlations allow us to predict the evolution of the pulsar’s X-ray spectral shape over all outbursts as a function of just one parameter, the source’s flux. If modified by an additional soft component, this prediction even holds during GRO J1008−57’s 2012 type II outburst.

No anticorrelation between cyclotron line energy and X-ray flux in 4U 0115+634

We report on an outburst of the high mass X-ray binary 4U 0115+634 with a pulse period of 3.6 s in 2008 March/April as observed with RXTE and INTEGRAL. During the outburst the neutron star’s luminosity varied by a factor of 10 in the 3–50 keV band. In agreement with earlier work we find evidence of five cyclotron resonance scattering features at ~10.7, 21.8, 35.5, 46.7, and 59.7 keV. Previous work had found an anticorrelation between the fundamental cyclotron line energy and the X-ray flux. We show that this apparent anticorrelation is probably due to the unphysical interplay of parameters of the cyclotron line with the continuum models used previously, e.g., the negative and positive exponent power law (NPEX). For this model, we show that cyclotron line modeling erroneously leads to describing part of the exponential cutoff and the continuum variability, and not the cyclotron lines. When the X-ray continuum is modeled with a simple exponentially cutoff power law modified by a Gaussian emission feature around 10 keV, the correlation between the line energy and the flux vanishes, and the line parameters remain virtually constant over the outburst. We therefore conclude that the previously reported anticorrelation is an artifact of the assumptions adopted in the modeling of the continuum.

Hot thermal X-ray emission from the Be star HD 119682

We present an analysis of a series of four consecutive Chandra high-resolution transmission gratings observations, amounting to a total of 150 ks, of the Be X-ray source HD 119682 (=1WGA J1346.5–6255), a member of the new class of γ Cas analogs. The Chandra light curve shows significant brightness variations on timescales of hours. However, the spectral distribution appears rather stable within each observation and during the whole campaign. A detailed analysis is not able to detect any coherent pulsation up to a frequency of 0.05 Hz. The Chandra High Energy Transmission Gratings spectrum seems to be devoid of any strong emission line, including Fe Kα fluorescence. The continuum is well described with the addition of two collisionally ionized plasmas of temperatures kT ≈ 15 keV and 0.2 keV, respectively, by the apec model. Models using photoionized plasma components (mekal) or non-thermal components (powerlaw) give poorer fits, providing support for the pure thermal scenario. These two components are absorbed by a single column with N H = (0.20+0.15 –0.03) × 1022 cm–2 compatible with the interstellar value. We conclude that HD 119682 can be regarded as a pole-on γ Cas analog.

Chandra and Suzaku observations of the Be/X-ray star HD110432

We present an analysis of a pointed 141 ks Chandra high-resolution transmission gratings observation of the Be X-ray emitting star HD110432, a prominent member of the γ Cas analogs. This observation represents the first high-resolution spectrum taken for this source as well as the longest uninterrupted observation of any γ Cas analog. The Chandra light curve shows a high variability but its analysis fails to detect any coherent periodicity up to a frequency of 0.05 Hz. Hardness ratio versus intensity analyses demonstrate that the relative contributions of the [1.5-3] Å, [3-6] Å, and [6-16] Å energy bands to the total flux change rapidly in the short term. The analysis of the Chandra High Energy Transmission Grating (HETG) spectrum shows that, to correctly describe the spectrum, three model components are needed. Two of those components are optically thin thermal plasmas of different temperatures (kT ≈ 8-9 and 0.2-0.3 keV, respectively) described by the models vmekal or bvapec. The Fe abundance in each of these two components appears equal within the errors and is slightly subsolar with Z ≈ 0.75 Z ☉. The bvapec model better describes the Fe L transitions, although it cannot fit well the Na XI Lyα line at 10.02 Å, which appears to be overabundant. Two different models seem to describe well the third component. One possibility is a third hot optically thin thermal plasma at kT = 16-21 keV with an Fe abundance Z ≈ 0.3 Z ☉, definitely smaller than for the other two thermal components. Furthermore, the bvapec model describes well the Fe K shell transitions because it accounts for the turbulence broadening of the Fe XXV and Fe XXVI lines with a v turb ≈ 1200 km s–1. These two lines, contributed mainly by the hot thermal plasma, are significantly wider than the Fe Kα line whose FWHM < 5 mÅ is not resolved by Chandra. Alternatively, the third component can be described by a power law with a photon index of Γ = 1.56. In either case, the Chandra HETG spectrum establishes that each one of these components must be modified by distinct absorption columns. The analysis of a noncontemporaneous 25 ks Suzaku observation shows the presence of a hard tail extending up to at least 33 keV. The Suzaku spectrum is described with the sum of two components: an optically thin thermal plasma at kT ≈ 9 keV and Z ≈ 0.74 Z ☉, and a very hot second plasma with kT ≈ 33 keV or, alternatively, a power law with photon index of Γ = 1.58. In either case, each one of the two components must be affected by different absorption columns. Therefore, the kT = 8-9 keV component is definitely needed while the nature of the harder emission cannot be unambiguously established with the present data sets. The analysis of the Si XIII and S XV He-like triplets present in the Chandra spectrum points to a very dense (ne ~ 1013 cm–3) plasma located either close to the stellar surface (r < 3R *) of the Be star or, alternatively, very close (r ~ 1.5R WD) to the surface of a (hypothetical) white dwarf companion. We argue, however, that the available data support the first scenario.

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.

A highly resistive layer within the crust of X-ray pulsars limits their spin periods

The lack of isolated X-ray pulsars with spin periods longer than 12 s raises the question of where the population of evolved high-magnetic-field neutron stars has gone. Unlike canonical radiopulsars, X-ray pulsars are not subject to physical limits to the emission mechanism nor observational biases against the detection of sources with longer periods. Here we show that a highly resistive layer in the innermost part of the crust of neutron stars naturally limits the spin period to a maximum value of about 10–20 s. This highly resistive layer is expected if the inner crust is amorphous and heterogeneous in nuclear charge, possibly owing to the existence of a nuclear ‘pasta’ phase. Our findings suggest that the maximum period of isolated X-ray pulsars may be the first observational evidence for an amorphous inner crust, whose properties can be further constrained by future X-ray timing missions combined with more detailed models.

A new low magnetic field magnetar: the 2011 outburst of Swift J1822.3–1606

We report on the long-term X-ray monitoring with Swift, RXTE, Suzaku, Chandra, and XMM-Newton of the outburst of the newly discovered magnetar Swift J1822.3–1606 (SGR 1822–1606), from the first observations soon after the detection of the short X-ray bursts which led to its discovery, through the first stages of its outburst decay (covering the time span from 2011 July until the end of 2012 April). We also report on archival ROSAT observations which detected the source during its likely quiescent state, and on upper limits on Swift J1822.3–1606's radio-pulsed and optical emission during outburst, with the Green Bank Telescope and the Gran Telescopio Canarias, respectively. Our X-ray timing analysis finds the source rotating with a period of P = 8.43772016(2) s and a period derivative P = 8.3(2)×10−14 s s−1, which implies an inferred dipolar surface magnetic field of B sime 2.7 × 1013 G at the equator. This measurement makes Swift J1822.3–1606 the second lowest magnetic field magnetar (after SGR 0418+5729). Following the flux and spectral evolution from the beginning of the outburst, we find that the flux decreased by about an order of magnitude, with a subtle softening of the spectrum, both typical of the outburst decay of magnetars. By modeling the secular thermal evolution of Swift J1822.3–1606, we find that the observed timing properties of the source, as well as its quiescent X-ray luminosity, can be reproduced if it was born with a poloidal and crustal toroidal fields of Bp ~ 1.5 × 1014 G and B tor ~ 7 × 1014 G, respectively, and if its current age is ~550 kyr.

Modeling magnetar outbursts: flux enhancements and the connection with short bursts and glitches

The availability of a large amount of observational data recently collected from magnetar outbursts is now calling for a complete theoretical study of outburst characteristics. In this Letter (the first of a series dedicated to modeling magnetar outbursts), we tackle the long-standing open issue of whether or not short bursts and glitches are always connected to long-term radiative outbursts. We show that the recent detection of short bursts and glitches seemingly unconnected to outbursts is only misleading our understanding of these events. We show that, in the framework of the starquake model, neutrino emission processes in the magnetar crust limit the temperature, and therefore the luminosity. This natural limit to the maximum luminosity makes outbursts associated with bright persistent magnetars barely detectable. These events are simply seen as a small luminosity increase over the already bright quiescent state, followed by a fast return to quiescence. In particular, this is the case for 1RXS J1708–4009, 1E 1841–045, SGR 1806–20, and other bright persistent magnetars. On the other hand, a similar event (with the same energetics) in a fainter source will drive a more extreme luminosity variation and longer cooling time, as for sources such as XTE J1810–197, 1E 1547–5408, and SGR 1627–41. We conclude that the non-detection of large radiative outbursts in connection with glitches and bursts from bright persistent magnetars is not surprising per se, nor does it need any revision of the glitches and burst mechanisms as explained by current theoretical models.