High-order effective-one-body tidal interactions and gravitational scattering

This paper enhances the tidal sector of four effective-one-body (EOB) formalisms by incorporating next-to-next-to-leading post-Minkowskian order adiabatic and post-adiabatic quadrupolar tidal effects derived from scattering results, demonstrating improved agreement with numerical-relativity data for neutron star scattering and establishing a foundation for more accurate tidal models in both scattering and bound orbits.

The Gist

Imagine two massive, dense stars (like neutron stars) hurtling toward each other in the vast emptiness of space. They don't crash immediately; instead, they swing past one another like two skaters gliding on ice, whipping around each other due to gravity before flying apart again. This is called gravitational scattering.

This paper is about building a better "instruction manual" (a mathematical model) to predict exactly how these stars will behave during this high-speed dance, specifically focusing on how they squish and stretch each other due to gravity.

Here is the breakdown of the paper using simple analogies:

1. The Problem: The "Squishy" Stars

In the old days, scientists treated stars like perfect, hard billiard balls. But neutron stars are actually more like giant, super-dense marshmallows. When they get close, the gravity of one star pulls on the other, stretching it out like taffy. This "squishing" (called tidal effects) changes the path they take.

If we want to understand what happens when these stars collide (which creates the gravitational waves we detect on Earth), we need to know exactly how they squish. If our math is wrong, we can't figure out what the stars are made of (their "equation of state").

2. The Tool: The "Effective One-Body" (EOB)

Calculating how two complex objects interact is incredibly hard. It's like trying to predict the path of two tangled jump ropes.

• The Solution: Physicists use a clever trick called the Effective One-Body (EOB) formalism. Instead of tracking two messy ropes, they pretend the whole system is just one single, magical particle moving through a warped landscape.
• The Analogy: Imagine you are trying to predict how a car drives over a bumpy road. Instead of modeling every single wheel and suspension part, you model the car as a single point moving over a "bump map" that represents all the bumps combined. The EOB is that "bump map."

3. The Innovation: Using "Scattering" to Fix the Map

The authors realized that the best way to fix their "bump map" is to look at scattering events (the fly-by scenario mentioned earlier).

• The Old Way: They used approximations that worked okay for slow, circular orbits (like planets), but broke down when things got fast and messy.
• The New Way: They took the most advanced, high-precision calculations available (from a field called Post-Minkowskian gravity) regarding how stars fly past each other. They then "transcribed" this high-speed data into their EOB map.
• The Metaphor: Think of it like calibrating a GPS. If you only drive in a city (slow, circular orbits), your GPS might get confused on a highway (fast, scattering orbits). The authors took data from a Formula 1 race (high-speed scattering) and used it to update the GPS software so it works perfectly everywhere.

4. The Four "Flavors" of the Model

The paper didn't just make one map; they made four different versions of the EOB model to see which one worked best:

1. The "Post-Schwarzschild" version: A standard approach.
2. The "Lagrange" version: A newer, more flexible method that handles the math differently.
3. The "w-EOB" version: Another variation using a specific coordinate system.
4. The "PN" version: The traditional method used for decades.

They tested all four against Numerical Relativity (NR) data.

• What is NR? Imagine a supercomputer simulation that solves Einstein's equations exactly, like a high-definition video game rendering of the universe. It's the "truth" we compare our math against.

5. The Results: Who Won?

When they compared their new, improved maps to the supercomputer simulations:

• The Lagrange models (LEOB) were the clear winners. They matched the "truth" (the supercomputer) much better than the old models, especially when the stars got very close.
• The "Post-Adiabatic" Twist: The authors also looked at a subtle effect called "post-adiabatic" tides.
  • Analogy: If you stretch a rubber band slowly, it snaps back perfectly (adiabatic). If you stretch it very fast, it heats up and behaves differently (post-adiabatic).
  • They found that including this "heating up" effect was crucial for accuracy, though the exact size of this effect is still a bit of a mystery that needs more study.

6. Why Does This Matter?

This work is like upgrading the engine of a race car before the big race.

• Better Predictions: With these improved models, when we detect gravitational waves from colliding neutron stars in the future, we will be able to decode the signal much more accurately.
• Unlocking Secrets: By understanding the "squish" better, we can finally figure out the Equation of State of neutron stars. This tells us what matter looks like when it is crushed to the density of a nuclear atom. It's essentially telling us what the universe is made of at its most extreme.

Summary

The authors took the most advanced math available for high-speed star encounters and used it to upgrade the "instruction manual" (EOB) that physicists use to predict gravitational waves. They found that a specific new version of the manual (the Lagrange model) is much more accurate, bringing us one step closer to understanding the mysterious, super-dense matter inside neutron stars.

Technical Summary

1. Problem Statement

The accurate modeling of gravitational waves (GWs) from binary neutron star (NS) mergers is crucial for constraining the Equation of State (EOS) of dense nuclear matter. Tidal interactions significantly alter the dynamics of the final orbits before merger. While the Effective-One-Body (EOB) formalism is the standard analytical framework for generating GW templates, its tidal sector has historically relied on Post-Newtonian (PN) expansions.

• Limitations: Existing PN-based EOB models often lack sufficient accuracy for next-generation detectors, particularly in the strong-field regime near merger.
• Gap: Recent advances in Post-Minkowskian (PM) gravity have provided high-order scattering results (up to 3PM/NNLO and 4PM) for tidal interactions. However, these results had not yet been systematically transcribed into the EOB framework for both unbound (scattering) and bound (inspiral) orbits.
• Goal: To integrate state-of-the-art PM tidal scattering results into four different flavors of the EOB formalism, validate them against numerical relativity (NR) scattering data, and assess their impact on bound orbit dynamics.

2. Methodology

The authors employed a "matching" procedure to transcribe high-order PM scattering results into EOB potentials. The core methodology involves:

• Gauge-Invariant Matching: The authors equate the analytical EOB scattering angle ($\chi_{\text{EOB}}$) with the real scattering angle ($\chi_{\text{real}}$) derived from PM theory and NR simulations. Since the scattering angle is gauge-invariant, this allows for the unique determination of EOB potential coefficients.
• Four EOB Variants: The study constructs tidal sectors for four distinct EOB formulations:
  1. PM EOB (PS Gauge): A Hamiltonian formulation in the post-Schwarzschild gauge, where higher-order PM effects are encoded in a non-geodesic $Q$-term.
  2. Lagrange-EOB (LJBL Gauge): A Lagrangian formulation in the Lagrange-Just-Boyer-Lindquist gauge, where tidal effects modify the energy-dependent $A$-potential.
  3. Lagrange-EOB (w-EOB Gauge): A formulation using the radial potential $w(\bar{u}, \gamma)$ in isotropic coordinates.
  4. PN EOB (DJS Gauge): The standard PN-based EOB in the Damour-Jaranowski-Schäfer gauge, extended to include high-order tidal terms for unbound orbits.
• Inclusion of Post-Adiabatic Effects: The models incorporate both adiabatic (instantaneous) and post-adiabatic (time-delayed/dynamical) tidal effects. The latter are parameterized by coefficients $\kappa_{\dot{E}2}$ and $\kappa_{\dot{B}2}$, which account for the finite response time of the neutron stars.
• Regularization: To handle divergent integrals arising from the PM expansion of the scattering angle, the authors utilize Hadamard's partie finie (Pf) regularization.
• Validation: The analytical models are compared against recent NR scattering data for equal-mass non-spinning neutron stars (using SLy and MS1b EOS) from Fontbuté et al. [30].

3. Key Contributions

The paper provides several major theoretical and practical advancements:

• High-Order Tidal Coefficients:
  • PM EOB: Determined three new coefficients for the tidal $Q$-term, completing the adiabatic tidal sector up to $O(G^8)$ (3PM) with octupolar contributions at $O(G^9)$.
  • LEOB (LJBL & w-EOB): Derived the tidal part of the $A$-potential (LJBL) and $w$-potential (w-EOB) up to 3PM, introducing four new coefficients for each gauge.
  • PN EOB: Determined twelve new coefficients for the $A$, $D$, and $Q$ potentials in the DJS gauge. This includes the $O(u^3)$ contribution to the gravitoelectric amplification factor, complete knowledge of the tidal $D$-potential up to $O(u^8)$, and the leading-order post-adiabatic contributions.
• Post-Adiabatic Formalism: The paper explicitly integrates post-adiabatic tidal terms (involving proper-time derivatives of tidal tensors) into the EOB potentials. It discusses the logarithmic running of these coefficients and provides estimates based on neutron star f-mode frequencies.
• Unbound to Bound Mapping: The authors demonstrate that the coefficients derived from unbound scattering data (valid for $\gamma > 1$) can be analytically continued to bound orbits ($\gamma < 1$) within the PN framework, as the scattering angle's PM expansion contains only local information up to the orders considered.

4. Results

• Comparison with Numerical Relativity (Scattering):

  • The new LEOB models (both LJBL and w-EOB gauges) show significantly improved agreement with NR scattering data compared to existing PN-EOB models and non-resummed PM expansions.
  • For angular momenta $J_{in}/M^2 \gtrsim 1.5$, the LEOB models match NR data within error bars.
  • Discrepancy at Low Angular Momentum: At lower angular momenta (closer to the plunge), all analytical models overestimate the tidal scattering angle. The authors attribute this partly to the need for improved resummation schemes and partly to hydrodynamic effects (mass exchange) present in NR but not in the point-particle EOB model.
  • Post-Adiabatic Impact: Including physically estimated positive post-adiabatic coefficients ($\kappa_{\dot{E}2} \sim 100-200$) increases the discrepancy with NR. However, fitting $\kappa_{\dot{E}2}$ as a free parameter yields negative values, suggesting that current resummation schemes or the definition of post-adiabatic Love numbers in the EFT context may need refinement.

• Impact on Bound Orbits (PN EOB):

  • The new $O(u^3)$ term in the gravitoelectric amplification factor is negative for equal-mass binaries, contrasting with lower-order positive terms. This suggests a sign flip in the strong-field regime, consistent with gravitational self-force (GSF) results.
  • Orbital Frequency: The inclusion of new tidal terms in the $A$ and $D$ potentials leads to a positive shift in the orbital frequency $\Omega$.
    • The $D$-potential contribution increases $\Omega$ by $\sim 1\%$.
    • The $Q$-term contribution is smaller ($\sim 0.01\%$).
    • The post-adiabatic contribution to the $A$-potential is expected to be the dominant effect, potentially improving the phase agreement with NR by making the potential more attractive.

5. Significance

• Foundation for PM-EOB: This work lays the groundwork for a fully PM-based EOB model that incorporates tidal effects to high orders, bridging the gap between scattering amplitude methods and waveform modeling.
• EOS Constraints: By improving the accuracy of the tidal sector, these models are essential for extracting precise tidal deformability parameters from future GW detections, which directly constrain the neutron star EOS.
• Resummation Challenges: The results highlight the limitations of current Taylor-expanded or simple Padé-resummed PM series in the strong-field regime. The paper calls for new resummation strategies that account for the singularity structure of the scattering angle near the plunge.
• Post-Adiabatic Physics: The study underscores the complexity of defining and matching post-adiabatic Love numbers in the EFT framework, noting that current estimates may be insufficient for the high-precision requirements of next-generation detectors.

In summary, this paper successfully transcribes cutting-edge PM scattering results into the EOB framework, significantly improving the description of tidal interactions in neutron star scattering and providing a robust foundation for future high-precision gravitational waveform models.