Thermodynamically consistent accreted crust of neutron stars: The role of proton shell effects

Abstract

Observations of accreting neutron stars are widely used to constrain the microphysical properties of superdense matter. A key ingredient in this analysis is the heating associated with nuclear reactions in the outer layers of the neutron star (crust), as well as the equation of state and composition of these layers. As recently shown, the neutron hydrostatic/diffusion (nHD) condition is valid in the inner part of the crust, where some of the neutrons are not bound to the nuclei, and this condition should be properly incorporated into crustal models. Here we construct models of the accreted crust of a neutron star, taking into account the nHD condition and proton shell effects in nuclei. For numerical illustration, we employ the recently proposed compressible liquid drop model, which incorporates shell effects. However, our approach is general and can also be used in future studies relying on more sophisticated nuclear physics models.