Learning a general class of admissible multi-species collision operators from molecular dynamics
Yue Zhao, Andrew Christlieb, Huan Lei
Abstract
We develop a structure-preserving, data-driven collision operator for spatially homogeneous multi-species kinetic systems from molecular dynamics (MD). The operator consists of diagonal self-collision blocks and ordered off-diagonal cross-species blocks to describe intra- and inter-species momentum and energy exchange. Within a local and point-wise identifiable kernel class, we develop the necessary and sufficient condition for the admissible kernel class satisfying the conservation laws, the H-theorem, and the frame indifference. Unlike the classical Landau operator, the off-diagonal kernels are not restricted to be symmetric under permutation of the two velocity variables. This unique structural freedom captures the distinct responses of different species to unresolved correlations and many-body effects arising from micro-scale particle interactions. The equivalent parameterizable kernel formalization enables us to learn a generalized data-driven collision operator directly from MD, where the low-rank tensor representations and random sampling are used to achieve efficient kernel training and numerical simulation. Numerical experiments show that the learned operator accurately predicts transport coefficients and the non-equilibrium relaxation, while retaining discrete conservation and entropy production. In particular, it captures plasma kinetics in the moderately coupled regime, where the predictions of both the Landau and the data-driven model restricted to velocity-permutation symmetry show significant discrepancies.
Create a lesson
Related papers
Learning Spectral-Like Mesh-Free Discretisations
Lucas Gerken Starepravo, Henry Broadley, Steven Lind et al.
RBMD 2.0: Random batch molecular dynamics package for large-scale simulations on multi-GPU architectures
Qi Zhou, Yongfa Guo, Jincheng Zhong et al.
Physics-informed Learning for Orbital Uncertainty Propagation with Error Bounds
Chun-Wei Kong, Morteza Lahijanian, Jay McMahon
A Non-Isothermal Viscoplastic Constitutive Model for Clay Slip Surfaces
Saeed Tourchi, Ehsan Badakhshan, Milad Jabbarzadeh et al.
Mixture of Polyconvex Neural Potentials for Parametric Hyperelasticity: Towards Foundation Material Models
Steven J. Yang, Govinda Anantha Padmanabha, D. Thomas Seidl et al.
OpenMUSTANC (MUltiple Scattering Theory At Nanoplasmonic Cavities): A MATLAB toolbox for the simulation of Plasmonic Sphere Aggregates
Xin Zheng, Christos Mystilidis, Christos Tserkezis et al.