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.
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 magnet
osphere 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 $P_0lesssim 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 $langlelog B_0{rm [G]}rangle approx 13.0-13.2$ with width $sigma (log B_0) = 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 $gamma$-rays).
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.
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.
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.
