arxiv: v1 [astro-ph] 11 Oct 2008

Transcription

1 The bifurcation periods in low-mass X-ray binaries: the effect of magnetic braking and mass loss Bo Ma and Xiang-Dong Li arxiv: v1 [astro-ph] 11 Oct 2008 Department of Astronomy, Nanjing University, Nanjing , China ABSTRACT The bifurcation period in low-mass X-ray binaries is the initial orbital period which separates the formation of converging systems (which evolve with decreasing orbital periods until the donor becomes degenerate) from the diverging systems (which evolve with increasing orbital periods until the donor star loses its envelope and a wide detached binary is formed). We calculate systematically the bifurcation periods of binary systems with a 1.4M neutron star and a 0.5 2M donor star, taking into account different kinds of magnetic braking and mass loss mechanisms. Our results show that the saturated magnetic braking can considerably decrease the values of bifurcation period compared to the traditional magnetic braking, while the influence of mass loss mechanisms on bifurcation periods is quite weak. We also develop a semi-analytical method to compute the bifurcation period, the result of which agrees well with the numerical method in the leading order. Subject headings: binaries: close stars: evolution X-rays: binaries 1. Introduction One interesting and important topic in the secular evolution of low-mass X-ray binaries (LMXBs) is the so-called bifurcation period P bif, the initial binary orbital period which separates the formation of converging systems (which evolve with decreasing orbital periods until the donor becomes degenerate) from the diverging systems (which evolve with increasing orbital periods until the donor star loses its envelope and a wide detached binary is formed) (Tutukov et al. 1985). The first systematic investigations on the bifurcation period were done by Pylyser & Savonije (1988, 1989). Neglecting mass loss from the binary system and assuming angular momentum loss due to magnetic braking (MB; Verbunt & Zwaan

2 2 1981) and gravitational radiation (GR; Landau & Lifshitz 1975), these authors found that the bifurcation period is in the range P bif day for LMXBs, and strongly depends on magnetic braking efficiency. Ergma et al. (1998) included mass loss from the binary system and re-calculated the bifurcation period for two mass configurations (M 1 /M, M 2 /M ) = (1.4, 1) and (1.4, 1.5) and two chemical compositions (Z = 0.003, 0.03). They pointed out that the mass loss from the binary system also plays an impotent role besides magnetic braking in determining the value of P bif, while the chemical composition could only cause small change in P bif. Their bifurcation periods are P bif day under conservative mass transfer, and times larger if moderate non-conservative mass transfer is assumed. Podsiadlowski et al. (2002) found a bifurcation period around 18 hr for a 1.4M NS and a 1M companion star, where they defined the bifurcation period as the orbital period when the Roche lobe overflow just began, instead of the initial orbital period. van der Sluys et al. (2005a,b) also investigated the bifurcation period in LMXBs focusing the formation of ultra-compact X-ray binaries (UCXBs), and specified the bifurcation period as the longest initial period that leads to UCXBs within a Hubble time (13.7 Gyr). UCXBs are bright X-ray sources with very short orbital periods (P 1 h). The donor has to be a compact source like a white dwarf or a compact core of an evolved giant star to fit in the small Roche lobe size. Such sources may be formed through dynamical processes including stellar collisions and common envelop evolution (Clark et al. 1975; Rasio et al. 2000; Lombardi et al. 2006). An alternative scenario for the formation of such sources is through stable mass transfer in X-ray binaries with a low- or intermediate-mass donor star, which may explain the negative derivative of the 11-min source in NGC 6624 (van der Klis et al. 1993; Chou & Grindlay 2001). It has been found that systems with initial orbital period just below the bifurcation period may form UCXBs (Nelson et al. 1986; Tutukov et al. 1987; Pylyser & Savonije 1988; Podsiadlowski et al. 2002; van der Sluys et al. 2005a). Podsiadlowski et al. (2002) showed that the closer the initial orbital period to the bifurcation period from below, the smaller the minimum orbital period will be achieved. So the value of bifurcation period is crucial to understanding the formation of UCXBs (van der Sluys et al. 2005b). In this paper we make a systematic investigation on the bifurcation period for binary systems containing an NS with a main-sequence (MS) companion of mass from 0.5M to 2M. This work was motivated by recent progress in studies on mass and angular momentum loss mechanisms in LMXB evolution. In previous works the MB law originally postulated by Verbunt & Zwaan (1981) and Rappaport et al. (1983) was usually adopted. However, this law predicts too fast spin-down of low-mass MS stars, contradicted with the observation of rapid rotators in young open clusters (Sills et al. 2000; Andronov et al. 2003). Obviously a modification of the MB law will have significant influence on the period evolu-

3 3 tion (van der Sluys et al. 2005b). Additionally, there is strong evidence that during LMXB evolution the mass transfer is highly non-conservative. Recent measurements of the masses of binary and millisecond pulsars indicate that a large fraction of the transferred mass may be lost from the systems rather accreted by the NS (Bassa et al. 2006; Steeghs & Jonker 2007, and references therein) 1. Theoretically, possible ways of mass loss have been suggested, including evaporation of the donor (Ruderman et al. 1989) or radio-ejection of the transferred material (Burderi et al. 2001, 2002; D Antona et al. 2006) due to the pulsar radiation/wind impinging on. In the latter case, the matter is lost from the system at the inner Lagrangian (L 1 ) point, carrying away angular momentum and altering the period evolution. This paper is organized as follows. 2 briefly describes the stellar evolution code, the binary models, and the physical assumptions, especially the MB laws and the mass loss mechanisms. Then we present the calculated results in 3. Our discussion and conclusions are given in Evolution code and binary model 2.1. The stellar evolution code We use an updated version of the stellar evolution code originally developed by Eggleton (1971, 1972, see also Han et al. 2004, Pols et al. 1995) to calculate the evolutions of binaries consisting of an NS (of mass M 1 ) and an MS secondary (of mass M 2 ). For the secondary star we assume a solar chemical composition (X = 0.70, Y = 0.28, and Z = 0.02), the ratio of mixing length to pressure scale height α = 2.0, and convective overshooting parameter to be The opacity table is from Rogers & Iglesias (1992), Alexander & Ferguson (1994) and Hubbard & Lampe (1969). The effective radius of the Roche lobe for the secondary is taken from Eggleton (1983), R L,2 = 0.49q 2/3 0.6q 2/3 + ln(1 + q 1/3 ) a (1) where q = M 2 /M 1 is the mass ratio, and a is the orbital separation. Mass transfer rate via Roche lobe overflow is evaluated as Ṁ2 = RMT max(0, (R 2 /R L,2 1) 3 )M yr 1, and we adopt RMT = 10 3 in the calculations. 1 The massive ( 1.9M ) NS discovered in the globular cluster M5 (Freire et al. 2008) may reflect a bimodal distribution of the initial masses of NSs (rather heavy accretion during the previous LMXB evolution), as already predicted by hydrodynamical core collapse simulations (Timmes et al. 1996).

4 Mass and Angular momentum loss mechanisms For LMXBs the timescale of tidal synchronization is much shorter than the characteristic evolutionary timescale of the binary, so we can assume that the spin of the secondary star and the binary orbital revolution are always synchronized. Assuming rigid body rotation of the secondary star and neglecting the spin angular momentum of the neutron star, the total angular momentum of the binary system can be expressed as J = I 2 ω + J orb = I 2 ω + G 2/3 M 1 M 2 (M 1 + M 2 ) 1/3 ω 1/3 (2) where I 2 is the moment of inertia of the secondary star, ω is the angular velocity of the binary. We consider three kinds of mechanisms of angular momentum loss. The first is the angular momentum loss due to gravitational radiation (Landau & Lifshitz 1975) dj GR dt = 32 G 7/2 M1M 2 2(M M 2 ) 1/2, (3) 5 c 5 a 7/2 where c is the light speed. This mechanism is important only in very short period binary systems. The second angular momentum loss mechanism is for non-conservative mass transfer. We assume that a fraction α of the transferred mass is accreted by the NS, and the remaining mass is ejected out of the binary as isotropic winds from the NS, carrying away the specific angular momentum of the NS, dj ML dt q = (1 α)ṁ2( 1 + q )2 a 2 ω. (4) In our numerical calculations we have set α = 0. Alternatively, if the NS is spun up to be a millisecond pulsar, its radiation pressure may be strong enough to halt the transferred matter at the L 1 point and quench the accretion. This radio ejection may cause almost all the matter from the secondary to be lost from the binary (Burderi et al. 2001, 2002). The corresponding rate of angular momentum loss is dj ML dt = Ṁ2a 2 L1 ω (5) where a L1 is the distance from the L 1 point to the center of mass of the binary system. The third angular momentum loss mechanism is MB. For a low-mass MS star with deep convection zone, stellar winds which are magnetically coupled with the star can decelerate the stellar spin efficiently, thus carrying away the orbital angular momentum because

5 5 of tidal synchronization. The widely used formula for such MB effect was postulated by Verbunt & Zwaan (1981) and Rappaport et al. (1983) as dj MB dt = M 2 R 4 2 ω3 dyn cm. (6) However, observations of rapid rotators in young open clusters suggest a modification of the MB law at high rotation rate (Sills et al. 2000), dj MB dt dj MB dt ( ) 0.5 ( ) 0.5 = Kω 3 R2 M2, ω ω crit, R M ( ) 0.5 ( ) 0.5 = Kωcritω 2 R2 M2, ω > ω crit, (7) R M where K = gcm 2 s (Andronov et al. 2003), ω crit is the critical angular velocity at which the angular momentum loss rate reaches a saturated state, given by (Krishnamurthi et al. 1997) ω crit (t) = ω crit τ t0, τ t, (8) where τ t0, and τ t are the global turnover timescales for the convective envelope of the Sun at its current age and of the secondary star at age t, respectively. They can be calculated by integrating the inverse local convective velocity over the surface convective envelope (Kim & Demarque 1996): R2 dr τ t = v, (9) R b where R b is the radial distance from the center of the star to the bottom of the surface convective envelope, and v is the local convective velocity from mixing-length theory (Böhm-Vitense 1958). Our calculation gives τ t0, 28.4 d, slightly larger than τ t0, 13.8 d in van der Sluys et al. (2005b), but consistent with the results of Kim & Demarque (1996) and Jung & Kim (2007). See Eggleton (2006, page 46) for the discussion of a possible reason for different values of τ t0, calculated. Following the suggestion of Podsiadlowski et al. (2002), we also add an ad hoc factor exp( 0.02/q conv + 1) if q conv < 0.02, in Eqs. (6) and (7), where q conv is the mass fraction of the surface convective envelop, to reduce the MB effect when the convective envelope becomes too small.

6 Binary models To examine the influence of mass and angular momentum loss mechanisms on the period evolution, we construct four models with various mass and angular momentum loss combinations: (1) model 1 - conservative mass transfer with traditional MB law (Eq. [6]); (2) model 2 - conservative mass transfer with saturated MB law (Eq. [7]); (3) model 3 - non-conservative mass transfer with mass loss from L 1 point (Eq. [5]) and saturated MB law (Eq. [7]); and (4) model 4 - non-conservative mass transfer with mass loss from the NS (Eq. [4]) and saturated MB law (Eq. [7]). In all the four models, the initial NS mass is set to be M 1,i = 1.4M, and the initial mass of the secondary M 2,i ranges from 0.5 to 2.0M. 3. Numeric Results 3.1. The bifurcation periods Throughout this paper we define the bifurcation period P bif as the initial binary orbital period P i with a zero-age main-sequence (ZAMS) companion star that separates converging from diverging systems. We use P f to denote the final orbital period after the mass transfer. Another definition of the bifurcation period used by Podsiadlowski et al. (2002) is the orbital period when the Roche lobe overflow just begins, which is expressed as P rlof in this paper. The results of the bifurcation periods for the four models described in 2.3 are summarized in Fig. 1 and Table 1. We also draw the minimum initial period P ZAMS that corresponds to a lobe-filling ZAMS donor star in Fig. 1. Several features are noted for the bifurcation periods in Fig. 1. First, the bifurcation periods for all the four models decrease with increasing initial secondary mass from 0.5M to 1.3M. Second, in models with saturated MB, there exists an upper limit of the initial secondary mass, beyond which no bifurcation period exists. This upper limit is in the range M for model 3, and M for models 2 and 4. Third, comparing the bifurcation periods of model 1 with those of models 2 4 indicates that the MB law plays the most important role in determining the values of the bifurcation periods: different MB laws can change the bifurcation periods by as much as 60%, compared to 14% (see Table 1) caused by different mass loss mechanisms. In Table 1 we also present P rlof following Podsiadlowski et al. (2002). In model 1, we get P rlof 18.3 hr for M 2,i = 1M, which is in close line with the result of Podsiadlowski et al. (2002) (17.7 hr with Z = 0.001, Y = 0.27), where the difference could be explained as the difference between the metallicities we used. When we use saturated MB, P rlof decreases to 11 hr.

7 7 According to the calculated orbital period evolutions, LMXBs can be classified into three categories: the diverging systems with P f P i, the converging systems with P f P i, and the parallel systems with P f P i. As an example, we present the calculated results for a 1.4M + 1.0M binary in model 3, to illustrate the three kinds of evolutionary sequences in Fig. 2. The corresponding bifurcation period is found to be 1.25 day, and the initial orbital periods are chosen to be P i = 1.20, 1.25, and 1.40 days, which represent the converging, parallel, and diverging systems respectively Effect of MB and mass loss Pylyser & Savonije (1988) have emphasized the effect of MB on the evolution of LMXBs. Comparing the results of models 1 and 2 presented in Table 1, we find that the bifurcation periods with traditional MB law are smaller (larger) than those with saturated MB law, when the initial secondary star mass M 2,i is less (larger) than 0.7M. Our results suggest that mass loss also influences the value of P bif, though in an less important way compared with MB. The bifurcation periods in non-conservative models 3 and 4 are lower than those in model 2, in which conservative mass transfer has been assumed. This result is consistent with van der Sluys et al. (2005b), but contradicted with Ergma et al. (1998). It is also interesting to see whether an UCXB can form with saturated MB. For an LMXB with an initial orbital period below the bifurcation period, mass transfer is mainly driven by the loss of angular momentum. The orbital period will decrease with the donor mass until a minimum period is reached. Paczynski & Sienkiewicz (1981) found a minimum period about 80 min without MB, while Podsiadlowski et al. (2002) showed that minimum orbital periods less than 11 min could be reached for binaries with an initial orbital period very close to the bifurcation period if traditional MB is included, but in a time longer than the age of the universe. van der Sluys et al. (2005b) further investigated this magnetic capture scenario for the formation of UCXBs. Our calculations show that when the initial orbital period is close to the bifurcation period, ultra-compact systems (P < 1 h) can indeed form with saturated MB, but also in a time longer than the age of the universe. For example, for an LMXB with M 2,i = 1.3M and P i = 0.46 day in model 4, a final period of P f = 22 min can be reached after 15 Gyr of mass transfer. All the works done by previous authors show that a more efficient angular momentum loss mechanism is required to produce UCXBs within 13.7 Gyr in this scenario.

8 Semi-analytical Method In this subsection we will try to use a semi-analytical method to understand our numerical results. First from Eq. (2), we have the following equation 3 J 2 J = 3 2 ( M 1 M + 2 ) 1 M 1 + M P M 1 M 2 2 M 1 + M 2 2 P. (10) If we assume a fraction α of the mass lost by the donor is accreted by the NS, i.e., M 1 = αm 2, we can write the period derivative as where P P = 3 J J A(M M 1, M 2, α) 2, (11) M 2 A(M 1, M 2, α) = 3M (1 α)m 1M 2 3αM 2 2 M 1 (M 1 + M 2 ). (12) Our analysis is limited to binary evolution with M 2 < M 1. In this case it is clearly seen that mass transfer increases the orbital period and angular momentum loss decreases the orbital period. The bifurcation period is decided by the balance of these two factors. Keeping the orbital period unchanged (i.e. P 0), we calculate the maximum mass transfer rates for orbital periods from 0.35 day to 0.95 day. This period interval covers the whole range of the bifurcation periods obtained in this work and in Podsiadlowski et al. (2002). We show the calculated mass transfer rates in Fig. 3, and find that they can be fitted by an approximate expression as M 2 (P) (P/day) 6.41(±0.11) M yr 1. (13) Then we calculate the mean mass transfer rates with constant P and find that they lie between Ṁ2(P)/2 and Ṁ2(P). This means that if we use Eq. (13) to calculate P rlof, it will deviate no more than 10% from the true value. The mean mass transfer rates here are calculated as follows. We fix the binary period P in a constant value in our code, and then evolve the donor from its initial mass M 2,i to the time when it loses half of its initial mass 0.5M 2,i. This mass transfer process takes a time of T 1/2. Then we use 0.5M 2,i /T 1/2 as the mean Ṁ2 for this period P. If we assume that angular momentum loss is dominated by saturated MB when P < 10 d, from Eqs. (2) and (7) we have J J orb = Kω 2 crit ( R 2 R ) 0.5 ( M 2 M ) 0.5 ( 4π2 G )2/3(M 1 + M 2 ) 1/3 M 1 M 2 P 4/ (P/day) 1 yr 1. (14)

9 9 Here we adopt K = gcm 2 s (Andronov et al. 2003), ω crit = Hz (Sills et al. 2000), M 1 = 1.4M, M 2 = 1M, and replace the radius of the secondary R 2 with its Roche lobe radius R L,2 from Eq. (1). Combining Eqs. (11)-(14) with P = 0 we obtain P rlof 12.4 hr for conservative mass transfer (α = 1), and P rlof 10.7 hr for non-conservative mass transfer (α = 0). These values agree well with our numerical results ( hr). If the traditional MB law is used, similarly, from Eq. (6) we get J J orb = M 2 R 4 2ω 3 G 2/3(M 1 + M 2 ) 1/3 M 1 M 2 ω 1/ (P/day) 2/3 yr 1, (15) for M 1 = 1.4M and M 2 = 1M. Combining Eqs. (11)-(13), (15) with P = 0 we get P rlof 22.3 hr for conservative mass transfer (α = 1), which is about 20% larger than 18.3 hr from our numerical calculations and 17.7 hr in Podsiadlowski et al. (2002). The main reason for this difference is that we use the constant value 1.4M, 1M for M 1, M 2 in Eq. (12), which should change with time to 2.2M, 0.2M. This will decrease the coefficient in Eq. (12) and increase the value of P rlof by 10 20%. For donor mass 1.2M, this will increase the value of P rlof by as much as 30%. So it is better to use α = 0 instead of α = 1 for donors mass 1.2M, which could yield more accurate results (from our numerical results we find that the deviation of P rlof between conservative and non-conservative mass transfer is smaller than 10%). Using Eq. (11)-(15), we also compute the semi-analytical results of P rlof for 0.5 2M donors, and compare them with our numerical results of models 1 and 4 in Fig. 4, where the semi 4 results are calculated with α = 0. When we calculate the semi 1 results in Fig. 4, for the reasons mentioned above and below, we use α = 1 for M 2,i 1.1M and α = 0 for M 2,i 1.2M, where α = 1 should be used. A few points need to be noted for the semi-analytical results in Fig. 4. Firstly with the above mentioned equations it is impossible to compute the P rlof for binaries with M 2,i 1.4M under conservative mass transfer, since both terms in the right side of Eq. (11) are negative when α = 1 and M 2,i 1.4 M, and there will be no solutions for P = 0. We instead adopt α = 0 when M 2,i > 1.4M (from our numerical results we find that the deviation of P rlof between conservative and nonconservative mass transfer is smaller than 10%). Secondly, Eq. (13) is derived only for 1M donor star rather donors in the whole mass range (0.5 2 M ), because in the latter case it is impossible to find a unified expression of the mass transfer rate like Eq. (13). As seen in Fig. 4, the difference between the semi-analytical and numerical results is generally smaller than 20% except for donors smaller than 0.7M. The reasons for the big discrepancies when M 2,i < 0.7M are discussed in 4.

10 10 4. Discussion and Conclusions Motivated by new ideas about MB and mass loss in LMXB evolution, we have made a systematic investigation on the bifurcation periods in binary models, taking into account different MB laws and mass loss mechanisms. We find that the strength of MB is the dominant factor in determining the value of bifurcation periods compared with mass loss. The stronger MB, the larger the bifurcation periods. This also results in an upper limit for the secondary masses beyond which no converging systems exist. In our calculations we assume either fully conservative (models 1 and 2) or non-conservative (models 3 and 4) mass transfer to constrain the bifurcation period distribution in different mass transfer modes. From the expression of A(M 1, M 2, α) we always have A(M 1, M 2, 1) = 3 3M 2 M 2 < A(M 1, M 2, 0) = 3, (16) M 1 M 1 + M 2 which means that non-conservative mass transfer contributes more to the increase of the orbital period than conservative mass transfer. This explains why we generally have a lower bifurcation period in non-conservative mass transfer models (models 3 and 4) than in conservative mass transfer model (model 2) under the same MB law. The real situation may lie between these two extreme cases. For binary systems with donors M 2,i M, it would take more that 13.7 Gyr before mass transfer begins via Roche lobe overflow. So the bifurcation period for these system seems meaningless, unless there exist some unknown mechanisms of loss of orbital angular momentum. In our semi-analytical analysis in 3.3 we use the condition P 0 to derive the values of P rlof. This expression seems different from P f P rlof, which is the original definition of bifurcation period. We argue here that these two expressions are roughly the same except for binaries with M 2,i 1.4M under conservative mass transfer (α = 1), the reason of which has been given in 3.3. For M 2,i < 1.4M, we find that P/P always scales with P from Eq. (11) and (13)-(15). This means that if initially P > 0 (< 0), P/P will become larger (smaller) during the evolution, leading to monotonic increase (decrease) of the period, as seen in Fig. 2. So for these systems P f P rlof is approximately equivalent with P 0. Several rough assumptions in this semi-analytical method contribute to the discrepancies between the semi-analytical results and the numerical results in Fig. 4, especially for M 2,i < 0.7M. First is the use of P 0 as the definition of P rlof, which may not work well sometimes. Second is the use of Eq. (13), which is most suitable for binaries with M 2,i = 1M as pointed out in 3.3. Third is the assumption we made that magnetic braking law is the dominated mechanism for the angular momentum loss, while the true case is that the MB may not work sometimes (for example when the convective envelop is too small). Fourth is that we use a constant initial value of M 1,i, M 2,i for M 1, M 2 in Eq. (12) and Eq.(14)-(15), while in the

11 11 true case M 1, M 2 should change with time. This will cause a big problem for α = 1 when M 2,i > 1M, which has been pointed out in 3.3. Fifth reason is the use of ω crit, in Eq. (14) as the value of ω crit for all the donors ranging from 0.5M to 1.3M. At last we conclude that (1) P is a fair definition of bifurcation period, and (2)the period evolution during the mass transfer phase is in first approximation sufficiently well described by the balance of mass transfer and angular momentum loss caused by MB. For these rough assumptions made in this semi-analytical method, its results agree with the numerical results only in the leading order. Our numerical calculations show that there is an upper limit for the donor mass beyond which no converging systems will form. Pylyser & Savonije (1988) found that, in the case of M 1,i = 4.0M, there is no converging system existing if M 2,i > 1.7M, and concluded that for any given initial accretor mass there exists a maximum initial secondary mass for the formation of converging systems. From our calculations with M 1,i = 1.4M, we find an upper limit for the initial secondary mass M 2,i between 1.2 and 1.4M under saturated MB. The reason is that for binaries with a MS donor of initial mass > 1.4M, the bifurcation period is shorter than the minimum ZAMS period, so that these systems will diverge. For traditional MB, this upper limit is > 2M, beyond the range of donor masses we adopt. We are grateful to an anonymous referee for helpful comments. This work was supported by the Natural Science Foundation of China under grant numbers and B.M. thank P. P. Eggleton, Xuefei Chen, and Xiaojie Xu for their help during this work. REFERENCES Alexander, D. R., & Ferguson, J. W. 1994, ApJ, 437, 879 Andronov, N., Pinsonneault, M., & Sills, A. 2003, ApJ, 582, 358 Bassa, C. G., van Kerkwijk, M. H., Koester, D. & Verbunt, F. 2006, A&A456, 295 Böhm-Vitense, E. 1958, Zeitschrift fur Astrophysik, 46, 108 Burderi, L., et al. 2001, ApJ, 560, L71 Burderi, L., D Antona, F., & Burgay, M. 2002, ApJ, 574, 325 Chou, Y., & Grindlay, J. E. 2001, ApJ, 563, 934 Clark, G. W., Markert, T. H., & Li, F. K. 1975, ApJ, 199, L93

12 12 D Antona, F., Ventura, P.,Burderi, L., Di Salvo, T., Lavagetto, G., Possenti, A. & Teodorescu, A. 2006, ApJ, 640, 950 Eggleton, P. P. 1971, MNRAS, 151, 351 Eggleton, P. P. 1972, MNRAS, 156, 361 Eggleton, P. P. 1983, ApJ, 268, 368 Eggleton, P. 2006, Evolutionary Processes in Binary and Multiple Stars, by Peter Eggleton, pp.. ISBN Cambridge, UK: Cambridge University Press, 2006., Ergma, E., & Sarna, M. J. 1996, MNRAS, 280, 1000 Ergma, E., Sarna, M. J., & Antipova, J. 1998, MNRAS, 300, 352 Frank, J., King, A. R., & Raine, D. 2002, Accretion Power in Astrophysics (Cambridge: Cambridge Univ. Press) Freire, P. C. C., Wolszczan, A., van den Berg, M., & Hessels, J. W. T. 2008, ApJ, 679, 1433 Jung, Y. K., & Kim, Y.-C. 2007, Journal of Astronomy and Space Sciences, 24, 1 Han, Z., Podsiadlowski, P., & Eggleton, P. P. 1994, MNRAS, 270, 121 Hubbard, W. B., & Lampe, M. 1969, ApJS, 18, 297 Kim, Y.-C., & Demarque, P. 1996, ApJ, 457, 340 Krishnamurthi, A., Pinsonneault, M. H., Barnes, S., & Sofia, S. 1997, ApJ, 480, 303 Landau, L. D., & Lifshitz, E. M. 1975, Course of theoretical physics - Pergamon International Library of Science, Technology, Engineering and Social Studies, Oxford: Pergamon Press, 1975, 4th rev.engl.ed., Lombardi, J. C., Jr., Proulx, Z. F., Dooley, K. L., Theriault, E. M., Ivanova, N., & Rasio, F. A. 2006, ApJ, 640, 441 Nelson, L. A., Rappaport, S. A., & Joss, P. C. 1986, ApJ, 311, 226 Paczynski, B., & Sienkiewicz, R. 1981, ApJ, 248, L27 Podsiadlowski, P., Rappaport, S., & Pfahl, E. D. 2002, ApJ, 565, 1107 Pols, O. R., Tout, C. A., Eggleton, P. P., & Han, Z. 1995, MNRAS, 274, 964

13 13 Pylyser, E., & Savonije, G. J. 1988, A&A, 191, 57 Pylyser, E. H. P., & Savonije, G. J. 1989, A&A, 208, 52 Rappaport, S., Verbunt, F., & Joss, P. C. 1983, ApJ, 275, 713 Rasio, F. A., Pfahl, E. D., & Rappaport, S. 2000, ApJ, 532, L47 Rogers, F. J., & Iglesias, C. A. 1992, ApJS, 79, 507 Ruderman, M., Shaham, J., & Tavani, M. 1989, ApJ, 336, 507 Sills, A., Pinsonneault, M. H., & Terndrup, D. M. 2000, ApJ, 534, 335 Steeghs, D. & Jonker, P. G. 2007, ApJ, 669, L85 Timmes, F. X., Woosley, S. E., & Weaver, T. A. 1996, ApJ, 457, 834 Tutukov, A. V., Fedorova, A. V., Ergma, E. V., & Yungelson, L. R. 1985, Soviet Astronomy Letters, 11, 52 Tutukov, A. V., Fedorova, A. V., Ergma, E. V., & Yungelson, L. R. 1987, Soviet Astronomy Letters, 13, 328 van der Klis, M., Hasinger, G., Verbunt, F., van Paradijs, J., Belloni, T., & Lewin, W. H. G. 1993, A&A, 279, L21 van der Sluys, M. V., Verbunt, F., & Pols, O. R. 2005, A&A, 431, 647 van der Sluys, M. V., Verbunt, F., & Pols, O. R. 2005, A&A, 440, 973 Verbunt, F., & Zwaan, C. 1981, A&A, 100, L7 Webbink, R. F., Rappaport, S., & Savonije, G. J. 1983, ApJ, 270, 678 This preprint was prepared with the AAS L A TEX macros v5.2.

14 14 Table 1. Calculated results of the bifurcation periods for different binary models P bif (day) P rlof (hr) M 2,i Model 1 Model 2 Model 3 Model 4 Model 1 Model 2 Model 3 Model 4 0.5M M M M M M M M M M M M M

15 Model 1 Model 2 Model 3 Model 4 Bifurcation Period (days) Secondary Mass (M O ) Fig. 1. Bifurcation periods as a function of the secondary mass in an LMXB for the four kinds of models described in 2.3. The dotted line shows the minimum initial period P ZAMS that corresponds to a Roche lobe filling zero-age main-sequence secondary star.

16 16 10 Period (days) Time (Gyr) Fig. 2. Period evolutions of an LMXB with M 2,i = 1.0M and P i = 1.20, 1.25, 1.40 day in model 3. Circle, triangle, and cross mark the onset of the Roche lobe overflow (RLOF), the end of the RLOF, and the end of the calculation, respectively.

17 Mass transfer rate (10-10 M O yr -1 ) Period (days) Fig. 3. Maximum mass transfer rates in an LMXB consisting of 1.4M neutron star and a 1M secondary at fixed orbital periods from 0.35 d to 0.95 d. Crosses marks the calculated data and the solid line represents a logarithmic fit.

18 Model 1 Semi 1 Model 4 Semi 4 P r l o f (hr) Secondary Mass (M O ) Fig. 4. Comparison of the semi-analytical results of P rlof with the numerical values for model 1 and model 4. The dotted line shows the minimum initial period P ZAMS that corresponds to a Roche lobe filling zero-age main-sequence secondary star. Here the semi 1 results are calculated with α = 1 for M 2,i 1.1M, and α = 0 for M 2,i 1.2M where α = 1 should be used. The reasons why we do this are given in the text.