Fuzzy Landau Equation

• The fuzzy Landau equation is a nonlocal kinetic model where spatial delocalisation of collisions is achieved via a convolution kernel, modifying classical Landau dynamics.
• Its GENERIC formulation and De Giorgi-type variational characterization enable rigorous proofs of global existence, entropy dissipation, and stability through explicit estimates.
• Recent research establishes strong well-posedness, uniqueness, and the grazing-collision limit from fuzzy Boltzmann to fuzzy Landau, highlighting advances in nonlocal transport-diffusion analysis.

Searching arXiv for papers on the fuzzy Landau equation and closely related variational/solvability results. The fuzzy Landau equation is a nonlocal kinetic equation in which the spatially inhomogeneous Landau collision operator is modified so that collisions are delocalised in the position variable through a nonnegative interaction kernel. For a one-particle density $f=f(t,x,v)$, the basic form is

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$

with

$$Q_{\sf fuz}(f,f) =\nabla_v\!\cdot\!\int_{\Omega\times\mathbb R^3} \kappa(x-x_*)\,a(v-v_*)\, \bigl[f_*\,\nabla_v f-f\,\nabla_{v_*}f_*\bigr]\, dx_*\,dv_*,$$

where $f_*=f(t,x_*,v_*)$ and $\Omega=\mathbb T^3$ or $\mathbb R^3$ in the three-dimensional formulation. The current theory places this equation in the GENERIC framework, derives a De Giorgi-type variational characterization, proves global $H$-solutions, establishes global-in-time existence and uniqueness of smooth solutions for moderately soft potentials in a specific regime, identifies new Fisher-information monotonicity, derives the grazing-collision limit from non-cutoff fuzzy Boltzmann equations, and proves conditional uniqueness of weak solutions through explicit $2$-Wasserstein stability estimates (Duong et al., 10 Apr 2025, Duong et al., 14 Jul 2025, Gualdani et al., 11 Jul 2025, Duong et al., 4 Dec 2025, Caja-Lopez, 29 Mar 2026).

1. Kinetic definition and the role of spatial delocalisation

The defining modification is the replacement of the spatially local collision rule of the classical inhomogeneous Landau equation by a delocalised interaction. In the classical operator, collisions are fully localized in $x$, corresponding formally to the limit $\kappa\to\delta_0$. In the fuzzy equation, the spatial “fuzziness” is encoded by a kernel $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$0 with $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$1; examples in the literature include smooth kernels and exponentially decaying kernels, while later formulations also allow bounded spatial interaction kernels and smooth Gaussian-type kernels (Duong et al., 10 Apr 2025, Duong et al., 4 Dec 2025).

The velocity interaction has the standard Landau tensor structure

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$2

For Coulomb-type interactions one takes $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$3 with $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$4 in the three-dimensional setting, while the solvability theory in general dimension treats either $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$5 with $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$6, $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$7, or the milder class

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$8

(Duong et al., 10 Apr 2025, Duong et al., 14 Jul 2025).

This spatial delocalisation changes the analytic character of the equation. Part II explicitly contrasts the fuzzy model with the classical inhomogeneous Landau equation, noting that fuzziness introduces a genuinely non-local transport–diffusion in $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$9 but improves compactness and coercivity in $$Q_{\sf fuz}(f,f) =\nabla_v\!\cdot\!\int_{\Omega\times\mathbb R^3} \kappa(x-x_*)\,a(v-v_*)\, \bigl[f_*\,\nabla_v f-f\,\nabla_{v_*}f_*\bigr]\, dx_*\,dv_*,$$0. The same source records two asymptotic regimes: $$Q_{\sf fuz}(f,f) =\nabla_v\!\cdot\!\int_{\Omega\times\mathbb R^3} \kappa(x-x_*)\,a(v-v_*)\, \bigl[f_*\,\nabla_v f-f\,\nabla_{v_*}f_*\bigr]\, dx_*\,dv_*,$$1 gives the classical inhomogeneous Landau equation, while $$Q_{\sf fuz}(f,f) =\nabla_v\!\cdot\!\int_{\Omega\times\mathbb R^3} \kappa(x-x_*)\,a(v-v_*)\, \bigl[f_*\,\nabla_v f-f\,\nabla_{v_*}f_*\bigr]\, dx_*\,dv_*,$$2 yields the homogeneous Landau equation (Duong et al., 14 Jul 2025).

2. GENERIC formulation and variational characterization

Duong and He recast the fuzzy Landau equation as a GENERIC system, with state variable $$Q_{\sf fuz}(f,f) =\nabla_v\!\cdot\!\int_{\Omega\times\mathbb R^3} \kappa(x-x_*)\,a(v-v_*)\, \bigl[f_*\,\nabla_v f-f\,\nabla_{v_*}f_*\bigr]\, dx_*\,dv_*,$$3 in

$$Q_{\sf fuz}(f,f) =\nabla_v\!\cdot\!\int_{\Omega\times\mathbb R^3} \kappa(x-x_*)\,a(v-v_*)\, \bigl[f_*\,\nabla_v f-f\,\nabla_{v_*}f_*\bigr]\, dx_*\,dv_*,$$4

The energy and entropy functionals are

$$Q_{\sf fuz}(f,f) =\nabla_v\!\cdot\!\int_{\Omega\times\mathbb R^3} \kappa(x-x_*)\,a(v-v_*)\, \bigl[f_*\,\nabla_v f-f\,\nabla_{v_*}f_*\bigr]\, dx_*\,dv_*,$$5

The reversible part is generated by the antisymmetric operator

$$Q_{\sf fuz}(f,f) =\nabla_v\!\cdot\!\int_{\Omega\times\mathbb R^3} \kappa(x-x_*)\,a(v-v_*)\, \bigl[f_*\,\nabla_v f-f\,\nabla_{v_*}f_*\bigr]\, dx_*\,dv_*,$$6

equivalently

$$Q_{\sf fuz}(f,f) =\nabla_v\!\cdot\!\int_{\Omega\times\mathbb R^3} \kappa(x-x_*)\,a(v-v_*)\, \bigl[f_*\,\nabla_v f-f\,\nabla_{v_*}f_*\bigr]\, dx_*\,dv_*,$$7

The irreversible part is encoded by a diffusion operator $$Q_{\sf fuz}(f,f) =\nabla_v\!\cdot\!\int_{\Omega\times\mathbb R^3} \kappa(x-x_*)\,a(v-v_*)\, \bigl[f_*\,\nabla_v f-f\,\nabla_{v_*}f_*\bigr]\, dx_*\,dv_*,$$8 built from the delocalised Landau tensor. Introducing the fuzzy Landau-gradient

$$Q_{\sf fuz}(f,f) =\nabla_v\!\cdot\!\int_{\Omega\times\mathbb R^3} \kappa(x-x_*)\,a(v-v_*)\, \bigl[f_*\,\nabla_v f-f\,\nabla_{v_*}f_*\bigr]\, dx_*\,dv_*,$$9

one has

$f_*=f(t,x_*,v_*)$0

The GENERIC equation

$f_*=f(t,x_*,v_*)$1

recovers exactly the fuzzy Landau evolution. The construction satisfies the standard degeneracies

$f_*=f(t,x_*,v_*)$2

together with antisymmetry of the Poisson bracket and formal Jacobi identity (Duong et al., 10 Apr 2025).

A second structural layer is a De Giorgi-type variational characterization. For a transport–grazing-rate pair $f_*=f(t,x_*,v_*)$3 satisfying

$f_*=f(t,x_*,v_*)$4

the action

$f_*=f(t,x_*,v_*)$5

and entropy dissipation

$f_*=f(t,x_*,v_*)$6

enter the functional

$f_*=f(t,x_*,v_*)$7

Equality holds if and only if

$f_*=f(t,x_*,v_*)$8

which is equivalent to $f_*=f(t,x_*,v_*)$9 being a weak, or “$\Omega=\mathbb T^3$0-”, solution of the fuzzy Landau equation. This formulation unifies the GENERIC flow and the $\Omega=\mathbb T^3$1-theorem in a single min-action principle (Duong et al., 10 Apr 2025).

3. Conservation laws, entropy production, and global $\Omega=\mathbb T^3$2-solutions

The GENERIC structure yields the fundamental balance laws. For smooth solutions,

$\Omega=\mathbb T^3$3

Entropy is nonincreasing: $\Omega=\mathbb T^3$4 Part I describes this as automatic enforcement of thermodynamic consistency (1st/2nd laws), correct conservation laws and entropy dissipation, together with a “Lyapunov structure” well suited to long-time asymptotics (Duong et al., 10 Apr 2025).

Part II develops the weak theory in a functional framework based on weighted $\Omega=\mathbb T^3$5 spaces, mixed $\Omega=\mathbb T^3$6 norms, entropy space $\Omega=\mathbb T^3$7, velocity moments $\Omega=\mathbb T^3$8, Fisher information

$\Omega=\mathbb T^3$9

and fuzzy entropy dissipation

$\mathbb R^3$0

For initial data

$\mathbb R^3$1

Duong and He define weak solutions and $\mathbb R^3$2-solutions through the integrated GENERIC identity

$\mathbb R^3$3

They prove existence of global-in-time $\mathbb R^3$4-solutions. For such solutions, mass and momentum are conserved, energy satisfies

$\mathbb R^3$5

the spatial second moment obeys

$\mathbb R^3$6

and the entropy inequality holds: $\mathbb R^3$7 (Duong et al., 14 Jul 2025).

The analytic mechanism behind this existence theory combines several ingredients recorded explicitly in Part II: the Desvillettes-type bound

$\mathbb R^3$8

a weighted Sobolev embedding

$\mathbb R^3$9

Hardy–Littlewood–Sobolev estimates for soft potentials, coercivity of the diffusion matrix

$H$0

through lower bounds of the form

$H$1

or

$H$2

and finally approximation, parabolic regularisation in $H$3, and a fixed-point argument leading to smooth approximants, a priori bounds, compactness, and passage to the limit (Duong et al., 14 Jul 2025).

4. Strong well-posedness, ellipticity, and Fisher information

A complementary line of work studies a fuzzy variant on $H$4 with a smooth, everywhere-positive mollifier $H$5 in place of the spatial delta kernel. In this formulation,

$H$6

with

$H$7

For $H$8 and integer $H$9, the main theorem gives a unique strong solution provided the initial datum satisfies the stated $2$0, $2$1, moment, and weighted Sobolev assumptions. The solution obeys the weighted bounds

$2$2

and

$2$3

for all multi-indices with $2$4; mass, momentum, and energy are conserved (Gualdani et al., 11 Jul 2025).

The central a priori mechanism is a global ellipticity effect created by delocalisation. Defining

$2$5

the paper proves that there exists $2$6 such that

$2$7

for all $2$8. The source emphasizes that because $2$9 everywhere, any two particles “feel” a weak collisional interaction regardless of spatial separation, and that this forces an immediate gain of ellipticity in $x$0 even where $x$1 would be vacuum (Gualdani et al., 11 Jul 2025).

The same paper establishes additional Lyapunov structure through Fisher information. For $x$2 and $x$3, any smooth solution satisfies

$x$4

and

$x$5

The proof uses a doubled-variable representation with

$x$6

which converts the Landau operator into a diffusion-type operator on $x$7. In addition, for Maxwell molecules $x$8 the mixed Fisher functional

$x$9

satisfies $\kappa\to\delta_0$0 (Gualdani et al., 11 Jul 2025).

5. Grazing-collision limit and asymptotic connections

The fuzzy Landau equation arises as the grazing-collision limit of non-cutoff fuzzy Boltzmann equations. In the Boltzmann model, particles at $\kappa\to\delta_0$1 and $\kappa\to\delta_0$2 undergo delocalised elastic collisions with collision operator

$\kappa\to\delta_0$3

where

$\kappa\to\delta_0$4

The non-cutoff regime is described by an angular singularity

$\kappa\to\delta_0$5

To isolate grazing events one rescales

$\kappa\to\delta_0$6

so that $\kappa\to\delta_0$7 and $\kappa\to\delta_0$8 is fixed. Formally, as $\kappa\to\delta_0$9, the Boltzmann operator converges to the fuzzy Landau operator with

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$00

(Duong et al., 4 Dec 2025).

The variational theory persists across this limit. Duong, Golubkov, and He formulate the fuzzy Boltzmann equation through a non-quadratic dual dissipation pair $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$01, while the fuzzy Landau equation corresponds to a quadratic dissipation pair. The associated action functionals satisfy

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$02

with equality characterizing solutions in each case. Their main theorem states that if $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$03 are $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$04-solutions of the fuzzy Boltzmann equation under the stated moment–entropy bounds and uniform controls on $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$05, then

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$06

and the limit $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$07 is an $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$08-solution of the fuzzy Landau equation. The proof proceeds through compactness, velocity averaging and fractional diffusion from the angular singularity, together with a $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$09-liminf argument for the actions and dissipation functionals (Duong et al., 4 Dec 2025).

This asymptotic result aligns with the broader program already indicated in the variational paper: the framework is stated to connect smoothly to the grazing-collision limit from Boltzmann to Landau and to the homogeneous fuz or $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$10-local limits to classical Landau (Duong et al., 10 Apr 2025).

6. Uniqueness, Wasserstein stability, and terminology

A later development addresses uniqueness of weak solutions in a singular-potential regime. For $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$11, the fuzzy Landau equation is written on full phase space $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$12 as

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$13

with

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$14

A weak solution is required to belong to $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$15 with finite second moment and

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$16

If, in addition,

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$17

and $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$18 is bounded and Lipschitz, then the main theorem yields

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$19

for a continuous increasing modulus $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$20 with $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$21. The proof gives the explicit Osgood-type differential inequality

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$22

where

$\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$23

This yields conditional uniqueness in the class of weak solutions with finite second moment and the stated $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$24 integrability (Caja-Lopez, 29 Mar 2026).

Two proof strategies are developed. The first builds on the stochastic coupling method of Fournier and Guerin for the homogeneous Landau equation, reformulated analytically through a coupling $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$25 and a corresponding linear PDE or McKean–Vlasov-type SDE. The second is an alternative argument based on the symmetrization technique of Guillen and Silvestre, producing comparable stability estimates after a doubled-variable integration-by-parts and commutator analysis. The same framework also covers a wider class of nonlinear Fokker–Planck equations with singular kernels, including the $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$26D incompressible Euler equations, the Vlasov–Poisson system, and the Patlak–Keller–Segel model (Caja-Lopez, 29 Mar 2026).

A common terminological confusion is that “fuzzy” in the fuzzy Landau equation refers to spatial delocalisation of collisions by the kernel $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$27. This usage is distinct from “fuzzy” in “fuzzy dark matter,” where the term appears in the study of a quantum Landau equation for bosonic particles in gravitational interaction, developed in an infinite and spatially homogeneous system or a local approximation (Chavanis, 2020). The two contexts both involve Landau-type kinetic operators, conservation laws, and an $\partial_t f+v\cdot\nabla_x f=Q_{\sf fuz}(f,f),$28-theorem, but they arise from different modelling assumptions and represent different physical regimes.