Abstract:Plasma-neutral interactions, including reactive kinetics, are often either studied in 0D using ODE based de- scriptions, or in multi-dimensional fluid or particle based plasma codes. The latter case involves a complex assembly of procedures that are not always necessary to test effects of underlying physical models and mech- anisms for particle-based descriptions. Here we present ThunderBoltz, a lightweight, publicly available 0D Direct Simulation Monte Carlo (DSMC) code designed to accommodate a generalized combination of species and arbitrary cross sections without the overhead of expensive field solves. It can produce electron, ion, and neutral velocity distributions in applied AC/DC E-field and/or static B-field scenarios. The code is built in the C++ standard library and includes a convenient Python interface that allows for input file generation (including compatibility with cross section data from the LXCat database), electron transport and reaction rate constant post-processing, input parameter constraint satisfaction, calculation scheduling, and diagnostic plotting. These codes can be accessed at the repository: https://github.com/lanl/ThunderBoltz. In this work we compare ThunderBoltz transport calculations against Bolsig+ calculations, benchmark test prob- lems, and swarm experiment data, finding good agreement with all three in the appropriate field regimes. In addition, we present example use cases where the electron, ion, and background neutral particle species are self consistently evolved to model the background kinetics, a feature that is absent in fixed background Monte Carlo and n-term Boltzmann solvers. The latter functionality allows for the possibility of particle- based chemical kinetics simulations of the plasma and neutral species as a new alternative to ODE based approaches.

Benchmark calculations for anisotropic scattering in kinetic models for low temperature plasma

A tutorial overview on the angular scattering models of Electron-Neutral, Ion-Neutral, Neutral-Neutral, and Coulomb Collisions in Monte Carlo collision modeling on low temperature plasma

Particle simulation methods for the Landau-Fokker-Planck equation with uncertain data

Monte Carlo Collision method for low temperature plasma simulation

A multiscale hybrid Maxwellian-Monte-Carlo Coulomb collision algorithm for particle simulations

Benchmarks for infinite medium, time dependent transport problems with isotropic scattering

Monte Carlo simulations of charged particle transport in plasmas: A fast single-scattering model

Numerical Investigation of Non-equilibrium Electron Effects on the Collisional Ionization Rate in the Collisional-Radiative Model

Computation of effective collision strengths for plasma applications using JAC

Monte Carlo Simulation Study of Quasi-Elastic Electron Scattering from an Overlayer/Substrate System

Monte Carlo Simulation Study of Electron Interaction with Solids and Surfaces

Multidisciplinary benchmarks of a conservative spectral solver for the nonlinear Boltzmann equation

Particle simulation of collision dynamics for ion beam injection into a rarefied gas

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker-Planck-Boltzmann Model

A fast solver for the spatially homogeneous electron Boltzmann equation

Electron Scattering in Solids

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker–Planck–Boltzmann Model

ThunderBoltz: An open-source direct simulation Monte Carlo Boltzmann solver for plasma transport, chemical kinetics, and 0D modeling

Historical development of electron swarm physics based on the Boltzmann equation towards in-depth understanding of a low-temperature collisional plasma