The Particle-in-Cell Method
The workhorse algorithm that solves the Vlasov-Maxwell system by following sampled particles on a field grid.
Bridging particles and fields
The particle-in-cell (PIC) method solves the kinetic Vlasov-Maxwell (or Vlasov-Poisson) system by representing the distribution function with many macro-particles while computing the fields on a grid. Each step cycles through four operations, closing the self-consistent loop between particles and fields.
- Deposit: scatter particle charge and current onto the grid
- Solve: advance Maxwell's equations (or Poisson) on the grid
- Gather: interpolate the grid fields back to particle positions
- Push: advance particle velocities and positions with the Boris algorithm
Why macro-particles work
Each macro-particle represents a large number of real particles sampled from the distribution. Because the fields are smoothed on the grid, the method captures collective behavior (waves, instabilities, Landau damping) with far fewer particles than reality contains, at the cost of statistical noise that falls as one over the square root of the particle count.
Numerical constraints
- The grid must resolve the Debye length or artificial grid heating occurs
- Explicit schemes must resolve the plasma period for stability
- The Boris pusher conserves phase-space volume, preventing long-term energy drift
Variants
Electromagnetic PIC uses a Yee-grid FDTD Maxwell solver; electrostatic PIC uses a Poisson solve. Implicit PIC relaxes the Debye and plasma-period constraints for larger scales. Gyrokinetic PIC follows gyrocenters instead of full orbits to reach confinement timescales. Monte-Carlo collision operators add collisions to the otherwise collisionless push.
Where it is used
PIC is the standard tool for kinetic problems: wave-plasma interaction, radio-frequency sheaths, laser-plasma, and (in gyrokinetic form) turbulence. Kinetic and gyrokinetic PIC simulations support the physics analysis of Kronos machines, including turbulence and fast-particle studies for the Hyperion breeder, which are design and simulation efforts.