Systematic bias due to eccentricity in parameter estimation for merging binary neutron stars : Spinning case

This study extends previous research on eccentricity-induced systematic biases in binary neutron star parameter estimation to include spin effects, revealing that while biases in chirp mass, mass ratio, and effective spin follow predictable quadratic trends with eccentricity, biases in effective tidal deformability exhibit a complex dependence on mass and spin parameters, potentially impacting neutron star property inference.

The Gist

Imagine you are a detective trying to solve a mystery by listening to a faint, distorted recording of a car crash. Your goal is to figure out exactly how heavy the cars were, how fast they were spinning, and what kind of metal they were made of.

This paper is about a specific problem with that detective work: what happens if you ignore the fact that the cars were swerving (eccentricity) and assume they were driving in perfect circles?

Here is the breakdown of the research in simple terms:

The Setup: The "Perfect Circle" Assumption

For a long time, scientists studying gravitational waves (the "sound" of colliding stars) assumed that when two neutron stars crash, they are moving in perfect, smooth circles. It's a convenient assumption, like assuming a car is driving straight down a highway.

However, in reality, some of these stars might be swerving or moving in slightly oval paths (eccentricity) before they crash. This paper asks: If we pretend they are driving in perfect circles when they are actually swerving, how much do we mess up our calculations?

The Experiment: The "Swerving" Simulation

The researchers, Eunjung Lee, Chang-Hwan Lee, and Hee-Suk Cho, ran thousands of computer simulations. They created "fake" signals of colliding neutron stars that had a little bit of swerving (eccentricity).

Then, they tried to analyze these signals using the standard tools that ignore the swerving. They wanted to see how much the "detective" would get the details wrong.

The Findings: The "Distorted Lens"

They found that even a tiny bit of swerving (which is hard to detect) acts like a distorted lens on a camera. It makes the picture look very different from reality. Here is what went wrong:

1. The Weight Mistake (Mass):

   • The Analogy: Imagine you see a heavy truck swerving. If you assume it's driving straight, you might think it's actually a tiny sports car that is just moving very fast, or a giant truck that is moving very slowly.
   • The Result: The scientists found that if they ignored the swerving, they could estimate a pair of normal-sized neutron stars (about the size of Earth but as heavy as the Sun) to be either much lighter or much heavier than they actually were. In some extreme cases, they might even think one of the stars is so heavy it shouldn't exist (falling into the "mass gap" between stars and black holes).

2. The Spin Mistake:

   • The Analogy: If a spinning top is wobbling, and you assume it's spinning perfectly straight, you might guess it's spinning faster or slower than it really is.
   • The Result: The estimated spin of the stars was also thrown off, though not as wildly as the mass.

3. The "Recipe" Mistake (Equation of State):

   • The Analogy: Neutron stars are made of incredibly dense matter. Scientists have different "recipes" (Equations of State) for how this matter behaves—some are "soft" (squishy) and some are "stiff" (hard).
   • The Result: This is the most dangerous error. Because the mass and spin estimates were wrong, the scientists' calculation of the star's "squishiness" was also wrong.
   • The Consequence: They could look at a star made of "Recipe A" and, because of the swerving, conclude it was actually made of "Recipe B." This means we could be completely wrong about the fundamental laws of physics governing these stars.

The "Magic Number"

The researchers found that you don't need a lot of swerving to cause these big mistakes. Even if the stars are only swerving a tiny bit (about 2% deviation from a perfect circle), the errors in our calculations become larger than the natural uncertainty of our instruments.

The Takeaway

The paper warns us that as we get better at listening to the universe (with future, super-sensitive detectors), we must stop assuming everything moves in perfect circles.

If we don't account for the "swerving" (eccentricity), we risk:

• Thinking normal stars are monsters or ghosts.
• Getting the "recipe" of the universe wrong.
• Misidentifying the nature of the objects we are studying.

In short: You can't solve the mystery of the crashing stars if you ignore the fact that they were driving in a circle of chaos. We need to update our "detective tools" to account for the wobble, or we'll keep getting the wrong answer.

Technical Summary

Here is a detailed technical summary of the paper "Systematic bias due to eccentricity in parameter estimation for merging binary neutron stars: Spinning case" by Eunjung Lee, Chang-Hwan Lee, and Hee-Suk Cho.

1. Problem Statement

Gravitational wave (GW) parameter estimation for merging binary neutron stars (BNS) typically assumes quasi-circular orbits, ignoring orbital eccentricity ($e$). While most BNS systems are expected to circularize before entering the detector's frequency band, dynamically formed binaries in dense stellar clusters may retain significant eccentricity ($e \gtrsim 0.01$) at 10 Hz.

• The Core Issue: If an eccentric signal is analyzed using non-eccentric waveform templates, it introduces systematic biases in the estimated parameters (mass, spin, tidal deformability).
• The Gap: Previous work by the authors (and others) addressed this for non-spinning BNS systems. However, realistic BNS systems possess spin. This paper investigates how the inclusion of spin parameters alters the nature and magnitude of these systematic biases, specifically focusing on how they affect the inference of the Neutron Star Equation of State (EoS).

2. Methodology

The authors employ a hybrid approach combining analytic approximations with rigorous numerical Bayesian inference.

• Waveform Model: They utilize the TaylorF2 post-Newtonian (PN) waveform model, which includes terms for point-particle dynamics, spin, tidal effects, and eccentricity (up to 3PN order for eccentricity).
• Bias Calculation (Analytic):
  • They use the Fisher Matrix (FM) approach to estimate measurement errors.
  • They employ the Cutler-Vallisneri (FCV) method to analytically calculate the systematic bias ($\Delta \theta$) introduced by model mismatch (using a non-eccentric template for an eccentric signal).
  • The bias formula used is: $\Delta \theta_i = 4A^2 \Sigma^P_{ij} \int \frac{f^{-7/3}}{S_n} (\Psi_{true} - \Psi_{approx}) \partial_j \Psi_{approx} df$.
• Validation (Numerical):
  • To verify the reliability of the FM/FCV methods, they perform Bayesian parameter estimation using the BILBY and DYNESTY libraries.
  • They compare the posterior distributions from Bayesian inference against the Gaussian approximations derived from the Fisher Matrix.
• Simulation Setup:
  • Detector: Advanced LIGO (aLIGO) O4 high-sensitivity curve.
  • Source Population: $10^4$ Monte Carlo BNS signals randomly distributed in the parameter space:
    • Masses: $1 M_\odot \le m_{1,2} \le 2 M_\odot$($m_2 \le m_1$).
    • Effective Spin: $-0.2 \le \chi_{eff} \le 0.2$.
    • Eccentricity: $0 \le e_0 \le 0.024$ (at 10 Hz).
    • Tidal Deformability: Determined by the APR4 EoS model.
  • Signal-to-Noise Ratio (SNR): Analyses assume a high SNR of $\rho = 300$ to ensure the validity of the high-SNR approximation used in the FCV method.

3. Key Contributions

1. Extension to Spinning Systems: This is the first comprehensive study of eccentricity-induced bias in BNS parameter estimation that explicitly includes spin parameters ($\chi_{eff}$).
2. Reversal of Bias Direction: A critical finding is that the inclusion of spin reverses the direction of the bias for certain parameters compared to the non-spinning case. Specifically, the bias direction for the effective tidal deformability ($\tilde{\lambda}$) flips from positive to negative.
3. EoS Inference Failure: The paper demonstrates concretely how these biases can lead to false predictions of the Neutron Star Equation of State, causing an observer to prefer a softer or stiffer EoS model than the true one.
4. Mass Range Distortion: The study shows that eccentricity can cause the estimated component masses to deviate significantly from the typical BNS mass range ($1-2 M_\odot$), potentially pushing estimates into the "mass gap" or below typical neutron star masses.

4. Key Results

A. Validation of Methods

The FCV method was validated against Bayesian results. For high SNRs ($\rho=300$), the Gaussian approximation of the posterior distribution from the Fisher Matrix closely matches the Bayesian results, confirming the reliability of the analytic bias calculation for small eccentricities ($e_0 < 0.024$).

B. Bias Trends for $M_c$, $\eta$, and $\chi_{eff}$

• Distribution: The fractional biases ($\Delta \theta / \sigma_\theta$) for Chirp Mass ($M_c$), Symmetric Mass Ratio ($\eta$), and Effective Spin ($\chi_{eff}$) form narrow bands.
• Dependence: These biases show a weak dependence on the specific mass and spin values of the source.
• Scaling: The magnitude of the bias increases quadratically with eccentricity ($e_0$).
  • $\Delta M_c / \sigma_{M_c} \propto +e_0^2$ (Positive bias)
  • $\Delta \eta / \sigma_{\eta} \propto -e_0^2$ (Negative bias)
  • $\Delta \chi_{eff} / \sigma_{\chi_{eff}} \propto +e_0^2$ (Positive bias)

C. Bias Trends for Tidal Deformability ($\tilde{\lambda}$)

• Wide Distribution: Unlike the other parameters, the bias in $\tilde{\lambda}$ is widely distributed and strongly dependent on the specific mass and spin configuration of the binary.
• Spin-Induced Reversal:
  • Non-spinning case (Previous work): Bias was positive (recovered $\tilde{\lambda}$ > true $\tilde{\lambda}$), favoring a stiffer EoS.
  • Spinning case (Current work): Bias becomes negative (recovered $\tilde{\lambda}$ < true $\tilde{\lambda}$), favoring a softer EoS.
  • Significance: This reversal implies that ignoring spin when analyzing eccentric signals could lead to incorrect conclusions about the stiffness of nuclear matter.

D. Impact on Mass Inference

• Asymmetry: Eccentricity biases tend to make the recovered component masses more asymmetric than the true values.
• Mass Gap/Range Violation:
  • A signal with true masses $(1.2 M_\odot, 1.2 M_\odot)$ can be recovered as $(1.9 M_\odot, 0.8 M_\odot)$.
  • A signal with true masses $(1.8 M_\odot, 1.8 M_\odot)$ can be recovered as $(2.8 M_\odot, 1.2 M_\odot)$.
  • In extreme cases (e.g., $m_1=m_2=2 M_\odot, \chi_{eff}=0.2, e_0=0.02$), the recovered mass can enter the mass gap ($>2.5 M_\odot$), leading to misidentification of the object type.

E. Impact on EoS Inference (Concrete Examples)

The authors provide two specific Bayesian examples demonstrating EoS misidentification:

1. APR4 vs. WFF1: An eccentric signal with true EoS APR4 (and masses $\sim 1.37 M_\odot$) is analyzed with non-eccentric templates. The recovered masses shift to $\sim (1.87, 1.01) M_\odot$, and the recovered tidal parameter favors the WFF1 EoS model over the true APR4 model.
2. MPA1 vs. APR4: For a $(1.4 M_\odot, 1.4 M_\odot)$ signal following the MPA1 EoS, an eccentricity as low as $e_0 \approx 0.0041$ causes the parameter estimation to prefer the APR4 model (a softer EoS) over the true MPA1 model.

5. Significance and Conclusion

• Critical for 3G Detectors: With the upcoming Einstein Telescope and Cosmic Explorer, the ability to detect signals with higher SNR and lower frequency sensitivity will make even tiny eccentricities ($e_0 \sim 10^{-3}$) detectable. If unmodeled, these will dominate statistical errors and lead to systematic errors in EoS constraints.
• Need for Eccentric Models: The results strongly suggest that for high-precision BNS science, especially regarding the EoS, eccentric waveform models must be included in parameter estimation pipelines, even for signals that appear nearly circular.
• Diagnostic Tool: The specific bias trends (e.g., the relationship between $\Delta \eta$ and $e_0$) can be used as a diagnostic tool. If a recovered signal shows an unexpectedly large asymmetry in masses or a specific bias in $\tilde{\lambda}$, it may indicate the presence of hidden eccentricity, allowing researchers to place upper limits on the source's eccentricity.

In summary, this paper establishes that spin and eccentricity are coupled in their impact on parameter estimation, and neglecting either can lead to severe misinterpretations of neutron star properties, including the fundamental nature of dense nuclear matter.