The Cost of Circularity: Quantifying Eccentricity-Induced Biases in Binary Black Hole Inference

This study demonstrates that unmodeled orbital eccentricity in binary black hole mergers significantly biases inferred parameters—such as mass, spin, and distance—when analyzed with standard circular waveform models, particularly for systems with eccentricities above $e \sim 0.2$, thereby establishing the critical thresholds where eccentric waveform models become essential for accurate astrophysical inference.

The Gist

Imagine you are a detective trying to solve a crime, but the only witness is a very specific type of musical instrument: a pair of black holes dancing together, creating ripples in the fabric of space-time called gravitational waves.

For years, scientists have been listening to these cosmic dances. Most of the time, they assume the black holes are dancing in a perfect, smooth circle, like ice skaters holding hands and spinning on a frozen lake. They use a "circular" map to figure out who the dancers are, how heavy they are, and how far away they are.

But what if the dancers aren't skating in a perfect circle?

What if they are actually doing a chaotic, wobbly dance, zooming in and out like a rollercoaster? This is called eccentricity.

This paper asks a critical question: If we assume the black holes are dancing in a perfect circle, but they are actually wobbling, how badly will we get the details of the dance wrong?

The "Round Peg in a Square Hole" Problem

The researchers decided to test this by creating fake signals (like playing a recording of a wobbly dance) and then trying to analyze them using the old, "perfect circle" maps.

Here is what they found, explained through some everyday analogies:

1. The "Wobble" Confuses the Scale
Imagine you are looking at a spinning top. If it's spinning perfectly straight, you can easily guess how heavy it is. But if it's wobbling wildly, your brain might trick you into thinking it's lighter or heavier than it really is.

• The Finding: When the black holes have even a little bit of wobble (eccentricity), the "circular" maps start to guess the wrong mass. It's like trying to weigh a person while they are jumping on a trampoline; the scale gives you a confusing number.

2. The "Fake Spin" Illusion
This is the most interesting part. The researchers found that when a black hole system is wobbly (eccentric), the circular maps try to "fix" the error by inventing a new feature: spin.

• The Analogy: Imagine you see a car swerving on the road. If you assume the car is driving straight, you might conclude, "The driver must be drunk and swerving!" But actually, the road itself was bumpy.
• The Reality: The circular maps think the black holes are spinning wildly (precessing) to explain the wobble. In reality, the black holes might not be spinning at all; they are just moving in an oval path. The map is creating a "ghost" spin to make the math work.

3. The "Heavy Dancers" are the Problem
The study found that this confusion gets much worse when the black holes are very heavy.

• The Analogy: If you are trying to guess the weight of a feather that is fluttering in the wind, it's hard. But if you are trying to guess the weight of a giant boulder that is rolling down a hill in a zig-zag pattern, and you assume it's rolling straight, your guess will be wildly off.
• The Finding: For massive black holes (like the ones that make up the heaviest systems we've ever seen), even a small amount of wobble causes the scientists to get the distance, the mass, and the spin completely wrong.

The "Cost" of Being Circular

The title of the paper is "The Cost of Circularity." This "cost" isn't money; it's truth.

If we keep using the old "perfect circle" maps for these wobbly, heavy systems, we will:

• Think the black holes are in different places than they really are.
• Think they are spinning when they aren't.
• Misunderstand how they formed. (Wobbly orbits usually mean the black holes met by accident in a crowded star cluster, while perfect circles mean they were born together. If we get the orbit wrong, we get their "family tree" wrong.)

The Solution: A New Map

The researchers tested a new, more advanced map (called TEOBResumS-Dalí) that knows how to handle wobbles.

• The Result: When they used this new map, it correctly identified the wobble and gave the right answers for mass, distance, and spin. It didn't invent fake spins to fix the math.

The Bottom Line

We are entering a new era of astronomy where we are finding heavier and heavier black holes. The authors warn us: We can no longer pretend the dance is a perfect circle.

If we want to understand the universe correctly, we need to upgrade our "detective tools" to handle the wobble. Otherwise, we are going to keep misidentifying the stars in our cosmic neighborhood, thinking they are doing a spin move when they are actually just tripping over their own feet.

In short: Don't force a square peg into a round hole. If the black holes are wobbly, we need a wobbly map to find them.

Technical Summary

1. Problem Statement

Dynamically assembled binary black hole (BBH) systems, formed in dense stellar environments (e.g., globular clusters, galactic nuclei, AGN), are expected to retain measurable orbital eccentricity ($e$) when entering the LIGO-Virgo-KAGRA (LVK) frequency band. However, standard parameter estimation (PE) pipelines in the LVK collaboration predominantly assume quasi-circular inspirals.

The core problem addressed is the systematic bias introduced when analyzing eccentric signals using circular waveform models. Specifically, the authors investigate:

• How unmodeled eccentricity distorts the inference of physical parameters (masses, spins, distance, inclination).
• The degeneracy between orbital eccentricity and spin precession, where circular precessing templates may mimic eccentric amplitude/phase modulations by introducing artificial precession.
• The threshold of eccentricity and mass beyond which circular models fail, potentially leading to incorrect astrophysical interpretations regarding formation channels.

2. Methodology

The study employs a rigorous injection-recovery framework to quantify these biases:

• Signal Injection (Ground Truth):

  • Waveform Model: TEOBResumS–Dalí, a state-of-the-art effective-one-body (EOB) model capable of handling eccentric inspirals.
  • Source Configuration: Equal-mass ($q=1$), non-spinning binary black holes.
  • Parameter Space:
    • Total Mass ($M$): $20 - 80 , M_\odot$(in steps of$10 , M_\odot$).
    • Initial Eccentricity ($e_{5Hz}$): $0.0$to$0.5$(in steps of$0.1$).
  • Environment: Signals injected into zero-noise data for a two-detector network (Hanford and Livingston) with design sensitivities.

• Signal Recovery (Inference):

  • Waveform Models: The injected eccentric signals were recovered using the IMRPhenomXPHM model, a phenomenological frequency-domain model that assumes circular orbits but includes spin precession.
  • Recovery Configurations: Three distinct setups were tested:
    1. Precessing spins (non-aligned).
    2. Aligned spins.
    3. Non-spinning.
  • Inference Engine: The Bilby Bayesian parameter estimation pipeline using nested sampling (Dynesty).

• Validation:

  • A specific case ($M=26.1 \, M_\odot, e_{20Hz}=0.1$) was recovered using the eccentric TEOBResumS–Dalí model to compare against the circular recovery and calculate Bayes factors.

3. Key Results

A. Parameter Biases and Mass Dependence

• Mass-Eccentricity Correlation: The bias in recovered parameters is strongly correlated with both total mass and eccentricity.
  • Low Mass ($<30 \, M_\odot$): Circular models remain robust for low eccentricities ($e_{10Hz} \lesssim 0.1$).
  • High Mass ($>40 \, M_\odot$): Significant systematic offsets appear even at relatively low eccentricities. For $M > 40 \, M_\odot$ and $e_{5Hz} > 0.3$ ($e_{10Hz} \approx 0.15$), the injected parameters often fall outside the 90% credible intervals.
• Specific Parameter Deviations:
  • Masses: Recovered component masses show significant deviations, often biased toward lower mass ratios.
  • Distance ($D_L$): Luminosity distance is frequently overestimated due to a loss in matched-filter efficiency.
  • Inclination ($\theta_{jn}$): High-eccentricity systems recovered with precessing models show a shift toward edge-on orientations.

B. The Eccentricity-Precession Degeneracy

• Artificial Precession: The study confirms that circular precessing templates attempt to compensate for the unmodeled eccentric phase/amplitude modulations by introducing artificial spin precession.
• Spin Parameter Recovery:
  • For $e_{5Hz} \leq 0.2$, the recovered effective precession spin ($\chi_p$) remains consistent with zero (the true injection value).
  • For $e_{5Hz} \geq 0.3$, the posterior distributions for $\chi_p$ shift significantly toward high values, falsely suggesting a precessing system.
  • The effective spin ($\chi_{eff}$) also shows non-negligible deviations for $e_{5Hz} \geq 0.4$.

C. Model Selection (Bayes Factors)

• The authors calculated the Bayes factor ($B$) between a precessing circular recovery and a non-spinning circular recovery.
• Result: As eccentricity increases, the Bayes factor favors the precessing model ($B > 1$), even though the true system is non-spinning. This indicates that the circular precessing model is "fitting" the eccentricity with spin parameters, reinforcing the degeneracy.

D. Validation with Eccentric Models

• When the same signal was recovered using the eccentric TEOBResumS–Dalí model, the parameters were recovered accurately.
• The Bayes factor strongly favored the eccentric model over the precessing circular model ($B_{PE} \approx 10^{-7}$), demonstrating that using the correct waveform family is essential for accurate inference.

4. Key Contributions

1. Quantification of Bias Boundaries: The paper establishes specific boundaries in the parameter space where circular models fail. It identifies that for high-mass systems ($M > 40 \, M_\odot$), even moderate eccentricity ($e_{10Hz} \sim 0.15$) causes significant biases.
2. Degeneracy Characterization: It explicitly demonstrates how circular precessing models mimic eccentricity, leading to false positives for spin precession and incorrect astrophysical conclusions about formation channels.
3. Methodological Framework: The study provides a systematic injection-recovery analysis using TEOBResumS–Dalí (eccentric) and IMRPhenomXPHM (circular precessing) within the Bilby framework, setting a standard for future bias studies.

5. Significance and Implications

• Astrophysical Interpretation: Ignoring eccentricity in high-mass BBH mergers can lead to misidentifying the formation channel. A dynamically formed (eccentric) binary might be incorrectly classified as an isolated binary with high spin precession.
• Population Studies: Systematic biases in mass, distance, and spin distributions will skew population synthesis models if not corrected, affecting estimates of merger rates and formation scenarios.
• Future Observing Runs: As LVK moves to O4 and O5, and with the advent of next-generation detectors (Einstein Telescope, Cosmic Explorer), the sensitivity to high-mass and potentially eccentric systems will increase. The authors conclude that eccentric waveform models must be integrated into PE pipelines to avoid systematic errors, particularly for high-mass events.

In summary, the paper argues that the "cost of circularity" is high for massive, moderately eccentric binaries, necessitating a shift toward eccentric waveform models to ensure the accuracy of gravitational-wave astrophysics.