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.
One of the most important predictions of any gap model for pulsar magnetospheres is the predicted $gamma$-ray spectra. In the outer gap model, the properties of the synchro-curvature radiation are sensitive to many parameters, whose realistic ranges
have been studied in detail in an accompanying paper. There we demonstrated that the uncertainty in the radius of curvature, the magnetic field geometry, and the X-ray surface flux may affect by orders of magnitude the predicted flux and spectral peak in the $gamma$-ray regime. Here, we present a systematic, numerical study of the impact of the different parameters on the particle dynamics along the gap and calculate the emitted synchro-curvature radiation along the trajectory. By integrating the emitted radiation along the gap and convolving it with a parametrized particle distribution, we discuss how the comparison with the wealth of {em Fermi}-LAT data can be used to constrain the applicability of the model. The resulting spectra show very different energy peaks, fluxes and shapes, qualitatively matching the great variety of the observed {em Fermi}-LAT pulsars. In particular, if we see a large fraction of photons emitted from the initial part of the trajectory, we show that the spectra will be flatter at the low-energy {it Fermi}-LAT regime (100 MeV -- 1 GeV). This provides a solution for such observed flat spectra, while still maintain synchro-curvature radiation as the origin of these photons.
The popular outer gap model of magnetospheric emission from pulsars has been widely applied to explain the properties observed in $gamma$-rays. However, its quantitative predictions rely on a number of approximations and assumptions that are usually
overlooked. Here we examine them, reviewing the main ingredients entering in the model, evaluating their range of uncertainties. Usually, in the quantitative applications of the model, key parameters like the radius of curvature and the energies of the interacting photons are taken to be a fixed, single value. Instead, here we explore their realistic ranges, and the impact of these on the consistency of the model itself. We conclude that the popular evaluation of the trans-field size of the gap as a function of period and period derivative, is unreliable and affected by a huge dispersion. Last, the exploration of the possible values for the radius of curvature, the local magnetic field and other quantities deserve more attention for quantitative applications of the outer gap model, like the calculation of $gamma$-ray spectra, which is the subject of an accompanying paper.
We consider the fundamental problem of charged particles moving along and around a curved magnetic field line, revising the synchro-curvature radiation formulae introduced by Cheng and Zhang (1996). We provide more compact expressions to evaluate the
spectrum emitted by a single particle, identifying the key parameter that controls the transition between the curvature-dominated and the synchrotron-dominated regime. This parameter depends on the local radius of curvature of the magnetic field line, the gyration radius, and the pitch angle. We numerically solve the equations of motion for the emitting particle by considering self-consistently the radiative losses, and provide the radiated spectrum produced by a particle when an electric acceleration is balanced by its radiative losses, as it is assumed to happen in the outer gaps of pulsars magnetospheres. We compute the average spectrum radiated throughout the particle trajectory finding that the slope of the spectrum before the peak depends on the location and size of the emission region. We show how this effect could then lead to a variety of synchro-curvature spectra. Our results reinforce the idea that the purely synchrotron or curvature losses are, in general, inadequate to describe the radiative reaction on the particle motion, and the spectrum of emitted photons. Finally, we discuss the applicability of these calculations to different astrophysical scenarios.
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).
The thermal X-ray spectra of several isolated neutron stars display deviations from a pure blackbody. The accurate physical interpretation of these spectral features bears profound implications for our understanding of the atmospheric composition, ma
gnetic field strength and topology, and equation of state of dense matter. With specific details varying from source to source, common explanations for the features have ranged from atomic transitions in the magnetized atmospheres or condensed surface, to cyclotron lines generated in a hot ionized layer near the surface. Here we quantitatively evaluate the X-ray spectral distortions induced by inhomogeneous temperature distributions of the neutron star surface. To this aim, we explore several surface temperature distributions, we simulate their corresponding general relativistic X-ray spectra (assuming an isotropic, blackbody emission), and fit the latter with a single blackbody model. We find that, in some cases, the presence of a spurious spectral line is required at a high significance level in order to obtain statistically acceptable fits, with central energy and equivalent width similar to the values typically observed. We also perform a fit to a specific object, RX J0806.4-4123, finding several surface temperature distributions able to model the observed spectrum. The explored effect is unlikely to work in all sources with detected lines, but in some cases it can indeed be responsible for the appearance of such lines. Our results enforce the idea that surface temperature anisotropy can be an important factor that should be considered and explored also in combination with more sophisticated emission models like atmospheres.
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.
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.
