Black-hole - neutron-star mergers: new numerical-relativity simulations and multipolar effective-one-body model with spin precession and eccentricity

This paper presents 52 new numerical-relativity simulations of black-hole-neutron-star mergers focusing on tidal disruption, which are used to develop and validate the improved, multipolar TEOBResumS-Dalí model capable of accurately describing precessing, eccentric systems with enhanced accuracy at merger.

The Gist

Imagine the universe as a giant, silent ocean. When two massive objects, like a black hole and a neutron star, dance toward each other, they create ripples in the fabric of space and time called gravitational waves. Detecting these ripples is like trying to hear a whisper in a hurricane; the signal is incredibly faint and complex.

This paper is about building a better "ear" to hear that whisper, specifically when a black hole swallows a neutron star. Here is the story of what the authors did, explained simply:

1. The Problem: The "Whisper" is Hard to Decode

For a long time, scientists have been great at predicting the sound of two black holes merging (like two heavy bowling balls colliding). But when a black hole meets a neutron star (a city-sized ball of super-dense matter), the physics gets messy. The black hole's gravity can stretch and tear the neutron star apart before it gets swallowed, creating a splash of matter and a different kind of "sound."

Current models were like a blurry photograph of this event. They didn't capture the details of the "splash" (tidal disruption) well enough to tell us exactly what happened.

2. The Solution: Running 52 "Cosmic Movies"

To fix this, the authors ran 52 new, high-definition computer simulations. Think of these as running 52 different movies of black holes eating neutron stars, changing the ingredients slightly each time:

• The Recipe: They used different types of "neutron star dough" (Equations of State) to see how stiff or squishy the star was.
• The Spin: They changed how fast the black hole was spinning and in which direction.
• The Dance: They simulated the stars spinning around each other, sometimes wobbling (precession) or moving in slightly oval paths (eccentricity).

These simulations were so detailed that they produced "convergent" waveforms—meaning the results were stable and reliable, not just noisy guesses.

3. The Discovery: Listening to the "Splash"

By watching these 52 movies, the authors found something new about the sound of the merger:

• The Tidal Signature: When the neutron star gets torn apart, it leaves a specific "fingerprint" in the gravitational wave. The authors found that certain "notes" in the sound (specifically the (2,0) and (3,0) modes) get much louder when the star is ripped apart. It's like hearing a distinct crunch in the sound of a car crash that tells you the metal was bent, not just broken.
• The Kick: When the black hole eats the star, it doesn't just sit there; it gets kicked backward, like a gun recoiling. The authors found that if the star is torn apart early (tidal disruption), the "kick" is actually smaller than expected because the splash of matter absorbs some of the momentum.

4. The New Tool: TEOBResumS-Dalí

Using the data from these 52 movies, the authors built a new mathematical model called TEOBResumS-Dalí.

• The Analogy: Imagine you have a recipe for a cake (the old model). It tastes okay, but it's not quite right. The authors took the 52 new movies, analyzed exactly how the cake rose and browned, and wrote a new, improved recipe.
• The Result: This new model is much more accurate. When they tested it against a new, complex simulation (a 12-orbit dance with a wobbling spin), the model predicted the sound almost perfectly, with less than a 0.5-radian error in the timing. It's like a GPS that finally tells you exactly when you'll arrive, rather than just "sometime in the afternoon."

5. Why This Matters Now

The authors used their new model to look at a real event detected by LIGO/Virgo called GW230529.

• They found that their new model, which accounts for the "splash" of the neutron star, matches the real data much better than older models that ignored the splash.
• They also used the model to predict what happens if the stars are moving in oval paths or wobbling wildly. They generated the first-ever theoretical waveforms for these messy, eccentric, wobbling black-hole-neutron-star dances.

6. The Roadmap for the Future

Finally, the authors used their new model to act as a "guidebook" for other scientists. They ran a computer search to figure out: "Which 200 specific black-hole-neutron-star dances should we simulate next to learn the most?"
They found that the most urgent simulations to run are those where the neutron star is very "squishy" (high tidal disruption) and the black hole is spinning fast. These are the scenarios where our current knowledge is weakest.

Summary

In short, this paper is a massive upgrade to our understanding of how black holes eat neutron stars.

1. They made 52 new, high-quality movies of these events.
2. They discovered new "sounds" that tell us when a star is ripped apart.
3. They built a new, sharper model (TEOBResumS-Dalí) that predicts these events with high precision.
4. They used this model to decode a real cosmic event and to map out exactly what simulations scientists need to run next to keep improving our cosmic hearing.

The data from these simulations is now public, allowing the whole scientific community to use these new "movies" to tune their instruments and listen more clearly to the universe.

Technical Summary

Technical Summary: Black-hole - neutron-star mergers: new numerical-relativity simulations and multipolar effective-one-body model with spin precession and eccentricity

Problem Statement
The detection of Black Hole-Neutron Star (BHNS) mergers by LIGO-Virgo-KAGRA (LVK) has highlighted a critical gap in gravitational wave (GW) modeling. While Binary Black Hole (BBH) models are highly sophisticated, current BHNS waveform models lack the necessary precision, particularly regarding tidal disruption effects, spin precession, and orbital eccentricity. Existing numerical-relativity (NR) simulations for BHNS are often limited in length (fewer than four orbits), lack full convergence, or do not cover the parameter space where significant tidal disruption occurs. This scarcity of high-fidelity NR data hinders the development of accurate Effective One Body (EOB) models required to interpret current and future GW observations, especially for events where the neutron star (NS) is tidally disrupted, potentially producing electromagnetic counterparts.

Methodology
The authors present a comprehensive study combining new NR simulations with the development of an updated EOB model, TEOBResumS-Dalí.

1. Numerical-Relativity Simulations:

   • Dataset: 52 new BHNS simulations were performed using the BAM code with the Z4c formulation of Einstein's equations coupled to relativistic hydrodynamics.
   • Parameters: The simulations target quasicircular mergers in the region of significant tidal disruption. They cover a range of mass ratios ($q$), black hole spins ($a_{BH}$), and three different equations of state (EoS): SLy, MS1b, and ALF2 (including deconfined quark matter).
   • Initial Data: Generated using the Elliptica pseudospectral code, capable of handling high spin magnitudes ($\sim 0.8$) and arbitrary orientations.
   • Validation: Extensive convergence tests and error budgeting were performed. One specific simulation (BAM:0223) was evolved for 12 orbits with precessing spins to serve as an independent validation case.
   • Analysis: The authors extracted multipolar waveforms up to $\ell=4$, analyzed the remnant black hole (BH) properties (mass and spin), and calculated the GW recoil (kick) velocity.

2. Model Development (TEOBResumS-Dalí):

   • The authors extended the TEOBResumS framework to include BHNS physics, incorporating precession and eccentricity.
   • Key Components:
     • Remnant Models: New NR-informed fits for the final mass and spin of the remnant BH, smoothly connecting to BBH limits.
     • Multipolar Structure: Explicit modeling of subdominant modes $(2,1), (3,2), (3,3), (4,4)$ and the identification of characteristic tidal signatures in $(2,0)$ and $(3,0)$ modes.
     • Ringdown & NQC: A new ringdown model based on Quasi-Normal Modes (QNMs) and updated Next-to-Quasicircular (NQC) corrections. The model classifies mergers into three types (I: tidal disruption, II: direct plunge, III: intermediate) based on the peak amplitude of the $(2,2)$ mode.
     • Recoil Model: A fitting formula for the GW kick velocity, accounting for tidal suppression effects.

Key Contributions and Results

• New NR Catalog: The paper releases 52 convergent BHNS waveforms with detailed error budgets to the CoRe database. These include the first publicly available precessing BHNS simulation (BAM:0223) with 12 orbits.
• Multipolar Hierarchy: The study identifies a distinct multipolar amplitude hierarchy for BHNS compared to BBH. Notably, the $(2,0)$ and $(3,0)$ modes (associated with nonlinear memory) show significant contributions in BHNS, comparable to higher-order modes like $(3,2)$ and $(4,4)$, a feature not seen in BBH.
• Remnant and Kick Models: Updated analytical models for remnant mass and spin are provided, showing how tidal disruption reduces the final spin in certain configurations. A new model for GW recoil velocity demonstrates that tidal effects can suppress the kick velocity compared to BBH predictions, particularly for anti-aligned spins and high tidal polarizability ($\Lambda$).
• TEOBResumS-Dalí Performance:
  • Amplitude Improvement: The new model shows an order-of-magnitude improvement in waveform amplitude accuracy at merger compared to the previous TEOBResumS-GIOTTO model.
  • Validation: Validated against the independent 12-orbit precessing simulation (BAM:0223), the model achieves phase differences $< 0.5$ rad throughout the inspiral and mismatches at the $\sim 1\%$ level for low inclinations.
  • Eccentricity and Precession: The model successfully generates robust waveforms for eccentric and precessing BHNS systems, a capability not previously available for BHNS.
• Event GW230529: The authors note that simulation BAM:0223 matches the parameters of the observed event GW230529. Comparing the NR waveform to the LVK analysis (using SEOBNRv5PHM) yields a mismatch of $\sim 0.3$, attributed largely to unmodeled tidal effects in the LVK analysis.
• Event GW200105: The model was tested against a post-Newtonian (PN) analysis of GW200105 (which claimed eccentricity and precession). The PN model accumulated significant dephasing compared to TEOBResumS-Dalí, suggesting that PN-based inference for such events may be unreliable.

Significance and Claims
The paper claims that the combination of high-fidelity NR data and the TEOBResumS-Dalí model represents a significant step forward in BHNS waveform modeling. The authors emphasize that:

1. Tidal Signatures: The distinct behavior of the $(2,0)$ and $(3,0)$ modes offers a potential new channel for distinguishing BHNS from BBH in future observations.
2. Guiding Future Simulations: Using a greedy algorithm, the authors identify that future NR efforts should prioritize systems with high tidal polarizability ($\Lambda > 3000$), mass ratios $q \le 2$, and large effective spins to maximize the information gain for waveform modeling.
3. Robustness: The model is the first to produce robust, physically consistent waveforms for BHNS systems with both eccentricity and precession, filling a critical gap for analyzing potential future detections with these complex orbital parameters.
4. Limitations: The authors acknowledge that while the model is robust, the primary challenge remains the limited suite of NR simulations suitable for calibration, particularly regarding the convergence and length required for high-fidelity modeling of the strong-field dynamics and the ringdown phase. They call for further quantitative assessment of how matter and ejecta suppress QNMs to better characterize the ringdown.