Abstract:Coulomb collisions in particle simulations for weakly coupled plasmas are modeled by the Landau-Fokker-Planck equation, which is typically solved by Monte-Carlo (MC) methods. One of the main disadvantages of MC is the timestep accuracy constraint {\nu}{\Delta}t << 1 to resolve the collision frequency {\nu}. The constraint becomes extremely stringent for self-collisions in the presence of high charge state species and for inter-species collisions with large mass disparities (such as present in Inertial Confinement Fusion hohlraums), rendering long-time-scale simulations prohibitively expensive or impractical. To overcome these difficulties, we explore a hybrid Maxwellian-MC (HMMC) model for particle simulations. Specifically, we devise a collisional algorithm that describes weakly collisional species with particles, and highly collisional species and fluid components with Maxwellians. We employ the Lemons method for particle-Maxwellian collisions, enhanced with a more careful treatment of low-relative-speed particles, and a five-moment model for Maxwellian-Maxwellian collisions. Particle-particle binary collisions are dealt with classic Takizuka-Abe MC, which we extend to accommodate arbitrary particle weights to deal with large density disparities without compromising conservation properties. HMMC is strictly conservative and significantly outperforms standard MC methods in situations with large mass disparities among species or large charge states, demonstrating orders of magnitude improvement in computational efficiency. We will substantiate the accuracy and performance of the proposed method with several examples of varying complexity, including both zero-dimensional relaxation and one-dimensional transport problems, the latter using a hybrid kinetic-ion/fluid-electron model.

First principles calculation of the effect of Coulomb collisions in partially ionized gases

A molecular dynamics simulation of charged particle transfer through metals: The role of multi-body collisions

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

First-Principles Determination of Electron-Ion Couplings in the Warm Dense Matter Regime

Mixed, quantum-classical description of electron density transfer in the collision process

Reduction of the Coulomb logarithm due to electron-neutral collisions

A Simple Electron-Density Based Force Field Model for High-Energy Interactions between Atoms and Molecules

Importance of accurate consideration of the electron inertia in hybrid-kinetic simulations of collisionless plasma turbulence: 1. The 2D limit

Demarcating the classical and quantum approaches for the Coulomb logarithm in plasmas

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

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

Numerical investigation of nonequilibrium electron effects on the collisional ionization rate in the collisional-radiative model

Experimental validation of particle-in-cell/Monte Carlo collisions simulations in low-pressure neon capacitively coupled plasmas

Monte Carlo Collision method for low temperature plasma simulation

A Monte Carlo simulation of ion transport at finite temperatures

A Kinetic Theory of Particle Charging in the Free-Molecule Regime

Auxiliary field simulation and Coulomb's Law

Collisional charging of dust particles by suprathermal particles. I -- Standard anisotropic Kappa distributions

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

Effects of fast atoms and energy-dependent secondary electron emission yields in PIC/MCC simulations of capacitively coupled plasmas

Ionisation Calculations using Classical Molecular Dynamics