SIMBAD references

We describe a new method to incorporate thermonuclear heating in the envelope of accreting neutron stars into long-term simulations of their thermal evolution. We obtain boundary conditions for the heat exchange between the envelope and the crust based on stationary models, which include nuclear burning, and validate these values by comparing them to the results of the time-dependent code MESA. These simple boundary conditions allow us to explore a large parameter space. We quantify the amount of heat flowing from the envelope into the crust, or vice versa, depending on the mass accretion rate, outburst duration, and duty cycle, and especially crust/core physical parameters such as impurities, crustal heating, and neutrino cooling rate.