,Draft version September 29, 2018 Typeset using LAT X twocolumn style in AASTeX61 , E THE INFLUENCE OF STELLAR SPIN ON IGNITION OF THERMONUCLEAR RUNAWAYS Duncan K. Galloway,1, 2 Jean J. M. in 't Zand,3 Jer´ ome^ Chenevez,4 Laurens Keek,5, 6 Celia Sanchez-Fernandez,7 Hauke Worpel,8 Nathanael Lampe,9 Erik Kuulkers,10 Anna Watts,11 and Laura Ootes11 The MINBAR collaboration 1School of Physics & Astronomy, Monash University, Clayton, VIC 3800, Australia 2Monash Centre for Astrophysics, Monash University, Clayton, VIC 3800, Australia 3SRON Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, the Netherlands 4National Space Institute, Technical University of Denmark, Elektrovej 327-328, DK-2800 Lyngby, Denmark 5X-ray Astrophysics Laboratory, Astrophysics Science Division, NASA/GSFC, Greenbelt, MD 20771 6CRESST and the Department of Astronomy, University of Maryland, College Park, MD 20742 7European Space Astronomy Centre (ESA/ESAC), Science Operations Department, E-28691, Villanueva de la Ca~nada,Madrid, Spain 8Leibniz-Institut f¨urAstrophysik Potsdam (AIP), An der Sternwarte 16, Potsdam 14482, Germany 9CNRS, IN2P3, CENBG, UMR 5797, F-33170 Gradignan, France 10ESA/ESTEC, Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands 11Anton Pannekoek Institute for Astronomy, University of Amsterdam, Postbus 94249, 1090GE Amsterdam, the Netherlands (Accepted April 9, 2018) Submitted to ApJL ABSTRACT Runaway thermonuclear burning of a layer of accumulated fuel on the surface of a compact star provides a brief but intense display of stellar nuclear processes. For neutron stars accreting from a binary companion, these events manifest as thermonuclear (type-I) X-ray bursts, and recur on typical timescales of hours to days. We measured the burst rate as a function of accretion rate, from seven neutron stars with known spin rates, using a burst sample accumulated over several decades. At the highest accretion rates, the burst rate is lower for faster spinning stars. The observations imply that fast (> 400 Hz) rotation encourages stabilization of nuclear burning, suggesting a dynamical dependence of nuclear ignition on the spin rate. This dependence is unexpected, because faster rotation entails less shear between the surrounding accretion disk and the star. Large-scale circulation in the fuel layer, leading to enhanced mixing of the burst ashes into the fuel layer, may explain this behavior; further numerical simulations are required to confirm. Keywords: stars: neutron | X-rays: bursts | nuclear reactions | stars: rotation arXiv:1804.03380v1 [astro-ph.HE] 10 Apr 2018 Corresponding author: Duncan Galloway [email protected] 2 Galloway et al. 1. INTRODUCTION at high accretion rates; Done et al. 2007) then the shear Thermonuclear runaways arise from unstable ignition between the surface and the disk material (orbiting at of accreted fuel on the surface of compact objects in bi- roughly Keplerian speeds) may contribute heating even nary systems. One such class of events has been known into the fuel layer (Inogamov & Sunyaev 2010). Because since antiquity as novae (Jos´e& Hernanz 2007), and the Keplerian motion of the disk is always faster than the occur on the surface of white dwarfs accreting from a rotation speed of the neutron star, the shear heating can stellar companion. A second class of runaway was dis- be expected to have a greater effect for slower neutron covered in the early 1970s with the advent of satellite- star spin; in contrast, rotation-induced mixing will have based X-ray astronomy (Grindlay et al. 1976; Clark et al. a greater effect at faster spin rates. 1976), and occurs on neutron stars, where the surface Observationally, variations in the persistent spectral temperature is sufficiently hot that the star emits pri- shape as a function of spin rate have been attributed to marily in X-rays. These X-ray bursts occur typically the effects of shear heating (Burke et al. 2018). The spin 106 times more frequently than novae, with recurrence rate also seems to determine in part the type of bursts times of hours rather than years. At typical accretion in which \burst oscillations" (Watts 2012) are observed. −8 −1 For rapidly spinning sources (νspin 500 Hz) burst oscil- rates (∼ 0:1 M_ Edd, where M_ Edd = 1:3×10 M yr is & lations are predominantly observed in bursts exhibiting the Eddington limit for a 1:4 M object accreting solar composition fuel), bursts occur quasi-periodically every so-called photospheric radius-expansion (PRE; Muno few hours, with durations of 10{100 s (e.g. Lewin et al. et al. 2004). The X-ray intensity during such events 1993; Galloway & Keek 2017). reaches the local Eddington limit, where the force due Few recurrent novae are known (e.g. Darnley et al. to radiation pressure from the burst emission balances 2016), although the population has expanded dramat- the (considerable) force due to gravity at the neutron ically in the last decade, thanks primarily to surveys star surface. In contrast, burst oscillations in slower- of other galaxies. The relative frequency of runaways rotating sources are prevalent equally likely in bursts on neutron stars enables much more detailed studies with or without PRE. This effect may arise partially than for novae. Most of the ≈ 100 known X-ray burst from the predominance in the slower-spinning sources sources1 within our Galaxy have shown multiple bursts, for (almost) pure He accretion (Galloway et al. 2008, with over a thousand events detected from the most pro- 2010). lific sources. In addition, many of the bursting neutron Here we present a study of the variation in the rate stars have been shown to be spinning with spin frequen- of thermonuclear runaways in response to the accretion cies in the range 11{620 Hz (Watts 2012), introducing rate in different sources, to determine the influence of the possibility of studying rotational effects that may the spin rate on thermonuclear ignition. influence the spreading of the thermonuclear flame. The influence of stellar spin on thermonuclear igni- 2. OBSERVATIONS AND ANALYSIS tion has been investigated theoretically via its effect We employ the largest available sample of thermonu- on mixing (Fujimoto 1993). At the depth of ther- clear bursts, incorporating more than 7000 events from monuclear burning on neutron stars (column depth y ∼ 85 bursting sources, observed by three long-duration 108 g cm−2) viscosity likely results in nearly uniform ro- satellite X-ray missions: the Dutch-Italian BeppoSAX, tation. Even so, turbulent mixing can occur, partic- NASA's Rossi X-ray Timing Explorer (RXTE), and ularly at high accretion rates (Piro & Bildsten 2007). ESA's INTEGRAL satellites. The sample, comprising Numerical simulations suggest that this mixing may sta- the Multi-INstrument Burst ARchive (MINBAR2), in- bilize the burning at accretion rates lower than would be cludes uniform analyses of the burst properties, as well expected in nonrotating layers (Keek et al. 2009). How- as the persistent (non-bursting) X-ray emission, from ever, most calculations have focussed on pure helium which the rate of accretion onto the neutron star can be layers, while many neutron stars accrete both hydrogen estimated. and helium. We carried out our analysis for seven sources for which Another factor that may influence the conditions in there are ≈ 100 or more bursts in MINBAR, and for the fuel layer is the shear with respect to the accretion which the spin frequency is known (Table1). The spin flow. If the accretion disk meets the neutron star surface frequency for these sources has been inferred from the forming a boundary layer (as is thought to be the case detection of burst oscillations (e.g. Watts 2012), which 1 http://burst.sci.monash.edu/sources 2 http://burst.sci.monash.edu/minbar The influence of stellar spin on ignition of thermonuclear runaways 3 3 are known to trace the neutron star spin (Chakrabarty a mean peak flux hFEddi for each source . These val- et al. 2003). ues are listed in Table1; note that the quoted uncer- We measured the (average) burst rate as a function of tainty includes the intrinsic variation of peak intensities accretion rate, as follows. We estimated the accretion of radius-expansion bursts from each source (cf. with rate for each observation from the intensity of the non- Galloway et al. 2008). bursting (persistent) emission, averaged over each obser- Rather than calculate the distance, accretion luminos- vation, largely following Galloway et al.(2008). We then ity and hence M_ , we can estimate the accretion rate as a grouped the bursts detected within each observation to fraction of the Eddington value M_ Edd more directly (fol- measure the burst rate at different accretion rates. lowing van Paradijs et al. 1988), by calculating the ratio The instantaneous accretion rate onto bursting neu- of bolometric persistent flux to the inferred Eddington tron stars can be inferred from measurements of the flux measured from the photospheric radius-expansion persistent (non-bursting) X-ray emission, and varies on bursts, i.e. timescales of minutes to decades. Some sources are tran- _ _ ξpFX cbol sient, meaning that they undergo intermittent periods M=MEdd = (1 + X) (3) ξb hFEddi of active accretion lasting weeks to months, punctuated where ξ is the equivalent anisotropy term specific to the by longer quiescent intervals lasting months to years. b burst emission (which is expected to be different from The cumulative effect of these variations for the current ξ ). This approach is equivalent to calculating the dis- study is that our observational data, which spans 16 p tance and luminosity, but omits the intermediate steps. years, samples most bursting sources over a wide range Here we include one additional correction factor (con- of intensities (or accretion rates). stant for all the sources) related to the hydrogen frac- The MINBAR observations were performed with one tion X of the material in the atmosphere.