Future global stability of Maxwell-Jüttner equilibria and vacuum for the massless Boltzmann equation on FLRW spacetimes

This paper establishes the future global-in-time existence and uniqueness of small perturbations for both Maxwell-Jüttner equilibria and vacuum solutions of the massless Boltzmann equation on decelerating FLRW spacetimes with $\mathbb{T}^3$ topology, covering hard ball interactions for all expansion rates $\mathfrak{q} \in [0,1]$ and vacuum stability for $\mathfrak{q} > 1/3$.

The Gist

Imagine the universe as a giant, expanding balloon. Inside this balloon, there are countless tiny, invisible particles zooming around, bouncing off one another like hyperactive billiard balls. This paper is a mathematical study of how these particles behave when the balloon is inflating, specifically focusing on two scenarios: when the particles are already in a calm, balanced state, and when there are almost no particles at all.

Here is a breakdown of the paper's findings using simple analogies:

The Setting: The Expanding Balloon

The authors are looking at a model of the universe called FLRW spacetime. Think of this as a 3D grid (like a video game world that wraps around on itself, called a torus) that is stretching over time.

• The Scale Factor ($t^q$): The universe isn't just expanding; it's expanding at different speeds depending on a number called $q$.
  • If $q$ is small, the universe expands slowly (decelerating).
  • If $q$ is large (up to 1), it expands faster (linearly).
  • The "time" in this story starts at the Big Bang ($t=0$) and moves forward.

The Particles: Massless Billiard Balls

The particles being studied are massless (like photons of light) and they collide with each other. The math used to describe their collisions is called the Boltzmann equation.

• The "Hard Ball" Rule: The authors assume these particles interact like hard spheres (or hard balls). When they hit, they bounce off instantly. This is a specific, simplified way of modeling how they crash into each other.

Scenario 1: The Calm State (Maxwell–Jüttner Equilibrium)

Imagine the particles are dancing in a very specific, organized pattern. In a static room, this pattern would stay the same forever. But because the universe (the balloon) is expanding, this "dance" has to change shape to keep up.

• The Equilibrium: The authors found a special, non-stationary dance routine (called a Maxwell–Jüttner equilibrium) that the particles naturally fall into as the universe expands. It's like a dance that slowly slows down and spreads out as the room gets bigger.
• The Stability Test: The big question was: If you nudge this dance slightly (add a little chaos), will it eventually settle back into the rhythm, or will it spiral out of control?
• The Result:
  • It's Stable: For small nudges, the system always returns to the rhythm. The particles don't go wild; they find their way back to the "equilibrium dance."
  • The Speed of Recovery: How fast they settle down depends on how fast the universe is expanding ($q$).
    • Slow Expansion ($q$ is small): The particles settle down very fast. In fact, they settle down faster than any standard polynomial speed (super-polynomial decay). It's like a shock absorber that works incredibly well.
    • Fast Expansion ($q$ is large): The universe is stretching so fast that it actually fights the particles' ability to calm down. The "friction" from collisions isn't strong enough to overcome the stretching. The particles still settle down, but much slower (polynomial decay).
    • The Tipping Point ($q = 1/3$): There is a magic number, $1/3$. Below this, the universe's expansion is slow enough that particle collisions act like a strong brake, quickly stopping the chaos. Above this, the expansion is so strong it weakens the braking effect of the collisions.

Scenario 2: The Empty Room (Vacuum Solution)

Now, imagine the room is almost empty. There are very few particles.

• The Question: If you start with just a few particles in this expanding universe, will they eventually disappear (decay to zero), or will they clump together and cause trouble?
• The Result:
  • If the universe is expanding fast enough ($q > 1/3$), the particles will naturally spread out and fade away until the room is effectively empty (the vacuum is stable). The expansion acts like a giant fan blowing the particles apart so they never collide enough to cause a problem.
  • If the expansion is too slow ($q \le 1/3$), the authors couldn't prove this stability with their current methods. The particles might hang around too long and interact in ways that are hard to predict.

The "Secret Sauce" of the Math

The authors had to invent new mathematical tools to solve this.

• The Problem: Standard math tools for particle physics assume the room is a fixed size. Here, the room is stretching.
• The Solution: They created a "time-normalized" view. Imagine watching the particles through a camera that zooms out at the exact same rate the universe expands. In this zoomed-out view, the particles look like they are in a normal, static room, making it possible to apply standard stability tests.
• The Energy Method: They tracked the "energy" of the chaos. They proved that even though the universe is stretching, the energy of the disturbance (the nudge) eventually drains away, either through the particles hitting each other (dissipation) or just being stretched out by the universe (dispersion).

Summary

In simple terms, this paper proves that:

1. Order wins: Even in an expanding universe, if particles are close to a calm state, they will stay calm.
2. Expansion matters: How fast the universe expands determines how quickly the particles calm down. If the universe expands too fast, it weakens the natural "braking" effect of particle collisions.
3. Empty is safe: If the universe is expanding fast enough, a nearly empty universe will stay empty and stable.

This is a theoretical proof about the long-term behavior of gas particles in a cosmological setting, ensuring that our mathematical models of the universe don't break down over time.

Technical Summary

Technical Summary: Future Global Stability of Maxwell–Jüttner Equilibria and Vacuum for the Massless Boltzmann Equation on FLRW Spacetimes

Problem Statement
This work investigates the long-time dynamics of the general relativistic massless Boltzmann equation on Friedmann–Lemaître–Robertson–Walker (FLRW) spacetimes with spatial topology $T^3$. The study focuses on the future global-in-time stability of two specific states:

1. Maxwell–Jüttner Equilibria: Non-stationary equilibrium solutions of the form $J(t^2 q p) = \exp(-|t^2 q p|)$, where the scale factor is $a(t) = t^q$ with $q \in [0, 1]$. These equilibria decay in time for $q > 0$.
2. The Vacuum Solution: The trivial solution $F \equiv 0$.

The analysis covers both linearly expanding ($q=1$) and decelerated expanding ($0 < q < 1$) regimes, as well as the non-expanding case ($q=0$). The collision operator is modeled using a hard-ball interaction kernel (constant scattering cross-section), though the authors discuss potential extensions to more general kernels.

Methodology
The authors employ a perturbative approach, analyzing small deviations from the equilibrium and vacuum states. The core methodology involves:

• Reformulation: The distribution function $F$ is decomposed into an equilibrium part and a perturbation $f$. For the equilibrium case, $F = J + \sqrt{J}f$, leading to a linearized equation involving a linear operator $L$ and a nonlinear term $\Gamma$. For the vacuum case, $F = \sqrt{J}f$, resulting in a purely nonlinear equation.
• Macro-Micro Decomposition: The perturbation $f$ is decomposed into a macroscopic part $Pf$ (projection onto the null space of the linearized operator $L$, corresponding to mass, momentum, and energy conservation) and a microscopic part $(I-P)f$.
• Time-Rescaled Orthonormal Basis: A critical innovation is the use of a time-normalized orthonormal basis for the null space of $L$. Unlike previous works on fixed backgrounds, the basis vectors here depend on time due to the cosmological expansion, requiring a new derivation of macroscopic equations.
• Energy Methods: The proofs rely on Guo's energy method combined with low-regularity Wiener algebra techniques (using Fourier transforms on the torus). The authors define specific energy norms weighted by powers of time, such as $\|t^{3q}f\|_{L^1_k L^\infty_T L^2_p}$, to capture the decay rates.
• Coercivity and Compactness: The linearized operator $L$ is shown to be coercive modulo its null space. This is achieved by decomposing the integral operator $K$ (part of $L = \nu - K$) into a compact part and a small part, utilizing the specific structure of the hard-ball kernel on FLRW backgrounds.
• Macroscopic Estimates: The authors derive a system of macroscopic equations for the coefficients of the projection $Pf$. By exploiting the vanishing of spatial averages (due to conservation laws) and the upper-triangular structure of the system, they obtain decay estimates for the macroscopic part.

Key Contributions and Results

1. Global Stability of Maxwell–Jüttner Equilibria ($0 \le q \le 1$):

   • The authors prove the future global-in-time existence and uniqueness of small perturbations around the Maxwell–Jüttner equilibrium for hard-ball interactions without symmetry assumptions.
   • Decay Rates: The decay of perturbations depends critically on the expansion rate $q$:
     • For $0 \le q < 1/3$: The perturbation decays at a super-polynomial rate of $t^{-3q} \exp(-t^{1-3q})$.
     • For $q = 1/3$: The decay is $t^{-3q - c}$ for a uniform constant $c > 0$.
     • For $1/3 < q \le 1$: The decay is polynomial, specifically $t^{-3q}$.
   • Dissipation Threshold: The value $q = 1/3$ is identified as a threshold. Below this threshold, the collisional dissipation dominates the cosmological expansion, leading to faster decay. Above this threshold, the expansion weakens the dissipative effect, and the decay rate is determined solely by the expansion factor.

2. Global Stability of the Vacuum Solution ($1/3 < q \le 1$):

   • The paper establishes the future global-in-time existence and uniqueness of small perturbations around the vacuum solution for the massless Boltzmann equation.
   • This result holds for $q > 1/3$. The proof relies on the time decay induced by cosmological expansion to control the nonlinear term. The restriction $q > 1/3$ arises because the decay rate $t^{-3q}$ must be integrable to close the energy estimates for the nonlinear term in the absence of a linearized collision operator (which is trivial near vacuum).

3. Weighted Norms and Regularity:

   • The results extend to weighted norms involving powers of momentum, demonstrating the propagation of spatial regularity and the preservation of moments.
   • In the non-expanding case ($q=0$), the authors recover exponential decay rates.

Significance and Claims
The authors claim this work provides the first future global-in-time stability result for the Boltzmann equation in the presence of both spatial dependence (on a torus) and a non-trivial gravitational field (FLRW spacetime).

• Novelty in Kinetic Theory: The work addresses the interplay between cosmological expansion and particle collisions. It demonstrates that expansion can fundamentally alter the decay rates of kinetic systems, acting as a "weakening" mechanism for dissipation when the expansion rate is sufficiently high ($q > 1/3$).
• Mathematical Innovation: The development of a time-rescaled orthonormal basis and the derivation of macroscopic equations adapted to the expanding background are presented as essential tools for obtaining precise large-time estimates in this geometric setting.
• Threshold Phenomenon: The identification of $q=1/3$ as a critical threshold for the dissipation mechanism in the massless case is a central finding. The authors conjecture that for more general collision kernels (e.g., soft potentials), this threshold shifts to $(3+a)q = 1$, where $a$ characterizes the kernel's growth.

The paper concludes that while the Einstein–Boltzmann system admits explicit FLRW solutions, the rigorous nonlinear stability of these solutions in the massless, non-homogeneous setting requires new analytical techniques to handle the time-dependent nature of the equilibrium and the competition between expansion and collisional dissipation.