A fast spectral particle method for the Landau equation

 

Giacomo Borghi, Lorenzo Pareschi (preprint arXiv:2609.23086)

We propose a fast deterministic particle method for the spatially homogeneous Landau equation. The method combines a particle representation of the solution with a Fourier approximation of the nonlinear collision flux. Nonuniform fast Fourier transforms are used to reconstruct the density from the particles and to evaluate the flux and density at the particle locations, while the convolutional structure of the flux enables its efficient computation by FFTs. For a fixed transform tolerance, the cost per time step is O(N+M^d logM), where N is the number of particles and M the number of Fourier modes per velocity dimension. We establish a consistency estimate for the reconstructed velocity field on regions where the reference density is bounded away from zero. The estimate separates the spectral truncation error from the particle-density reconstruction error, showing how the latter can dominate for sufficiently smooth densities.

Control variates with neural surrogates for uncertainty quantification in kinetic equations

Wei Chen, Giacomo Dimarco, Lorenzo Pareschi (preprint arXiv:2608.01360)

Efficient uncertainty quantification for kinetic equations with random inputs is challenging because it requires repeated simulations of high-dimensional models, such as the Boltzmann, Landau, and related collisional equations, whose computational cost can quickly become prohibitive. Multifidelity control variates address this difficulty by coupling a small number of high-fidelity simulations with many evaluations of lower-complexity reduced models. In this work, we analyze the case in which the reduced model is replaced by a neural surrogate rather than evaluated through a classical numerical scheme.

High-Order Asymptotic-Preserving Schemes for Kinetic Equations from Rarefied to Incompressible Regimes

Giacomo Dimarco, Axel Klar, Theresa Köfler, Lorenzo Pareschi, Tiwari Sudarshan (preprint arXiv:2606.28040)

This work introduces a novel high-order numerical framework for solving kinetic equations, designed to remain uniformly valid across all regimes of the mean free path, spanning from the rarefied kinetic scale to the incompressible hydrodynamic limit. The method is built upon a micro-macro decomposition, which reformulates the underlying kinetic equation into a coupled system consisting of a macroscopic part, representing the fluid-dynamic evolution, and a microscopic part, describing the non-equilibrium deviations. The proposed framework ensures high-order temporal accuracy through the use of Implicit-Explicit Runge-Kutta methods, which provide stability and efficiency in stiff regimes, while spatial resolution is enhanced by combining finite-difference WENO reconstructions with high-order central difference approximations.

How opinions shape epidemics: a graphon-based kinetic approach

Abu Safyan Ali, Elisa Calzola, Giacomo Dimarco, Lorenzo Pareschi, Thomas Rey (preprint arXiv:2605.14139)

Understanding the mutual influence between social behavior and physical health is crucial for designing effective epidemic mitigation strategies. Individual interactions drive the evolution of opinions, which in turn shape how infectious diseases are perceived and consequently how they spread within a population, for instance through the adoption or rejection of preventive measures. At the same time, the distribution and dynamics of physical contacts play a fundamental role in determining transmission patterns.

Two-Time-Scale Learning Dynamics: A Population View of Neural Network Training

Giacomo Borghi, Hyesung Im, Lorenzo Pareschi (preprint arXiv:2603.19808)

Population-based learning paradigms, including evolutionary strategies, Population-Based Training (PBT), and recent model-merging methods, combine fast within-model optimisation with slower population-level adaptation. Despite their empirical success, a general mathematical description of the resulting collective training dynamics remains incomplete. We introduce a theoretical framework for neural network training based on two-time-scale population dynamics. We model a population of neural networks as an interacting agent system in which network parameters evolve through fast noisy gradient updates of SGD/Langevin type, while hyperparameters evolve through slower selection--mutation dynamics.

Micro-Macro Tensor Neural Surrogates for Uncertainty Quantification in Collisional Plasma

Wei Chen, Giacomo Dimarco, Lorenzo Pareschi (preprint arXiv:2512.24205)

Plasma kinetic equations exhibit pronounced sensitivity to microscopic perturbations in model parameters and data, making reliable and efficient uncertainty quantification (UQ) essential for predictive simulations. However, the cost of uncertainty sampling, the high-dimensional phase space, and multiscale stiffness pose severe challenges to both computational efficiency and error control in traditional numerical methods. These aspects are further emphasized in presence of collisions where the high-dimensional nonlocal collision integrations and conservation properties pose severe constraints. To overcome this, we present a variance-reduced Monte Carlo framework for UQ in the Vlasov--Poisson--Landau (VPL) system, in which neural network surrogates replace the multiple costly evaluations of the Landau collision term. The method couples a high-fidelity, asymptotic-preserving VPL solver with inexpensive, strongly correlated surrogates based on the Vlasov--Poisson--Fokker--Planck (VPFP) and Euler--Poisson (EP) equations.

Asymptotic preserving methods for the low Mach limit in discrete velocity models approximating kinetic equations

Giacomo Dimarco, Axel Klar, Theresa Köfler, Lorenzo Pareschi (preprint arXiv:2512.19847)

We consider a Lattice Boltzmann type discrete velocity model in the low Mach number scaling and develop a corresponding numerical scheme that remains uniformly valid across all regimes of the mean free path, from the kinetic to the hydrodynamic scale. The proposed framework ensures high order temporal accuracy through the use of Implicit Explicit Runge Kutta methods, which provide stability and efficiency in stiff regimes, while spatial resolution is enhanced by combining finite difference WENO reconstructions with high order central difference approximations.

Collective Annealing by Switching Temperatures: a Boltzmann-type description

Frédéric Blondeel, Lorenzo Pareschi, Giovanni Samaey (preprint arXiv:2512.13522)

The design of effective cooling strategies is a crucial component in simulated annealing algorithms based on the Metropolis method. Traditionally, this is achieved through inverse logarithmic decays of the temperature to ensure convergence to global minima. In this work, we propose Collective Annealing by Switching Temperatures (CAST), a novel collective simulated annealing dynamic in which agents interact to learn an adaptive cooling schedule.

High-Order Asymptotic-Preserving IMEX schemes for an ES-BGK model for Gas Mixtures

Domenico Caparello, Lorenzo Pareschi, Thomas Rey (preprint arXiv:2511.22304)

In this work we construct a high-order Asymptotic-Preserving (AP) Implicit-Explicit (IMEX) scheme for the ES-BGK model for gas mixtures introduced in [Brull, Commun. Math. Sci., 2015]. The time discretization is based on the IMEX strategy proposed in [Filbet, Jin, J. Sci. Comput., 2011] for the single-species BGK model and is here extended to the multi-species ES-BGK setting. The resulting method is fully explicit, uniformly stable with respect to the Knudsen number and, in the fluid regime, it reduces to a consistent and high-order accurate solver for the limiting macroscopic equations of the mixture.

A DSMC-PIC coupling method for the Vlasov-Maxwell-Landau system

Andrea Medaglia, Lorenzo Pareschi, Mattia Zanella (preprint arXiv:2510.22226)

We present a numerical framework for the simulation of collisional plasma dynamics, based on a coupling between Direct Simulation Monte Carlo (DSMC) and Particle-in-Cell (PIC) methods for the Vlasov-Maxwell-Landau system. The approach extends previously developed DSMC techniques for the homogeneous Landau equation to the fully inhomogeneous, electromagnetic regime. The Landau collision operator is treated through a stochastic particle formulation inspired by the grazing-collision limit of the Boltzmann equation, which enables an efficient and physically consistent representation of Coulomb interactions without relying on the full Boltzmann structure.

Augmented data and neural networks for robust epidemic forecasting: application to COVID-19 in Italy

Giacomo Dimarco, Federica Ferrarese, Lorenzo Pareschi (preprint arXiv:2510.09192)

In this work, we propose a data augmentation strategy aimed at improving the training phase of neural networks and, consequently, the accuracy of their predictions. Our approach relies on generating synthetic data through a suitable compartmental model combined with the incorporation of uncertainty. The available data are then used to calibrate the model, which is further integrated with deep learning techniques to produce additional synthetic data for training.

Multi-Order Monte Carlo IMEX hierarchies for uncertainty quantification in multiscale hyperbolic systems

Giulia Bertaglia, Walter Boscheri, Lorenzo Pareschi (preprint arXiv:2508.20187)

We introduce a novel Multi-Order Monte Carlo approach for uncertainty quantification in the context of multiscale time-dependent partial differential equations. The new framework leverages Implicit-Explicit Runge-Kutta time integrators to satisfy the asymptotic-preserving property across different discretization orders of accuracy. In contrast to traditional Multi-Level Monte Carlo methods, which require costly hierarchical re-meshing, our method constructs a multi-order hierarchy by varying both spatial and temporal discretization orders within the Monte Carlo framework. This enables efficient variance reduction while naturally adapting to the multiple scales inherent in the problem.

Swarm-based optimization with jumps: a kinetic BGK framework and convergence analysis

 

Giacomo Borghi, Hyesung Im, Lorenzo Pareschi (to appear in Comm. Pure Appl. Analysis. Preprint arXiv:2507.00871)

Metaheuristic algorithms are powerful tools for global optimization, particularly for non-convex and non-differentiable problems where exact methods are often impractical. Particle-based optimization methods, inspired by swarm intelligence principles, have shown effectiveness due to their ability to balance exploration and exploitation within the search space. In this work, we introduce a novel particle-based optimization algorithm where velocities are updated via random jumps, a strategy commonly used to enhance stochastic exploration.

Structure and asymptotic preserving deep neural surrogates for uncertainty quantification in multiscale kinetic equations

Wei Chen, Giacomo Dimarco, Lorenzo Pareschi (to appear in SIAM MMS. Preprint arXiv:2506.10636)

The high dimensionality of kinetic equations with stochastic parameters poses major computational challenges for uncertainty quantification (UQ). Traditional Monte Carlo (MC) sampling methods, while widely used, suffer from slow convergence and high variance, which become increasingly severe as the dimensionality of the parameter space grows. To accelerate MC sampling, we adopt a multiscale control variates strategy that leverages low-fidelity solutions from simplified kinetic models to reduce variance.

Superlinear Drift in Consensus-Based Optimization with Condensation Phenomena

Jonathan Franceschi, Lorenzo Pareschi, Mattia Zanella (preprint arXiv:2506.09001)

Consensus-based optimization (CBO) is a class of metaheuristic algorithms designed for global optimization problems. In the many-particle limit, classical CBO dynamics can be rigorously connected to mean-field equations that ensure convergence toward global minimizers under suitable conditions. In this work, we draw inspiration from recent extensions of the Kaniadakis--Quarati model for indistinguishable bosons to develop a novel CBO method governed by a system of SDEs with superlinear drift and nonconstant diffusion. The resulting mean-field formulation in one dimension exhibits condensation-like phenomena, including finite-time blow-up and loss of L2-regularity.

Robust feedback control of collisional plasma dynamics in presence of uncertainties

Giacomo Albi, Giacomo Dimarco, Federica Ferrarese, Lorenzo Pareschi (preprint arXiv:2505.19992)

Magnetic fusion aims to confine high-temperature plasma within a device, enabling the fusion of deuterium and tritium nuclei to release energy. Due to the very large temperatures involved, it is essential to isolate the plasma from the device walls to prevent structural damage and the external magnetic fields play a fundamental role in achieving this confinement. In realistic settings, the physical mechanisms governing plasma behavior are highly complex, involving numerous uncertain parameters and intricate particle interactions, such as collisions, that significantly affect both confinement efficiency and overall stability.

Hierarchical dynamic domain decomposition for the multiscale Boltzmann equation

Domenico Caparello, Lorenzo Pareschi, Thomas Rey (preprint arXiv:2505.03360)

In this work, we present a hierarchical domain decomposition method for the multi-scale Boltzmann equation based on moment realizability matrices, a concept introduced by Levermore, Morokoff, and Nadiga. This criterion is used to dynamically partition the two-dimensional spatial domain into three regimes: the Euler regime, an intermediate kinetic regime governed by the ES-BGK model, and the full Boltzmann regime. The key advantage of this approach lies in the use of Euler equations in regions where the flow is near hydrodynamic equilibrium, the ES-BGK model in moderately non-equilibrium regions where a fluid description is insufficient but full kinetic resolution is not yet necessary, and the full Boltzmann solver where strong non-equilibrium effects dominate, such as near shocks and boundary layers. This allows for both high accuracy and significant computational savings, as the Euler solver and the ES-BGK models are considerably cheaper than the full kinetic Boltzmann model.

Wasserstein convergence rates for stochastic particle approximation of Boltzmann models

 G. Borghi, L. Pareschi (SIAM J. Numerical Analysis, to appear. Preprint arXiv: 2504.10091). 

We establish quantitative convergence rates for stochastic particle approximation based on Nanbu-type Monte Carlo schemes applied to a broad class of collisional kinetic models. Using coupling techniques and stability estimates in the Wasserstein-1 (Kantorovich-Rubinstein) metric, we derive sharp error bounds that reflect the nonlinear interaction structure of the models. Our framework includes classical Nanbu Monte Carlo method and more recent developments as Time Relaxed Monte Carlo methods. The results bridge the gap between probabilistic particle approximations and deterministic numerical error analysis, and provide a unified perspective for the convergence theory of Monte Carlo methods for Boltzmann-type equations. As a by-product, we also obtain existence and uniqueness of solutions to a large class of Boltzmann-type equations.

A data augmentation strategy for deep neural networks with application to epidemic modelling

M. Awais, A.S. Ali, G. Dimarco, F. Ferrarese, L. Pareschi (Bollettino dell'Unione Matematica Italiana
10.1007/s40574-025-00486-3, preprint arXiv:2502.21033)

In this work, we integrate the predictive capabilities of compartmental disease dynamics models with machine learning ability to analyze complex, high-dimensional data and uncover patterns that conventional models may overlook. Specifically, we present a proof of concept demonstrating the application of data-driven methods and deep neural networks to a recently introduced SIR-type model with social features, including a saturated incidence rate, to improve epidemic prediction and forecasting. Our results show that a robust data augmentation strategy trough suitable data-driven models can improve the reliability of Feed-Forward Neural Networks (FNNs) and Nonlinear Autoregressive Networks (NARs), providing a complementary strategy to Physics-Informed Neural Networks, particularly in settings where data augmentation from mechanistic models can enhance learning.

Multi-fidelity and multi-level Monte Carlo methods for kinetic models of traffic flow

Elisa Iacomini, Lorenzo Pareschi (preprint arXiv:2501.15967)

In traffic flow modeling, incorporating uncertainty is crucial for accurately capturing the complexities of real-world scenarios. In this work we focus on kinetic models of traffic flow, where a key step is to design effective numerical tools for analyzing uncertainties in vehicles interactions. To this end we discuss space-homogeneous Boltzmann-type equations, employing a non intrusive Monte Carlo approach both on the physical space, to solve the kinetic equation, and on the stochastic space, to investigate the uncertainty. To address the high dimensional challenges posed by this coupling, control variate approaches such as multi-fidelity and multi-level Monte Carlo methods are particularly effective.