If the thermal evolution of the hot young neutron star in the supernova remnant HESS J1731-347 is driven by neutrino emission, it provides a stringent constraint on the coupling of light (mass 10 keV) axion-like particles to neutrons. Using Markov-Chain Monte Carlo we find that for the values of axion-neutron coupling g 2 ann > 7.7 × 10 −20 (90% c.l.) the axion cooling from the bremsstrahlung reaction n + n → n + n + a is too rapid to account for the high observed surface temperature. This implies that the Pecci-Quinn scale or axion decay constant fa > 6.7 × 10 7 GeV for KSVZ axions and fa > 1.7 × 10 9 GeV for DFSZ axions. The high temperature of this neutron star also allows us to tighten constraints on the size of the nucleon pairing gaps.

Ages and thermal luminosities of neutron stars, inferred from observations, can be interpreted with the aid of the neutron star cooling theory to gain information on the properties of superdense matter in neutron-star interiors. We present a survey of estimated ages, surface temperatures, and thermal luminosities of middle-aged neutron stars with relatively weak or moderately strong magnetic fields, which can be useful for these purposes. The catalogue includes results selected from the literature, supplemented with new results of spectral analysis of a few cooling neutron stars. The data are compared with the theory. We show that overall agreement of theoretical cooling curves with observations improves substantially for models where neutron superfluidity in stellar core is weak.

Thermal evolution of neutron stars is known to depend on the properties of superdense matter in neutron star cores. We suggest a statistical analysis of isolated cooling middle-aged neutron stars and old transiently accreting quasi-stationary neutron stars warmed up by deep crustal heating in low-mass X-ray binaries. The method is based on simulations of the evolution of stars of different masses and on averaging the results over respective mass distributions. This gives theoretical distributions of isolated neutron stars in the surface temperature-age plane and of accreting stars in the photon thermal luminosity-mean mass accretion rate plane to be compared with observations. This approach permits to explore not only superdense matter but also the mass distributions of isolated and accreting neutron stars. We show that the observations of these stars can be reasonably well explained by assuming the presence of the powerful direct Urca process of neutrino emission in the inner cores of massive stars, introducing a slight broadening of the direct Urca threshold (for instance, by proton superfluidity), and by tuning mass distributions of isolated and accreted neutron stars.

The possible detection of a compact object in the remnant of SN 1987A presents an unprecedented opportunity to follow its early evolution. The suspected detection stems from an excess of infrared emission from a dust blob near the compact object’s predicted position. The infrared excess could be due to the decay of isotopes like 44Ti, accretion luminosity from a neutron star or black hole, magnetospheric emission or a wind originating from the spin down of a pulsar, or to thermal emission from an embedded, cooling neutron star (NS 1987A). It is shown that the last possibility is the most plausible as the other explanations are disfavored by other observations and/or require fine-tuning of parameters. Not only are there indications that the dust blob overlaps the predicted location of a kicked compact remnant, but its excess luminosity also matches the expected thermal power of a 30 yr old neutron star. Furthermore, models of cooling neutron stars within the minimal cooling paradigm readily fit both NS 1987A and Cas A, the next-youngest known neutron star. If correct, a long heat transport timescale in the crust and a large effective stellar temperature are favored, implying relatively limited crustal n-1S0 superfluidity and an envelope with a thick layer of light elements, respectively. If the locations do not overlap, then pulsar spin down or accretion might be more likely, but the pulsar’s period and magnetic field or the accretion rate must be rather finely tuned. In this case, NS 1987A may have enhanced cooling and/or a heavy-element envelope.

A critical relation in the study of neutron star cooling is the one between surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be affected by the process of diffusive nuclear burning (DNB), which occurs when elements diffuse to depths where the density and temperature are sufficiently high to ignite nuclear burning. We calculate models of H-He and He-C envelopes that include DNB and obtain analytic temperature relations that can be used in neutron star cooling simulations. We find that DNB can lead to a rapidly changing envelope composition and prevents the build-up of thermally stable hydrogen columns y H 10 7 g cm −2 , while DNB can make helium envelopes more transparent to heat flux for surface temperatures T s 2× 10 6 K. We perform neutron star cooling simulations in which we evolve temperature and envelope composition, with the latter due to DNB and accretion from the interstellar medium. We find that a time-dependent envelope composition can be relevant for understanding the longterm cooling behaviour of isolated neutron stars. We also report on the latest Chandra observations of the young neutron star in the Cassiopeia A supernova remnant; the resulting 13 temperature measurements over more than 18 years yield a ten-year cooling rate of ≈ 2%. Finally, we fit the observed cooling trend of the Cassiopeia A neutron star with a model that includes DNB in the envelope.