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.

The role of nuclear physics in supernovae and the evolution of neutron stars, Neutrino Opacities, Equation of State, Transport Coefficients, and Dark Matter Production

Thesis (Ph.D.)--University of Washington, 2016-08

On neutron stars in ƒ (R) theories: Small radii, large masses and large energy emitted in a merger

In the context of ƒ (R) gravity theories, we show that the apparent mass of a neutron star as seen from an observer at infinity is numerically calculable but requires careful matching, first at the star’s edge, between interior and exterior solutions, none of them being totally Schwarzschild-like but presenting instead small oscillations of the curvature scalar R; and second at large radii, where the Newtonian potential is used to identify the mass of the neutron star. We find that for the same equation of state, this mass definition is always larger than its general relativistic counterpart. We exemplify this with quadratic R^2 and Hu-Sawicki-like modifications of the standard General Relativity action. Therefore, the finding of two-solar mass neutron stars basically imposes no constraint on stable ƒ (R) theories. However, star radii are in general smaller than in General Relativity, which can give an observational handle on such classes of models at the astrophysical level. Both larger masses and smaller matter radii are due to much of the apparent effective energy residing in the outer metric for scalar-tensor theories. Finally, because the ƒ (R) neutron star masses can be much larger than General Relativity counterparts, the total energy available for radiating gravitational waves could be of order several solar masses, and thus a merger of these stars constitutes an interesting wave source.

Self-bound models of compact stars and recent mass-radius measurements

The exact composition of a specific class of compact stars, historically referred to as ""neutron stars,'' is still quite unknown. Possibilities ranging from hadronic to quark degrees of freedom, including self-bound versions of the latter, have been proposed. We specifically address the suitability of strange star models (including pairing interactions) in this work, in the light of new measurements available for four compact stars. The analysis shows that these data might be explained by such an exotic equation of state, actually selecting a small window in parameter space, but still new precise measurements and also further theoretical developments are needed to settle the subject.

The isolated neutron star candidate 2XMM J104608.7-594306

Over the last decade, X-ray observations have revealed the existence of several classes of isolated neutron stars (INSs) which are radio-quiet or exhibit radio emission with properties much at variance with those of ordinary radio pulsars. The identification of new sources is crucial in order to understand the relations among the different classes and to compare observational constraints with theoretical expectations. A recent analysis of the 2XMMp catalogue provided fewer than 30 new thermally emitting INS candidates. Among these, the source 2XMM J104608.7-594306 appears particularly interesting because of the softness of its X-ray spectrum, kT = 117 +/- 14 eV and N(H) = (3.5 +/- 1.1) x 10(21) cm(-2) (3 sigma), and of the present upper limits in the optical, m(B) greater than or similar to 26, m(V) greater than or similar to 25.5 and m(R) greater than or similar to 25 (98.76% confidence level), which imply a logarithmic X-ray-to-optical flux ratio log(F(X)/F(V)) greater than or similar to 3.1, corrected for absorption. We present the X-ray and optical properties of 2XMM J104608.7-594306 and discuss its nature in the light of two possible scenarios invoked to explain the X-ray thermal emission from INSs: the release of residual heat in a cooling neutron star, as in the seven radio-quiet ROSAT-discovered INSs, and accretion from the interstellar medium. We find that the present observational picture of 2XMM J104608.7-594306 is consistent with a distant cooling INS with properties in agreement with the most up-to-date expectations of population synthesis models: it is fainter, hotter and more absorbed than the seven ROSAT sources and possibly located in the Carina Nebula, a region likely to harbour unidentified cooling neutron stars. The accretion scenario, although not entirely ruled out by observations, would require a very slow (similar to 10 km s(-1)) INS accreting at the Bondi-Hoyle rate.

Compact Stars and Magnetized CFL Matter

The stability of the color flavor locked phase in the presence of a strong magnetic field is investigated within the phenomenological MIT bag model. It is found that the minimum value of the energy per baryon in a color flavor locked state at vanishing pressure is lower than the corresponding one for unpaired magnetized strange quark matter and, as the magnetic field increases, the energy per baryon decreases. This implies that magnetized color flavor locked matter is more stable and could become the ground state inside neutron stars. The anisotropy of the pressures is discussed. The mass-radius relation for such stars is also studied.

Magnetized color flavor locked state and compact stars

The stability of the color flavor locked phase in the presence of a strong magnetic field is investigated within the phenomenological MIT bag model, taking into account the variation of the strange quark mass, the baryon density, the magnetic field, as well as the bag and gap parameters. It is found that the minimum value of the energy per baryon in a color flavor locked state at vanishing pressure is lower than the corresponding one for unpaired magnetized strange quark matter and, as the magnetic field increases, the energy per baryon decreases. This implies that magnetized colorflavor locked matter is more stable and could become the ground state inside neutron stars. The mass-radius relation for such stars is also studied.

Neutron skin of 208Pb, nuclear symmetry energy, and the parity radius experiment

A precise determination of the neutron skin thickness of a heavy nucleus sets a basic constraint on the nuclear symmetry energy (the neutron skin thickness is the difference of the neutron and proton rms radii of the nucleus). The parity radius experiment (PREX) may achieve it by electroweak parity-violating electron scattering (PVES) on 208Pb. We investigate PVES in nuclear mean field approach to allow the accurate extraction of the neutron skin thickness of 208Pb from the parity-violating asymmetry probed in the experiment. We demonstrate a high linear correlation between the parity-violating asymmetry and the neutron skin thickness in successful mean field forces as the best means to constrain the neutron skin of 208Pb from PREX, without assumptions on the neutron density shape. Continuation of the experiment with higher precision in the parity-violating asymmetry is motivated since the present method can support it to constrain the density slope of the nuclear symmetry energy to new accuracy.

Search for semiclassical gravity effects in relativistic stars

We discuss the possible influence of gravity in the neutronization process p+e-→νe, which is particularly important as a cooling mechanism of neutron stars. Our approach is semiclassical in the sense that leptonic fields are quantized on a classical background spacetime, while neutrons and protons are treated as excited and unexcited nucleon states, respectively. We expect gravity to have some influence wherever the energy content carried by the in state is barely above the neutron mass. In this case the emitted neutrinos would be soft enough to have a wavelength of the same order as the space curvature radius. ©2000 The American Physical Society.

Medium effects in the pion-pole mechanism [gamma gamma ->pi(0)->nu(R)(nu)over-bar(L)(nu(L)(nu)over-bar(R))] of neutron star cooling

Nuclear medium effects in the neutrino cooling of neutron stars through the reaction channel γγ→π0 →ν Rν̄L(νLν̄R) are incorporated. Throughout the paper we discuss different possibilities of right-handed neutrinos, massive left-handed neutrinos, and standard massless left-handed neutrinos (reaction is then allowed only with medium modified vertices). It is demonstrated that multiparticle effects suppress the rate of this reaction channel in the dense hadron matter by 6-7 orders of magnitude that does not allow to decrease existing experimental upper limit on the corresponding π0νν̄ coupling. Other possibilities of the manifestation of the given reaction channel in different physical situations, e.g., in the quark color superconducting cores of the most massive neutron stars, are also discussed. We demonstrate that in the color-flavor-locked superconducting phase for temperatures T≲ 0.1-10 MeV (depending on the effective pion mass and the decay width) the process is feasibly the most efficient neutrino cooling process, although the absolute value of the reaction rate is rather small.

Entropy, complexity and disequilibrium in compact stars

We used the statistical measurements of information entropy, disequilibrium and complexity to infer a hierarchy of equations of state for two types of compact stars from the broad class of neutron stars, namely, with hadronic composition and with strange quark composition. Our results show that, since order costs energy. Nature would favor the exotic strange stars even though the question of how to form the strange stars cannot be answered within this approach. (C) 2012 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.