Constraining the neutron star-black hole merger rate

By incorporating spin-induced orbital precession into gravitational-wave searches for the first time, this study significantly improves detection sensitivity for neutron star-black hole binaries, leading to a 16% downward revision of their estimated merger rate and the identification of four new subthreshold candidates likely of terrestrial origin.

The Gist

Imagine the universe is a giant, dark ocean, and gravitational waves are the ripples created when massive objects, like neutron stars and black holes, crash into each other. For years, scientists have been trying to "hear" these ripples using giant ears called detectors (LIGO, Virgo, and KAGRA).

To find these ripples, they use a method called "matched filtering." Think of this like trying to find a specific song in a noisy room. You have a playlist of known songs (called templates), and you compare the noise in the room against your playlist to see if a match pops up.

The Problem: The Playlist Was Missing a Key Feature
Until now, the scientists' playlist had a major blind spot. They assumed that when a black hole and a neutron star dance together, they spin perfectly in sync, like a figure skater spinning straight up.

However, in reality, these cosmic dancers often wobble. If the black hole is spinning at a weird angle compared to the orbit, the whole system starts to precess (wobble like a spinning top that's about to fall over). The old playlist didn't include songs with this "wobble." So, if a wobbly pair crashed, the scientists' ears might have missed it entirely because the sound didn't match their rigid playlist.

The authors of this paper realized that because black holes and neutron stars have very different masses, this wobble is actually quite common and creates a very distinct "sound." By ignoring it, they might have been missing up to 85% of the crashes happening in our local universe.

The Solution: A New, Smarter Playlist
The researchers created a new search method that includes these "wobbly" signals in their playlist for the first time. They tested this new method on data from the third major observing run of the gravitational wave detectors.

Here is what they found:

• Super Sensitive Ears: For systems that wobble strongly, their new method is up to 100% more sensitive than the old way. It's like upgrading from a tin can telephone to a high-tech radar; they can hear the same signal from twice as far away.
• Fewer Crashes Than We Thought: Because they can now hear these signals from much further away, they realized that the "volume" of space they are listening to is much bigger than before. When you listen to a bigger volume of space and still only hear a few crashes, it means the actual rate of crashes in the universe is likely lower than we previously calculated. Specifically, they found the overall rate of these mergers is about 16% smaller than earlier estimates.

The "Wobbly" Sub-Group
They also looked specifically at the "wobbly" (precessing) pairs. Even with their super-sensitive new ears, they didn't find any confirmed wobbly crashes in the data. This allows them to set a strict limit: there are likely no more than 79 of these specific wobbly crashes happening per cubic billion light-years every year.

The "Almost" Finds
The new search also picked up four "marginal" candidates—signals that were just barely too quiet to be confirmed as real crashes. Interestingly, all four of these faint signals showed strong signs of wobbling. However, the scientists are cautious: they believe these are likely just "static" or noise from Earth (terrestrial origin) rather than real cosmic events, so they didn't count them in their final numbers.

Why This Matters
By fixing the "playlist" to include wobbling spins, the scientists aren't just finding more signals; they are getting a more accurate picture of how often these cosmic collisions happen. This helps us understand how these pairs form in the first place—whether they are born from stars that evolved together peacefully (which usually don't wobble much) or from chaotic crowds of stars in dense clusters (which often do wobble).

In short: They built a better hearing aid, realized the universe is actually quieter than we thought, and learned that the cosmic dancers might be wobbling more than we expected, even if we haven't heard them crash yet.

Technical Summary

Technical Summary: Constraining the Neutron Star-Black Hole Merger Rate

Problem Statement
Current template-based gravitational-wave (GW) searches for compact binary mergers, specifically those conducted by the LIGO–Virgo–KAGRA (LVK) collaboration, rely on matched filtering against pre-generated waveform banks. A critical limitation of these existing analyses is the neglect of general relativistic spin-induced orbital precession. Templates are typically restricted to scenarios where component spins are aligned or anti-aligned with the orbital angular momentum.

This restriction is problematic for neutron star-black hole (NSBH) binaries due to their asymmetric masses, which make them prime candidates for exhibiting strong spin-precession imprints. If the astrophysical population includes binaries with significant spin misalignment (high effective perpendicular spin, $\chi_p$), current search pipelines may fail to detect a substantial fraction of these events. Consequently, the detected NSBH population could be suppressed or biased, leading to inaccurate inferences regarding the NSBH merger rate and the underlying formation channels (e.g., isolated binary evolution vs. dynamical capture in dense stellar environments).

Methodology
The authors present the first dedicated template-based matched filter search for NSBH binaries that explicitly incorporates spin-precession. The study utilizes data from the entirety of the third LIGO–Virgo–KAGRA observing run (O3).

• Search Technique: The authors apply a novel method (previously proposed in Ref. [22]) to incorporate misaligned spins into the search. Unlike standard LVK searches that reduce the parameter space dimensionality by assuming spin alignment (reducing the problem to 4 dimensions: two masses and two spin magnitudes), this "generic-spin" search accounts for the full 15-dimensional parameter space required to describe NSBH mergers on circular orbits without neglecting precession terms.
• Sensitivity Estimation: Search sensitivity is estimated by injecting simulated signals drawn from a reference population model into real GW detector noise and calculating the recovery rate. The authors assume an agnostic population model due to large uncertainties in the underlying NSBH distribution.
• Rate Calculation: The merger rate density is calculated by combining the estimated search sensitivity with the number of confident GW signals detected. The analysis employs a Poisson likelihood over the astrophysical rate, a Jeffreys prior, and marginalizes over a 15% uncertainty arising from detector calibration.
• Comparison: The authors generated a reduced version of the standard LVK search (restricting spins to be aligned) covering the same parameter space to serve as a baseline for comparison. All searches in this study were restricted to coincident signals in the two most sensitive detectors, LIGO-Hanford and LIGO-Livingston.

Key Results

1. Increased Sensitivity: The inclusion of spin-precession significantly enhances search sensitivity. On average, the generic-spin search is 44% more sensitive than current LVK techniques. For systems demonstrating strong precessional effects (effective perpendicular spin $> 0.75$) and high total mass ($> 14 M_\odot$), the sensitive volume increases by up to 100%.
2. Detection Statistics: The search identified only one confident detection ($FAR < 10^{-2} \text{ yr}^{-1}$): GW200115 042309. Bayesian inference confirmed this event has minimal precession ($\chi_p \lesssim 0.4$). The search did not recover GW200105 162426, as it was a single-detector candidate when Virgo data was excluded.
3. Merger Rate Constraints:
   • Precessing Subpopulation: For a theoretical subpopulation of NSBH binaries with $\chi_p > 0.4$, the authors place an upper limit on the merger rate of $R_{90} = 79 \text{ Gpc}^{-3}\text{yr}^{-1}$ (at 90% confidence). This represents a 40% reduction in the upper limit compared to searches using only aligned-spin waveforms.
   • Overall Rate: When accounting for the full parameter space, the inferred overall NSBH merger rate in the local Universe is $44^{+102}_{-37} \text{ Gpc}^{-3}\text{yr}^{-1}$. This is, on average, 16% smaller than rates estimated using search techniques that neglect spin-precession (which yielded $52^{+130}_{-44} \text{ Gpc}^{-3}\text{yr}^{-1}$).
4. Population Fraction: By comparing the rate of precessing binaries (assuming zero observations) with the expected total rate, the authors estimate that no more than 30% of NSBH binaries in the Universe undergo spin-precession for a subpopulation with $\chi_p > 0.4$. This estimate is approximately 40% lower than what would be inferred from an aligned-spin search.
5. Subthreshold Candidates: The search identified four new subthreshold candidates (GW190421 203205, GW190916 134302, GW191225 091958, GW200310 030555). All exhibit strong imprints of spin precession but have low probabilities of astrophysical origin ($p_{astro} < 1\%$), suggesting they are likely of terrestrial origin and were excluded from the merger rate calculation.

Significance and Claims
The paper claims that by accounting for spin-precession, the search for NSBH binaries becomes significantly more sensitive, particularly for systems formed through dynamical channels which are expected to have isotropic spin distributions. The primary significance of this work is the ability to place tighter, more accurate constraints on the NSBH merger rate in the local Universe.

The authors assert that previous rate estimates were likely overestimated by approximately 16% because they failed to account for the increased sensitive volume gained by including precession. Furthermore, the improved sensitivity allows for a more rigorous test of NSBH formation mechanisms; the lack of observed highly precessing candidates allows the authors to place stringent upper limits on the fraction of the NSBH population that undergoes spin-precession. This capability is crucial for distinguishing between isolated binary evolution and dynamical formation channels. The paper notes that these results are indicative and do not account for uncertainties in the rate estimates, but they demonstrate that neglecting spin-precession introduces a bias in population inference.