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.

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 spindown of a pulsar, or 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 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 year 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 don't overlap, then pulsar spindown 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.

Valuable information about the neutron star interior can be obtained by comparing observations of thermal radiation from a cooling neutron star crust with theoretical models. Nuclear burning of lighter elements that diffuse to deeper layers of the envelope can alter the relation between surface and interior temperatures and can change the chemical composition over time. We calculate new temperature relations and consider two effects of diffusive nuclear burning (DNB) for H-C envelopes. First, we consider the effect of a changing envelope composition and find that hydrogen is consumed on short timescales and our temperature evolution simulations correspond to those of a hydrogen-poor envelope within ~100 days. The transition from a hydrogen-rich to a hydrogen-poor envelope is potentially observable in accreting NS systems as an additional initial decline in surface temperature at early times after the outburst. Second, we find that DNB can produce a non-negligible heat flux, such that the total luminosity can be dominated by DNB in the envelope rather than heat from the deep interior. However, without continual accretion, heating by DNB in H-C envelopes is only relevant for <1-80 days after the end of an accretion outburst, as the amount of light elements is rapidly depleted. Comparison to crust cooling data shows that DNB does not remove the need for an additional shallow heating source. We conclude that solving the time-dependent equations of the burning region in the envelope self-consistently in thermal evolution models instead of using static temperature relations would be valuable in future cooling studies.

A neo-neutron star is a hot neutron star that has just become transparent to neutrinos. In a core collapse supernova or accretion induced collapse of a white dwarf the neo-neutron star phase directly follows the proto-neutron star phase, about 30 to 60 seconds after the initial collapse. It will also be present in a binary neutron star merger in the case the "born-again" hot massive compact star does not immediately collapse into a black hole. Eddington or even super-Eddington luminosities are present for some time. A neo-neutron star produced in a core collapse supernova is not directly observable but the one produced by a binary merger, likely associated with an off-axis short gamma-ray burst, may be observable for some time as well as when produced in the accretion induced collapse of a white dwarf. We present a first step in the study of this neo-neutron star phase in a spherically symmetric configuration, thus neglecting fast rotation, and also neglecting the effect of strong magnetic fields. We put particular emphasis on determining how long the star can sustain a near-Eddington luminosity and also show the importance of positrons and contraction energy during neo-neutron star phase. We finally discuss the observational prospects for neutron star mergers triggered by LIGO and for accretion-induced collapse transients.

The study of how neutron stars cool over time can provide invaluable insights into fundamental physics such as the nuclear equation of state and superconductivity and superfluidity. A critical relation in neutron star cooling is the one between observed surface temperature and interior temperature. This relation is determined by the composition of the neutron star envelope and can be influenced by the process of diffusive nuclear burning (DNB). We calculate models of envelopes that include DNB and find that DNB can lead to a rapidly changing envelope composition which can be relevant for understanding the long-term cooling behavior of neutron stars. We also report on analysis of the latest temperature measurements of the young neutron star in the Cassiopeia A supernova remnant. The 13 Chandra observations over 18 years show that the neutron star's temperature is decreasing at a rate of 2-3 percent per decade, and this rapid cooling can be explained by the presence of a proton superconductor and neutron superfluid in the core of the star.

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 $\times 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 long-term 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 $\approx$ 2%. Finally, we fit the observed cooling trend of the Cassiopeia A neutron star with a model that includes DNB in the envelope.

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 $\ll 10$ keV) axion-like particles to neutrons. Using Markov-Chain Monte Carlo we find that for the values of axion-neutron coupling $g_{ann}^2 > 7.7 \times 10^{-20}$ (90% c.l.) the axion cooling from the bremsstrahlung reaction $n+n\rightarrow 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 $f_a > 6.7 \times 10^7$ GeV for KSVZ axions and $f_a > 1.7 \times 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.

We study the effects of heat blanketing envelopes of neutron stars on their cooling. To this aim, we perform cooling simulations using newly constructed models of the envelopes composed of binary ion mixtures (H--He, He--C, C--Fe) varying the mass of lighter ions (H, He or C) in the envelope. The results are compared with those calculated using the standard models of the envelopes which contain the layers of lighter (accreted) elements (H, He and C) on top of the Fe layer, varying the mass of accreted elements. The main effect is that the chemical composition of the envelopes influences their thermal conductivity and, hence, thermal insulation of the star. For illustration, we apply these results to estimate the internal temperature of the Vela pulsar and to study the cooling of neutron stars of ages of 0.1 - 1 Myr at the photon cooling stage. The uncertainties of the cooling models associated with our poor knowledge of chemical composition of the heat insulating envelopes strongly complicate theoretical reconstruction of the internal structure of cooling neutron stars from observations of their thermal surface emission.

We construct new models of outer heat blanketing envelopes of neutron stars composed of binary ion mixtures (H - He, He - C, C - Fe) in and out of diffusive equilibrium. To this aim, we generalize our previous work on diffusion of ions in isothermal gaseous or Coulomb liquid plasmas to handle non-isothermal systems. We calculate the relations between the effective surface temperature Ts and the temperature Tb at the bottom of heat blanketing envelopes (at a density rhob= 1e8 -- 1e10 g/cc) for diffusively equilibrated and non-equilibrated distributions of ion species at different masses DeltaM of lighter ions in the envelope. Our principal result is that the Ts - Tb relations are fairly insensitive to detailed distribution of ion fractions over the envelope (diffusively equilibrated or not) and depend almost solely on DeltaM. The obtained relations are approximated by analytic expressions which are convenient for modeling the evolution of neutron stars.

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.

We apply our recently suggested statistical approach to thermal evolution of isolated neutron stars and accreting quasistationary neutron stars in X-ray transients for constraining the position and relative broadening alpha of the direct Urca threshold of powerful neutrino emission in neutron star cores. We show that most likely explanation of observations corresponds to alpha = 0.08 - 0.10 and to the neutron star mass, at which the direct Urca process is open, M_D = (1.6 - 1.8) M_sun.

We calculate interdiffusion coefficients in a two-component, weakly or strongly coupled ion plasma (gas or liquid, composed of two ion species immersed into a neutralizing electron background). We use an effective potential method proposed recently by Baalrud and Daligaut [PRL, 110, 235001, (2013)]. It allows us to extend the standard Chapman-Enskog procedure of calculating the interdiffusion coefficients to the case of strong Coulomb coupling. We compute binary diffusion coefficients for several ionic mixtures and fit them by convenient expressions in terms of the generalized Coulomb logarithm. These fits cover a wide range of plasma parameters spanning from weak to strong Coulomb couplings. They can be used to simulate diffusion of ions in ordinary stars as well as in white dwarfs and neutron stars.

We analyze diffusion equations in strongly coupled Coulomb mixtures of ions in dense stellar matter. Strong coupling of ions in the presence of gravitational forces and electric fields (induced by plasma polarization in the presence of gravity) produces a specific diffusion current which can separate ions with the same A/Z (mass to charge number) ratios but different Z. This Coulomb separation of ions can be important for the evolution of white dwarfs and neutron stars.