A Neural Score-Based Particle Method for the Vlasov-Maxwell-Landau System
arXiv:2603.25832v1 Announce Type: cross Abstract: Plasma modeling is central to the design of nuclear fusion reactors, yet simulating collisional plasma kinetics from first principles remains a formidable computational challenge: the Vlasov-Maxwell-Landau (VML) system describes six-dimensional phase-space transport under self-consistent electromagnetic fields together with the nonlinear, nonlocal Landau collision operator. A recent deterministic particle method for the full VML system estimates the velocity score function via the blob method, a kernel-based approximation with $O(n^2)$ cost. In — Vasily Ilin, Jingwei Hu
View PDF HTML (experimental)
Abstract:Plasma modeling is central to the design of nuclear fusion reactors, yet simulating collisional plasma kinetics from first principles remains a formidable computational challenge: the Vlasov-Maxwell-Landau (VML) system describes six-dimensional phase-space transport under self-consistent electromagnetic fields together with the nonlinear, nonlocal Landau collision operator. A recent deterministic particle method for the full VML system estimates the velocity score function via the blob method, a kernel-based approximation with $O(n^2)$ cost. In this work, we replace the blob score estimator with score-based transport modeling (SBTM), in which a neural network is trained on-the-fly via implicit score matching at $O(n)$ cost. We prove that the approximated collision operator preserves momentum and kinetic energy, and dissipates an estimated entropy. We also characterize the unique global steady state of the VML system and its electrostatic reduction, providing the ground truth for numerical validation. On three canonical benchmarks -- Landau damping, two-stream instability, and Weibel instability -- SBTM is more accurate than the blob method, achieves correct long-time relaxation to Maxwellian equilibrium where the blob method fails, and delivers $50%$ faster runtime with $4\times$ lower peak memory.
Comments: presented at ICLR AI&PDE workshop
Subjects:
Numerical Analysis (math.NA); Machine Learning (cs.LG); Analysis of PDEs (math.AP)
MSC classes: 35Q83, 65M75, 82D10
ACM classes: G.1.8; I.2.6
Cite as: arXiv:2603.25832 [math.NA]
(or arXiv:2603.25832v1 [math.NA] for this version)
https://doi.org/10.48550/arXiv.2603.25832
arXiv-issued DOI via DataCite
Submission history
From: Vasily Ilin [view email] [v1] Thu, 26 Mar 2026 18:55:44 UTC (24,149 KB)
Sign in to highlight and annotate this article

Conversation starters
Daily AI Digest
Get the top 5 AI stories delivered to your inbox every morning.
More about
researchpaperarxiv
Google launches Gemma 4: four open-weight models from smartphones to workstations
Built from the same research as Gemini 3, the new family spans a 2B edge model that runs on a Raspberry Pi to a 31B dense model currently ranked third on the Arena AI open-model leaderboard. The Apache 2.0 licence is a significant shift from previous Gemma releases. Google has released Gemma 4, the latest [ ] This story continues at The Next Web
Knowledge Map
Connected Articles — Knowledge Graph
This article is connected to other articles through shared AI topics and tags.






Discussion
Sign in to join the discussion
No comments yet — be the first to share your thoughts!