Reduced-Order Analysis and Prediction of Kinetic Plasma Behavior using Dynamic Mode Decomposition Indranil Nayak a , Mrinal Kumar b , Fernando L. Teixeira a a ElectroScience Laboratory and Department of Electrical and Computer Engineering, The Ohio State University, Columbus, Ohio 43212, USA b Laboratory for Autonomy in Data-Driven and Complex Systems, Department of Mechanical and Aerospace Engineering, The Ohio State University, Columbus, Ohio 43210, USA Abstract A dynamic mode decomposition (DMD) based reduced-order model (ROM) is developed for analysis and prediction of kinetic plasma behavior. DMD is applied to the high-fidelity kinetic plasma model based on the electromagnetic particle-in-cell (EMPIC) algorithm to extract the underlying dynamics and key features of the model. In particular, the ability of DMD to reconstruct the self electric field distribution from high-fidelity data and the effect of DMD extrapolated self-fields on charged particle dynamics are investigated. An in- line sliding-window DMD method is presented for identifying the transition from transient to equilibrium state based on the loci of DMD eigenvalues in the complex plane. The in-line detection of equilibrium state combined with time extrapolation ability of DMD has the potential to effectively expedite the simulation. Case studies involving electron beams and plasma ball are presented to assess the strengths and limitations of the proposed method. Keywords: Equilibrium detection, kinetic plasma, limit cycle detection, particle-in-cell, reduced-order models, dynamic mode decomposition. 1. Introduction Kinetic plasma simulations are important for a wide range of applications, including but not limited to the design and analysis of high-power microwave sources, particle accelerators, laser ignited devices, and ionosphere and magneto- sphere problems [1–7]. Electromagnetic particle-in-cell (EMPIC) algorithms are typically used for simulating kinetic collisionless plasmas governed by Maxwell- Vlasov equations. EMPIC algorithms compute the electromagnetic field on the spatial mesh based on a discretized form of Maxwell’s equations while simul- taneously updating, via a kinetic model based on the Lorentz force equation, * Corresponding author Email address: nayak.77@osu.edu (Indranil Nayak) Preprint submitted to Elsevier 8th November 2021 arXiv:2010.09613v1 [physics.plasm-ph] 19 Oct 2020