  We consider a particle system in 1D, interacting via repulsive or attractive
Coulomb forces. We prove the trajectorial propagation of molecular chaos
towards a nonlinear SDE associated to the Vlasov-Poisson-Fokker-Planck
equation. We obtain a quantitative estimate of convergence in expectation, with
an optimal convergence rate of order $N^{-1/2}$. We also prove some exponential
concentration inequalities of the associated empirical measures. A key argument
is a weak-strong stability estimate on the (nonlinear) VPFP equation, that we
are able to adapt for the particle system in some sense.
