Fixing the center-of-mass frame of numerical relativity waveforms using the post-Newtonian center-of-mass charge

This paper improves the robustness of center-of-mass frame fixing for numerical relativity waveforms by incorporating post-Newtonian boosted center-of-mass charge calculations to capture physical out-spiraling oscillations, thereby significantly reducing the sensitivity of fit parameters to the choice of fitting window and integrating the method into the `scri` Python package.

The Gist

The Big Picture: Tuning the "Gravity Radio"

Imagine the universe is a giant radio station broadcasting gravitational waves (ripples in space-time) from colliding black holes. Scientists use detectors like LIGO to "listen" to these signals. To understand what they hear, they need a perfect "map" or "script" of what the signal should look like.

This paper is about fixing the static and distortion on that radio signal so the map is accurate.

The Problem: The "Drifting Camera"

When scientists run supercomputer simulations of two black holes crashing into each other, they get a waveform (the signal). However, because of how the computer code is set up, the simulation often starts with a hidden flaw: the camera is moving.

Think of it like filming a dance with a shaky, drifting camera.

• The Drift: The camera slowly drifts away from the dancers (the black holes).
• The Wobble: The camera also wobbles slightly as the dancers spin.

Because the camera is moving, the video looks weird. The main dancers (the dominant signal) seem to leak energy into the background dancers (faint, higher-order ripples). In the real world, this "leakage" is just an illusion caused by the moving camera, not a real physical effect. If scientists use these "drifting" videos to build their maps, their maps will be wrong, and they might misinterpret the real signals from space.

The Old Solution: Drawing a Straight Line

Previously, scientists tried to fix this "moving camera" problem by looking at the center of the dance floor. They noticed the center was drifting in a straight line. So, they drew a straight line through the data and said, "Okay, let's just subtract that straight line to steady the camera."

The Flaw: This was like trying to fix a bumpy ride by only accounting for the straight road, ignoring the actual bumps. The center of mass of black holes doesn't just drift in a straight line; it also spirals and oscillates (wiggles) as they orbit each other. The old "straight line" method missed these wiggles. This made the fix very sensitive: if you picked a slightly different chunk of time to analyze, your "straight line" would change, and your fix would be different. It was unstable.

The New Solution: The "Physics GPS"

The authors of this paper came up with a better way. Instead of just drawing a straight line, they used Post-Newtonian (PN) theory.

Think of PN theory as a physics-based GPS. It knows the exact laws of gravity and can predict exactly how the center of the dance floor should move, including the straight drift and the spiraling wiggles.

1. The Prediction: They calculated a mathematical formula that predicts the "perfect" path of the center of mass, including all the wiggles.
2. The Comparison: They compared the messy, drifting computer simulation data against this perfect GPS prediction.
3. The Fix: They adjusted the simulation's "camera" (the frame) until the messy data matched the perfect prediction.

Why This is a Game Changer

The paper tested this new method against the old "straight line" method using 20 different black hole collision simulations.

• Stability: The old method was like trying to balance a broom on your finger; if you moved your hand slightly (changed the time window), the broom fell. The new method is like a gyroscopic stabilizer. Even if you change the time window you are looking at, the result stays almost exactly the same.
• The Numbers: The new method made the results 25 times more stable for the "boost" (the speed of the drift) and 20 times more stable for the "translation" (the position of the drift).

The "Sweet Spot"

The researchers also found the best place to apply this fix. Imagine the black hole collision is a song.

• The Start: The beginning of the song has some "static" (junk radiation) from the computer starting up.
• The End: The end of the song is the crash, where things get chaotic.
• The Middle: The middle of the song is the steady, rhythmic orbit.

They found that applying their fix in the middle of the song (the center of the inspiral) works best. It avoids the static at the start and the chaos at the end, giving the cleanest, most accurate result.

The Bottom Line

This paper provides a super-stable, physics-based tool to clean up gravitational wave simulations. By using a detailed mathematical model (the GPS) instead of a simple guess (the straight line), scientists can now create much more accurate maps of black hole collisions. This means when the real gravitational waves hit Earth's detectors, we will be able to decode them with much higher precision, helping us understand the universe better.

In short: They stopped guessing the camera's movement and started using a physics-based GPS to keep the camera perfectly steady.

Technical Summary

1. Problem Statement

Numerical Relativity (NR) simulations of binary black hole mergers produce gravitational waveforms that are extracted at future null infinity ($\mathscr{I}^+$). However, these waveforms are initially defined in an arbitrary Bondi–van der Burg–Metzner–Sachs (BMS) frame due to gauge choices made during the simulation setup.

• The Issue: An arbitrary BMS frame introduces unphysical gauge artifacts. Specifically, if the center of mass (CoM) of the binary is not fixed, the waveform exhibits spurious linear drifts and oscillations. This causes power from the dominant $(2, \pm 2)$ modes to leak into subdominant higher-order modes (mode mixing), degrading the accuracy of waveform models used for gravitational wave detection and parameter estimation.
• Limitations of Previous Methods: Previous approaches to fix the CoM frame (e.g., in the SXS catalog) relied on fitting the numerically computed CoM charge vector $\vec{G}(u)$ to a linear function of time. This method only accounted for the linear drift caused by an initial boost but ignored the physical, oscillatory "out-spiraling" motion of the CoM driven by linear momentum conservation. Consequently, the resulting boost and translation parameters were highly sensitive to the choice of the fitting window's size and location.

2. Methodology

The authors propose an improved frame-fixing procedure that incorporates analytical Post-Newtonian (PN) theory to model the full behavior of the CoM charge.

• Theoretical Derivation:

  • The authors derive the PN expression for the CoM charge $\vec{G}$ for quasicircular, non-precessing binaries.
  • They utilize the CoM balance law (Compère et al.) and PN equations of motion to express the time derivative of the CoM charge ($d\vec{G}/dt$) in terms of source multipole moments.
  • They derive the leading-order PN expression for $\vec{G}$, which includes both a linear drift term and a physical oscillatory term proportional to the orbital phase.
  • They derive the transformation rules for $\vec{G}$ under small boosts and translations to construct a fitting function that relates the numerical data to the desired "PNBMS" frame (Post-Newtonian BMS frame).

• The Fitting Function:
  Instead of a simple linear fit, the new method fits the numerical CoM charge $\vec{G}_{num}$ to an analytical model:
  $$\vec{G}' = \frac{1}{P^0} \left[ (\alpha_1 \hat{\lambda} + \alpha_2 \hat{n})|\vec{G}| - \vec{\beta} \times \vec{J} - u P^0 \vec{\beta} \right] + \vec{\Delta}$$
  Where:

  • $\vec{\beta}$ and $\vec{\Delta}$ are the boost and translation parameters to be determined.
  • $\hat{\lambda}$ and $\hat{n}$ are unit vectors defining the orbital plane.
  • $\alpha_1, \alpha_2$ are nuisance parameters accounting for amplitude/direction mismatches between PN and NR.
  • The term $u P^0 \vec{\beta}$ captures the linear drift, while the oscillatory terms capture the physical out-spiral.

• Implementation:

  • The method is applied to 20 quasicircular, non-precessing binary black hole simulations from the SXS catalog.
  • Waveforms are extracted using the Cauchy-characteristic evolution (CCE) method (via the SpECTRE code), which provides the necessary Weyl scalars ($\Psi_0-\Psi_4$) to compute the full BMS charges.
  • The authors perform a sensitivity analysis by varying the fitting window size (from $\sim 500M$ to $4000M$) and its location (start-fixed, center-fixed, end-fixed).

3. Key Contributions

1. Analytical Modeling of CoM Oscillations: The paper provides the first analytical PN derivation of the CoM charge that captures both the linear drift and the physical oscillatory behavior, moving beyond the purely linear approximation used in prior work.
2. Robust Frame Fixing: The authors demonstrate that using the PN-based fitting function significantly reduces the sensitivity of the derived boost and translation parameters to the choice of the fitting window.
3. Software Integration: The new method has been implemented into the scri Python package, a standard tool for analyzing NR waveforms, making it accessible for the broader community.
4. Quantitative Improvement: The study quantifies the improvement in robustness, showing that the new method yields stable parameters even with smaller fitting windows.

4. Results

• Robustness Metrics: The authors quantified the improvement by calculating the ratio of variances of the fit parameters (boost $\vec{\beta}$ and translation $\vec{\Delta}$) when the window size is varied.
  • Boost Vector ($\vec{\beta}$): The new method improved robustness by a factor of $\sim 25$ compared to the linear fit when the window was centered.
  • Translation Vector ($\vec{\Delta}$): The improvement was a factor of $\sim 20$ under the same conditions.
• Window Sensitivity:
  • The center-fixed window placement yielded the best results, as it minimizes the impact of "junk radiation" (early inspiral artifacts) and higher-order PN effects (late inspiral).
  • The authors recommend a minimum window size of $1500M$ (where $M$ is the total mass) for optimal stability.
• Visual Confirmation: Plots of the CoM charge show that the PN analytical curve (dashed orange) fits the numerical data (solid blue) significantly better than a linear fit, accurately capturing the oscillatory "out-spiral" pattern.

5. Significance

• Improved Waveform Accuracy: By fixing the BMS frame more robustly, the resulting waveforms have reduced mode mixing. This is critical for constructing accurate surrogate models and phenomenological templates (like Phenom and EOB) used by the LIGO-Virgo-KAGRA (LVK) collaboration.
• Standardization: The method provides a more rigorous, physics-based standard for defining the "rest frame" of numerical simulations, reducing arbitrary gauge choices that could introduce systematic errors in parameter estimation.
• Future Directions: While this work is limited to quasicircular, non-precessing systems, it lays the groundwork for fixing the entire BMS gauge (including supertranslations and rotations) for generic systems, including eccentric and precessing binaries, which is essential for future gravitational wave detectors like LISA and the Einstein Telescope.

In summary, this paper replaces a heuristic linear approximation with a physically motivated Post-Newtonian model to fix the center-of-mass frame of numerical relativity waveforms, resulting in a substantial increase in the reliability and precision of gravitational waveform models.