A new method for estimates of binary pulsar coalescence rates

Estimates of the Galactic coalescence rate (R) of close binaries with two neutron stars (NS-NS) are known to be uncertain by large factors (about two orders of magnitude) mainly due to the small number of systems detected as binary radio pulsars. We present an analysis method that allows us to estimate the Galactic NS-NS coalescence using the current observed sample and, importantly, to assign a statistical significance to these estimates and to calculate the allowed ranges of values at various confidence levels. The method involves the simulation of selection effects inherent in all relevant radio pulsar surveys and a Bayesian statistical analysis for the probability distribution of R. The most likely values for the total Galactic coalescence rate (R_tot) lie in the range 2-60 Myr^-1 depending on different pulsar population models. For our reference model 1, where the most likely estimates of pulsar population properties are adopted, we obtain R_tot = 8_-5⁺⁹ Myr^-1 at a 68% statistical confidence level. The corresponding range of expected detections rates of NS-NS inspiral are 3_-2⁺⁴ × 10^-3 yr^-1 for the initial LIGO and 18_-11⁺²¹ yr_-1 for the advanced LIGO.