Abstract In this paper we develop a new martingale method to show the convergence of the regularized empirical measure of many particle systems in probability under a Sobolev norm to the corresponding mean field PDE. Our method works well for the simple case of Fokker Planck equation and we can estimate a lower bound of the rate of convergence. This method can be generalized to more complicated systems with interactions.