Abstract:We present the parallel particle-in-cell (PIC) code Apar-T and, more importantly, address the fundamental question of the relations between the PIC model, the Vlasov-Maxwell theory, and real plasmas. First, we present four validation tests: spectra from simulations of thermal plasmas, linear growth rates of the relativistic tearing instability and of the filamentation instability, and nonlinear filamentation merging phase. For the filamentation instability we show that the effective growth rates measured on the total energy can differ by more than 50% from the linear cold predictions and from the fastest modes of the simulation. We link these discrepancies to the superparticle number per cell and to the level of field fluctuations. Second, we detail a new method for initial loading of Maxwell-Jüttner particle distributions with relativistic bulk velocity and relativistic temperature, and explain why the traditional method with individual particle boosting fails. The formulation of the relativistic Harris equilibrium is generalized to arbitrary temperature and mass ratios. Both are required for the tearing instability setup. Third, we turn to the key point of this paper and scrutinize the question of what description of (weakly coupled) physical plasmas is obtained by PIC models. These models rely on two building blocks: coarse-graining, i.e., grouping of the order of p ~ 1010 real particles into a single computer superparticle, and field storage on a grid with its subsequent finite superparticle size. We introduce the notion of coarse-graining dependent quantities, i.e., quantities depending on p. They derive from the PIC plasma parameter ΛPIC, which we show to behave as ΛPIC ∝ 1/p. We explore two important implications. One is that PIC collision- and fluctuation-induced thermalization times are expected to scale with the number of superparticles per grid cell, and thus to be a factor p ~ 1010 smaller than in real plasmas, a fact that we confirm with simulations. The other is that the level of electric field fluctuations scales as 1/ΛPIC ∝ p. We provide a corresponding exact expression, taking into account the finite superparticle size. We confirm both expectations with simulations. Fourth, we compare the Vlasov-Maxwell theory, often used for code benchmarking, to the PIC model. The former describes a phase-space fluid with Λ = + ∞ and no correlations, while the PIC plasma features a small Λ and a high level of correlations when compared to a real plasma. These differences have to be kept in mind when interpreting and validating PIC results against the Vlasov-Maxwell theory and when modeling real physical plasmas.

Diff-PIC: Revolutionizing Particle-In-Cell Nuclear Fusion Simulation with Diffusion Models

A Pairwise Nuclear Fusion Algorithm for Particle-in-cell Simulations: Weighted Particles at Relativistic Energies

Toward a First-Principles Integrated Simulation of Tokamak Edge Plasmas

QuickPIC: a highly efficient fully parallelized PIC code for plasma-based acceleration

PIFE-PIC: Parallel Immersed-Finite-Element Particle-In-Cell For 3-D Kinetic Simulations of Plasma-Material Interactions

Physics-informed Meta-instrument for eXperiments (PiMiX) with applications to fusion energy

Modern Gyrokinetic Particle-In-Cell Simulation of Fusion Plasmas on Top Supercomputers

JefiPIC: A 3-D Full Electromagnetic Particle-in-Cell Simulator Based on Jefimenko's Equations on GPU

Recent development of fully kinetic particle-in-cell method and its application to fusion plasma instability study

A Machine-Learning Augmented Mixed Method for Simulating Α Particle Transport in Inertial Confinement Fusion

Study on Extreme Plasma Dynamics by Quantum Electrodynamic Particle-in-Cell Simulations

A Comparative Study on the Acceleration Techniques for Solving Finite Difference Discretization Poisson's Equation in the PIC/MCC Method

Unraveling Diffusion in Fusion Plasma: A Case Study of In Situ Processing and Particle Sorting

Explicit energy-conserving modification of relativistic PIC method

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Controlling the Numerical Cerenkov Instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction

PM2D: A parallel GPU-based code for the kinetic simulation of laser plasma instabilities in large scale plasmas

Apar-T: code, validation, and physical interpretation of particle-in-cell results

PolyPIC: The Polymorphic-Particle-in-Cell Method for Fluid-Kinetic Coupling

Characterizing the Performance of the Implicit Massively Parallel Particle-in-Cell iPIC3D Code