Constructing a gravitational wave analysis pipeline for extremely large mass ratio inspirals

TL;DR

XMRI signals from a brown dwarf orbiting Sgr A* offer a unique probe of strong-field gravity in the millihertz band. The authors introduce a three-stage hierarchical semi-coherent search that combines a multi-harmonic F-statistic, particle swarm optimization, and a novel multi-harmonic eccentricity estimator, tailored to TianQin. Validation with simulated data shows precise recovery of the orbital frequency, eccentricity, and black-hole parameters, with 90-day observations yielding sub-percent to per-mille level accuracy and a BD mass estimate at the $2\times10^{-3}$ level. This framework provides a practical, scalable path to XMRI discovery and high-fidelity tests of the Kerr spacetime with future space-based GW detectors.

Abstract

Extremely large mass-ratio inspirals (XMRIs), consisting of a brown dwarf orbiting a supermassive black hole, emit long-lived and nearly monochromatic gravitational waves in the millihertz band and constitute a promising probe of strong-field gravity and black-hole properties. However, dedicated data-analysis pipelines for XMRI signals have not yet been established. In this work, we develop, for the first time, a hierarchical semi-coherent search pipeline for XMRIs tailored to space-based gravitational-wave detectors, with a particular focus on the TianQin mission. The pipeline combines a semi-coherent multi-harmonic $\mathcal{F}$-statistic with particle swarm optimization, and incorporates a novel eccentricity estimation method based on the relative power distribution among harmonics. We validate the performance of the pipeline using simulated TianQin data for a Galactic center XMRI composed of a brown dwarf and Sgr A*. For a three-month observation, the pipeline successfully recovers the signal and achieves high-precision parameter estimation, including fractional uncertainties of $<10^{-6}$ in the orbital frequency, $\lesssim10^{-3}$ in the eccentricity, $\lesssim2\times10^{-3}$ in the black-hole mass, and $\lesssim10^{-3}$ in the black-hole spin. Our framework establishes a practical foundation for future XMRI searches with space-based detectors and highlights the potential of XMRIs as precision probes of stellar dynamics and strong-field gravity in the vicinity of supermassive black holes.

Figures (9)

• Figure 1: Geometry of an eccentric XMRI system in the black-hole spin–aligned frame $(X,Y,Z)$ and the orbital frame $(X',Y',Z')$. $Z$-axis is parallel to the spin of Sgr A* ($\vec{S}$), while $Z'$-axis aligned with the orbital angular momentum $\vec{L}$, where the $X'$-axis points towards the pericenter $P$. The misalignment between $\vec{L}$ and $\vec{S}$ is characterized by the inclination angle $\iota$. The unit vector $\hat{k}$ denotes the propagation direction of the GW toward Earth, with $(\theta, \phi)$ being its corresponding polar and azimuthal angles in the orbital frame, and $\delta$ being its angle relative to $\vec{S}$. The geometry is further specified by the longitude of the ascending node (point $A$) $\alpha$ and the argument of pericenter $\gamma$.
• Figure 2: The detection statistic $2\mathcal{F}_n$ versus harmonic index $n$. Noiseless results (red dots) coincide with the optimal SNR squared $\rho_n^2$ (blue open circles), while noisy realizations are shown as green crosses.
• Figure 3: Comparison between the GWSpace template and the analytic template, obtained via direct Fourier transform of the time-domain signals. Upper panel: The strain of the $n=3$ harmonic for the GWSpace template (blue), the analytic template (grey), and their resulting strain residual (red). The vertical green dashed line marks the frequency $3f_0$, where $f_0$ denotes the orbital frequency of the XMRI system, while the frequency sidebands surrounding the peak result from the combined effects of the XMRI orbital precession and the periodic motion of TianQin. Lower panel: The phase evolution of the two templates and their corresponding phase residual (red). The close agreement between the templates, especially near the peak frequency, demonstrates the reliability of the analytic approximation used in our search pipeline.
• Figure 4: Schematic of the hierarchical semi-coherent search pipeline for XMRI signals. The top panels illustrate the foundational components: (left) the generation of the $\mathcal{F}$-statistic using downsampled data and analytic TDI templates; (center) the evolution of the $2\mathcal{F}$ peak as the coherence time $T_c$ increases from 15 to 90 days. The broader peaks at shorter $T_c$ provide a larger capture range that facilitates the initial global search, while the narrower peaks at longer $T_c$ yield significantly higher parameter precision, motivating our hierarchical semi-coherent approach; and (right) the use of the multi-harmonic power distribution $2\mathcal{F}_n$ to provide an initial probe of the initial eccentricity $e_0$. The bottom panel details the three-stage hierarchical PSO workflow. Stage I ($T_c = 15$ d) focuses on the initial localization of orbital frequency $f_0$ and initial eccentricity $e_0$. Stage II ($T_c = 30$ d) performs an intermediate refinement to break parameter degeneracies among the intrinsic parameters. Stage III ($T_c = 90$ d) executes the final fully coherent search to obtain high-precision estimates for the complete parameter set, including the analytically inferred extrinsic parameters. At each transition, the parameter space is refined based on a $t$-distribution analysis of the best-performing PSO runs.
• Figure 5: Recovered values of the orbital inclination angle $\iota$ (blue circles) and Sgr A* spin $s$ (red squares) from ten independent PSO runs for Stage I (open) and Stage II (filled). The results obtained in Stage II exhibit a substantially reduced scatter and cluster more closely around the injected parameters (vertical dashed lines), although a residual spread remains.