Asymptotic behavior of small solutions to the Vlasov--Klein--Gordon system in high dimensions

1. Problem Statement

The paper investigates the global existence and asymptotic behavior of small solutions to the Vlasov–Klein–Gordon (VKG) system in spatial dimensions $n \ge 4$. The system describes a collisionless ensemble of relativistic particles interacting via a self-consistent massive scalar field.

The system is given by:
$$\begin{aligned} \partial_t f + \hat{v} \cdot \nabla_x f - \nabla_x \phi \cdot \nabla_v f &= 0, \\ (-\partial_t^2 + \Delta_x - 1)\phi &= \int_{\mathbb{R}^n} f \, dv, \end{aligned}$$
where $f(t,x,v)$ is the distribution function, $\phi(t,x)$ is the Klein–Gordon field, and $\hat{v} = v/\sqrt{1+|v|^2}$.

The Core Difficulty:
While the Vlasov–Maxwell (VM) system (massless field) admits global small solutions via the standard argument of Glassey and Strauss, this approach fails for the VKG system.

• Massiveness: The Klein–Gordon field has mass, violating the strong Huygens' principle. Consequently, decay estimates based on the $(t-r)^{-1}$ factor (crucial for the VM system inside the light cone) are irrelevant.
• Decay Rates: Although the massive field decays faster ($t^{-3/2}$) than the massless field ($t^{-1}$), this rate alone is insufficient to close the energy estimates for the Vlasov equation without additional structural tools. Previous works required artificial damping terms to achieve necessary decay; this paper addresses the undamped system.

2. Methodology

The author employs the vector field method, originally introduced by Klainerman for wave equations and later adapted to kinetic equations (Vlasov–Poisson, Vlasov–Nordström, Vlasov–Maxwell).

Key Technical Components:

1. Hyperboloidal Foliation:
   Instead of standard constant-time slices, the analysis is conducted on the hyperboloidal foliation of Minkowski spacetime defined by $\tau = \sqrt{t^2 - |x|^2} \ge 1$. This geometry is particularly effective for capturing the decay of massive fields and handling the interaction between the transport and wave equations.

2. Vector Fields and Commutation:
   The method utilizes the set of Killing vector fields $\mathcal{P}$ (translations, rotations, and boosts) and their complete lifts $\hat{\mathcal{P}}$ to the phase space $(x, v)$.

   • The commutators $[T_\phi, \hat{Z}]$ (where $T_\phi$ is the perturbed transport operator) are computed explicitly.
   • The commutation formulas generate source terms involving polynomials of velocity components ($v/v_0$) and derivatives of the field $\phi$.
   • A crucial step involves decomposing the error terms into "good" and "bad" parts, utilizing the specific algebraic structure of the commutators to control the growth of high-order derivatives.

3. Energy Norms:
   The paper defines specific energy norms adapted to the hyperboloidal foliation:

   • $E_N(\phi)$: Energy for the Klein–Gordon field involving $L^2$ norms of derivatives.
   • $\hat{E}_{N,1}(f)$ and $\hat{E}_{N+n}(f)$: Energy norms for the distribution function, utilizing $L^1$-type estimates in velocity space weighted by $v_0$.
   • The distinction between low-order and high-order derivatives is managed by splitting the energy norms (e.g., derivatives up to $\lfloor N/2 \rfloor$ are weighted by $v_0$).

4. Bootstrap Argument:
   The proof relies on a continuity (bootstrap) argument:

   • Assume the solution remains small in the defined energy norms up to time $T$.
   • Derive improved estimates for the field $\phi$ and the distribution $f$ using the commuted equations and Klainerman–Sobolev inequalities.
   • Show that the improved bounds are strictly better than the bootstrap assumptions, allowing the extension of the solution to $T = \infty$.

3. Key Contributions

• Resolution of the Undamped Case: The paper successfully establishes the asymptotic stability of small solutions for the undamped VKG system in dimensions $n \ge 4$, removing the need for artificial damping terms used in prior studies (e.g., Ha and Lee).
• Adaptation of Vector Field Method: It extends the vector field method to the VKG system, overcoming the obstruction caused by the violation of the strong Huygens' principle. This is achieved by leveraging the hyperboloidal foliation to derive necessary decay rates without relying on $(t-r)^{-1}$ estimates.
• Handling High Dimensions: The analysis is restricted to $n \ge 4$. The author notes that the case $n=3$ would require "modified vector fields" (a technique used in Vlasov–Maxwell studies) due to weaker decay rates, which is left as future work.

4. Main Results (Theorem 1.1)

For initial data on the unit hyperboloid $H_1$ that is sufficiently small (measured by specific energy norms), the system admits a global solution with the following asymptotic properties for $\tau \ge 1$:

1. Boundedness of Field Energy: The energy of the Klein–Gordon field remains uniformly bounded:
   $$E_N(\phi)(\tau) \lesssim \varepsilon.$$
2. Growth of Distribution Energy: The energy of the distribution function grows at most polynomially (or remains bounded):
   $$\hat{E}_{N,1}(f)(\tau) \lesssim \varepsilon \tau^\delta,$$
   where $\delta = 0$ for $n \ge 5$ and $\delta = \varepsilon^{1/4}$ for $n=4$.
3. Pointwise Decay Estimates:
   • Distribution Function: For high-order derivatives ($|A| \le N-n$), the velocity integral decays as:
     $$\int_{\mathbb{R}^n} |\hat{Z}^A f| \, dv \lesssim \frac{\varepsilon}{t^{n-1+\delta/2}(t-|x|)^{(1-\delta)/2}}.$$
   • Klein–Gordon Field: The field and its derivatives decay as:
     $$|Z^A \phi| \lesssim \frac{\varepsilon^{1/2}}{t^{n/2}}, \quad |\partial Z^A \phi| \lesssim \frac{\varepsilon^{1/2}}{t^{(n-1)/2}(t-|x|)^{1/2}}.$$

5. Significance

• Mathematical Rigor: This work provides a rigorous proof of global existence for small data in the VKG system, a problem that has been challenging due to the interplay between the massive wave equation and the transport equation.
• Methodological Advancement: It demonstrates the robustness of the vector field method combined with hyperboloidal foliation for kinetic-wave systems where the strong Huygens' principle fails.
• Foundation for Future Work: By solving the problem for $n \ge 4$, it sets the stage for tackling the more difficult 3-dimensional case, which is physically most relevant but requires more sophisticated modifications to the vector field method.
• Physical Insight: The results confirm that the mass of the Klein–Gordon field, while complicating the decay structure, does not prevent global stability for small perturbations, provided the dimension is high enough to ensure sufficient dispersion.