Asymptotic expansions for semilinear waves on asymptotically flat spacetimes

Technical Summary: Asymptotic Expansions for Semilinear Waves on Asymptotically Flat Spacetimes

Problem Statement
This paper investigates the precise late-time asymptotic behavior of solutions to semilinear wave equations with power-type nonlinearities on stationary, asymptotically flat spacetimes. The primary goal is to move beyond standard decay rate estimates to establish full asymptotic expansions, explicitly identifying the leading-order tail terms and their coefficients. The study focuses on two distinct regimes of nonlinearity: cubic nonlinearities ($p=3$) and higher-order power nonlinearities ($p \geq 4$). The analysis is conducted on the exterior region of subextremal black holes (modeled by metrics satisfying specific spectral admissibility conditions), utilizing a resolved spacetime framework to capture behavior at spatial infinity, null infinity, and timelike infinity simultaneously.

Methodology
The authors employ a synthesis of geometric microlocal analysis (specifically the $b$-calculus) and classical physical-space techniques. The core methodology involves:

1. Radiation Field Expansions: The authors derive high-order radiation field expansions at future null infinity ($\mathscr{I}^+$). By utilizing available pointwise decay estimates for semilinear waves, they construct expansions for the solution $\phi$ in terms of powers of the radial coordinate $\rho = |x|^{-1}$ and time $t_*$. Crucially, they establish the existence of second-order (for $p=3$) and third-order (for $p \geq 4$) radiation fields, which are necessary to characterize the nonlinear forcing terms with sufficient precision.
2. Low-Energy Resolvent Analysis: The initial value problem is converted into a forcing problem via a time cutoff. The solution is then analyzed in the frequency domain using the Fourier transform. The authors extend the low-energy resolvent expansion techniques developed by Hintz for linear waves to inputs with conormal asymptotics. This involves analyzing the resolvent operator $\widehat{\square}_g(\sigma)^{-1}$ near zero energy ($\sigma \to 0$).
3. Conormal Spaces and Resolved Geometry: The analysis is performed on a resolved spacetime manifold $M^+$, obtained by compactifying the causal future and blowing up the corner where time and radius tend to infinity. This geometric framework allows for a unified description of asymptotics across different boundary hypersurfaces (spatial infinity $K^+$, null infinity $\mathscr{I}^+$, and timelike infinity $I^+$). The authors utilize conormal spaces $A(E, \alpha)$ to handle the singular behavior of the resolvent and the forcing terms.
4. Separation of Frequency Regimes: The solution is decomposed into low-frequency ($\phi_0$), intermediate-frequency ($\phi_1$), and high-frequency ($\phi_2$) components. The high-frequency part is shown to decay rapidly. The low-frequency part is analyzed using the resolvent expansion, where the singular terms (logarithmic in frequency) dictate the leading-order polynomial decay in time.

Key Contributions and Results

• Cubic Nonlinearities ($p=3$):

  • The authors prove that for the equation $\square_g \phi = a(t_*, x)\phi^3$, the solution admits a full asymptotic expansion globally on the resolved spacetime.
  • On any fixed compact spatial set $K$, the solution behaves as:
    $$\phi(t_*, x) = 2c_0 t_*^{-2} + O(t_*^{-3+})$$
  • The coefficient $c_0$ is explicitly computable in terms of the initial data, the background geometry (mass $m$), and the first radiation field. Specifically, $c_0$ is given by an integral involving the radiation field $\phi_{rad,1}$ and the angular coefficients of the initial data and metric.
  • This result extends the classical linear Price's law (which predicts $t^{-3}$ decay for linear waves with non-zero angular momentum or specific initial data conditions, though here the nonlinear interaction generates a $t^{-2}$ tail) to the semilinear cubic case.

• Higher-Order Nonlinearities ($p \geq 4$):

  • For equations of the form $\square_g \phi = b(t_*, x)\phi^p$ with $p \geq 4$, the leading-order behavior on compact sets is shown to be:
    $$\phi(t_*, x) = d t_*^{-3} + O(t_*^{-4+})$$
  • The constant $d$ is explicitly determined by the background metric, the nonlinearity, and integrals of the radiation fields.
  • This confirms that for $p \geq 4$, the nonlinear contribution does not alter the leading decay rate of the linear Price's law ($t^{-3}$), but it does modify the coefficient.
  • The paper notes that in the Minkowski limit ($m=0$) for $p \geq 5$, the leading coefficient vanishes trivially, and the leading term becomes $O(t_*^{-4})$, requiring higher-order expansions of both radiation fields and the resolvent.

• Global Asymptotics:

  • Theorem 1.3 provides a global description of the solution on the resolved spacetime $M^+$. It characterizes the leading terms at the three distinct boundaries:
    • At spatial infinity ($K^+$): $2c_0 t_*^{-2}$.
    • At timelike infinity ($I^+$): $2c_0 t_*^{-2} \frac{\rho t_*}{\rho t_* + 2}$.
    • At null infinity ($\mathscr{I}^+$): $\rho \phi_{rad,1}(t_*, \omega)$, where $\phi_{rad,1}$ decays as $c_0 t_*^{-1}$.

Significance and Claims
The authors claim that these results "sharpen previous decay estimates" found in the literature (such as those by Angelopoulos, Aretakis, Gajic, and others) by providing precise asymptotic expansions with explicit leading-order terms and computable coefficients, rather than just decay rates.

The significance of the work lies in:

1. Extension of Price's Law: It rigorously extends the classical linear Price's law to semilinear settings, demonstrating how nonlinearities affect the late-time tails.
2. Methodological Advancement: The paper successfully adapts Hintz's linear low-energy resolvent expansion approach to the semilinear setting. This is achieved by developing a systematic method to handle nonlinear forcing terms with conormal remainders and deriving high-order radiation field expansions.
3. Relevance to Strong Cosmic Censorship (SCC): The authors note that quantitative information about late-time tails, particularly along the event horizon, is a critical input for analyzing the Strong Cosmic Censorship conjecture. While this paper restricts itself to the exterior region, the precise tail coefficients ($c_0$ and $d$) provide a refined "exterior ingredient" for models relevant to SCC, offering a more detailed profile of the exterior field than decay rates alone.

The paper explicitly states that it assumes the global existence of solutions (a standard assumption in this context, often guaranteed for small data or defocusing nonlinearities) and does not address the existence theory itself, focusing instead on the asymptotic structure of existing global solutions. It also notes that while the current work focuses on power-type nonlinearities, the core methodology is robust and potentially applicable to more general nonlinear structures, such as those with derivative nonlinearities.