On the mass distribution of neutron stars

The distribution of masses for neutron stars is analysed using the Bayesian statistical inference, evaluating the likelihood of the proposed Gaussian peaks by using 54 measured points obtained in a variety of systems. The results strongly suggest the existence of a bimodal distribution of the masses, with the first peak around 1.37 M(circle dot) and a much wider second peak at 1.73 M(circle dot). The results support earlier views related to the different evolutionary histories of the members for the first two peaks, which produces a natural separation (even if no attempt to `label` the systems has been made here). They also accommodate the recent findings of similar to M(circle dot) masses quite naturally. Finally, we explore the existence of a subgroup around 1.25 M(circle dot), finding weak, if any, evidence for it. This recently claimed low-mass subgroup, possibly related to the O-Mg-Ne core collapse events, has a monotonically decreasing likelihood and does not stand out clearly from the rest of the sample.

Evolution of magnetic field and spin period in accreting neutron stars

Based on the accretion-induced magnetic field decay model, in which a frozen field and an incompressible fluid are assumed, we obtain the following results: (1) an analytic relation between the magnetic field and spin period, if the fastness parameter of the accretion disk is neglected (The evolutionary tracks of accreting neutron stars in the P-B diagram in our model are different from the equilibrium period lines when the influence of the fastness parameter is taken into account.); (2) the theoretical minimum spin period of an accreting neutron star is max(1.1ms (DeltaM/M(circle dot))(-1)R(6)(-5/14) I(45)(M/M(circle dot))(-1/2),1.1ms (M/M(circle dot))(-1/2) R(6)(17/14)), independent of the accretion rate (X-ray luminosity) but dependent on the total accretion mass, DeltaM; however, the minimum magnetic field depends on the accretion rate; (3) the magnetic field strength decreases faster with time than does the period.

Force-free twisted magnetospheres of neutron stars

Context. 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 magnetosphere, which has motivated an increasing number of efforts to improve and generalize existing magnetosphere models. Aims. 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. Methods. 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. Results. 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. Conclusions. 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.

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 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 imprint of the crustal magnetic field on the thermal spectra and pulse profiles of isolated neutron stars

Isolated neutron stars (NSs) show a bewildering variety of astrophysical manifestations, presumably shaped by the magnetic field strength and topology at birth. Here, using state-of-the-art calculations of the coupled magnetic and thermal evolution of NSs, we compute the thermal spectra and pulse profiles expected for a variety of initial magnetic field configurations. In particular, we contrast models with purely poloidal magnetic fields to models dominated by a strong internal toroidal component. We find that, while the former displays double-peaked profiles and very low pulsed fractions, in the latter, the anisotropy in the surface temperature produced by the toroidal field often results in a single pulse profile, with pulsed fractions that can exceed the 50–60 per cent level even for perfectly isotropic local emission. We further use our theoretical results to generate simulated ‘observed’ spectra, and show that blackbody (BB) fits result in inferred radii that can be significantly smaller than the actual NS radius, even as low as ∼1–2 km for old NSs with strong internal toroidal fields and a high absorption column density along their line of sight. We compute the size of the inferred BB radius for a few representative magnetic field configurations, NS ages and magnitudes of the column density. Our theoretical results are of direct relevance to the interpretation of X-ray observations of isolated NSs, as well as to the constraints on the equation of state of dense matter through radius measurements.

Spectral features in isolated neutron stars induced by inhomogeneous surface temperatures

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, magnetic 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.

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 〈log B0[G]〉 ≈ 13.0–13.2 with width σ(log B0) = 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).

The many lives of magnetized neutron stars

The magnetic field strength at birth is arguably one of the most important properties to determine the evolutionary path of a neutron star. Objects with very high fields, collectively known as magnetars, are characterized by high X-ray quiescent luminosities, occurrence of outbursts, and, for some of them, sporadic giant flares. While the magnetic field strength is believed to drive their collective behaviour, however, the diversity of their properties, and, especially, the observation of magnetar-like bursts from “low-field” pulsars, has been a theoretical puzzle. In this review, we discuss results of long-term simulations following the coupled evolution of the X-ray luminosity and the timing properties for a large, homogeneous sample of X-ray emitting isolated neutron stars, accounting for a range of initial magnetic field strengths, envelope compositions, and neutron star masses. In addition, by following the evolution of magnetic stresses within the neutron star crust, we can also relate the observed magnetar phenomenology to the physical properties of neutron stars, and in particular to their age and magnetic field strength and topology. The dichotomy of “high-B” field pulsars versus magnetars is naturally explained, and occasional outbursts from old, low B-field neutron stars are predicted. We conclude by speculating on the fate of old magnetars, and by presenting observational diagnostics of the neutron star crustal field topology.

Quiescent thermal emission from neutron stars in low-mass X-ray binaries

Context. We monitored the quiescent thermal emission from neutron stars in low-mass X-ray binaries after active periods of intense activity in X-rays (outbursts). Aims. The theoretical modeling of the thermal relaxation of the neutron star crust may be used to establish constraints on the crust composition and transport properties, depending on the astrophysical scenarios assumed. Methods. We numerically simulated the thermal evolution of the neutron star crust and compared them with inferred surface temperatures for five sources: MXB 1659−29, KS 1731−260, XTE J1701−462, EXO 0748−676  and IGR J17480−2446. Results. We find that the evolution of MXB 1659−29, KS 1731−260 and EXO 0748−676 can be well described within a deep crustal cooling scenario. Conversely, we find that the other two sources can only be explained with models beyond crustal cooling. For the peculiar emission of XTE J1701−462 we propose alternative scenarios such as residual accretion during quiescence, additional heat sources in the outer crust, and/or thermal isolation of the inner crust due to a buried magnetic field. We also explain the very recent reported temperature of IGR J17480−2446 with an additional heat deposition in the outer crust from shallow sources.

Population synthesis of isolated neutron stars with magneto-rotational evolution – II. From radio-pulsars to magnetars

Population synthesis studies constitute a powerful method to reconstruct the birth distribution of periods and magnetic fields of the pulsar population. When this method is applied to populations in different wavelengths, it can break the degeneracy in the inferred properties of initial distributions that arises from single-band studies. In this context, we extend previous works to include X-ray thermal emitting pulsars within the same evolutionary model as radio-pulsars. We find that the cumulative distribution of the number of X-ray pulsars can be well reproduced by several models that, simultaneously, reproduce the characteristics of the radio-pulsar distribution. However, even considering the most favourable magneto-thermal evolution models with fast field decay, lognormal distributions of the initial magnetic field overpredict the number of visible sources with periods longer than 12 s. We then show that the problem can be solved with different distributions of magnetic field, such as a truncated lognormal distribution, or a binormal distribution with two distinct populations. We use the observational lack of isolated neutron stars (NSs) with spin periods P > 12 s to establish an upper limit to the fraction of magnetars born with B > 1015 G (less than 1 per cent). As future detections keep increasing the magnetar and high-B pulsar statistics, our approach can be used to establish a severe constraint on the maximum magnetic field at birth of NSs.

Neutron-proton bremsstrahlung from intermediate energy heavy-ion reactions as a probe of the nuclear symmetry energy?

Hard photons from neutron-proton bremsstrahlung in intermediate energy heavy-ion reactions are examined as a potential probe of the nuclear symmetry energy within a transport model. Effects of the symmetry energy on the yields and spectra of hard photons are found to be generally smaller than those due to the currently existing uncertainties of both the in-medium nucleon-nucleon cross sections and the photon production probability in the elementary process pn -> pn gamma. Very interestingly, nevertheless, the ratio of hard photon spectra R-1/2(gamma) from two reactions using isotopes of the same element is not only approximately independent of these uncertainties but also quite sensitive to the symmetry energy. For the head-on reactions of Sn-132 + Sn-124 and Sn-112 + Sn-112 at E-beam/A = 50 MeV, for example, the R-1/2(gamma) displays a rise up to 15% when the symmetry energy is reduced by about 20% at rho = 1.3 rho(0) which is the maximum density reached in these reactions. (C) 2008 Elsevier B.V. All rights reserved.

The peculiar isolated neutron star in the Carina Nebula Deep XMM-Newton and ESO-VLT observations of 2XMM J104608.7-594306

While fewer in number than the dominant rotation-powered radio pulsar population, peculiar classes of isolated neutron stars (INSs) which include magnetars, the ROSAT-discovered "Magnificent Seven" (M7), rotating radio transients (RRATs), and central compact objects in supernova remnants (CCOs) - represent a key element in understanding the neutron star phenomenology. We report the results of an observational campaign to study the properties of the source 2XMM J104608.7-594306, a newly discovered thermally emitting INS. The evolutionary state of the neutron star is investigated by means of deep dedicated observations obtained with the XMM-Newton Observatory, the ESO Very Large Telescope, as well as publicly available gamma-ray data from the Fermi Space Telescope and the AGILE Mission. The observations confirm previous expectations and reveal a unique type of object. The source, which is likely within the Carina Nebula (N-H = 2.6x10(21) cm(-2)), has a spectrum that is both thermal and soft, with kT(infinity) = 135 eV. Non-thermal (magnetospheric) emission is not detected down to 1% (3 sigma, 0.1-12 keV) of the source luminosity. Significant deviations (absorption features) from a simple blackbody model are identified in the spectrum of the source around energies 0.6 keV and 1.35 keV. While the former deviation is likely related to a local oxygen overabundance in the Carina Nebula, the latter can only be accounted for by an additional spectral component, which is modelled as a Gaussian line in absorption with EW = 91 eV and sigma = 0.14 keV (1 sigma). Furthermore, the optical counterpart is fainter than m(V) = 27 (2 sigma) and no gamma-ray emission is significantly detected by either the Fermi or AGILE missions. Very interestingly, while these characteristics are remarkably similar to those of the M7 or the only RRAT so far detected in X-rays, which all have spin periods of a few seconds, we found intriguing evidence of very rapid rotation, P = 18.6ms, at the 4 sigma confidence level. We interpret these new results in the light of the observed properties of the currently known neutron star population, in particular those of standard rotation-powered pulsars, recycled objects, and CCOs. We find that none of these scenarios can satisfactorily explain the collective properties of 2XMM J104608.7-594306, although it may be related to the still poorly known class of Galactic anti-magnetars. Future XMM-Newton data, granted for the next cycle of observations (AO11), will help us to improve our current observational interpretation of the source, enabling us to significantly constrain the rate of pulsar spin down.

Magnetic color-flavor-locked stars

Using the solutions of the gap equations of the magnetic-color-flavor-locked (MCFL) phase of paired quark matter in a magnetic field, and taking into consideration the separation between the longitudinal and transverse pressures due to the field-induced breaking of the spatial rotational symmetry, the equation of state (EoS) of the MCFL phase is self-consistently determined. Implications for stellar models of magnetized (self-bound) strange stars and hybrid (MCFL core) stars are discussed.

The masses of the neutron and donor star in the high-mass X-ray binary IGR J18027-2016

Aims. We report near-infrared observations of the supergiant donor to the eclipsing high mass X-ray binary pulsar IGR J18027-2016. We aim to determine its spectral type and measure its radial velocity curve and hence determine the stellar masses of the components. Methods. ESO/VLT observations of the donor utilising the NIR spectrograph ISAAC were obtained in the H and K bands. The multi-epoch H band spectra were cross-correlated with RV templates in order to determine a radial solution for the system. Results. The spectral type of the donor was confirmed as B0-1 I. The radial velocity curve constructed has a semi-amplitude of 23.8 ± 3.1 km s-1. Combined with other measured system parameters, a dynamically determined neutron star mass of 1.4  ±  0.2–1.6  ±  0.3 M⊙ is found. The mass range of the B0-B1 I donor was 18.6  ±  0.8–21.8  ±  2.4 M⊙. These lower and upper limits were obtained under the assumption that the system is viewed edge-on (i = 90° with β = 0.89) for the lower limit and the donor fills its Roche lobe (β = 1 with i = 73.1°) for the upper limit respectively.

Simulated magnetic field expulsion in neutron star cores

The study of long-term evolution of neutron star (NS) magnetic fields is key to understanding the rich diversity of NS observations, and to unifying their nature despite the different emission mechanisms and observed properties. Such studies in principle permit a deeper understanding of the most important parameters driving their apparent variety, e.g. radio pulsars, magnetars, X-ray dim isolated NSs, gamma-ray pulsars. We describe, for the first time, the results from self-consistent magnetothermal simulations considering not only the effects of the Hall-driven field dissipation in the crust, but also adding a complete set of proposed driving forces in a superconducting core. We emphasize how each of these core-field processes drive magnetic evolution and affect observables, and show that when all forces are considered together in vectorial form, the net expulsion of core magnetic flux is negligible, and will have no observable effect in the crust (consequently in the observed surface emission) on megayear time-scales. Our new simulations suggest that strong magnetic fields in NS cores (and the signatures on the NS surface) will persist long after the crustal magnetic field has evolved and decayed, due to the weak combined effects of dissipation and expulsion in the stellar core.