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

C. Kim,V. Kalogera,D. R. Lorimer

DOI: https://doi.org/10.1086/345740

2002-10-22

Abstract: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 rate 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 the rate. The most likely values for the total Galactic coalescence rate (R_peak) lie in the range 2-60 per Myr 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}^{+9} per Myr at a 68% statistical confidence level. The corresponding range of expected detection rates of NS-NS inspiral are 3_{-2}^{+4}x10^{-3} per yr for the initial LIGO and 18_{-11}^{+21} per yr for the advanced LIGO.

Astrophysics,General Relativity and Quantum Cosmology

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to estimate the merger rate ($R$) of double - neutron - star (NS - NS) systems in the Milky Way, assign statistical significance to these estimates, and calculate the range of allowed values at different confidence levels. Due to the small number of detected binary radio pulsars, there is currently a large uncertainty (about two orders of magnitude) in the estimate of the merger rate. To improve this estimate, the authors propose an analytical method by simulating selection effects in all relevant radio pulsar surveys and using Bayesian statistical analysis to derive the probability distribution of the merger rate.

### Core problems of the paper

1. **Estimating the merger rate of double - neutron - star systems in the Milky Way**:
- Traditional estimation methods rely on observational data, but due to the small sample size, there is a large uncertainty.
- This paper proposes a new method, combining Monte Carlo simulation and Bayesian statistical analysis, to estimate the merger rate more accurately.

2. **Assigning statistical significance**:
- Through Bayesian statistical analysis, calculate the probability distribution of the merger rate and determine the value ranges at different confidence levels.

3. **Considering selection effects**:
- Simulate the selection effects in radio pulsar surveys to ensure more reliable estimation results.

### Method overview

- **Model construction**: Construct models of the spatial distribution and luminosity function of pulsars in the Milky Way.
- **Monte Carlo simulation**: Generate a large number of simulated pulsar samples, consider different total numbers ($N_{tot}$), and simulate observational selection effects.
- **Bayesian statistical analysis**: Calculate the probability distribution of the merger rate according to the observational data, and determine the most likely merger rate and its confidence interval.

### Results

For the reference model (assuming the most likely pulsar population properties), the most likely value of the total merger rate ($R_{tot}$) in the Milky Way is $8^{+9}_{-5}\, \text{Myr}^{-1}$ at the 68% confidence level. The corresponding LIGO detection rates are respectively:
- Initial LIGO: $3^{+4}_{-2}\times 10^{-3}\, \text{yr}^{-1}$
- Advanced LIGO: $18^{+21}_{-11}\, \text{yr}^{-1}$

These results provide an important reference for understanding the origin, evolution of double - neutron - star systems and their importance as gravitational - wave sources.