No Arabic abstract
Time-centered, hence second-order, methods for integrating the relativistic momentum of charged particles in an electromagnetic field are derived. A new method is found by averaging the momentum before use in the magnetic rotation term, and an implementation is presented that differs from the relativistic Boris Push [1] only in the method for calculating the Lorentz factor. This is shown to have the same second-order accuracy in time as that (Boris Push) [1] found by splitting the electric acceleration and magnetic rotation and that [2] found by averaging the velocity in the magnetic rotation term. All three methods are shown to conserve energy when there is no electric field. The Boris method and the current method are shown to be volume-preserving, while the method of [2] and the current method preserve the $vec{E} times vec{B}$ velocity. Thus, of these second-order relativistic momentum integrations, only the integrator introduced here both preserves volume and gives the correct $vec{E} times vec{B}$ velocity. While all methods have error that is second-order in time, they deviate from each other by terms that increase as the motion becomes relativistic. Numerical results show that [2] develops energy errors near resonant orbits of a test problem that neither volume-preserving integrator does. [1] J. Boris, Relativistic plasma simulation-optimization of a hybrid code, in: Proc. Fourth Conf. Num. Sim. Plasmas, Naval Res. Lab, Wash. DC, 1970, pp. 3-67. [2] J.-L. Vay, Simulation of beams or plasmas crossing at relativistic velocity, Physics of Plasmas (1994-present) 15 (5) (2008) 056701.
The dynamics of charged particles in electromagnetic fields is an essential component of understanding the most extreme environments in our Universe. In electromagnetic fields of sufficient magnitude, radiation emission dominates the particle motion and effects of nonlinear quantum electrodynamics (QED) are crucial, which triggers electron-positron pair cascades and counterintuitive particle-trapping phenomena. As a result of recent progress in laser technology, high-power lasers provide a platform to create and probe such fields in the laboratory. With new large-scale laser facilities on the horizon and the prospect of investigating these hitherto unexplored regimes, we review the basic physical processes of radiation reaction and QED in strong fields, how they are treated theoretically and in simulation, the new collective dynamics they unlock, recent experimental progress and plans, as well as possible applications for high-flux particle and radiation sources.
A class of explicit numerical schemes is developed to solve for the relativistic dynamics and spin of particles in electromagnetic fields, using the Lorentz-BMT equation formulated in the Clifford algebra representation of Baylis. It is demonstrated that these numerical methods, reminiscent of the leapfrog and Verlet methods, share a number of important properties: they are energy-conserving, volume-conserving and second order convergent. These properties are analysed empirically by benchmarking against known analytical solutions in constant uniform electrodynamic fields. It is demonstrated that the numerical error in a constant magnetic field remains bounded for long time simulations in contrast to the Boris pusher, whose angular error increases linearly with time. Finally, the intricate spin dynamics of a particle is investigated in a plane wave field configuration.
When two collisionless plasma shells collide, they interpenetrate and the overlapping region may turn Weibel unstable for some values of the collision parameters. This instability grows magnetic filaments which, at saturation, have to block the incoming flow if a Weibel shock is to form. In a recent paper [J. Plasma Phys. (2016), vol. 82, 905820403], it was found implementing a toy model for the incoming particles trajectories in the filaments, that a strong enough external magnetic field $mathbf{B}_0$ can prevent the filaments to block the flow if it is aligned with. Denoting $B_f$ the peak value of the field in the magnetic filaments, all test particles stream through them if $alpha=B_0/B_f > 1/2$. Here, this result is extended to the case of an oblique external field $B_0$ making an angle $theta$ with the flow. The result, numerically found, is simply $alpha > kappa(theta)/costheta$, where $kappa(theta)$ is of order unity. Noteworthily, test particles exhibit chaotic trajectories.
Electron evaporation plays an important role in the electron temperature evolution and thus expansion rate in low-density ultracold plasmas. In addition, evaporation is useful as a potential tool for obtaining colder electron temperatures and characterizing plasma parameters. Evaporation theory has been developed for atomic gases and has been applied to a one-component plasma system. We numerically investigate whether such an adapted theory is applicable to ultracold neutral plasmas. We find that it is not due to the violation of fundamental assumptions of the model. The details of our calculations are presented as well as a discussion of the implications for a simple description of the electron evaporation rate in ultracold plasmas.
Since the invention of chirped pulse amplification, which was recognized by a Nobel prize in physics in 2018, there has been a continuing increase in available laser intensity. Combined with advances in our understanding of the kinetics of relativistic plasma, studies of laser-plasma interactions are entering a new regime where the physics of relativistic plasmas is strongly affected by strong-field quantum electrodynamics (QED) processes, including hard photon emission and electron-positron ($e^+$-$e^-$) pair production. This coupling of quantum emission processes and relativistic collective particle dynamics can result in dramatically new plasma physics phenomena, such as the generation of dense $e^+$-$e^-$ pair plasma from near vacuum, complete laser energy absorption by QED processes or the stopping of an ultrarelativistic electron beam, which could penetrate a cm of lead, by a hairs breadth of laser light. In addition to being of fundamental interest, it is crucial to study this new regime to understand the next generation of ultra-high intensity laser-matter experiments and their resulting applications, such as high energy ion, electron, positron, and photon sources for fundamental physics studies, medical radiotherapy, and next generation radiography for homeland security and industry.