Go to previous page Go up Go to next page

3.4 The population of relativistic binaries

Although no radio pulsar has so far been observed in orbit around a black hole companion, we now know of several double neutron star and neutron star-white dwarf binaries which will merge due to gravitational wave emission within a reasonable timescale. The current sample of objects is shown as a function of orbital period and eccentricity in Figure 18View Image. Isochrones showing various coalescence times tg are calculated using the expression
( )8/3 ( ) -2/3( )- 1 t -~ 9.83 × 106 yr Pb- m1--+-m2- -m-- (1 - e2)7/2 , (7) g hr Mo . Mo .
where m1 and m2 are the masses of the two stars, m = m1m2/(m1 + m2) is the so-called “reduced mass”, Pb is the binary period and e is the eccentricity. This formula is a good analytic approximation (within a few percent) to the numerical solution of the exact equations for tg [245246].
View Image

Figure 18: The relativistic binary merging plane. Top: Orbital eccentricity versus period for eccentric binary systems involving neutron stars. Bottom: Orbital period distribution for the massive white dwarf-pulsar binaries.
In addition to testing general relativity through observations of these systems (see Section 4.4), estimates of their Galactic population and merger rate are of great interest as one of the prime sources for current gravitational wave detectors such as GEO600 [112], LIGO [46], VIRGO [131] and TAMA [225]. In the following, we review recent empirical determination of the population sizes and merging rates of binaries where at least one component is visible as a radio pulsar.

3.4.1 Double neutron star binaries

As discussed in Section 2.2, double neutron star (DNS) binaries are expected to be rare. This is certainly the case; despite extensive searches, only five certain DNS binaries are currently known: PSRs J0737-3039 [44Jump To The Next Citation Point], B1534+12 [343Jump To The Next Citation Point], J1756-2251 [93Jump To The Next Citation Point], B1913+16 [129Jump To The Next Citation Point] and B2127+11C [257Jump To The Next Citation Point]. Although we only see both neutron stars as pulsars in J0737-3039 [198Jump To The Next Citation Point], we are “certain” of the identification in the other four systems from precise component mass measurements from pulsar timing observations (see Section 4.4). The spin and orbital parameters, timescales and mass constraints for these systems are listed in Table 1, along with three other binaries with eccentric orbits, mass functions and periastron advance measurements that are consistent with a DNS identification, but for which there is presently not sufficient component mass information to confirm their nature.












J0737-3039 J1518+4904 B1534+12 J1756-2251





         
P [ms] 22.7/2770 40.9 37.9 28.5
Pb [d] 0.102 8.6 0.4 0.32
e 0.088 0.25 0.27 0.18
tc [108 yr] 2.1/0.5 200 2.5 4.4
t g [108 yr] 0.87 24000 27 169
Masses measured? Yes No Yes Yes










J1811-1736 J1829+2456 B1913+16 B2127+11C





         
P [ms] 104.2 41.0 59.0 30.5
Pb [d] 18.8 1.18 0.3 0.3
e 0.83 0.14 0.62 0.68
tc [108 yr] 9.2 130 1.1 0.97
tg [108 yr] 10000 600 3.0 2.2
Masses measured? No No Yes Yes











Table 1: Known and likely DNS binaries. Listed are the pulse period P, orbital period P b, orbital eccentricity e, characteristic age tc, and expected binary coalescence time-scale tg due to gravitational wave emission calculated from Equation (7View Equation). To distinguish between definite and candidate DNS systems, we also list whether the masses of both components have been determined via the measurement of two or more post Keplerian parameters as described in Section 4.4.

Update Despite the uncertainties in identifying DNS binaries, for the purposes of determining the Galactic merger rate, the systems2 for which t g is less than t u (i.e. PSRs J0737-3039, B1534+12, J1756-2251, B1913+16 and B2127+11C) are primarily of interest. Of these PSR B2127+11C is in the process of being ejected from the globular cluster M15 [257252] and is thought to make only a negligible contribution to the merger rate [248Jump To The Next Citation Point]. The general approach with the remaining systems is to derive scale factors for each object, construct the probability density function of their total population (as outlined in Section 3.2.1) and then divide these by a reasonable estimate for the lifetime. It is generally accepted [151] that the observable lifetimes for these systems are determined by the timescale on which the current orbital period is reduced by a factor of two [9]. Below this point, the orbital smearing selection effect discussed in Section 3.1.3 will render the binary undetectable by current surveys. The results of the most recent study of this kind [147Jump To The Next Citation Point148Jump To The Next Citation Point], which take into account the discovery of PSR J0737-3039, are summarised in the left panel of Figure 19View Image. The combined Galactic merger rate, dominated by the double pulsar, is found to be 80+-27100 Myr -1, where the uncertainties reflect the 95% confidence level using the techniques summarised in Section 3.2.2. Extrapolating this number to include DNS binaries detectable by LIGO in other galaxies [248], the expected event rate is 35+90 × 10-3 yr-1 -30 for initial LIGO and +470 - 1 190 -150 yr for advanced LIGO. Future prospects for detecting gravitational wave emission from binary neutron star inspirals are therefore very encouraging. Since much of the uncertainty in the rate estimates is due to our ignorance of the underlying distribution of double neutron star systems, future detection rates could ultimately constrain the properties of this exciting binary species.
View Image

Figure 19: The current best empirical estimates of the coalescence rates of relativistic binaries involving neutron stars. The individual contributions from each known binary system are shown as dashed lines, while the solid line shows the total probability density function on a logarithmic and (inset) linear scale. The left panel shows DNS binaries [147, 148Jump To The Next Citation Point], while the right panel shows the equivalent results for NS-WD binaries [157Jump To The Next Citation Point].
Although the double pulsar system J0737-3039 will not be important for ground-based detectors until its final coalescence in another 85 Myr, it may be a useful calibration source for the future space-based detector LISA [224Jump To The Next Citation Point]. Recent calculations [146Jump To The Next Citation Point] show that a 1-yr observation with LISA would detect (albeit with S/N ~ 2) the continuous emission at a frequency of 0.2 mHz based on the current orbital parameters. Although there is the prospect of using LISA to detect similar systems systems through their continuous emission, current calculations [146] suggest that significant (S/N > 5) detections are not likely. Despite these limitations, it is likely that LISA observations will be able to place independent constraints on the Galactic DNS binary population after several years of operation.

3.4.2 White dwarf-neutron star binaries

Although the population of white dwarf-neutron star (WDNS) binaries in general is substantial, the fraction which will merge due to gravitational wave emission is small. Like the DNS binaries, the observed WDNS sample suffers from small-number statistics. From Figure 18View Image, we note that only three WDNS systems are currently known that will merge within tu, PSRs J0751+1807 [192Jump To The Next Citation Point], J1757-5322 [90Jump To The Next Citation Point] and J1141-6545 [154Jump To The Next Citation Point]. Applying the same techniques as used for the DNS population, the merging rate contributions of the three systems can be calculated [157Jump To The Next Citation Point] and are shown in Figure 19View Image. The combined Galactic coalescence rate is 4+5 Myr -1 -3 (at 68% confidence interval). Although the orbital frequencies of these objects at coalescence are too low to be detected by LIGO, they do fall within the band planned for LISA [224]. Unfortunately, an extrapolation of the Galactic event rate out to distances at which such events would be detectable by LISA does not suggest that these systems will be a major source of detection [157Jump To The Next Citation Point]. Similar conclusions were reached by considering the statistics of low-mass X-ray binaries [69].


  Go to previous page Go up Go to next page