No Arabic abstract
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 strong magnetic field of neutron stars is intimately coupled to the observed temperature and spectral properties, as well as to the observed timing properties (distribution of spin periods and period derivatives). Thus, a proper theoretical and numerical study of the magnetic field evolution equations, supplemented with detailed calculations of microphysical properties (heat and electrical conductivity, neutrino emission rates) is crucial to understand how the strength and topology of the magnetic field vary as a function of age, which in turn is the key to decipher the physical processes behind the varied neutron star phenomenology. In this review, we go through the basic theory describing the magneto-thermal evolution models of neutron stars, focusing on numerical techniques, and providing a battery of benchmark tests to be used as a reference for present and future code developments. We summarize well-known results from axisymmetric cases, give a new look at the latest 3D advances, and present an overview of the expectations for the field in the coming years.
We study the mutual influence of thermal and magnetic evolution in a neutron stars crust in axial symmetry. Taking into account realistic microphysical inputs, we find the heat released by Joule effect consistent with the circulation of currents in the crust, and we incorporate its effects in 2D cooling calculations. We solve the induction equation numerically using a hybrid method (spectral in angles, but a finite--differences scheme in the radial direction), coupled to the thermal diffusion equation. We present the first long term 2D simulations of the coupled magneto-thermal evolution of neutron stars. This substantially improves previous works in which a very crude approximation in at least one of the parts (thermal or magnetic diffusion) has been adopted. Our results show that the feedback between Joule heating and magnetic diffusion is strong, resulting in a faster dissipation of the stronger fields during the first million years of a NSs life. As a consequence, all neutron stars born with fields larger than a critical value (about 5 10^13 G) reach similar field strengths (approximately 2-3 10^{13} G) at late times. Irrespectively of the initial magnetic field strength, after $10^6$ years the temperature becomes so low that the magnetic diffusion timescale becomes longer than the typical ages of radio--pulsars, thus resulting in apparently no dissipation of the field in old NS. We also confirm the strong correlation between the magnetic field and the surface temperature of relatively young NSs discussed in preliminary works. The effective temperature of models with strong internal toroidal components are systematically higher than those of models with purely poloidal fields, due to the additional energy reservoir stored in the toroidal field that is gradually released as the field dissipates.
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 $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).
The possibility for direct investigation of thermal emission from isolated neutron stars was opened about 25 years ago with the launch of the first X-ray observatory, Einstein. A significant contribution to this study was provided by ROSAT in 1990s. The outstanding capabilities of the currently operating observatories, Chandra and XMM-Newton, have greatly increased the potential to observe and analyze thermal radiation from the neutron star surfaces. Confronting observational data with theoretical models of thermal emission, presumably formed in neutron star atmospheres, allows one to infer the surface temperatures, magnetic fields, chemical composition, and neutron star masses and radii. This information, supplemented with model equations of state and neutron star cooling models, provides an opportunity to understand the fundamental properties of the superdense matter in the neutron star interiors. I review the current status and most important results obtained from modeling neutron star thermal emission and present selected Chandra and XMM-Newton results on thermal radiation from various types of these objects: ordinary radio pulsars with ages ranging from about 2 kyr to 20 Myr (J1119-6127, Vela, B1706-44, J0538+2817, B2334+61, B0656+14, B1055-52, Geminga, B0950+08, J2043+2740), millisecond pulsars (J0030+0451, J2124-3358, J1024-0719, J0437-4715), putative pulsars (CXOU J061705.3+222127, RX J0007.0+7302), central compact objects in supernova remnats (in particular, 1E 1207.4-5209), and isolated radio-quiet neutron stars.
Neutron stars harbour extremely strong magnetic fields within their solid outer crust. The topology of this field strongly influences the surface temperature distribution, and hence the stars observational properties. In this work, we present the first realistic simulations of the coupled crustal magneto-thermal evolution of isolated neutron stars in three dimensions with account for neutrino emission, obtained with the pseudo-spectral code Parody. We investigate both the secular evolution, especially in connection with the onset of instabilities during the Hall phase, and the short-term evolution following episodes of localised energy injection. Simulations show that a resistive tearing instability develops in about a Hall time if the initial toroidal field exceeds ~$10^{15}$ G. This leads to crustal failures because of the huge magnetic stresses coupled with the local temperature enhancement produced by dissipation. Localised heat deposition in the crust results in the appearance of hot spots on the star surface which can exhibit a variety of patterns. Since the transport properties are strongly influenced by the magnetic field, the hot regions tend to drift away and get deformed following the magnetic field lines while cooling. The shapes obtained with our simulations are reminiscent of those recently derived from NICER X-ray observations of the millisecond pulsar PSR J0030+0451.