PRECESSION 2.1: black-hole binary spin precession on eccentric orbits

This paper introduces version 2.1 of the public Python code PRECESSION, which extends its capabilities to model the post-Newtonian dynamics of precessing black hole binaries on eccentric orbits through a semi-automatic adaptation method and new averaged evolutionary equations.

The Gist

Imagine two massive black holes dancing around each other in the vast darkness of space. For a long time, scientists have had a very good way to predict how this dance moves, but they mostly assumed the dancers were moving in perfect, smooth circles.

This paper introduces an upgrade to a computer program called precession (specifically version 2.1). Think of this program as a "dance simulator" for black holes. The new update allows the simulator to handle a much more chaotic and realistic scenario: elliptical orbits, where the black holes swoop in close and then swing far out, like a comet around the sun, rather than staying in a perfect circle.

Here is a breakdown of what the authors did, using some everyday analogies:

1. The "Magic Translator" (The Python Decorator)

The biggest challenge was that the old code was built for circular orbits, and the new code needs to handle oval (eccentric) ones. Rewriting the whole thing from scratch would be like rebuilding a car engine just to add a new seat.

Instead, the authors built a clever "translator" tool (a Python decorator named eccentricize).

• The Analogy: Imagine you have a recipe for a perfect round cake. Now you want to bake an oval cake. Instead of rewriting the whole recipe, you use a magic kitchen gadget that automatically adjusts the measurements. If the recipe says "use a 10-inch round pan," the gadget says, "Okay, for an oval shape, we need a 10-inch semi-latus rectum (a specific oval measurement) and we need to bake it for a slightly different amount of time."
• The Result: The code can now take the old, trusted logic for circular orbits and automatically adapt it to work for oval orbits without breaking anything.

2. The "Wobbly Top" (Spin Precession)

Black holes aren't just heavy balls; they spin like tops. As they orbit each other, their spins wobble and tilt, a phenomenon called precession.

• The Analogy: Think of a spinning top on a table. If the table is perfectly flat (a circular orbit), the top spins steadily. But if the table is tilted and bumpy (an elliptical orbit), the top wobbles wildly.
• The New Feature: The update tracks how these "wobbles" change when the orbit is oval. It also tracks a new angle called the periastron angle (the point where the black holes get closest). It's like tracking not just how the top spins, but also where the "closest point" of its wobble is moving around the table.

3. The "Radio Station" (Gravitational Waves)

When black holes dance, they send out ripples in space-time called gravitational waves. In a perfect circle, these waves are like a steady hum at one specific pitch.

• The Analogy: In an oval orbit, the dance speeds up and slows down. This is like a siren on a police car: as it speeds up, the pitch goes up; as it slows down, the pitch goes down. The sound isn't just one note anymore; it's a complex chord with many different frequencies (harmonics).
• The New Feature: The code now knows how to calculate this complex "chord." It can figure out exactly what frequency the gravitational waves will be, even when the black holes are speeding up and slowing down in their oval path. It also fixed a small error in the old code regarding how the black holes' own spin affects the "shape" of the space around them.

4. Why Does This Matter?

Detecting these gravitational waves is like listening to the universe's history.

• The Detective Work: If we see a black hole dance in a perfect circle, it likely formed in a quiet, calm environment. If we see it dancing in a wild, oval path, it might have formed in a crowded, chaotic star cluster where black holes bumped into each other.
• The Goal: By updating the simulator to handle these wild, oval dances, scientists can better understand how these black hole pairs were born in the first place. It helps us distinguish between different "origins stories" for these cosmic giants.

Summary

In short, the authors took a powerful tool for simulating black hole dances and gave it a "wild card" mode. They created a smart system that automatically translates "circular orbit" rules into "oval orbit" rules, added new math to track the complex wobbles and sound frequencies of oval dances, and fixed a few tiny bugs. This allows scientists to listen to the universe more clearly and understand the chaotic history of black holes.

Technical Summary

1. Problem Statement

Gravitational wave (GW) detections of black hole (BH) binaries provide crucial data regarding the properties of these systems. While some parameters (masses, spin magnitudes) remain constant, others (spin orientations, orbital eccentricity) evolve significantly over time. To accurately infer the formation channels of BH binaries, it is necessary to connect measured GW parameters back to their values at formation.

Previous versions of the numerical code precession were designed exclusively for quasi-circular orbits. However, many astrophysical formation channels (e.g., dynamical capture in dense clusters) produce binaries with significant orbital eccentricity when they enter the GW detector band. Existing circular-only models cannot accurately track the evolution of spin vectors and orbital parameters for these eccentric systems, creating a gap in population synthesis and parameter estimation capabilities.

2. Methodology

The authors present version 2.1 of the precession Python module, extending its capabilities to handle eccentric binaries using a multi-timescale post-Newtonian (PN) formalism. The methodology involves two primary technical approaches:

A. Mapping to Effective Circular Sources

To leverage existing numerical infrastructure, the authors map the dynamics of an eccentric binary onto an effective circular source on the precessional timescale. This mapping replaces the orbital separation $r$ and time $t$ with eccentricity-dependent variables:

• Semi-latus rectum: $a \to a(1-e^2)$
• Rescaled time: $t \to t(1-e^2)^{3/2}$

Implementation: A custom Python decorator named eccentricize was developed. This utility dynamically wraps existing functions, modifying their user-facing interface to accept semi-major axis ($a$) and eccentricity ($e$) instead of separation ($r$). Internally, the functions compute $r$ from $a$ and $e$ using the mapping equations, preserving the original logic of the circular code while extending compatibility.

B. Specialized Treatment for Eccentric-Only Functions

Certain functions require manual modification because they have no circular counterpart or involve complex non-linearities:

• Evolutionary Equations: The code implements 2PN spin-precession equations coupled with 3PN evolutionary equations for the semi-major axis and eccentricity (including spin terms up to 2PN).
• New Angular Variables: Three new variables are introduced to track the system's geometry:
  • $\Psi_1, \Psi_2$: Angles between spin vector projections and the periastron line.
  • $\delta\lambda$: The periastron angle (precession of the orbit).
• Regularization: A new regularized version of the Peters' equation for $da/de$ is derived to ensure numerical stability as eccentricity approaches zero ($e \to 0$).
• GW Frequency Conversion: The relationship between orbital separation and GW emission frequency ($f_{GW}$) is revised. Unlike circular binaries (dominated by the second harmonic), eccentric binaries emit across multiple harmonics. The code now:
  1. Calculates orbital frequency $\omega$ from $a$ and $e$.
  2. Maps $\omega$ to $f_{GW}$ using fitting formulae to identify the dominant harmonic (either the second harmonic or the one with maximum power).
  3. Analytically inverts the frequency-separation relation to include spin-induced quadrupole terms (entering at 2PN order), correcting previous circular-only approximations.

3. Key Contributions

1. Code Release (v2.1): Public release of precession with full support for eccentric BH binaries, available via pip and GitHub.
2. Semi-Automatic Extension: Introduction of the eccentricize decorator, allowing for the efficient adaptation of circular-orbit functions to eccentric cases without rewriting core logic.
3. High-Precision Evolutionary Equations: Implementation of 3PN orbital evolution and 2PN spin precession specifically for eccentric orbits, including the evolution of the periastron angle $\delta\lambda$.
4. Improved Frequency-Separation Mapping:
   • Derivation of analytical inversion formulas (Eq. 2–5) for converting between semi-major axis and GW frequency, including mass ratio ($q$), spin magnitudes ($\chi$), and spin orientations ($\theta, \Delta\Phi$) dependencies.
   • Correction of circular-case formulas to include the previously missing spin-induced quadrupole term.
5. Numerical Stability: Development of regularized equations for the $e \to 0$ limit to prevent numerical singularities during integration.

4. Results

• Validation: The authors present simulations of BH binaries evolving from $a = 100M$ to $a = 10M$.
  • Comparison: They compare binaries with initial eccentricity $e_0 = 0.5$ against $e_0 = 0$.
  • Observation: The spin tilt angles ($\theta_1, \theta_2$) and periastron angle ($\delta\lambda$) evolve differently depending on eccentricity.
  • Continuity Check: The authors note that setting $e=0$ in the eccentric equations does not perfectly coincide with the circular evolution due to the non-zero derivative of eccentricity ($\dot{e} \neq 0$). This is a known feature of orbit-averaged equations that encodes the influence of spins on eccentricity.
• Performance: The code successfully handles the complex interplay between spin precession and orbital eccentricity, providing a seamless transition between circular and eccentric regimes.

5. Significance

• Formation Channel Discrimination: By accurately modeling the evolution of eccentric binaries, this tool allows researchers to better distinguish between different BH formation channels (e.g., isolated binary evolution vs. dynamical capture), as eccentricity is a key discriminant.
• Population Synthesis: The code enables the evolution of population synthesis predictions from formation to the GW detection regime, bridging the gap between theoretical models and observational data.
• Data Analysis Readiness: As GW detectors (like LIGO/Virgo/KAGRA and future space-based detectors like LISA) become more sensitive, the ability to model eccentric precession is critical for parameter estimation and waveform modeling.
• Community Resource: By making the code public and modular (separating eccentric functionality into precession.eccentricity), the authors facilitate broader adoption and further development within the gravitational wave community.

Note on Limitations: The current implementation focuses on the inspiral phase. Post-merger analytic fits (for remnant mass, spin, and recoil) remain restricted to quasi-circular binaries and have not yet been ported to the eccentric version.