Accretion heated atmospheres of X-ray bursting neutron stars

Some thermonuclear (type I) X-ray bursts at the neutron star surfaces in low-mass X-ray binaries take place during hard persistent states of the systems. Spectral evolution of these bursts is well described by the atmosphere model of a passively cooling neutron star when the burst luminosity is high enough. The observed spectral evolution deviates from the model predictions when the burst luminosity drops below a critical value of 20-70% of the maximum luminosity. We suggest that these deviations are induced by the additional heating of the accreted particles. We present a method for computation of the neutron star atmosphere models heated by accreted particles assuming that their energy is released via Coulomb interactions with electrons. We compute the temperature structures and the emergent spectra of the atmospheres of various chemical compositions and investigate the dependence of the results on the other model parameters. We show that the heated atmosphere develops the hot (20--100 keV) corona-like surface layer cooled by Compton scattering, and the deeper, almost isothermal optically thick region with a temperature of a few keV. The emergent spectra deviate strongly from those of undisturbed neutron star atmospheres, with the main differences being the presence of a high-energy tail and a strong excess in the low-energy part of the spectrum. They also lack the iron absorption edge, which is visible in the spectra of undisturbed low-luminosity atmospheres with solar chemical composition. Using the computed spectra, we obtained the dependences of the dilution and color-correction factors as functions of relative luminosities for pure helium and solar abundance atmospheres. We show that the helium model atmosphere heated by accretion corresponding to 5% of the Eddington luminosity describes well the late stages of the X-ray bursts in 4U 1820-30.

A window into the neutron star: Modelling the cooling of accretion heated neutron star crusts

In accreting neutron star X-ray transients, the neutron star crust can be substantially heated out of thermal equilibrium with the core during an accretion outburst. The observed subsequent cooling in quiescence (when accretion has halted) offers a unique opportunity to study the structure and thermal properties of the crust. Initially crust cooling modelling studies focussed on transient X-ray binaries with prolonged accretion outbursts (> 1 year) such that the crust would be significantly heated for the cooling to be detectable. Here we present the results of applying a theoretical model to the observed cooling curve after a short accretion outburst of only ~10 weeks. In our study we use the 2010 outburst of the transiently accreting 11 Hz X-ray pulsar in the globular cluster Terzan 5. Observationally it was found that the crust in this source was still hot more than 4 years after the end of its short accretion outburst. From our modelling we found that such a long-lived hot crust implies some unusual crustal properties such as a very low thermal conductivity (> 10 times lower than determined for the other crust cooling sources). In addition, we present our preliminary results of the modelling of the ongoing cooling of the neutron star in MXB 1659-298. This transient X-ray source went back into quiescence in March 2017 after an accretion phase of ~1.8 years. We compare our predictions for the cooling curve after this outburst with the cooling curve of the same source obtained after its previous outburst which ended in 2001.

Potential cooling of an accretion-heated neutron star crust in the low-mass X-ray binary 1RXS J180408.9-342058

We have monitored the transient neutron star low-mass X-ray binary 1RXS J180408.9-342058 in quiescence after its ~4.5 month outburst in 2015. The source has been observed using Swift and XMM-Newton. Its X-ray spectra were dominated by a thermal component. The thermal evolution showed a gradual X-ray luminosity decay from ~18 x 10^32 to ~4 x 10^32 (D/5.8 kpc)^2 erg s^{-1} between ~8 and ~379 d in quiescence, and the inferred neutron star surface temperature (for an observer at infinity; using a neutron star atmosphere model) decreased from ~100 to ~71 eV. This can be interpreted as cooling of an accretion-heated neutron star crust. Modelling the observed temperature curve (using NSCOOL) indicated that the source required ~1.9 MeV per accreted nucleon of shallow heating in addition to the standard deep crustal heating to explain its thermal evolution. Alternatively, the decay could also be modelled without the presence of deep crustal heating, only having a shallow heat source (again ~1.9 MeV per accreted nucleon was required). However, the XMM-Newton data statistically required an additional power-law component. This component contributed ~30 per cent of the total unabsorbed flux in 0.5-10 keV energy range. The physical origin of this component is unknown. One possibility is that it arises from low-level accretion. The presence of this component in the spectrum complicates our cooling crust interpretation because it might indicate that the smooth luminosity and temperature decay curves we observed may not be due to crust cooling but due to some other process.

Neutron stars - cooling and transport

Observations of thermal radiation from neutron stars can potentially provide information about the states of supranuclear matter in the interiors of these stars with the aid of the theory of neutron-star thermal evolution. We review the basics of this theory for isolated neutron stars with strong magnetic fields, including most relevant thermodynamic and kinetic properties in the stellar core, crust, and blanketing envelopes.

Pairing and the Cooling of Neutron Stars

In this review, I present a brief summary of the impact of nucleon pairing at supra-nuclear densities on the cooling of neutron stars. I also describe how the recent observation of the cooling of the neutron star in the supernova remnant Cassiopeia A may provide us with the first direct evidence for the occurrence of such pairing. It also implies a size of the neutron 3P-F2 energy gap of the order of 0.1 MeV.