KEMPIC-3D: A transparent and extensible electromagnetic Particle-in-Cell framework for kinetic plasma simulations
Abstract
KEMPIC-3D, a fully electromagnetic three-dimensional particle-in-cell (PIC) code, has been developed in C/C++ and systematically verified for the self-consistent simulation of kinetic plasma phenomena. The numerical framework combines a staggered Yee finite-difference time-domain solver for Maxwell's equations with a relativistic Boris particle pusher, trilinear field interpolation, and charge-conserving current deposition. Its modular architecture emphasizes algorithmic transparency and extensibility, allowing the principal numerical routines to be directly inspected and modified to incorporate additional physical models and computational capabilities. The implementation is assessed through a hierarchy of benchmark problems addressing both the individual numerical components and the self-consistent particle--field coupling. The particle pusher is verified against analytical solutions for relativistic charged-particle motion, while the electromagnetic field solver is evaluated through guided-wave propagation and grid-convergence analysis. The complete PIC algorithm is subsequently verified through plasma oscillations, global energy conservation, and discrete charge continuity, demonstrating the numerical accuracy, conservation properties, and self-consistency of the coupled formulation. The capabilities of KEMPIC-3D are further demonstrated through a representative simulation of laser-driven plasma wakefield generation in an underdense plasma. The simulations reproduce the characteristic nonlinear wakefield structure and exhibit stable numerical behavior across different particle discretizations without additional spatial filtering. Together, these results establish KEMPIC-3D as a reliable, transparent, and extensible computational framework for multidimensional kinetic plasma simulations.
Related
- An Extensible Agentic Framework for Expert-Level Automation of Atomistic Simulations
- OPALX: A performance-portable framework for particle accelerator simulations
- NEOPIC: A Neural Operator Framework for Particle-based Kinetic Plasma Simulations
- GWPTools: An extensible Wolfram Language framework for time-dependent Gaussian wave packet dynamics
- X-dispersionless solver for electromagnetic and axion fields in a cylindrical particle-in-cell code
BibTeX
@misc{lopez2026kempic3d,
title = {{KEMPIC-3D: A transparent and extensible electromagnetic Particle-in-Cell framework for kinetic plasma simulations}},
author = {Jesús E. López and Keren C. Vanegas and Eduardo A. Orozco-Ospino},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.16575},
eprint = {2609.16575},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.16575},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - KEMPIC-3D: A transparent and extensible electromagnetic Particle-in-Cell framework for kinetic plasma simulations
AU - Jesús E. López
AU - Keren C. Vanegas
AU - Eduardo A. Orozco-Ospino
PY - 2026
DA - 2026/09/15
PB - arXiv
DO - 10.48550/arxiv.2609.16575
UR - https://arxiv.org/abs/2609.16575
AB - KEMPIC-3D, a fully electromagnetic three-dimensional particle-in-cell (PIC) code, has been developed in C/C++ and systematically verified for the self-consistent simulation of kinetic plasma phenomena. The numerical framework combines a staggered Yee finite-difference time-domain solver for Maxwell's equations with a relativistic Boris particle pusher, trilinear field interpolation, and charge-conserving current deposition. Its modular architecture emphasizes algorithmic transparency and extensibility, allowing the principal numerical routines to be directly inspected and modified to incorporate additional physical models and computational capabilities. The implementation is assessed through a hierarchy of benchmark problems addressing both the individual numerical components and the self-consistent particle--field coupling. The particle pusher is verified against analytical solutions for relativistic charged-particle motion, while the electromagnetic field solver is evaluated through guided-wave propagation and grid-convergence analysis. The complete PIC algorithm is subsequently verified through plasma oscillations, global energy conservation, and discrete charge continuity, demonstrating the numerical accuracy, conservation properties, and self-consistency of the coupled formulation. The capabilities of KEMPIC-3D are further demonstrated through a representative simulation of laser-driven plasma wakefield generation in an underdense plasma. The simulations reproduce the characteristic nonlinear wakefield structure and exhibit stable numerical behavior across different particle discretizations without additional spatial filtering. Together, these results establish KEMPIC-3D as a reliable, transparent, and extensible computational framework for multidimensional kinetic plasma simulations.
ER -