In this talk, I will discuss fine modified scattering properties of small data solutions to the Vlasov-Poisson system. First, I will show that the spatial density and the force field satisfy polyhomogeneous expansions of any order. As a result, we will obtain an enhanced modified scattering result for this non-linear system. I will show that the distribution function converges, with an arbitrary rate, to a regular distribution function along high order modifications to the characteristics of the linearised problem. This is joint work with Léo Bigorgne (Université de Rennes).