Isolated neutron stars show a diversity in timing and spectral properties, which has historically led to a classification in different sub-classes. The magnetic field plays a key role in many aspects of the neutron star phenomenology: it regulates th
e braking torque responsible for their timing properties and, for magnetars, it provides the energy budget for the outburst activity and high quiescent luminosities (usually well above the rotational energy budget). We aim at unifying this observational variety by linking the results of the state-of-the-art 2D magneto-thermal simulations with observational data. The comparison between theory and observations allows to place two strong constraints on the physical properties of the inner crust. First, strong electrical currents must circulate in the crust, rather than in the star core. Second, the innermost part of the crust must be highly resistive, which is in principle in agreement with the presence of a novel phase of matter so-called nuclear pasta phase.
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 em
ission. 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.
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.
High magnetic fields are a distinguishing feature of neutron stars and the existence of sources (the soft gamma repeaters and the anomalous X-ray pulsars) hosting an ultra-magnetized neutron star (or magnetar) has been recognized in the past few deca
des. Magnetars are believed to be powered by magnetic energy and not by rotation, as with normal radio pulsars. Until recently, the radio quietness and magnetic fields typically above the quantum critical value (Bq~4.4x10^{13} G), were among the characterizing properties of magnetars. The recent discovery of radio pulsed emission from a few of them, and of a low dipolar magnetic field soft gamma repeater, weakened further the idea of a clean separation between normal pulsars and magnetars. In this Letter we show that radio emission from magnetars might be powered by rotational energy, similarly to what occurs in normal radio pulsars. The peculiar characteristics of magnetars radio emission should be traced in the complex magnetic geometry of these sources. Furthermore, we propose that magnetar radio activity or inactivity can be predicted from the knowledge of the stars rotational period, its time derivative and the quiescent X-ray luminosity.
The X-ray spectra observed in the persistent emission of magnetars are evidence for the existence of a magnetosphere. The high-energy part of the spectra is explained by resonant cyclotron upscattering of soft thermal photons in a twisted magnetosphe
re, which has motivated an increasing number of efforts to improve and generalize existing magnetosphere models. We want to build more general configurations of twisted, force-free magnetospheres as a first step to understanding the role played by the magnetic field geometry in the observed spectra. First we reviewed and extended previous analytical works to assess the viability and limitations of semi-analytical approaches. Second, we built a numerical code able to relax an initial configuration of a nonrotating magnetosphere to a force-free geometry, provided any arbitrary form of the magnetic field at the star surface. The numerical code is based on a finite-difference time-domain, divergence-free, and conservative scheme, based of the magneto-frictional method used in other scenarios. We obtain new numerical configurations of twisted magnetospheres, with distributions of twist and currents that differ from previous analytical solutions. The range of global twist of the new family of solutions is similar to the existing semi-analytical models (up to some radians), but the achieved geometry may be quite different. The geometry of twisted, force-free magnetospheres shows a wider variety of possibilities than previously considered. This has implications for the observed spectra and opens the possibility of implementing alternative models in simulations of radiative transfer aiming at providing spectra to be compared with observations.
We report on a new mechanism for heat conduction in the neutron star crust. We find that collective modes of superfluid neutron matter, called superfluid phonons (sPhs), can influence heat conduction in magnetized neutron stars. They can dominate the
heat conduction transverse to magnetic field when the magnetic field $B gsim 10^{13}$ G. At density $rho simeq 10^{12}-10^{14} $ g/cm$^3$ the conductivity due to sPhs is significantly larger than that due to lattice phonons and is comparable to electron conductivity when temperature $simeq 10^8$ K. This new mode of heat conduction can limit the surface anisotropy in highly magnetized neutron stars. Cooling curves of magnetized neutron stars with and without superfluid heat conduction could show observationally discernible differences.
Context: Many thermally emitting isolated neutron stars have magnetic fields larger than 10^13 G. A realistic cooling model that includes the presence of high magnetic fields should be reconsidered. Aims: We investigate the effects of anisotropic tem
perature distribution and Joule heating on the cooling of magnetized neutron stars. Methods: The 2D heat transfer equation with anisotropic thermal conductivity tensor and including all relevant neutrino emission processes is solved for realistic models of the neutron star interior and crust. Results: The presence of the magnetic field affects significantly the thermal surface distribution and the cooling history during both, the early neutrino cooling era and the late photon cooling era. Conclusions: There is a large effect of the Joule heating on the thermal evolution of strongly magnetized neutron stars. Both magnetic fields and Joule heating play a key role in keeping magnetars warm for a long time. Moreover, this effect is important for intermediate field neutron stars and should be considered in radio-quiet isolated neutron stars or high magnetic field radio-pulsars.
We present two-dimensional simulations for the cooling of neutron stars with strong magnetic fields (B > 1e13 Gauss). We study how the cooling curves are influenced by magnetic field decay. We show that the Joule heating effects are very large and in
some cases control the thermal evolution. We characterize the temperature anisotropy induced by the magnetic field and predict the surface temperature distribution for the early and late stages of the evolution of isolated neutron stars, comparing our results with available observational data of isolated neutron stars.
The impact of strong magnetic fields B>10e13 G on the thermal evolution of neutron stars is investigated, including crustal heating by magnetic field decay. For this purpose, we perform 2D cooling simulations with anisotropic thermal conductivity con
sidering all relevant neutrino emission processes for realistic neutron stars. The standard cooling models of neutron stars are called into question by showing that the magnetic field has relevant (and in many cases dominant) effects on the thermal evolution. The presence of the magnetic field significantly affects the thermal surface distribution and the cooling history of these objects during both, the early neutrino cooling era and the late photon cooling era. The minimal cooling scenario is thus more complex than generally assumed. A consistent magneto-thermal evolution of magnetized neutron stars is needed to explain the observations.
We present 2D simulations of the cooling of neutron stars with strong magnetic fields (B geq 10^{13} G). We solve the diffusion equation in axial symmetry including the state of the art microphysics that controls the cooling such as slow/fast neutrin
o processes, superfluidity, as well as possible heating mechanisms. We study how the cooling curves depend on the the magnetic field strength and geometry. Special attention is given to discuss the influence of magnetic field decay. We show that Joule heating effects are very large and in some cases control the thermal evolution. We characterize the temperature anisotropy induced by the magnetic field for the early and late stages of the evolution of isolated neutron stars.
