The Probability Distribution of Binary Pulsar Coalescence Rates. I. Double Neutron Star Systems in the Galactic Field

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ₚeak 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ₜot = 8-₅^+9 Myr^-1 at a 68% statistical confidence level. The corresponding range of expected detection rates of NS--NS inspiral are 3-₂^+4x 10^-3 yr^-1 for the initial LIGO and 18-₁₁^+21 yr^-1 for the advanced LIGO.