Daniele Viganò

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 radiopulsars, 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

N.R. is supported by a Ramon y Cajal Research Fellowship, and by grants AYA2009-07391, AYA2012-39303, SGR2009-811, TW2010005, and iLINK 2011-0303. J.A.P. and D.V. acknowledge support from the grants AYA 2010-21097-C03-02 and Prometeo/2009/103. R.T. and S.M. are partially funded through an INAF 2011 PRIN grant. A.P. is supported by a JAE-Doc CSIC fellowship co-funded with the European Social Fund under the program “Junta para la Ampliacion de Estudios,” by the Spanish MICINN grant AYA2011-30228-C03-02 (co-funded with FEDER funds), and by the AGAUR grant 2009SGR1172 (Catalonia).

A Strongly Magnetized Pulsar within the Grasp of the Milky Way's Supermassive Black Hole

We acknowledge support by grants AYA 2012-39303, SGR2009-811, iLINK 2011-0303, AYA 2010-21097-C03-02, Prometeo 2009/103, AYA2010-17631, P08-TIC-4075, INAF 2010 PRIN grant, Chandra Awards GO2-13076X, G03-14060X, GO3-14099X and G03-14121X, and an EU Marie Curie IEF (FP7-PEOPLE-2012-IEF-331095).

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.

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

A pulsar is a rotating neutron star that beams out electromagnetic waves. The absence of isolated X-ray pulsars with periods longer than 12 s could be a clue to the structural composition of a neutron star’s crust, as simulations show that an amorphous layer would prevent a pulsar from spinning down.

CCOs and the hidden magnetic field scenario

CCOs are X-ray sources lying close the center of supernova remnants, with inferred values of the surface magnetic fields significantly lower (less than about 1e11 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 reemergence of the magnetic field in the crust of a neutron star. A post-supernova accretion stage of about 1e-4-1e-3 solar masses over a vast region of the surface is required to bury the magnetic field into the inner crust. When accretion stops, the field reemerges on a typical timescale 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 alternative to the anti-magnetar scenario, and it could provide the missing link between CCOs and the other classes of isolated neutron stars.

3XMM J185246.6+003317: ANOTHER LOW MAGNETIC FIELD MAGNETAR

N.R. is supported by a Ramon y Cajal fellowship and by an NWO Vidi award. N.R., D.V., and D.F.T. acknowledge support from grants AYA2012-39303, SGR 2009-811, and iLINK 2011-0303. J.A.P. acknowledges support from the grants AYA 2010-21097-C03-02 and Prometeo/2009/103.

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 timing irregularities and the imprint of magnetic field evolution

Context. The rotational evolution of isolated neutron stars is dominated by the magnetic field anchored to the solid crust of the star. Assuming that the core field evolves on much longer timescales, the crustal field evolves mainly though Ohmic dissipation and the Hall drift, and it may be subject to relatively rapid changes with remarkable effects on the observed timing properties. Aims. We investigate whether changes of the magnetic field structure and strength during the star evolution may have observable consequences in the braking index n. This is the most sensitive quantity to reflect small variations of the timing properties that are caused by magnetic field rearrangements. Methods. We performed axisymmetric, long-term simulations of the magneto-thermal evolution of neutron stars with state-of-the-art microphysical inputs to calculate the evolution of the braking index. Relatively rapid magnetic field modifications can be expected only in the crust of neutron stars, where we focus our study. Results. We find that the effect of the magnetic field evolution on the braking index can be divided into three qualitatively different stages depending on the age and the internal temperature: a first stage that may be different for standard pulsars (with n ∼ 3) or low field neutron stars that accreted fallback matter during the supernova explosion (systematically n 3 is expected; in the third stage, at late times, when the interior temperature has dropped to very low values, Hall oscillatory modes in the neutron star crust result in braking indices of a high absolute value and both positive and negative signs. Conclusions. Current magneto-thermal evolution models predict a large contribution to the timing noise and, in particular, to the braking index, from temporal variations of the magnetic field. Models with strong (10 14 G) multipolar or toroidal components, even with a weak (∼10 12 G) dipolar field are consistent with the observed trend of the timing properties.

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 hlogB0[G]i ≈13.0−13.2 with width �(logB0) = 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).