The Hidden Simplicity of Neutron Star and Black Hole Scattering ✛

Original: "Hidden simplicity in the scattering for neutron stars and black holes"
· Rafael Aoude, Andreas Helset
arXiv:2509.04425 · 2025-09-04 · CC BY 4.0 · 3 min
Kerr generating functions reveal a simple structure in multiloop scattering amplitudes of a neutron star off a rotating black hole.

Context

Scattering of compact objects such as neutron stars and black holes underpins the description of gravitational waves from mergers. Traditional post-Minkowskian calculations quickly become complicated as the number of loops and the spin order increase. This work sits alongside studies that use quantum field theory methods for scattering amplitudes, including resummation of spin effects and calculations in Schwarzschild backgrounds. However, here for the first time, the rotation of the black hole is systematically accounted for to all orders in spin for multiloop tidal processes.

Methods

The authors use the heavy particle effective theory (HPET) to describe a rotating black hole as a heavy particle with higher spin. The key tool is the Kerr generating functions, which are multiloop integrals with an exponential factor depending on the spin vector. These functions are computed recursively and expressed in terms of hypergeometric series. Tensor reduction is achieved by differentiating with respect to the spin vector, replacing cumbersome integration-by-parts identities. The neutron star is modeled by tidal operators built from the electric and magnetic components of the Weyl tensor.

Results

They obtain compact expressions for the scattering amplitudes of a neutron star off a Kerr background up to four loops (order $G^5$) for the leading nonlinear tidal operators. The spin dependence of the black hole is summed to all orders and expressed in terms of hyperbolic cosine and generalized hypergeometric functions. The amplitudes are organized by the helicities of the exchanged gravitons (MHV, NMHV, $N^{2MHV}$), and in each sector the same structure emerges, generated by the same generating functions. In the spinless limit, the results agree with known calculations in a Schwarzschild background. They also show that a special choice of Wilson coefficients leads to boost-invariant operators (generalized Kretschmann scalars), for which the spin dependence reduces to a Newman–Janis shift.

Implications

This work stands out from neighboring studies by pushing spin resummation to four loops for tidal interactions. Previously, such results were limited to one or two loops, or only considered spinless backgrounds. The Kerr generating functions provide a universal method that can be applied to any process in a rotating black hole background, including the scattering of two black holes. This significantly simplifies calculations and opens the way to higher orders in perturbation theory.

Future development

The generating function method can be extended to the full scattering of two rotating black holes, where both particles carry spin. It is also natural to include gravitational wave emission and self-force corrections. Further development may connect these results to exact solutions of the equations of motion in a Kerr background and to calculations within self-force theory.

Impact

The results will impact the accuracy of gravitational wave models for neutron star–black hole systems, especially when accounting for tidal effects and rotation. The method is also useful for high-energy theoretical physics, where similar techniques are applied in QCD.

Next steps

The next step is to apply the Kerr generating functions to the scattering of two rotating black holes and to include radiative corrections. It is also important to explore the connection with classical solutions and to verify the results numerically.

Key open problems

This work is directly related to the unsolved problem of accurately describing the two-body dynamics in general relativity, especially with spins and tidal deformations. The question of complete resummation of all loop corrections and the behavior of the perturbative series at high energies remains open.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

This work introduces Kerr generating functions for the tensor reduction of multiloop integrals in the background of a rotating black hole. Using them, compact scattering amplitudes for a neutron star are obtained up to four loops, with the spin resummed to all orders. A universal structure is uncovered that connects different loop orders and helicity sectors.

Where its strength is

The Kerr generating function method allows the black hole spin to be systematically incorporated at all orders, replacing cumbersome tensor reduction with differentiation. This opens a direct route to computations at higher orders in perturbation theory.

Where it could go next

Worth trying

  • Test whether the hidden simplicity persists when the finite size of the neutron star is taken into account via higher-order tidal operators.
  • Use Kerr generating functions to compute radiative corrections to the impulse in the scattering of two spinning bodies.
  • Compare predictions for spin-dependent tidal effects with numerical simulations of neutron star–black hole mergers.

Why this may matter

Such computations could improve the accuracy of gravitational-wave templates for next-generation detectors, which is important for extracting information about the nature of compact objects and testing general relativity.

The Kerr generating functions vanish when acted upon by the Laplacian with respect to the spin vector, reflecting the special structure of the interaction of a rotating black hole with the gravitational field.

The idea of hidden simplicity in complex calculations echoes a motif from Carl Sagan's novel "Contact," where a message from an extraterrestrial civilization is encoded in seemingly random data, but with the right approach reveals a deep structure.

G^{(L,k)}(a) = \int \hat{d}^D\ell_1 \cdots \hat{d}^D\ell_L \frac{1}{(2v_1\cdot\ell_1)\cdots(2v_1\cdot\ell_{1\cdots L})} \frac{1}{\ell_1^2\cdots\ell_L^2(q-\ell_{1\cdots L})^2} \exp[\ell_{1\cdots k}\cdot a]
Integral over L loop momenta with an exponential factor depending on the spin vector a; differentiation with respect to a yields tensor integrals.
c^{(L,k)}_{i,j} = \frac{2^{2i+j} (d(L+1-k))_i (dk)_{i+j}}{(dL)_i (d(L+1))_{2i+j}} \, {}_3F_2\left(\genfrac{}{}{0pt}{}{d(L+1-k)+i, \frac{1-j}{2}, -\frac{j}{2}}{dL+i, 1-dk-i-j}; 1\right)
The coefficients are expressed in terms of the hypergeometric function 3F2 and Pochhammer symbols; they determine the contribution of each order in spin.
\chi_{\mathrm{Kr}} = \frac{1}{2} \sum_{\pm} \chi^{(k=0,s=0)}_{\phi^2 R^{L+1}}(b \mp i a)
The phase is expressed in terms of the spinless eikonal phase with a complex shift of the impact parameter (Newman–Janis shift).

Key numbers

  • maximum loop order: 4 (order G^5)
  • order in spin: all orders (resummed)
  • number of independent tidal operators at 3-4 loops: 5
  • spacetime dimension in calculations: D (general, with limit D=4 for spin vector)
  • coefficient of c_L at 4 loops: 1/42
Scientists
Stephen HawkingJacob BekensteinFritz ZwickyBernhard RiemannJoseph WeberKarl Schwarzschild
Concepts
kerr black hole Penrose process effective spin black hole neutrino source kerr effect accretion disk wind torus spin rate neutron scattering Hawking radiation hawking radiation mechanism ringdown active galaxy gravastar kerr metric rotating black hole wald entropy neutron star dark matter dark matter lensing