Incompressible Navier-Stokes-Fourier limit from a nonlinear quantum Fokker-Planck equation

Technical Summary: Incompressible Navier–Stokes–Fourier Limit from a Nonlinear Quantum Fokker–Planck Equation

Problem Statement
This paper addresses the rigorous derivation of the incompressible Navier–Stokes–Fourier (NSF) system from a nonlinear quantum Fokker–Planck equation. The governing kinetic equation describes the evolution of a distribution function $f(t, x, p)$ incorporating Bose–Einstein or Fermi–Dirac statistics through a nonlinear collision term of the form $f(1 + \hbar\kappa f)$, where $\kappa \in \{1, -1\}$ and $\hbar > 0$ is a fixed quantum parameter.

The specific challenge lies in the diffusive scaling limit where the Knudsen number $\epsilon \to 0$, while the quantum parameter $\hbar$ remains fixed. Unlike classical hydrodynamic limits where the equilibrium state approaches a Maxwellian distribution, this regime requires identifying a limiting fluid system that retains the full effects of quantum statistics (Bose enhancement or Pauli exclusion) in its transport coefficients and normalization constants. The model features a self-consistent collision structure where the local density, bulk velocity, and temperature are determined by nonlinear quantum-weighted moments of the distribution, rather than being prescribed.

Methodology
The authors employ a perturbative macro–micro decomposition around a global quantum equilibrium $F_\hbar$. The analysis proceeds through the following stages:

1. Scaling and Perturbation: The distribution function is expanded as $f^\epsilon = F_\hbar + \epsilon\sqrt{\mu_\hbar}g^\epsilon$, where $\mu_\hbar$ is a quantum weight derived from $F_\hbar$. This transforms the nonlinear quantum Fokker–Planck equation into a rescaled equation for the perturbation $g^\epsilon$:
   $$\epsilon\partial_t g^\epsilon + p \cdot \nabla_x g^\epsilon = \frac{1}{\epsilon} \mathcal{L}g^\epsilon + \Gamma^\epsilon(g^\epsilon)$$
   Here, $\mathcal{L}$ is the linearized collision operator, and $\Gamma^\epsilon$ contains scaled nonlinear remainders arising from the self-consistent macroscopic fields.

2. Uniform Estimates: The authors establish global energy estimates uniform in $\epsilon$. By adapting the perturbative well-posedness theory from previous work (Choi et al., [6]) to the singular diffusive scaling, they prove that the microscopic component $(I-\mathcal{P})g^\epsilon$ (where $\mathcal{P}$ is the projection onto the null space of $\mathcal{L}$) relaxes strongly to zero at a rate of order $\epsilon$. This ensures that any weak limit of the solution belongs to the quantum macroscopic space $\mathcal{N} = \ker \mathcal{L}$.

3. Macroscopic Constraints and Compactness:

   • Conservation Laws: By taking velocity moments of the rescaled equation, the authors derive local conservation laws for mass, momentum, and energy.
   • Incompressibility and Boussinesq Relation: Passing to the limit $\epsilon \to 0$ in the momentum and mass equations yields the incompressibility condition $\nabla_x \cdot u = 0$ and the Boussinesq relation $\varrho + \vartheta = 0$, linking density and temperature fluctuations.
   • Strong Compactness: To handle nonlinear terms, strong compactness is required. The authors identify a "quantum-adapted" thermal variable $q^\epsilon = \vartheta^\epsilon - \frac{\beta}{1+\beta}\varrho^\epsilon$ (where $\beta$ is a constant depending on quantum moments) and the solenoidal velocity component. Using the Aubin–Lions lemma on these variables, they establish strong convergence in local Sobolev spaces.
   • Acoustic Modes: The acoustic variables (density/temperature sum and irrotational velocity) satisfy a fast wave system. The authors employ a local dispersive estimate for the acoustic group to show these modes vanish strongly in local spaces, preventing them from contributing to the limiting nonlinear convection.

4. Identification of Fluxes: The limiting viscous stress tensor and heat flux are identified by solving auxiliary microscopic equations $\mathcal{L}\tilde{A} = A$ and $\mathcal{L}\tilde{B} = B$ on the orthogonal complement of the null space. The transport coefficients (quantum viscosity $\nu_\hbar$ and thermal diffusivity $\kappa_\hbar$) are defined via quadratic forms involving these auxiliary functions and the quantum equilibrium. The nonlinear convection terms are derived by expanding the local quantum equilibrium manifold to second order.

Key Results
The main theorem (Theorem 1.1) establishes that, for small initial perturbations in high-order Sobolev spaces, the global solution $g^\epsilon$ converges (up to a subsequence) to a limit $g$ of the form:
$$g(t, x, p) = \left( \frac{u(t, x) \cdot p}{m_2} + \vartheta(t, x) \left[ \frac{1}{m_0\beta}\left(\frac{|p|^2}{\phi} - 1\right) - \frac{1}{m_0} \right] \right) \sqrt{\mu_\hbar}(p)$$
where $u$ and $\vartheta$ satisfy the incompressible Navier–Stokes–Fourier system:
$$\begin{cases} \partial_t u + \frac{1}{m_2} u \cdot \nabla_x u + \nabla_x p = \nu_\hbar \Delta_x u, \\ \nabla_x \cdot u = 0, \\ \partial_t \vartheta + \frac{1}{m_2} u \cdot \nabla_x \vartheta = \kappa_\hbar \Delta_x \vartheta. \end{cases}$$
The initial data for the limit are determined by the projection of the initial kinetic data. The coefficients $\nu_\hbar$ and $\kappa_\hbar$ are strictly positive and depend explicitly on the fixed quantum parameter $\hbar$ and the equilibrium $F_\hbar$.

Significance and Claims
The paper claims to provide the first rigorous incompressible Navier–Stokes–Fourier limit derived from a nonlinear quantum Fokker–Planck equation with self-consistent macroscopic fields.

• Quantum vs. Semiclassical: The authors emphasize that their limit is distinct from a semiclassical limit. By keeping $\hbar$ fixed as $\epsilon \to 0$, the resulting fluid system retains the full Bose–Einstein or Fermi–Dirac corrections in its transport coefficients and normalization constants, rather than approximating them with classical Maxwellian counterparts.
• Structural Complexity: The work addresses the specific difficulties introduced by the nonlinear, self-consistent nature of the collision frequency and macroscopic fields, which generate additional terms in the macroscopic fluxes that must be carefully separated into transport and convection parts.
• Extension of Theory: This result extends the known hydrodynamic limits for classical nonlinear Vlasov–Fokker–Planck models to the quantum regime, providing a counterpart that accounts for quantum statistical effects at the level of the limiting fluid equations.

The paper does not propose new applications or experimental setups but focuses on the mathematical justification of the hydrodynamic regime for this specific class of quantum kinetic models.