In the Vlasov variant of Boltzmann equation, the Coulomb pair forces are included via the force term that appears as the prefactor of the derivative of the density with respect to velocity. In that case, one is effectively calculating a mean force or a "force field" by averaging over all the pair forces. The average is obtained by an integral over the pair potential weighted with the density.
Since this Vlasov approach allows to have good analytical insight into systems mainly governed by Coulomb forces (i.e Landau damping and other), is it possible to apply the same approach for Lennard Jones pair forces ?
I could not find references and I do not understand why this is the case.