Singular mean-field limits via a multiscale mollification metric

We consider a general class of first order ODE systems for the evolution of $N$ interacting particles (in Euclidean space $\mathbb{R}^d$) in a mean-field regime. The class of interactions treated includes singular interactions of inverse power type up to power $d+1$, attractive or repulsive, and not necessarily deriving from a potential -- unlike, for instance, the modulated energy method. We introduce a new method to prove quantitative convergence of the discrete system to solutions of the mean-field equation. It relies on studying the evolution of a metric encoding a multiscale control of the difference between the empirical measure and its limit, via mollification by heat kernels. We prove that the desired convergence holds (i) up to the maximal time of existence of the smooth solution to the limiting equation if the singularity is sub-coulombic in any dimension, or coulombic in dimensions 1 and 2 (where, to do so, we introduce a notion of weak solution to the ODE system), or (ii) for short time in the case of Coulomb singularity in dimension 3 and above and (iii) up to a short $N$-dependent timescale for super-coulombic interactions in all dimensions. The latter two results are demonstrated to be optimal as we prove that collisions occur within the same timescale for a class of attractive interactions.