Measuring Eccentricity and Addressing Waveform Systematics in GW231123

This paper reanalyzes the heavy binary black hole event GW231123 using a complete physical model to demonstrate that the event exhibits no significant eccentricity and that previous discrepancies in parameter estimates stem from waveform model disagreements regarding strong spin precession rather than orbital eccentricity.

The Gist

The Big Picture: A Cosmic Mystery

Imagine the LIGO and Virgo detectors as giant, ultra-sensitive ears listening to the universe. On November 23, 2023, they heard a very loud "thump" from space, an event called GW231123.

This wasn't just any thump; it was the heaviest collision of two black holes ever recorded. It was so massive that the black holes involved were heavier than what standard physics says should be possible (like finding a 100-pound goldfish in a pond where the biggest fish usually weighs 10 pounds). This suggests these black holes were formed in a very unusual way, perhaps by merging with other black holes first, rather than just dying stars.

However, when scientists tried to measure the details of this crash, they hit a wall. Different computer models gave them different answers. It was like asking three different weather forecasters to predict the storm, and one said "rain," one said "snow," and one said "hail."

This paper is the team's attempt to figure out why the models disagreed and what actually happened.

---

The Two Suspects: "Eccentricity" and "Spinning"

To understand the mystery, we need to know two things about how black holes dance before they crash:

1. Spin Precession (The Wobbly Top): Imagine a spinning top. If it's spinning perfectly upright, it's stable. But if it's tilted, it wobbles as it spins. Black holes can wobble too. In this event, the black holes were spinning so fast and tilted so much that they were wobbling wildly.
2. Eccentricity (The Oval Track): Usually, planets orbit the sun in a perfect circle. But sometimes, they orbit in an oval (an ellipse). If black holes orbit in an oval, that's called "eccentricity."

The Problem: The scientists suspected the black holes might have been orbiting in an oval (eccentric). If they were, the standard computer models (which assume perfect circles) would get the math wrong, leading to the conflicting results.

---

The Investigation: A New, Better Model

The authors of this paper decided to build a "super-model" called TEOBResumS-Dalí. Think of this model as a high-definition, 3D simulator that can handle both the wild wobbling (spin) and the oval tracks (eccentricity) at the same time.

They ran the data through this new model and found two major things:

1. The "Oval Track" Theory is Unlikely

They asked: "Was the black hole orbit actually oval-shaped?"

• The Result: No. The data didn't show strong evidence of an oval track. The best guess was that the orbit was actually very close to a perfect circle.
• The Analogy: Imagine trying to tell if a race car is driving on a slightly bumpy oval track or a perfectly smooth circle. The authors found that even if the track was slightly bumpy (up to a certain limit), the car's speed and position looked so similar to a smooth circle that you couldn't confidently say, "It's definitely an oval!"
• Conclusion: The disagreement between the old models wasn't caused by the black holes having an oval orbit.

2. The Real Culprit: The "Wobbly Top" Confusion

If it wasn't the oval track, why did the models disagree?

• The Result: The disagreement came from how the models handled the wobbling (spin precession).
• The Analogy: Imagine trying to describe a dancer who is spinning, wobbling, and jumping all at once.
  • Model A is good at describing the jump but bad at the wobble.
  • Model B is good at the wobble but assumes the dancer is standing still.
  • Because the black holes were wobbling so much (a "strong spin precession"), the models got confused. They tried to compensate for the missing wobble physics by changing the numbers for mass and distance. It's like a photographer trying to take a picture of a fast-moving, spinning object with a blurry camera; they might guess the object is bigger or further away just to make the blur make sense.

The authors proved this by creating a "fake" signal in a computer (a zero-noise injection) that had the exact same wild wobble as the real event. When they ran this fake signal through the old models, they got the exact same wrong answers as they did with the real data. This confirmed that the models themselves are the problem, not the data.

---

The "Double-Blind" Trap

There was one more twist. The authors tested a specific scenario: What if we force the model to assume the orbit is oval, but we ignore the wobbling?

• The Result: The model would confidently say, "Yes! It's definitely an oval orbit!"
• The Trap: This is a degeneracy. Because the signal was so short (the black holes crashed quickly), the "wobble" looked exactly like an "oval track" to the computer.
• The Analogy: It's like hearing a sound that could be a drum beating in a rhythm (wobble) or a drum being hit on an uneven surface (oval). If you only listen for a split second, you can't tell the difference. If you guess it's the uneven surface, you get the answer "oval," but you're actually wrong.
• The Fix: When the authors used their "super-model" that accounted for both wobble and oval tracks, the "oval" guess disappeared, and the "wobble" explanation won.

---

The Takeaway: What Does This Mean for Us?

1. The Black Holes: They are likely not on an oval track. They are probably on a near-perfect circle, but they are wobbling like crazy.
2. The Models: Our current computer models for black holes are struggling when the black holes spin very fast and wobble wildly. They are "breaking" under the pressure of these extreme conditions.
3. The Future: To understand these cosmic monsters, we need better "maps" (waveform models). We need to teach our computers how to handle the wildest, wobbliest dances in the universe. Until we do, we might keep misinterpreting the size and shape of these collisions.

In short: The universe sent us a confusing message because our translators (the computer models) aren't fluent enough in the language of "wobbly, spinning black holes." This paper didn't find a new type of black hole, but it fixed the translation errors so we can finally read the message correctly.

Technical Summary

1. Problem Statement

The gravitational wave event GW231123 is the most massive binary black hole (BBH) system observed by the LIGO–Virgo–KAGRA (LVK) Collaboration to date. Initial LVK analysis suggested component masses ($137 M_\odot$ and $103 M_\odot$) and high spins ($\chi \approx 0.9$ and $0.8$) that place the system within or above the pair-instability mass gap, challenging standard stellar evolution models.

However, a critical issue emerged: different waveform models yielded significantly different parameter estimates (e.g., total mass, mass ratio, spin magnitudes, and luminosity distance). The LVK analysis relied on quasi-circular models, potentially missing physics such as:

• Orbital Eccentricity: Possible if the system formed via dynamical channels.
• Strong Spin Precession: The inferred spins exceed the calibration limits ($\chi > 0.8$) of many standard waveform models.
• Degeneracies: The potential confusion between signatures of eccentricity and spin precession in short-duration, high-mass signals.

The authors aim to reanalyze GW231123 to determine if eccentricity is present and to identify the source of the systematic discrepancies between models.

2. Methodology

The authors employed a rigorous Bayesian analysis framework using the RIFT parameter estimation algorithm.

• Waveform Models:
  • Primary Model: TEOBResumS-Dalí, a physically complete model capable of handling eccentric and spin-precessing binaries simultaneously.
  • Comparison Models: NRSur7dq4 (numerical relativity surrogate) and SEOBNRv5PHM (effective-one-body), both of which are quasi-circular but spin-precessing.
  • Hypothesis Testing: They analyzed the data under three distinct hypotheses:
    1. Eccentric + Spin Precessing (TEOB).
    2. Quasi-circular + Spin Precessing (TEOB-P, NRSur, SEOB).
    3. Eccentric + Aligned-Spin (TEOB-E) to test degeneracies.
• Data Processing:
  • Used 8 seconds of data from LIGO Hanford and Livingston (20–448 Hz).
  • Priors on eccentricity ($e_{10}$) and anomaly were uniform; other priors matched the original LVK analysis.
• Validation Studies (Zero-Noise Injection-Recovery):
  • Eccentricity Threshold: Injected signals with $e_{10} \in \{0.0, 0.1, 0.15\}$ to determine the minimum eccentricity detectable for this system.
  • Systematic Error Source: Injected the maximum-likelihood point from the TEOB analysis (high spin, strong precession) into zero noise and recovered it with different models to see if model disagreement reproduces the real-data discrepancies.
  • Degeneracy Test: Injected a quasi-circular, strongly precessing signal and analyzed it with an eccentric, aligned-spin model to see if it falsely infers eccentricity.

3. Key Contributions & Results

A. Eccentricity Measurement

• Result: The analysis using the full eccentric, spin-precessing model (TEOB) yielded a median eccentricity at 10 Hz of $e_{10} = 0.062^{+0.063}_{-0.062}$.
• Conclusion: The 90% highest-density interval (HDI) includes zero. The data does not provide strong evidence for eccentricity.
• Detection Limit: Injection studies showed that even for injected eccentricities as high as $e_{10} = 0.15$, the posterior distribution still retained substantial support for zero eccentricity. A confident non-zero measurement is only possible if $e_{10} > 0.15$ for systems like GW231123.

B. Impact of Excluding Eccentricity

• Comparison: Comparing the full TEOB model (eccentric) with TEOB-P (quasi-circular) showed negligible differences in the posterior distributions for mass, spin, and distance.
• Bayes Factor: The log Bayes factor ($\ln \text{BF} \approx 0.7$) slightly favored the quasi-circular hypothesis, indicating that excluding eccentricity does not explain the parameter discrepancies observed in the LVK analysis.

C. Source of Waveform Systematics

• High Spin Calibration: Testing models at high spins ($\chi > 0.8$) with aligned spins showed only modest discrepancies.
• Strong Spin Precession: When the injection included strong spin precession ($\chi_p = 0.77$), significant discrepancies emerged between models (TEOB, NRSur, SEOB).
  • The bias patterns in the zero-noise recovery (e.g., NRSur/SEOB inferring lower primary mass and higher secondary spin compared to TEOB) closely matched the discrepancies seen in the real GW231123 data.
  • Conclusion: The primary source of systematic error is the disagreement between waveform models in the regime of strong spin precession, not the omission of eccentricity.

D. Eccentricity-Spin Precession Degeneracy

• The Trap: Analyzing the event with an eccentric, aligned-spin model (TEOB-E) yielded a "confident" non-zero eccentricity ($e_{10} \approx 0.55$).
• The Reality: This was a false positive caused by the degeneracy between eccentricity and spin precession in short-duration signals. The aligned-spin model incorrectly attributed precession-induced modulations to eccentricity.
• Resolution: Bayesian model selection strongly favored the eccentric, spin-precessing hypothesis ($\ln \text{BF} \sim -26$) over the eccentric, aligned-spin hypothesis. The full model correctly identified that the signal is consistent with zero eccentricity.

4. Significance and Implications

• Astrophysical Interpretation: GW231123 is likely a high-mass, high-spin binary formed via channels beyond standard stellar collapse (e.g., hierarchical mergers or active galactic nucleus disks), but not necessarily via a dynamical channel requiring high eccentricity.
• Modeling Challenges: The study highlights a critical limitation in current gravitational wave astronomy: waveform models disagree significantly in the high-spin, strong-precession regime. This leads to biased parameter estimates even at moderate signal-to-noise ratios (SNR $\approx 21.8$).
• Future Directions:
  • Development of waveform models that self-consistently couple spin precession and eccentricity.
  • Need for high-accuracy numerical relativity simulations covering generic spins ($\chi > 0.8$) and eccentric orbits to calibrate future models.
  • Caution is required when interpreting "confident" eccentricity measurements in short, high-mass signals, as they may be artifacts of spin-precession degeneracy.

In summary, the paper resolves the mystery of GW231123 by demonstrating that the event is consistent with a quasi-circular, highly precessing binary, and that the previously observed parameter discrepancies were caused by waveform model limitations in handling strong spin precession, rather than unmodeled orbital eccentricity.