Impact of LIGO-Virgo black hole binaries on gravitational wave background searches

We study the impact of the black hole binary population currently probed by LIGO-Virgo on future searches for the primordial gravitational wave background. We estimate the foreground generated by the binaries using the observed event rate and a simple modeling of the black hole population. We subtract individually resolvable binaries from the foreground and utilize Fisher analysis to derive sensitivity curves for power-law signals including these astrophysical foregrounds. Even with optimistic assumptions, we find that the reach of future experiments will be severely reduced.

Introduction -Various cosmological processes may have generated gravitational waves (GWs) that contribute to the stochastic GW background [1][2][3][4][5].If sufficiently strong, this background can be probed with GW detectors.The non-observation of such background with LIGO-Virgo detectors [6] translates into mild constraints on the possible cosmological GW sources [7][8][9].As shown in Fig. 1, many new experimental programs, built on the success of currently running LIGO-Virgo network, are expected to enter GW search in the next decades.With tremendous sensitivity improvements, these experiments will further probe the early Universe processes that can generate GWs.
The LIGO-Virgo detectors have observed dozens of GW signals from compact object mergers [10][11][12].From these observations we now have a good understanding of the population of O(10M ) black holes (BHs) and their present merger rate [13,14].We can therefore estimate how strong astrophysical GW foreground these binary BHs (BBHs) generate.This foreground will impact the detectability of the primordial GW background with future GW experiments.We will present a simple analysis quantifying the effect of the BBH foreground on searches of the primordial GW background.
Other astrophysical GW foregrounds arise from neutron star (NS) binaries, NS-BH binaries, and from binary white dwarfs (BWDs).The latter has been extensively studied in the literature [15,16], and we include its impact on the reach of GW experiments in our analysis as modeled for the LISA experiment in [17,18].We neglect the contributions on the foreground arising from NS-NS and NS-BH binaries, as their population characteristics and merger rate is currently subject to large uncertainties and their contribution to the total GW foreground is subdominant at low frequencies to that arising from BBHs [14].However, they may still have an impact on the reach of the GW experiments as it is typically more difficult to resolve GW signals from NS-NS and NS-BH binaries than from BBHs.
Our analysis complements the earlier studies by subtracting the resolvable BBHs from the total foreground and systematically analysing the detectability of a primordial GW background. 1 We account for the uncertainties in the BBH merger rate and perform a Fisher analysis to study the impact the BBH foreground on detection of primordial GW backgrounds modeled in the simplest case as power-laws.Moreover, we develop a new graphical presentation of the sensitivity of a given detector that accounts for the foreground.We find that the sensitivity of all experiments will be severely impacted by the BBH foreground, limiting their reach in probing primor- 1 Subtracting the resolvable binaries from the foreground has been studied in the earlier literature [20,24,29,30], however without systematical analysis of the detectability of the primordial GW background.Ref. [28] presented a Bayesian implementation that estimates the astrophysical foreground and the detectability of the primordial GW background, but they demonstrate their method only in one benchmark case.
dial GW sources.
BBH merger rate -Assuming that the BH population that LIGO-Virgo is probing is astrophysical, we use a truncated power-law mass function and a merger rate that follows the star formation rate with a power-law delay time distribution distribution. 2The differential BBH merger rate then is of the form where the normalization factor Z ψ is defined such that R = R 0 at z = 0, M = m 1 + m 2 and η = m 1 m 2 /M are the total mass and the symmetric mass ratio of the binary, and the mass dependence of the binary population is parametrized by α and β.The function P b (z) is given by the star formation rate [49], and for the time delay distribution we use P d (t) ∝ t −1 at t > 50 Myr [50].The mass function is with the normalization ψ(m)d ln m = 1.We take m min = 3.0M and m max = 55M corresponding to the estimates for the maximal mass of neutron stars and the beginning of the pair instability mass gap.We use the values α = 0, β = 6, ζ = −1.5 and R 0 = 10 +6 −5 Gpc −3 yr −1 , following the fit performed in [51] to the LIGO-Virgo observations. 3For indicative errors or the resulting stochastic GW foreground, we include the uncertainties only for the present merger rate R 0 and use the central values for the other parameters.This greatly simplifies our analysis.

Detectability of individual BBHs -The amplitude
| h(f )| of the Fourier transform of the inspiral-mergerringdown GW signal from a BBH can be approximated as [52] where D L is the luminosity distance of the BBH, and the frequencies f merg , f ring , f cut and σ are of the form with the coefficients a j , b j and c j given in Table I of Ref. [52].
The signal-to-noise ratio, SNR(f ), of a BBH signal depends on the orientation of the detector with respect to the sky location of the binary and on the inclination of the binary.For an optimally oriented binarydetector system, neglecting the rotation of the detector, the signal-to-noise ratio of the BBH signal that in the end of the experiment is at frequency f , or the BBH has merged, is given by where S n is the noise power spectrum of the GW detector, that includes the instrumental noise and the GW foreground.The lower integration limit depends on the duration of the experiment, T , and the coalescence time of the binary, We assume a duration of T = 4 yr for all experiments.We estimate how large fraction of the BBH signals can be resolved by calculating the probability [53,54] where r(f ) = SNR c /SNR(f ), that accounts for the antenna patterns F +,× of the detector and averages over the binary sky location and inclination, and the polarization of the signal.The parameter ω is defined as and its probability distribution p(ω) is calculated assuming uniform distributions for the binary inclination c i ∈ (−1, 1), sky location cos θ ∈ (−1, 1) and φ ∈ (0, 2π), and the polarization angle ψ ∈ (0, 2π).We use SNR c = 8 as the threshold signal-to-noise ratio.
GW foreground from LIGO-Virgo binaries -To calculate the GW foreground generated by the unresolvable BBH, we subtract from the total GW energy density emitted by all binaries the contribution from those that are sufficiently loud to be seen individually.The GW foreground then is different for different detectors, and its dimensionless energy density is (see e.g.[55,56]) where dV c is the differential comoving volume element at redshift z, ρ c is the critical energy density of the Universe, dρ GW is the GW energy density emitted by a binary in the frequency range (f, f + df ) [56], dR is given by Eq. ( 1), and the factor 1−p det (f ) removes the binaries that are individually resolvable within the duration of the given experiment.Note that, in the analysis of actual data the estimated gravitational waveform for every detected event is subtracted from the time series data and the subtraction is affected by the errors in the parameter estimation.By comparing to the results obtained in [29] through the latter method, we find that our analysis is conservative regarding to how much the subtraction reduces the noise.In Fig. 2 we show the GW foreground from unresolvable BBHs for different GW experiments 6 .The solid curves in Fig. 2 indicate the foreground for the central value of R 0 .Since Ω BBH ∝ R 0 the boundaries of the 1σ band on the BBH foreground can be obtained by multiplying the central value result by 1.6 and 0.5 for the upper and lower limit.The gray dotted curve in Fig. 2 indicates the total GW foreground before any binaries are resolved, obtained by taking p det → 0 in Eq. (10).At f 100 Hz this curve is simply a powerlaw, Ω BBH ≈ 3.2 × 10 −11 (f /Hz) 2/3 .In Fig. 1 the BBH strain curve corresponds to the total BBH foreground for the central value of R 0 .This is what we use as the BBH contribution to the instantaneous noise.
Solid curves show the central value of the GW foreground from unresolvable BH binaries for different observatories, while the bands indicate their uncertainties stemming from uncertainty on BBH merger rate.Dashed curves show the reach of the corresponding observatories in terms of PI sensitivities, and the gray dotted curve indicates the total GW energy density emitted by the BH population.We assume T = 4 yr.
The dashed curves in Fig. 2 show the standard powerlaw integrated (PI) sensitivities [65] of the experiments, which account for the foregrounds only in the instantaneous detector noise.The signal-to-noise ratio of the GW background is given by [66] and the PI curves are obtained as the envelope of the power-laws that give SNR BG = SNR c = 8.The detector noise is the sum of the instrumental noise, and the GW foregrounds from BWDs (which we include as in Refs.[17,18]) and BBHs, Ω n (f ) = Ω instr (f )+Ω BWD (f )+ Ω BBH (f ).Notice that here Ω BBH (f ) is the total BBH foreground before any binaries are resolved.We note that the GW foreground generated by BH binaries is not exactly smooth but includes fluctuations [67,68].At low frequencies, f 0.01 Hz, the foreground from solar mass BBHs consists of a large number of nearly-monochromatic signals and we expect the fluctuations to follow a narrow Gaussian distribution.At higher frequencies the fluctuations become larger and constitute so called 'popcorn' signal as the number of binaries contributing to the signal in a given frequency bin decreases.These fluctuations may help distinguishing this foreground from the primordial GW background.In this work we use the average GW foreground (10) and leave these improvements for future work.

Detectability of a primordial GW background -
To study the impact of the BBH foreground on the detectability of the primordial GW background, we estimate how accurately the background can be measured using Fisher analysis [69] (see also e.g.[20,22,29,30,32,33,[70][71][72][73] for applications to GW data analysis).The total stochastic noise in the GW detector is given by the sum of the primordial GW background, which we model as a power-law, the WDB and the BBH foregrounds, and the instrumental noise, where Ω is the amplitude of the primordial GW background at a reference frequency f ref , α is the spectral index of the primordial GW background, and A Ω BBH (f ) is the contribution from unresolvable BBHs.We calculate Ω BBH (f ) using the central value for the present merger rate, R 0 = 10 Gpc −3 yr −1 , and we have introduced the parameter A = R 0 / R 0 to account for the uncertainties in the amplitude of the BBH foreground.
As our main aim is to demonstrate the impact of the BBH foreground on the detectability of a primordial GW background, for simplicity we don't include the uncertainties in the BWD foreground or in the instrumental noise in our analysis. 7Moreover, simultaneously utilizing data from multiple experiments would in principle improve the prospects of detecting cosmological GW backgrounds [89].However, in practice the relevant signals can be modeled as power-laws only locally [1][2][3][4][5] and in particular signals whose detectability in a given experiment depends on the impact of the foreground are not guaranteed to be visible in other experiments.For this reason we consider each experiment individually.
The variances of the model parameters can be estimated from the inverse of the Fisher matrix, whose ij component is given by as Here the indices i and j label the parameters of the model, i, j ∈ {Ω, α, A}.Keeping α and A fixed (that is, calculating the Fisher matrix only in Ω) we get which shows that the amplitude of any power-law that touches the PI sensitivity curve can be measured with 7 The uncertainties in the BWD foreground and in the instrumental noise, and their effects are considered previously in [27,[74][75][76][77][78][79][80][81][82] and [83][84][85][86][87][88].
3. Solid and dot-dashed curves show the sensitivities on primordial GW background accounting for the BBH foreground with the prior σA,0 = 0.6 and σA,0 = 0.06, respectively.Dashed curves show the PI sensitivities as in Fig. 2.
accuracy σ Ω /Ω = SNR −1 c , assuming that the spectral index α and the BBH foreground are known with negligible uncertainties.
Next we aim to find a curve similar to the PI one, which takes into account the uncertainties in the BBH foreground.We perform a Fisher analysis in Ω and A for different fixed values of α which we assume to know exactly as in the PI prescription.We account for the measured uncertainties in the present BBH merger rate R 0 by setting a prior on A by the replacement We use σ A,0 = 0.6, which arises from the 1σ upper bound R 0 < 16 Gpc −3 yr −1 .The solid curves in Fig. 3 show the sensitivities of different experiments obtained as described above by requiring that σ Ω /Ω < SNR −1 c .In the same way as the PI curve, the curves shown in Fig. 3 depict the envelope of different power-laws that are detectable.As expected, these curves lie above the PI curves.The effect is particularly strong for LIGO, LISA and AEDGE, whereas ET and BBO can probe the primordial GW background well below the BBH foreground.The reason is that ET and BBO will be able to measure the BBH foreground very accurately even without an accurate prior knowledge of the BBH merger rate.Therefore also changing the prior of A does not significantly change the sensitivity prospects for these experiments.Instead, for LISA and LIGO more accurate measurement of the BBH merger rate will greatly improve the prospects for probing primordial GW background.We show this by the dot-dashed curves that are calculated assuming an order of magnitude improvement in the accuracy of the BBH merger rate, which would imply σ A,0 = 0.06.For LIGO this curve overlaps with the PI sensitivity, whereas Minimal detectable amplitude of the power-law primordial GW background as a function of its spectral index α.The solid and dot-dashed curves include the BBH foreground with the current and an order of magnitude improved accuracy of the BBH merger rate, respectively.The dashed curves correspond to the usual PI sensitivities.We assume T = 4 yr.
for BBO and ET it overlaps the curve obtained with the current accuracy on the rate, σ A,0 = 0.6.Fig. 3, however, does not give a complete picture of the detectability of the primordial GW background.While any power-law that touches the sensitivity curves shown in Fig. 3 can be detected, for example with LISA assuming σ A = 0.6 a flat primordial GW background can be detected even if it doesn't touch the solid green curve.As the foreground for LISA is a power-law with spectral index 2/3, detecting a primordial GW background with α = 2/3 requires a larger amplitude, and therefore that power-law cuts the sensitivity curve for LISA in Fig. 3.We show the minimal detectable amplitude for a given value of α at the peak sensitivity of each experiment in Fig. 4. For example, we see that for the primordial GW background with α = 2/3 to be detectable by LISA, its amplitude needs to be higher than 3 × 10 −12 at 3 mHz, assuming the current uncertainties in the BBH merger rate.
Finally, for an easy comparison of different cases, in Table I we show the values of minimal detectable amplitude of a flat, α = 0, primordial GW background, including also the case that the uncertainties in α are accounted for in the Fisher analysis.The latter case denoted with σ α,0 → ∞, and the case where α is fixed is denoted by σ α,0 → 0. When α is included in the Fisher analysis, the result depends also on the reference frequency f ref .In Table I the reference frequencies for each experiment are fixed in the same way as in Fig. 4. In all cases the minimal detectable amplitude of the primordial GW background is determined by requiring that σ Ω /Ω < SNR −1 c .For the first row in the table we assume that the foreground is known and calculate the Fisher matrix only in Ω.As described above, this corresponds to the usual PI result.
Conclusions -We have analysed the impact of the BH population currently probed by LIGO-Virgo on future searches for a primordial GW background.The GW foreground associated with the current observations has not been probed directly yet.However, the rate of events inferred by the LIGO-Virgo observations suggests that it will have a large amplitude, and will cut significantly into the integrated sensitivities of upcoming experiments, limiting their reach.
We have estimated the BBH foreground by subtracting from the total GW energy density emitted by BBHs the contribution from the binaries that can be resolved individually.Then, using Fisher analysis, we have assessed the impact of this foreground on future searches of a primordial GW background.Describing the primordial background as a power-law with a given slope, we have checked what amplitudes could be probed in the presence of the foreground with each experiment.We have found the impact to be much smaller than one would expect from simply looking at the amplitude of the foreground, especially for detectors that can themselves probe the foreground accurately.Instead for detectors unable to significantly improve the foreground measurement the final result depends on priors.We have illustrated this by showing also results assuming an order of magnitude improvement in the accuracy of the BBH merger rate.In the case of LISA and LIGO this significantly alleviates the problem, while giving almost negligible change in AEDGE, BBO and ET.However, while subtraction of the loud identifiable binaries and improvements in the measurement of the BBH merger rate help, the sensitivity of all future experiments in terms of the abundance of the primordial background they can reach is lowered in each case.

5 ×
10 −9TABLE I. Minimal detectable amplitude Ω of a flat (α = 0) primordial GW background.The first two columns show the priors used for the relative amplitude A of the BBH foreground and the spectral index α of the primordial GW background.