Hierarchical search method for gravitational waves from stellar-mass binary black holes in noisy space-based detector data

Yao Fu,Yan Wang,Soumya D. Mohanty

2024-07-15

Abstract:Future space-based laser interferometric detectors, such as LISA, will be able to detect gravitational waves (GWs) generated during the inspiral phase of stellar-mass binary black holes (SmBBHs). The detection and characterization of GWs from SmBBHs poses a formidable data analysis challenge, arising from the large number of wave cycles that make the search extremely sensitive to mismatches in signal and template parameters in a likelihood-based approach. This makes the search for the maximum of the likelihood function over the signal parameter space an extremely difficult task. We present a data analysis method that addresses this problem using both algorithmic innovations and hardware acceleration driven by GPUs. The method follows a hierarchical approach in which a semi-coherent $\mathcal{F}$-statistic is computed with different numbers of frequency domain partitions at different stages, with multiple particle swarm optimization (PSO) runs used in each stage for global optimization. An important step in the method is the judicious partitioning of the parameter space at each stage to improve the convergence probability of PSO and avoid premature convergence to noise-induced secondary maxima. The hierarchy of stages confines the semi-coherent searches to progressively smaller parameter ranges, with the final stage performing a search for the global maximum of the fully-coherent $\mathcal{F}$-statistic. We test our method on 2.5 years of a single LISA TDI combination and find that for an injected SmBBH signal with a SNR between $\approx 11$ and $\approx 14$, the method can estimate (i) the chirp mass with a relative error of $\lesssim 0.01\%$, (ii) the time of coalescence within $\approx 100$ sec, (iii) the sky location within $\approx 0.2$ ${\rm deg}^2$, and (iv) orbital eccentricity at a fiducial signal frequency of 10 mHz with a relative error of $\lesssim 1\%$. (abr.)

Instrumentation and Methods for Astrophysics,High Energy Astrophysical Phenomena,General Relativity and Quantum Cosmology

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is: in the noise data of space - baseline gravitational - wave detectors (such as LISA), detect and estimate the gravitational - wave signals from stellar - mass binary black holes (SmBBHs). Specifically, since the SmBBH signal persists for several years throughout the entire sensitive frequency band and contains more than \(10^5\) wave cycles, this makes the search method based on matched filtering extremely sensitive to small mismatches between the signal and template parameters. Therefore, the traditional grid - based global optimization method becomes computationally infeasible.

To solve this problem, the author proposes a hierarchical search method, which combines algorithmic innovation and hardware acceleration driven by Graphics Processing Units (GPUs). The main features of this method are as follows:

1. **Hierarchical search strategy**:
- Gradually narrow the search range of the parameter space in a phased manner.
- Use different numbers of frequency - domain partitions to calculate the semi - coherent F - statistic at each stage, and perform multiple Particle Swarm Optimization (PSO) operations at each stage for global optimization.

2. **Parameter space partitioning**:
- Reasonably partition the parameter space at each stage to increase the probability of PSO converging to the global maximum and avoid premature convergence to secondary maxima caused by noise.

3. **Multi - stage PSO runs**:
- Perform multiple PSO runs at each stage and in each partition to further improve performance.

4. **Fully coherent search in the final stage**:
- In the last stage, perform a search for the global maximum of the fully coherent F - statistic to obtain more accurate results.

Through this method, the author can effectively search for SmBBH signals in a noisy environment and estimate their parameters. Specifically, for the injected SmBBH signals (with a signal - to - noise ratio SNR between approximately 11 and 14), this method can:
- Keep the relative error of the chirp mass within 0.01%,
- Estimate the time of coalescence within approximately 100 seconds,
- Estimate the sky location within approximately 0.2 square degrees,
- Keep the relative error of the orbital eccentricity at a reference signal frequency of 10 mHz within 1%.

In summary, the method proposed in this paper significantly improves the ability to detect and estimate SmBBH gravitational - wave signals in noisy data, thus promoting the research progress in this field.