Asymptotic Expansion of the Gravitational Field

We work in Bondi coordinates $(u, r, x^A)$ where $u$ is retarded time, $r$ is the luminosity distance, and $x^A = (\theta, \phi)$ parametrize the celestial 2-sphere. The Bondi gauge is defined by three conditions: (i) $g_{rr} = 0$, (ii) $g_{rA} = 0$, and (iii) $\det(g_{AB})/r^4 = \det(\gamma_{AB})$ where $\gamma_{AB}$ is the unit round metric on $S^2$. The most general metric satisfying $g_{rr} = 0$ and $g_{rA} = 0$ is:

Here $V(u,r,x^A)$, $\beta(u,r,x^A)$, $U^A(u,r,x^A)$, and $g_{AB}(u,r,x^A)$ are four metric functions. Reading off the individual components:

All remaining components $g_{rr} = g_{rA} = 0$ vanish by the gauge choice. This leaves the angular metric $g_{AB}$ as a $2 \times 2$ symmetric matrix on the sphere with one constraint from the determinant condition.

The third Bondi gauge condition states:

This means $\det(g_{AB}) = r^4 \det(\gamma_{AB})$ for all $r$. Now write $g_{AB}$ as a large-$r$ expansion:

Computing the determinant to leading orders using $\det(M + \epsilon N) = \det(M)(1 + \epsilon\,\text{tr}(M^{-1}N) + \ldots)$:

Demanding this equals $r^4\det(\gamma_{AB})$ order by order forces: at $\mathcal{O}(r^3)$: $\gamma^{AB}C_{AB} = 0$ (traceless), and at $\mathcal{O}(r^2)$: $\gamma^{AB}D_{AB} = \frac{1}{2}C_{AB}C^{AB}$. The shear tensor $C_{AB}$ is therefore symmetric and trace-free (STF), the fundamental radiative degree of freedom.

The vacuum Einstein equation $R_{\mu\nu} = 0$ in Bondi gauge gives a hierarchy of equations. The $rr$-component is a radial ODE for $\beta$. Computing $R_{rr}$ from the Christoffel symbols:

The leading-order part simplifies to:

Substituting $g_{AB} = r^2\gamma_{AB} + rC_{AB} + \ldots$ gives $\partial_r g_{AB} = 2r\gamma_{AB} + C_{AB} + \ldots$ and contracting:

The boundary condition $\beta \to 0$ as $r \to \infty$ is built in. This result shows $\beta$ is fully determined by the shear at leading order.

The mixed $rA$-component of the Einstein equation constrains the angular shift vector $U^A$. From $R_{rA} = 0$ one obtains:

Expanding $U^A = U_0^A + \frac{U_1^A}{r} + \frac{U_2^A}{r^2} + \ldots$ and matching powers of $r$, the boundary condition $U^A \to 0$ as $r \to \infty$ forces $U_0^A = 0$. At next order:

Since $g_{uA} = -g_{AB}U^B = -r^2\gamma_{AB}U^B + \ldots$, we recover the subleading $g_{uA}$ expansion. The angular momentum aspect $N_A$ first appears at order $r^{-1}$ in $g_{uA}$.

The news tensor is defined as the retarded-time derivative of the shear:

To verify gauge invariance, consider a supertranslation $u \to u + f(x^A)$. Under this transformation $C_{AB} \to C_{AB} + 2D_{\langle A}D_{B\rangle}f$, where $D_{\langle A}D_{B\rangle}f = D_A D_B f - \frac{1}{2}\gamma_{AB}D^2 f$ is the STF part. Since $f$ depends only on angles, not on $u$:

Therefore $N_{AB}$ is invariant under supertranslations. It is the gravitational analogue of the electromagnetic field strength $F_{\mu\nu}$: gauge-invariant, directly measurable, and encoding the radiative content. When $N_{AB} = 0$, no gravitational radiation passes through $\mathscr{I}^+$.

The constraint equation $G_{uA} = 0$ expanded to order $r^{-2}$ yields the evolution equation for the angular momentum aspect $N_A$. We start from the full Einstein tensor component:

Setting this to zero and collecting terms, we use the identity $D^B N_{AB} = \frac{1}{2}D_A(N_{BC}C^{BC})$ at the linearized level. Including the covariant derivative terms on $S^2$, the full evolution equation reads:

Each term has a physical interpretation: $-\partial_A m_B$ is the gradient of mass loss driving angular momentum change; the $N_{AC}C^{BC}$ terms represent nonlinear coupling between the news and shear; and the $D^B D^2 C_{AB}$ terms arise from the angular structure of the radiation field on $S^2$. This equation is the foundation of the spin memory effect.

On the 2-sphere, introduce a complex null dyad $(m^A, \bar{m}^A)$ satisfying $\gamma_{AB} = m_A \bar{m}_B + \bar{m}_A m_B$ and $m^A \bar{m}_A = 1$. In standard coordinates $m^A = \frac{1}{\sqrt{2}r}(\partial_\theta + \frac{i}{\sin\theta}\partial_\phi)^A$. The shear tensor, being STF, has two independent components captured by a single complex scalar of spin weight $s = -2$:

The Weyl scalar $\Psi_4^0$ at $\mathscr{I}^+$ is related to the second time derivative of the shear:

Identifying $h = h_+ - i h_\times$ with $\bar{\sigma}^0$, we have $\Psi_4^0 = \ddot{h}_+ - i\ddot{h}_\times$. Expanding in spin-weighted spherical harmonics $ {}_{-2}Y_{\ell m}$:

where $Y_{AB}^{E,\ell m}$ and $Y_{AB}^{B,\ell m}$ are the electric and magnetic parity tensor harmonics. The displacement memory resides in the E-mode sector, while the spin memory is encoded in the B-mode sector.

We model a gravitational-wave chirp burst and compute the shear $C_{AB}(u)$, news tensor $N_{AB} = \partial_u C_{AB}$, accumulated shear, and angular momentum aspect evolution: