Issue 
A&A
Volume 537, January 2012



Article Number  A104  
Number of page(s)  18  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201117880  
Published online  16 January 2012 
The evolution of ultracompact Xray binaries
^{1} Department of Astrophysics/ IMAPP, Radboud University Nijmegen, PO Box 9010, 6500 GL Nijmegen, The Netherlands
email: L.vanHaaften@astro.ru.nl
^{2} Institute for Astronomy, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
^{3} Department of Physics and Space Sciences, Florida Institute of Technology, 150 W. University Blvd., Melbourne, FL 32901, USA
Received: 12 August 2011
Accepted: 19 November 2011
Context. Ultracompact Xray binaries (UCXBs) typically consist of a white dwarf donor and a neutron star or black hole accretor. The evolution of UCXBs and very low mass ratio binaries in general is poorly understood. In particular, the dynamical behavior of an accretion disk extending to a large radius (relative to the orbit) is unclear.
Aims. We investigate the evolution of UCXBs in order to learn for which mass ratios and accretor types these systems can exist, and if they do, what are their orbital and neutron star spin periods, mass transfer rates and evolutionary timescales.
Methods. We compute tracks of a binary containing a Rochelobe overflowing helium white dwarf in which mass transfer is driven by gravitational wave emission. For different assumptions concerning accretion disk behavior we calculate for which system parameters dynamical instability, thermalviscous disk instability or the propeller effect emerge. The significance of these processes during the evolution of an UCXB is considered.
Results. At the onset of mass transfer, the survival of the UCXB is determined by how efficiently the accretor can eject matter in the case of a superEddington mass transfer rate. At later times, the evolution of systems strongly depends on the binary’s capacity to return angular momentum from the disk to the orbit. We find that this feedback mechanism most likely remains effective even at very low mass ratio. In the case of steady mass transfer, the propeller effect can stop accretion onto recycled neutron stars completely at a sufficiently low mass transfer rate, based on energy considerations. However, mass transfer will likely be nonsteady because disk instability allows for accretion of some of the transferred matter. Together, the propeller effect and disk instability cause the low mass ratio UCXBs to be visible a small fraction of the time at most, thereby explaining the lack of observations of such systems.
Conclusions. Most likely UCXBs avoid latetime dynamically unstable mass loss from the donor and continue to evolve as the age of the Universe allows. This implies the existence of a large population of low mass ratio binaries with orbital periods ~70–80 min, unless some other mechanism has destroyed these binaries. Even though none have been discovered yet, black hole UCXBs could also exist, at orbital periods of typically 100–110 min.
Key words: binaries: close / accretion, accretion disks / Xrays: binaries
© ESO, 2012
1. Introduction
Ultracompact Xray binaries (UCXBs) are a subclass of lowmass Xray binaries and consist of a white dwarf or helium star losing mass to a neutron star or black hole, at a subhour orbital period (Savonije et al. 1986). The short orbital periods point to white dwarf or helium burning star (sdB star) donors, as those are the only known stellar types that have the same size as the donor Roche lobe corresponding to these orbital periods (Nelson et al. 1986). The observed helium and carbonoxygen composition of the transferred matter confirms this picture (Schulz et al. 2001; Nelemans et al. 2004, 2006). The accretor type can be identified via its inferred mass, magnetic field, spin period, explosive nuclear fusion on its surface or inner accretion disk behavior, and can be either a neutron star or a black hole (the latter have not yet been observed in UCXBs). To date about 30 UCXBs and candidates have been identified (in’t Zand et al. 2007).
UCXBs are important objects to study because of the absence of hydrogen in the accretion disk and Xray bursts (in’t Zand et al. 2005), their gravitational wave signal (Nelemans 2009), and them being tests for the commonenvelope phase (Nelemans & Jonker 2010). Furthermore, UCXBs are candidate progenitors of radio millisecond pulsars (Alpar et al. 1982) because of the high number of accreting millisecond pulsars found in UCXBs (Wijnands 2010). The ratio between the number of black hole and neutron star accretors constrains the highmass end of the binary initial mass function.
1.1. Evolutionary history
The formation of ultracompact Xray binaries starting from a zeroage main sequence binary involves one or two common envelope stages and a supernova. The initially more massive star evolves off the main sequence and starts losing mass to its companion. Depending on the mass ratio and stage of donor evolution, this may happen in an unstable way, leading to a common envelope. The orbit decays significantly due to energy and angular momentum loss via friction. If the system does not merge, it emerges as a much shorterperiod system consisting of the core of the giant and a rejuvenated main sequence companion. There are several scenarios for the subsequent evolution; the core of the giant may develop into a corecollapse supernova, or the mainsequence companion may evolve into a (sub)giant, causing a second common envelope before it becomes a supernova itself (Tutukov & Yungelson 1993; Podsiadlowski et al. 2002). If precisely one of the stars becomes a supernova and the system is not unbound during this event, the system can develop into a neutron star/black hole – white dwarf/helium star binary. If neither of the components is massive enough to become a corecollapse supernova, a second common envelope can lead to a double white dwarf or a white dwarf – helium star binary. An accretioninduced collapse of a white dwarf may still yield an UCXB, once mass transfer has resumed (van den Heuvel 1984).
More models for UCXB formation exist. In globular clusters, home to about onethird of the presently known UCXB population, formation may be dominated by dynamical interactions (Verbunt 1987; Ivanova et al. 2005; Voss & Gilfanov 2007). After a system with a neutron star or black hole component has formed, UCXBs can evolve through 3 scenarios: via a white dwarf donor (Yungelson et al. 2002), heliumstar donor (Savonije et al. 1986; Yungelson 2008) or a mainsequence donor (Podsiadlowski et al. 2002; van der Sluys et al. 2005; Ma & Li 2009). However, after a relatively short time of mass transfer, each donor becomes degenerate and the subsequent evolution is very similar for each scenario. Therefore, we consider only the white dwarf donor scenario.
As suggested by the rarity of UCXBs, the initial stellar and binary parameters must be finely tuned in order to arrive at a stable ultracompact configuration. Several things can go wrong: binaries can either evolve too slowly to have two postmain sequence components, be too wide to have interacting components, merge during a common envelope, possibly explode as a single degenerate type Ia supernova, become unbound by a supernova, or experience two supernovae.
1.2. Present research
The objective of this paper is to investigate the occurrence of a dynamical instability in low mass ratio UCXBs, and the influence on UCXB evolution of magnetospheredisk interactions and the thermalviscous disk instability. By improving our understanding of all evolutionary stages of UCXBs, we will be able to better estimate their observational properties, such as Xray luminosity and orbital period, fraction of time visible and lifetime. In a forthcoming paper, we will use the results to predict the presentday UCXB population, which can be compared to observations by the Galactic Bulge Survey (Jonker et al. 2011).
In particular the behavior of old, low mass ratio (q < 0.01) and low mass transfer rate systems is unclear. None have been discovered, even though their long evolutionary timescales suggest that many should exist if no disruptive process emerges at some point. Are they invisible to our instruments (most of the time), or do they no longer exist?
Section 2 gives an overview of regular UCXB evolution and the main angular momentum flows in the absence of instabilities. Several complications that may play an important role during the evolution of an UCXB are mentioned in Sect. 2.6. In the method, Sect. 3, we describe the evolutionary tracks as a function of the degree of feedback of angular momentum from accretion disk to orbit, the propeller effect and the disk instability model. In the results, Sect. 4, we analyze the tracks, consider how likely limited feedback is, and look into the magnetic field, both for constant mass transfer and nonsteady behavior such as outbursts. In Sect. 5 we discuss the results and give a conclusion. The appendix contains fits of evolutionary tracks and some analytical approximations.
2. Overview of UCXB orbital mechanics
The evolution of a binary consisting of only degenerate stars or black holes cannot be driven by single star evolution, since such components hardly evolve, if at all, and certainly do not increase their radius in the absence of interaction with other stars. Since we assume white dwarfs to be cold, i.e. thermally relaxed, angular momentum loss and redistribution are the only processes that can alter the system.
2.1. Angular momentum flows
The angular momentum flows in an UCXB are given by (1)where is the change in orbital angular momentum, the loss to gravitational wave emission, angular momentum advected from the donor to the accretion disk along with the transferred matter, is the return of angular momentum from the disk to the orbit by means of a tidal torque between the outer disk and the donor (Lin & Papaloizou 1979; Frank et al. 2002) and is the angular momentum carried by matter ejected from the system. The high surface gravity of degenerate stars prohibits significant wind mass loss from the donor, so spin angular momentum loss via magnetic braking can be neglected. Change in accretor spin can present a sink or source of angular momentum (Marsh et al. 2004), but contrary to white dwarf accretors, this is negligible compared to the angular momentum transport to and from the disk (Priedhorsky & Verbunt 1988) for neutron star and black hole accretors.
When the system has an accretion disk, is precisely balanced by the feedback flow in the opposite direction, (Priedhorsky & Verbunt 1988). In that case, the accretion disk does not gain or lose angular momentum, which leaves orbital angular momentum loss via gravitational wave radiation as the sole process driving binary evolution, giving the usual expression (2)However, below we will consider the more general case. The orbital period of two stars must be rather short for gravitational wave emission to be able to significantly shrink the orbit within the age of the Universe, e.g. shorter than 14 hr in the case of two 1 M_{⊙} components. Once the orbit has shrunk to the point where the Roche lobe of a white dwarf component is smaller than the associated star’s volume, the white dwarf starts to lose surface matter.
2.2. Stable mass transfer
To calculate the mass transfer rate, we assume that the donor precisely fills its Roche lobe at all times t during stable mass loss (in reality the donor slightly overfills its Roche lobe, but this distance is small compared to the Rochelobe radius and can be neglected): (3)where R_{d} is the donor radius and R_{L} is the effective radius of the donor Roche lobe. Equation 3 implies that the ratio of the component masses and the donor radius determines the semimajor axis and orbital period, which increase with time. The change in R_{d} and R_{L} can be separated in a masstransfer dependent and independent part as (Soberman et al. 1997) where ζ_{i} ≡ ∂lnR_{i}/∂lnM_{d} (Hut & Paczyński 1984) represents the change in R_{i} resulting from mass transfer or loss, i.e. a change in donor mass M_{d} (). The subscript evo denotes single star evolution and aml denotes angular momentum loss from the system (in the absence of mass transfer and loss). The advantage of the above separation of the change in radii is that it allows for more easily solving the equilibrium mass loss rate. From Roche lobe geometry and the orbital angular momentum equation follows that (∂lnR_{L})_{aml} = (∂lna)_{aml} = 2(∂lnJ_{orb})_{aml} in the case of fixed component masses (a is the semimajor axis). Combining the above and taking (∂lnR_{d})_{evo} = 0 for the fully cooled, nonevolving, white dwarf case, leads to the (negative) donor mass loss rate^{1}(6)where orbital angular momentum (unrelated to mass transfer or mass loss) is lost exclusively by gravitational wave emission (Eq. (2)), given by Landau & Lifshitz (1975). The effect of angular momentum loss via mass ejection is contained in ζ_{L} and will be discussed in Sect. 3.2. In Fig. 1 we show the mass transfer rate given by Eq. (6) for ζ_{d} given in Sect. 3.3.
Fig. 1
Tracks for a helium white dwarf donor with an initially 1.4 M_{⊙}, 12 km radius neutron star (solid) and 10 M_{⊙} black hole (dashed) accretor. The tracks start at the respective dynamical instability limits of 0.83 M_{⊙} and 1.04 M_{⊙}, since systems can only survive if they start with a lower donor mass (Sect. 2.3). The isotropic reemission limit is ignored in the first part of the track. The numbers correspond to the two circles on their righthand side and represent the logarithm of the system age (yr), which is the time since the onset of mass transfer. The relative location of both tracks will be discussed in Sect. 4.2.2. 
Figure 1 shows how the mass transfer rate decreases with increasing orbital period and time. The white dwarf grows in size as it loses mass, and since it is assumed to fill its Roche lobe at all times (Eq. (3)), the orbital separation and hence orbital period must increase as well. Since gravitational wave radiation becomes significantly weaker in a wider orbit, the mass transfer rate decreases. This can be understood as follows: mass transfer from a less massive to a more massive component in itself tends to widen the orbit, and the weaker the orbitshrinking gravitational waves, the less mass transfer is needed to keep the orbit wide enough to accommodate the donor.
If the binary transfers too much matter, the donor detaches and mass transfer stops. Conversely, if too little matter is transferred, the donor overfills its Roche lobe more and more, increasing mass loss. An equilibrium is reached; the system naturally arrives at a mass transfer rate for which the donor overfills its Roche lobe by the right amount.
2.3. Dynamical stability
From Eq. (6) it can be seen that the mass transfer rate becomes unbounded when ζ_{L} and ζ_{d} approach each other. Dynamical stability of the donor requires that, upon mass loss, the change in R_{L} exceeds the change in R_{d}, i.e.^{2}(7)Dividing both sides by R_{d} = R_{L} (Eq. (3)) and by Ṁ_{d}/M_{d} < 0 yields the well known criterion (e.g. Webbink 1985; Pols & Marinus 1994)^{3}(8)Physically this means that when a donor that does not obey this criterion starts transferring matter via Rochelobe overflow, its volume grows faster than the volume of the Roche lobe it is contained in. The donor will overflow its Roche lobe more and more, leading to runaway mass loss on the dynamical timescale of the donor (~5 s for an 0.6 M_{⊙} white dwarf), and followed by disruption of the donor.
Near the dynamical instability limit, the mass transfer rate becomes in principle arbitrarily high, as shown by the upper left part of the tracks in Fig. 1. This means that in order to find the critical mass ratio at which this instability occurs, we have to consider ζ_{L} for the case in which a fraction 1 of transferred mass is ejected. Figure 2 shows that for a system with an accretor mass M_{a} = 1.4 M_{⊙} (dotted curve), mass transfer is dynamically unstable if M_{d} > 0.83 M_{⊙}. In the case of a 10 M_{⊙} accretor (dashdotted) this becomes M_{d} > 1.04 M_{⊙}.
Fig. 2
Massradius relation exponent ζ_{L} (Sect. 3.2) for a 1.4 M_{⊙} accretor (dotted) and a 10 M_{⊙} accretor (dashdotted) and the zerotemperature helium white dwarf ζ_{d} (dashed, Sect. 3.3). Nearly all matter is assumed to be lost via the isotropic reemission mechanism (Sect. 2.4). 
2.4. Isotropic reemission
Even when the donor mass is low enough to avoid a dynamical instability, the high amount of mass it needs to lose to avoid a merger of the binary components shortly after the onset of mass transfer (van den Heuvel & Bonsema 1984) may be too much for the accretor to accrete or eject. The ultimate condition for survival is how efficiently the accretor can eject matter in the case of a superEddington mass transfer rate. The upper limit of mass that can be removed from the system follows from the energy balance. By assuming that precisely the Eddington limit is accreted, the large amount of gravitational energy liberated can in principle be employed to unbind additional transferred matter (with energy corresponding to the first Lagrangian point) from the system. This mechanism is called isotropic reemission. Soberman et al. (1997) and Tauris & Savonije (1999) described the mechanism as matter transferred via an accretion disk to the vicinity of the accretor, from where it is subsequently ejected as a fast, isotropic wind.
2.5. The deciding stability criterion at the onset of mass transfer
The onset of Rochelobe overflow is a very important event in the evolution of a binary. White dwarf donors that are too massive lose mass in a dynamically unstable way (Sect. 2.3) and the systems containing these do not survive. Figure 3 (upper solid curve) shows that systems with a 1.4 M_{⊙} neutron star accretor are dynamically unstable if the white dwarf donor is more massive than 0.83 M_{⊙}. Since there is no accretion disk present yet when mass transfer first starts, there is no tidal torque transferring back angular momentum to the orbit (to be discussed in Sect. 2.6.1), and the criterion becomes stricter; the system must lie within the lower solid curve in Fig. 3. However, if the donor is less massive than ~ 0.50 M_{⊙} (dotted curve), the disk forms on a timescale short enough to start transferring back enough angular momentum to stop the dynamical instability (Verbunt & Rappaport 1988). Moreover, the 0.50 M_{⊙} disk formation criterion itself is replaced by an even stricter criterion: the isotropic reemission limit of M_{d} ≈ 0.37 M_{⊙} (dashed curve)^{4}. In this paper we assume that isotropic reemission is the criterion which in practice determines whether a system survives the onset of mass transfer, or merges^{5}.
If the binary system survives the onset of mass transfer, it will settle in a state of stable mass transfer, as described in Sect. 2.2.
2.6. Complications
The perpetually expanding binary with declining mass transfer rate described above is idealized. In reality, several processes can potentially disturb continuous accretion or even disrupt the binary, and they may become important in view of the low mass ratios and mass transfer rates expected for compact systems driven by gravitational wave radiation.
2.6.1. Limited feedback of angular momentum at low mass ratio
In Sect. 2.3 the possibility of a dynamical instability at the onset of mass transfer was discussed. According to Fig. 2 such an instability can only occur at a high donor mass of M_{d} ≳ 0.8 M_{⊙}, depending on accretor mass. This is no longer necessarily true if an extra angular momentum sink is introduced. Here we look into possible reasons for limited feedback of angular momentum from disk to orbit.
Angular momentum is advected from the donor to the accretion disk along with transferred matter. Initially, when a disk has yet to form, matter accumulates in a ring around the accretor and, due to viscous friction, spreads out radially to form an accretion disk. The disk keeps expanding until it reaches a radius that is so large that all added angular momentum can be returned to the orbit via the tidal torque between the donor and the outer disk, and no further (see Frank et al. 2002, chap. 5). A torque exists because the disk is slightly elongated due to the tidal force of the donor, and also slightly ahead of the donor because of the shorter orbital period and a delay in adjusting its shape. At the same time, the outwards transport of angular momentum through the disk allows for matter to move inwards, and eventually to be accreted onto the central object.
Fig. 3
Stability regimes in terms of zerotemperature helium white dwarf donor mass versus accretor mass. The upper solid curve is the upper limit for dynamical stability in the presence of an accretion disk that feeds back angular momentum, where all transferred matter is lost via isotropic reemission. The dashed curve is the upper limit for isotropic reemission (i.e. if the system lies above this curve, the accretor will be engulfed by transferred matter). The lower solid curve encloses the parameter space where systems are dynamically stable in the absence of angular momentum feedback. The dotted curve indicates a very rough estimate of the upper donor mass for which there is enough time to form an accretion disk and suppress the nofeedback instability, based on Verbunt & Rappaport (1988). 
Possible causes for reduced feedback include a viscosity that is too low and tidal torque that is too weak (Lin & Papaloizou 1979; Ruderman & Shaham 1983; Hut & Paczyński 1984), and the resumption of mass transfer after a detached stage caused by orbital widening, which in turn is due to mass loss caused by an accretioninduced collapse of a white dwarf accretor (van den Heuvel 1984). Recently, Yungelson et al. (2006); Lasota (2007) stated that it remains unknown what happens to mass transfer and the accretion disk in systems with a mass ratio q = M_{d}/M_{a} ≲ 0.01. When q ≲ 0.02, the circularization radius exceeds the estimates of the outer radius by Paczyński (1977) and Papaloizou & Pringle (1977), causing the disk to be truncated, and matter to circularize onto unstable orbits, thereby preventing accretion (Yungelson et al. 2006; Lasota 2008). The donor may cause gaps in the disk at certain radii via destabilizing orbital resonances, causing disk particles to be confined to certain radius ranges. This process is similar to the formation of gaps in the ring system of Saturn by Mimas and other moons (Franklin & Colombo 1970), and the Kirkwood gaps in the asteroid belt due to Jupiter’s gravitational influence. Matter may be unable to pass these gaps, which would prevent both the outward expansion of the disk to the radius required for full angular momentum feedback to the orbit, as well as the corresponding accretion on the inside.
2.6.2. Occurrence of a dynamical instability at low mass ratio
A dynamical instability at low mass ratio may happen only if feedback would be reduced for any of the reasons given in Sect. 2.6.1. In the absence of feedback and Eq. (1) becomes (9)The disk would turn into an extra angular momentum sink and the orbit loses more angular momentum than in the case of feedback, and has a stronger tendency to shrink. More mass must be transferred to compensate for this and keep the donor within its Roche lobe, and thereby a new equilibrium is established. For systems with a low mass ratio this effect is stronger because of the geometry of the system: mass entering the disk carries a high specific angular momentum because of the large distance (relative to the semimajor axis) of the first Lagrangian point to the center of mass. At very low mass ratio, the orbit loses so much angular momentum to the disk that no amount of mass transfer can compensate for this; the donor will be disrupted^{6}.
Ruderman & Shaham (1983) were the first to identify this instability, and Ruderman & Shaham (1985) showed that in the case of a helium white dwarf donor, the dynamical instability can be reached after 9 Gyr, with M_{d} = 4 × 10^{3} M_{⊙} for M_{a} = 1.4 M_{⊙}. Bonsema & van den Heuvel (1985) and Verbunt & Rappaport (1988) came to a similar conclusion.
On the other hand, Jeffrey (1986) reported that the fast spin period of the single millisecond pulsar PSR 1937+214 contradicts a history in a binary followed by donor disruption, because the spin would have been much slower at the evolutionary stage such a disruption would take place.
Bildsten (2002) considered finitetemperature white dwarfs and found that those would reach the instability at higher donor mass than zerotemperature white dwarfs, for instance at M_{d} = 0.01 in the case of a 10^{6} K helium white dwarf donor and a M_{a} = 1.4 M_{⊙} accretor.
The questions we want to answer are, 1) could feedback be reduced at some stage in the evolution of an UCXB; and 2) if so, does this lead to a dynamical instability within the age of the Universe?
2.6.3. Propeller effect
More processes can potentially limit accretion. In the case of a neutron star accretor, accretion may be disturbed by the magnetic field. This is the propeller effect, in which the fastrotating magnetosphere accelerates transferred matter in the azimuthal direction and stops it from accreting (Davidson & Ostriker 1973; Illarionov & Sunyaev 1975).
2.6.4. Thermalviscous disk instability
At low mass transfer rates, the disk instability model describes a thermalviscous instability resulting from a relatively large and sudden local increase in opacity and viscosity of the disk material (Osaki 1974; Lasota 2001) (note, this disk instability is unrelated to the dynamical instability discussed in Sects. 2.3 and 2.6.2). The high viscosity causes a much higher accretion rate, called an outburst. Outbursts are alternated by lowviscosity stages during which the disk builds up again. Stable behavior persists if the entire disk has a rather homogeneous degree of ionization; UCXBs in particular can have a stable disk if the mass transfer rate is sufficiently high to keep the entire disk ionized through Xray irradiation (in’t Zand et al. 2007).
3. Method
In this section we describe the details of our treatment of the evolutionary stages of UCXBs. The concepts that will play a role during the evolution of UCXBs are discussed in order of appearance.
3.1. Isotropic reemission
As discussed in Sect. 2.4, the energy released by accreting some matter, can unbind additional matter via radiation pressure. The highest spherical accretion rate possible without the arriving matter stopped and blown outwards by radiation pressure is the Eddington accretion limit Ṁ_{Edd} = 4πcR_{a}/κ_{es}, with c the speed of light, R_{a} the accretor radius and κ_{es} ≈ 0.199 cm^{2} g^{1} the electron Thomson scattering mean opacity for hydrogendeficient matter. The maximum mass transfer rate that can be survived via isotropic reemission is given by the ratio of the potential well at the accretor surface Φ_{a} and the potential in the first Lagrangian point (Begelman 1979; King & Begelman 1999; Yungelson et al. 2002). The potential in L1 by definition is equal to the corotatingbinary potential at the Roche lobe surface Φ_{L}, which can be approximated to within 0.15% by (10)with 0 < q < ∞, M_{tot} = M_{a} + M_{d} the total binary mass and G the gravitational constant. If we define φ(q) as the righthand side of Eq. (10), the maximum survivable mass transfer rate is (11)and is independent of R_{a}. Since the matter is ejected from the vicinity of the accretor (Soberman et al. 1997; Tauris & Savonije 1999), it has approximately the same specific orbital angular momentum as the accretor.
3.2. Response of Rochelobe radius to mass transfer
Any effect of change in angular momentum because of mass transfer and ejection from the system is contained in ζ_{L} = ∂lnR_{L}/∂lnM_{d} (Sect. 2.2). The effective radius of the donor Roche lobe R_{L} can be approximated to within 1% by Eggleton (1983)(12)Given q = M_{d}/M_{a} and R_{d}, Eqs. (3) (first part) and (12) yield a. The Rochelobe geometry only applies when the spin period of the Rochelobe filling component is equal to the orbital period, and when the orbit is circular. These conditions are normally met in compact Rochelobe filling systems.
Logarithmic differentiation of R_{L} yields (13)where ζ_{a} ≡ ∂lna/∂lnM_{d} and (14)with (15)where β is the fraction of mass lost by the donor that leaves the system via isotropic reemission, see Eq. (21). ζ_{rL} with β = 0 is the same as in Marsh et al. (2004), who studied mass transfer in binary white dwarfs.
3.2.1. Mass flows
In order to determine ζ_{a}, we consider three mass flows in the binary system: mass lost by the donor can either be accreted, stay in the disk or be unbound from the system.
Accreted mass
The accretion efficiency ϵ (where 0 ≤ ϵ ≤ 1) is the fraction of the matter lost by the donor that is accreted by the primary (16)When we assume full accretion in the case of a subEddington mass transfer rate (in the absence of magnetic fields or instabilities) (17)Ṁ_{Edd} ≈ 3 × 10^{8} M_{⊙} yr^{1} for a neutron star that is accreting hydrogendeficient matter. The evolution of ϵ for two different accretor masses is illustrated in Fig. 4.
Fig. 4
Accretion efficiency ϵ (Eq. (17)) for a helium white dwarf donor and an initially 1.4 M_{⊙} neutron star (solid) and a 10 M_{⊙} black hole (dashed) accretor, in the case of full feedback of angular momentum. Both curves start at their respective dynamical stability limit (Sect. 2.3). 
Mass lost from the system
In the case of nonconservative mass transfer the system loses angular momentum along with the ejected matter. The specific angular momenta of the ejected matter and the binary are and J_{orb}/M_{tot}, respectively. The parameter γ is defined as the ratio between these (18)(note that γ in Soberman et al. 1997 has a different meaning, unlike β and ϵ). It can easily be derived that γ = q in the case of isotropic reemission, because matter has the specific orbital angular momentum of the accretor (Sect. 3.1).
Mass added to accretion disk
In the absence of an accretion disk that can efficiently transfer angular momentum back to the orbit, the stream and torque terms in Eq. (1) do not cancel out, and the disk becomes a net sink of angular momentum. In the extreme case of no feedback at all , leaving (Hut & Paczyński 1984). The circularization radius R_{h} is the orbital distance at which the specific angular momentum is the same as in the incoming stream from the donor (Osaki 1989), minus the amount of angular momentum removed by the donor, approximated by (Verbunt & Rappaport 1988) (19)adjusted for our definition of q = M_{d}/M_{a}, and valid for 10^{3} < q < 1.
In the general case of limited feedback, we introduce a factor η (where 0 ≤ η ≤ 1) describing the fraction of mass lost from the donor that is (net) added to the disk because it cannot be accreted, defined analogous to ϵ in Eq. (16) (20)Conservation of mass is represented by β + ϵ + η = 1. In the case of full feedback (η = 0), β may be estimated by assuming that all nonaccreted matter leaves the system, unless the maximum isotropic reemission rate Ṁ_{mir} > 0 (Eq. (11)) is too low to accomplish this, so (21)
3.2.2. Response of semimajor axis to mass transfer
Combining the above equations, the change in angular momentum of the system because of mass leaving the system or entering the disk is (22)with r_{h} = R_{h}/a.
To find ζ_{a}, we first differentiate the logarithm of the equation for orbital angular momentum to time. We introduce a disk mass term to allow for transferred matter to stay in the disk instead of being accreted, corresponding to possible limited feedback to find (23)where M_{disk} ≪ M_{a} has been used. Mass added to a disk around the accretor has the same effect on J_{orb} as mass that is accreted, because the accretordisk subsystem is orbiting the center of mass of the binary. Hence the generalization of accretor mass to include disk mass.
Equating Eqs. (22) and (23) eliminates . Inserting Eq. (16), Ṁ_{tot} = βṀ_{d} and Ṁ_{a} + Ṁ_{disk} = (β − 1)Ṁ_{d} gives (24)In the case of full feedback, η = 0 and β + ϵ = 1. At the onset of mass transfer from a massive donor, when a disk has not yet formed, η = 1 and , identical to the equation in Verbunt & Rappaport (1988); Marsh et al. (2004). In the case of lowmass donors where accretion is made impossible by limited feedback, accretion has ceased (ϵ = 0) and isotropic reemission has stopped accordingly (β = 0), hence η = 1. The actual value of η during the evolution of an UCXB will be considered in Sect. 4.1.
The ζ_{L} relation as a function of q for ϵ = 0 is shown in Fig. 2 for two accretor masses and for more general mass flow parameters in Fig. 7, for only a neutron star accretor. The Rochelobe radius depends on both component masses and angular momentum losses via its dependence on a, and on the mass ratio via R_{L}/a.
Note that the variables ϵ, Ṁ_{d}, ζ_{L}, , M_{d}, M_{a} and Ṁ_{Edd}, and in the case of a black hole accretor also R_{a}, are circularly related to each other through Eq. (6), the gravitational wave equation, and Eqs. (17) and (13). Hence, iteration at each mass step is necessary to arrive at the equilibrium values.
3.3. White dwarf radius
Now ζ_{L} is known, only the response of the white dwarf donor radius to mass loss, ζ_{d}, is still needed to evaluate Eq. (6). We justify the use of a zerotemperature white dwarf massradius relation by noting that thermal pressure is negligible for significantly degenerate objects such as white dwarfs with a mass above 0.01 M_{⊙}, and still quite small for lower masses. From the massradius relations derived by Deloye & Bildsten (2003) it follows that the effect of ideal gas pressure on the radius, additional to degeneracy pressure and Coulomb attraction, is less than 10% for M_{d} > 10^{2} M_{⊙} when the temperature T = 10^{6} K, less than 2.3% for M_{d} > 4 × 10^{3} M_{⊙} when T = 10^{5} K and less than 0.1% for M_{d} > 10^{3} M_{⊙} when T = 10^{4} K. Temperature is important only for very lowmass or very hot donors. For a mass ratio below ~ 0.01, donor radii could be 10 − 30% higher than for the zerotemperature case (Nelson & Rappaport 2003). Furthermore, we consider only degenerate donors, since any fusion will be extinguished on a relatively short timescale after some mass has been lost (Savonije et al. 1986; Yungelson 2008). Also, a nondegenerate outer layer can have a significant effect on the radius response to mass loss, but again this is temporary – once the layer has been lost, a degenerate object remains.
The radius R_{d} of a zerotemperature white dwarf (the ratio of atomic number to atomic weight Z/A = 1/2) of mass M_{d} by Eggleton (Rappaport et al. 1987) is (25)with M_{Ch} = 1.44 M_{⊙}. For a helium white dwarf, M_{p} = 5.66 × 10^{4} M_{⊙} and for a carbonoxygen white dwarf (assuming 50% of each element) M_{p} = 1.44 × 10^{3} M_{⊙}, valid for 0 < M_{d} < M_{Ch} (Marsh et al. 2004). ζ_{d} ≡ ∂lnR_{d}/∂lnM_{d} is determined by logarithmic differentiation (26)and is illustrated by the dashed curve in Fig. 2. Throughout this paper we consider helium white dwarfs, though all results except quantitative details and the disk instability details also apply to carbonoxygen composition.
3.4. Propeller effect
3.4.1. Magnetic field of neutron star accretors
Neutron stars typically possess a strong magnetic field. During accretion, the magnetic field strength is reduced (Bhattacharya & van den Heuvel 1991), but even for old millisecond pulsars with a spin period below 30 ms, it can still be 10^{8−9} G (Wang et al. 2011). The magnetic field which is locked to the neutron star is supposed to penetrate the accretion disk, and in the case of a rapidly spinning neutron star, the outer magnetosphere forces the orbiting matter to velocities exceeding the local Kepler velocity by exerting a torque. As a result of the angular momentum gained, disk material in this region moves outwards, which reduces accretion and can lead to unbinding the matter from the system.
3.4.2. Critical mass transfer rate
The Alfvén radius (also known as magnetosphere radius) is the distance from the neutron star below which the magnetic energy density (i.e. magnetic pressure) dominates the ram pressure in the accretion disk (Lamb et al. 1973). (27)with B = B_{NS} (R_{NS}/R)^{3} the magnetic induction of the neutron star (NS) in vacuum, B_{NS} the stellar surface field at the magnetic equator, ρ the average density at radius R and the Keplerian orbital velocity.
We use the spherical continuity equation (see justification in Sect. 3.4.3) (28)with the freefall velocity, to eliminate ρ and solve for the Alfvén radius (Davidson & Ostriker 1973; Elsner & Lamb 1977) (29)where is the magnetic dipole moment of the neutron star.
Since the magnetic field lines are dragged along with the neutron star’s rotation, their orbital angular frequency equals the spin angular frequency ω_{s} of the neutron star. At a distance R from the neutron star, the corotating field lines have an orbital velocity ω_{s}R. This equals the Keplerian velocity at the corotation radius (30)The propeller effect acts only if the magnetic field energy dominates the kinetic energy of the disk beyond the corotation radius, and when the orbital velocity of the field lines exceeds the escape velocity at the Alfvén radius (Rappaport et al. 2004), i.e. , where v_{K,μ} is the Kepler velocity at R_{μ}, so which implies R_{μ} > 2^{1/3}R_{co}. Inserting Eqs. (29) and (30) yields the propeller criterion (31)
3.4.3. Disk accretion
By using the thin disk description by Dunkel et al. (2006) it follows that the magnetic field is far too weak to dominate the gas flow in the disk even for the lowest expected accretion rates. Even for the slightly unrealistic case of a neutron star with a 10^{9} G equatorial magnetic field and a 1 ms spin period, the magnetic field only dominates the kinetic energy of the Kepler flow at the corotation radius if the mass transfer rate is lower than 10^{16} M_{⊙} yr^{1}, ~ 10^{10} times lower than for the spherically infalling case. High density blobs of gas can also dominate the magnetic field and be accreted (Aly & Kuijpers 1990). Such blobs could form if matter piles up near the Alfvén radius, held back by the magnetic field, and enter once enough matter has accumulated.
The conclusion is that Eq. (29) is an upper limit for the Alfvén radius since the accretion is not spherical as assumed in Eq. (28) but occurs in a disk. However, since the density and the geometry of the interacting diskmagnetosphere region is poorly known, we choose to allow for disk accretion in the equations by also considering an effective magnetic field strength at the stellar surface which is smaller than the real field strength B_{NS}.
3.5. Disk instability model
in’t Zand et al. (2007) make a rough estimate for the lowest stable mass transfer rate for helium composition based on Lasota (2001) and Menou et al. (2002)(32)where P_{orb} is the orbital period. Once the mass transfer rate decreases below this value, a thermalviscous instability appears (Sect. 2.6.4) and the disk experiences outbursts. The timescale on which an accretion disk rebuilds after a collapse is equal to M_{disk}/Ṁ_{d}, where the disk mass is given by (33)with Σ the surface density for zone C in Dunkel et al. (2006). R_{in} and R_{out} are the inner and outer disk radii and will be estimated later.
The disk instability model is still not fully understood, as illustrated by the case of SS Cygni, a dwarf nova that according to the disk instability model is too bright and hot to have an unstable disk (Schreiber & Lasota 2007). This increases the uncertainty in the mass transfer rate given by Eq. (32). On the other hand, the intermediate state in AM CVn systems seems to be predicted quite well (Tsugawa & Osaki 1997).
4. Results
4.1. Feedback of angular momentum
Since η is the most uncertain mass flow parameter in Sect. 3.2.1, here we first look into the likelihood of limited feedback of angular momentum at some low mass ratio. If the torque between outer disk and donor is too weak, a dynamical instability may result (Sect. 2.6.1).
Fig. 5
Important radii within the accretion disk; the distance between the accretor and the L1 point d_{L1} (solid), the average equatorial Roche lobe radius (i.e. in the orbital plane) (dashed), the circularization radius R_{circ} (dashdotted) and R_{circ} corrected for donor influence during freefall, R_{h} (dotted, Eq. (19)). Furthermore, the pluses (+) indicate the outer tidally stable radius from Papaloizou & Pringle (1977). The stars ( ⋆ ) indicate the largest orbit that does not intersect others, according to Paczyński (1977) (which lies near the 3:1 resonance). The two open circles are outer disk estimates by Priedhorsky & Verbunt (1988), and the three filled triangles are outer disk radii from our SPH simulations. 
Also mentioned in Sect. 2.6.1 was that Yungelson et al. (2006) and Lasota (2008) have suggested that a resonance between the orbital period and the (mean) period of particles in the disk may prevent the accretion disk from redistributing angular momentum. Also, when q ≲ 0.02, the circularization radius exceeds the estimates of the outer radius by Paczyński (1977) and Papaloizou & Pringle (1977) (~0.48a, see the stars and pluses in Fig. 5) which could hinder accretion, and at the same time reduce the feedback of angular momentum, making the binary less stable. Their value of q = 0.02 can be refined by using the adjusted circularization radius equation by Verbunt & Rappaport (1988) (dotted curve in Fig. 5), who correct the circularization radius R_{circ}/a = (1 + q)(d_{L1}/a)^{4} (Frank et al. 2002) (dashdotted) for the effect of the donor on the stream of matter. During this process, angular momentum is extracted from the stream, so R_{h} < R_{circ}. Using R_{h} (Eq. (19)) instead of R_{circ} lowers the critical mass ratio to q = 0.01 (Fig. 5). Figure 5 also shows that the location of the first Lagrangian point (solid curve) is very close to the donor at low mass ratio, and that the equatorial Roche lobe radius (dashed) also approaches the L1 point, because the equatorial accretor Roche lobe becomes more circular rather than teardropshaped at low mass ratio^{7}.
Fig. 6
Numerical smoothed particle hydrodynamics simulation of an accretion disk in dynamical equilibrium for mass ratio 10^{3} after 400 orbits. The solid curve represents the equatorial Roche lobes. The accretor has been marked by a small open circle. The dashed circles indicate the 3:1 (inner) and 2:1 (outer) resonances with the orbital period. 
We ran numerical smoothed particle hydrodynamics (SPH) simulations such as that in Fig. 6 which show that disks function properly even at a very low mass ratio. The figure shows a disk in dynamical equilibrium after 400 orbits (see Simpson & Wood 1998 and Wood et al. 2009 for details of the numerics). The simulations show that the disk fits well inside the Roche lobe (triangles in Fig. 5). This confirms the conclusion of Priedhorsky & Verbunt (1988), who used an accretion disk model to find that the disk expansion is not stopped at 0.48a and that the disk stays within the accretor Roche lobe (circles in Fig. 5) and tidal torques are strong enough to transfer back all angular momentum. Also, nothing special appears to happen in the disk near the most important resonant radii, the circles in Fig. 6. Hence, η remains zero at very low mass ratio. In the remainder of the paper we will consider the fullfeedback case, except for Sect. 4.2.3, where we will show evolutionary tracks for the case in which feedback somehow stops at a low mass ratio.
Observations of binaries cannot rule out the existence of a dynamical instability since no systems with an upper mass ratio below 0.01 have been discovered, although a couple have minimum mass ratios below this limit (Galloway et al. 2002; Krimm et al. 2007; Altamirano et al. 2010).
4.2. The evolution of ultracompact Xray binaries
The key elements in understanding the evolution of UCXBs are represented by Eq. (6). Mass transfer is driven by gravitational wave radiation, and also depends on the relative responses of donor and Roche lobe size to mass transfer. In Sect. 4.1 we found that advected angular momentum will be returned from disk to orbit for any mass ratio, which constrains ζ_{L} and allows us to solve all relevant binary parameters during the entire evolution when assuming stable accretion, rather than in outbursts. A system can survive the onset of mass transfer only if the accretor can eject all transferred matter exceeding the Eddington limit via isotropic reemission. For any realistic stellarmass accretor, a merger can be avoided if the white dwarf donor has a mass below ~0.3 M_{⊙}. Mass transfer is dynamically stable if, upon mass loss, the donor expands less than its Roche lobe (Eq. (8)), that is, when the solid curve in Fig. 7 lies below the dashed curve. UCXBs with black hole accretors can theoretically exist, but may be rarer than neutron star UCXBs due to the steep slope of the initial mass function. However, the binary survival rate during the supernova event can be different for black hole formation than it is for neutron star formation due to a different kick velocity distribution.
4.2.1. Evolution in the case of full feedback
Fig. 7
Massradius relation exponent ζ_{L} (Sect. 3.2) for various modes of mass transfer from a zerotemperature helium white dwarf to a neutron star of initially 1.4 M_{⊙} with a 12 km radius. The lower dotted curve shows the case of accretion efficiency ϵ = 0 and full feedback (η = 0) (almost all transferred matter lost via isotropic reemission, hardly any accretion). The middle dotted curve is for ϵ = 1 and η = 0 (conservative mass transfer). The upper dotted curve is for no feedback (η = 1). The zerotemperature white dwarf exponent ζ_{d} (Sect. 3.3) is represented by the dashed curve. The solid curves starting at q = 0.59 (M_{d} = 0.83 M_{⊙}, Sect. 2.3) indicate the actual value of ζ_{L} (for ϵ shown in Fig. 4) for both the full feedback and nofeedback case. Until q = 0.01 there is full feedback in either case. Below q = 0.01 the full feedback case is shown by the lower solid branch, where the open circle denotes an age of 10 Gyr (since the onset of mass transfer). The nofeedback case has a jump in ζ_{L} at q = 0.01, after which it follows the upper dotted curve until the intersection with the dashed curve (filled circle) where donor disruption occurs at a system age of 300 Myr. 
In Fig. 8, the mass transfer rate starts very high because of the massive donor and especially because of the short orbit. The first part of the tracks lies above the isotropic reemission limit (dashed, Sect. 2.4), which means that a system cannot survive here, this part of the tracks is shown just for illustration. Isotropic reemission of superEddington mass transfer is possible for donor mass below 0.37 M_{⊙} and orbital period above 1.8 min. For lower donor mass (0.08 M_{⊙}, with a corresponding orbital period of 8 min), the mass transfer rate decreases below the Eddington limit (dotted), and from here on conservative mass transfer is possible (Fig. 4). At later times, the mass transfer rate reaches the thermalviscous disk instability rate (dashdotted, Sect. 3.5). The donor mass is about 0.02 M_{⊙} at this point, and the orbital period 28 min in the case of a helium disk. Systems with a longer period than this do no longer have stable accretion disks. After 10 Gyr of mass transfer, a system with a 1.4 M_{⊙} neutron star accretor would have a donor mass of 4.2 × 10^{3} M_{⊙} and an orbital period of 85 min. In the case of a 10 M_{⊙} black hole accretor, the donor mass would be 2.6 × 10^{3} M_{⊙} and the orbital period 111 min. (The 10 Gyr given here, as well as ages given elsewhere in the paper, is the time measured after the onset of mass transfer to the neutron star or black hole.) These systems technically do not fall under the present (observational) definition of UCXBs (an orbital period ≲ 60 min). Gravitational wave radiation is stronger in systems with a more massive accretor, therefore the black hole track in Fig. 8 lies above the neutron star track (to be explained in Sect. 4.2.2).
Fig. 8
Mass transfer rates for a zerotemperature helium white dwarf donor with an initially 1.4 M_{⊙}, 12 km radius neutron star accretor (lower solid curve) and a 10 M_{⊙} black hole accretor (upper solid curve), starting at the dynamical instability limit (i.e. ignoring possible merging), for full feedback. The dashed curves branching away from these represent the evolution assuming no feedback below q = 0.01 (Sect. 4.1). The numbers correspond to the two circles on their lefthand side and represent the logarithm of the system age (yr), which is the time since the onset of mass transfer. Indicated are the isotropic reemission regime (dashed, Eq. (11)), the Eddington limit (dotted) and the accretion disk instability limit by in’t Zand et al. (2007) (dashdotted, Sect. 3.5) Note that the dashed neutron star curve at its local maximum “touching” the upper dashdotted curve is a coincidence. 
Figure 9 is a linear graph showing more clearly how the orbital period increases with time, initially the period increases very rapidly, and slowly later on. Systems with a massive accretor reach a given orbital period earlier on, because of their shorter evolutionary timescale.
Fig. 9
Orbital period against time for a zerotemperature helium white dwarf donor with an initially 1.4 M_{⊙}, 12 km radius neutron star (lower solid) and 10 M_{⊙} black hole (upper solid) accretor. 
The change in accretor mass during evolution
For young systems with massive donors, the orbit is so compact that the donor is forced to lose mass at a high rate in order to avoid a merger. For these systems ϵ ≈ 0, because most mass is expelled from the system by isotropic reemission. Once the transfer rate drops below the Eddington limit, ϵ = 1 from Eq. (17) as all transferred mass can be accreted (Fig. 4). This happens at a donor mass of 0.08 M_{⊙}. Because ϵ goes from negligible to 1 quickly (in terms of donor mass), the total mass increase of a neutron star ΔM_{a} will not exceed − ^{∫}ϵṀ_{d} dt = 0.08 M_{⊙} by Eq. (16). The angular momentum associated with ΔM_{a} is relevant when considering spin up of the accretor, in Sect. 4.3.
Observed potential longperiod UCXBs
The binary millisecond pulsar HETE J1900.1–2455 has an orbital period of 83.3 min, but a donor much more massive than systems of the same period following the evolution described here. The donor is probably a brown dwarf of 0.016 − 0.07 M_{⊙} (Kaaret et al. 2006).
Because of their shorter evolutionary timescale, white dwarf – black hole systems can potentially reach an orbital period of 110 min within the age of the Universe, if the system is not disrupted at some low mass ratio and mass transfer starts relatively quickly after the Big Bang. There are 2 lowmass Xray binaries known with 90 − 115 min periods, viz. 1E 1603.6+2600 at 111 min (Morris et al. 1990) and XTE 1748–361 at 97 min (Bhattacharyya et al. 2006). Another 2 lowmass Xray binaries have periods just outside this range, SAX J1808.4–3658 at 120 min (Chakrabarty & Morgan 1998) and GS 1826–238 at 126 min (Homer et al. 1998), although the last period is uncertain (Liu et al. 2007). All these systems exhibit type I Xray bursts so they must harbor a neutron star. The donors are not expected to be cold white dwarfs since these long periods can not have been reached by evolution on the gravitational wave timescale.
Why is conservative mass transfer less stable than nonconservative mass transfer?
ζ_{L} for conservative mass transfer (ϵ = 1) is higher than for fully nonconservative mass transfer (ϵ = 0) (the lower two dotted curves in Fig. 7) implying dynamically less stable mass transfer for ϵ = 1 by Eq. (6), when we assume isotropic reemission and full feedback of angular momentum from disk to orbit (η = 0). This may seem counterintuitive since less mass ejection makes a system more stable because less angular momentum is lost. The explanation is that the latter effect (represented by the γterm in ζ_{a}, Eq. (24)) is outweighed by the effect of the (1 − β)q ( = ϵq here) term for all ϵ and q due to the unfavorable effect on the accretor mass: in the case of nonconservative mass transfer the accretor gains less mass, implying a must increase more to compensate for this () than in the conservative case. This effect is enhanced by the increase of ζ_{rL} with ϵ (Eq. 14). For ϵ = 1, the mass ratio q = M_{d}/M_{a} decreases faster than for ϵ = 0. The more extreme q causes the donor Roche lobe to shrink more (or expand less) upon mass loss relative to a.
4.2.2. Why do massive accretors speed up evolution?
Figures 1 and 8 show that UCXBs evolve faster if they have a more massive accretor. This is because evolution is driven by gravitational wave radiation. Here we give the explanation and look into some consequences.

The effect of accretor mass on mass transfer rate:
For q ≲ 1, the orbital period of a Rochelobe filling system is only weakly dependent on the accretor mass (which can be seen from the Rochelobe approximation by Paczyński 1971), so the semimajor axis for fixed donor mass, according to Kepler’s third law. From Eq. (6) and (Landau & Lifshitz 1975) follows that , so the mass transfer rate increases with accretor mass for a fixed donor, even though the orbital separation increases as well. This explains why the black hole track in Fig. 8 lies above the neutron star track. The correction from ζ_{d} − ζ_{L} is small unless ζ_{d} − ζ_{L} approaches zero, as evidenced by the intersection of the tracks in Fig. 8, where the black hole system is still further away from its dynamical stability limit (1.04 M_{⊙}). In reality, this intersection is not important since systems with an orbital period below ~2 min cannot survive.

Systems of the same age after the onset of mass transfer: When we compare systems of the same age, mass transfer is lower in the case of a more massive accretor (again, except when they are extremely young). The influence of the lower donor mass outweighs the influence of the more massive accretor in the gravitational wave equation, causing a lower mass transfer rate for more massive accretors of a fixed age. An old 1.4 M_{⊙} neutron star accretor system has a mass transfer rate that is ~50% higher than a 10 M_{⊙} black hole accretor system of the same age, as shown by the fits in Appendix A. (If the systems did not start with identical donors, the more massive donor will lose mass fast enough to (almost) catch up with the lower mass donor, on a timescale much shorter than the typical age of the systems.) A system with a lower accretor mass consistently experiences a higher mass transfer rate, but still its donor remains more massive. It can never catch up, because its mass transfer rate would drop below that of the black hole system in the case of nearly equal donor masses. Lastly, the gravitational wave timescale for a system with a neutron star accretor is ~7% longer than that of a 10 M_{⊙} accretor system of the same age, and the orbital period is ~25% shorter.
To summarize, in the case of an identical donor, Ṁ_{d}/M_{d} and Ṁ_{d} are higher for systems with a higher accretor mass. For systems of the same age, however, M_{d} and Ṁ_{d} are lower for a higher accretor mass, whereas Ṁ_{d}/M_{d} is marginally higher (i.e. mass is transferred on a slightly shorter timescale).

Consequences for AM CVn systems: The above also implies that UCXBs evolve on a shorter timescale than AM CVn systems (which have white dwarfs both as donor and accretor) with an identical donor. The higher value of ϵ in AM CVns makes the mass transfer dynamically less stable and the mass transfer rate higher, but the accretor mass has a larger impact, unless the donor mass is close to the dynamical stability limit. For a low donor mass, ζ_{d} − ζ_{L} is almost independent on the mass ratio and ϵ, so the mass transfer rate only depends on the accretor mass. Also, AM CVns are expected to experience a higher mass transfer rate than UCXBs of the same age after onset of mass transfer, except again for very young systems. Old AM CVns have a more massive donor than equally old UCXBs.
4.2.3. Evolution in the case of no feedback at low mass ratio
In Sect. 2.6.1 we mentioned several proposed mechanisms for reducing feedback of angular momentum from the disk to the orbit, but we think that none of these will actually work (Sect. 4.1). In the case that after all a mechanism exists that leads to a reduction in or absence of feedback in systems below a certain low mass ratio, here we briefly discuss the evolution of such systems, including the dynamical instability that may follow. We consider the case in which all feedback stops once q < 0.01, i.e. η = 1 and β = 0 in Eqs. (24) and (15) for these mass ratios.
Under this condition, the evolution of ζ for a neutron star accretor UCXB is shown by the upper solid curve in Fig. 7. The sudden jump in ζ leads to an increased mass transfer rate per Eq. (6), which speeds up evolution, therefore the low mass ratio where ζ_{d} and ζ_{L} meet can be reached well within the age of the Universe. Two evolutionary tracks are shown as the dashed curves branching away at q = 0.01 from the solid curves in Fig. 8.
Dynamically stable parameter ranges
Fig. 10
Stability regimes in terms of zerotemperature helium white dwarf donor mass versus accretor mass. The upper solid curve is the upper limit for dynamical stability in the presence of an accretion disk that feeds back angular momentum, where almost all transferred matter is lost via isotropic reemission. The upper dashed curve is the upper stability limit for isotropic reemission and the lower dashed curve indicates the proposed q = 0.01 lower mass ratio limit for feedback (Sect. 4.1). The lower solid curve encloses the parameter space where systems are dynamically stable in the absence of angular momentum feedback; all matter enters and stays in the disk. The dashdotted lines indicate when the propeller effect comes into play according to Eq. (31) for several neutron star spin periods and an equatorial magnetic field of 10^{8.5} G. The dotted isochrones indicate the donor masses for system ages (since the onset of mass transfer) of 0.1,1,10 Gyr based on full feedback. The dynamical stability and isotropic reemission rates are independent of accretor radius and type, hence these curves are continuous. Moreover, the dynamical stability rates are independent of accretor type and therefore also hold for white dwarf accretors. 
Figure 10 shows which combinations of donor and accretor masses are dynamically stable, both in the case of feedback and the case of no feedback. The accretor mass has a large influence on when an instability occurs, if at all. The bottom left region enclosed by the solid line is dynamically stable for mass transfer, even in the absence of feedback. Initial mass transfer is dynamically unstable for every M_{d} if M_{a} > 11 M_{⊙}, because of high specific angular momentum of the stream in low mass ratio systems. The upper border of this region shows that the highest donor mass that can possibly avoid dynamical instability in the case of no feedback is 0.42 M_{⊙}, for an accretor mass of 5.3 M_{⊙}^{8}. As explained in Sect. 2.5, in practice the upper border will be replaced by a limit of ~0.5 M_{⊙}. Note that the donor mass at which the low mass ratio dynamical instability could occur is very sensitive to accretor mass near M_{a} = 1.4 M_{⊙}. At exactly which donor mass this instability occurs depends on the massradius relation of the donor, which in turn depends on its temperature and composition.
The upper limit for dynamical stability in the case of full feedback is represented by the upper solid curve in Fig. 10. Because ζ_{L} (Eq. (13)) has a minimum of −5/3 for q → 0 regardless of ϵ and γ, any system with a zerotemperature white dwarf donor mass exceeding 1.07 M_{⊙} is dynamically unstable, even for arbitrarily high accretor mass^{9}.
Age of disruption
If feedback would suddenly stop at q = 0.01, a dynamical instability would happen immediately if M_{a} > 4.7 M_{⊙}, and later for less massive accretors (shown by the intersection between the lower solid and dashed curves in Fig. 10). For an UCXB with a 1.4 M_{⊙} accretor, the instability would happen at an age of 300 Myr, for M_{d} = 3.2 × 10^{3} M_{⊙} (the upperleft intersection at q = 2.3 × 10^{3} in Fig. 7). The reason such a low donor mass could be reached well within 10 Gyr is the increased mass transfer rate caused by the extra angular momentum sink, about one order of magnitude for a 1.4 M_{⊙} accretor (Fig. 8). A 10 M_{⊙} black hole system would reach q = 0.01 after already 80 kyr.
Disruption of lowmass donors would happen later than 10 Gyr if M_{a} < 0.6 M_{⊙}, this could apply to AM CVn systems. UCXBs with any realistic accretor mass would suffer disruption well within 10 Gyr.
4.3. System evolution during the propeller phase
Even though most likely no dynamical instability occurs at low mass ratios, the regular accretion rates as described in Sect. 4.2.1 are idealized. Here we investigate the UCXB evolution taking a magnetized neutron star accretor into account, using angular momentum and energy considerations. This gives insight into the feasibility of the propeller effect (Sect. 3.4) removing matter from the system and hence suppressing accretion onto and Xray emission from UCXBs.
4.3.1. Steady mass transfer
First we assume the absence of disk instability (Sect. 3.5). This means that matter approaches the magnetosphere at a constant rate (that is, only varying on the evolutionary timescale). From a magnetic perspective, the evolution can be divided in a spin up and a spin down stage.
Spin up
Initially the mass transfer rate is so high that the kinetic energy in the disk completely dominates the magnetic field energy; the Alfvén radius lies either at the surface or just outside the neutron star. Moreover, the relatively low initial accretor spin frequency (and correspondingly large corotation radius) ensures that any interaction between disk and magnetic field that is still present, causes the neutron star to spin up rather than matter to be expelled, because the fast rotating disk applies a torque on the field lines, which leads to transfer of angular momentum from disk to neutron star. We assume that as long as the Alfvén radius is smaller than 2^{1/3} ≈ 1.26 times the corotation radius (Sect. 3.4.2), all arriving matter is accreted, and the amount of angular momentum added to the accretor corresponds to the specific angular momentum of disk matter at the Alfvén radius. Specific angular momentum of matter in a Keplerian orbit is proportional to the square root of the distance to the accretor, and we will find that the typical ratio of the outer disk radius to the Alfvén radius exceeds 10^{3} during the spinup stage. Therefore, the fraction of angular momentum that is not transferred back via the outer disk is at most a few percent, and our initial assumption of dynamically stable mass transfer holds.
Fig. 11
Radii in UCXB accretion disk versus orbital period for a 1.4 M_{⊙} neutron star accretor with a residual equatorial magnetic field of 10^{8.5} G. The shaded layer at the bottom represents the neutron star surface, at a radius of 12 km. The curve annotated L1 shows the distance of the first Lagrangian point to the neutron star, which is an upper limit on the outer disk radius. The circularization radius by Verbunt & Rappaport (1988) is discussed in Sect. 4.1. The two upper dashed curves indicate the radii where the temperature of a helium disk (Dunkel et al. 2006) equals the He I ionization temperature, for two values of the viscosity parameter α, as discussed in Sect. 2.6.4. The solid corotation radius curve applies to the full propeller case, the dashed one to the hypothetical full accretion case. The speedoflight cylinder radius and Alfvén radius are also shown. 
As the neutron star spins up, the corotation radius approaches the neutron star. At the same time the Alfvén radius moves outwards as a result of decreasing mass transfer. At an orbital period of 10 min (and a system age of 1 Myr) these two meet. The surface magnetic field strength at this time is determined by how much mass the accretor has gained up to this point (Sect. 4.2.1). We adopt the empirical relation for the accretioninduced field decay by Shibazaki et al. (1989), with a residual field strength of B_{NS} = 10^{8.5} G (Zhang & Kojima 2006). We also assume that the neutron star does not spin when the mass transfer starts, i.e. J_{NS} = 0, since the earlier commonenvelope stage is too shortlived to allow for significant accretion. Also, the strong initial magnetic field suggests a very low initial spin frequency relative to the eventual spin frequency.
The amount of angular momentum added to the neutron star by accretion is (34)where R_{in} is the inner disk radius, which determines the specific angular momentum of matter in a Keplerian orbit at this radius. We assume that matter stops losing angular momentum outwards once it moves inside this radius. R_{in} = R_{μ} (Li & Wang 1999) except when the magnetic field does not dominate the Kepler flow anywhere, in which case the stellar surface is taken, or when R_{μ} lies outside the speedoflight cylinder radius R_{c} = c/ω_{s}, where corotating field lines would travel at the speed of light. There, field lines are open and particles can escape by following them outwards. The value of the inner radius is summarized as (35)From J_{NS} the spin angular frequency ω_{s} is found using J_{NS} = I_{NS}ω_{s}, where is the moment of inertia of a 1.4 M_{⊙}, 12 km neutron star (Lattimer & Schutz 2005). We consider angular momentum accretion rather than energy accretion, because the latter is complicated by radiation losses and orbit coupling.
Spin down
Once the Alfvén radius moves beyond the corotation radius, the torque is reversed, and matter in Keplerian orbits between the corotation radius and the Alfvén radius may be accelerated by the superKeplerian velocity of the field lines, reducing the accretion rate. Matter near R_{μ} (or R_{c}, once smaller) that reaches the escape velocity then leaves the system on trajectories outside the orbital plane, the angle between the outflow and the plane being smaller for a stronger propeller effect (Lovelace et al. 1999; Ustyugova et al. 2006).
Due to its accretion history, the neutron star has gained a large amount of rotational kinetic energy at this stage. The evolution of spin frequency ω_{s} once R_{μ} > 2^{1/3}R_{co} depends on whether there is enough energy, and whether it can be employed to unbind all or some fraction of arriving matter. At this point we assume that there is sufficient energy to unbind all transferred matter when R_{μ} > 2^{1/3}R_{co}, and will later find that this is indeed the case. Then, the energy of the neutron star changes as (36)where dE_{NS} and Ṁ_{d} are negative, and the factor 2 appears because the total energy in a Keplerian orbit is half the potential energy. The rotational kinetic energy of the neutron star is , which yields ω_{s}.^{10}
Spin evolution
Figure 11 illustrates the evolution of the accretion disk, magnetosphere and corotation radius for an UCXB. The dotted part of the Alfvén radius curve (Eq. (29)) lies beyond the speedoflight cylinder radius, and therefore the field lines cannot corotate in this region. The corotation radius is given by Eq. (30) and is proportional to , where the accretor spin period P_{s} = 2π/ω_{s}. The solid corotation curve shows the evolution when no matter is accreted once R_{μ} > 2^{1/3}R_{co}, but instead all transferred matter is unbound from radius R_{in}, extracting energy from the rotation of the neutron star. For comparison, the dashed corotation curve shows the spin evolution in the hypothetical case in which, once R_{μ} > 2^{1/3}R_{co}, all matter would continue to be accreted with specific angular momentum corresponding to R_{in} instead of stopped by the magnetosphere. Before R_{μ} = 2^{1/3}R_{co}, both assume full accretion so there these two are equal.
Because the corotation radius stays within the Alfvén radius, apparently the early accretion at a high rate has sufficiently spun up the neutron star to provide enough energy to unbind any fraction of transferred matter during the whole remainder of the evolution, when the mass transfer rate is low. This justifies our assumption above Eq. (36).
By using a lower effective equatorial surface magnetic field strength to accommodate for disk accretion, as discussed in Sect. 3.4.3, less angular momentum is accreted per unit mass due to the smaller Alfvén radius, but also accretion can continue for a longer time because the propeller effect appears at a later stage. The result is that the neutron star is spun up to a lower spin period, while subsequently less energy is needed to unbind arriving matter because the propeller effect starts at a lower mass transfer rate. The eventual spin period is lower in the case of a lower effective magnetic field.
Figure 12 shows the spin period directly, rather than the related corotation radius. The hypothetical case in which the propeller effect does not work, so all matter is accreted at the stellar surface carrying the local Keplerian specific angular momentum, is shown by the dashdotted curve (Sect. 3.4.3). For comparison we include the observed UCXB neutron star spin periods, which are too small in number and have too scattered spin periods to have any constraining value to our model.
Fig. 12
Spin period versus orbital period for a 1.4 M_{⊙} neutron star accretor with a residual equatorial magnetic field of 10^{8.5} G, assuming all transferred matter is propelled out of the system once R_{μ} > 2^{1/3}R_{co} (solid, Sect. 3.4.2), or hypothetically keeps being accreted (dashed), or hypothetically is accreted at the neutron star surface all the time (dashdotted). Observed millisecond pulsarUCXB parameters are represented by filled circles; XTE J1807–294 at an orbital period of 40.07 min (Markwardt et al. 2003a) and a spin period of 5.25 ms (Markwardt et al. 2003b), XTE J1751–305 at 42.42 min, 2.30 ms (Markwardt et al. 2002), XTE J0929–314 at 43.58 min, 5.40 ms (Galloway et al. 2002), SWIFT J1756.9–2508 at 54.70 min, 5.49 ms (Krimm et al. 2007) and NGC 6440 X2 at 57.3 min, 4.86 ms (Altamirano et al. 2010). The open circle represents 4U 0614+091, which harbors a 2.41 ms pulsar (Strohmayer et al. 2008) and has a tentative orbital period of 51.3 min (Shahbaz et al. 2008; Hakala et al. 2011). 
Now we know that full propeller is energetically possible for UCXBs, does it actually happen? The observed neutron star spin periods, typically 5 ms, suggest that at least some spin down is plausible. On the other hand, the faint transient Xray binary SAX J1808.4–3658 proves that at least some accretion is still possible onto a very fast spinning (2.5 ms) neutron star, since it shows type I Xray bursts (in’t Zand et al. 1998). Spruit & Taam (1993); Rappaport et al. (2004) showed that some accretion likely persists assuming the inner disk radius comes quite close to the corotation radius. They also mentioned that this disk structure probably breaks down at very low mass transfer rate, and that full propelling probably happens then. Simulations by D’Angelo & Spruit (2011) showed that accelerated matter does not necessarily leave the system, because the inner disk radius may remain trapped at some short distance outside the corotation radius. This scenario, however, appears more likely for stronger magnetic fields (~10^{12} G) than the ones considered here, and for stars with a much shorter spindown timescale than millisecond pulsars. Therefore, a nonaccreting state seems more likely for UCXBs containing a recycled neutron star, and the propeller effect may work depending on whether the inner disk radius exceeds the corotation radius by a sufficiently large ratio.
Whether accretion at very low mass transfer rate is possible or not in a system with a stable accretion disk perhaps is not the right question to ask, since accretion most likely occurs irregularly. Figure 11 suggests that a helium disk (and also a carbonoxygen disk, which is not shown) remains divided in a neutral outer part and an ionized inner part, where the inner part at all times is too large for the magnetosphere to disrupt. This implies that the disk instability remains relevant and will come into play once the mass transfer rate becomes too low (Sect. 2.6.4).
4.3.2. Nonsteady mass transfer
Varying mass flow rates and densities near the accretor have a large impact on the interaction between the disk and the magnetic field. For instance, Spruit & Taam (1993) described cycles of accretion, during which matter piles up outside the Alfvén radius until the local density is such that the inner disk forces the Alfvén radius inwards, allowing accretion. When the disk density decreases, the Alfvén radius moves outwards again. In Fig. 11 it can be seen that the Alfvén radius at Eddingtonlimit mass transfer is almost equal to the latetime corotation radius, but taking into account the decreasing magnetic field, the Alfvén radius will always lie well inside the corotation radius at the maximum accretion rate, so a nearEddington outburst can always be accreted.
Fig. 13
Recurrence time for a 1.4 M_{⊙} neutron star accretor with a residual equatorial magnetic field of 10^{8.5} G (solid) and for a 10 M_{⊙} black hole accretor (dashed) with a helium disk described by Dunkel et al. (2006), for two values of the viscosity parameter α, of which 0.02 probably best describes the disk viscosity in quiescence. Disk instability becomes relevant to the right of the vertical line at 28 min. The circles indicate an age of 10 Gyr. The discontinuity in the slopes of the neutron star curves at an orbital period of 29 min is because the speedoflight cylinder defines the inner disk radius beyond this period (Fig. 11). 
In Sect. 3.5 we have discussed under which circumstances accretion could become irregular due to thermalviscous disk instability. When the mass transfer rate has become sufficiently low, the disk periodically collapses and builds, which leads to high mass transfer rates near the accretor every cycle. The recurrence time is the time it takes until the empty disk is refilled. Figure 13 shows the recurrence time for our UCXB model. The inner disk radius R_{in} is given by Eq. (35) in the case of a neutron star accretor and 3 Schwarzschild radii in the case of a black hole accretor. The outer radius is estimated by R_{out} ≈ 0.9d_{L1} based on Fig. 5. The recurrence time is increasing because the mass transfer rate decreases faster than the disk mass. During quiescence the disk viscosity is low, of the order of α = 0.02 (Lasota 2001). Systems of age 10^{9 − 10} yr with a neutron star accretor have a recurrence time of hundreds of years. In the case of an equally old system with black hole accretor, this time is a bit shorter. The long recurrence time suggests that systems are only rarely observable as bright Xray sources, so that at any given time, only a small fraction of a population is visible.
5. Discussion and conclusions
We have studied the evolution of ultracompact Xray binaries, and in particular their characteristics at an old age, when they reach a very low mass ratio and low mass transfer rate. The aim is to investigate how likely the occurrence of a dynamical instability is, and how important the propeller effect and disk instability are.
If the feedback of angular momentum from disk to orbit were to be reduced, the lifetime of an UCXB would be dramatically shortened. Based on accretion disk simulations (Sect. 4.1), we conclude that an accretion disk in a system of any realistically achievable mass ratio is neither significantly affected by resonant orbits, nor does it expand outside the accretor Roche lobe. The latter is in agreement with Priedhorsky & Verbunt (1988). Most likely the disk is still capable of full feedback of angular momentum to the orbit so that the donor will never be disrupted. Consequently, unless a new mechanism is discovered that limits or stops feedback, or pulsar radiation evaporates the donor by heating and stripping off the outer layers (Kluzniak et al. 1988; Ruderman et al. 1989; Shaham 1992), or UCXBs become temporarily detached at some stage, the vast majority of UCXBs today will have orbital periods around 70 − 80 min and experience low mass transfer rates of 10^{13}−10^{12} M_{⊙} yr^{1}, depending on the precise massradius relation of the donor. (Systems that started mass transfer more recently will have shorter orbital periods, but since the evolutionary timescale increases rapidly with increasing orbital period, most systems are expected to have a long orbital period at the present. The time between the zeroage main sequence and the onset of mass transfer is in general short, as will be shown in a forthcoming paper.)
However, all ~15 known UCXBs have much shorter orbital periods, below 60 min. Since these shortperiod systems evolve on a much shorter timescale than longperiod systems, there must exist a much larger population of apparently invisible, old, longperiod UCXBs in the Galaxy. A strong selection effect hides the large population of relatively common systems so that we see only a couple of much rarer systems during their shortlived high mass transfer phase. The bolometric luminosity corresponding to a constant mass transfer rate of ~10^{13} M_{⊙} yr^{1} is ~5 × 10^{32} erg s^{1}, within the observational limits of the best presentday Xray telescopes even at a distance of 8 kpc (the distance to the Galactic bulge) if these systems emit predominantly in Xrays. Then, the only consistent resolutions are that these longperiod systems are even fainter than expected, or visible only a small fraction of the time. These resolutions are supported by the propeller effect and the disk instability model, respectively.
The propeller effect probably strongly influences accretion in UCXBs, especially at very low mass transfer rates. In the event of a stable disk, rotational kinetic energy stored in the neutron star during the accreting stage is sufficient to expel all matter arriving at the magnetosphere from the system once the propeller effect starts to work. Therefore, systems with a stable accretion disk could have a significantly reduced Xray luminosity.
However, the disk instability model, which describes how disks periodically build and collapse, still applies at the very low mass transfer rates expected in UCXBs. During a collapse, matter is capable of penetrating the magnetosphere and reach the stellar surface, which means that observationally the system appears as a transient Xray source. Based on estimates of the disk structure, the recurrence time of old (10 Gyr) systems is of the order of a few 100 yr, depending on at what mass transfer rate the disk becomes unstable, among other things. Since black hole accretor systems do not suffer from the propeller effect, these would be expected to be better visible, but none have been discovered so far. Apparently they are much rarer than their neutron star accretor counterparts, or disk instability is the dominant process in making low mass ratio UCXBs nearly invisible.
A third mechanism that can inhibit the Xray visibility of a black hole UCXB is the inefficiently radiating disk (ADAF – advection dominated accretion flow) (Narayan & Yi 1994; Lasota 2008). In ADAFs, most of the thermal energy stored in the inner disk is accreted into the event horizon. For a neutron star UCXB this does not apply as the energy would still be released upon accretion.
Lastly, UCXBs with a low mass transfer rate may also become fainter in Xrays if they emit a significant fraction of their energy as UV or even optical radiation. The Dunkel et al. (2006) central disk temperature for helium composition is given by T_{c} ∝ Ṁ^{3/10} or Ṁ^{2/5} depending on the choice of opacity, while a thick disk has T ∝ Ṁ^{1/4}, so the disk spectrum shifts to longer wavelengths as the mass transfer rate decreases. The hot inner part of the disk, origin of most radiation, still emits mainly Xrays, with the result that the disk as a whole emits more than half of its energy as Xray radiation, even at the lowest mass transfer rates we find. The inner disk, however, is affected by the magnetic field in the case of neutron star accretors (in particular at low mass transfer rate), and can become an ADAF in the case of black hole accretors. In these cases the outer disk spectrum becomes more important and the system may emit most radiation in the UV, which is difficult to observe due to absorption (Campana & Stella 2000).
The consequences of the evolutionary behavior of UCXBs for the population in the Galactic bulge will be the topic of a forthcoming paper.
A stability criterion expressed in values, i.e. logarithmic derivatives, is to be preferred over linear derivatives of radius to donor mass because the expressions for and , Eqs. (13) and (26) respectively, are simpler than the expression for as the factors and introduce several extra variables.
Yungelson et al. (2002) found a value of , for a neutron star accretor, due to a different adopted isotropic reemission prescription combined with the relatively weak sensitivity to donor mass of both mass transfer rate and isotropic reemission rate.
The dynamical instability for high mass white dwarf donors on the other hand is a result of the steep massradius relation of the donor, and also because the orbit expands less easily when the mass ratio is closer to (Fig. 2).
Note that the equatorial Roche lobe radius is significantly larger than the average radius over the whole volume (Eq. (12)), similar to an equatorial bulge.
Around an intermediate mass or supermassive black hole, an accretion disk can only exist if M_{a} ≲ 10^{5} M_{⊙} in the case of a 0.6 M_{⊙} white dwarf donor and M_{a} ≲ 10^{6} M_{⊙} in the case of a 0.1 M_{⊙} white dwarf donor, since otherwise the innermost stable circular orbit exceeds the semimajor axis.
Acknowledgments
We thank E. Körding for helpful discussions on accretion disks. LMvH is supported by the Netherlands Organisation for Scientific Research (NWO). G.N. and R.V. are supported by NWO Vidi grant 016.093.305.
References
 Alpar, M. A., Cheng, A. F., Ruderman, M. A., & Shaham, J. 1982, Nature, 300, 728 [NASA ADS] [CrossRef] [Google Scholar]
 Altamirano, D., Patruno, A., Heinke, C. O., et al. 2010, ApJ, 712, L58 [NASA ADS] [CrossRef] [Google Scholar]
 Aly, J. J., & Kuijpers, J. 1990, A&A, 227, 473 [NASA ADS] [Google Scholar]
 Begelman, M. C. 1979, MNRAS, 187, 237 [NASA ADS] [Google Scholar]
 Bhattacharya, D., & van den Heuvel, E. P. J. 1991, Phys. Rep., 203, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Bhattacharyya, S., Strohmayer, T. E., Markwardt, C. B., & Swank, J. H. 2006, ApJ, 639, L31 [NASA ADS] [CrossRef] [Google Scholar]
 Bildsten, L. 2002, ApJ, 577, L27 [NASA ADS] [CrossRef] [Google Scholar]
 Bonsema, P. F. J., & van den Heuvel, E. P. J. 1985, A&A, 146, L3 [NASA ADS] [Google Scholar]
 Campana, S., & Stella, L. 2000, ApJ, 541, 849 [NASA ADS] [CrossRef] [Google Scholar]
 Chakrabarty, D., & Morgan, E. H. 1998, Nature, 394, 346 [Google Scholar]
 D’Angelo, C. R., & Spruit, H. C. 2011, MNRAS, 416, 893 [NASA ADS] [CrossRef] [Google Scholar]
 Davidson, K., & Ostriker, J. P. 1973, ApJ, 179, 585 [Google Scholar]
 Deloye, C. J., & Bildsten, L. 2003, ApJ, 598, 1217 [NASA ADS] [CrossRef] [Google Scholar]
 Dunkel, J., Chluba, J., & Sunyaev, R. A. 2006, Astron. Lett., 32, 257 [NASA ADS] [CrossRef] [Google Scholar]
 Eggleton, P. P. 1983, ApJ, 268, 368 [NASA ADS] [CrossRef] [Google Scholar]
 Elsner, R. F., & Lamb, F. K. 1977, ApJ, 215, 897 [NASA ADS] [CrossRef] [Google Scholar]
 Frank, J., King, A., & Raine, D. J. 2002, Accretion Power in Astrophysics: 3rd ed. [Google Scholar]
 Franklin, F. A., & Colombo, G. 1970, Icarus, 12, 338 [NASA ADS] [CrossRef] [Google Scholar]
 Galloway, D. K., Chakrabarty, D., Morgan, E. H., & Remillard, R. A. 2002, ApJ, 576, L137 [NASA ADS] [CrossRef] [Google Scholar]
 Hakala, P. J., Charles, P. A., & Muhli, P. 2011, MNRAS, 991 [Google Scholar]
 Homer, L., Charles, P. A., & O’Donoghue, D. 1998, MNRAS, 298, 497 [NASA ADS] [CrossRef] [Google Scholar]
 Hut, P., & Paczyński, B. 1984, ApJ, 284, 675 [NASA ADS] [CrossRef] [Google Scholar]
 Illarionov, A. F., & Sunyaev, R. A. 1975, A&A, 39, 185 [NASA ADS] [Google Scholar]
 in’t Zand, J. J. M., Heise, J., Muller, J. M., et al. 1998, A&A, 331, L25 [NASA ADS] [Google Scholar]
 in’t Zand, J. J. M., Cumming, A., van der Sluys, M. V., Verbunt, F., & Pols, O. R. 2005, A&A, 441, 675 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 in’t Zand, J. J. M., Jonker, P. G., & Markwardt, C. B. 2007, A&A, 465, 953 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ivanova, N., Rasio, F. A., Lombardi, Jr., J. C., Dooley, K. L., & Proulx, Z. F. 2005, ApJ, 621, L109 [NASA ADS] [CrossRef] [Google Scholar]
 Jeffrey, L. C. 1986, Nature, 319, 384 [NASA ADS] [CrossRef] [Google Scholar]
 Jonker, P. G., Bassa, C. G., Nelemans, G., et al. 2011, ApJS, 194, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Kaaret, P., Morgan, E. H., Vanderspek, R., & Tomsick, J. A. 2006, ApJ, 638, 963 [Google Scholar]
 King, A. R., & Begelman, M. C. 1999, ApJ, 519, L169 [NASA ADS] [CrossRef] [Google Scholar]
 Kluzniak, W., Ruderman, M., Shaham, J., & Tavani, M. 1988, Nature, 334, 225 [NASA ADS] [CrossRef] [Google Scholar]
 Krimm, H. A., Markwardt, C. B., Deloye, C. J., et al. 2007, ApJ, 668, L147 [NASA ADS] [CrossRef] [Google Scholar]
 Lamb, F. K., Pethick, C. J., & Pines, D. 1973, ApJ, 184, 271 [NASA ADS] [CrossRef] [Google Scholar]
 Landau, L. D., & Lifshitz, E. M. 1975, The classical theory of fields (Oxford: Pergamon Press) [Google Scholar]
 Lasota, J. 2001, New Astron. Rev., 45, 449 [NASA ADS] [CrossRef] [Google Scholar]
 Lasota, J. 2007, Comptes Rendus Physique, 8, 45 [NASA ADS] [CrossRef] [Google Scholar]
 Lasota, J. 2008, New Astron. Rev., 51, 752 [NASA ADS] [CrossRef] [Google Scholar]
 Lattimer, J. M., & Schutz, B. F. 2005, ApJ, 629, 979 [NASA ADS] [CrossRef] [Google Scholar]
 Li, X.D., & Wang, Z.R. 1999, ApJ, 513, 845 [Google Scholar]
 Lin, D. N. C., & Papaloizou, J. 1979, MNRAS, 186, 799 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, Q. Z., van Paradijs, J., & van den Heuvel, E. P. J. 2007, A&A, 469, 807 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lovelace, R. V. E., Romanova, M. M., & BisnovatyiKogan, G. S. 1999, ApJ, 514, 368 [NASA ADS] [CrossRef] [Google Scholar]
 Ma, B., & Li, X.D. 2009, ApJ, 698, 1907 [NASA ADS] [CrossRef] [Google Scholar]
 Markwardt, C. B., Juda, M., & Swank, J. H. 2003a, ATel, 127, 1 [NASA ADS] [Google Scholar]
 Markwardt, C. B., Smith, E., & Swank, J. H. 2003b, ATel, 122, 1 [Google Scholar]
 Markwardt, C. B., Swank, J. H., Strohmayer, T. E., in’t Zand, J. J. M., & Marshall, F. E. 2002, ApJ, 575, L21 [Google Scholar]
 Marsh, T. R., Nelemans, G., & Steeghs, D. 2004, MNRAS, 350, 113 [NASA ADS] [CrossRef] [Google Scholar]
 Menou, K., Perna, R., & Hernquist, L. 2002, ApJ, 564, L81 [NASA ADS] [CrossRef] [Google Scholar]
 Morris, S. L., Liebert, J., Stocke, J. T., et al. 1990, ApJ, 365, 686 [NASA ADS] [CrossRef] [Google Scholar]
 Narayan, R., & Yi, I. 1994, ApJ, 428, L13 [NASA ADS] [CrossRef] [Google Scholar]
 Nelemans, G. 2009, Classical and Quantum Gravity, 26, 094030 [Google Scholar]
 Nelemans, G., & Jonker, P. G. 2010, New Astron. Rev., 54, 87 [NASA ADS] [CrossRef] [Google Scholar]
 Nelemans, G., Jonker, P. G., Marsh, T. R., & van der Klis, M. 2004, MNRAS, 348, L7 [NASA ADS] [CrossRef] [Google Scholar]
 Nelemans, G., Jonker, P. G., & Steeghs, D. 2006, MNRAS, 370, 255 [NASA ADS] [CrossRef] [Google Scholar]
 Nelson, L. A., & Rappaport, S. 2003, ApJ, 598, 431 [NASA ADS] [CrossRef] [Google Scholar]
 Nelson, L. A., Rappaport, S. A., & Joss, P. C. 1986, ApJ, 304, 231 [NASA ADS] [CrossRef] [Google Scholar]
 Osaki, Y. 1974, PASJ, 26, 429 [NASA ADS] [Google Scholar]
 Osaki, Y. 1989, PASJ, 41, 1005 [NASA ADS] [Google Scholar]
 Paczyński, B. 1971, ARA&A, 9, 183 [NASA ADS] [CrossRef] [Google Scholar]
 Paczyński, B. 1977, ApJ, 216, 822 [NASA ADS] [CrossRef] [Google Scholar]
 Papaloizou, J., & Pringle, J. E. 1977, MNRAS, 181, 441 [NASA ADS] [CrossRef] [Google Scholar]
 Podsiadlowski, P., Rappaport, S., & Pfahl, E. D. 2002, ApJ, 565, 1107 [NASA ADS] [CrossRef] [Google Scholar]
 Pols, O. R., & Marinus, M. 1994, A&A, 288, 475 [NASA ADS] [Google Scholar]
 Priedhorsky, W. C., & Verbunt, F. 1988, ApJ, 333, 895 [NASA ADS] [CrossRef] [Google Scholar]
 Rappaport, S., Ma, C. P., Joss, P. C., & Nelson, L. A. 1987, ApJ, 322, 842 [NASA ADS] [CrossRef] [Google Scholar]
 Rappaport, S. A., Fregeau, J. M., & Spruit, H. 2004, ApJ, 606, 436 [NASA ADS] [CrossRef] [Google Scholar]
 Ruderman, M., Shaham, J., & Tavani, M. 1989, ApJ, 336, 507 [NASA ADS] [CrossRef] [Google Scholar]
 Ruderman, M. A., & Shaham, J. 1983, Nature, 304, 425 [NASA ADS] [CrossRef] [Google Scholar]
 Ruderman, M. A., & Shaham, J. 1985, ApJ, 289, 244 [NASA ADS] [CrossRef] [Google Scholar]
 Savonije, G. J., de Kool, M., & van den Heuvel, E. P. J. 1986, A&A, 155, 51 [NASA ADS] [Google Scholar]
 Schreiber, M. R., & Lasota, J.P. 2007, A&A, 473, 897 [Google Scholar]
 Schulz, N. S., Chakrabarty, D., Marshall, H. L., et al. 2001, ApJ, 563, 941 [NASA ADS] [CrossRef] [Google Scholar]
 Shaham, J. 1992, in XRay Binaries and the Formation of Binary and Millisecond Radio Pulsars, 375 [Google Scholar]
 Shahbaz, T., Watson, C. A., Zurita, C., Villaver, E., & HernandezPeralta, H. 2008, PASP, 120, 848 [NASA ADS] [CrossRef] [Google Scholar]
 Shibazaki, N., Murakami, T., Shaham, J., & Nomoto, K. 1989, Nature, 342, 656 [NASA ADS] [CrossRef] [Google Scholar]
 Simpson, J. C., & Wood, M. A. 1998, ApJ, 506, 360 [NASA ADS] [CrossRef] [Google Scholar]
 Soberman, G. E., Phinney, E. S., & van den Heuvel, E. P. J. 1997, A&A, 327, 620 [NASA ADS] [Google Scholar]
 Spruit, H. C., & Taam, R. E. 1993, ApJ, 402, 593 [Google Scholar]
 Strohmayer, T. E., Markwardt, C. B., & Kuulkers, E. 2008, ApJ, 672, L37 [NASA ADS] [CrossRef] [Google Scholar]
 Tauris, T. M., & Savonije, G. J. 1999, A&A, 350, 928 [NASA ADS] [Google Scholar]
 Tsugawa, M., & Osaki, Y. 1997, PASJ, 49, 75 [NASA ADS] [Google Scholar]
 Tutukov, A. V., & Yungelson, L. R. 1993, Astron. Rep., 37, 411 [NASA ADS] [Google Scholar]
 Ustyugova, G. V., Koldoba, A. V., Romanova, M. M., & Lovelace, R. V. E. 2006, ApJ, 646, 304 [NASA ADS] [CrossRef] [Google Scholar]
 van den Heuvel, E. P. J. 1984, JA&A, 5, 209 [NASA ADS] [CrossRef] [Google Scholar]
 van den Heuvel, E. P. J., & Bonsema, P. F. J. 1984, A&A, 139, L16 [NASA ADS] [Google Scholar]
 van der Sluys, M. V., Verbunt, F., & Pols, O. R. 2005, A&A, 431, 647 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Verbunt, F. 1987, ApJ, 312, L23 [NASA ADS] [CrossRef] [Google Scholar]
 Verbunt, F., & Rappaport, S. 1988, ApJ, 332, 193 [NASA ADS] [CrossRef] [Google Scholar]
 Voss, R., & Gilfanov, M. 2007, MNRAS, 380, 1685 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, J., Zhang, C. M., Zhao, Y. H., et al. 2011, A&A, 526, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Warner, B. 2003, Cataclysmic Variable Stars, ed. B. Warner [Google Scholar]
 Webbink, R. F. 1985, Stellar evolution and binaries, ed. Pringle, J. E., & Wade, R. A., 39 [Google Scholar]
 Wijnands, R. 2010, Highlights Astron., 15, 121 [Google Scholar]
 Wood, M. A., Thomas, D. M., & Simpson, J. C. 2009, MNRAS, 398, 2110 [NASA ADS] [CrossRef] [Google Scholar]
 Yungelson, L. R. 2008, Astron. Lett., 34, 620 [NASA ADS] [CrossRef] [Google Scholar]
 Yungelson, L. R., Nelemans, G., & van den Heuvel, E. P. J. 2002, A&A, 388, 546 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Yungelson, L. R., Lasota, J., Nelemans, G., et al. 2006, A&A, 454, 559 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zhang, C. M., & Kojima, Y. 2006, MNRAS, 366, 137 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Fitted tracks
Sometimes it is useful to have simple analytic estimates of the evolution of the quantities involved. For this purpose we provide powerlaw approximations of the tracks with full feedback of angular momentum, optimized for system age (time since the onset of mass transfer) t ∈ [1 Myr,20 Gyr] , (A.1)with t in yr, the mass transfer rate Ṁ_{d} in M_{⊙} yr^{1}, the orbital period P_{orb} in min and the donor mass M_{d} in M_{⊙}. The corresponding parameters f_{i} and e_{i} for accretor masses M_{a} = 1.4 M_{⊙} and 10 M_{⊙} are listed in Table A.1. The exponent in the Ṁ_{d}(t) relation (e_{1}) must be roughly 1 lower than in the M_{d}(t) relation (e_{3}), but not necessarily precisely 1 because both tracks have been fitted separately.
For an initially 1.4 M_{⊙} neutron star accretor, these relations fit the modeled helium white dwarf tracks (shown in Figs. 1, 8 and 9) with an average error of 4.3%, in the tracks t = 1 Myr corresponds to P_{orb} = 9.0 min, Ṁ_{d} = 2.1 × 10^{8} M_{⊙} yr^{1} (lower than the Eddington limit; mass transfer is conservative for the whole range) and M_{d} = 0.074 M_{⊙}.
Similarly for an initially 10 M_{⊙} black hole accretor, where the average error is 5.3%. For a 10 M_{⊙} black hole, t = 1 Myr corresponds to P_{orb} = 12.3 min, Ṁ_{d} = 1.4 × 10^{8} M_{⊙} yr^{1} and M_{d} = 0.050 M_{⊙}.
From the fitted parameters as well as in Fig. 1, it can be seen that all slopes are nearly independent of accretor mass. Also, the M_{d}(P_{orb}) relations are rather similar for both accretor masses, as expected.
Appendix B: Approximate analytic solution
By making a few mild assumptions it is possible to solve almost the entire evolution (assuming full feedback of angular momentum to the orbit) of a binary containing a zerotemperature helium white dwarf donor analytically, except for the early stages when the donor mass is high. These derivations only serve to gain physical insight into the evolution of UCXBs – throughout the paper numerical solutions of more precise equations have been used.
First we approximate the lowmass white dwarf massradius relation from Eq. (25) by the nonrelativistic relation (B.1)which is accurate around M_{d} = 0.01 M_{⊙} (and within 20% of the radius given by Eq. (25) for 4 × 10^{3} < M_{d}/M_{⊙} < 1), but still represents the largest simplification made here. By substituting the white dwarf radius from Eq. (B.1) into the Rochelobe radius in the wellknown periodmean density relation we find the orbital period of a system with a white dwarf donor, (B.2)with the constant (B.3)The orbital period is inversely proportional to donor mass (and proportional to donor volume). The mass transfer rate is given by Eq. (6), in which ζ_{d} = −1/3 follows from Eq. (B.1) and ζ_{L} = −5/3 is the low mass ratio limit of Eq. (13) for η = 0. Furthermore, we can use the gravitational wave equation and Kepler’s third law to replace a by P_{orb} to obtain (B.4)Lastly, M_{d} ≪ M_{a} implies M_{tot} ≈ M_{a}, and we substitute Eq. (B.2) for P_{orb} to express the mass transfer rate in only the component masses (B.5)with the constant (B.6)Equation (B.5) shows that the mass transfer rate depends strongly on donor mass. 8/3 out of the 14/3 power is due to massive donors having a shorter semimajor axis when they fill their Roche lobe, which results in stronger gravitational wave radiation. The direct effect of the high donor mass is slightly smaller with 2 out of 14/3. Alternatively, the mass transfer rate can be expressed as a function of accretor mass and orbital period (B.7)By integrating the mass transfer rate given by Eq. (B.5), we obtain the system age t as a function of component masses. (B.8)where integrating from ∞ is justified by the negligibly short timescale associated with a high donor mass, as evidenced by the − 11/3 power of the donor mass. The donor mass as a function of time follows, (B.9)and by substituting Eq. (B.2) (B.10)The evolutionary timescale is , even steeper than the gravitational wave timescale in a detached binary () due to the effect of the donor mass. Finally, differentiating Eq. (B.9) to time yields (B.11)As mentioned in Sect. 4.2.2 it shows that when comparing for systems with the same age, a less massive accretor corresponds to a higher donor mass and mass transfer rate. The mass transfer timescale of systems of fixed age t is given by . This value is 2/3 times the gravitational wave timescale, and is independent of accretor mass in this approximation. As we saw in Sect. 4.2.2, however, this timescale is very slightly (~7% for a factor 7 in accretor mass) shorter for highmass accretor systems when we use the more precise massradius relation.
Above relations are valid for M_{d} ≲ 0.2 M_{⊙}. Table A.1 lists numerical values for the parameters in several equations given above.
B.1. Luminosity
By assuming an accretor radius and taking the timeaveraged bolometric luminosity of an UCXB as (B.12)where Ṁ_{d} is still defined as negative, and by using Eq. (B.7), we can express the following quantities in the observables L and P_{orb}: (B.13)and (B.14)By combining Eqs. (B.10) and (B.13) the system age follows (B.15)
Appendix C: First Lagrangian point location approximations
The distance between the first Lagrangian point (L1) and binary component 2 is given by (C.1)where the relative errors in d(L1,2) and d(L1,1) are less than 0.015% for 7 × 10^{4} < q = M_{2}/M_{1} < 1.5 × 10^{3} and less than 0.1% for 10^{6} < q < 10^{6}.
The next two expressions are only valid for q ≤ 1. The distance between L1 and the least massive binary component (2), valid within 0.02% (again, this is the relative error in the distance to either component) for 7 × 10^{5} < q = M_{2}/M_{1} ≤ 1, is given by (C.2)A simpler but less accurate (<0.5%) expression, inspired by the form of the Eggleton Roche lobe Eq. (12), is^{11}(C.3)where 0 < q = M_{2}/M_{1} ≤ 1.
Warner (2003, lists three approximations of the L1 location, but those are valid over a much smaller range in mass ratio and have much larger errors than our equations.
Appendix D: Roche lobe potential approximations
Two simpler but less accurate versions of Eq. (10) are given here. For a binary with 0 < q = M_{2}/M_{1} ≤ 1, the potential at the Roche lobe surface Φ_{L} is approximated to within 0.4% by (D.1)and to within 1% by (D.2)
All Tables
All Figures
Fig. 1
Tracks for a helium white dwarf donor with an initially 1.4 M_{⊙}, 12 km radius neutron star (solid) and 10 M_{⊙} black hole (dashed) accretor. The tracks start at the respective dynamical instability limits of 0.83 M_{⊙} and 1.04 M_{⊙}, since systems can only survive if they start with a lower donor mass (Sect. 2.3). The isotropic reemission limit is ignored in the first part of the track. The numbers correspond to the two circles on their righthand side and represent the logarithm of the system age (yr), which is the time since the onset of mass transfer. The relative location of both tracks will be discussed in Sect. 4.2.2. 

In the text 
Fig. 2
Massradius relation exponent ζ_{L} (Sect. 3.2) for a 1.4 M_{⊙} accretor (dotted) and a 10 M_{⊙} accretor (dashdotted) and the zerotemperature helium white dwarf ζ_{d} (dashed, Sect. 3.3). Nearly all matter is assumed to be lost via the isotropic reemission mechanism (Sect. 2.4). 

In the text 
Fig. 3
Stability regimes in terms of zerotemperature helium white dwarf donor mass versus accretor mass. The upper solid curve is the upper limit for dynamical stability in the presence of an accretion disk that feeds back angular momentum, where all transferred matter is lost via isotropic reemission. The dashed curve is the upper limit for isotropic reemission (i.e. if the system lies above this curve, the accretor will be engulfed by transferred matter). The lower solid curve encloses the parameter space where systems are dynamically stable in the absence of angular momentum feedback. The dotted curve indicates a very rough estimate of the upper donor mass for which there is enough time to form an accretion disk and suppress the nofeedback instability, based on Verbunt & Rappaport (1988). 

In the text 
Fig. 4
Accretion efficiency ϵ (Eq. (17)) for a helium white dwarf donor and an initially 1.4 M_{⊙} neutron star (solid) and a 10 M_{⊙} black hole (dashed) accretor, in the case of full feedback of angular momentum. Both curves start at their respective dynamical stability limit (Sect. 2.3). 

In the text 
Fig. 5
Important radii within the accretion disk; the distance between the accretor and the L1 point d_{L1} (solid), the average equatorial Roche lobe radius (i.e. in the orbital plane) (dashed), the circularization radius R_{circ} (dashdotted) and R_{circ} corrected for donor influence during freefall, R_{h} (dotted, Eq. (19)). Furthermore, the pluses (+) indicate the outer tidally stable radius from Papaloizou & Pringle (1977). The stars ( ⋆ ) indicate the largest orbit that does not intersect others, according to Paczyński (1977) (which lies near the 3:1 resonance). The two open circles are outer disk estimates by Priedhorsky & Verbunt (1988), and the three filled triangles are outer disk radii from our SPH simulations. 

In the text 
Fig. 6
Numerical smoothed particle hydrodynamics simulation of an accretion disk in dynamical equilibrium for mass ratio 10^{3} after 400 orbits. The solid curve represents the equatorial Roche lobes. The accretor has been marked by a small open circle. The dashed circles indicate the 3:1 (inner) and 2:1 (outer) resonances with the orbital period. 

In the text 
Fig. 7
Massradius relation exponent ζ_{L} (Sect. 3.2) for various modes of mass transfer from a zerotemperature helium white dwarf to a neutron star of initially 1.4 M_{⊙} with a 12 km radius. The lower dotted curve shows the case of accretion efficiency ϵ = 0 and full feedback (η = 0) (almost all transferred matter lost via isotropic reemission, hardly any accretion). The middle dotted curve is for ϵ = 1 and η = 0 (conservative mass transfer). The upper dotted curve is for no feedback (η = 1). The zerotemperature white dwarf exponent ζ_{d} (Sect. 3.3) is represented by the dashed curve. The solid curves starting at q = 0.59 (M_{d} = 0.83 M_{⊙}, Sect. 2.3) indicate the actual value of ζ_{L} (for ϵ shown in Fig. 4) for both the full feedback and nofeedback case. Until q = 0.01 there is full feedback in either case. Below q = 0.01 the full feedback case is shown by the lower solid branch, where the open circle denotes an age of 10 Gyr (since the onset of mass transfer). The nofeedback case has a jump in ζ_{L} at q = 0.01, after which it follows the upper dotted curve until the intersection with the dashed curve (filled circle) where donor disruption occurs at a system age of 300 Myr. 

In the text 
Fig. 8
Mass transfer rates for a zerotemperature helium white dwarf donor with an initially 1.4 M_{⊙}, 12 km radius neutron star accretor (lower solid curve) and a 10 M_{⊙} black hole accretor (upper solid curve), starting at the dynamical instability limit (i.e. ignoring possible merging), for full feedback. The dashed curves branching away from these represent the evolution assuming no feedback below q = 0.01 (Sect. 4.1). The numbers correspond to the two circles on their lefthand side and represent the logarithm of the system age (yr), which is the time since the onset of mass transfer. Indicated are the isotropic reemission regime (dashed, Eq. (11)), the Eddington limit (dotted) and the accretion disk instability limit by in’t Zand et al. (2007) (dashdotted, Sect. 3.5) Note that the dashed neutron star curve at its local maximum “touching” the upper dashdotted curve is a coincidence. 

In the text 
Fig. 9
Orbital period against time for a zerotemperature helium white dwarf donor with an initially 1.4 M_{⊙}, 12 km radius neutron star (lower solid) and 10 M_{⊙} black hole (upper solid) accretor. 

In the text 
Fig. 10
Stability regimes in terms of zerotemperature helium white dwarf donor mass versus accretor mass. The upper solid curve is the upper limit for dynamical stability in the presence of an accretion disk that feeds back angular momentum, where almost all transferred matter is lost via isotropic reemission. The upper dashed curve is the upper stability limit for isotropic reemission and the lower dashed curve indicates the proposed q = 0.01 lower mass ratio limit for feedback (Sect. 4.1). The lower solid curve encloses the parameter space where systems are dynamically stable in the absence of angular momentum feedback; all matter enters and stays in the disk. The dashdotted lines indicate when the propeller effect comes into play according to Eq. (31) for several neutron star spin periods and an equatorial magnetic field of 10^{8.5} G. The dotted isochrones indicate the donor masses for system ages (since the onset of mass transfer) of 0.1,1,10 Gyr based on full feedback. The dynamical stability and isotropic reemission rates are independent of accretor radius and type, hence these curves are continuous. Moreover, the dynamical stability rates are independent of accretor type and therefore also hold for white dwarf accretors. 

In the text 
Fig. 11
Radii in UCXB accretion disk versus orbital period for a 1.4 M_{⊙} neutron star accretor with a residual equatorial magnetic field of 10^{8.5} G. The shaded layer at the bottom represents the neutron star surface, at a radius of 12 km. The curve annotated L1 shows the distance of the first Lagrangian point to the neutron star, which is an upper limit on the outer disk radius. The circularization radius by Verbunt & Rappaport (1988) is discussed in Sect. 4.1. The two upper dashed curves indicate the radii where the temperature of a helium disk (Dunkel et al. 2006) equals the He I ionization temperature, for two values of the viscosity parameter α, as discussed in Sect. 2.6.4. The solid corotation radius curve applies to the full propeller case, the dashed one to the hypothetical full accretion case. The speedoflight cylinder radius and Alfvén radius are also shown. 

In the text 
Fig. 12
Spin period versus orbital period for a 1.4 M_{⊙} neutron star accretor with a residual equatorial magnetic field of 10^{8.5} G, assuming all transferred matter is propelled out of the system once R_{μ} > 2^{1/3}R_{co} (solid, Sect. 3.4.2), or hypothetically keeps being accreted (dashed), or hypothetically is accreted at the neutron star surface all the time (dashdotted). Observed millisecond pulsarUCXB parameters are represented by filled circles; XTE J1807–294 at an orbital period of 40.07 min (Markwardt et al. 2003a) and a spin period of 5.25 ms (Markwardt et al. 2003b), XTE J1751–305 at 42.42 min, 2.30 ms (Markwardt et al. 2002), XTE J0929–314 at 43.58 min, 5.40 ms (Galloway et al. 2002), SWIFT J1756.9–2508 at 54.70 min, 5.49 ms (Krimm et al. 2007) and NGC 6440 X2 at 57.3 min, 4.86 ms (Altamirano et al. 2010). The open circle represents 4U 0614+091, which harbors a 2.41 ms pulsar (Strohmayer et al. 2008) and has a tentative orbital period of 51.3 min (Shahbaz et al. 2008; Hakala et al. 2011). 

In the text 
Fig. 13
Recurrence time for a 1.4 M_{⊙} neutron star accretor with a residual equatorial magnetic field of 10^{8.5} G (solid) and for a 10 M_{⊙} black hole accretor (dashed) with a helium disk described by Dunkel et al. (2006), for two values of the viscosity parameter α, of which 0.02 probably best describes the disk viscosity in quiescence. Disk instability becomes relevant to the right of the vertical line at 28 min. The circles indicate an age of 10 Gyr. The discontinuity in the slopes of the neutron star curves at an orbital period of 29 min is because the speedoflight cylinder defines the inner disk radius beyond this period (Fig. 11). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.