🤖 AI Summary
To address the challenge of numerically solving the homogeneous Landau equation for plasmas—specifically, the nonlinear dependence of the velocity field on the distribution function—this paper proposes a structure-preserving Lagrangian particle method. The core method introduces fractional matching, a novel technique that explicitly models this nonlinearity via a learnable, dynamic fractional function, thereby eliminating the need for kernel density estimation. We theoretically prove that the fractional matching loss is equivalent to KL-divergence minimization, ensuring exact conservation of mass and momentum. Coupled with normalizing flow-based transport mapping, the approach achieves scalability to high dimensions. In benchmark physical settings—including Coulomb interactions—the method attains high accuracy (reducing error by one to two orders of magnitude versus conventional particle methods) and computational efficiency (nearly linear complexity), substantially enhancing the feasibility and reliability of high-dimensional phase-space simulations.
📝 Abstract
We propose a novel score-based particle method for solving the Landau equation in plasmas, that seamlessly integrates learning with structure-preserving particle methods [arXiv:1910.03080]. Building upon the Lagrangian viewpoint of the Landau equation, a central challenge stems from the nonlinear dependence of the velocity field on the density. Our primary innovation lies in recognizing that this nonlinearity is in the form of the score function, which can be approximated dynamically via techniques from score-matching. The resulting method inherits the conservation properties of the deterministic particle method while sidestepping the necessity for kernel density estimation in [arXiv:1910.03080]. This streamlines computation and enhances scalability with dimensionality. Furthermore, we provide a theoretical estimate by demonstrating that the KL divergence between our approximation and the true solution can be effectively controlled by the score-matching loss. Additionally, by adopting the flow map viewpoint, we derive an update formula for exact density computation. Extensive examples have been provided to show the efficiency of the method, including a physically relevant case of Coulomb interaction.