arxiv: v2 [astro-ph.sr] 23 Aug 2013

Transcription

1 Astronomy & Astrophysics manuscript no. aic c ESO 2018 July 18, 2018 Evolution towards and beyond accretion-induced collapse of massive white dwarfs and formation of millisecond pulsars T. M. Tauris 1,2,, D. Sanyal 1,3, S.-C. Yoon 1,4, and N. Langer 1 1 Argelander-Insitut für Astronomie, Universität Bonn, Auf dem Hügel 71, Bonn, Germany 2 Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, Bonn, Germany 3 Department of Physical Sciences, IISER Kolkata, Mohanpur , Nadia, West Bengal, India 4 Department of Physics & Astronomy, Seoul National University, Gwanak-gu, Gwanak-ro 1, Seoul, , Korea arxiv: v2 [astro-ph.sr] 23 Aug 2013 Received April 8, 2013; Accepted July 30, 2013 ABSTRACT Context. Millisecond pulsars (MSPs) are generally believed to be old neutron stars (NSs), formed via type Ib/c core-collapse supernovae (SNe), which have been spun up to high rotation rates via accretion from a companion star in a low-mass X-ray binary (LMXB). In an alternative formation channel, NSs are produced via the accretion-induced collapse (AIC) of a massive white dwarf (WD) in a close binary. Aims. Here we investigate binary evolution leading to AIC and examine if NSs formed in this way can subsequently be recycled to form MSPs and, if so, how they can observationally be distinguished from pulsars formed via the standard core-collapse SN channel in terms of their masses, spins, orbital periods and space velocities. Methods. Numerical calculations with a detailed stellar evolution code were used for the first time to study the combined pre- and post-aic evolution of close binaries. We investigated the mass transfer onto a massive WD (treated as a point mass) in 240 systems with three different types of non-degenerate donor stars: main-sequence stars, red giants, and helium stars. When the WD is able to accrete sufficient mass (depending on the mass-transfer rate and the duration of the accretion phase) we assumed it collapses to form a NS and we studied the dynamical effects of this implosion on the binary orbit. Subsequently, we followed the mass-transfer epoch which resumes once the donor star refills its Roche lobe and calculated the continued LMXB evolution until the end. Results. We show that recycled pulsars may form via AIC from all three types of progenitor systems investigated and find that the final properties of the resulting MSPs are, in general, remarkably similar to those of MSPs formed via the standard core-collapse SN channel. However, as a consequence of the finetuned mass-transfer rate necessary to make the WD grow in mass, the resultant MSPs created via the AIC channel preferentially form in certain orbital period intervals. In addition, their predicted small space velocities can also be used to identify them observationally. The production time of NSs formed via AIC can exceed 10 Gyr which can therefore explain the existence of relatively young NSs in globular clusters. Our calculations are also applicable to progenitor binaries of SNe Ia under certain conditions. Key words. binaries: close accretion white dwarfs supernovae stars: neutron X-rays: binaries 1. Introduction The final outcome of close binary stellar evolution is a pair of compact objects if the system avoids disruption following a supernova (SN) explosion or a merger event in a common envelope (CE). The nature of the compact objects formed can be either black holes, neutron stars (NSs), or white dwarfs (WDs), depending primarily on the initial stellar masses and orbital period of the zero-age main-sequence (ZAMS) binary. Systems in which the primary star is not massive enough to end its life as a NS may leave behind an oxygen-neon-magnesium white dwarf (ONeMg WD) of mass M WD M, which may stem from primaries with initial masses in the interval of M ZAMS 6 16 M, depending on the mass of the secondary star and, in particular, on the orbital period which affects both the important carbon/oxygen ratio at the depletion of core helium burning and the occurrence of a second dredge-up phase at the beginning of the asymptotic giant branch (see e.g. Wellstein & Langer 1999; Podsiadlowski et al. 2004). In systems where an ONeMg WD forms, the WD may later accrete sufficient material, when the secondary star sub- tauris@astro.uni-bonn.de sequently evolves and fills its Roche lobe, such that it reaches the Chandrasekhar-mass limit and implodes via an accretioninduced collapse (AIC) to form a NS (e.g. Nomoto et al. 1979a; Taam & van den Heuvel 1986; Michel 1987; Canal et al. 1990; Nomoto & Kondo 1991). The donor star in these systems can either be a main-sequence star, a low-mass giant, or a helium star. The first aim of this paper is to investigate in which binaries AIC can occur. Millisecond pulsars (MSPs) are traditionally believed to be old NSs which have been spun up to high rotation rates via accretion of mass and angular momentum from a companion star in a low-mass X-ray binary, LMXB (e.g. Bhattacharya & van den Heuvel 1991; Tauris & van den Heuvel 2006, and references therein). It is important to investigate whether or not this standard recycling scenario is the sole formation channel of MSPs. Three key questions arise in the context of AIC: 1. Could the implosion of the ONeMg WD lead directly to the formation of an MSP? 2. When an additional accretion phase is necessary for NSs formed via AIC in order to explain the observed rapid spins 1

2 and low B-fields of MSPs, is it possible at all to form MSPs via post-aic accretion from the same donor star? 3. If so, can MSPs produced this way be distinguished observationally from MSPs that formed via the standard corecollapse SN channel? In the context of direct and indirect formation of MSPs via AIC (i.e. prompt or after additional accretion, respectively), we point out that the origin of pulsar B-fields is not well understood. In direct MSP formation, B-fields could be created from conservation of magnetic flux of the collapsing core, as originally hypothesized for NSs by Woltjer (1964), or as suggested by theoretical work on thermomagnetic effects during or shortly after the NS is formed (Reisenegger 2003; Spruit 2008). In principle, a Chandrasekhar mass WD with an equatorial radius of about 3000 km and a B-field of 10 3 G could, assuming flux conservation, undergo AIC and directly produce an MSP with B 10 8 G and a spin period of a few ms, equivalent to typical values of observed recycled radio MSPs (Tauris & van den Heuvel 2006) and accreting X-ray MSPs (Patruno & Watts 2012) near the spin-up line in the PṖ diagram. Producing the spin of the MSP via AIC is not a problem if just a small fraction of the spin angular momentum of the WD is conserved during the collapse (Dessart et al. 2006). Until recently, it was thought that the distribution of WD B-fields is strongly bimodal, with a large majority of WDs being non-magnetic and a smaller fraction ( 15%, primarily in binaries) having larger fields of typically 10 7 G (Wickramasinghe & Ferrario 2000; Liebert et al. 2003). However, by using more sensitive instruments it has been demonstrated that some 15 20% of WDs could have weak B- fields of the order of 10 3 G (e.g. Jordan et al. 2007). Hence, AIC may lead to the formation of NSs with a potential large range of possible B-field values, if flux conservation is at work. One the other hand, regardless of its formation mechanism any newborn NS is extremely hot and liquid, and therefore differentially rotating, which may amplify any seed magnetic field such that the B-fields of NSs formed via AIC could be similar to those formed via an iron-core collapse. In this respect, it is important that none of the more than 40 known young NSs associated with SN remnants are observed with the characteristic properties of MSPs: low B-fields and fast spin. The B-fields of MSPs ( G) are typically lower than that of young, normal pulsars ( G) by five orders of magnitude. It is therefore questionable if an MSP would form directly from any type of collapse: iron-core collapse, electron capture SN, or AIC 1. For the rest of this work, we therefore focus our attention exclusively on the possibility of indirect formation of MSPs via AIC, i.e. following the scenario where an AIC leads to the formation of a normal NS which then subsequently accretes matter, resulting in a weaker B-field and a faster spin. The AIC route of forming MSPs has three main advantages: 1) it can explain the low space velocities of many recycled pulsars and the large fraction of NSs retained in globular clusters (Bailyn & Grindlay 1990) due to both the small amount of mass lost and the small momentum kick expected to be associated with the implosion (for details of simulations of the implosion, see e.g. Kitaura et al. 2006; Dessart et al. 2006); 2) it can explain the existence of apparently young NSs in globular clusters (Lyne et al. 1996; Boyles et al. 2011); and finally 3) it may explain the presence of a number of peculiar high B-field and slowly spinning Galactic disk NSs in close binaries with 1 Although the argument for the AIC is less certain since it may not give rise to the formation of an observable remnant. (semi)degenerate companions (e.g. Yungelson et al. 2002). In addition, direct formation of MSPs via AIC could possibly help explain the postulated birthrate problem (Kulkarni & Narayan 1988) between the small number of LMXB progenitor systems and the large observed number of MSPs. A recent, detailed population synthesis study by Hurley et al. (2010) concluded that one cannot ignore the AIC route to MSP formation and that some binary MSPs in wide orbits are best explained by an AIC scenario (see also Ivanova et al. 2008, for a specialized study on NS formation in globular clusters). There are, however, many uncertainties involved in even the best population synthesis studies and in particular in the applied physical conditions for making the ONeMg WD mass grow sufficiently. The weaknesses of the AIC formation channel are that it lacks direct observational evidence of the AIC event itself and, as already mentioned, the difficulty in predicting the spin rate and the surface B-field associated with a newborn NS formed via AIC (e.g. Kitaura et al. 2006; Dessart et al. 2006, 2007). It is not clear whether or not the AIC would be an observable transient event. According to Dessart et al. (2006), during the AIC only a few M of material is ejected (of which 25% is 56 Ni, decaying into 56 Fe via 56 CO) which quickly becomes optically thin. Hence, AIC events are most likely underluminous and very short-lived. The studies by Kitaura et al. (2006); Metzger et al. (2009); Darbha et al. (2010) also yield somewhat small amounts (< M ) of 56 Ni ejected in the AIC process, possibly synthesized in a disk, which may result in a radioactively powered SN-like transient that peaks after 1 day with a bolometric luminosity erg s 1. It is also possible that a transient radio source may appear, lasting for a few months, following the AIC event (Piro & Kulkarni 2013). In any case, these amounts are small enough to justify our assumption that the whole WD mass gets converted into the mass of the newborn NS (see Sect. 3). Another issue is that any NS formed via AIC may shortly afterwards begin to accrete additional material from its companion star, once this donor star re-fills its Roche lobe when recovering from the dynamical effects of the implosion (partly caused by the released gravitational binding energy in the transition from a WD to a more compact NS). Therefore, regardless of its initial properties, any NS formed via AIC could in principle subsequently be spun up to become an MSP, as suggested by Helfand et al. (1983). The problem is, as pointed out by Hurley et al. (2010), that this post-aic accretion phase should then resemble the conditions under which normal, old NSs are spun up to become MSPs via the conventional channel and, consequently, one cannot easily distinguish the outcome of this formation path from the standard scenario. In this work, we therefore concentrate on answering the second and the third questions raised above, i.e. if MSPs can be produced via AIC events which are immediately followed by subsequent mass transfer and, if so, how they can be distinguished observationally from those MSPs formed via the standard SN channel. We aim at investigating which progenitor binaries lead to AIC in the first place and we present detailed modelling of both pre- and post-aic evolution to predict the properties of MSPs formed via the indirect AIC channel. The structure of our paper is as follows: In Sect. 2 we review the suggested observational evidence for NS production via AIC. The computer code and our assumptions governing the pre- and post-aic mass-transfer processes are given in Sect. 3. In Sects. 4, 5, and 6 we present those of our calculated systems which successfully lead to AIC with main-sequence star, giant star, and helium star donors, respectively, and review the properties of the binary pulsars formed. 2

3 We discuss our results in a broader context in Sect. 7 and summarize our conclusions in Sect Observational evidence for AIC The question of the origin of NSs is closely related to many of their observable parameters: spin, B-field, age, space velocity, and the nature of their companion star. As already pointed out, the observational evidence suggested in the literature for NSs formed via AIC can be categorized into three groups. We now review this evidence in more detail The role of NS kicks It has been well established from observations of radio pulsar velocities that most NSs receive a momentum kick at birth (Lyne & Lorimer 1994; Hobbs et al. 2005). These kicks are possibly associated with SN explosion asymmetries and may arise from non-radial hydrodynamic instabilities (neutrinodriven convection and the standing accretion-shock instability) in the collapsing stellar core. These instabilities lead to large-scale anisotropies of the innermost SN ejecta, which interact gravitationally with the proto-ns and accelerate the nascent NS on a timescale of several seconds (e.g. Janka 2012; Wongwathanarat et al. 2013). For the entire population of NSs, the range of required kick velocity magnitudes extends basically from a few 10 km s 1 to more than 1000 km s 1, in order to explain both the existence of NSs residing inside globular clusters (which have small escape velocities, v esc < 50 km s 1 ) as well as bow shocks and SN remnants associated with highvelocity pulsars. However, when considering only young NSs in the Galactic disk the study by Hobbs et al. (2005) is interesting; it revealed that the velocities of young (< 3 Myr) radio pulsars are well described by a single Maxwellian distribution with a three-dimensional mean speed of 400 km s 1. Furthermore, there are no detections of any low-velocity (v < 60 kms 1 ) single radio pulsars with a characteristic age,τ<1myr. These facts indicate that NSs which formed recently in young stellar environments (the Galactic disk) received large kicks and that iron core-collapse SNe of type II and type Ib/c therefore, in general, result in these large kicks. On the other hand, about half of the approximately 300 known MSPs are detected in globular clusters (Ransom et al. 2005; Freire et al. 2008). Obviously, pulsars retained in globular clusters (GCs) cannot have formed with large kicks since these clusters have small escape velocities, except in a few rare cases where an isolated low-velocity pulsar could form in a disrupted binary involving a large kick with a finetuned direction (cf. Fig. 5 in Tauris & Takens 1998). It is therefore tempting to believe that many of these MSPs in GCs were not formed by iron core-collapse SNe Electron capture SNe It seems clear that the lowest mass SN progenitors may not evolve all the way to form iron cores (see Langer 2012, for a recent review on pre-sn evolution of massive single and binary stars). The final fate of these stars with ONeMg cores is an electron-capture SN (EC SN), i.e. a collapse triggered by loss of pressure support owing to the sudden capture of electrons by neon and/or magnesium nuclei (e.g. Nomoto 1984; Wheeler et al. 1998). Work by Poelarends et al. (2008) shows that the initial mass range for EC SNe is quite narrow, only about Table 1. Neutron stars that are candidates for being formed via AIC in a globular cluster (a d) or in the Galactic disk (e h), respectively. See text for explanations and discussion. Object P B P orb Mcomp Ref. ms G days M PSR B a PSR J A b PSR J B c PSR J C d GRO J e PSR J f PSR B g 4U h B-field values calculated from Eq.(5) in Tauris et al. (2012) which includes a spin-down torque due to a plasma-filled magnetosphere. Median masses calculated for i=60 and M NS = 1.35 M. a) Lyne et al. (1993); b) Lyne et al. (1996); c) Biggs et al. (1994); d) Boyles et al. (2011). e) van Paradijs et al. (1997); f) Breton et al. (2007); g) Sutantyo & Li (2000); h) Yungelson et al. (2002); 0.25 M wide, which would imply that some 4% of all singlestar SNe would be of this type. However, it has been suggested by Podsiadlowski et al. (2004) that EC SNe could occur in close binaries for stars with masses between 8 11 M since these stars lose their envelopes via mass transfer before entering the AGB phase and thus avoid the dredge-up and the consequent erosion of the CO core by this process. Therefore, these stars undergo EC SNe rather than becoming ONeMg WDs, the likely outcome of most single stars of the same mass. Furthermore, these authors argue that EC SNe lead to prompt explosions (rather than slow, delayed neutrino-driven explosions) that naturally produce NSs with low-velocity kicks (see also van den Heuvel 2004, who proposed similar ideas). The idea of different NS kick magnitudes comes from the discovery of two classes of Be/X-ray binaries with significantly different orbital eccentricities (Pfahl et al. 2002). Furthermore, the low eccentricities and the low masses ( 1.25 M ) of second-born NSs in double NS systems supports this picture (Schwab et al. 2010, and references therein) The role of young NSs in GCs In Table 1 we list a number of apparently young NSs (characterized by slow spin and relatively high B-fields) that are found in GCs. The lifetime as an observable radio source is of the order of 100 Myr for a young (i.e. non-recycled) pulsar. Therefore, if these NSs had formed via iron core-collapse SNe their existence in GCs would not only be unlikely for kinematic reasons (as explained above), it would simply be impossible given that the stellar progenitor lifetimes of SNe II and SNe Ib/c are less than a few 10 Myr, much shorter than the age of the many Gyr old stellar populations in GCs. Similarly, the nuclear evolution timescales of stars undergoing EC SNe is of the order of Myr, which is still short compared to the age of GCs, and for this reason also EC SNe cannot explain the existence of young NSs in GCs today. It is therefore clear that these NSs in GCs, if they are truly young 2, are formed via a different channel. 2 For an alternative view, see Verbunt & Freire (2013) who argue that these NSs that appear to be young, are not necessarily young. 3

4 A strong link to AIC first piece of evidence An AIC event is not very different from an EC SN and it is therefore expected that also NSs formed via AIC will receive a small kick (if any significant kick at all). For this reason formation via AIC could explain the many NSs in GCs and, more importantly, also the young ages of some of these NSs. As we shall demonstrate in this work, pre-aic binaries may, in some cases, reach ages exceeding 10 Gyr before a low-mass giant companion star initiates Roche-lobe overflow (RLO) leading to the AIC event. Therefore, we would expect ongoing AIC events, and thus formation of newborn NSs, in GCs even today NS formation via the merger of two massive WDs? For the sake of completeness, we also mention that the merger of two massive WDs may also produce a pulsar (Saio & Nomoto 1985, 2004). This could be an alternative way of producing young pulsars in an old stellar population like a GC. Even binary pulsars may be produced this way in a GC, since in a dense environment a single produced NS can capture a companion star later on (Ransom et al. 2005). It is more difficult for this scenario to produce MSPs in binaries in the Galactic disk. This scenario would not only require an initial triple system origin, which is somewhat rare in the Galactic disk (although recent work by Rappaport et al. (2013) suggests that at least 20% of all close binaries have tertiary companions), it would also require the remaining binary orbit to survive the dynamical effects of the merger event. However, it should be noted that the merging double CO WD event is also a key scenario (the socalled double-degenerate scenario) suggested as a progenitor of SNe Ia (cf. Webbink 1984; Iben & Tutukov 1984; Yoon et al. 2007; van Kerkwijk et al. 2010; Pakmor et al. 2012) AIC candidates in the Galactic disk The evidence for AIC is found not only in GCs. In Table 1 we list a number of Galactic disk binary NS systems which are postulated candidates for having formed via AIC. A common feature of these NSs is a slow spin and a relatively high B-field and an ultra-light ( 0.10 M ) companion star in a close orbit. The idea that the origin of some high B-field, slow spinning NSs (e.g. 4U , Her X 1, and PSR B ) is associated with AIC was originally suggested by Taam & van den Heuvel (1986). Although it was believed at that time that B-fields decay spontaneously on a timescale of only 50 Myr (and therefore these NSs could not have much larger ages), many of these sources remain good candidates for AIC today even though it has been demonstrated that pulsar B-fields can remain high on much longer timescales (Kulkarni 1986; Bhattacharya et al. 1992). One reason why these NS systems remain good AIC candidates is the very small masses of their companion stars which indicate that a significant amount of material ( M ) was transfered towards the compact object 3. The paradox is therefore that these NSs still have high B-fields and slow spins even though a significant mass transfer has occurred (see below). 3 A small companion-star mass suggests that the previous (or ongoing) mass-transfer episode was (is) dynamically stable (Tauris & Savonije 1999; Podsiadlowski et al. 2002). However, even observed radio pulsars which are thought to have evolved via a CE phase have B< G (Tauris et al. 2012), in general much smaller that the B-fields of the NSs listed in Table The role of accretion-induced B-field decay in NSs There is solid observational evidence that the surface B-field strengths of NSs decrease with accretion (Taam & van den Heuvel 1986; Shibazaki et al. 1989; van den Heuvel & Bitzaraki 1994, 1995). The exact mechanism for this process is still unknown. It may be related to decay of crustal fields by ohmic dissipation and diffusion from heating via nuclear processing of accreted material (Romani 1990; Geppert & Urpin 1994; Konar & Bhattacharya 1997), burial (screening) of the field (Zhang 1998; Cumming et al. 2001; Payne & Melatos 2007), or decay of core fields due to flux tube expulsion from the superfluid interior induced by rotational slow-down in the initial phases of mass accretion (Srinivasan et al. 1990); see also review by Bhattacharya (2002). Even a small amount of material accreted may lead to significant B-field decay, contradicting the observational evidence that these binary NSs have accreted large amounts of material A strong link to AIC second piece of evidence A possible solution to the above-mentioned paradox of observing close binaries with ultra-light companions orbiting high B- field NSs, would be if these NSs were formed recently via AIC during the very final stages of the mass-transfer process in a binary. The associated required finetuning of the AIC event to occur near the termination of the mass-transfer phase is important for preventing accretion of significant amounts of matter after the formation of the NS (resulting in low B-fields and fast spin). This finetuning problem will be investigated further in this paper. We note that the slow spins are expected from efficient loss of rotational energy due to the emission of magnetodipole waves from these high B-field NSs. In Sect we return to additional discussions of observational evidence for AIC in view of our theoretical calculations, and also comment on a handful of recycled radio pulsars in the Galactic disk with puzzling characteristics. 3. Numerical methods and physical assumptions of AIC The numerical models presented in this work are divided into two parts: 1) evolution prior to AIC and 2) post-aic LMXB evolution. Both parts are computed with a binary stellar evolution code originally developed by Braun (1997) on the basis of a single-star code (Langer 1998, and references therein). It is a one-dimensional implicit Lagrangian code which solves the hydrodynamic form of the stellar structure and evolution equations (Kippenhahn & Weigert 1990). The evolution of the donor star, the mass-transfer rate, and the orbital separation are computed simultaneously through an implicit coupling scheme (see also Wellstein & Langer 1999) using the Roche-approximation in the formulation of Eggleton (1983). To compute the mass-transfer rate, we use the prescription of Ritter (1988). In Sect. 5.3 we discuss the limitations of this description in wide-orbit LMXBs with giant donor stars. We employ the radiative opacities of Iglesias & Rogers (1996), which we interpolated in tables as function of density, temperature, and chemical element mass fractions, including carbon and oxygen. For the electron conduction opacity, we follow Hubbard & Lampe (1969) in the nonrelativistic case, and Canuto (1970) in the relativistic case. The stellar models are computed using extended nuclear networks including the PP I, II, and III chains and the four CNO-cycles. 4

5 In our default models we assumed a mixing-length parameter of α=l/h p = 1.5 (Langer 1991) and a core convective overshooting parameter ofδ ov = We tested several models using α=l/h p = 2.0 which resulted in only slightly larger final WD masses ( 1%) orbiting recycled pulsars in somewhat larger orbits (up to 3% increase in P orb ). If the accreting ONeMg WD reached the limiting Chandrasekhar mass for a rigidly rotating WD (i.e. M Chan = 1.48 M, e.g. Yoon & Langer 2005) we assumed it collapsed to form a NS. Differential rotation can persist if the timescale of angular momentum transport is smaller than the accretion timescale in accreting WDs, and leads to a critical mass that is significantly higher than the canonical Chandrasekhar mass (e.g. Yoon & Langer 2004). However, magnetic torques resulting from the Spruit-Tayler dynamo may enforce nearly rigid rotation in accreting WDs with the considered accretion rates in this study. Further research is needed to verify this. The WD collapse was modelled both with and without a momentum kick imparted to the newborn NS. Although AIC is generally believed to result in no kick, or possibly a small kick, we applied three different kick values (w=0, 50, 450 km s 1 ) to the newborn NS. Given that the physics behind the kick mechanism is still unclear, we have included a few extra models with a large kick of 450 km s 1 to probe the extreme boundary conditions of our calculations. We solved for the combined effects of sudden mass loss and imparted momentum kick on the orbital dynamics following Hills (1983). In all cases we assumed the mass equivalent of 0.20 M was lost as released gravitational binding energy yielding a post-aic NS gravitational mass of 1.28 M. The post-aic LMXB evolution was followed using the same computer code (see also Tauris et al. 2011, 2012, for additional details). For a more general discussion of MSP formation from LMXBs we refer to e.g. Ergma et al. (1998); Tauris & Savonije (1999); Podsiadlowski et al. (2002); Deloye (2008) Accretion onto a white dwarf The mass-transfer process and the physics of accretion onto a WD has been described in a large number of papers, e.g. Whelan & Iben (1973); Nomoto & Sugimoto (1977); Nomoto (1982); Prialnik & Kovetz (1995); Iben & Tutukov (1996); Li & van den Heuvel (1997); Hachisu et al. (1999); Langer et al. (2000); Livio (2000); Han & Podsiadlowski (2004); Yoon & Langer (2003, 2004); Nomoto et al. (2007), and more recently in Hachisu et al. (2012); Wheeler (2012); Idan et al. (2012); Starrfield et al. (2012); Newsham et al. (2013); Denissenkov et al. (2013); Ma et al. (2013). Most of these papers aim at investigating progenitors of type Ia supernovae (SN Ia) which are important for cosmology studies of the accelerating Universe (e.g. Riess et al. 1998; Perlmutter et al. 1999). In these progenitor systems an accreting CO WD reaches the Chandrasekhar limit and explodes. Observationally, these systems manifest themselves as cataclysmic variables (Hellier 2001), symbiotic systems (Kenyon 1986), and supersoft X-ray sources (van den Heuvel et al. 1992; Kahabka & van den Heuvel 1997), depending on the companion star mass, its evolutionary status, and the mass-transfer rate. It should be noted that the classification scheme overlaps somewhat on certain aspects. The cataclysmic variables (CVs) can be further subdivided into several classes: (classical) novae, recurrent novae, dwarf novae, polars, AM CNn stars etc., again depending on the nature of the companion star, the orbital separation, as well as the magnetic field strength of the accreting WD and accretion disk morphology. The classical novae have only been reported to erupt once whereas the recurrent novae usually erupt every one or two decades. The symbiotic variable star binaries have low-mass red giant donors which transfer mass at a fairly low rate via (beginning atmospheric) RLO. They undergo nova-like outbursts which last for a few decades before decaying back to their original luminosity. The luminous persistent supersoft X- ray sources, on the other hand, display luminosities between erg s 1 revealing high mass-transfer rates of the order of 10 7 M yr 1. These systems often have a more massive mainsequence donor star undergoing thermal timescale mass transfer. It is generally believed that the accreting WDs in supersoft X-ray sources can grow in mass by accretion since the thermonuclear fusion of accreted hydrogen can be fairly stable with these high accretion rates, in contrast to the case of nova systems where nova explosions (caused by a violent ignition in a thin shell of degenerate hydrogen) may erode the WDs (Wiescher et al. 1986; Patterson et al. 2013). The supersoft X-ray sources are therefore believed to represent progenitors of SNe Ia, because many of these systems may produce a Chandrasekhar mass CO WD. However, existing multicycle computations of hydrogen accretion onto massive WDs at a high rate are still somewhat controversial. In a recent study Idan et al. (2012) showed that the accumulated helium is completely lost in strong helium flashes, thereby making SN Ia and AIC impossible, whereas another study by Newsham et al. (2013) concluded that WDs continue to grow toward the Chandrasekhar limit. A number of numerical simulations show that when a CO WD approaches the Chandrasekhar limit with accretion rates of the order of Ṁ WD M yr 1, a thermonuclear-runaway caused by carbon burning can occur at central densities of aboutρ c = g cm 3 (e.g. Nomoto et al. 1984; Yoon & Langer 2003; Lesaffre et al. 2006), which may result in a SN Ia explosion. In case an ONeMg WD mass grows to the Chandrasekhar limit, electron-captures onto 24 Mg and 20 Ne make the central density increase to about g cm 3 before oxygen ignites at the center (Miyaji et al. 1980; Miyaji & Nomoto 1987; Nomoto 1987). The consequent oxygen deflagration with this high density cannot lead to a thermonuclear explosion because of very rapid electroncaptures onto heavy elements produced by the oxygen burning (Miyaji et al. 1980; Timmes & Woosley 1992). Therefore, a NS is the most likely outcome in the case of the collapse of an ONeMg WD. Although it is generally believed that accreting CO WDs reaching the Chandrasekhar mass limit lead to a SNe Ia, and that accreting ONeMg WDs reaching this limit undergo AIC and produce NSs, the outcome could in some cases be the opposite see Sect. 7.2 for further discussion Dependence on mass-transfer rates The response of the WD to mass transfer from the donor star depends on the mass-transfer rate, Ṁ 2. In this work we assume that the WD receives (but does not necessarily accrete) mass at the same rate as it is lost from the donor star, Ṁ 2. White dwarfs which accrete faster than the rate at which their cores can grow (Ṁ up ) will puff up their envelope to giant star dimensions (see below). If the WD envelope does expand to such huge radii it will engulf the donor star and the system is likely to evolve through a common envelope (CE) and spiral-in phase (Paczyński 5

6 Fig. 1. The accretion window of the WD, the corresponding critical mass-transfer rates, and our assumed mass-accumulation fractions,η H, here calculated with a hydrogen abundance of X= 0.70 of the accreted matter (see text). 1976; Iben & Livio 1993). The relevant critical mass-transfer rates for our work are (see Fig. 1): Ṁ CE Ṁ up Ṁ Edd,WD Ṁ steady Ṁ accu upper limit before the formation of a giant and a CE upper limit for steady hydrogen burning Eddington limit for spherical accretion lower limit for steady hydrogen burning lower limit for WD mass accumulation where we have adopted the following ad hoc condition: Ṁ CE = 3 Ṁ Edd,WD. (1) The Eddington accretion limit is found by equating the outward radiation pressure to the gravitational force per unit area acting on the nucleons of the accreted plasma. The radiation pressure (from photons that scatter on plasma electrons) is generated from both nuclear burning at the WD surface and from the release of gravitational binding energy of the accreted material. The total energy production is given by: L = (ǫ nuc +ǫ acc ) Ṁ WD, where Ṁ WD is the accretion rate of the WD andǫ nuc andǫ acc denote the specific energy production from nuclear burning at the WD surface and release of gravitational binding energy, respectively. Hence, the Eddington accretion limit depends on the chemical composition of the accreted material and on the mass of the WD and is roughly given by: Ṁ Edd,WD ( ) 10 7 M yr 1 (2) for the WDs studied in this work. The value of Ṁ CE is probably the largest uncertainty in our modelling. The rationale of even considering an accretion rate Ṁ 2 > Ṁ Edd,WD is that the WD may drive a strong wind in a bipolar outflow which may prevent the WD envelope from otherwise expanding into giant star dimensions (Nomoto et al. 1979b). Many previous studies of accreting WDs have adapted the optically thick wind model of Kato & Iben (1992) and Kato & Hachisu (1994), without any restrictions on Ṁ 2 (see Cassisi et al. 1998; Langer et al. 2000, for a critique of this assumption). Given the uncertainties regarding the validity of this model, we adopt a maximum allowed mass-transfer rate limit of Ṁ 2 = Ṁ CE = 3 Ṁ Edd,WD. For most of our models we stopped the calculations if Ṁ 2 > Ṁ CE (assuming the system evolved into a CE and merged). To test the dependence of this limit we also computed some models by allowing Ṁ CE = 10 Ṁ Edd,WD for comparison. As we shall see, the role of adapting the optically thick wind model or not has important consequences for the progenitor parameter space leading to AIC. Constraints of the wind mass loss from an accreting WD can be determined directly from radio observations of SN Ia remnants (Chomiuk et al. 2012). Hence, there is some hope that future observations can clarify the situation. If a wind is driven from the WD (when Ṁ 2 > Ṁ up ) we calculate its mass-accumulation fraction using:η H = Ṁ up / Ṁ 2 to restrict the accretion rate to a maximum of Ṁ up. This value represents the upper limit for steady shell hydrogen burning of a WD and can be estimated from the growth rate of the degenerate core in red giant stars undergoing hydrogen shell burning by applying the relation between core mass and luminosity (Iben & Tutukov 1989). As an example, Hachisu & Kato (2001) found: Ṁ up = M yr 1 ( MWD M 0.40 )( 1.7 X X ), (3) which is valid for a hydrogen mass fraction X There is a variety of similar expressions for this critical limit in the literature (e.g. Nomoto 1982; Hachisu et al. 1996, 1999; Nomoto et al. 2007). However, they do not differ by much and our results are stable against these minor variations. (We note that Ṁ up Ṁ Edd,WD ). The minimum value for steady hydrogen shell burning is given by Nomoto (1982): Ṁ steady = 0.4 Ṁ up. (4) Although the hydrogen burning is not steady below this limit, the shell flashes are found to be weak for accretion rates just slightly below the limit. Following Hachisu et al. (1999) we therefore assumeη H = 1 in the entire interval: Ṁ accu < Ṁ 2 < Ṁ up, where Ṁ accu = 1/8 Ṁ up. (5) If the mass-accretion rate is below Ṁ accu violent shell flashes and nova outbursts cannot be avoided and thus the WD is prevented from increasing its mass (i.e.η H = 0), or may even erode. Following hydrogen burning the helium is processed into carbon and oxygen. The mass accumulation efficiency in helium shell flashes was studied in detail by Kato & Hachisu (2004). We have adapted their mass accumulation efficiencies for helium burning,η He into our code. The long-term effective massaccretion rate of the WD is therefore given by: Ṁ WD =η H η He Ṁ 2. (6) For accretion of pure helium we used Ṁ WD =η He Ṁ 2 and ( ) Ṁ up,he = M yr 1 MCO 0.6, (7) M where M CO is the mass of the CO core of the helium donor star (Nomoto 1982). According to Jose et al. (1993),η He might be somewhat smaller for direct accretion of helium, which leads to stronger shell flashes, compared to the case where helium is accumulated via multiple cycles of hydrogen burning (i.e. double shell burning). For recent discussions on the WD growth rate and the dependence on the WD mass and the mixing of the accreted material, see e.g. Denissenkov et al. (2013) and Newsham et al. (2013). 6

7 3.2. Orbital dynamics We consider close interacting binary systems which consist of a non-degenerate (evolved) donor star and a compact object, in our case initially a massive WD and later on, in the case of an AIC event, a NS. When the donor star fills its Roche lobe, any exchange and loss of mass from such an X-ray binary will also lead to alterations of the orbital dynamics, via modifications in the orbital angular momentum, and hence changes in the size of the critical Roche-lobe radius of the donor star. The stability of the mass-transfer process therefore depends on how these two radii evolve (i.e. the radius of the star and its Roche-lobe radius). The various possible modes of mass exchange and mass loss include, for example, fast wind mass loss (Jeans mode), Roche-lobe overflow (with or without isotropic re-emission), and common envelope evolution (e.g. van den Heuvel 1994; Soberman et al. 1997, and references therein). The RLO mass transfer can be initiated while the donor star is still on the main sequence (Case A RLO), during hydrogen shell burning (Case B RLO), or during helium shell burning (Case C RLO). The corresponding evolutionary timescales for these different cases will in general proceed on a nuclear, thermal, or dynamical timescale, respectively, or a combination thereof. This timescale is important for the amount of mass that can be accreted and for the extent to which the NS produced in the AIC can be recycled after its formation. The dynamical evolution of a binary system can be found by solving for the changes in the orbital separation, a. The orbital angular momentum of a circular binary system is given by: J orb =µωa 2, whereµis the reduced mass and the orbital angular velocity is:ω= GM/a 3. A simple logarithmic differentiation of the orbital angular momentum equation yields the rate of change in orbital separation: ȧ a = 2 J orb 2 Ṁ1 2 Ṁ2 + Ṁ1+Ṁ 2 J orb M 1 M 2 M, (8) where the two stellar masses are given by M 1 and M 2, the total mass is M=M 1 +M 2, and the total change in orbital angular momentum per unit time is given by: J orb = J gwr + J mb + J ls + J ml. These four terms represent gravitational wave radiation, magnetic braking, other spin-orbit couplings, and mass loss, respectively (e.g. Tauris & van den Heuvel 2006, and references therein). In this work we adopt the so-called isotropic re-emission mode 4 for modelling the mass-transfer and the mass loss from the binary (e.g. Bhattacharya & van den Heuvel 1991; Soberman et al. 1997). We let the mass ratio between the donor star (M 2 ) and the accretor (M 1 ) be denoted by q = M 2 /M 1. Assuming that the direct wind mass loss of the donor star is negligible compared to the RLO mass-transfer rate from the donor star, Ṁ 2, and ignoring mass stored in a circumbinary torus, one can show that a binary system always widens (ȧ>0) as a result of mass transfer if q < 1. Similarly, a binary always decreases (ȧ<0) if q>(1+ 17)/4 1.28, irrespective of the amount of mass ejected from the vicinity of the accretor (e.g. see Sect in Tauris & van den Heuvel 2006). This fact is worth remembering when we analyse our results. 4 Although at present there is not much observational evidence behind this model (not surprising given that systems with extremely high mass-transfer rates are short lived) there is some evidence of excess mass loss in the case of Cyg X-2, as demonstrated by King & Ritter (1999) and Podsiadlowski & Rappaport (2000). This mass loss could have been be caused in the past by a relativistic jet (as observed in SS433) or by coronal winds from the outer parts of the accretion disk (Blandford & Begelman 1999). In our calculations we have ignored any changes in J orb due to magnetic braking or any other tidal spin-orbit interactions. Magnetic braking is only known to operate efficiently in lowmass stars (< 1.5 M ) which have convective envelopes. As we shall see, our calculations show that only main-sequence donor stars with masses in the range M (depending on metallicity) transfer mass at the rates needed for AIC events to take place. Therefore, to be consistent with all our calculations, we have not included magnetic braking in the very few borderline cases. (These low-mass donors all have low-metallicities resulting in some stability against convective envelopes.) Our code includes gravitational wave radiation by calculating J gwr and correcting for it. However, for these AIC progenitor systems the RLO mass-transfer timescales are always significantly shorter than the timescales on which gravitational wave radiation is important. In all our calculations we assumed an initial (ONeMg) WD mass of M WD = 1.20 M prior to accretion. As mentioned earlier, for these calculations the WD is simply treated as a point mass The dynamical effect of the AIC on the orbit If the WD mass reached 1.48 M we assumed that the WD was subject to instantaneous AIC. The effect of sudden mass loss in a binary has been studied in detail by Hills (1983) for bound systems, and by Tauris & Takens (1998) for disrupted systems. Assuming a circular pre-aic orbit we here follow Hills (1983) to find the changes of the binary orbital parameters. The change of the binary semi-major axis, as a result of an asymmetric AIC, is given by: a a 0 = [ ] 1 ( M/M 0 ), (9) 1 2( M/M 0 ) (w/v c ) 2 2 cosθ (w/v c ) where a 0 is the pre-aic semi-major axis (radius), a the post-aic semi-major axis, M= 0.20 M the effective mass loss during the AIC when the 1.48 M WD is compressed to a NS with a gravitational mass of 1.28 M (Zeldovich & Novikov 1971), M 0 = M WD +M 2 the pre-aic total mass, v c = GM 0 /a 0 the pre- AIC orbital velocity of the collapsing WD in a reference fixed on the companions star, w the magnitude of the kick velocity, and θ the angle between the kick velocity vector, w and the pre-aic orbital velocity vector, v c. The post-aic eccentricity is given by: e= 1+ 2E orbjorb 2, (10) µ G 2 MNS 2 M2 2 where the post-aic orbital energy of the system is given by: E orb = GM NS M 2 /2a, and the orbital angular momentum is given by: J orb = a 0 µ (v c + w cosθ) 2 + (w sinθ sinφ) 2, (11) whereµis the post-aic reduced mass andφis the angle between the projection of the kick velocity vector onto a plane perpendicular to the pre-aic velocity vector of the WD and the pre-aic orbital plane. We neglected any shell impact effects on the companion star since the amount of material ejected in an AIC event is expected to be negligible. Even for SN Ib/c where a significant shell is ejected, the impact effect on the orbital dynamics is small if the pre-sn separation is larger than a few R (Tauris & Takens 1998). 7

8 Following the AIC we checked if the post-aic periastron separation, a(1 e) is smaller than the radius of the companion star, R 2. In that case we assumed that the system merges. The post-aic orbit is expected to be circularized with time and we therefore assumed that the tidal interactions reduced the semimajor axis by a factor of (1 e 2 ) in order to conserve J orb. When no momentum kick was added (w=0) to the newborn NS, the relation between post-aic orbital separation (including the subsequent effect of tidal circularization), a circ and the pre- AIC orbital radius, a 0 is simply given by (Verbunt et al. 1990; Bhattacharya & van den Heuvel 1991): M 0 a circ = a 0 M = a 1.48 M + M 2 0. (12) 1.28 M + M 2 Since post-aic evolution calculations also require a significant amount of computational time, and each AIC event can lead to a large set of possible parameter outcomes for (w,θ,φ), we have restricted ourselves to those cases which best probe the extreme cases for the post-aic evolution with respect to orbital periods and systemic recoil velocities resulting from the AIC event Post-AIC LMXB evolution The evolution of post-aic binaries is, in principle, similar to normal LMXB evolution (i.e. a donor star which transfers matter and angular momentum to an accreting NS). For this recycling process we follow Tauris & Savonije (1999) and Tauris et al. (2012) for tracking the evolution of the LMXB (see also Sect. 3.2 below). The accretion rate onto the NS is assumed to be Eddington limited and is given by: Ṁ NS = ( Ṁ 2 max [ Ṁ 2 Ṁ Edd, 0 ]) e acc k def, (13) where e acc is the fraction of matter transfered to the NS which actually ends up being accreted and remains on the NS, and k def is a factor that expresses the ratio of gravitational mass to rest mass of the accreted matter (depending on the equation-of-state of supranuclear matter k def ; e.g. Lattimer & Prakash 2007). Here we assumed e acc k def = Our motivation for this value is the increasing evidence of inefficient accretion in LMXBs, even in close systems where the mass-transfer rate is expected to be sub-eddington ( Ṁ 2 < Ṁ Edd ) at all times (e.g. Jacoby et al. 2005; Antoniadis et al. 2012). Possible mechanisms for inefficient accretion include propeller effects, accretion disc instabilities, and direct irradiation of the donor s atmosphere from the pulsar (Illarionov & Sunyaev 1975; van Paradijs 1996; Dubus et al. 1999). For the post-aic NS, we calculated the Eddington mass-accretion rate using: Ṁ Edd = M yr 1 M 1/3 2 NS 1+X. (14) In Sect we specify a relation between the amount of mass accreted and the final pulsar spin period Complete evolution in the (P orb, M 2 )-plane To demonstrate the orbital evolution: 1) during pre-aic mass transfer, followed by 2) post-aic mass transfer in an LMXB, and to show the effect of possible kicks associated with the AIC, we have plotted complete evolutionary tracks in the (P orb, M 2 )- plane in Fig. 2. The two donor stars in these examples (M 2 = 2.6 M and M 2 = 2.2 M ) are both more massive than the accreting ONeMg WD (initially with mass ratio, q 2.) This explains why both systems decrease in P orb prior to the AIC event. At the moment of the AIC the orbits widen instantaneously as a consequence of the sudden mass loss (c.f. Sect ). The larger the kick, w, the larger the post-aic P orb becomes at which the donor star refills its Roche lobe and continues mass transfer to the newborn NS in the LMXB source. The final products are binary MSPs with He WDs. During the LMXB phase the orbit changes from a converging system to a diverging system when the mass ratio inverses (the exact value depends on the masstransfer rate, c.f. Sect. 3.2). It is interesting that for the original 2.2 M donor star, the post-aic mass transfer is not dynamically stable if the AIC was asymmetric (i.e. if w = 50 km s 1 or w = 450 km s 1 ). In this case the post-aic binary becomes wide enough that the donor star develops a deep convective envelope before refilling its Roche lobe. The result is that the subsequent mass-transfer stage leads to excessive mass-transfer rates and thus to the formation of a CE. We did not follow the evolution of these systems further, although it is possible, in principle, that some of the donor star envelopes would be loosely enough bound to allow for ejection during spiral-in and thereby leave behind a mildly recycled pulsar orbiting a WD in a tight orbit. We return to this possibility in Sect AIC in systems with main-sequence star donors In the following three sections we present our results of AIC calculations in systems with main-sequence, giant, and helium star donors, respectively. We evolved a total of 240 binary systems. Key parameters from 36 examples of complete calculations (i.e. both pre-aic and post-aic LMXB computations leading to a recycled binary pulsar) are given in Table 2. Model names beginning with MS and MZ refer to main-sequence donor stars at solar metallicity (Z = 0.02) and Z = 0.001, respectively; giant donor star models are denoted by GS (or GSZ for Z=0.001), and helium donor star models are denoted by He. The given parameters for each model are the following: M ZAMS 2 and P ZAMS orb refer to initial donor star mass and orbital period; t RLO is the age of the donor star when it initiates RLO; and X c is its central hydrogen content at that time; t CV is the duration of the masstransfer phase until the AIC event; and M2 AIC and P AIC orb are the donor star mass and the orbital period at the moment of the AIC. When a momentum kick (w>0) is added to the newborn NS, its magnitude and direction are given. The resulting systemic recoil velocity of each post-aic system is given by v sys. As a result of the AIC the system temporarily detaches. The time it takes until the donor star refills its Roche lobe is denoted by t detach, and P circ orb is the orbital period at that time (after circularization). The duration of the subsequent post-aic LMXB phase is given by t LMXB. Finally, the parameters M WD, P MSP orb, M NS, M NS, P spin, and t total denote the WD mass, the orbital period, the NS mass, the amount of mass accreted by the NS during the post-aic mass transfer, the final spin period of the recycled pulsar, and the total age of the binary system at this end point. For hydrogen-rich donors we have initiated our calculations assuming a ZAMS star orbiting an ONeMg WD. The error in placing the companion star on the ZAMS, and neglecting the evolution of this star while the ONeMg WD forms, is not very significant. As we shall see in a moment, our main-sequence companions have maximum masses of M 2 = 2.6 M, and even these stars evolve on a much longer timescale (at least by a factor of 10) compared to the typical 6 8 M progenitor stars of ONeMg WDs. 8