On a re-examination of neutron star cooling in transient sources -- No shallow heating required?

Abstract

Context: For the typical modeling of neutron stars cooling after an accretion episode in Low-Mass X-ray Binaries, an extra heating source of unknown physical origin, \textit{the shallow heating}, is invoked in order to account for the inferred high effective temperatures of the star up to hundreds of days after the end of accretion. The amount of the shallow heating generated in the crust is usually taken to be proportional to the accretion rate, although the proportionality constant may change from source to source. Aims: In this paper, we intend to model the effective temperature data of eight outburst episodes from seven different sources (\mxb, \xte, \exo, \igr, \swift, \rxs\ and \ks) without {\it ad hoc} shallow heating but accounting for the presence of thermonuclear heating due to the burning of the accreted H/He. Methods: We employ the fully relativistic code \texttt{nscool}, which simulates both the crust and core of a neutron star, equipped with a new boundary condition at the envelope/crust transition which considers thermonuclear heating leakage from the envelope into the crust and depends on the mass accretion rate. Results: We find that the neutron star cooling for seven out of eight of these outbursts can be well explained with this new boundary condition and without the requirement of {\it ad hoc} shallow heating. While the qualitative features of \exo\ cooling curve can be explained, a good fitting still requires additional physics.