Relatively non-degenerate integrated decay estimates for massless Vlasov fields on Schwarzschild spacetimes

Technical Summary: Relatively Non-Degenerate Integrated Decay Estimates for Massless Vlasov Fields on Schwarzschild Spacetimes

Problem Statement
This paper addresses the asymptotic behavior of massless collisionless systems (described by the massless Vlasov equation) on the exterior of a Schwarzschild black hole. The primary objective is to establish integrated local energy decay (ILED) estimates that are relatively non-degenerate. Previous works on this topic, such as [ABJ18] and [Big23], established ILEDs where the bulk decay term degenerated relative to the boundary energy flux, particularly at the photon sphere ($r=3M$) and the future event horizon ($H^+$). This degeneracy prevented the derivation of sharp decay rates and limited the control of derivatives of the energy-momentum tensor $T[f]$, specifically the radial derivative $\partial_r T[f]$. The authors aim to close this gap by constructing a commutation framework that yields ILEDs without such relative degeneracy, thereby enabling the control of all first-order derivatives of $T[f]$ and establishing sharp time decay.

Methodology
The authors develop a commutation vector field approach tailored to the kinetic nature of the Vlasov equation, distinct from the scalar wave equation methods used in [HK23b, Mav23]. The strategy involves several key technical components:

1. Weight Functions and Trapping: The authors utilize a weight function $\phi_-$ associated with the vector field $G$, which captures the concentration phenomenon of future-trapped null geodesics near the photon sphere. This weight is shown to be a Lyapunov function for the geodesic flow.
2. Symplectic Gradients and Projections: A central innovation is the identification of a vector field $V_+$ on the null-shell $P$. $V_+$ is constructed as the projection (parallel to $\partial_{p_t}$) of the symplectic gradient of $\phi_-$. While $V_+$ commutes well with the geodesic spray $X_g$, its commutator contains a problematic error term proportional to $\phi_- \partial_{p_r} f$.
3. Modified Commutator: To eliminate the problematic error term, the authors introduce a modified vector field $V_+^{\text{mod}} = V_+ + \Omega \Phi r \phi_- \partial_{p_r}$. The correction term $\Phi$ is determined by a transport equation along the geodesic flow. This modification ensures that the commutator $[X_g, V_+^{\text{mod}}]$ allows for the control of the bulk integral without relative degeneracy.
4. Rescaling and Redshift: To handle the behavior near the event horizon, the authors work with rescaled vector fields ($V_+^+ = \Omega^{-1} V_+$ and $V_+^{\text{mod},+} = \Omega^{-1} V_+^{\text{mod}}$) and a rescaled weight $\phi_- = \Omega^{-1} \phi_-$. This captures the redshift effect, ensuring non-degenerate control near $H^+$.
5. Hierarchical Energy Estimates: The proof involves splitting the null-shell into regions (far-away incoming, far-away outgoing, and bounded regions) and constructing a specific $W^{1,1}_{x,p}$ weighted energy norm $E[f]$. This norm includes terms involving $p_N f$, $\phi_- \partial_t f$, $\phi_- X_{|/p|} f$, $p_N V_+^{\text{mod}} f$, and weighted derivatives of $f$.
6. $r^p$-Weighted Method: The non-degenerate ILED is combined with the $r^p$-weighted energy method of Dafermos–Rodnianski to derive time decay rates.

Key Contributions and Results

• Theorem 1.1 (ILED without Relative Degeneration): The paper proves that for sufficiently regular solutions $f$ with finite initial energy, the energy flux $E[f]$ satisfies an ILED where the bulk term is equivalent to the boundary term (up to a factor vanishing at infinity, which is handled by the $r^p$ method). Specifically:
  $$\sup_{\tau \ge 0} E[f](\tau) + \int_0^\infty E[r^{-1} \log^{-2}(2+r) f](\tau) d\tau \le C E[f](0).$$
  This is the first such result for massless Vlasov fields on black hole exteriors.

• Control of Derivatives: The construction of $V_+^{\text{mod}}$ allows for the control of the radial derivative of the energy-momentum tensor, $\partial_r T[f]$. This is achieved because the projection of $V_+^{\text{mod}}$ onto the base manifold spans the tangent space of the Schwarzschild exterior, allowing the estimation of all derivatives of $T[f]$.

• Decay Estimates:

  • Inverse Polynomial Decay: Using the $r^p$-weighted method, the authors establish that the energy flux decays as $\langle \tau \rangle^{-p}$ for any $p \in \mathbb{N}$, provided the initial data has sufficient regularity (Corollary 1.1).
  • Exponential Decay: If the initial data decays exponentially, the energy flux decays exponentially (Corollary 1.4).
  • Pointwise Decay: The paper derives sharp pointwise decay estimates for the components of $T[f]$ (Corollary 1.1 and Corollary 8.8).

• Compatibility with Non-Linear Stability Programs: The authors note that their method is compatible with the approaches used in [Mav24, DHRT22] for quasi-linear wave equations. The vector fields used commute with the Vlasov operator in a way that allows the energy norms to control the source term $T[f]$ for the Einstein equations, a crucial step for studying the non-linear stability of the Schwarzschild spacetime for the massless Einstein–Vlasov system.

Significance
The paper claims that these results represent a significant improvement over previous works [ABJ18, Big23, Wei24] by removing the relative degeneracy in the energy decay estimates. This removal is essential for:

1. Deriving exponential decay via purely energy-based methods, rather than relying on super-polynomial decay or more complex hierarchies.
2. Controlling the radial derivatives of the energy-momentum tensor, which was previously an open difficulty due to the lack of a radial symmetry generator in Schwarzschild spacetime.
3. Providing a framework that is directly applicable to the study of the non-linear stability of the Schwarzschild spacetime for the massless Einstein–Vlasov system, bridging the gap between linear kinetic theory and non-linear geometric analysis.

The authors emphasize that while the results are established for the linear massless Vlasov equation on a fixed Schwarzschild background, the techniques are designed to be robust enough to be integrated into future non-linear stability proofs.