Time-independent Simulations of Steady-State Accretion with Nuclear Burning

Abstract: We construct a new formulation that allows efficient exploration of steady-state accretion processes onto compact objects. Accretion onto compact objects is a common scenario in astronomy. These systems serve as laboratories to probe the nuclear burning of the accreted matter. Conventional stellar evolution codes have been developed to simulate in detail the nuclear reactions on the compact objects. In order to follow the case of steady burning, however, using these codes can be very expensive as they are designed to follow a time-dependent problem. Here we introduce our new code $\textsc{StarShot}$, which resolves the structure of the compact objects for the case of stable thermonuclear burning, and is able to follow all nuclear species using an adaptive nuclear reaction network and adaptive zoning. Compared to dynamical codes, the governing equations can be reduced to time-independent forms under the assumption of steady-state accretion. We show an application to accreting low mass X-ray binaries (LMXBs) with accretion onto a neutron-star as compact object. The computational efficiency of $\textsc{StarShot}$ allows us to explore the parameter space for stable burning regimes, and can be used to generate initial conditions for time-dependent evolution models.