Gravitational wave background from mergers of large primordial black
  holes
  The Peters formula, which tells how the coalescence time of a binary system
emitting gravitational radiation is determined by the initial size and shape of
the elliptic orbit, is often used in estimating the merger rate of primordial
black holes and the gravitational wave background from the mergers. Valid as it
is in some interesting scenarios, such as the analysis of the LIGO-Virgo
events, the Peters formula fails to describe the coalescence time if the
orbital period of the binary exceeds the value given by the formula. This could
underestimate the event rate of mergers that occur before the cosmic time
$t\sim 10^{13}\ \text{s}$. As a result, the energy density spectrum of the
gravitational wave background could develop a peak, which is from mergers
occurring at either $t\sim 10^{13}\ \text{s}$ (for black holes with mass
$M\gtrsim 10^8 M_\odot$) or $t\sim 10^{26}(M/M_\odot)^{-5/3}\ \text{s}$ (for
$10^5 M_\odot \lesssim M\lesssim 10^8 M_\odot$). This can be used to constrain
the fraction of dark matter in primordial black holes (denoted by $f$) if
potential probes (such as SKA and U-DECIGO) do not discover such a background,
with the result $f\lesssim 10^{-6}\text{-}10^{-4}$ for the mass range
$10\text{-} 10^9M_\odot$. We then consider the effect of mass accretion onto
primordial black holes at redshift $z\sim 10$, and find that the merger rate
could drop significantly at low redshifts. The spectrum of the gravitational
wave background thus gets suppressed at the high-frequency end. This feature
might be captured by future detectors such as ET and CE for initial mass $M=
\mathcal{O}(10\text{-}100) M_\odot$ with $f\gtrsim 10^{-4}$.
