1 From Newton to the Vlasov equation
The derivation opens with \(N\) identical particles in \(\mathbb {R}^d\) under the mean-field force \(-\frac{1}{N}\sum _j \nabla W\). The empirical measure of the particle system is not merely close to a solution of the limit equation: under the standing evenness assumption it is one, exactly, for every \(N\). The chapter records the standing assumption, the particle system, the computation that exposes the \(O(1/N)\) diagonal remainder, the cancellation that kills it, and the two solution concepts — weak and Lagrangian — whose relation the rest of the development settles.
Let \(W : \mathbb {R}^d \to \mathbb {R}\) be an even function of class \(C^{1,1}\): differentiable, with a globally Lipschitz gradient of constant \(L := \mathrm{Lip}(\nabla W) {\lt} \infty \). This assumption is in force throughout. Evenness gives \(\nabla W(0) = 0\) — the cancellation that makes the empirical measure an exact weak solution below.
The mean-field Hamiltonian of \(N\) identical unit-mass particles,
The \(1/N\) scaling keeps the kinetic and potential energies of the same order as \(N \to \infty \).
Newton’s equations for \(N\) particles under the mean-field force:
A solution is a curve \((X,V) : \mathbb {R} \to (\mathbb {R}^d \times \mathbb {R}^d)^N\) satisfying both systems at every time.
The empirical measure of a configuration \((X,V) \in (\mathbb {R}^d \times \mathbb {R}^d)^N\):
a probability measure on phase space.
The curve of empirical measures \(t \mapsto \mu ^N_t\) carried by a solution of Newton’s equations.
Along any Newton solution, for every test function \(\varphi \in C_c^\infty (\mathbb {R}^d \times \mathbb {R}^d)\) the pairing \(t \mapsto \langle \mu ^N_t, \varphi \rangle \) is differentiable and
where the remainder is the diagonal correction \(R_N = \frac{1}{N^2} \sum _i \nabla W(0) \cdot \nabla _v \varphi (x_i, v_i)\), of size \(|R_N| \le \frac{1}{N} \| \nabla W\| _\infty \| \nabla _v \varphi \| _\infty \).
Differentiate the pairing along the trajectories and apply the chain rule with Newton’s equations. Completing the interaction sum \(\sum _{j \ne i}\) to \(\sum _j\) produces exactly the diagonal term \(\nabla W(0)\), and the stated bound follows from the sup-norms of \(\nabla W\) and \(\nabla _v \varphi \).
Under Assumption 1 the remainder vanishes identically: \(W\) even forces \(\nabla W(0) = 0\), so the empirical measure of every Newton solution satisfies the distributional Vlasov equation exactly — not merely up to an \(O(1/N)\) error. This is the sense in which the particle system already solves the limit equation.
Evenness of \(W\) gives \(\nabla W(0) = 0\), which kills the diagonal correction in the weak-evolution identity.
The convolution of a vector field \(k : \mathbb {R}^d \to \mathbb {R}^d\) with a finite measure \(\rho \): \((k * \rho )(x) := \int k(x - y)\, \mathrm{d}\rho (y)\). The Vlasov force is \(\nabla W * \rho _t\), with \(\rho _t\) the spatial marginal of \(f_t\).
A weak solution of the Vlasov equation: a curve \(t \mapsto f_t\) of measures such that for every \(\varphi \in C_c^\infty (\mathbb {R}^d \times \mathbb {R}^d)\) the map \(t \mapsto \langle f_t, \varphi \rangle \) is differentiable at every \(t \in \mathbb {R}\) with
The predicate states exactly this identity; membership in \(\mathcal{P}_1\), windows, and continuity in time enter as explicit hypotheses of the theorems that consume it (cf. the definition in Section 1.3 of the paper, whose prose imposes them up front).
A characteristic flow for the field driven by \(\rho \): a map \(\Phi = (X,V)\) solving the mean-field ODE
the system along which the Vlasov equation transports mass.
A Lagrangian solution is a weak solution \(f\) for which there exists a characteristic flow \(\Phi \) driven by its own marginal \(\rho ^f\) with \(f_t = (\Phi _t)_\# f_0\) for every \(t\) — the solution is transported by the flow it generates. Every Lagrangian solution is weak by definition; the converse is the superposition principle (Theorem 1.7 of the paper, the final chapter).