Stochastic numerical approximation for nonlinear Fokker-Planck equations with singular kernels

Abstract: This paper studies the convergence rate of the Euler-Maruyama scheme for systems of interacting particles used to approximate solutions of nonlinear Fokker-Planck equations with singular interaction kernels, such as the Keller-Segel model. We derive explicit error estimates in the large-particle limit for two objects: the empirical measure of the interacting particle system and the density distribution of a single particle. Specifically, under certain assumptions on the interaction kernel and initial conditions, we show that the convergence rate of both objects towards solutions of the corresponding nonlinear FokkerPlanck equation depends polynomially on N (the number of particles) and on h (the discretization step). The analysis shows that the scheme converges despite singularities in the drift term. To the best of our knowledge, there are no existing results in the literature of such kind for the singular kernels considered in this work.