Radial fall: the gravitational waveform up to the second-and-half Post-Newtonian order

This paper applies the Multipolar Post-Minkowskian formalism to compute the gravitational waveform, energy, and linear momentum emissions for a two-body system in radial fall up to the 2.5 Post-Newtonian order, while also evaluating inertial force contributions at the 4.5PN level to facilitate future high-precision calculations.

The Gist

Imagine two heavy objects, like two bowling balls, floating in space. Usually, when we talk about black holes and gravity, we imagine these balls orbiting each other like dancers, slowly spiraling inward until they crash. But this paper is about a much more dramatic, "head-on" scenario: What happens if one object simply drops straight down into a black hole, with zero sideways motion?

The authors, Donato Bini and Giorgio Di Russo, are essentially trying to write the "script" for the gravitational waves (ripples in spacetime) that would be created during this straight-line plunge.

Here is a breakdown of their work using simple analogies:

1. The Problem: The "Too Fast" Fall

In physics, we have a powerful tool called the Post-Newtonian (PN) approximation. Think of this as a set of rules that work great when things are moving slowly and gravity is weak (like planets orbiting the Sun).

However, when an object falls straight into a black hole, it starts far away (where the rules work) but speeds up incredibly fast as it gets close to the "event horizon" (the point of no return).

• The Analogy: Imagine trying to predict the path of a car using only the laws of driving on a flat highway. That works fine on the highway. But as soon as the car hits a steep, twisting mountain road and starts accelerating dangerously, your highway rules break down.
• The Paper's Goal: The authors wanted to use their "highway rules" (PN approximation) to describe the fall as far as they possibly could before the car hits the mountain. They stopped the calculation just before the object gets too close to the black hole, where the math gets too messy for their current tools.

2. The "Brake" That Pushes Back (Radiation Reaction)

One of the key findings is about radiation reaction.

• The Analogy: Imagine a swimmer diving into a pool. As they move, they create waves. Those waves carry away energy. Because energy is leaving the swimmer, the water pushes back against them, slowing them down slightly.
• In the Paper: As the object falls, it creates gravitational waves (ripples in space). These ripples carry energy away. This loss of energy acts like a tiny "brake" or a drag force on the falling object. The authors calculated exactly how this "brake" changes the object's fall at a very specific level of precision (called 2.5PN). They found that this force modifies the fall, creating a specific type of burst of radiation (called bremsstrahlung, or "braking radiation").

3. The "Ghost" Forces (Inertial Forces)

The paper also looks at what happens to the "center of mass" (the average center point between the two objects).

• The Analogy: Imagine you are standing on a skateboard holding a heavy ball. If you throw the ball hard to the right, you will roll to the left to conserve momentum.
• In the Paper: As the falling object shoots out gravitational waves in one direction, it pushes the black hole (and the system's center) slightly in the opposite direction. This creates "inertial forces" or "dragging effects." The authors calculated these forces, which are like the "recoil" of the system. They noted that while these forces are tiny now, they become important for future, ultra-precise calculations.

4. The "Soundtrack" of the Fall (The Waveform)

The main output of the paper is the gravitational waveform.

• The Analogy: If the falling object were a singer, the gravitational waves would be the sound recording of their voice. The authors didn't just guess the song; they wrote out the sheet music note-by-note, calculating the pitch and volume up to a very high level of detail.
• The Result: They found that the "song" (the energy emitted) has a specific shape. Interestingly, they found that the loudest part of the "song" happens while the object is still far away, in the "weak gravity" zone, which is exactly where their math works best. This gives them confidence that their calculations are accurate for that part of the fall.

5. Why Does This Matter?

You might ask, "Why study a straight drop when real black holes usually have spinning, orbiting partners?"

• The "Control Group" Analogy: In science, to understand a complex system, you often start with the simplest possible version. A straight drop is the "simplest" version of a black hole collision.
• The Future: By perfecting the math for this simple "straight drop," the authors are building a foundation. Once they have this "simple script" perfect, they can add complexity (like spin or orbits) to understand the more chaotic, real-world collisions that detectors like LIGO actually see.

Summary

This paper is a high-precision mathematical simulation of a heavy object falling straight into a black hole. The authors used advanced physics tools to predict:

1. How the object falls (its path).
2. How the "braking" from gravitational waves changes that path.
3. How the system recoils (pushes back) as it emits energy.
4. The exact "sound" (waveform) of the event.

They stopped the simulation just before the object hit the black hole's "point of no return," acknowledging that their tools work best in the "slow and far" zone, but they successfully mapped out the entire approach up to that critical moment. This work serves as a crucial building block for future, even more accurate models of cosmic collisions.

Technical Summary

1. Problem Statement

The paper addresses the gravitational radiation emitted by a two-body system undergoing radial infall (a "head-on" collision) into a black hole. While this problem has a long history involving numerical relativity and perturbation theory (e.g., the Zerilli and Teukolsky formalisms), there is a gap in analytical Post-Newtonian (PN) treatments for this specific trajectory.

• The Challenge: Unlike inspiraling binaries (which spend many cycles in the weak-field regime) or scattering events, a radial fall involves a particle starting from rest at infinity and rapidly entering the strong-field regime near the event horizon.
• The Limitation: Standard PN expansions break down as the particle approaches the horizon (strong-field region).
• The Goal: The authors aim to perform a partial analytical study valid in the weak-field regime (from infinity down to a cutoff time $t_{max}$ before the horizon) to compute the gravitational waveform, energy loss, and momentum losses up to 2.5PN order. This serves as a benchmark for future numerical studies and a "computational block" for higher-order analytical work.

2. Methodology

The authors employ the Multipolar Post-Minkowskian (MPM) formalism, a rigorous framework for computing gravitational waveforms in the post-Newtonian approximation.

• Formalism: They utilize the MPM approach to decompose the asymptotic waveform ($h_{ij}^{TT}$) into radiative multipole moments ($U_\ell, V_\ell$). These are expressed as nonlinear functionals of "source moments" (mass and current multipoles derived from the matter distribution) and "gauge moments."
• Coordinate System: The calculations are performed in modified harmonic coordinates within the Center-of-Mass (CM) frame.
• Orbit Modeling:
  • The relative motion is modeled as a particle released from rest at infinity.
  • The trajectory is computed including radiation-reaction forces up to 2.5PN order.
  • The authors derive the position $x_p(t)$ and velocity as a series expansion in the PN parameter $\eta = 1/c$.
• Multipole Calculation:
  • They explicitly calculate the source mass quadrupole $I_{ij}$ and higher-order moments ($I_{ijk}, I_{ijkl}$, etc.) by substituting the PN-expanded orbit into the definitions.
  • They account for the specific symmetries of radial motion (e.g., vanishing angular momentum $J_a=0$, diagonal quadrupole tensor structure).
  • They compute the radiative moments $U_\ell$ by applying the MPM formulas, which include hereditary terms (tail effects) and memory terms.

3. Key Contributions

The paper provides the first complete analytical derivation of the gravitational waveform for radial infall up to 2.5PN accuracy.

• Waveform Derivation: Explicit expressions for the complex waveform $h_c$ (and its frequency domain counterpart) are derived, decomposed into spin-weighted spherical harmonics ($h_{\ell m}$).
• Radiation Reaction: The authors incorporate the 2.5PN radiation-reaction force into the equations of motion, showing how it modifies the trajectory compared to the conservative (geodesic) case.
• Loss Quantification: They calculate the fluxes of:
  • Energy ($dE/dt$): Up to 2.5PN.
  • Angular Momentum: Proven to vanish identically due to the radial nature of the motion.
  • Linear Momentum: Computed up to 2.5PN, revealing a non-zero recoil (kick) effect.
• Inertial Forces (4.5PN): Although the main waveform is 2.5PN, the authors compute the inertial forces (dragging effects) appearing in the CM frame at 4.5PN order. These arise because the system loses linear momentum, causing the CM frame to become non-inertial. This is a preparatory step for future 4.5PN waveform calculations.

4. Key Results

A. The Waveform and Multipole Expansion

The waveform is expanded in powers of $\eta$ (PN order). The dominant contribution comes from the mass quadrupole ($U_2$), but higher-order multipoles ($U_3, U_4, \dots$) are non-zero due to the non-linearities of the MPM formalism, even in radial motion.

• Time Domain: The waveform components are expressed as functions of the dimensionless time $T = t/M$.
• Frequency Domain: The Fourier transform yields the energy spectrum $dE/d\omega$. The spectrum exhibits a peak at $M\omega \approx 0.36$, which aligns with previous numerical results (e.g., Detweiler's work using the Teukolsky formalism).

B. Energy and Momentum Losses

• Energy Loss: The power radiated ($dE/dt$) is derived explicitly. In the limit of equal masses ($\nu \to 0$), the results match the test-particle limit found in previous literature (Ref. [35]).
• Linear Momentum Kick: The system emits a net linear momentum, leading to a recoil velocity of the center of mass. The rate of change of linear momentum is given by:
  $$\frac{dP_1}{dt} \propto \eta^4 \nu^2 \Delta \frac{1}{T^{13/3}}$$
  where $\Delta = \sqrt{1-4\nu}$ is the mass difference parameter. This implies that asymmetric mass ratios ($\nu \neq 0$) are required for a kick.
• Inertial Forces: The authors calculate the acceleration of the CM frame due to the momentum flux. This "dragging" effect, while appearing at 4.5PN, is crucial for high-precision modeling of the center-of-mass dynamics.

C. Comparison with Numerical Results

The paper validates its analytical approach by comparing the frequency of the energy spectrum peak ($M\omega \sim 0.36$) with numerical relativity and perturbation theory results. The agreement confirms that the PN expansion captures the essential physics of the "burst" of radiation emitted as the particle approaches the strong-field region, even before the horizon is reached.

5. Significance and Future Outlook

• Analytical Benchmark: This work provides a crucial analytical benchmark for numerical relativity simulations of head-on collisions, particularly in the early inspiral/weak-field phase.
• Bridge to Strong Field: While PN theory cannot describe the plunge through the horizon, this work defines the "initial conditions" and the radiation emitted up to the point where PN breaks down.
• Foundation for Higher Orders: By computing the 4.5PN inertial forces and organizing the multipole moments, the authors lay the groundwork for future computations at 3PN, 3.5PN, and 4.5PN orders.
• Contextual Relevance: The study is motivated by recent gravitational wave events (like GW250114 mentioned in the text) and the need to understand quasi-normal modes (QNMs) and tidal deformations. It also connects to broader theoretical frameworks like Conformal Field Theory (CFT) and scattering amplitudes.

In summary, the paper successfully extends the Multipolar Post-Minkowskian formalism to the challenging case of radial infall, delivering a comprehensive 2.5PN description of the gravitational waveform and associated losses, while preparing the theoretical infrastructure for future high-precision analytical studies.