We consider a Vlasov-Fokker-Planck equation governing the evolution of the density of interacting and diffusive matter in the space of positions and velocities. We use a probabilistic interpretation to obtain convergence towards equilibrium in Wasserstein distance with an explicit exponential rate. We also prove a propagation of chaos property for an associated particle system, and give rates on the approximation of the solution by the particle system. Finally, a transportation inequality for the distribution of the particle system leads to quantitative deviation bounds on the approximation of the equilibrium solution of the equation by an empirical mean of the particles at given time.Here (W t ) t 0 is a Brownian motion in the velocity space R d and f t is the law of (X t , V t ) in R 2d , so that ρ[f t ] is the law of X t in R d .Space homogeneous models of diffusive and interacting granular media (see [4]) have been studied by P. Cattiaux and the last two authors in particular [14], [19,20], by means of a stochastic interpretation analogous to (2) and a particle approximation analogous to (4) below. They were interpreted as gradient flows in the space of probability measures by J. A. Carrillo, R. J. McCann and C. Villani [11,12] (see also [26]), both approaches leading to explicit exponential (or algebraic for non uniformly convex potentials) rates of convergence to equilibrium. Also possibly time-uniform propagation of chaos was proven for the associated particle system.Obtaining rates of convergence to equilibrium for (1) is much more complex, as the equation simultaneously presents hamiltonian and gradient flows aspects. Much attention has recently been called to the linear noninteracting case of (1), when C = 0, also known as the kinetic Fokker-Planck equation. First of all a probabilistic approach based on Lyapunov functionals, and thus easy to check conditions, lead D. Talay [24], L. Wu [27] or D. Bakry, P. Cattiaux and the second author [2] to exponential or subexponential convergence to equilibrium in total variation distance. The case when A(v) = v and B(x) = ∇Ψ(x), and when the equilibrium solution is explicitely given by f ∞ (x, v) = e −Ψ(x)−|v| 2 /2 is studied in [16], [17] and [25, Chapter 7]: hypocoercivity analytic techniques are developed which, applied to this situation, give sufficient conditions, in terms of Poincaré or logarithmic Sobolev inequalities for the measure e −Ψ , to L 2 or entropic convergence with an explicit exponential rate. We also refer to [23] for the evolution of two species, modelled by two coupled Vlasov-Fokker-Planck equations.C. Villani's approach extends to the selfconsistent situation when C derives from a nonzero potential U (see [25, Chapter 17]): replacing the confinement force B(x) by a periodic boundary condition, and for small and smooth potential U , he obtains an explicit exponential rate of convergence of all solutions toward the unique normalized equilibrium solution e −|v| 2 /2 .In this work we consider the case when the equation is set on the whole R d , ...