An improved description for nonlinear plasma wakefields with phase velocities near the speed of light is presented and compared against fully kinetic particle-in-cell simulations. These wakefields are excited by intense particle beams or lasers pushi
ng plasma electrons radially outward, creating an ion bubble surrounded by a sheath of electrons characterized by the source term $S equiv -frac{1}{en_p}(rho-J_z/c)$ where $rho$ and $J_z$ are the charge and axial current densities. Previously, the sheath source term was described phenomenologically with a positive-definite function, resulting in a positive definite wake potential. In reality, the wake potential is negative at the rear of the ion column which is important for self-injection and accurate beam loading models. To account for this, we introduce a multi-sheath model in which the source term, $S$, of the plasma wake can be negative in regions outside the ion bubble. Using this model, we obtain a new expression for the wake potential and a modified differential equation for the bubble radius. Numerical results obtained from these equations are validated against particle-in-cell simulations for unloaded and loaded wakes. The new model provides accurate predictions of the shape and duration of trailing bunch current profiles that flatten plasma wakefields. It is also used to design a trailing bunch for a desired longitudinally varying loaded wakefield. We present beam loading results for laser wakefields and discuss how the model can be improved for laser drivers in future work. Finally, we discuss differences between the predictions of the multi- and single-sheath models for beam loading.
We show using particle-in-cell (PIC) simulations and theoretical analysis that a high-quality electron beam whose density is modulated at angstrom scales can be generated directly using density downramp injection in a periodically modulated density i
n nonlinear plasma wave wakefields. The density modulation turns on and off the injection of electrons at the period of the modulation. Due to the unique longitudinal mapping between the electrons initial positions and their final trapped positions inside the wake, this results in an electron beam with density modulation at a wavelength orders of magnitude shorter than the plasma density modulation. The ponderomotive force of two counter propagating lasers of the same frequency can generate a density modulation at half the laser wavelength. Assuming a laser wavelength of $0.8micrometer$, fully self-consistent OSIRIS PIC simulations show that this scheme can generate high quality beams modulated at wavelengths between 10s and 100 angstroms. Such beams could produce fully coherent, stable, hundreds of GW X-rays by going through a resonant undulator.
When a fs duration and hundreds of kA peak current electron beam traverses the vacuum and high-density plasma interface a new process, that we call relativistic transition radiation (R-TR) generates an intense $sim100$ as pulse containing $sim$ TW po
wer of coherent VUV radiation accompanied by several smaller fs duration satellite pulses. This pulse inherits the radial polarization of the incident beam field and has a ring intensity distribution. This R-TR is emitted when the beam density is comparable to the plasma density and the spot size much larger than the plasma skin depth. Physically, it arises from the return current or backward relativistic motion of electrons starting just inside the plasma that Doppler up-shifts the emitted photons. The number of R-TR pulses is determined by the number of groups of plasma electrons that originate at different depths within the first plasma wake period and emit coherently before phase mixing.
A customized finite-difference field solver for the particle-in-cell (PIC) algorithm that provides higher fidelity for wave-particle interactions in intense electromagnetic waves is presented. In many problems of interest, particles with relativistic
energies interact with intense electromagnetic fields that have phase velocities near the speed of light. Numerical errors can arise due to (1) dispersion errors in the phase velocity of the wave, (2) the staggering in time between the electric and magnetic fields and between particle velocity and position and (3) errors in the time derivative in the momentum advance. Errors of the first two kinds are analyzed in detail. It is shown that by using field solvers with different $mathbf{k}$-space operators in Faradays and Amperes law, the dispersion errors and magnetic field time-staggering errors in the particle pusher can be simultaneously removed for electromagnetic waves moving primarily in a specific direction. The new algorithm was implemented into OSIRIS by using customized higher-order finite-difference operators. Schemes using the proposed solver in combination with different particle pushers are compared through PIC simulation. It is shown that the use of the new algorithm, together with an analytic particle pusher (assuming constant fields over a time step), can lead to accurate modeling of the motion of a single electron in an intense laser field with normalized vector potentials, $eA/mc^2$, exceeding $10^4$ for typical cell sizes and time steps.
The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, i
nertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles are generated and interact with plasmas. However, when relativistic particles stream across the grid both in vacuum and in plasma there are many numerical issues that may arise which can lead to incorrect physics. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented. Two types of errors to these expressions are identified. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. The second arises from errors to the denominator of the integrand and lead to Cerenkov like radiation in vacuum. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS. Results from the customized solver are also presented. This solution will be useful for a beam of particles that all move in one direction with a small angular divergence.
A new method of controllable injection to generate high quality electron bunches in the nonlinear blowout regime driven by electron beams is proposed and demonstrated using particle-in-cell simulations. Injection is facilitated by decreasing the wake
phase velocity through varying the spot size of the drive beam and can be tuned through the Courant-Snyder (CS) parameters. Two regimes are examined. In the first, the spot size is focused according to the vacuum CS beta function, while in the second, it is focused by the plasma ion column. The effects of the driver intensity and vacuum CS parameters on the wake velocity and injected beam parameters are examined via theory and simulations. For plasma densities of $sim 10^{19} ~text{cm}^{-3}$, particle-in-cell (PIC) simulations demonstrate that peak normalized brightnesses $gtrsim 10^{20}~text{A}/text{m}^2/text{rad}^2$ can be obtained with projected energy spreads of $lesssim 1%$ within the middle section of the injected beam, and with normalized slice emittances as low as $sim 10 ~text{nm}$.
In this paper we present a customized finite-difference-time-domain (FDTD) Maxwell solver for the particle-in-cell (PIC) algorithm. The solver is customized to effectively eliminate the numerical Cerenkov instability (NCI) which arises when a plasma
(neutral or non-neutral) relativistically drifts on a grid when using the PIC algorithm. We control the EM dispersion curve in the direction of the plasma drift of a FDTD Maxwell solver by using a customized higher order finite difference operator for the spatial derivative along the direction of the drift ($hat 1$ direction). We show that this eliminates the main NCI modes with moderate $vert k_1 vert$, while keeps additional main NCI modes well outside the range of physical interest with higher $vert k_1 vert$. These main NCI modes can be easily filtered out along with first spatial aliasing NCI modes which are also at the edge of the fundamental Brillouin zone. The customized solver has the possible advantage of improved parallel scalability because it can be easily partitioned along $hat 1$ which typically has many more cells than other directions for the problems of interest. We show that FFTs can be performed locally to current on each partition to filter out the main and first spatial aliasing NCI modes, and to correct the current so that it satisfies the continuity equation for the customized spatial derivative. This ensures that Gauss Law is satisfied. We present simulation examples of one relativistically drifting plasmas, of two colliding relativistically drifting plasmas, and of nonlinear laser wakefield acceleration (LWFA) in a Lorentz boosted frame that show no evidence of the NCI can be observed when using this customized Maxwell solver together with its NCI elimination scheme.
When modeling laser wakefield acceleration (LWFA) using the particle-in-cell (PIC) algorithm in a Lorentz boosted frame, the plasma is drifting relativistically at $beta_b c$ towards the laser, which can lead to a computational speedup of $sim gamma_
b^2=(1-beta_b^2)^{-1}$. Meanwhile, when LWFA is modeled in the quasi-3D geometry in which the electromagnetic fields and current are decomposed into a limited number of azimuthal harmonics, speedups are achieved by modeling three dimensional problems with the computation load on the order of two dimensional $r-z$ simulations. Here, we describe how to combine the speed ups from the Lorentz boosted frame and quasi-3D algorithms. The key to the combination is the use of a hybrid Yee-FFT solver in the quasi-3D geometry that can be used to effectively eliminate the Numerical Cerenkov Instability (NCI) that inevitably arises in a Lorentz boosted frame due to the unphysical coupling of Langmuir modes and EM modes of the relativistically drifting plasma in these simulations. In addition, based on the space-time distribution of the LWFA data in the lab and boosted frame, we propose to use a moving window to follow the drifting plasma to further reduce the computational load. We describe the details of how the NCI is eliminated for the quasi-3D geometry, the setups for simulations which combine the Lorentz boosted frame and quasi-3D geometry, the use of a moving window, and compare the results from these simulations against their corresponding lab frame cases. Good agreement is obtained, particularly when there is no self-trapping, which demonstrates it is possible to combine the Lorentz boosted frame and the quasi-3D algorithms when modeling LWFA to achieve unprecedented speedups.
A hybrid Maxwell solver for fully relativistic and electromagnetic (EM) particle-in-cell (PIC) codes is described. In this solver, the EM fields are solved in $k$ space by performing an FFT in one direction, while using finite difference operators in
the other direction(s). This solver eliminates the numerical Cerenkov radiation for particles moving in the preferred direction. Moreover, the numerical Cerenkov instability (NCI) induced by the relativistically drifting plasma and beam can be eliminated using this hybrid solver by applying strategies that are similar to those recently developed for pure FFT solvers. A current correction is applied for the charge conserving current deposit to correctly account for the EM calculation in hybrid Yee-FFT solver. A theoretical analysis of the dispersion properties in vacuum and in a drifting plasma for the hybrid solver is presented, and compared with PIC simulations with good agreement obtained. This hybrid solver is applied to both 2D and 3D Cartesian and quasi-3D (in which the fields and current are decomposed into azimuthal harmonics) geometries. Illustrative results for laser wakefield accelerator simulation in a Lorentz boosted frame using the hybrid solver in the 2D Cartesian geometry are presented, and compared against results from 2D UPIC-EMMA simulation which uses a pure spectral Maxwell solver, and from OSIRIS 2D lab frame simulation using the standard Yee solver. Very good agreement is obtained which demonstrates the feasibility of using the hybrid solver for high fidelity simulation of relativistically drifting plasma with no evidence of the numerical Cerenkov instability.
When using an electromagnetic particle-in-cell (EM-PIC) code to simulate a relativistically drifting plasma, a violent numerical instability known as the numerical Cerenkov instability (NCI) occurs. The NCI is due to the unphysical coupling of electr
omagnetic waves on a grid to wave-particle resonances, including aliased resonances, i.e., $omega + 2pimu/Delta t=(k_1+ 2pi u_1/Delta x_1)v_0$, where $mu$ and $ u_1$ refer to the time and space aliases and the plasma is drifting relativistically at velocity $v_0$ in the $hat{1}$-direction. Recent studies have shown that an EM-PIC code which uses a spectral field solver and a low pass filter can eliminate the fastest growing modes of the NCI. Based on these studies a new spectral PIC code for studying laser wakefield acceleration (LWFA) in the Lorentz boosted frame was developed. However, we show that for parameters of relevance for LWFA simulations in the boosted frame, a relativistically drifting plasma is susceptible to a host of additional unstable modes with lower growth rates, and that these modes appear when the fastest growing unstable modes are filtered out. We show that these modes are most easily identified as the coupling between modes which are purely transverse (EM) and purely longitudinal (Langmuir) in the rest frame of the plasma for specific time and space aliases. We rewrite the dispersion relation of the drifting plasma for a general field solver and obtain analytic expressions for the location and growth rate for each unstable mode, i.e, for each time and space aliased resonances. We show for the spectral solver that when the fastest growing mode is eliminated a new mode at the fundamental resonance ($mu= u_1=0$) can be seen. (Please check the whole abstract in the paper).
