Consider the Vlasov-Poisson-Landau system with Coulomb potential in the

weakly collisional regime on a $3$-torus, i.e. $$\begin{aligned} \partial_t

F(t,x,v) + v_i \partial_{x_i} F(t,x,v) + E_i(t,x) \partial_{v_i} F(t,x,v) = \nu

Q(F,F)(t,x,v),\\ E(t,x) = \nabla \Delta^{-1} (\int_{\mathbb R^3} F(t,x,v)\,

\mathrm{d} v – \frac{1}{(2\pi)^3}\int_{\mathbb T^3} \int_{\mathbb R^3}

F(t,x,v)\, \mathrm{d} v \, \mathrm{d} x), \end{aligned}$$ with $\nu\ll 1$. We

prove that for $\epsilon>0$ sufficiently small (but independent of $\nu$),

initial data which are $O(\epsilon \nu^{1/3})$-Sobolev space perturbations from

the global Maxwellians lead to global-in-time solutions which converge to the

global Maxwellians as $t\to \infty$. The solutions exhibit uniform-in-$\nu$

Landau damping and enhanced dissipation.

Our main result is analogous to an earlier result of Bedrossian for the

Vlasov-Poisson-Fokker-Planck equation with the same threshold. However, unlike

in the Fokker-Planck case, the linear operator cannot be inverted explicitly

due to the complexity of the Landau collision operator. For this reason, we

develop an energy-based framework, which combines Guo’s weighted energy method

with the hypocoercive energy method and the commuting vector field method. The

proof also relies on pointwise resolvent estimates for the linearized density

equation.