Skip to content

Scattering Angle and Radial Action

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

The conservative scattering angle is a compact meeting point for otherwise different descriptions of the relativistic two-body problem. It can be read from asymptotic momenta, derived from a Hamiltonian trajectory, obtained by differentiating a radial action, or extracted from an impact-parameter eikonal. Coordinate positions and canonical momenta differ among those descriptions, but the angle at fixed invariant energy and angular momentum does not. That makes it a stringent comparison quantity—provided all calculations use the same masses, asymptotic frame, boundary conditions, and definition of the conservative sector.

Required background. Events, Intervals, and Minkowski Geometry supplies asymptotic Lorentz invariants, and Weak Metric Perturbations and Gauge Freedom supplies the weak-field domain in which PM coefficients are defined.

Helpful background. The Einstein Equation fixes the underlying interaction; Presymplectic Potentials and Currents provides the action and canonical language; and Gravitational Radiation and Waveforms is needed to distinguish a conservative angle from a radiation-reacted one.

Angle from a central Hamilton–Jacobi problem

Section titled “Angle from a central Hamilton–Jacobi problem”

Assume first a time-independent, rotationally invariant, local conservative Hamiltonian in the center-of-mass frame. At fixed total energy E=sE=\sqrt{s} and angular momentum JJ, solve

H(r,pr,J)=EH\bigl(r,p_r,J\bigr)=E

for the radial momentum pr(r;E,J)p_r(r;E,J). Let rmin⁡r_{\min} be the largest positive turning point, chosen by continuing the incoming branch from r=∞r=\infty. Hamilton’s characteristic function is

S=−Et+Jϕ+∫pr(r;E,J) dr.S=-Et+J\phi+\int p_r(r;E,J)\,dr .

Differentiating SS with respect to JJ at fixed EE gives the azimuth accumulated along the orbit. The incoming and outgoing halves yield

χ(E,J)+π=−2∫rmin⁡∞∂pr(r;E,J)∂J dr.\chi(E,J)+\pi =-2\int_{r_{\min}}^\infty \frac{\partial p_r(r;E,J)}{\partial J}\,dr .

The added π\pi removes the straight-line change in polar angle. Although the integral of prp_r itself diverges at infinity, its JJ derivative is finite after the free contribution is treated consistently. Damour uses precisely this Hamilton–Jacobi relation to connect the real two-body and effective-one-body scattering functions Damour 2016, §§ III and VI.

Free subtraction removes the linear large-distance divergence but leaves a logarithm for a 1/r1/r force. Suppose the outgoing radial momentum has the asymptotic form

pr(r;E,J)−pr(0)(r;E,J)=β(E)r+O(r−2),p_r(r;E,J)-p_r^{(0)}(r;E,J) =\frac{\beta(E)}r+O(r^{-2}),

where the leading long-range coefficient is independent of JJ. A finite hyperbolic radial action therefore needs an infrared subtraction as well:

Irscatt(E,J)=lim⁡R→∞[2∫rmin⁡Rpr dr−2∫rmin⁡(0)Rpr(0) dr−2β(E)log⁡ ⁣(RℓIR)],I_r^{\rm scatt}(E,J) = \lim_{R\to\infty} \left[ 2\int_{r_{\min}}^R p_r\,dr -2\int_{r_{\min}^{(0)}}^R p_r^{(0)}\,dr -2\beta(E)\log\!\left(\frac R{\ell_{\rm IR}}\right) \right],

where pr(0)p_r^{(0)} and rmin⁡(0)=J/p∞r_{\min}^{(0)}=J/p_\infty belong to free motion at the same EE and J>0J\gt0, and ℓIR\ell_{\rm IR} fixes the infrared convention. For example, Newtonian scattering has pr2=p∞2+2GMμ2/r−J2/r2p_r^2=p_\infty^2+2GM\mu^2/r-J^2/r^2, so β=GMμ2/p∞\beta=GM\mu^2/p_\infty: the remaining logarithm can be checked without any relativistic machinery. Its subtraction is JJ independent and therefore cannot affect the angle. Since the derivative of the free radial action is −π-\pi, the free subtraction removes precisely the straight-line angle, giving

  χ(E,J)=−∂Irscatt(E,J)∂J  \boxed{\; \chi(E,J) =-\frac{\partial I_r^{\rm scatt}(E,J)}{\partial J} \;}

at fixed EE and masses. The action is more informative than the angle: its energy derivative also encodes a regulated time delay once the asymptotic clock and infrared subtraction are fixed. For a long-range gravitational force there is no finite absolute time delay relative to an infinitely distant free trajectory without such a prescription.

At first PM order, write

A(E)=2Gm1m2(2γ2−1)γ2−1.A(E) = \frac{2Gm_1m_2(2\gamma^2-1)} {\sqrt{\gamma^2-1}} .

The leading angle is χ1PM=A/J\chi_{\rm 1PM}=A/J. Integration with respect to JJ gives

Ir,1PMscatt=−A(E)log⁡ ⁣(JJ0(E)).I_{r,{\rm 1PM}}^{\rm scatt} =-A(E)\log\!\left(\frac{J}{J_0(E)}\right).

The arbitrary J0J_0 reflects the infrared subtraction; it disappears from χ\chi. This logarithm is the action-space counterpart of the transverse Fourier transform of a massless 1/q21/q^2 exchange. The same relationship between a radial action and the PM angle persists at higher orders, though iterations and tail regions must be included with their correct boundary conditions Bern et al. 2022, main text around eqs. (5)–(8).

A canonical transformation changes the one-form by an exact differential,

pi dqi=Pi dQi+dF.p_i\,dq^i=P_i\,dQ^i+dF .

Consequently, an action integral changes by endpoint terms. If the transformation approaches the identity at the incoming and outgoing boundaries, or if its endpoint terms are independent of JJ, then −∂Irscatt/∂J-\partial I_r^{\rm scatt}/\partial J is unchanged. This is the practical canonical-invariance statement behind comparisons among harmonic-coordinate, ADM, isotropic, amplitude-derived, and EOB Hamiltonians.

The qualifications matter. A transformation that rotates the asymptotic frame, changes the identification of JJ, or carries a JJ-dependent boundary term is not innocuous. Likewise, quoting a coordinate turning radius is not equivalent to quoting χ(E,J)\chi(E,J). A clean comparison specifies the invariant energy, the total canonical angular momentum, the masses and spins, and whether radiative or time-symmetric boundary conditions were used.

The angle also connects directly to the impulse. For elastic planar scattering,

∣Δp1∣=2p∞sin⁡χ2.\lvert\Delta\boldsymbol p_1\rvert =2p_\infty\sin\frac{\chi}{2}.

Thus Impulse and Momentum Transfer supplies a vector check, while the radial action supplies a scalar generating function. Agreement of both exposes sign errors or missing longitudinal iterations that an angle-only comparison can hide.

Take m1=μ≪m2=Mm_1=\mu\ll m_2=M and hold the test body’s specific energy γ\gamma fixed. Since J=μj+O(μ2)J=\mu j+O(\mu^2), the general 1PM result becomes

χ1PMtest=2GMj2γ2−1γ2−1.\chi_{\rm 1PM}^{\rm test} = \frac{2GM}{j} \frac{2\gamma^2-1}{\sqrt{\gamma^2-1}} .

This is the weak-deflection expansion of timelike Schwarzschild geodesic scattering. For γ=1+v∞2/2+O(v∞4)\gamma=1+v_\infty^2/2+O(v_\infty^4) and j=bv∞+O(v∞3)j=bv_\infty+O(v_\infty^3), it gives 2GM/(bv∞2)2GM/(b v_\infty^2). For γ→∞\gamma\to\infty and j=bγ2−1j=b\sqrt{\gamma^2-1}, it gives 4GM/b4GM/b, the leading null deflection. One formula therefore checks both the Newtonian and ultrarelativistic limits without a change of PM order.

The radial-action representation is valid while a real scattering turning point exists and the weak-field series can be continued from infinity to it. It breaks down at capture, when the turning point merges with a separatrix, and it becomes nonuniform near such a merger. For spinning motion that is not planar or integrable, more actions and scattering angles are required. For hereditary conservative tails, the full on-shell functional is nonlocal in time; an “effectively local” pr(r)p_r(r) reconstructed on hyperbolic orbits need not continue unchanged to bound motion. This limitation is central to Bound-Orbit Information from Scattering Data.

Suppose an independently computed scattering angle is χ(E,J)=A(E)/J+O(G2)\chi(E,J)=A(E)/J+O(G^2). Determine the JJ-dependent part of the radial action and explain the integration ambiguity.

Solution

Using χ=−∂JIrscatt\chi=-\partial_J I_r^{\rm scatt},

∂Irscatt∂J=−A(E)J+O(G2).\frac{\partial I_r^{\rm scatt}}{\partial J} =-\frac{A(E)}{J}+O(G^2).

Integration at fixed EE gives

Irscatt=−A(E)log⁡J+C(E)+O(G2)=−A(E)log⁡ ⁣(JJ0(E))+O(G2).I_r^{\rm scatt} =-A(E)\log J+C(E)+O(G^2) =-A(E)\log\!\left(\frac{J}{J_0(E)}\right)+O(G^2).

The function C(E)=A(E)log⁡J0(E)C(E)=A(E)\log J_0(E) depends on the infrared subtraction and affects a convention-dependent time delay through its energy derivative. It cannot affect the scattering angle because it is independent of JJ.

Use the test-body formula above to recover both Newtonian massive deflection and null deflection.

Solution

At small velocity, γ2−1=v∞+O(v∞3)\sqrt{\gamma^2-1}=v_\infty+O(v_\infty^3), 2γ2−1=1+O(v∞2)2\gamma^2-1=1+O(v_\infty^2), and j=bv∞+O(v∞3)j=bv_\infty+O(v_\infty^3). Therefore

χ1PMtest=2GMbv∞2+O ⁣(GMb).\chi_{\rm 1PM}^{\rm test} =\frac{2GM}{b v_\infty^2} +O\!\left(\frac{GM}{b}\right).

At large γ\gamma, (2γ2−1)/γ2−1=2γ+O(γ−1)(2\gamma^2-1)/\sqrt{\gamma^2-1}=2\gamma+O(\gamma^{-1}) and j=bγ+O(γ−1)j=b\gamma+O(\gamma^{-1}), so

χ1PMtest=4GMb+O(γ−2).\chi_{\rm 1PM}^{\rm test} =\frac{4GM}{b}+O(\gamma^{-2}).

The first limit tests the long interaction time of a slow particle; the second tests the energy-independent leading bending of a null ray by a fixed Schwarzschild mass.

  • Bern, Zvi, Julio Parra-Martinez, Radu Roiban, Michael S. Ruf, Chia-Hsien Shen, Mikhail P. Solon, and Mao Zeng. “Scattering Amplitudes, the Tail Effect, and Conservative Binary Dynamics at O(G4)\mathcal O(G^4).” Physical Review Letters 128, 161103 (2022). DOI. Open PDF.
  • Damour, Thibault. “Gravitational Scattering, Post-Minkowskian Approximation and Effective One-Body Theory.” Physical Review D 94, 104015 (2016). DOI. Open PDF.
  • Kälin, Gregor, and Rafael A. Porto. “From Boundary Data to Bound States II: Scattering Angle to Dynamical Invariants (with Twist).” Journal of High Energy Physics 2020, 120 (2020). DOI. Open PDF.