The merger rate of primordial-black-hole binaries

Primordial black holes (PBHs) have long been a candidate for the elusive dark matter (DM), and remain poorly constrained in the ~20-100 Msun mass range. PBH binaries were recently suggested as the possible source of LIGO's first detections. In this paper, we thoroughly revisit existing estimates of the merger rate of PBH binaries. We compute the probability distribution of orbital parameters for PBH binaries formed in the early Universe, accounting for tidal torquing by all other PBHs, as well as standard large-scale adiabatic perturbations. We then check whether the orbital parameters of PBH binaries formed in the early Universe can be significantly affected between formation and merger. Our analytic estimates indicate that the tidal field of halos and interactions with other PBHs, as well as dynamical friction by unbound standard DM particles, do not do significant work on nor torque PBH binaries. We estimate the torque due to baryon accretion to be much weaker than previous calculations, albeit possibly large enough to significantly affect the eccentricity of typical PBH binaries. We also revisit the PBH-binary merger rate resulting from gravitational capture in present-day halos, accounting for Poisson fluctuations. If binaries formed in the early Universe survive to the present time, as suggested by our analytic estimates, they dominate the total PBH merger rate. Moreover, this merger rate would be orders of magnitude larger than LIGO's current upper limits if PBHs make a significant fraction of the dark matter. As a consequence, LIGO would constrain ~10-300 Msun PBHs to constitute no more than ~1% of the dark matter. To make this conclusion fully robust, though, numerical study of several complex astrophysical processes - such as the formation of the first PBH halos and how they may affect PBH binaries, as well as the accretion of gas onto an extremely eccentric binary - is needed.