Including higher-order modes in a quadrupolar eccentric numerical relativity surrogate using universal eccentric modulation functions

This paper introduces the \texttt{gwNRHME} framework, which leverages universal eccentric modulation functions to extend quasi-circular numerical relativity surrogates into multi-modal, non-spinning eccentric waveform models (specifically \model{}) that achieve high accuracy against numerical relativity simulations while also providing analytical tools for eccentricity evolution.

The Gist

Imagine you are trying to listen to a symphony played by two black holes spiraling into each other. For years, scientists have built "sheet music" (waveform models) to predict what this music sounds like, but they mostly assumed the black holes were dancing in a perfect circle.

However, nature is messy. Sometimes, these black holes are caught in a chaotic dance, swinging in and out on eccentric (oval-shaped) orbits, like a planet swinging wildly around a star. The problem is that the old sheet music doesn't work for these wild dancers. If you try to listen for them using circular music, you might miss them entirely or misunderstand the story they are telling.

This paper introduces a clever new toolkit called gwNRHME (a mouthful of a name, but think of it as a "Universal Eccentricity Translator"). Here is how it works, broken down into simple concepts:

1. The Problem: The "Circular" Bias

Most of our current models for black hole collisions are like a recipe for a perfect, round cake. They work great if the ingredients (the black holes) are moving in a circle. But if the black holes are moving in a weird, stretched-out oval (eccentric), the recipe fails.

Until now, scientists had to build a completely new recipe from scratch for every type of oval orbit. This is slow, expensive, and requires massive computer simulations that take months to run.

2. The Solution: The "Universal Modulation" Magic

The authors discovered a hidden pattern in the chaos. They found that no matter how weird the orbit is, the "wobble" (eccentricity) affects all the different notes (modes) of the gravitational wave in the exact same way.

Think of it like this: Imagine a drummer playing a steady beat (the circular orbit). Now, imagine the drummer gets excited and starts speeding up and slowing down rhythmically (the eccentricity).

• The bass drum (the main, loudest note) wobbles in a specific pattern.
• The snare drum (a quieter, higher note) wobbles in the exact same pattern, just slightly quieter.
• The cymbals (the highest, faintest notes) also wobble in that same pattern.

The paper proves that if you know how the bass drum wobbles, you can mathematically predict exactly how the snare and cymbals will wobble. You don't need to simulate the whole orchestra again; you just need to apply the "wobble pattern" to the existing music.

3. The Toolkit: Mixing and Matching

The authors built a framework that acts like a modular Lego set:

• Piece A: A highly accurate model of the "bass drum" (the main gravitational wave) for eccentric orbits.
• Piece B: A library of "sheet music" for circular orbits that includes all the complex, high-pitched notes (higher-order modes).
• The Glue: The "Universal Modulation" function.

They take Piece B (the complex circular music) and use the Glue to stretch and squeeze it, turning it into the correct eccentric music. They then attach Piece A to ensure the main beat is perfect.

4. The Results: A New Super-Model

By using this method, they created a new model called gwNRHME_NRSur_q4.

• What it does: It predicts the sound of colliding black holes that are spinning (or not spinning) and moving in oval orbits, capturing nine different "notes" (modes) at once.
• How good is it? It is incredibly accurate. When they tested it against the "gold standard" (supercomputer simulations), the error was tiny—like trying to hear a pin drop in a library and being off by a fraction of a whisper.
• Why it matters: It allows astronomers to listen for these eccentric black holes without needing to run a new supercomputer simulation for every single event. It's like having a universal translator that can instantly turn a circular song into an eccentric one.

5. The Bigger Picture

The paper also shows that this "translator" is modular. You can swap out the "sheet music" part.

• They tested it with a standard circular model (NRHybSur3dq8) and got great results.
• They also tested it with two other famous models (SEOBNRv5HM and TEOBResumS-Dali) and it worked there too!

This means the framework is flexible. If scientists invent a better circular model in the future, they can just plug it into this translator, and poof—they instantly have a better eccentric model.

Summary

In short, this paper solves the problem of "how do we model messy, oval-shaped black hole collisions?" by realizing that the messiness follows a simple, universal rule. Instead of building a new house for every different shape of orbit, they built a machine that can take a standard house and instantly remodel it to fit any shape, saving time and opening the door to discovering a whole new class of cosmic events.

Technical Summary

1. Problem Statement

Detecting and characterizing gravitational waves (GWs) from eccentric binary black hole (BBH) mergers is a critical objective in GW astronomy. While most detected events are well-described by quasi-circular templates, astrophysical scenarios (e.g., dynamical captures in dense clusters) suggest a population of eccentric binaries may exist.
Current challenges include:

• Lack of Accurate Models: Existing eccentric waveform models often rely on semi-analytical approximations (Post-Newtonian or Effective-One-Body) that assume circularization by merger or lack calibration to high-fidelity Numerical Relativity (NR) simulations.
• Missing Higher-Order Modes (HOMs): Most eccentric models are restricted to the dominant quadrupolar mode $(2,2)$. However, accurate parameter estimation and detection of high-mass-ratio or high-eccentricity systems require higher-order spherical harmonic modes.
• Computational Cost: Generating full NR simulations for every eccentric configuration with HOMs is computationally prohibitive.

2. Methodology

The authors propose and utilize the gwNRHME (gravitational-wave Numerical Relativity Higher-Order Modes with Eccentricity) framework. This framework leverages the universality of eccentric modulation functions to convert a quasi-circular, multi-modal waveform model into an eccentric, multi-modal model, provided a quadrupolar eccentric model is available.

Key Methodological Steps:

• Universal Modulation Functions: The framework exploits the empirical finding that eccentricity induces characteristic, mode-independent modulations in both waveform amplitudes ($A_{\ell m}$) and instantaneous frequencies ($\omega_{\ell m}$). These modulations are quantified by comparing eccentric waveforms to their quasi-circular counterparts.
  • Amplitude modulation: $\xi^A_{\ell m} \approx B \cdot \xi^\omega_{\ell m}$, where $B \approx 0.9$.
  • A common modulation function $\xi(t)$ is defined using the robust $(2,2)$ mode.
• Model Construction:
  • Base Eccentric Model: The quadrupolar, non-spinning eccentric surrogate NRSurE_q4NoSpin_22 (trained on ~156 NR simulations) provides the eccentric modulation $\xi(t)$ and the dominant $(2,2)$ mode.
  • Base Quasi-Circular Model: The multi-modal, quasi-circular NR surrogate NRHybSur3dq8 provides the baseline higher-order modes.
  • Synthesis: The framework maps the quasi-circular HOMs to eccentric ones using the modulation function:
    $$A^{gwNRHME}_{\ell m}(t) = A_{\ell m}(t; \lambda_0) \left[ 1 + \frac{\ell}{2} \xi(t) \right]$$
    $$\omega^{gwNRHME}_{\ell m}(t) = \omega_{\ell m}(t; \lambda_0) \left[ 1 + \frac{\xi(t)}{B} \right]$$
• Eccentricity Evolution: A new surrogate model (gwEccEvolve_q4NoSpin_Sur) and an analytical model (gwEccEvNSv2) are constructed to describe the evolution of eccentricity $e_\xi(t)$ up to merger, based on the universal modulation functions.

3. Key Contributions

1. gwNRHME_NRSur_q4: A new, accurate, multi-modal, non-spinning eccentric surrogate model.
   • Modes: Includes 9 spherical harmonic modes: $(2, \{1,2\})$, $(3, \{1,2,3\})$, $(4, \{2,3,4\})$, and $(5,5)$.
   • Parameter Space: Valid for mass ratios $q \in [1, 4]$ and reference eccentricities $e_{ref} \in [0.001, 0.43]$.
2. Modularity Demonstration: The framework successfully combines the eccentric quadrupolar model with other state-of-the-art quasi-circular models, specifically:
   • SEOBNRv5HM (Effective-One-Body), yielding gwNRHME_SEOB_q4.
   • TEOBResumS-Dali, yielding gwNRHME_TEOB_q4.
3. Eccentricity Evolution Models:
   • gwEccEvolve_q4NoSpin_Sur: A data-driven surrogate for eccentricity evolution.
   • gwEccEvNSv2: An updated analytical phenomenological model for eccentricity evolution.
4. Public Availability: The framework is released via the gwModels package, and the waveform models via gwsurrogate.

4. Results and Performance

The models were validated against 156 eccentric SXS NR waveforms.

• Accuracy (gwNRHME_NRSur_q4):
  • Achieves median frequency-domain mismatches of $\sim 9 \times 10^{-5}$ (using Advanced LIGO design sensitivity) with a standard deviation of $\sim 2 \times 10^{-4}$.
  • The time-domain relative $L_2$-norm errors are $\sim 10^{-4}$ for the $(2,2)$ mode and comparable to NR resolution limits for higher-order modes.
  • The model correctly reproduces physical features such as mode mixing during the ringdown phase and the vanishing of odd-$m$ modes in the equal-mass limit.
• Impact of Higher-Order Modes:
  • Using only the quadrupolar mode (NRSurE_q4NoSpin_22) against multi-modal NR data results in mismatches of $\sim 10^{-3}$ to $10^{-2}$, which is insufficient for high-mass-ratio systems.
  • Including HOMs via gwNRHME reduces mismatches by orders of magnitude, making the model suitable for GW data analysis.
• Modularity Validation:
  • When combined with EOB models, the framework yields median mismatches of $\sim 2 \times 10^{-4}$ (SEOBNRv5HM) and $\sim 10^{-3}$ (TEOBResumS-Dali). While slightly less accurate than the NR-based hybrid, these results confirm the framework's robustness in extending analytical models.
• Eccentricity Evolution:
  • The new eccentricity evolution models reproduce NR data with fractional errors bounded within 5% for 95% of systems (for $t \le -100M$).

5. Significance

• Detection Capability: By providing accurate multi-modal eccentric waveforms, this work enables the detection and characterization of eccentric BBH mergers that would otherwise be missed or poorly characterized by quasi-circular templates.
• Astrophysical Insights: Accurate models for eccentric systems are crucial for probing dynamical formation channels (e.g., globular clusters, galactic nuclei) where binaries may retain measurable eccentricity near merger.
• Efficiency and Flexibility: The gwNRHME framework offers a computationally efficient "plug-and-play" solution. It allows researchers to upgrade existing quasi-circular models (NR or EOB) to eccentric versions without re-running expensive NR simulations for every new configuration.
• Foundation for Future Work: The paper sets the stage for extending this framework to spinning and precessing eccentric binaries, further expanding the reach of GW astronomy.

In conclusion, this paper presents a robust, modular, and highly accurate framework for generating eccentric gravitational waveforms with higher-order modes, significantly advancing the toolkit available for analyzing eccentric binary black hole mergers.