Accurately simulating nine-dimensional phase space of relativistic particles in strong fields
arXiv:2007.07556 · doi:10.1016/j.jcp.2021.110367
Abstract
Next-generation high-power lasers that can be focused to intensities exceeding 10^23 W/cm^2 are enabling new physics and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In strong fields, the motion of charged particles and their spin is affected by radiation reaction. Standard PIC codes usually use Boris or its variants to advance the particles, which requires very small time steps in the strong-field regime to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a 9D particle phase space (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime are desired. We present a new 9D phase space particle pusher based on analytical solutions to the position, momentum and spin advance from the Lorentz force, together with the semi-classical form of RR in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step for each particle from the simulation time step. Due to the analytical integration, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. We present single-particle simulations and full PIC simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 9D phase space for given laser fields. A discussion on the numerical efficiency of the proposed pusher is also provided.
References in corpus (7)
- Generating multi-GeV electron bunches using single stage laser wakefield acceleration in a 3D nonlinear regime
- Radiation Reaction Effects on Radiation Pressure Acceleration
- Ultrarelativistic electron beam polarization in single-shot interaction with an ultraintense laser pulse
- Temporal resolution criterion for correctly simulating relativistic electron motion in a high-intensity laser field
- Particle-In-Cell Modeling of Plasma-Based Accelerators in Two and Three Dimensions
- Elimination of the numerical Cerenkov instability for spectral EM-PIC codes
- Special unitary particle pusher for extreme fields
Cited by in corpus (11)
- Polarized QED cascades
- Kinetic theory for spin-polarized relativistic plasmas
- Special unitary particle pusher for extreme fields
- RaDiO: an efficient spatiotemporal radiation diagnostic for particle-in-cell codes
- Particle acceleration and radiation reaction in a strongly magnetized rotating dipole
- An Analytic Boris Pusher for Plasma Simulation
- Particle motion in neutron stars ultra-strong electromagnetic field: the influence of radiation reaction
- Ultrafast dynamics of a spin-polarized electron plasma with magnetic ions
- Plasma acceleration of polarized particle beams
- General-relativistic and non-ideal radiative cooling in neutron star magnetospheres
- Radiation Reaction Enhancement in Flying Focus Pulses