pith. sign in

arxiv: 2410.10612 · v2 · pith:2542XXISnew · submitted 2024-10-14 · 🧮 math.AP · math-ph· math.MP

A Probabilistic Mean-Field Limit for the Vlasov-Poisson System for Ions

Megan Griffin-Pickering This is my paper

Pith reviewed 2026-05-23 18:44 UTC · model grok-4.3

classification 🧮 math.AP math-phmath.MP
keywords mean-field limitVlasov-Poisson systemionsprobabilistic methodslaw of large numbersplasma kinetic equationsCoulomb interactionmean-field derivation
0
0 comments X

The pith

The Vlasov-Poisson system for ions is the mean-field limit of N point ions interacting with thermalized electrons when the Coulomb force is truncated at scale N to the power minus beta with beta less than one third.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper proves that the ionic Vlasov-Poisson equation arises as the continuum limit of a microscopic system of N ions whose mutual interactions are cut off at spatial scale N to the minus beta. The allowed range is beta less than one third, which relaxes the earlier restriction of beta less than one fifteenth. The argument adapts a probabilistic mean-field technique to the nonlinear electron-density coupling that appears in the Poisson equation for ions. A quantitative uniform law of large numbers controls the fluctuations between the empirical particle density and its expected value for locally Lipschitz test functions. This yields convergence in probability of the microscopic trajectories to the solution of the kinetic equation.

Core claim

The Vlasov-Poisson system for ions is derived from the microscopic dynamics of N ions with interaction truncated at N^{-β} for β<1/3 by generalizing the probabilistic mean-field approach to handle the nonlinear Poisson coupling with the electron density.

What carries the argument

A quantitative uniform law of large numbers for convolutions between empirical measures of i.i.d. random variables and locally Lipschitz functions, applied to the nonlinear coupling.

If this is right

  • The mean-field limit holds for a wider range of truncation parameters than the previous derivation.
  • The probabilistic method extends directly to other kinetic equations whose force depends nonlinearly on the particle density.
  • Convergence in probability is obtained from the microscopic N-particle system to the ionic Vlasov-Poisson solution on a fixed time interval.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the uniform law of large numbers can be sharpened further, the same argument might remove the cutoff entirely in the limit N to infinity.
  • The technique may adapt to other plasma models that retain a nonlinear electron response, such as certain Vlasov-Maxwell systems.
  • Quantitative rates from the law of large numbers could be used to obtain explicit error bounds between particle simulations and the kinetic equation.

Load-bearing premise

The quantitative uniform law of large numbers for convolutions between empirical measures and locally Lipschitz functions holds and can be applied to the nonlinear ionic coupling.

What would settle it

An explicit counterexample in which the error term in the law of large numbers grows too fast to close the estimate when the truncation exponent beta reaches or exceeds one third.

Figures

Figures reproduced from arXiv: 2410.10612 by Megan Griffin-Pickering.

Figure 1
Figure 1. Figure 1: Straight line paths in T 2 between two points x and y contained in Br(0). (a) When r is large, the shortest path (solid line), which has length |x − y|T2 , may pass outside of Br(0). (b) When r is small, the shortest path is contained in Br(0). Lemma 3.6. Let x, y ∈ Br(0) for some r < 1 4 . Let y˜ ∈ [PITH_FULL_IMAGE:figures/full_fig_p013_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The union Br(x) ∪ Br(y). If an optimal z ∗ y lies in Br(x) ∩ Br(y) (shaded region), then h(y) ≤ h(x). Otherwise, we compare |∇g(z ∗ y )| with |∇g(u ∗ y )|, by using a mean value theorem argument along the solid line segment connecting u ∗ y and z ∗ y , whose length is no greater than |x−y|Td . Then, we estimate |z ∗ y −u ∗ y |Td . Using the definition of u ∗ y , |z ∗ y −u ∗ y |Td = (1−λ ∗ y )|z ∗ y −x|Td =… view at source ↗
read the original abstract

The Vlasov-Poisson system for ions is a kinetic equation for dilute, unmagnetised plasma. It describes the evolution of the ions in a plasma under the assumption that the electrons are thermalized. Consequently, the Poisson coupling for the electrostatic potential contains an additional exponential nonlinearity not present in the electron Vlasov-Poisson system. The system can be formally derived through a mean-field limit from a microscopic system of ions interacting with a thermalized electron distribution. However, it is an open problem to justify this limit rigorously for ions modelled as point charges. Existing results on the derivation of the three-dimensional ionic Vlasov-Poisson system, obtained by the author and Iacobelli [J. Math. Pures Appl. 135 (2020), pp. 199-255], require a truncation of the singularity in the Coulomb interaction at spatial scales of order $N^{-\beta}$ with $\beta<1/15$, which is more restrictive than the available results for the electron Vlasov-Poisson system. In this article, we prove that the Vlasov-Poisson system for ions can be derived from a microscopic system of ions and thermalized electrons with interaction truncated at scale $N^{-\beta}$ with $\beta<1/3$. We develop a generalisation of the probabilistic approach to mean-field limits developed in the works of Boers and Pickl [J. Stat. Phys. 164(1) (2016), pp. 1-16] and Lazarovici and Pickl [Arch. Ration. Mech. Anal. 225(3) (2017), pp. 1201-1231] that is applicable to interaction forces defined through a nonlinear coupling with the particle density. The proof is based on a quantitative uniform law of large numbers for convolutions between empirical measures of independent, identically distributed random variables and locally Lipschitz functions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper proves that the three-dimensional Vlasov-Poisson system for ions can be derived as a mean-field limit from a microscopic system of N ions interacting via a Coulomb potential truncated at scale N^{-β} with β < 1/3, where electrons are assumed thermalized. This improves upon the previous result requiring β < 1/15 by generalizing the probabilistic approach of Boers-Pickl and Lazarovici-Pickl to handle interaction forces defined through nonlinear Poisson coupling with the particle density, relying on a quantitative uniform law of large numbers for convolutions of empirical measures with locally Lipschitz functions.

Significance. If the result holds, this work substantially advances the rigorous justification of mean-field limits for the ionic Vlasov-Poisson system by relaxing the truncation requirement to a wider range of β, approaching the regime available for the standard electron Vlasov-Poisson system. The generalization of the probabilistic method to nonlinear couplings provides a useful framework that could apply to other kinetic models with density-dependent forces. The use of quantitative LLN is a strength in making the estimates explicit.

major comments (2)
  1. [Main proof (generalization of Boers-Pickl/Lazarovici-Pickl method)] The application of the quantitative uniform law of large numbers to the force field F[μ] obtained from the nonlinear elliptic problem (ρ_ions - exp(φ)) requires uniform control on the local Lipschitz modulus of F with respect to the empirical measure μ^N. It is not evident from the argument whether fluctuations in μ^N at scale N^{-β} are controlled sufficiently to prevent the Lipschitz constant from growing with N, which could invalidate the Gronwall estimate for β close to 1/3.
  2. [Theorem 1.1 and its proof] The claim that the method extends the previous β<1/15 to β<1/3 relies on the LLN holding for the nonlinear map; however, the dependence of the Poisson solution on the global measure may introduce additional error terms not fully accounted for in the uniform bound.
minor comments (2)
  1. [Introduction] The comparison with the electron Vlasov-Poisson results could include a brief remark on why the ionic case is more challenging due to the exponential nonlinearity.
  2. [Notation and setup] Ensure that the definition of the truncated interaction potential is clearly stated with the exact form of the cutoff.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation of the significance of our work and for the detailed comments on the proof. We address the major comments point by point below, providing clarifications based on the estimates in the manuscript.

read point-by-point responses

  1. Referee: [Main proof (generalization of Boers-Pickl/Lazarovici-Pickl method)] The application of the quantitative uniform law of large numbers to the force field F[μ] obtained from the nonlinear elliptic problem (ρ_ions - exp(φ)) requires uniform control on the local Lipschitz modulus of F with respect to the empirical measure μ^N. It is not evident from the argument whether fluctuations in μ^N at scale N^{-β} are controlled sufficiently to prevent the Lipschitz constant from growing with N, which could invalidate the Gronwall estimate for β close to 1/3.

    Authors: We appreciate the referee raising this point on the control of the Lipschitz modulus. The manuscript establishes uniform control in Section 3 via the quantitative LLN (Proposition 3.4) combined with elliptic regularity for the nonlinear Poisson problem. The truncation at scale N^{-β} (β<1/3) ensures density fluctuations remain small enough in the relevant Sobolev norms that the local Lipschitz constant of F[μ] stays bounded independently of N (see the estimates after (3.15) and the application in (4.8)). This bound is independent of the empirical fluctuations at the truncation scale and closes the Gronwall argument. We are prepared to add an explicit remark or auxiliary lemma in a revision to make this uniformity more transparent. revision: partial

  2. Referee: [Theorem 1.1 and its proof] The claim that the method extends the previous β<1/15 to β<1/3 relies on the LLN holding for the nonlinear map; however, the dependence of the Poisson solution on the global measure may introduce additional error terms not fully accounted for in the uniform bound.

    Authors: The dependence of the Poisson solution on the global empirical measure is incorporated directly into the formulation of the quantitative uniform LLN for the nonlinear map F[μ]. In the proof of Theorem 1.1 (Section 4), the total error is split into the LLN approximation error for the force and the stability error of the elliptic map; both are controlled uniformly by the same truncation scale and the a priori L^∞ bounds on the potential (Lemma 2.3 and (4.15)–(4.20)). No additional unaccounted terms arise beyond those already bounded in the argument. If the referee prefers, we can expand the error decomposition explicitly in the revised manuscript. revision: partial

Circularity Check

0 steps flagged

No circularity: proof derives mean-field limit via independent generalization of external probabilistic methods

full rationale

The paper's derivation chain consists of a rigorous mathematical proof that generalizes the probabilistic mean-field approach of Boers-Pickl and Lazarovici-Pickl (distinct authors) to handle the nonlinear Poisson coupling for ions. The quantitative uniform LLN for convolutions with locally Lipschitz functions is invoked as the core tool and applied directly to the ionic force map; no step reduces by construction to a fitted parameter, self-defined quantity, or load-bearing self-citation. The author's prior joint work with Iacobelli is cited only to contextualize the improved exponent β<1/3 relative to the earlier β<1/15 result and does not supply the new argument. The central claim therefore remains self-contained against external benchmarks and does not exhibit any of the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

No free parameters or invented entities; the work relies on standard assumptions from probability and analysis in kinetic theory.

axioms (2)
  • standard math Empirical measures of i.i.d. particles satisfy a quantitative uniform law of large numbers when convolved with locally Lipschitz functions.
    This is the key tool invoked for the mean-field limit in the abstract.
  • domain assumption The truncation of the Coulomb interaction at scale N^{-β} preserves the necessary properties for the limit to hold.
    Assumed in the setup of the microscopic system.

pith-pipeline@v0.9.0 · 5891 in / 1220 out tokens · 56862 ms · 2026-05-23T18:44:18.799988+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Mean-field derivation of a two-dimensional signal-dependent parabolic-elliptic Keller-Segel system in algebraic scaling

    math.PR 2026-05 unverdicted novelty 6.0

    Proves weak propagation of chaos and short-time L1 density convergence for a particle system with algebraic scaling where moderate interactions occur in the diffusive term, yielding a signal-dependent 2D Keller-Segel system.

Reference graph

Works this paper leans on

51 extracted references · 51 canonical work pages · cited by 1 Pith paper

  1. [1]

    Ambrosio

    L. Ambrosio. Transport equation and Cauchy problem for Non-smooth Vector Fields. In B. Dacorogna and P. Marcellini, editors, Calculus of Variations and Nonlinear Partial Differential Equations: Lectures given at the C.I.M.E. Summer School held in Cetraro, Italy, June 27 - July 2, 2005. Springer-Verlag Berlin Heidelberg, 2008

    work page 2005

  2. [2]

    Ambrosio, N

    L. Ambrosio, N. Gigli, and G. Savaré. Gradient Flows In Metric Spaces and in the Space of Probability Measures. Lectures in Mathematics. ETH Zürich. Birkhäuser Basel, 2008

    work page 2008

  3. [3]

    Bardos, F

    C. Bardos, F. Golse, T. T. Nguyen, and R. Sentis. The Maxwell-Boltzmann approximation for ion kinetic modeling. Phys. D, 376/377:94–107, 2018

    work page 2018

  4. [4]

    P. M. Bellan. Fundamentals of Plasma Physics. Cambridge University Press, 2008

    work page 2008

  5. [5]

    Bernstein

    S. Bernstein. The Theory of Probabilities. Gastehizdat Publishing House, Moscow., 1946

    work page 1946

  6. [6]

    Boers and P

    N. Boers and P. Pickl. On mean field limits for dynamical systems. J. Stat. Phys., 164(1):1–16, 2016

    work page 2016

  7. [7]

    Bonhomme, T

    G. Bonhomme, T. Pierre, G. Leclert, and J. Trulsen. Ion phase space vortices in ion beam-plasma systems and their relation with the ion acoustic instability: numerical and experimental results. Plasma Physics and Controlled Fusion, 33(5):507–520, 1991

    work page 1991

  8. [8]

    Boucheron, G

    S. Boucheron, G. Lugosi, and P. Massart. Concentration Inequalities: A Nonasymptotic Theory of Inde- pendence. Oxford University Press, 2013

    work page 2013

  9. [9]

    F. Bouchut. Global weak solution of the Vlasov-Poisson system for small electrons mass. Comm. Partial Differential Equations, 16(8-9):1337–1365, 1991

    work page 1991

  10. [10]

    Bouchut and J

    F. Bouchut and J. Dolbeault. On long time asymptotics of the Vlasov-Fokker-Planck equation and of the Vlasov-Poisson-Fokker-Planck system with Coulombic and Newtonian potentials. Differential Integral Equations, 8(3):487–514, 1995. 32

    work page 1995

  11. [11]

    Braun and K

    W. Braun and K. Hepp. The Vlasov dynamics and its fluctuations in the 1 /N limit of interacting classical particles. Comm. Math. Phys., 56(2):101–113, 1977

    work page 1977

  12. [12]

    Bresch, P.-E

    D. Bresch, P.-E. Jabin, and J. Soler. A new approach to the mean-field limit of Vlasov-Fokker-Planck equations, 2023. Preprint, arXiv:2203.15747

    work page arXiv 2023

  13. [13]

    Cesbron and M

    L. Cesbron and M. Iacobelli. Global well-posedness of Vlasov-Poisson-type systems in bounded domains. Anal. PDE, 16(10):2465–2494, 2023

    work page 2023

  14. [14]

    Chaintron and A

    L.-P. Chaintron and A. Diez. Propagation of chaos: a review of models, methods and applications. I. Models and methods. Kinet. Relat. Models, 15(6):895–1015, 2022

    work page 2022

  15. [15]

    Chaintron and A

    L.-P. Chaintron and A. Diez. Propagation of chaos: a review of models, methods and applications. II. Applications. Kinet. Relat. Models, 15(6):1017–1173, 2022

    work page 2022

  16. [16]

    R. J. DiPerna and P. L. Lions. Ordinary differential equations, transport theory and Sobolev spaces. Inventiones mathematicae, 98(3):511–547, Oct 1989

    work page 1989

  17. [17]

    R. L. Dobrushin. Vlasov Equations. Funktsional. Anal. i Prilozhen., 13(2):48–58, 1979

    work page 1979

  18. [18]

    P. Flynn. Local well-posedness of the Vlasov-Poisson-Landau system and the massless electron system,

    work page

  19. [19]

    Preprint, arXiv:2310.00777

    work page arXiv

  20. [20]

    Flynn and Y

    P. Flynn and Y . Guo. The Massless Electron Limit of the Vlasov–Poisson–Landau System. Communica- tions in Mathematical Physics, 405(2):27, 2024

    work page 2024

  21. [21]

    Fournier and A

    N. Fournier and A. Guillin. On the rate of convergence in Wasserstein distance of the empirical measure. Probability Theory and Related Fields, 162(3):707–738, 2015

    work page 2015

  22. [22]

    F. Golse. On the dynamics of large particle systems in the mean field limit. In A. Muntean, J. Rademacher, and A. Zagaris, editors, Macroscopic and Large Scale Phenomena: Coarse Graining, Mean Field Limits and Ergodicity, volume 3 of Lect. Notes Appl. Math. Mech., pages 1–144. Springer, 2016

    work page 2016

  23. [23]

    P. Graß. Microscopic derivation of Vlasov equations with singular potentials . PhD thesis, Ludwig- Maximilians-Universität München, 2019

    work page 2019

  24. [24]

    Griffin-Pickering and M

    M. Griffin-Pickering and M. Iacobelli. Singular limits for plasmas with thermalised electrons. J. Math. Pures Appl. (9), 135:199–255, 2020

    work page 2020

  25. [25]

    Griffin-Pickering and M

    M. Griffin-Pickering and M. Iacobelli. Global strong solutions in R3 for ionic Vlasov-Poisson systems. Kinet. Relat. Models, 14(4):571–597, 2021

    work page 2021

  26. [26]

    Griffin-Pickering and M

    M. Griffin-Pickering and M. Iacobelli. Global well-posedness for the Vlasov-Poisson system with mass- less electrons in the 3-dimensional torus.Comm. Partial Differential Equations, 46(10):1892–1939, 2021

    work page 1939

  27. [27]

    Griffin-Pickering and M

    M. Griffin-Pickering and M. Iacobelli. Recent developments on the well-posedness theory for Vlasov- type equations. In From particle systems to partial differential equations, volume 352 of Springer Proc. Math. Stat., pages 301–319. Springer, Cham, [2021] ©2021

    work page 2021

  28. [28]

    Griffin-Pickering and M

    M. Griffin-Pickering and M. Iacobelli. Stability in quasineutral plasmas with thermalized electrons, 2023. Preprint, arXiv:2307.07561

    work page arXiv 2023

  29. [29]

    A. V . Gurevich and L. P. Pitaevsky. Non-linear dynamics of a rarefied ionized gas.Progress in Aerospace Sciences, 16(3):227 – 272, 1975

    work page 1975

  30. [30]

    Han-Kwan and M

    D. Han-Kwan and M. Iacobelli. The quasineutral limit of the Vlasov-Poisson equation in Wasserstein metric. Commun. Math. Sci., 15(2):481–509, 2017

    work page 2017

  31. [31]

    M. Hauray. Mean field limit for the one dimensional Vlasov-Poisson equation. In Séminaire Laurent Schwartz—Équations aux dérivées partielles et applications. Année 2012–2013 , Exp. No. XXI, Sémin. Équ. Dériv. Partielles. École Polytech., Palaiseau, 2014. 33

    work page 2012

  32. [32]

    Hauray and P.-E

    M. Hauray and P.-E. Jabin. N-particles approximation of the Vlasov equations with singular potential. Arch. Ration. Mech. Anal., 183(3):489–524, 2007

    work page 2007

  33. [33]

    Hauray and P.-E

    M. Hauray and P.-E. Jabin. Particle approximation of Vlasov equations with singular forces: propagation of chaos. Ann. Sci. Éc. Norm. Supér. (4), 48(4):891–940, 2015

    work page 2015

  34. [34]

    M. Herda. On massless electron limit for a multispecies kinetic system with external magnetic field. Journal of Differential Equations, 260(11):7861 – 7891, 2016

    work page 2016

  35. [35]

    Iacobelli

    M. Iacobelli. A New Perspective on Wasserstein Distances for Kinetic Problems. Archive for Rational Mechanics and Analysis, 244(1):27–50, 2022

    work page 2022

  36. [36]

    P.-E. Jabin. A review of the mean field limits for Vlasov equations. Kinet. Relat. Models, 7(4):661–711, 2014

    work page 2014

  37. [37]

    Jabin and Z

    P.-E. Jabin and Z. Wang. Mean field limit and propagation of chaos for Vlasov systems with bounded forces. J. Funct. Anal., 271(12):3588–3627, 2016

    work page 2016

  38. [38]

    Jabin and Z

    P.-E. Jabin and Z. Wang. Mean field limit for stochastic particle systems. In N. Bellomo, P. Degond, and E. Tadmor, editors, Active Particles, Volume 1 : Advances in Theory, Models, and Applications , pages 379–402. Springer International Publishing, Cham, 2017

    work page 2017

  39. [39]

    J. H. Jeans. On the theory of star-streaming and the structure of the universe. Monthly Notices of the Royal Astronomical Society, 76:70–84, 1915

    work page 1915

  40. [40]

    M. K. H. Kiessling. The Microscopic Foundations of Vlasov Theory for Jellium-Like Newtonian N-Body Systems. Journal of Statistical Physics, 155(6):1299–1328, 2014

    work page 2014

  41. [41]

    Lazarovici

    D. Lazarovici. The Vlasov-Poisson dynamics as the mean field limit of extended charges. Comm. Math. Phys., 347(1):271–289, 2016

    work page 2016

  42. [42]

    Lazarovici and P

    D. Lazarovici and P. Pickl. A mean field limit for the Vlasov-Poisson system. Arch. Ration. Mech. Anal., 225(3):1201–1231, 2017

    work page 2017

  43. [43]

    Ledoux and M

    M. Ledoux and M. Talagrand. Probability in Banach Spaces: Isoperimetry and Processes . Springer Berlin, Heidelberg, 1991

    work page 1991

  44. [44]

    R. J. Mason. Computer simulation of ion-acoustic shocks. The diaphragm problem. The Physics of Fluids, 14(9):1943–1958, 1971

    work page 1943

  45. [45]

    Y . V . Medvedev. Ion front in an expanding collisionless plasma.Plasma Physics and Controlled Fusion, 53(12):125007, nov 2011

    work page 2011

  46. [46]

    Neunzert and J

    H. Neunzert and J. Wick. Die Approximation der Lösung von Integro-Differentialgleichungen durch endliche Punktmengen. In Numerische Behandlung nichtlinearer Integrodifferential-und Differentialgle- ichungen., volume 395 of Lecture Notes in Math., pages 275–290, Berlin, Heidelberg, 1974. Springer

    work page 1974

  47. [47]

    Sakanaka, C

    P. Sakanaka, C. Chu, and T. Marshall. Formation of ion-acoustic collisionless shocks. The Physics of Fluids, 14(611), 1971

    work page 1971

  48. [48]

    S. Serfaty. Mean field limit for Coulomb-type flows. Duke Math. J., 169(15):2887–2935, 2020. With an appendix by M. Duerinckx and S. Serfaty

    work page 2020

  49. [49]

    Sznitman

    A.-S. Sznitman. Topics in propagation of chaos. InÉcole d’Été de Probabilités de Saint-Flour XIX—1989, volume 1464 of Lecture Notes in Math., pages 165–251. Springer, Berlin, 1991

    work page 1989

  50. [50]

    Titchmarsh

    E. Titchmarsh. Eigenfunction Expansions Associated with Second-Order Differential Equations, Part II. Oxford University Press, 1958

    work page 1958

  51. [51]

    A. A. Vlasov. On the vibration properties of the electron gas. Zh. Eksper. Teor. Fiz., 8(3):291, 1938. 34

    work page 1938