Local Existence and Finite-Time Singularity Formation in the Vlasov-Poisson-Isotropic Landau System

1. Problem Statement

The authors investigate the Vlasov-Poisson-Isotropic Landau (VPiL) system in $\mathbb{R}^3$, which models a collisional plasma or self-gravitating system where particles interact via long-range forces and collisions. The system is given by:
$$\begin{aligned} \partial_t f + v \cdot \nabla_x f + E_f \cdot \nabla_v f &= Q_{\text{iso}}(f, f), \\ E_f &= (\mp) \frac{x}{4\pi|x|^3} \star_x \rho_f, \\ \rho_f(t, x) &= \int_{\mathbb{R}^3} f(t, x, v) \, dv, \end{aligned}$$
where the collision operator is the Isotropic Landau operator introduced by Krieger and Strain:
$$Q_{\text{iso}}(f, f) = (-\Delta_v)^{-1}f \Delta_v f + f^2.$$
Here, $(-\Delta_v)^{-1}f$ represents the Newtonian potential of $f$ in velocity space. The sign $(\mp)$ corresponds to the plasma (repulsive, $-$) and gravitational (attractive, $+$) cases, respectively.

Key Distinction: Unlike the classical Vlasov-Poisson-Landau system, the total energy in the VPiL system is not conserved. In fact, the energy is non-decreasing (and strictly increasing if $f \geq 0$), which suggests that the system may amplify concentration rather than smooth out singularities. This lack of a dissipative energy law poses significant challenges for global existence theory.

2. Methodology

The paper employs a combination of parabolic PDE theory, iterative approximation schemes, and moment analysis to address two main objectives: local existence and finite-time singularity formation.

A. Local Existence Theory

To establish the existence of weak solutions for small initial data, the authors utilize a method of continuity applied to a suitably designed approximation scheme.

1. Regularization: They introduce artificial diffusion ($\varepsilon \Delta_{x,v}$), spatial/velocity cutoffs ($\chi(\varepsilon x), \chi(\varepsilon v)$), and time mollification to ensure uniform ellipticity and coefficient regularity.
2. Function Spaces: The analysis is conducted in weighted parabolic Sobolev spaces ($W^{1,2}_{p,q,w}$) and weighted $L^\infty$ spaces. A crucial feature is the use of the transport-invariant weight $\langle x - vt \rangle^m \langle v \rangle^m$ to control the decay of solutions in phase space, overcoming the non-invariance of standard weights under the transport operator $v \cdot \nabla_x$.
3. Iterative Scheme: They construct a sequence of approximate solutions $\{f_n\}$ solving linearized problems.
4. A Priori Estimates: Using Schauder estimates and $L^p$ estimates for parabolic equations with Muckenhoupt weights, they derive uniform bounds on the sequence. A key step involves proving that the sequence remains bounded in a weighted $L^\infty$ norm, ensuring the nonlinear terms (electric field and collision operator) are well-defined.
5. Limit Passage: By passing to the limit using weak-* compactness and dominated convergence, they establish the existence of a non-negative weak solution.

B. Finite-Time Singularity Formation

For the gravitationally attractive case, the authors investigate conditions under which solutions blow up in finite time.

1. Mixed-Moment Argument: They analyze the evolution of the second spatial moment (moment of inertia):
   $$I(t) = \frac{1}{2} \iint_{\mathbb{R}^6} |x|^2 f(t, x, v) \, dv dx.$$
2. Derivatives of $I(t)$:
   • $I'(t) = \iint (x \cdot v) f \, dv dx$.
   • $I''(t) = 2KE(t) - EE(t)$, where $KE$ is kinetic energy and $EE$ is gravitational field energy.
   • $I'''(t) = \frac{d}{dt}(2KE - EE) = 4 \iint (-\Delta_v)^{-1}f f \, dv dx \geq 0$.
     Unlike the classical Landau equation where $I'''(t)=0$, the VPiL system has a positive third derivative due to the energy growth.
3. Polynomial Bound: They derive an upper bound for $I(t)$ in the form of a cubic polynomial:
   $$I(t) \leq I(0) + I'(0)t + \frac{1}{2}(2KE(0) - EE(0))t^2 + C t^3.$$
4. Blow-up Mechanism: By assuming the initial gravitational field energy significantly exceeds the kinetic energy ($EE(0) > 2KE(0)$) and the initial mixed moment is negative ($I'(0) < 0$), they show that the negative quadratic term dominates the small cubic term for a finite time interval. This forces the upper bound of $I(t)$ to become negative.
5. Contradiction: Since $f \geq 0$ implies $I(t) \geq 0$, the solution cannot exist globally while maintaining non-negativity and regularity. This proves finite-time singularity formation.

3. Key Contributions and Results

Result 1: Local Existence of Weak Solutions

• Theorem 2.3: For initial data $f_{\text{in}}$ in a weighted $C^2$ space ($C^2_{m,0}$) satisfying a smallness condition, there exists a time $T > 0$ (depending on the size of the data) such that a non-negative weak solution exists in the space $X^m_T$.
• The solution exhibits polynomial decay in both position and velocity variables.
• This result fills a gap in the mathematical theory of collisional plasmas, specifically for the near-vacuum regime where standard conditional regularity theories (requiring uniform lower bounds on density) fail.

Result 2: Finite-Time Singularity Formation

• Theorem 2.7: In the gravitational case, if the initial data satisfies:
  1. $I'(0) < 0$ (negative initial mixed moment),
  2. $EE(0) - 2KE(0)$ is sufficiently large (gravitational attraction dominates kinetic dispersion),
  3. The energy gap is large relative to the diffusive effect of the collision operator.
     Then, the maximal existence time $t_{\text{max}}$ is finite. The solution loses non-negativity or regularity at $t_{\text{max}}$.
• Corollary 2.9 (Mass Collapse): If the solution can be continued as a measure-valued solution up to $t_{\text{max}}$, the mass collapses to a single point in physical space ($x=0$) at that time. The limiting measure is singular with respect to the Lebesgue measure.

Result 3: Energy Monotonicity

• The paper rigorously establishes that the total energy of the VPiL system is strictly increasing (unless $f=0$), contrasting sharply with the energy conservation in the classical Vlasov-Poisson-Landau system. This property is central to the singularity mechanism.

4. Significance

• Novelty of the Model: This is the first mathematical study of the Vlasov-Poisson-Isotropic Landau system. The isotropic Landau operator is a simplified model of the full Landau operator, retaining the key feature of velocity-dependent diffusion governed by an inverse Laplacian but neglecting directional degeneracy.
• Breakdown of Global Existence: The paper provides a rigorous mechanism for finite-time blow-up in a kinetic equation with collisions. While collisions usually promote smoothing, the specific structure of the isotropic operator combined with strong gravitational attraction leads to a "runaway" concentration effect where diffusion cannot prevent collapse.
• Methodological Advances: The work introduces robust techniques for handling transport-invariant weights in parabolic Sobolev spaces and extends the mixed-moment argument to systems with non-conservative energy laws (where the third derivative of the moment is non-zero).
• Physical Implications: The results suggest that in certain regimes of self-gravitating systems with specific collisional models, the formation of singularities (black hole-like collapse in a kinetic context) is inevitable, even in the presence of collisional diffusion, provided the initial gravitational energy is sufficiently high.

In summary, the paper establishes a local well-posedness theory for a new kinetic system and demonstrates that, under attractive gravitational forces and specific energy conditions, this system inevitably develops finite-time singularities characterized by the collapse of mass to a point.