Peeling and Polyhomogeneous Asymptotics

Peeling is a hierarchy of curvature decay along outgoing null directions, not merely the assertion that the metric approaches Minkowski space. In four vacuum dimensions, the classical hierarchy follows when the conformally rescaled curvature is regular enough at null infinity. Generic initial data or an unfortunate asymptotic gauge may instead produce logarithms. The resulting polyhomogeneous expansion can still describe radiation, but it need not have the differentiability required for full peeling.

Required background. Analysis for Mathematical Relativity supplies weighted estimates, and exterior algebra and differential forms supplies null-frame curvature decompositions.

Helpful background. The Einstein equation supplies the vacuum Bianchi identities; Causality and Global Structure supplies null infinity; PDEs and Green Functions for Gravitation supplies asymptotic transport; and Gravitational Radiation and Waveforms supplies the waveform interpretation.

The classical hierarchy

Choose an outgoing Newman–Penrose tetrad $(l,n,m,\bar m)$ with $l$ tangent to outgoing rays, affine radius $r$ in a Bondi frame, and normalization $g(l,n)=-1$, $g(m,\bar m)=1$, all other inner products zero. On this page the Weyl scalars are

$$\begin{aligned} \Psi_0&=-C(l,m,l,m),& \Psi_1&=-C(l,n,l,m),\\ \Psi_2&=-C(l,m,\bar m,n),& \Psi_3&=-C(l,n,\bar m,n),\\ \Psi_4&=-C(n,\bar m,n,\bar m). \end{aligned}$$

With these sign and orientation choices, the classical orders are

$$\Psi_0=O(r^{-5}),\quad \Psi_1=O(r^{-4}),\quad \Psi_2=O(r^{-3}),\quad \Psi_3=O(r^{-2}),\quad \Psi_4=O(r^{-1}).$$

The leading $r^{-1}$ term in $\Psi_4$ is outgoing radiation; the middle terms encode longitudinal and Coulombic information; the fastest-decaying $\Psi_0$ is the incoming-radiation component in this tetrad. A boost of the tetrad or a change of radial normalization changes component representatives, so a peeling statement must name its asymptotic frame.

One precise sufficient class is a smooth asymptotically simple vacuum spacetime near $\mathscr I^+$. Thus the physical metric $\widetilde g$ admits a $C^\infty$ unphysical metric $g=\Omega^2\widetilde g$ on a manifold with boundary, with $\Omega=0$ and $d\Omega\ne0$ there, and the rescaled Weyl tensor $K^a{}_{bcd}=\Omega^{-1}C^a{}_{bcd}$ extends smoothly to the boundary. The unphysical Bianchi equation $\nabla_aK^a{}_{bcd}=0$ is regular at $\mathscr I^+$. Transport along its generators and conversion from a regular unphysical tetrad to the physical affine tetrad then supply successive powers of $\Omega\sim r^{-1}$. This is the smooth conformal proof of the displayed orders Penrose 1963, pp. 66–68; the tensor rescaling and boundary argument are written explicitly in Frauendiener 2004, §2.3, eqs. (10)–(11). It is a sufficient route, not a necessary-and-sufficient characterization of all peeling data.

The powers can be checked without memorizing the five falloffs. Use a regular unphysical tetrad $(L,N,M,\bar M)$ near the chosen boundary cut. A physical tetrad adapted to outgoing affine rays scales as $\widetilde l=\Omega^2L$, $\widetilde n=N$, and $\widetilde m=\Omega M$, up to smooth finite normalizations. Conformal invariance of the mixed Weyl tensor gives the all-lowered physical tensor $\widetilde C_{abcd}=\Omega^{-1}K_{abcd}$, where indices on $K$ are lowered with the unphysical metric. Its contraction with $(l,m,l,m)$ therefore carries $\Omega^{-1}\Omega^2\Omega\Omega^2\Omega=\Omega^5$; $(n,\bar m,n,\bar m)$ carries $\Omega^{-1}\Omega\Omega=\Omega$. The three intermediate contractions give the remaining powers. Regularity of $K$, and a tetrad with these controlled asymptotic weights, are the essential inputs.

Polyhomogeneous expansions

A polyhomogeneous field admits an expansion of the form

$$F(u,r,\vartheta) \sim \sum_{j\ge j_0}\sum_{k=0}^{K_j} r^{-j}[\log(r/r_0)]^k F_{jk}(u,\vartheta).$$

Here $r_0\gt0$ is a fixed reference length, so the logarithm is dimensionless. Changing $r_0$ mixes coefficients at a given power but does not remove a nonzero leading logarithmic obstruction.

Such an expansion is organized and differentiable to controlled orders, but it is not smooth in $x=1/r$ at $x=0$ when a nonzero logarithmic coefficient is present. Radial transport equations explain the mechanism: if a source at one order resonates with the homogeneous power, integrating $r^{-1}$ produces $\log(r/r_0)$. For example, the exact radial equation $r\partial_r F+jF=A r^{-j}$ becomes $\partial_r(r^jF)=A/r$, so $F=r^{-j}[B+A\log(r/r_0)]$. Constraint equations can then pass such a term to curvature.

In generalized harmonic formulations this issue has two distinct sources. Some logarithms reflect the asymptotic geometry or initial data; others are introduced by gauge source functions and constraint additions. Duarte, Feng, Gasperin, and Hilditch construct formal polyhomogeneous solutions of a model class and of the generalized-harmonic Einstein system, find peeling- obstructing logarithms in harmonic gauge, and exhibit choices suppressing that particular mechanism to arbitrarily high order Duarte et al. 2022, §§2–5 and eqs. (68)–(76). Their analysis is formal asymptotics; it is not a global existence theorem for the expanded solutions.

Data at spatial infinity

For a Cauchy problem, outgoing null hypersurfaces meet spatial infinity at their past end. Smoothness there imposes compatibility conditions not visible in a finite weighted Sobolev norm. In the special class of asymptotically Euclidean, time-symmetric data that are conformally flat near infinity, Valiente Kroon’s cylinder-at-spatial-infinity calculation produces obstructions through order six. Nonvanishing Newman–Penrose constants and higher quantities generate nonsmooth terms Valiente Kroon 2004, §§4–6, pp. 143–153.

The limitation matters: that calculation establishes the displayed finite-order obstructions in a restricted data class and motivates a rigidity conjecture; it does not prove that every smooth null infinity forces arbitrary data to be Schwarzschild near spatial infinity. Conversely, a global stability theorem with weak curvature decay does not automatically supply the extra regularity needed to upgrade its conclusion to full classical peeling.

What survives without full peeling

Failure of the $r^{-5}$ order for $\Psi_0$ does not by itself eliminate a well-defined leading radiative field or finite Bondi flux. One can often retain the leading shear and news while accepting logarithms in subleading metric or curvature coefficients. The exact outcome depends on three axes:

• the differentiability or polyhomogeneous index set at $\mathscr I$;
• the tetrad and asymptotic gauge used to define the components;
• the initial-data compatibility imposed at spatial infinity.

Thus “asymptotically flat,” “has radiation,” “admits a conformal boundary,” and “peels” are different conclusions. A theorem should never replace one with another without an explicit regularity argument.

Diagnose a logarithmic obstruction

Suppose a curvature component has $\Psi_0=r^{-5}[A(u,\vartheta)\log(r/r_0)+B(u,\vartheta)] +O(r^{-6}[\log(r/r_0)]^q)$. Does it satisfy classical peeling, and what conformal regularity fails?

Solution

It has the expected power in the weak estimate $\lvert\Psi_0\rvert=O(r^{-5}\log(r/r_0))$, but not the classical expansion with a finite smooth coefficient $\lim_{r\to\infty}r^5\Psi_0$. With $x=1/r$, the coefficient is $-A\log(xr_0)+B$, which is unbounded when $A\ne0$. The rescaled curvature therefore fails the continuity, and hence the differentiability, needed by the smooth-conformal-boundary proof of peeling. If $A=0$, this particular obstruction disappears, though other components or higher orders must still be checked.

References

• Duarte, Miguel, Justin C. Feng, Edgar Gasperin, and David Hilditch. “Peeling in Generalized Harmonic Gauge.” Classical and Quantum Gravity 39 (2022): 215002. DOI.
• Frauendiener, Jörg. “Conformal Infinity.” Living Reviews in Relativity 7 (2004): 1. DOI.
• Penrose, Roger. “Asymptotic Properties of Fields and Space-Times.” Physical Review Letters 10 (1963): 66–68. DOI.
• Valiente Kroon, Juan Antonio. “A New Class of Obstructions to the Smoothness of Null Infinity.” Communications in Mathematical Physics 244 (2004): 133–156. DOI.