Adding equatorial-asymmetric effects for spin-precessing binaries into the SEOBNRv5PHM waveform model

This paper introduces SEOBNRv5PHM_w/asym, an upgraded gravitational waveform model that incorporates equatorial-asymmetric effects to significantly improve agreement with numerical relativity simulations, enhance recoil velocity predictions, and refine parameter estimation for spin-precessing binary systems.

The Gist

Here is an explanation of the paper, translated into everyday language with some creative analogies.

The Big Picture: Fixing the "Spin" in the Dance

Imagine two black holes orbiting each other like a pair of ice skaters holding hands and spinning faster and faster until they crash together. This crash creates ripples in space-time called gravitational waves.

For years, scientists have built "models" (like mathematical blueprints) to predict exactly what these waves sound like. These models are crucial because when our detectors (like LIGO) hear a "chirp," they compare it to these blueprints to figure out: How heavy were the black holes? How fast were they spinning? Where did they come from?

However, there was a missing piece in the puzzle. Most models assumed the ice skaters were spinning perfectly upright. But in reality, they often lean over, wobble, and spin on their sides. This "leaning" creates a weird, lopsided effect that the old models couldn't quite capture.

This paper introduces a major upgrade to one of the best models out there, called SEOBNRv5PHM. They added a new feature called "equatorial asymmetry" (a fancy way of saying "lopsidedness caused by leaning spins"). They call the new, improved model SEOBNRv5PHMw/asym.

---

The Problem: The "Perfectly Symmetrical" Lie

In the old models, the scientists assumed that if you looked at the black hole collision from the top or the bottom, the waves would look like mirror images.

The Analogy: Imagine a perfectly symmetrical spinning top. If you spin it, the air pushes out evenly in all directions. It goes straight up.

The Reality: But black holes are more like a wobbly, leaning spinning top. When they spin, they don't just push air out evenly; they push harder on one side than the other. This creates a "kick."

• The Missing Effect: This "leaning" (spin precession) breaks the symmetry. It causes the black holes to shoot out a burst of gravitational waves that is stronger on one side of the orbital plane than the other.
• The Consequence: Because the waves push harder in one direction, the leftover black hole (the remnant) gets kicked in the opposite direction, like a rocket firing a thruster. The old models often got the speed of this "kick" wrong, sometimes missing it entirely.

The Solution: Adding the "Wobble" to the Model

The authors took their existing blueprint and added the physics of this "wobble."

1. The Math: They used complex equations (Post-Newtonian theory) to calculate how the leaning spins create these lopsided waves during the approach (inspiral).
2. The Calibration: They didn't just guess. They took 1,500+ supercomputer simulations (Numerical Relativity) of black holes crashing and "tuned" their new model to match these simulations perfectly. Think of it like a musician tuning a guitar by listening to a perfect reference tone.
3. The Result: The new model, SEOBNRv5PHMw/asym, now includes these "lopsided" waves in its calculations.

Why Does This Matter? (The Payoff)

1. Better Accuracy (The "Unfaithfulness" Score)

Scientists measure how bad a model is by something called "unfaithfulness" (basically, how much the model lies compared to reality).

• The Result: The new model is up to 50% more accurate than the old one, especially when we are looking at the collision "face-on" (from the top).
• The Metaphor: Imagine trying to recognize a friend in a crowd. The old model was like a blurry photo; you might guess it's them, but you aren't sure. The new model is a high-definition 4K photo. You can see the details clearly.

2. Predicting the "Kick" (The Recoil)

When the black holes merge, the leftover black hole can get kicked away at incredible speeds (thousands of miles per second).

• The Old Model: Often predicted the kick was tiny or zero.
• The New Model: Correctly predicts the massive "superkicks" that can happen.
• Why it matters: If a black hole gets kicked too hard, it might get ejected from its galaxy entirely. This changes how we understand where black holes live and how they grow.

3. Solving a Mystery: The Case of GW200129

There was a real gravitational wave event detected in 2020 (called GW200129) that was confusing. Some models said the black holes were spinning wildly (precessing), while others said they were just spinning straight up.

• The Test: The authors re-analyzed this event with their new model.
• The Verdict: The new model confirmed that yes, these black holes were definitely spinning wildly. It made the evidence for this "wobble" three times stronger than before. This helps us understand how these black holes formed in the first place.

Summary

Think of the old model as a flat, 2D drawing of a spinning dancer. It looked okay from the front, but it missed the depth and the wobble.

This paper adds the 3D depth and the wobble.

• It makes our "maps" of the universe more accurate.
• It helps us predict where the "remnant" black hole will fly off to.
• It helps us solve mysteries about real cosmic events we've already detected.

By fixing this "lopsided" effect, the scientists have given us a sharper, more truthful view of the most violent collisions in the universe.

Technical Summary

Here is a detailed technical summary of the paper "Adding equatorial-asymmetric effects for spin-precessing binaries into the SEOBNRv5PHM waveform model."

1. Problem Statement

Gravitational wave (GW) signals from compact binary coalescences with misaligned spins (spin-precessing systems) exhibit equatorial asymmetries that are absent in non-precessing (aligned-spin) systems. These asymmetries arise from in-plane spin components, breaking the symmetry of GW emission across the orbital plane.

• Consequences: This asymmetry leads to a net emission of linear momentum perpendicular to the orbital plane, resulting in a significant recoil (kick) velocity of the remnant black hole. In extreme "superkick" configurations, these kicks can eject the remnant from its host galaxy.
• Current Limitations: Most existing waveform models (including the previous version of SEOBNRv5PHM) construct spin-precessing signals by "twisting up" aligned-spin models. This approach inherently assumes equatorial symmetry in the co-precessing frame, thereby neglecting antisymmetric mode contributions. This omission leads to:
  • Systematic underestimation of recoil velocities.
  • Biases in parameter estimation (particularly spin orientation and mass ratio).
  • Reduced faithfulness (accuracy) when compared to Numerical Relativity (NR) simulations, especially for face-on/face-off viewing angles where antisymmetric modes are most prominent.

2. Methodology

The authors developed SEOBNRv5PHMw/asym, an extension of the SEOBNRv5PHM model that explicitly incorporates equatorial-asymmetric contributions into the waveform modes.

A. Theoretical Framework

• Mode Decomposition: The waveform modes in the co-precessing frame are decomposed into symmetric ($h^{sym}_{\ell m}$) and antisymmetric ($h^{asym}_{\ell m}$) parts based on their behavior under $z$-parity (reflection through the orbital plane).
• Inspiral-Plunge Modeling:
  • The antisymmetric contributions are modeled using Post-Newtonian (PN) theory.
  • The dominant modes ($\ell=m=2, 3, 4$) include Next-to-Leading Order (NLO) PN terms, while subdominant modes use Leading Order (LO) terms.
  • Corrections: To match NR data near the merger, the authors applied:
    • Non-Quasi-Circular (NQC) corrections: A phase-only correction to align frequencies at the matching time.
    • Amplitude corrections: A phenomenological correction for the dominant $(2,2)$ mode amplitude near merger.
  • Note: The model restricts explicit PN antisymmetric contributions to $\ell=m$ modes ($2,2; 3,3; 4,4$) because PN predictions for$\ell \neq m$modes (e.g.,$2,1$) failed to reproduce NR morphology accurately.
• Merger-Ringdown Modeling:
  • The antisymmetric contributions in the merger-ringdown phase are modeled using a phenomenological ansatz similar to the symmetric part (sum of QNMs).
  • The dominant $(2,2)$ antisymmetric mode peaks later than the symmetric mode; thus, its amplitude coefficients were recalibrated, while phase coefficients were fixed to match the symmetric case to ensure stability.

B. Calibration Strategy

The model was calibrated against a massive dataset to determine free parameters (amplitude and phase coefficients):

• Datasets:
  • 1,523 spin-precessing NR simulations from the SXS catalog (including the second catalog and specific simulations from Ref. [52]).
  • ~6,000 test-body (Teukolsky) waveforms for plunging geodesics in the large mass-ratio limit.
• Fitting Technique:
  • A sparse polynomial regression using Orthogonal Matching Pursuit (OMP) was employed to fit the target quantities (amplitudes, frequencies, and ringdown coefficients) as functions of the 7-dimensional parameter space (mass ratio, spin magnitudes, and angles).
  • Input features included effective spins ($\chi_{eff}, \chi_p$), mass ratio ($\nu$), and in-plane spin projections.

3. Key Contributions

1. First EOB Implementation of Antisymmetry: This is the first implementation of equatorial-asymmetric effects within the Effective-One-Body (EOB) formalism for spin-precessing binaries.
2. Comprehensive Calibration: The model is calibrated against the largest set of precessing NR simulations to date, including both SXS and BAM datasets, and test-body limits.
3. Recoil Velocity Prediction: The model explicitly targets the accurate prediction of gravitational recoil, a critical observable for astrophysical population studies.
4. Real-Event Re-analysis: The model was applied to the real GW event GW200129, a strong candidate for a spin-precessing system, to demonstrate its impact on parameter inference.

4. Results

A. Faithfulness and Accuracy

• Improvement over SEOBNRv5PHM: The new model reduces the median unfaithfulness (mismatch) with NR waveforms by up to 50% for face-on systems compared to the previous SEOBNRv5PHM.
• Comparison with Other Models:
  • Outperforms IMRPhenomXPNR by 30–60% in median unfaithfulness across various inclinations.
  • Outperforms TEOBResumS Dali by 76–80%.
  • Reduces the number of cases exceeding the 1% unfaithfulness threshold significantly, particularly for low-inclination (face-on) binaries.
• Generalization: Validation on independent datasets (BAM catalog) and training/testing splits confirms the model generalizes well without overfitting.

B. Gravitational Recoil (Kick)

• Accuracy: SEOBNRv5PHMw/asym reproduces the distribution of NR kick velocities with high fidelity, including the maximum kick magnitude.
• Error Reduction: The median relative error in predicting the recoil velocity dropped from 70% (SEOBNRv5PHM) to 1% (SEOBNRv5PHMw/asym).
• Limitation: A small percentage of cases still show >10% error, likely due to the omission of the antisymmetric $(3,2)$ mode.

C. Parameter Estimation (Bayesian Inference)

• Synthetic Injections: Using zero-noise injections of NR waveforms, the model improved the recovery of spin orientations ($\theta_1, \phi_{12}$) and mass ratios.
• GW200129 Re-analysis:
  • Re-analyzing GW200129 with SEOBNRv5PHMw/asym yielded a threefold increase in the Bayes factor favoring the spin-precessing hypothesis compared to the non-asymmetric SEOBNRv5PHM.
  • The results showed closer agreement with the NRSur7dq4 surrogate model, confirming that antisymmetric modes are crucial for correctly interpreting the spin-precession in this event.

5. Significance

• Astrophysics: Accurate recoil predictions are essential for understanding black hole retention in globular clusters and the formation of hierarchical mergers.
• Data Analysis: As detector sensitivity improves (e.g., with the Einstein Telescope and Cosmic Explorer), systematic errors from waveform models will dominate statistical errors. Including antisymmetric effects reduces these systematics, preventing biased inference of spin and mass parameters.
• Modeling Standard: This work sets a new standard for EOB models, demonstrating that incorporating specific physical asymmetries (beyond simple "twisting up") is necessary for high-precision GW astronomy.

The model is publicly available in the pySEOBNR Python package (version > v0.3.3) and can be activated via the enable_antisymmetric_modes=True flag.