Variability in the Thermal Emission from Accreting Neutron Star   Transients
  The composition of the outer 100 m of a neutron star sets the heat flux that flows outwards from the core. For an accreting neutron star in an X-ray transient, the thermal quiescent flux depends sensitively on the amount of hydrogen and helium remaining on the surface after an accretion outburst and on the composition of the underlying ashes of previous H/He burning. Because H/He has a higher thermal conductivity, a larger mass of H/He implies a shallower thermal gradient through the low density envelope and hence a higher effective temperature for a given core temperature. The mass of residual H and He varies from outburst to outburst, so the thermal quiescent flux is variable even though the core temperature is constant for timescales < 10 000 yr. Heavy elements settle from a H/He envelope in a few hours; we therefore model the quiescent envelope as two distinct layers, H/He over heavier elements, and treat the mass of H/He as a free parameter. We find that the emergent thermal quiescent flux can vary by a factor of 2 to 3 between different quiescent epochs. The variation is more pronounced at lower interior temperatures, making systems with low quiescent luminosities and frequent outbursts, such as SAX J1808.4-3658, ideal candidates from which to observe this effect. We compute, for different ash compositions, the interior temperatures of Cen X-4, Aql X-1, and SAX J1808.4-3658. In the case of Aql X-1, the inferred high interior temperature suggests that neutrino cooling contributes to the neutron star's thermal balance. 
Edward F. Brown (1), Lars Bildsten (2) and Philip Chang (2) ((1)   University of Chicago, (2) University of California, Sant...