Orbital eccentricity in a neutron star - black hole merger

This paper reports the first confident measurement of orbital eccentricity in the neutron star–black hole merger GW200105 using a novel waveform model that incorporates both precession and eccentricity, suggesting a formation channel driven by dynamical interactions rather than isolated binary evolution.

The Gist

Imagine two cosmic dancers: a neutron star (an incredibly dense, city-sized corpse of a star) and a black hole (a gravity monster so heavy it swallows light). Usually, when these two meet, they spiral toward each other in a perfect, smooth circle, like a figure skater gliding on ice, before finally colliding in a spectacular crash.

For years, scientists assumed that all these cosmic couples danced this way. They believed that by the time they got close enough for our detectors to hear them, any initial wobble or oval shape in their dance would have been smoothed out by the friction of spacetime itself.

But this paper says: "Not so fast."

The researchers, led by Gonzalo Morras, Geraint Pratten, and Patricia Schmidt, have found a "smoking gun" in the gravitational waves from an event called GW200105. They discovered that this specific pair wasn't dancing in a perfect circle. They were dancing in an oval.

Here is the story of that discovery, broken down into simple concepts:

1. The "Eccentric" Dance

In physics, we call an oval orbit eccentricity. Think of it like a race track.

• Circular (Low Eccentricity): A perfect circle. The runners stay the same distance from the center.
• Eccentric (High Eccentricity): An oval. The runners get very close to the center, then swing far out, then get close again.

The team found that GW200105 had an eccentricity of about 0.145. While that sounds small, in the world of black holes, it's huge. It means the two objects were swinging in and out of each other's gravity like a pendulum, rather than spinning in a tight, neat circle.

2. The New "Microscope"

Why didn't anyone see this before? Because looking for an oval orbit in a gravitational wave signal is like trying to find a specific needle in a haystack while wearing foggy glasses.

Previous models used to analyze these crashes assumed the orbits were perfect circles. If you try to fit a round peg into a square hole, you get a bad fit. The researchers built a new, super-advanced model (called pyEFPE) that acts like a high-definition microscope. This model can see both the "wobble" of the spin and the "ovalness" of the orbit at the same time. It's the first time they've been able to look for both clues simultaneously without one confusing the other.

3. The "Cosmic Dating" Mystery

This discovery changes the story of how these couples meet.

• The Old Story (Isolated Evolution): Imagine two stars born together, holding hands, and slowly aging. Over billions of years, they lose energy and their orbit becomes a perfect circle. This is the "standard" story scientists told for years.
• The New Story (Dynamical Interaction): The oval orbit of GW200105 suggests a much more chaotic meeting. It's like two strangers bumping into each other at a crowded party and deciding to dance.
  • The Analogy: Imagine a pool table. If two balls are gently rolling toward each other, they might hit and spin. But if a third ball hits them, or if they are in a crowded room where other balls are crashing into them, they can get knocked into a weird, oval path.
  • The Conclusion: This binary system likely formed in a crowded stellar neighborhood (like a dense star cluster) where other stars bumped into them, or perhaps they were part of a trio (a triple system) where a third partner messed up their dance before one was kicked out.

4. The "Heavier" Partner

The new model also changed the weight of the dancers.

• Old Estimate: The black hole was about 9 times the mass of our Sun, and the neutron star was about 2 times the mass.
• New Estimate: The black hole is actually heavier (about 11.5 Suns), and the neutron star is lighter (about 1.5 Suns).

Why does this matter?
When you ignore the "ovalness" of the orbit, your math gets confused. It's like trying to weigh a spinning top while it's wobbling; you might think it's lighter or heavier than it really is. By correcting for the wobble, the scientists realized the neutron star is actually a very typical, "normal" size, which fits better with what we know about neutron stars in our galaxy.

5. The "Noise" Check

You might ask, "Could this just be a glitch in the data? A random bump in the noise?"

The team did a massive "lie detector" test. They took the data and asked: "If this were just random static noise, how often would we see an oval orbit this big?"
The answer: Almost never. The odds of this being a random glitch are less than 1 in 4,000. They are confident this is a real, physical oval orbit.

The Big Picture

This paper is a game-changer because it proves that not all black hole and neutron star couples are born in quiet, isolated nurseries. Some are the result of chaotic, crowded cosmic brawls.

It's like finding out that while most people meet their partners at a quiet coffee shop, some actually meet in the middle of a mosh pit at a rock concert. Both happen, but for a long time, we only knew about the coffee shop. Now, thanks to this "oval dance," we know the mosh pit is real too.

In short: We found a cosmic couple that was swinging in an oval, proving they met in a chaotic crowd, and our new tools are finally good enough to see the shape of their dance.

Technical Summary

1. Problem Statement

The formation channels of neutron star–black hole (NSBH) binaries are a major open question in astrophysics. Two primary scenarios exist:

• Isolated Binary Evolution: Massive binary stars evolve in isolation, interacting via mass transfer and common envelope phases. Standard models predict that gravitational wave (GW) emission will rapidly circularize the orbit long before the binary enters the sensitivity band of ground-based detectors.
• Dynamical Assembly: Binaries form or are hardened through dynamical interactions in dense stellar environments (e.g., globular clusters, nuclear star clusters) or hierarchical triple systems. These channels can produce binaries with significant orbital eccentricity ($e$) and misaligned spins even at small separations.

While spin precession and orbital eccentricity are powerful tracers to distinguish these channels, neither has been confidently detected in NSBH mergers to date. Previous analyses of events like GW200105 relied on waveform models that assumed circular orbits, potentially introducing systematic biases in parameter estimation (e.g., mass and spin measurements) and missing evidence for dynamical formation.

2. Methodology

The authors re-analyzed the GW event GW200105 (detected on January 5, 2020) using a novel theoretical framework designed to simultaneously model spin-induced orbital precession and orbital eccentricity.

• Waveform Models:
  • Primary Model: pyEFPE (G. Morras et al. 2025a), a post-Newtonian (PN) inspiral-only waveform model that incorporates both spin precession and eccentricity. It uses higher-order Fourier harmonics but excludes higher-order spherical harmonic modes (HOMs) and the merger-ringdown phase.
  • Comparative Models: The analysis was cross-checked against TaylorF2Ecc (eccentric, aligned-spin), IMRPhenomXP (precessing, circular), and IMRPhenomXPHM (precessing, circular with HOMs).
• Inference Framework:
  • Bayesian Inference: Performed using the nested sampling algorithm Dynesty within the Bilby framework.
  • Data: Utilized public strain data from LIGO Livingston and Virgo. The analysis focused on the inspiral phase, with a high-frequency cutoff ($f_{high}$) of 340 Hz (conservative MECO bound) because the pyEFPE model does not cover the merger-ringdown.
  • Priors: A uniform prior was adopted for eccentricity at 20 Hz ($e_{20} \in [0, 0.4]$). Alternative priors (log-uniform, inverse linear/quadratic) were tested for robustness.
• Robustness Checks:
  • Waveform Systematics: Verified that the eccentricity signal is not an artifact of missing HOMs or the merger phase by comparing with IMRPhenomXPHM and injecting synthetic signals.
  • Noise Artifacts: Conducted hierarchical empirical Bayesian analysis on 20 realizations of Gaussian noise to quantify the probability that noise alone could mimic the observed eccentricity.
  • Prior Sensitivity: Tested various prior choices to ensure the result was data-driven rather than prior-dominated.

3. Key Contributions

• First Joint Measurement: This is the first study to perform a joint Bayesian inference of orbital eccentricity and spin precession in an NSBH system.
• Novel Waveform Application: It represents the first application of the pyEFPE model to real GW data, enabling the simultaneous extraction of precession and eccentricity parameters.
• Systematic Bias Correction: The study demonstrates that neglecting eccentricity in NSBH analyses introduces significant biases in mass and chirp mass estimation, correcting previous interpretations of GW200105.

4. Key Results

• Orbital Eccentricity:
  • The analysis infers a median orbital eccentricity at a GW frequency of 20 Hz of $e_{20} = 0.145^{+0.007}_{-0.063}$.
  • Eccentricities smaller than 0.028 are ruled out with 99.5% confidence.
  • The posterior distribution strongly disfavors a circular binary ($e=0$).
  • No other NSBH or BNS events in the analyzed sample (GW170817, GW190425, GW200115, etc.) showed statistically significant evidence for eccentricity.
• Mass Parameters:
  • Primary Mass ($m_1$): $11.5^{+0.6}_{-1.6} M_\odot$ (Black Hole).
  • Secondary Mass ($m_2$): $1.5^{+0.2}_{-0.1} M_\odot$ (Neutron Star).
  • Chirp Mass ($\mathcal{M}_c$): $3.33^{+0.09}_{-0.07} M_\odot$.
  • Significance: These values differ from the original LVK analysis (which assumed circularity), which reported $m_1 \approx 8.9 M_\odot$ and $m_2 \approx 1.9 M_\odot$. The authors attribute this discrepancy to the bias introduced by neglecting eccentricity, which leads to an overestimation of the chirp mass. The new secondary mass is more consistent with observed galactic millisecond pulsar masses.
• Spin Properties:
  • Effective Spin ($\chi_{eff}$): $0.005^{+0.071}_{-0.095}$ (consistent with zero).
  • Effective Precession Spin ($\chi_p$): $< 0.19$ (95% upper limit).
  • The binary is consistent with very low, at most mildly precessing spins.
  • Correlation: No correlation was found between eccentricity and precession in this event.
• Periastron Advance: The inferred rate of periastron advance is $\dot{\omega} \approx 270^\circ \text{s}^{-1}$, roughly eight orders of magnitude higher than that of the eccentric double neutron star PSR J1949+2052.

5. Significance and Astrophysical Implications

• Evidence for Dynamical Formation: The detection of non-zero eccentricity ($e_{20} \sim 0.145$) in a system that has not yet merged provides compelling evidence that GW200105 did not form via standard isolated binary evolution. Instead, it likely formed through dynamical interactions (e.g., in dense star clusters) or hierarchical triple systems (via the Zeipel-Lidov-Kozai mechanism), which can sustain eccentricity against gravitational wave circularization.
• Formation Channels: The result suggests that a fraction of NSBH binaries may retain eccentricity even at small separations, challenging the assumption that all NSBHs are circular by the time they enter the detector band.
• Future Observations: The authors project that if this eccentricity evolves solely via GW emission, the system would have been highly eccentric in the mHz band (LISA). Future observations with next-generation ground-based detectors and space-based observatories (LISA) will be crucial for characterizing the population of eccentric NSBHs and disentangling the relative merger rates of different formation channels.
• Methodological Impact: The study highlights the necessity of including eccentricity in waveform models to avoid systematic biases in mass and spin measurements, which are critical for testing fundamental physics and nuclear equations of state.

Conclusion

The paper reports the first confident measurement of orbital eccentricity in an NSBH merger (GW200105). By utilizing a waveform model that accounts for both precession and eccentricity, the authors infer a significant eccentricity ($e_{20} \approx 0.145$) and a more unequal mass ratio than previously thought. This finding provides strong evidence for a dynamical formation channel for this specific event and underscores the importance of eccentric waveform models for accurate astrophysical inference in the era of gravitational wave astronomy.