Formation and evolution of primordial black hole binaries in the early universe

The abundance of primordial black holes (PBHs) in the mass range 0.1-10(3) M-Theta can potentially be tested by gravitational wave observations due to the large merger rate of PBH binaries formed in the early universe. To put the estimates of the latter on a fi rmer footing, we fi rst derive analytical PBH merger rate for general PBH mass functions while imposing a minimal initial comoving distance between the binary and the PBH nearest to it, in order to pick only initial con fi gurations where the binary would not get disrupted. We then study the formation and evolution of PBH binaries before recombination by performing N-body simulations. We fi nd that the analytical estimate based on the tidally perturbed 2-body system strongly overestimates the present merger rate when PBHs comprise all dark matter, as most initial binaries are disrupted by the surrounding PBHs. This is mostly due to the formation of compact N-body systems at matter-radiation equality. However, if PBHs make up a small fraction of the dark matter, f PBH less than or similar to 10%, these estimates become more reliable. In that case, the merger rate observed by LIGO imposes the strongest constraint on the PBH abundance in the mass range 2-160M(Theta). Finally, we argue that, even if most initial PBH binaries are perturbed, the present BH-BH merger rate of binaries formed in the early universe is larger than O (10) Gpc(-3) yr(-1) f(3)PBH.