The Universal Eccentricity Distribution for Dynamical Gravitational-Wave Merger Channels

1. Problem Statement

The detection of gravitational waves (GWs) from binary black hole (BBH) mergers has revealed a population of compact objects, yet the formation channels of these binaries remain difficult to disentangle. While mass and spin distributions often overlap between different formation scenarios (e.g., isolated binary evolution vs. dynamical formation in dense clusters), orbital eccentricity is considered a "smoking gun" for dynamical formation.

However, current observational analyses typically assume a uniform or log-uniform prior for eccentricity ($P(e) \propto 1/e$) when searching for eccentric sources. The authors argue that this assumption is not physically motivated. The central problem addressed is: What is the intrinsic eccentricity distribution of BBHs formed via dynamical mechanisms at the frequencies relevant for LIGO/Virgo/KAGRA (LVK) detectors ($f \geq 10$ Hz), and does this distribution depend on the specific astrophysical environment?

2. Methodology

The authors develop an analytical framework based on the separation of scales between the astrophysical formation environment and the GW-driven merger phase.

• The "Pinhole Regime" Concept:
  The authors argue that the pericentre distance ($r$) required for a BBH to merge with measurable eccentricity in the LVK band is orders of magnitude smaller than the characteristic scales of the formation environment (e.g., cluster size, binary separation). Consequently, the incoming black holes must pass through a tiny "pinhole" in phase space.

  • Because the "hole" is so small relative to the environment, the spatial distribution of incoming black holes is effectively stochastic and isotropic.
  • This implies that the distribution of impact parameters ($b$) is uniform over the area, leading to a specific distribution of pericentre distances.

• Analytical Derivation:

  1. Initial Conditions: Assuming an isotropic distribution of incoming black holes, the impact parameter distribution scales as $P(b) \propto b$. Using Keplerian relations ($b \propto \sqrt{r}$), the distribution of pericentre distances is derived as $P(r) \approx \text{constant}$.
  2. Mapping to Eccentricity: The authors map the constant $P(r)$ distribution to the eccentricity distribution $P(e)$ at a specific GW frequency ($f$). They utilize the post-Newtonian (PN) evolution of eccentricity and semi-major axis ($a$) as described by Peters (1964).
  3. Frequency Definitions: They derive distributions for two frequency definitions:
     • $f_p$: The peak GW frequency.
     • $f_T$: The orbital (Keplerian) frequency.
  4. Validation: The analytical results are compared against numerical $N$-body simulations of:
     • Chaotic binary-single interactions in dense stellar clusters.
     • Hierarchical triple systems undergoing Zeipel-Lidov-Kozai (ZLK) oscillations.

3. Key Contributions

• Universal Distribution Discovery: The paper demonstrates that all dynamical merger channels (binary-single, binary-binary, hierarchical triples, gas-assisted) converge to the same eccentricity distribution in the high-eccentricity limit ($e \gtrsim 0.05$) relevant for LVK.
• Analytical Solution: They provide a closed-form analytical solution for the probability density function (PDF) of eccentricity, which differs significantly from the commonly assumed log-flat prior.
• Scale Separation Argument: The work establishes that the "memory" of the specific formation environment is lost during the transition to the GW-dominated regime due to the scale separation, making the eccentricity distribution universal.
• Correction of Priors: The paper explicitly challenges the use of $P(e) \propto 1/e$ as a prior for eccentric searches, offering a physically motivated alternative.

4. Key Results

• The Universal PDF:
  In the high-eccentricity limit, the distribution of eccentricity $e_f$ at frequency $f$ follows a power law:
  $$P(e_f) \propto e_f^{-31/19}$$
  Or, expressed in terms of the logarithm of eccentricity:
  $$P(\log e_f) \propto e_f^{-12/19}$$
  This result holds regardless of the specific dynamical channel (e.g., cluster vs. triple) or the total mass of the system, provided finite-size effects are negligible.

• Comparison with Simulations:
  The analytical solution shows exceptional agreement with numerical simulations of:

  • Binary-single scatterings in globular clusters (tested for initial semi-major axes $a_0 = 0.1$ AU and $0.01$ AU).
  • Hierarchical triples (using data from Antonini et al.).
    The simulated distributions converge to the predicted $e^{-31/19}$ scaling at high eccentricities.

• Breakdown Limits:
  The distribution is valid down to a critical eccentricity $e_c$, below which the astrophysical environment (e.g., velocity dispersion $v_{env}$) influences the outcome. The critical eccentricity is estimated as:
  $$e_c \approx 0.01 \times \left(\frac{M}{50 M_\odot}\right)^{-19/18} \left(\frac{v_{env}}{50 \text{ km s}^{-1}}\right)^{19/21} \left(\frac{f}{10 \text{ Hz}}\right)^{-19/18}$$
  Below $e_c$, the distribution depends on the specific environment.

• Frequency Definition Nuance:
  The authors note a discrepancy when using the orbital frequency ($f_T$) versus the peak frequency ($f_p$). While $f_p$ yields a monotonically decreasing distribution, $f_T$ shows a significant upturn at high eccentricities. They argue that using $f_T$ for modest eccentricities maps to highly relativistic encounters ($r \lesssim 10 R_s$), suggesting that full Numerical Relativity (NR) simulations are needed to accurately describe the distribution in this specific limit.

5. Significance and Implications

• Optimized GW Searches: The derived universal distribution ($P(e) \propto e^{-31/19}$) should replace the log-flat prior in LVK template banks and search algorithms. This will improve the sensitivity to eccentric sources and reduce false negatives.
• Parameter Estimation: Using the correct prior will lead to more accurate parameter estimation for eccentric events, reducing biases in mass and spin measurements.
• Limits on Distinguishability: The paper highlights a fundamental limitation: in the high-eccentricity regime, different dynamical formation channels are indistinguishable based on eccentricity alone. The specific formation mechanism (e.g., cluster vs. triple) is encoded in:
  1. The fraction of sources that enter the "pinhole regime" (probability of eccentric merger).
  2. The critical eccentricity $e_c$ (set by the environment's velocity dispersion).
  3. Correlations with other parameters (mass, redshift, spin).
• Future Directions: To distinguish between formation channels, future work must focus on measuring lower eccentricities (where environmental signatures re-emerge) or correlating eccentricity with mass and spin distributions.

In conclusion, this paper provides a foundational theoretical framework for understanding eccentric BBH mergers, asserting that nature produces a universal eccentricity distribution for dynamical mergers in the high-eccentricity limit, driven by the stochastic nature of close encounters in dense environments.