Hydrodynamic limit of the Vlasov-Poisson-Boltzmann system for gas mixture

In this paper, we study the hydrodynamic limit of the Vlasov-Poisson-Boltzmann system for a gas mixture in the whole space $(x \in \mathbb{R}^3)$ with the potential range of $γ\in\left(-3, 1\right]$. Using the method of Hilbert expansion, we first derive a bi-Maxwellian determined by the Euler-Poisson system of two fluids. To justify the convergence of the solution rigorously as the Knudsen number tends to zero, we sequentially calculate the first $2k-1$ terms of the expansion series $(k \geq 6)$, and then truncate it, and express the solution as the sum of these first $2k-1$ terms and a remainder term. Within the framework of the $L_{x,v}^2-W_{x,v}^{1,\infty}$ interplay established by Guo and Jang \cite{[ininp]Guo2010CMP}, we construct a new weight function to estimate the remainder term in four different cases regarding the potential $γ$. Here, the particle masses $m^A, m^B > 0$ and their charges $e^A, e^B$ can be given arbitrarily. This causes the collision operator to exhibit asymmetric effects ($m^A \neq m^B$), rendering the system of equations impossible to decouple. So, it adds difficulties to both $L^2$, $L^{\infty}$ estimates for the remainder. Therefore, we adopt the framework of vector-valued functions and analyze the velocity decay rate of the operator $K_{M,2,w}^{α,c}$ to eliminate the singularity induced by small parameters in characteristic line iterations. Our results show that the validity time of the solution is $O(\varepsilon^{-y})$, where $y$ is $-\frac{2k-3}{2(2k-1)}$ when $-1 \leq γ\leq 1$, and it becomes $-\frac{2k-3}{(1-γ)(2k-1)}$, when $-3 < γ< -1$. These results possess strong physical realism and can be applied to analyze gas flow dynamics in the daytime ionosphere at high altitudes above the Earth.

💡 Research Summary

The paper addresses the rigorous derivation of the hydrodynamic limit for a Vlasov‑Poisson‑Boltzmann (VPB) system describing a binary gas mixture with arbitrary particle masses and charges. The authors consider the whole space ℝ³ and a collision kernel whose kinetic part scales like |v−v*|^γ with γ∈(‑3,1]. The Knudsen number is denoted by ε (ε≪1), and the goal is to prove that as ε→0 the kinetic description converges to a macroscopic two‑fluid Euler‑Poisson system.

The analysis begins with a Hilbert expansion of the distribution functions F^α_ε (α=A,B) and the electric potential ϕ_ε up to order 2k‑1, where k≥6. The leading order term (ε^{‑1}) forces the collision operator to vanish, yielding a bi‑Maxwellian equilibrium µ_Aµ_B. The macroscopic fields (densities n_α, common bulk velocity u, and temperature θ) satisfy the Euler‑Poisson equations for two interacting fluids. The authors then derive recursive equations for the higher‑order coefficients F_i^α (i=0,…,2k‑1) and the remainder (F_R^α, ϕ_R).

A major difficulty stems from the asymmetry of the collision operator when the masses differ (m^A≠m^B). This prevents the decoupling of the system into scalar equations and invalidates standard L²‑based energy methods. To overcome this, the authors work in a vector‑valued function setting and employ the L²_x,v–W^{1,∞}_x,v framework introduced by Guo and Jang (2010). Within this framework they construct a new velocity‑dependent weight function w_γ(v) that adapts to four regimes of the potential exponent γ: hard sphere (γ=1), hard potentials/Maxwell molecules (0≤γ<1), moderately soft potentials (‑1≤γ<0), and very soft potentials (‑3<γ<‑1). The weight is designed to compensate for the growth of the collision frequency ν(v)∼⟨v⟩^γ and to control the L^∞ norm of the remainder.

A crucial technical component is the analysis of the linearized operator K_{M,2,w}^{α,c}. Lemma 4.3 establishes a decay rate in velocity for this operator even when the masses are unequal, which eliminates singularities that would otherwise arise in characteristic‑line iterations involving small parameters. This decay is essential for closing the L^∞ estimates.

The remainder equations are treated by combining an L²_ν energy estimate (with ν(v)≈⟨v⟩^γ) and a weighted L^∞ estimate. The authors obtain an inequality of the form

d/dt‖F_R‖{2,ν}² + λ‖(I−P)F_R‖{2,ν}² ≤ C(ε^{k‑1}‖F_R‖_{2,ν} + ε^{2k‑2}),

where P projects onto the collision invariants. Grönwall’s lemma yields uniform bounds for both the weighted L² and L^∞ norms of the remainder, showing that the remainder stays of order ε^{k‑1} over a time interval that depends on γ.

The final result is a quantitative validity time for the hydrodynamic approximation:

• For –1 ≤ γ ≤ 1, the solution exists up to T_γ = O(ε^{‑(2k‑3)/