Formation of Millisecond Pulsars from Intermediate- and Low-Mass X-ray Binaries
We present a systematic study of the evolution of intermediate- and low-mass X-ray binaries consisting of an accreting neutron star of mass and a donor star of mass . In our calculations we take into account physical processes such as unstable disk accretion, radio ejection, bump-induced detachment, and outflow from the point. Comparing the calculated results with the observations of binary radio pulsars, we report the following results. (1) The allowed parameter space for forming binary pulsars in the initial orbital period - donor mass plane increases with increasing neutron star mass. This may help explain why some MSPs with orbital periods longer than days seem to have less massive white dwarfs than expected. Alternatively, some of these wide binary pulsars may be formed through mass transfer driven by planet/brown dwarf-involved common envelope evolution. (2) Some of the pulsars in compact binaries might have evolved from intermediate-mass X-ray binaries with anomalous magnetic braking. (3) The equilibrium spin periods of neutron stars in low-mass X-ray binaries are in general shorter than the observed spin periods of binary pulsars by more than one order of magnitude, suggesting that either the simple equilibrium spin model does not apply, or there are other mechanisms/processes spinning down the neutron stars.
Millisecond pulsars (MSPs) are neutron stars (NSs) characterized by short spin periods ( ms) and weak surface magnetic fields ( G), which are often found in binaries with a white dwarf (WD) companion. It is generally agreed that MSPs are old NSs, recycled by accretion of mass and angular momentum from the donor stars through Roche-lobe overflow (RLOF) during the previous low-mass X-ray binary (LMXB) evolution. The NS was spun up, due to mass accretion, into a MSP, while the donor evolved to be a He or CO WD (see e.g., Bhattacharya & van den Heuvel, 1991; Tauris & van den Heuvel, 2006, for reviews).
The stability of mass transfer in an X-ray binary depends on the ratio of the masses of the donor star and the NS, and the initial orbital period of the system. Traditionally, it was thought that, if the mass ratio is large enough (), the mass transfer is likely to be unstable, resulting in a common envelope (CE) evolution (Paczynski et al., 1976; Webbink, 1984; Iben & Livio, 1993), where the NS spirals into the envelope of the donor on a very short timescale ( yr). More recent investigations (e.g., Tauris et al., 2000; Podsiadlowski & Rappaport, 2000; Kolb et al., 2000; Podsiadlowski et al., 2002; Pfahl et al., 2003) show that X-ray binary systems with intermediate-mass (up to ) donor stars (i.e., IMXBs) can avoid the spiral-in phase and experience rapid mass transfer on a thermal timescale, successfully evolve to become LMXBs.
For LMXBs, a CE evolution is even unavoidable if the separation of binary components is large enough, so that the donor reaches the asymptotic giant branch (AGB) phase with a deep convective envelope before RLOF (i.e. Case C RLOF). Stable mass transfer in LMXBs usually occurs on a timescale yr when the donor star is on the main sequence or (sub)giant branch at the onset of RLOF (i.e., Case A or B RLOF), and forming MSPs seems to be feasible in this process. It was found by Pylyser & Savonije (1988, 1989) that there exists a critical bifurcation orbital period (), the initial orbital period that separates the formation of converging LMXBs (which evolve with decreasing orbital periods until the donor star becomes degenerate or an ultra-compact binary is formed) from diverging LMXBs. The value of is found to be d, but depends heavily on the processes of tidal interactions and the mechanisms and efficiency of orbital angular momentum loss (Ergma et al., 1998; Podsiadlowski et al., 2002; van der Sluys, Verbunt, & Pols, 2005; Ma & Li, 2009).
The final products of LMXB evolution are binary pulsars. The distributions of the spin periods of the pulsars, the orbital periods and the WD masses can be used to testify the models of I/LMXB evolution. Deloye (2008) compared the theoretical expectations of I/LMXB evolution to the populations of Galactic binary pulsars. He showed that a significant population of binary pulsars with 1 d d are generally consistent with being the descendants of long-period LMXBs or IMXBs. However, there remain quite a few unresolved puzzles. For example, binary pulsars with 60 d seem to have WD companions less massive than predicted by theory, as pointed out previously by Tauris & Savonije (1999), and those with 0.1 d 1 d are inconsistent with any I/LMXB evolution.
In the previous studies on I/LMXB evolution, a canonical NS (of mass ) was usually adopted. This seems to be supported by the finding that the NS mass distribution is consistent with a narrow Gaussian at (Thorsett & Chakrabarty, 1999). However, both observations (see Zhang et al., 2011; Kiziltan, Kottas, & Thorsett, 2010; Schwab, Podsiadlowski, & Rappaport, 2010; Özel et al., 2012, and references therein) and theories (Nomoto, 1984; Timmes et al., 1996; Heger et al., 2003; Woosley & Janka, 2005) suggest that the initial masses of NSs may occupy a large range, probably originating from two different mechanisms of forming NS: iron-core collapse supernovae and electron-capture supernovae. NSs in high-mass X-ray binaries (HMXBs) have experienced very little accretion because of their young ages, so their masses should be very close to those at birth. The measured masses of NSs in HMXBs range from for SMC X1 (van der Meeret al., 2007) to for Vela X1 (Barziv et al., 2001; Quaintrell et al., 2003). Recently Rawls et al. (2011) present an improved method for determining the mass of NSs in eclipsing X-ray pulsar binaries and apply it to six systems. They find that the NS masses range from (eccentric orbit) or (circular orbit) for 4U153852 to for Vela X1. So in a proper investigation, the influence of the NS masses should be included. Recent evolutionary calculations by De Vito & Benvenuto (2010) have provided evidence that the evolution of I/LMXBs depends upon the NS mass. In this work we perform systematic calculations of I/LMXB evolution and discuss the properties of the produced binary and millisecond pulsars (BMSPs), taking into account different initial NS masses. We adopt the initial donor masses to be , and two different initial masses (1.0 and 1.8 ) for the NS.
During the mass transfer processes, part of the transferred mass from the donor star may escape from the binary system, carrying away the orbital angular momentum. The formation of BMSPs is closely related to the mechanisms of mass and angular momentum loss, and the stability of mass transfer is also dependent on the angular mementum loss rate (Soberman et al., 1997). Generally, it is assumed that all the mass transferred from the donor will accrete onto the NS unless in a super-Eddington mass transfer phase, during which the NS will accrete at the Eddington rate ( yr for a NS), and the residual mass escapes from the binary system, carrying the NS’s specific orbital angular momentum. This is the so-called “isotropic reemission model”. Actually, there may also be mass loss even during the sub-Eddington mass transfer phase. For example, the accretion disk can become thermally and viscously unstable when the orbital period is larger than a critical value (van Paradijs, 1996; Dubus et al., 1999), leading to limit cycle behavior of the mass transfer rate. The NS can accrete mass only during outbursts, while most matter may be ejected out of the binary systems during quiescence by the radiation and magnetic pressure of the rapidly rotating NS (Ruderman et al., 1989; Burderi et al., 2002). In our calculations, mass loss due to super-Eddington mass transfer, radio ejection caused by an unstable disk (Burderi et al., 2002) or bump-induced detachment (D’Antona et al., 2006), and in some cases outflow from the Lagrangian point are included.
This paper is organized as follows. In Section 2 we describe the stellar evolution code and the binary model used in this paper. We present the calculated results in Section 3 and discuss their possible applications in the formation of BMSPs in Section 4, and summarize in Section 5.
2 Binary Evolutionary Calculations
2.1 The Stellar Evolution Code
All calculations were carried out with an updated version of the stellar evolution code developed by Eggleton (1971, 1972) (see also Han et al., 1994; Pols et al., 1995). We set initial solar chemical compositions (i.e., , , and ) for the donor star, and take the ratio of the mixing length to the pressure scale height to be 2.0, and the general convective overshooting parameter to be 0.12. In our calculations, we have considered a number of binary interactions in order to follow the details of the mass transfer process, including orbital angular momentum loss due to gravitational radiation (GR), magnetic braking (MB), mass loss, and the effect of disk instability.
2.2 The Input Physics
Each of the binaries initially consists of an NS of mass ( or 1.8 ) and a zero-age main-sequence (ZAMS) donor of mass (). The efective radius of the Roche lobe of the donor is given by the following formula (Eggleton, 1983),
where is the orbital separation of the binary, and is the mass ratio of the binary components. As usual we assume that tides keep the binary orbit circular (King, 1988), so that the orbital angular momentum is
where , and is the orbital angular velocity. Logarithmic differentiation of Eq. (2) with time gives the rate of change in the orbital separation
where the total changing rate in the orbital angular momentum is
The three terms on the right-hand-side of Eq. (4) represent angular momentum loss due to GR, MB, and mass loss, respectively. The rate of angular momentum loss due to GR is calculated according to the standard formula (Landau & Lifshitz, 1959; Faulkner, 1971)
where and are the gravitational constant and the speed of light, respectively. The prescription of Verbunt & Zwaan (1981) is adopted to calculate the angular momentum loss due to MB,
where is the radius of the donor.
For IMXBs and wide LMXBs, the mass transfer rate (note that ) can be larger than the Eddington accretion rate . Thus the mass loss rate from the binary system is . Here, we adopt the isotropic reemission model, assuming that the extra material leaves the binary in the form of isotropic wind from the NS, carrying off the NS’s specific orbital angular momentum . The related angular momentum loss rate can be derived to be
In some cases part of the material lost from the donor star may escape the system through the Lagrangian point rather accrete onto the NS. Assuming that a fraction of the matter flow escapes from the binary with the specific orbital angular momentum , the angular momentum loss rate due to the point outflow is given by (Shao & Li, 2012)
where is the distance between the mass center of binary and the point.
The transferred material from the donor will form an accretion disk surrounding the NS. If the effective temperature in the accretion disk is below K (the hydrogen ionization temperature), the accretion disk is likely to be thermally and viscously unstable (Lasota, 2001). In LMXBs irradiation from the NS may help stabilize the disk to some extent (van Paradijs, 1996; King et al., 1997; Lasota, 2001; Ritter, 2008). The critical mass transfer rate for the disk instability is given by (Dubus et al., 1999)
If the mass transfer rate is less than , the X-ray binary is assumed to become transient, experiencing short outbursts separated by long quiescent intervals. Material accumulates in the disk during the quiescent phase while the NS accretes mostly during outbursts. Define the duty cycle (, see King et al., 2003) to be the ratio of the outburst timescale to the recurrence time, the accretion rate during outbursts can be estimated as . Of course, the mass accretion rate in this case is also limited by the Eddington accretion rate.
3 Results of I/LMXB Evolution Calculations
3.1 The Initial Parameter Space for Successful Evolution
In Fig. 1 we outline the results of our calculations by showing the fate of the I/LMXB evolution in the initial diagram, similar as in Tauris et al. (2000). The left and right panels correspond to and NS, respectively. It is seen that the allowed space for successful evolution into binary pulsars (i.e., without CE evolution) is larger for than for , since the lower mass ratio in the former case can stabilize mass transfer during the IMXB evolution. When , systems with a low mass companion () can form binary radio pulsars with a He/CO WD companion if the initial orbital period is d, while the lower limit of increases to when . IMXBs can avoid spiral-in and CE evolution when and d for , and and d for , respectively. In both cases, if the initial is too short, the binary systems will experience either a CE phase, or become X-ray binaries with degenerate hydrogen stars (Ergma et al., 1998), possibly forming ultra-compact binaries. On the other hand, if the initial is too long, the donor will develop a deep convective envelope at the onset of RLOF, and a runaway mass transfer is initiated, leading to CE evolution. In Tauris et al. (2000) systems with , and d can survive unstable mass transfer. These systems are right between the limiting cases presented in Fig. 1.
To see in more detail how the NS mass affects the binary evolution, we show in Fig. 2 the evolutionary paths of two IMXBs with the same donor mass and orbital period = 10 d, but different NS mass (1.0 and 1.8 ). At the age of Myr, the donor star starts to transfer mass via RLOF. In the upper two panels with a 1.0 NS, mass transfer occurs rapidly, rising up to yr within yr. Meanwhile, the orbital period reduces to d. This dynamically unstable mass transfer will result in CE evolution and probably merging of the NS and the He core of the donor. In the lower two panels with a 1.8 NS, mass transfer initially proceeds at a rate of yr for yr, during which most of the hydrogen-rich envelope () is removed from the donor star. The mass ratio inverts during this phase, and subsequently the mass transfer rate decreases to be several yr for yr. During this time the NS accretes about mass, and this is the primary process to spin up the NS. The final product is a binary consisting of a (mildly) recycled NS and a CO WD in an orbital period of d.
3.2 The Diagram
Figure 3 shows the calculated correlation between the final orbital period and the mass of the WD (the remnant of the donor) , for two different initial NS masses, (left panel) and (right panel). As mentioned above, low-mass donor stars ( or for or NS, respectively) may evolve to be He WDs. However, when the mass of He is accumulated to exceed in the core of the donor, He flash will occur, giving rise to the formation of a CO core. The orbital period will reach days. Intermediate-mass donors can avoid He flash due to their higher temperature, forming CO WDs with the mass . Their distribution deviates from that of the low-mass branch obviously.
3.3 Radio Ejection during the “Bump-related” Detachment?
When its hydrogen shell reaches a discontinuity in the hydrogen content at the time of the first dredge-up, a low-mass star will suffer a temporary contraction, thus producing a “bump” in the luminosity function of the red giants. D’Antona et al. (2006) suggested that the orbital period gap ( days) of BMSPs may be related to the bump-induced detachment of the donor star from its RL. If the NS has already accreted sufficient mass, it may turn on as a MSP during the detachment phase. Material transferred from the donor, once it expands again to refill its Roche lobe, may be inhibited by the pulsar’s radiation pressure and ejected at the inner lagrangian point (Ruderman et al., 1989; Burderi et al., 2002), so that no further mass accretion would occur. Similarly, this so-called “radio ejection” process may also occur in LMXBs with an unstable accretion disk (Burderi et al., 2002). In both cases the mass loss will influence the orbital evolution of the binary and the spin evolution of the NS. D’Antona et al. (2006) considered the radio ejection only during the bump-related detachment. Our calculations show that the disk instability generally appears earlier than the bump phase. Two examples are presented in Fig. 4. In the upper two panels, the binary system initially consists of a 1.0 NS and a 1.0 ZAMS donor star in an orbital period of 5.0 d. RLOF initiates at the age yr, and the bump occurs yr later, at which the donor mass has decreased to be . However, the disk becomes thermally unstable when the donor mass is , much earlier than the bump phase. If the NS’s spin was accelerated during the prior mass transfer, the accumulated mass in the disk during quiescence will be ejected out of the binary by the radiation pressure of the rapidly spinning NS. We find a similar situation in the case of a binary containing a 1.8 NS and a 1.5 donor star shown in the lower two panels. Compared with the effect of unstable disk accretion, it seems that the bump-related detachment might play a less important role in the evolution of LMXBs.
We finally note that the process of pulsar-driven mass ejection is highly uncertain, since its condition and efficiency depend on several unknown parameters (Fu & Li, 2011), such as the value of the equilibrium period that a NS will reach during the mass transfer. However, as seen below, there is controversy on the estimate of the equilibrium period for NSs in LMXBs, and one should be cautious when considering the effect of radio ejection on the binary evolution.
3.4 The effect of outflow from the point
In the binary evolution the stability of mass transfer strongly depends on mass loss and related angular momentum loss. Although we have taken into account various ways for mass loss, including super-Eddington mass transfer, unstable mass transfer due to disk instability, and mass ejection due to the propeller effect and the pulsar’s radiation, the detail processes are complicated and uncertain, and simplified treatment might not reveal the realistic situation. In the following, we show an example of how outflows from the binary systems through the point influence the evolution of mass transfer, and change the conditions of forming BMSPs.
MSP J16142230 is a 3.15 ms pulsar of mass 1.97 0.04 (Demorest, 2010). The mass of its WD companion () and the orbital period (8.7 d) were also measured accurately. Lin et al. (2011) and Tauris et al. (2011) systematically investigated the formation channels of this pulsar from an IMXB. They showed that NS this massive are not easy to produce in spite of the initially high mass of the donor star, unless they were already born as a relatively massive NS. However, Tauris et al. (2011) found that, for the system with a 1.8 NS, the final orbital period will be always larger than 10 d, inconsistent with the observation of the pulsar. This conclusion will not hold if we let a small fraction of the transferred mass leave the system from the point (e.g., Bailyn et al., 1989). As an illustration, in Fig. 5 we compare the evolution of an IMXB containing a 1.8 NS and a 4.5 companion star in an orbit of 2.4 d, with (upper two panels) and (lower two panels), respectively. In the former case the primary angular momentum loss mechanism is the isotropic reemission around the NS during the super-Eddington accretion phase. The final system consists of a 2 NS and 0.5 CO WD, with a orbital period 15 d. In the latter case angular momentum loss due to the point outflow plays a role when the system evolves into an LMXB, and reduces the orbital period to days.
4 Comparison with BMSPs
4.1 The relation for low-mass BMSPs
In wide LMXBs, the donor will climb to the red-giant branch (RGB) in the HR diagram before RLOF. For low-mass stars () on the RGB, there is a well known relationship between the mass of the degenerate He core and the radius of the giant star, which is almost entirely independent of the mass of the hydrogen-rich envelope (Refsdal & Weigert, 1971; Webbink et al., 1983). Based on this relationship, a specific correlation between the orbital period and the WD mass is obtained (Joss et al., 1987; Rappaport et al., 1995; Tauris & Savonije, 1999). Comparison with the observations shows that a significant population of BMSPs with He WD companion is generally consistent with this relation. However, there seems to be a systematic deviation from the correlation for pulsars with d, which seem to have WD companions lighter than expected (Tauris, 1996).
Both systematic small values of the orbital inclination and large NS mass can increase for the given observed mass functions. Since there does not seem to be any observational selection effect favoring small inclination angle (Tauris, 1996), we first examine whether the correlation can be accounted for if the long-period BMSPs have CO WD companions. For example, Stairs et al. (2005) noticed that the Tauris & Savonije (1999) relation is incompatible at the level with a uniform distribution of if the pulsar masses are drawn from a Gaussian distribution centered on with width , and better agreement with uniformity in can be reached if the pulsar masses are large on average (e.g., ). An extreme example is PSR B082002, which has a CO WD companion in a very wide orbit with d (Koester & Reimers, 2000). However, in such wide binaries it is difficult for the NS to accrete enough matter (Li & Wang, 1998; Tauris & Savonije, 1999), and the NS must be born heavy if this interpretation is correct. In Fig. 6 we compare the relations between and in the cases that the NS has an initial mass of 1.0 (left panel) and 1.8 (right panel). Also plotted are binary pulsars with measured and (90% probability mass range for randomly oriented orbits) for a fixed NS mass of 1.2 and 2.0 , respectively (i.e., in each case the NS is assumed to have accreted during the mass transfer). It is seen that some binary pulsars with 100 d can fairly match the relation if they have massive NSs () and heavy CO WDs (although in some cases small orbital inclination angles may be required), while for those with 20 d, statistically lighter NSs () seem to follow the relation better.
It is expected that the NS has to accrete at least a few mass to reach a millisecond period, but this is difficult to achieve for NSs in wide binaries (see also Liu & Chen, 2011), because the mass transfer rate (which increases with increasing orbital period) is likely to be super-Eddington, and the accretion disk is likely to be unstable. In Fig. 7 we show the mass transfer rate , and the accreted mass of the NS as a function of the final orbital period . Note that here we plot the mass transfer rate only for the stable mass transfer phase - when the accretion disk becomes unstable, we take it to be the critical value at the onset of disk instability. Hence both the mass accretion rate and are limited by the value of : when is less than the , unstable disk accretion occurs, and the NS is assumed to accrete mass only during outbursts, and part of the matter will be ejected out of the binary systems if is super-Eddington.
The left panel of Fig. 7 shows the results for a NS with a companion star of initial mass , 2.0 and 3.0 , respectively. Generally more massive donor stars result in higher mass transfer rate, which also increases with longer orbital period . It is noted that (1) is always , and (2) , which decreases with , because in wider systems the mass accretion rates (with both stable and unstable accretion disks) are more likely to be super-Eddington. The right panel is for the systems with a NS. In the case of light companion star (), the initial orbital period should be larger than d, so that the critical mass transfer rate is always larger than the mass transfer rate, since increases during the mass transfer process (with increasing orbit period), while the mass transfer rate decreases all the way (Webbink et al., 1983). The accreted mass by the NS is thus significantly lower than those with more massive donors. So only results with 2.0 and 3.0 donor stars are presented, in which , and .
Figure 7 shows that in general , and . The distribution of the mass transfer rate is considerably higher than that ( yr) obtained by Podsiadlowski et al. (2002), but more compatible with the observations of persistent LMXBs. The main reason is that we have taken into account the effect of disk instability. The distribution is roughly in line with Liu & Chen (2011), who found that is generally less than 0.6 in their calculations for systems with a 1.4 NS and a donor star. Moreover, we find that systems with initially massive NSs () may accrete enough mass to evolve into BMSPs with 60 d d, , while light NSs in wide binaries are more likely to be partially recycled. This seems to support our conjecture that some of the wide BMSPs might be born massive. Hopefully accurate measurements of both the pulsar and the WD masses will help settle this issue.
Löhmer et al. (2005) performed timing observations of the BMSP J16402224 (with an orbital period of 175 d), and constrained the WD mass to be ( uncertainties), which indicates that the companion is very likely to be a low-mass He WD. In this case the massive NS + CO WD model obviously does not work, and one has to explore other possibilities. Evaporation of the companion star from a wind of relativistic particles after the pulsar turns on may decrease the companion mass significantly, but it is unlikely for wide BMSPs, since the evaporation timescale would be longer than the Hubble time (Tauris, 1996).
Here we suggest another possible solution.
It is interesting to note that solar-type stars are
usually found to be surrounded by sub-stellar companions
(usually planets and/or brown dwarfs) (Cassan et al., 2012).
One may expect that in some relatively wide
LMXBs the companion star had
possessed substellar companion(s) in close orbits
like “hot Jupiters”.
When the star evolved on the giant branch it would become
big enough to capture its planet/brown dwarf.
The planet/brown dwarf spiraled into the envelope of the
giant to initiate a CE phase. The
frictional drag arising from its motion through the
CE would lead to loss of its orbital angular
momentum and deposit of orbital energy in the envelope.
If there was enough orbital energy, the spiral-in process would
expel the envelope of the giant, leaving a WD remnant
(Nelemans & Tauris, 1998)
During the planet/brown dwarf-involved CE phase, the companion’s envelope expanded rapidly and filled its Roche-lobe, leading to mass transfer onto the NS. Since this phase was very short ( yr), the mass transfer rate would be much higher than the Eddington limit rate for the NS, so that the NS accreted very small mass (unless the CE phase lasted much longer time). If this is the case, we have to require that some MSPs were born this way, rather recycled during the LMXB evolution. In the literature, this idea has already been discussed by Miller & Hamilton (2001), who showed that, the existence of the innermost, moon-sized planet in the PSR 125712 system suggests that the pulsar was born with approximately its current spin frequency and magnetic field. A schematic view of the formation of BMSPs with planet/brown dwarf-involved CE evolution is shown in Fig. 8.
Not only having an impact on the formation of MSPs, the substellar objects, if really exist in low-mass binaries, may also play a role in the CE evolution, especially influence the estimate of the CE efficiency parameter . Recently Davis et al. (2012) reconstructed the CE phase for the current sample of post-CE binaries (PCEBs) with observationally determined component masses and orbital periods. Searching for correlations between and the binary parameters, they found that, when the internal energy of the progenitor primary envelope is taken into account, decreases with increasing mass of the primary (i.e., the progenitor of the WD) (see however De Marco et al., 2011), and for , which seems to be in contrast with derived by Ricker & Taam (2012) from numerical simulations. If there are planets/brown dwarfs around these low-mass primaries, they can contribute extra orbital energy to help expel the primary’s envelope during the RGB/AGB phase, making the CE efficiency parameter within the canonical range .
4.2 The Rebirth Periods of BMSPs
Accretion onto the NS in I/LMXBs changes both the mass and spin of the NS. The spin evolution of an accreting NS depends on the interaction between the magnetosphere and the accretion disk, as described as follows (e.g., Bhattacharya & van den Heuvel, 1991): the accretion disk is truncated at the magnetospheric radius , which is close to the Alfvén radius of the disk
where and are the surface magnetic field and the radius of the NS, respectively. Material is channeled to the NS at the inner radius, producing a spin-up torque. If the NS is spinning very rapidly, the disk may be truncated outside the corotation radius , and the NS experiences a centrifugal barrier that can inhibit accretion (i.e., the so-called “propeller” effect; Illarionov & Sunyaev, 1975). So the NS will eventually reach the equilibrium spin period such that , or
where G, and cm. This period can be regarded as the beginning or rebirth period of recycled pulsars when accretion terminates.
The rebirth period of a MSP can be derived only when its actual spin-down time is known. Bhalerao & Kulkarni (2011) recently reported the optical discovery of the companion to the MSP J16142230. The optical colors show that the companion is a 2.2 Gyr old CO WD. From the age of the WD, Bhalerao & Kulkarni (2011) calculated the period of the pulsar at birth, and found that pulsar should be born with a spin close to its current value, implying that the final accretion rate was . This value is two orders of magnitude smaller than the estimate from theoretical calculations by Lin et al. (2011) and Tauris et al. (2011). These authors suggested that the system began as an IMXB consisting of a NS and a main-sequence secondary, which evolved to be a CO WD with He envelope. The NS gained the most mass during the final LMXB phase lasting Myr at near-Eddington rates.
We present a systematic view on the relation
and derived from theory and observations in Fig. 9.
The solid lines outline the theoretically
expected distribution of
(, here is
taken to be 1) from binary evolution calculations, and
the symbols represent the observations of recycled pulsars
(data are taken from the ATNF Pulsar Catalogue; Manchester et al., 2005).
Comparison between observations and theory shows
that, except in a few cases,
the observational is generally
larger than by orders of
magnitude, challenging the simple recycling theory
(1) There is no definite spin equilibrium since GR may remove the angular momentum from the NS. This idea was first suggested by Papaloizou & Pringle (1978) to account for the cutoff in the spin distribution of NSs in LMXBs. The main emission mechanisms involve crustal mountains (Bildsten, 1998), magnetic deformations (Cutler, 2002), and unstable -modes (Andersson, 1998) in the NSs. All these processes can produce a substantial mass quadrupole moment and thus a spin-down torque due to GR. The problem with this interpretation is that, recent observational results on some millisecond X-ray pulsars have shown that the efficiency of GR induced spin-down might be too low to be responsible for balancing the spin-up process during outbursts (see Haskell & Patruno, 2011; Patruno et al., 2012, and references therein). So the following models focus on modifications of the equilibrium period (Eq. ) that the NS will finally reach.
(2) If the NS magnetic field lines can thread the accretion disk, an extra magnetic torque will be exerted on the NS (Ghosh & Lamb, 1979a), and the equilibrium period will be times that in Eq. (10) (Ghosh & Lamb, 1979b; Wang, 1995; Li & Wang, 1996). Andersson et al. (2005) further suggested that the inner disk region may be geometrically thick and sub-Keplerian, and dominated by radiation pressure, if the mass transfer rate is above a few percent Eddington accretion rate. The coupling between the disk and magnetic field can reduce the amount of angular momentum deposited onto the NS from accretion by a factor (the ratio of orbital angular velocity to Keplerian rotational speed in the disk, ), compared with the case of thin disks (), thus a decreased spin-up torque results. Patruno et al. (2012) show that the existence of spin equilibrium as set by the diskmagnetosphere interaction is able to explain the observations of millisecond X-ray pulsars, if the spin-down torque coming from the interaction between the disk and field outside corotation is sufficiently large. However, the radial extent of the coupling between the star and the disk is highly controversial (e.g. Matt & Pudritz, 2005; Ghosh, 2007). Stellar and disk winds may also take away the angular momentum of the NS (Goodson et al., 1997, 1999; Ustyugova et al., 2006; Romanova et al., 2009; Zhang & Li, 2010). In addition, the abrupt torque reversals observed in disk-fed X-ray pulsars Her X-1, Cen X-3, and 4U 162667 (Bildsten et al., 1997) suggest that current models of diskmagnetosphere interaction may have severe limitations.
(3) If Eq. (10) does apply, a decaying mass transfer rate at the end of the mass transfer process will result in relatively large rebirth period of BMSPs, provided that the evolution time is longer than the NS spin evolution timescale (Jeffrey, 1986; Ruderman et al., 1989). More recently Tauris (2012) suggested that this equilibrium will even be broken when decreases substantially, so that the magnetospheric radius becomes larger than the corotation radius , causing the NS to enter the propeller phase. A centrifugal barrier arises to expel matter entering the magnetosphere with a braking torque
(where is the Keplerian angular velocity at ) to act to slow down the NS. An alternatively appropriate form for the accretion torque is (e.g. Menou et al., 1999)
where is the angular velocity of the NS.
With Eqs. (11) and (12) we recalculate the spin evolution of the NS when the
donor star gradually decouples from its RL as in
Tauris (2012). In Fig. 10 the blue and red
lines denote the results calculated with Eqs. (11) and (12),
(4) A decaying magnetic field may also change the rebirth period of the BMSPs, if the field decay timescale is shorter than the spin evolution timescale, so that the final spin period deviates from the equilibrium period. During the decay of the NS magnetic field induced by accretion, the spin-up timescale increases as decreases because the accreting matter carries less specific angular momentum, while the field decay timescale decreases as decreases (field decay stops when ), resulting in a significant departure from the spin-up line (Burderi et al., 1996; Konar & Bhattacharya, 1999).
(5) The rebirth periods of BMSPs may be related to the long-term spin equilibrium with varying mass accretion rate. A significant fraction of LMXBs are likely to be transient systems subject to disk instability (Coriat et al., 2012). For example, an NS LMXB in a 10 hr orbit would be transient if the average mass transfer rate is lower than gs. In addition, in narrow LMXBs, X-ray irradiation of the donor star could destabilize the mass transfer, and lead to mass transfer cycles (Hameury et al., 1993) - mass transfer is spasmodic with phases of high mass transfer driven by the thermal expansion of the convective envelope of the irradiated donor alternating with phases with low or no mass transfer, during which the donor readjusts towards thermal equilibrium of the un-irradiated star. The final spin period is thus determined by the spin-up during outbursts/high mass transfer phase balanced by the spin-down during quiescence/low mass transfer phase, which may be considerably longer than that attained via stable mass accretion (Li et al., 1998; Patruno et al., 2012).
Currently our knowledge about the evolution of LMXBs and accretion disks is not sufficient to tell which one is the dominant factor in determining the rebirth periods of BMSPs. Perhaps they are influenced by a combination of (at least some of) the above-mentioned mechanisms.
4.3 Origin of Intermediate-Mass BMSPs with 1 d
It is well known that a large fraction of IMXBs will produce binary pulsars with a CO WD companion. As shown in Fig. 3, the orbital periods of these intermediate-mass binary pulsars (IMBPs) are generally d. This leaves the puzzle that the binary pulsars with 1 d and a CO WD companion are hard to explain reasonably in any binary population (Deloye, 2008). In Table 1 we list the parameters of IMBPs with 1 d. Except PSR B065564, other pulsars share the characteristics of rapid spin and low magnetic field.
The small orbital periods combined with massive WDs suggest that this type of systems may have evolved through a CE and spiral-in phase (van den Heuvel & Taam, 1984). Similar to that of the Double NSs, their direct progenitors could be binaries consisting of a He star and a NS star in a close and circular orbit, as a result of spiral-in of a wide binary in late Case B or Case C mass transfer. The core of the companion was a He star which could already have gone through quite some He burning. After the spiral-in one then had a He star with already CO in its core. During He-shell burning, the envelopes of the He stars, if their masses are , will expand and slowly transfer mass to the NS, which will have caused the spin-up and recycling of the NS. This interpretation seem to be responsible for PSR B065564 (van den Heuvel & Taam, 1984), which has the longest spin period (195.6 ms) and highest magnetic field ( G) among the pulsars in Table 1, implying that it has been partially recycled with little mass accreted. For other pulsars, the short periods and low magnetic fields ( G) require that at least some material has been accreted by the NSs. Recently Chen et al. (2011) presented the detailed binary evolution calculation for a binary consisting of a NS (of mass 1.3 ) and a low-mass He star (of mass 1.0 ) with an initial orbital period of 0.5 d. They showed that the mass transfer seems to be able to spin up the NS’s spin to milliseconds, producing BMSPs like PSR J18022124, which has a short orbital period 1 d.
Another possible way to form compact IMBPs invokes evolution of IMXBs with anomalous MB. Traditionally MB is thought to work in low-mass, main-sequence stars. However, There are many intermediate-mass stars which have anomalously high magnetic fields ( 100 G), i.e., Ap/Bp stars. In IMXBs, the irradiation-driven stellar wind from the surface of the Ap/Bp donor star by the NSs may couple with the high magnetic fields, resulting in angular momentum loss at a rate (Justham et al., 2006),
where is the magnetic field strength at the surface of the donor star, and is a parameter with typical value , combining the wind-driving energy conversion efficiency and irradiation geometry. Justham et al. (2006) explored the effect of anomalous MB in the formation of compact BH LMXBs, in which the orbital energy of the companion star, if of low-mass initially, would be unable to expel the envelope of the BH progenitor. For the black hole LMXB XTE J1118480, González Hernández et al. (2012) measured a decay rate of the orbital period ms yr. This value is much larger than predicted by most conventional MB and mass loss models for the binary parameters of XTE J1118480, unless the companion star has a surface magnetic field kG to enhance MB.
In Fig. 11 we show the evolution of a binary initially consisting of a 1.8 NS and a 3.5 donor with d. We assume that the donor possesses a surface magnetic field G. RLOF starts at the age yr with a mass transfer rate yr. Because the donor is initially more massive than the NS, the orbital period decreases rapidly. Most of the transferred mass is assumed to be ejected out of the binary in the vicinity of the NS, taking away both mass and orbital angular momentum. The orbit starts to widen when . When the donor mass decreases to be , the mass transfer rate reduces to be yr. In the subsequent yr, MB begins to dominate the evolution. The orbital period keeps decreasing, the NS accretes most mass in this phase and is spun up to milliseconds. The final outcome is a tight binary system consisting of a NS and a WD with 0.4 d.
The anomalous MB model will not work if the NS is initially of low-mass (), since the allowed parameter space for forming binary pulsars will be very small (see Fig. 1). Additionally we note that models involving accretion-induced collapse of an ONeMg WD cannot account for the distribution of these BMSPs (Sutantyo & Li, 2000).
This work is motivated by the fact the evolution of I/LMXBs with canonical NSs seems to meet difficulties in explaining some of the observational characteristics of BMSPs, and the measurements of the NS masses indicate a wide distribution . We have preformed numerical calculations of the evolution of I/LMXBs consisting of a 1.0 or NS and a donor star, to investigate its dependence on the initial NS mass, and the properties of the descendent BMSPs. The main results can be summarized as follows.
1. The allowed parameter space in the initial diagram for forming recycled pulsars increases with increasing NS mass. This may help explain the formation of BMSPs with d and their distribution in the diagram. Alternatively, some of these wide binary pulsars may be formed through mass transfer driven by planet/brown dwarf-involved CE evolution.
2. The equilibrium spin periods of accreting NSs in LMXBs derived from the standard magnetosphere-accretion disk interaction model are in general shorter than the observed spin periods of BMSPs by more than one order of magnitude. This implies that either the simple equilibrium spin model does not apply for the spin evolution in accreting NSs, or there are other mechanisms/processes to spin down the NSs when forming BMSPs.
3. Some of the compact IMBPs might have evolved from IMXBs in which the companion star were strongly magnetized Ap/Bp stars with enhanced MB.
4. Our calculations doubt the suggestion that the orbital period gap ( d) of BMSPs is related to the occurrence of the bump-related detached phase in the LMXB evolution, since the accretion disk would become unstable at earlier time.
- In addition, Bear & Soke (2010) suggest that the binary systems may reach stable synchronized orbits before the onset of the CE phase. Such stable synchronized orbits allow the RGB star to lose mass prior to the onset of the CE phase. Even after the secondary enters the giant envelope, the rotational velocity is high enough to cause an enhanced mass-loss rate.
- Here we assume that the current spin periods of MSPs are not far from their initial ones. This holds if their characteristic ages are longer than the cooling ages of the WDs (e.g. Lorimer et al., 1995) or the Hubble time.
- The discontinuity at the point where the two lines meets is due to the fact that we use Eq. (12) to calculate the spin evolution even before the propeller phase in both cases.
- footnotetext: Data are taken from the ATNF pulsar catalogue.
- footnotetext: Median mass for the inclination angle .
- Andersson, N. 1998, ApJ, 502, 708
- Andersson, N., Glampedakis, K., Haskell, B., & Watts, A. L. 2005, MNRAS, 361, 1153
- Bailyn, C. D., Garcia, M. R., & Grindlay, J. E. 1989, ApJ, 344, 786
- Barziv, O., Kaper, L., Van Kerkwijk, M. H., Telting, J. H., & van Paradijs, J. 2001, A&A, 377, 925
- Bear, E. & Soker, N. 2010, New Astronomy, 15, 483
- Bhalerao, V. B. & Kulkarni, S. R. 2011, ApJ, 737, L1
- Bhattacharya, D. & van den Heuvel, E. P. J. 1991, Phys. Rev., 203, 1
- Bildsten, L. 1998, ApJ, 501, L89
- Bildsten, L., Chakrabarty, D., Chiu, J., Finger, M. H., & Koh, D. T. et al. 1997, ApJS, 113, 367
- Burderi, L., D’Antona, F., & Burgay, M. 2002, ApJ, 574, 325
- Burderi, L., King, A. R., & Wynn, G. A. 1996, MNRAS, 283, L63
- Cassan, A., Kubas, D., Beaulieu, J.-P., Dominik, M., Horne, K. et al. 2012, Nature, 481, 167
- Chen, W.-C., Li, X.-D., & Xu, R.-X., 2011, A&A, 530, 104
- Coriat, M., Fender, R. P., & Dubus, G. 2012, arXiv:1205.5038
- Cutler C., 2002, Phys. Rev. D, 66, 084025
- D’Antona F., Ventura P., Burderi L., Di Salvo T., & Lavagetto G. , et al. 2006, ApJ, 640, 950
- Davis, P. J., Kolb, U., & Knigge, C. 2012, MNRAS, 419, 287
- Deloye C. J. 2008, in 40 YEARS OF PULSARS: Millisecond Pulsars, Magnetars and More. AIP Conference Proceedings, Eds. C. Bassa, Z. Wang, A. Cumming, & V. M. Kaspi, Vol. 983, p. 501
- De Marco, O., Passy, J., Moe, M., Herwig, F., & Mac Low, M. et al. 2011, MNRAS, 411, 2277
- Demorest P. B., Pennucci T., Ransom S. M., Roberts M. S. E., & Hessels J. W. T., 2010, Nature, 467, 1081
- De Vito, M. A. & Benvenuto, O. G., 2010, MNRAS, 401, 2552
- Dubus, G., Lasota, J.-P., Hameury, J.-M., & Charles, P. 1999, MNRAS, 303, 139
- Eggleton, P. P. 1971, MNRAS, 151, 351
- Eggleton, P. P. 1972, MNRAS, 156, 361
- Eggleton, P. P. 1983, ApJ, 268, 368
- Ergma, E., Sarna, M. J., & Antipova, J. 1998, MNRAS, 300, 352
- Faulkner, J. 1971, ApJ, 170, L99
- Ferdman, R. D., Stairs, I. H., Kramer, M., McLaughlin, M. A., & Lorimer, D. R. et al. 2010, ApJ, 711, 764
- Fu, L. & Li, X.-D. 2011, RAA, 11, 1457
- Ghosh, P. 2007, Rotation and Accretion Powered Pulsars (World Scientific, London), p. 772
- Ghosh, P. & Lamb, F. K. 1979a, ApJ, 232, 259
- Ghosh, P. & Lamb, F. K. 1979b, ApJ, 234, 296
- González Hernández, J., Rafael Rebolo, R. & Casares, J. 2012, ApJ, 744L, 25
- Goodson, A. P., Böhm, K. H., & Winglee, R. M. 1999, ApJ, 524, 142
- Goodson, A. P., Winglee, R. M., & Böhm, K. H. 1997, ApJ, 489, 199
- Hameury, J. M., King, A. R., Lasota, J. P., & Raison, F. 1993, A&A, 277, 81
- Han, Z., Podsiadlowski, P., & Eggleton, P.P, 1994, MNRAS, 270, 121
- Haskell, B. & Patruno, A. 2011, ApJ, 738, 14
- Heger, A., Fryer, C. L., Woosley, S. E., Langer, N. & Hartmann, D. H. 2003, ApJ, 591, 288
- Iben, Jr. & Livio, M., 1993, PASP, 105, 1373
- Jeffrey, L. C. 1986, Nature, 319, 384
- Illarionov, A. F. & Sunyaev, R. A. 1975, A&A, 39, 185
- Joss, P. C., Rappaport, S. A., & Lewis, W. 1987, ApJ, 319, 180
- Justham, S., Rappaport, S., & Podsiadlowski, P. 2006, MNRAS, 366, 1415
- King, A. R. 1988, QJRAS, 29, 1
- King, A. R., Kolb, U., & Sienkiewicz, E. 1997, ApJ, 488, 89
- King A. R., Rolfe D. J., Kolb U., & Sshenker K., 2003, MNRAS, 341, L35
- Kiziltan, B., Kottas, A., & Thorsett, S. E. 2010, arXiv:1011.4291
- Koester, D. & Reimers, D. 2000, A&A, 364, L66
- Kolb, U., Davies, M. B., King, A., & Ritter, H. 2000, MNRAS, 317, 438
- Konar, S. & Bhattacharya, D. 1999, MNRAS, 303, 588
- Landau, L. D. & Lifshitz, E. M. 1959, The Classical Theory of Fields (Oxford: Pergamon Press)
- Lasota, J.-P. 2001, New Astron. Rev., 45, 449
- Li, X.-D., van den Heuvel, E. P. J., & Wang, Z.-R. 1998, ApJ, 497, 865
- Li, X.-D. & Wang, Z.-R. 1996, A&A, 307, L5
- Li, X.-D. & Wang, Z.-R. 1998, ApJ, 500, 935
- Lin, J., Rappaport, S., Podsiadlowski, P., Nelson, L., Paxton, B., & Todorov, P. 2011, ApJ, 732, 70
- Liu, W.-M. & Chen, W.-C, 2011, MNRAS, 416, 2285
- Löhmer, O., Lewandowski, W., Wolszczan, A., & Wielebinski, R. 2005, ApJ, 621, 388
- Lorimer, D. R., Lyne, A. G., Festin, L., & Nicastro, L. 1995, Nature, 376, 393
- Ma, B. & Li, X.-D. 2009, ApJ, 691, 1611
- Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 1993
- Matt, S. & Pudritz, R. E. 2005, ApJ, 632, L135
- Maxted, R. F. L., Napiwotzki, R., Dobbie, P. D., & Burleigh, M. R. 2006, Nature, 442, 543
- Menou, K., Esin, A. A., Narayan, R., Garcia, M. R., Lasota, J.-P., & McClintock, J. E. 1999, ApJ, 520, 276
- Miller, M. C. & Hamilton, D. P. 2001, ApJ, 550, 863
- Nelemans, G. & Tauris, T. M. 1998, A&A, 335, L85
- Nomoto, K. 1984, ApJ, 277, 791
- Özel, F., Psaltis, D., Narayan, R., & Villarreal, A. S. 2012, ApJ submitted (arXiv:1201.1006)
- Paczynski, B. 1976, In Structure and Evolution in Close Binary Systems. Proc. IAU Symp. 73, Eds., Eggleton, P. P., Mitton, S., Whealan, J. (Reidel, Dordrecht), p. 75
- Papaloizou, J. & Pringle, J. E. 1978, MNRAS, 184, 501
- Patruno, A., Haskell, B., & D’Angelo, C. 2012, ApJ, 746, 9
- Pfahl, E., Rappaport, S., & Podsiadlowski, P. 2003, ApJ, 597, 1036
- Podsiadlowski, Ph. & Rappaport, S. 2000, ApJ, 529, 946
- Podsiadlowski, Ph., Rappaport, S., & Pfahl, E. D. 2002, ApJ, 565, 1107
- Pols, O.R., Tout, C.A., Eggleton, P.P., & Han Z., 1995, MNRAS, 274, 964
- Pylyser, E. & Savonije, G. J. 1988, A&A, 191, 57
- Pylyser, E. & Savonije, G. J. 1989, A&A, 208, 52
- Quaintrell, H., Norton, A. J., Ash, T. D. C., Roche, P., & Willems, B. et al. 2003, A&A, 401, 313
- Rappaport, S.A., Podsiadlowski, P., Joss, P.C., Di Stefano, R., & Han, Z. 1995, MNRAS, 273, 731
- Rawls, M. L., Orosz, J. A., McClintock, J. E., Torres, M. A. P., & Bailyn, C. D. et al. 2011, ApJ, 730, 25
- Refsdal, S. & Weigert, A. 1971, A&A, 13, 367
- Ricker, P. M. & Taam, R. E. 2012, ApJ, 746, 74
- Ritter, A. 2008, New A Rev., 51, 869
- Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., & Lovelace, R. V. E. 2009, MNRAS, 399, 1802
- Ruderman, M., Shaham, J., & Tavani, M. 1989, ApJ, 336, 507
- Schwab, J., Podsiadlowski, P., & Rappaport, S. 2010, ApJ, 719, 722
- Shao, Y. & Li, X.-D. 2012, ApJ, 745, 165
- Soberman, G. E., Phinney, E. S., & van den Heuvel, E. P. J. 1997, A&A, 327, 620
- Stairs, I. H., Faulkner, A. J., Lyne, A. G., Kramer, M., & Lorimer, D. R. et al. 2005, ApJ, 632,1060
- Sutantyo, W. & Li, X.-D. 2000, A&A, 360, 633
- Tauris, T. M. 1996. A&A, 315, 453
- Tauris, T. M. 2012, Science, 335, 561
- Tauris, T. M., Langer N., & Kramer M. 2011, MNRAS, 416, 2130
- Tauris, T. M., & Savonije, G. J. 1999, A&A, 350, 928
- Tauris, T. M. & van den Heuvel E. P. J., 2006, in Compact Stellar X-Ray Sources, Eds., W. H G. Lewin, M. van der Klis (Cambridge: Cambridge University Press), 623
- Tauris, T. M., van den Heuvel, E. P. J., & Savonije, G. J. 2000, ApJ, 530, L93
- The Fermi LAT Collaboration, 2011, Science, 334, 1107
- Thorsett, S. E. & Chakrabarty, D. 1999, ApJ, 512, 288
- Timmes, F. X., Woosley, S. E., & Weaver, Thomas A. 1996, ApJ, 457, 834
- Ustyugova, G. V., Koleoba, A. V., Romanova, M. M., & Lovelace, R. V. E. 2006, ApJ, 646, 304
- van den Heuvel, E. P. J. & Taam, R. 1984, Nature, 309, 235
- van der Meer, A., Kaper, L., van Kerkwijk, M. H., Heemskerk, M. H. M., & van den Heuvel, E. P. J. 2007, A&A, 473, 523
- van der Sluys, M. V., Verbunt, F., & Pols, O. R. 2005, A&A, 440, 973
- van Paradijs, J. 1996, ApJ, 464, L139
- Verbunt, F. & Zwaan, C. 1981, A&A, 100, L7
- Wang, Y.-M. 1995, ApJ, 449, L153
- Webbink, R. F. 1984, ApJ, 277, 355
- Webbink, R. F., Rappaport, S. A., & Savonije, G. J. 1983, ApJ, 270, 678
- Woosley, S. & Janka, T. 2005, Nat. Phys., 1, 147
- Zhang, C. M., Wang, J., Zhao, Y. H., Yin, H. X., & Song, L. M. et al. 2011, A&A, 527, A83
- Zhang, Z. & Li, X.-D. 2010, A&A, 518, A19