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 similar to 10(13) 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 similar to 10(13) s (for black holes with mass M >= 10(8) M-circle dot) or t similar to 10(26)(M/M-circle dot)(-5/3) s (for 10(5)M(circle dot) less than or similar to M-circle dot 10(8)M(circle dot)). This can be used to constrain the fraction of dark matter in primordial black holes (denoted by integral) if potential probes (such as SKA and U-DECIGO) do not discover such a background, with the result f less than or similar to 10(-6)-10(-4) for the mass range 10-10(9) M-circle dot. We then consider the effect of mass accretion onto primordial black holes at redshift z similar to 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 = O(10-100) M-circle dot with f greater than or similar to 10(-4).