-
A fully plasma based electron injector for a linear collider or XFEL
Authors:
Thamine N. Dalichaouch,
Xinlu L. Xu,
Fei Li,
Frank S. Tsung,
Warren B. Mori
Abstract:
We demonstrate through high-fidelity particle-in-cell simulations a simple approach for efficiently generating 20+ GeV electron beams with the necessary charge, energy spread, and emittance for use as the injector for an electron arm of a future linear collider or a next generation XFEL. The self-focusing of an unmatched, relatively low quality, drive beam results in self-injection by elongating t…
▽ More
We demonstrate through high-fidelity particle-in-cell simulations a simple approach for efficiently generating 20+ GeV electron beams with the necessary charge, energy spread, and emittance for use as the injector for an electron arm of a future linear collider or a next generation XFEL. The self-focusing of an unmatched, relatively low quality, drive beam results in self-injection by elongating the wakefield excited in the nonlinear blowout regime. Over pump depletion distances, the drive beam dynamics and self-loading from the injected beam leads to extremely high quality and high energy output beams. For plasma densities of $10^{18} \ \text{cm}^{-3}$, PIC simulation results indicate that self-injected beams with $0.52 \ \text{nC}$ of charge can be accelerated to $\sim 20$ GeV energies with projected energy spreads, $\lesssim 1\%$ within the beam core, slice normalized emittances as low as $110 \ \text{nm}$, a peak normalized brightness $\gtrsim 10^{19} \ \text{A}/\text{m}^2/\text{rad}^2$, and energy transfer efficiencies $\gtrsim 54\%$.
△ Less
Submitted 6 June, 2024;
originally announced June 2024.
-
Accurate simulation of direct laser acceleration in a laser wakefield accelerator
Authors:
Kyle G. Miller,
John P. Palastro,
Jessica L. Shaw,
Fei Li,
Frank S. Tsung,
Viktor K. Decyk,
Warren B. Mori
Abstract:
In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The particle-in-cell (PIC) algorithm, althou…
▽ More
In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The particle-in-cell (PIC) algorithm, although often the tool of choice to study DLA, contains inherent errors due to numerical dispersion and the time staggering of the electric and magnetic fields. Further, conventional PIC implementations cannot reliably disentangle the fields of the plasma wave and laser pulse, which obscures interpretation of the dominant acceleration mechanism. Here, a customized field solver that reduces errors from both numerical dispersion and time staggering is used in conjunction with a field decomposition into azimuthal modes to perform PIC simulations of DLA in an LWFA. Comparisons with traditional PIC methods, model equations, and experimental data show improved accuracy with the customized solver and convergence with an order-of-magnitude fewer cells. The azimuthal-mode decomposition reveals that the most energetic electrons receive comparable energy from DLA and LWFA.
△ Less
Submitted 22 March, 2023;
originally announced March 2023.
-
Ultra-compact attosecond X-ray free-electron lasers utilizing unique beams from plasma-based acceleration and an optical undulator
Authors:
Xinlu Xu,
Jiaxin Liu,
Thamine Dalichaouch,
Frank S. Tsung,
Zhen Zhang,
Zhirong Huang,
Mark J. Hogan,
Xueqing Yan,
Chan Joshi,
Warren B. Mori
Abstract:
Accelerator-based X-ray free-electron lasers (XFELs) are the latest addition to the revolutionary tools of discovery for the 21st century. The two major components of an XFEL are an accelerator-produced electron beam and a magnetic undulator which tend to be kilometer-scale long and expensive. Here, we present an ultra-compact scheme to produce 10s of attosecond X-ray pulses with several GW peak p…
▽ More
Accelerator-based X-ray free-electron lasers (XFELs) are the latest addition to the revolutionary tools of discovery for the 21st century. The two major components of an XFEL are an accelerator-produced electron beam and a magnetic undulator which tend to be kilometer-scale long and expensive. Here, we present an ultra-compact scheme to produce 10s of attosecond X-ray pulses with several GW peak power utilizing a novel aspect of the FEL instability using a highly chirped, pre-bunched and ultra-bright electron beam from a plasma-based accelerator interacting with an optical undulator. The self-selection of electrons from the combination of a highly chirped and pre-bunched beam leads to the stable generation of attosecond X-ray pulses. Furthermore, two-color attosecond pulses with sub-femtosecond separation can be produced by adjusting the energy distribution of the electron beam so that multiple FEL resonances occur at different locations within the beam. Such a tunable coherent attosecond X-ray sources may open up a new area of attosecond science enabled by X-ray attosecond pump/probe techniques.
△ Less
Submitted 17 February, 2023;
originally announced February 2023.
-
Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition
Authors:
Fei Li,
Weiming An,
Frank S. Tsung,
Viktor K. Decyk,
Warren B. Mori
Abstract:
High-fidelity modeling of plasma-based acceleration (PBA) requires the use of 3D fully nonlinear and kinetic descriptions based on the particle-in-cell (PIC) method. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the beam-plasma interaction. In a QSA-PIC algorithm, the plasma response to a charged particle beam or…
▽ More
High-fidelity modeling of plasma-based acceleration (PBA) requires the use of 3D fully nonlinear and kinetic descriptions based on the particle-in-cell (PIC) method. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the beam-plasma interaction. In a QSA-PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA-PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA-PIC algorithm that combines another speedup technique known as azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can significantly reduce the algorithmic complexity. Modeling the laser-plasma interaction in a full 3D PIC algorithm is very computationally expensive due to the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible with the QSA PIC algorithms based on the azimuthal mode expansion. This algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described.
△ Less
Submitted 7 September, 2022; v1 submitted 22 March, 2022;
originally announced March 2022.
-
Ultra-Bright Electron Bunch Injection in a Plasma Wakefield Driven by a Superluminal Flying Focus Electron Beam
Authors:
Fei Li,
Thamine N. Dalichaouch,
Jacob R. Pierce,
Xinlu Xu,
Frank S. Tsung,
Wei Lu,
Chan Joshi,
Warren B. Mori
Abstract:
We propose a new method for self-injection of high-quality electron bunches in the plasma wakefield structure in the blowout regime utilizing a "flying focus" produced by a drive beam with an energy chirp. In a flying focus the speed of the density centroid of the drive bunch can be superluminal or subluminal by utilizing the chromatic dependence of the focusing optics. We first derive the focal v…
▽ More
We propose a new method for self-injection of high-quality electron bunches in the plasma wakefield structure in the blowout regime utilizing a "flying focus" produced by a drive beam with an energy chirp. In a flying focus the speed of the density centroid of the drive bunch can be superluminal or subluminal by utilizing the chromatic dependence of the focusing optics. We first derive the focal velocity and the characteristic length of the focal spot in terms of the focal length and an energy chirp. We then demonstrate using multidimensional particle-in-cell simulations that a wake driven by a superluminally propagating flying focus of an electron beam can generate GeV-level electron bunches with ultralow normalized slice emittance ($\sim$30 nm rad), high current ($\sim$ 17 kA), low slice energy-spread ($\sim$0.1%) and therefore high normalized brightness ($>10^{19}$ A/rad$^2$/m$^2$) in a plasma of density $\sim10^{19}$ cm$^{-3}$. The injection process is highly controllable and tunable by changing the focal velocity and shaping the drive beam current. Near-term experiments at FACET II where the capabilities to generate tens of kA, <10 fs drivers are planned, could potentially produce beams with brightness near $10^{20}$ A/rad$^2$/m$^2$.
△ Less
Submitted 23 April, 2022; v1 submitted 30 July, 2021;
originally announced August 2021.
-
A multi-sheath model for highly nonlinear plasma wakefields
Authors:
Thamine N. Dalichaouch,
Xinlu Xu,
Adam Tableman,
Fei Li,
Frank S. Tsung,
Warren B. Mori
Abstract:
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 pushing plasma electrons radially outward, creating an ion bubble surrounded by a sheath of electrons characterized by the source term…
▽ More
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 pushing 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}(ρ-J_z/c)$ where $ρ$ 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.
△ Less
Submitted 24 March, 2021;
originally announced March 2021.
-
Generation and acceleration of high brightness electrons beams bunched at X-ray wavelengths using plasma-based acceleration
Authors:
Xinlu Xu,
Fei Li,
Frank S. Tsung,
Kyle Miller,
Vitaly Yakimenko,
Mark J. Hogan,
Chan Joshi,
Warren B. Mori
Abstract:
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 in nonlinear plasma wave wakefields. The density modulation turns on and off the injection of electrons at the period of the modulation. Due to the uni…
▽ More
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 in 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.8\micro\meter$, 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.
△ Less
Submitted 30 October, 2020;
originally announced October 2020.
-
Accurately simulating nine-dimensional phase space of relativistic particles in strong fields
Authors:
Fei Li,
Viktor K. Decyk,
Kyle G. Miller,
Adam Tableman,
Frank S. Tsung,
Marija Vranic,
Ricardo A. Fonseca,
Warren B. Mori
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 ch…
▽ More
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.
△ Less
Submitted 21 April, 2021; v1 submitted 15 July, 2020;
originally announced July 2020.
-
A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm
Authors:
Fei Li,
Kyle G. Miller,
Xinlu Xu,
Frank S. Tsung,
Viktor K. Decyk,
Weiming An,
Ricardo A. Fonseca,
Warren B. Mori
Abstract:
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) dispe…
▽ More
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 Faraday's and Ampere's 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.
△ Less
Submitted 7 April, 2020;
originally announced April 2020.
-
On numerical errors to the fields surrounding a relativistically moving particle in PIC codes
Authors:
Xinlu Xu,
Fei Li,
Frank S. Tsung,
Thamine N. Dalichaouch,
Weiming An,
Han Wen,
Viktor K. Decyk,
Ricardo A. Fonseca,
Mark J. Hogan,
Warren B. Mori
Abstract:
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, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves…
▽ More
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, inertial 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.
△ Less
Submitted 29 October, 2019;
originally announced October 2019.
-
Three-dimensional particle-in-cell modeling of parametric instabilities near the quarter-critical density in plasmas
Authors:
H. Wen,
A. V. Maximov,
R. Yan,
J. Li,
C. Ren,
F. S. Tsung
Abstract:
The nonlinear regime of laser-plasma interaction including both two-plasmon--decay (TPD) and stimulated Raman scattering (SRS) instabilities has been studied in three-dimensional (3-D) particle-in-cell simulations with parameters relevant to the inertial confinement fusion (ICF) experiments. SRS and TPD develop in the same region in plasmas, and the generation of fast electrons can be described ac…
▽ More
The nonlinear regime of laser-plasma interaction including both two-plasmon--decay (TPD) and stimulated Raman scattering (SRS) instabilities has been studied in three-dimensional (3-D) particle-in-cell simulations with parameters relevant to the inertial confinement fusion (ICF) experiments. SRS and TPD develop in the same region in plasmas, and the generation of fast electrons can be described accurately with only the full model including both SRS and TPD. The growth of instabilities in the linear stage is found to be in good agreement with analytical theories. In the saturation stage the enhanced low-frequency density perturbations driven by the daughter waves of the SRS sidescattering can saturate the TPD and consequently inhibit the fast-electron generation. The fast-electron flux in 3-D modeling is up to an order of magnitude smaller than previously reported in 2-D TPD simulations, bringing it close to the results of ICF experiments.
△ Less
Submitted 15 April, 2019; v1 submitted 19 June, 2018;
originally announced June 2018.
-
Optimizing Laser Wakefield Acceleration in the Nonlinear Self-Guided Regime for Fixed Laser Energy
Authors:
A. Davidson,
A. Tableman,
P. Yu,
W. An,
F. S. Tsung,
W. Lu,
R. A. Fonseca,
W. B. Mori
Abstract:
The scaling laws for laser wakefield acceleration in the nonlinear, self-guided regime [Lu et al. Phys. Rev. Spec. Top. Accel. Beams 10, 061301 (2007)] are examined in detail using the quasi-3D version of the particle-in-cell code OSIRIS. We find that the scaling laws continue to work well as the plasma density is reduced while the normalized laser amplitude is kept fixed. For fixed laser energy,…
▽ More
The scaling laws for laser wakefield acceleration in the nonlinear, self-guided regime [Lu et al. Phys. Rev. Spec. Top. Accel. Beams 10, 061301 (2007)] are examined in detail using the quasi-3D version of the particle-in-cell code OSIRIS. We find that the scaling laws continue to work well as the plasma density is reduced while the normalized laser amplitude is kept fixed. For fixed laser energy, the energy gain of an isolated bunch of electrons can be improved with some loss in the bunch charge by shortening the normalized pulse length until self-guiding no longer occurs, and through the use of asymmetric longitudinal profiles with rapid rise times. For example, without any external guiding a 15 J,.8μm laser with a pulse length of 46 fs (39 fs) is found to generate a quasi-mono-energetic bunch of 355pC (227pC) with a max energy of 3.25 (4.04) GeV with an acceleration distance of 2.43 cm (3.08cm). Furthermore, a bunch with 39.4pC and a maximum energy 4.6 GeV is produced for an asymmetric laser with a rapid rise time. Studies for 30 J and 100 J lasers are also presented.
△ Less
Submitted 1 May, 2019; v1 submitted 22 May, 2018;
originally announced May 2018.
-
Quenching of stimulated Raman scattering in the kinetic regime by external magnetic fields
Authors:
B. J. Winjum,
F. S. Tsung,
W. B. Mori
Abstract:
We show via particle-in-cell simulations that small normalized magnetic fields ($ω_c/ω_p \ll 1$) applied perpendicularly to a light wave can significantly modify the evolution of backward stimulated Raman scattering (SRS) in the kinetic regime. The presence of the magnetic field increases the threshold for kinetic inflation and decreases the amount of reflectivity when SRS is driven significantly…
▽ More
We show via particle-in-cell simulations that small normalized magnetic fields ($ω_c/ω_p \ll 1$) applied perpendicularly to a light wave can significantly modify the evolution of backward stimulated Raman scattering (SRS) in the kinetic regime. The presence of the magnetic field increases the threshold for kinetic inflation and decreases the amount of reflectivity when SRS is driven significantly above threshold. Analysis indicates this arises because trapped electrons are accelerated as they surf across the wave, leading to the continual dissipation of the electron plasma waves over a wider range of wave amplitudes. The simulation parameters are directly relevant for SRS in inertial confinement fusion devices and indicate that approximately 30 Tesla magnetic fields might significantly reduce SRS backscatter.
△ Less
Submitted 16 February, 2018;
originally announced February 2018.
-
Controlling the Numerical Cerenkov Instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction
Authors:
Fei Li,
Peicheng Yu,
Xinlu Xu,
Frederico Fiuza,
Viktor K. Decyk,
Thamine Dalichaouch,
Asher Davidson,
Adam Tableman,
Weiming An,
Frank S. Tsung,
Ricardo A. Fonseca,
Wei Lu,
Warren B. Mori
Abstract:
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 pl…
▽ More
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.
△ Less
Submitted 5 May, 2016;
originally announced May 2016.
-
Lorentz boosted frame simulation of Laser wakefield acceleration in quasi-3D geometry
Authors:
Peicheng Yu,
Xinlu Xu,
Asher Davidson,
Adam Tableman,
Thamine Dalichaouch,
Michael D. Meyers,
Frank S. Tsung,
Viktor K. Decyk,
Frederico Fiuza,
Jorge Vieira,
Ricardo A. Fonseca,
Wei Lu,
Luis O. Silva,
Warren B. Mori
Abstract:
When modeling laser wakefield acceleration (LWFA) using the particle-in-cell (PIC) algorithm in a Lorentz boosted frame, the plasma is drifting relativistically at $β_b c$ towards the laser, which can lead to a computational speedup of $\sim γ_b^2=(1-β_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…
▽ More
When modeling laser wakefield acceleration (LWFA) using the particle-in-cell (PIC) algorithm in a Lorentz boosted frame, the plasma is drifting relativistically at $β_b c$ towards the laser, which can lead to a computational speedup of $\sim γ_b^2=(1-β_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.
△ Less
Submitted 24 October, 2015;
originally announced October 2015.
-
Self-modulated laser wakefield accelerators as x-ray sources
Authors:
N. Lemos,
J. L. Martins,
F. S. Tsung,
J. L. Shaw,
K. A. Marsh,
F. Albert,
B. B. Pollock,
C. Joshi
Abstract:
The development of a directional, small-divergence, and short-duration picosecond x-ray probe beam with an energy greater than 50 keV is desirable for high energy density science experiments. We therefore explore through particle-in-cell (PIC) computer simulations the possibility of using x-rays radiated by betatron-like motion of electrons from a self-modulated laser wakefield accelerator as a po…
▽ More
The development of a directional, small-divergence, and short-duration picosecond x-ray probe beam with an energy greater than 50 keV is desirable for high energy density science experiments. We therefore explore through particle-in-cell (PIC) computer simulations the possibility of using x-rays radiated by betatron-like motion of electrons from a self-modulated laser wakefield accelerator as a possible candidate to meet this need. Two OSIRIS 2D PIC simulations with mobile ions are presented, one with a normalized vector potential a0 = 1.5 and the other with an a0 = 3. We find that in both cases direct laser acceleration (DLA) is an important additional acceleration mechanism in addition to the longitudinal electric field of the plasma wave. Together these mechanisms produce electrons with a continuous energy spectrum with a maximum energy of 300 MeV for a0 = 3 case and 180 MeV in the a0 = 1.5 case. Forward-directed x-ray radiation with a photon energy up to 100 keV was calculated for the a0 = 3 case and up to 12 keV for the a0 = 1.5 case. The x-ray spectrum can be fitted with a sum of two synchrotron spectra with critical photon energy of 13 and 45 keV for the a0 of 3 and critical photon energy of 0.3 and 1.4 keV for a0 of 1.5 in the plane of polarization of the laser. The full width at half maximum divergence angle of the x-rays was 62 x 1.9 mrad for a0 = 3 and 77 x 3.8 mrad for a0 = 1.5.
△ Less
Submitted 5 October, 2015;
originally announced October 2015.
-
Estimation of direct laser acceleration in laser wakefield accelerators using particle-in-cell simulations
Authors:
J. L. Shaw,
N. Lemos,
K. A. Marsh,
F. S. Tsung,
W. B. Mori,
C. Joshi
Abstract:
Many current laser wakefield acceleration (LWFA) experiments are carried out in a regime where the laser pulse length is on the order of or longer than the wake wavelength and where ionization injection is employed to inject electrons into the wake. In these experiments, the trapped electrons will co-propagate with the longitudinal wakefield and the transverse laser field. In this scenario, the el…
▽ More
Many current laser wakefield acceleration (LWFA) experiments are carried out in a regime where the laser pulse length is on the order of or longer than the wake wavelength and where ionization injection is employed to inject electrons into the wake. In these experiments, the trapped electrons will co-propagate with the longitudinal wakefield and the transverse laser field. In this scenario, the electrons can gain a significant amount of energy from both the direct laser acceleration (DLA) mechanism as well as the usual LWFA mechanism. Particle-in-cell (PIC) codes are frequently used to discern the relative contribution of these two mechanisms. However, if the longitudinal resolution used in the PIC simulations is inadequate, it can produce numerical heating that can overestimate the transverse motion, which is important in determining the energy gain due to DLA. We have therefore carried out a systematic study of this LWFA regime by varying the longitudinal resolution of PIC simulations from the standard, best-practice resolution of 30 points per laser wavelength to four times that value and then examining the energy gain characteristics of both the highest-energy electrons and the bulk electrons. By calculating the contribution of DLA to the final energies of the electrons produced from the LWFA, we find that although the transverse momentum and oscillation radii are over-estimated in the lower-resolution simulations, this over-estimation does not lead to artificial energy gain by DLA. Rather, the DLA contribution to the highest-energy electrons is larger in the higher-resolution cases because the DLA resonance is better maintained. Thus, even at the highest longitudinal resolutions, DLA contributes a significant portion of the energy gained by the highest-energy electrons and also contributes to accelerating the bulk of the charge in the electron beam produced by the LWFA.
△ Less
Submitted 25 September, 2015;
originally announced September 2015.
-
Satisfying the Direct Laser Acceleration Resonance Condition in a Laser Wakefield Accelerator
Authors:
J. L. Shaw,
N. Vafaei-Najafabadi,
K. A. Marsh,
N. Lemos,
F. S. Tsung,
W. B Mori,
C. Joshi
Abstract:
In this proceeding, we show that when the drive laser pulse overlaps the trapped electrons in a laser wakefield accelerator (LWFA), those electrons can gain energy from direct laser acceleration (DLA) over extended distances despite the evolution of both the laser and the wake. Through simulations, the evolution of the properties of both the laser and the electron beam is quantified, and then the…
▽ More
In this proceeding, we show that when the drive laser pulse overlaps the trapped electrons in a laser wakefield accelerator (LWFA), those electrons can gain energy from direct laser acceleration (DLA) over extended distances despite the evolution of both the laser and the wake. Through simulations, the evolution of the properties of both the laser and the electron beam is quantified, and then the resonance condition for DLA is examined in the context of this change. We find that although the electrons produced from the LWFA cannot continuously satisfy the DLA resonance condition, they nevertheless can gain a significant amount of energy from DLA.
△ Less
Submitted 26 February, 2015;
originally announced February 2015.
-
Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit
Authors:
Peicheng Yu,
Xinlu Xu,
Adam Tableman,
Viktor K. Decyk,
Frank S. Tsung,
Frederico Fiuza,
Asher Davidson,
Jorge Vieira,
Ricardo A. Fonseca,
Wei Lu,
Luis O. Silva,
Warren B. Mori
Abstract:
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 numeri…
▽ More
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.
△ Less
Submitted 4 February, 2015;
originally announced February 2015.
-
Elimination of the numerical Cerenkov instability for spectral EM-PIC codes
Authors:
Peicheng Yu,
Xinlu Xu,
Viktor K. Decyk,
F. Fiuza,
Jorge Vieira,
Frank S. Tsung,
Ricardo A. Fonseca,
Wei Lu,
Luis O. Silva,
Warren B. Mori
Abstract:
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 electromagnetic waves on a grid to wave-particle resonances, including aliased resonances, i.e., $ω+ 2πμ/Δt=(k_1+ 2πν_1/Δx_1)v_0$, where $μ$ and $ν_1$ refer…
▽ More
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 electromagnetic waves on a grid to wave-particle resonances, including aliased resonances, i.e., $ω+ 2πμ/Δt=(k_1+ 2πν_1/Δx_1)v_0$, where $μ$ and $ν_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 ($μ=ν_1=0$) can be seen. (Please check the whole abstract in the paper).
△ Less
Submitted 31 January, 2015; v1 submitted 1 July, 2014;
originally announced July 2014.
-
Proton acceleration by a relativistic laser frequency-chirp driven plasma snowplow
Authors:
Aakash A. Sahai,
T. C. Katsouleas,
R. A. Bingham,
F. S. Tsung,
A. R. Tableman,
M. Tzoufras,
W. B. Mori
Abstract:
We analyze the use of a relativistic laser pulse with a controlled frequency chirp incident on a rising plasma density gradient to drive an acceleration structure for proton and light-ion acceleration. The Chirp Induced Transparency Acceleration (ChITA) scheme is described with an analytical model of the velocity of the snowplow at critical density on a pre-formed rising plasma density gradient th…
▽ More
We analyze the use of a relativistic laser pulse with a controlled frequency chirp incident on a rising plasma density gradient to drive an acceleration structure for proton and light-ion acceleration. The Chirp Induced Transparency Acceleration (ChITA) scheme is described with an analytical model of the velocity of the snowplow at critical density on a pre-formed rising plasma density gradient that is driven by a positive-chirp in the frequency of a relativistic laser pulse. The velocity of the ChITA-snowplow is shown to depend upon rate of rise of the frequency of the relativistic laser pulse represented by $\frac{ε_0}θ$ where, $ε_0 = \frac{Δω_0}{ω_0}$ and chirping spatial scale-length, $θ$, the normalized magnetic vector potential of the laser pulse $a_0$ and the plasma density gradient scale-length, $α$. We observe using 1-D OSIRIS simulations the formation and forward propagation of ChITA-snowplow, being continuously pushed by the chirping laser at a velocity in accordance with the analytical results. The trace protons reflect off of this propagating snowplow structure and accelerate mono-energetically. The control over ChITA-snowplow velocity allows the tuning of accelerated proton energies.
△ Less
Submitted 16 May, 2014;
originally announced May 2014.
-
Long Term Evolution of Plasma Wakefields
Authors:
Aakash A. Sahai,
T. C. Katsouleas,
F. S. Tsung,
W. B. Mori
Abstract:
We study the long-term evolution (LTE) of plasma wakefields over multiple plasma-electron periods and few plasma-ion periods, much less than a recombination time. The evolution and relaxation of such a wakefield-perturbed plasma over these timescales has important implications for the upper limits of repetition-rates in plasma colliders. Intense fields in relativistic lasers (or intense beams) cre…
▽ More
We study the long-term evolution (LTE) of plasma wakefields over multiple plasma-electron periods and few plasma-ion periods, much less than a recombination time. The evolution and relaxation of such a wakefield-perturbed plasma over these timescales has important implications for the upper limits of repetition-rates in plasma colliders. Intense fields in relativistic lasers (or intense beams) create plasma wakefields (modes around ωpe) by transferring energy to the plasma electrons. Charged-particle beams in the right phase may be accelerated with acceleration/focusing gradients of tens of GeV/m. However, wakefields leave behind a plasma not in equilibrium, with a relaxation time of multiple plasma-electron periods. Ion motion over ion timescales, caused by energy transfer from the driven plasma-electrons to the plasma-ions can create interesting plasma states. Eventually during LTE, the dynamics of plasma de-coheres (multiple modes through instability driven mixing), thermalizing into random motion (second law of thermodynamics), dissipating energy away from the wakefields. Wakefield-drivers interacting with such a relativistically hot-plasma lead to plasma wakefields that differ from the wakefields in a cold-plasma.
△ Less
Submitted 16 May, 2014;
originally announced May 2014.
-
Role of direct laser acceleration in energy gained by electrons in a laser wakefield accelerator with ionization injection
Authors:
J L Shaw,
F S Tsung,
N Vafaei-Najafabadi,
K A Marsh,
N Lemos,
W B Mori,
C Joshi
Abstract:
We have investigated the role that the transverse electric field of the laser plays in the acceleration of electrons in a laser wakefield accelerator (LWFA) operating in the quasi-blowout regime through particle-in-cell code simulations. In order to ensure that longitudinal compression and/or transverse focusing of the laser pulse is not needed before the wake can self-trap the plasma electrons, w…
▽ More
We have investigated the role that the transverse electric field of the laser plays in the acceleration of electrons in a laser wakefield accelerator (LWFA) operating in the quasi-blowout regime through particle-in-cell code simulations. In order to ensure that longitudinal compression and/or transverse focusing of the laser pulse is not needed before the wake can self-trap the plasma electrons, we have employed the ionization injection technique. Furthermore, the plasma density is varied such that at the lowest densities, the laser pulse occupies only a fraction of the first wavelength of the wake oscillation (the accelerating bucket), whereas at the highest density, the same duration laser pulse fills the entire first bucket. Although the trapped electrons execute betatron oscillations due to the ion column in all cases, at the lowest plasma density they do not interact with the laser field and the energy gain is all due to the longitudinal wakefield. However, as the density is increased, there can be a significant contribution to the maximum energy due to direct laser acceleration (DLA) of those electrons that undergo betatron motion in the plane of the polarization of the laser pulse. Eventually, DLA can be the dominant energy gain mechanism over acceleration due to the longitudinal field at the highest densities.
△ Less
Submitted 19 April, 2014;
originally announced April 2014.
-
Implementation of a hybrid particle code with a PIC description in r-z and a gridless description in $φ$ into OSIRIS
Authors:
A. Davidson,
A. Tableman,
W. An,
F. S. Tsung,
W. Lu,
J. Vieira,
R. A. Fonseca,
L. O. Silva,
W. B. Mori
Abstract:
For many plasma physics problems, three-dimensional and kinetic effects are very important. However, such simulations are very computationally intensive. Fortunately, there is a class of problems for which there is nearly azimuthal symmetry and the dominant three-dimensional physics is captured by the inclusion of only a few azimuthal harmonics. Recently, it was proposed [A. Lifschitz et al., J. C…
▽ More
For many plasma physics problems, three-dimensional and kinetic effects are very important. However, such simulations are very computationally intensive. Fortunately, there is a class of problems for which there is nearly azimuthal symmetry and the dominant three-dimensional physics is captured by the inclusion of only a few azimuthal harmonics. Recently, it was proposed [A. Lifschitz et al., J. Comp. Phys. 228 (5) (2009) 1803-1814] to model one such problem, laser wakefield acceleration, by expanding the fields and currents in azimuthal harmonics and truncating the expansion after only the first harmonic. The complex amplitudes of the fundamental and first harmonic for the fields were solved on an r-z grid and a procedure for calculating the complex current amplitudes for each particle based on its motion in Cartesian geometry was presented using a Marder's correction to maintain the validity of Gauss's law. In this paper, we describe an implementation of this algorithm into OSIRIS using a rigorous charge conserving current deposition method to maintain the validity of Gauss's law. We show that this algorithm is a hybrid method which uses a particles-in-cell description in r-z and a gridless description in $φ$. We include the ability to keep an arbitrary number of harmonics and higher order particle shapes. Examples, for laser wakefield acceleration, plasma wakefield acceleration, and beam loading are also presented and directions for future work are discussed.
△ Less
Submitted 26 March, 2014;
originally announced March 2014.
-
Relativistically Induced Transparency Acceleration (RITA) of Protons and Light-ions with Ultrashort Laser Interaction with Heavy-ion Plasma Density Gradient
Authors:
Aakash A. Sahai,
F. S. Tsung,
A. R. Tableman,
W. B. Mori,
T. C. Katsouleas
Abstract:
The relativistically induced transparency acceleration (RITA) scheme of proton and ion acceleration using laser-plasma interactions is introduced, modeled, and compared to the existing schemes. Protons are accelerated with femtosecond relativistic pulses to produce quasimonoenergetic bunches with controllable peak energy. The RITA scheme works by a relativistic laser inducing transparency to densi…
▽ More
The relativistically induced transparency acceleration (RITA) scheme of proton and ion acceleration using laser-plasma interactions is introduced, modeled, and compared to the existing schemes. Protons are accelerated with femtosecond relativistic pulses to produce quasimonoenergetic bunches with controllable peak energy. The RITA scheme works by a relativistic laser inducing transparency to densities higher than the cold-electron critical density, while the background heavy ions are stationary. The rising laser pulse creates a traveling acceleration structure at the relativistic critical density by ponderomotively driving a local electron density inflation, creating an electron snowplow and a co-propagating electrostatic potential. The snowplow advances with a velocity determined by the rate of the rise of the laser's intensity envelope and the heavy-ion-plasma density gradient scale length. The rising laser is incrementally rendered transparent to higher densities such that the relativistic-electron plasma frequency is resonant with the laser frequency. In the snowplow frame, trace density protons reflect off the electrostatic potential and get snowplowed, while the heavier background ions are relatively unperturbed. Quasimonoenergetic bunches of velocity equal to twice the snowplow velocity can be obtained and tuned by controlling the snowplow velocity using laser-plasma parameters. An analytical model for the proton energy as a function of laser intensity, rise time, and plasma density gradient is developed and compared to 1D and 2D PIC OSIRIS simulations. We model the acceleration of protons to GeV energies with tens-of-femtoseconds laser pulses of a few petawatts. The scaling of proton energy with laser power compares favorably to other mechanisms for ultrashort pulses.
△ Less
Submitted 21 February, 2014;
originally announced February 2014.
-
Modeling of Laser wakefield acceleration in Lorentz boosted frame using EM-PIC code with spectral solver
Authors:
Peicheng Yu,
Xinlu Xu,
Viktor K. Decyk,
Weiming An,
Jorge Vieira,
Frank S. Tsung,
Ricardo A. Fonseca,
Wei Lu,
Luis O. Silva,
Warren B. Mori
Abstract:
Simulating laser wakefield acceleration (LWFA) in a Lorentz boosted frame in which the plasma drifts towards the laser with $v_b$ can speedup the simulation by factors of $γ^2_b=(1-v^2_b/c^2)^{-1}$. In these simulations the relativistic drifting plasma inevitably induces a high frequency numerical instability that contaminates the interested physics. Various approaches have been proposed to mitiga…
▽ More
Simulating laser wakefield acceleration (LWFA) in a Lorentz boosted frame in which the plasma drifts towards the laser with $v_b$ can speedup the simulation by factors of $γ^2_b=(1-v^2_b/c^2)^{-1}$. In these simulations the relativistic drifting plasma inevitably induces a high frequency numerical instability that contaminates the interested physics. Various approaches have been proposed to mitigate this instability. One approach is to solve Maxwell equations in Fourier space (a spectral solver) as this has been shown to suppress the fastest growing modes of this instability in simple test problems using a simple low pass, ring (in two dimensions), or shell (in three dimensions) filter in Fourier space. We describe the development of a fully parallelized, multi-dimensional, particle-in-cell code that uses a spectral solver to solve Maxwell's equations and that includes the ability to launch a laser using a moving antenna. This new EM-PIC code is called UPIC-EMMA and it is based on the components of the UCLA PIC framework (UPIC). We show that by using UPIC-EMMA, LWFA simulations in the boosted frames with arbitrary $γ_b$ can be conducted without the presence of the numerical instability. We also compare the results of a few LWFA cases for several values of $γ_b$, including lab frame simulations using OSIRIS, a EM-PIC code with a finite difference time domain (FDTD) Maxwell solver. These comparisons include cases in both linear, and nonlinear regimes. We also investigate some issues associated with numerical dispersion in lab and boosted frame simulations and between FDTD and spectral solvers.
△ Less
Submitted 9 October, 2013;
originally announced October 2013.
-
Snowmass 2013 Computing Frontier: Accelerator Science
Authors:
P. Spentzouris,
E. Cormier-Michel,
C. Joshi,
J. Amundson,
W. An,
D. L. Bruhwiler,
J. R. Cary,
B. Cowan,
V. K. Decyk,
E. Esarey,
R. A. Fonseca,
A. Friedman,
C. G. R. Geddes,
D. P. Grote,
I. Kourbanis,
W. P. Leemans,
W. Lu,
W. B. Mori,
C. Ng,
Ji Qiang,
T. Roberts,
R. D. Ryne,
C. B. Schroeder,
L. O. Silva,
F. S. Tsung
, et al. (2 additional authors not shown)
Abstract:
This is the working summary of the Accelerator Science working group of the Computing Frontier of the Snowmass meeting 2013. It summarizes the computing requirements to support accelerator technology in both Energy and Intensity Frontiers.
This is the working summary of the Accelerator Science working group of the Computing Frontier of the Snowmass meeting 2013. It summarizes the computing requirements to support accelerator technology in both Energy and Intensity Frontiers.
△ Less
Submitted 8 October, 2013;
originally announced October 2013.
-
Exploiting multi-scale parallelism for large scale numerical modelling of laser wakefield accelerators
Authors:
Ricardo A. Fonseca,
Jorge Vieira,
Frederico Fiúza,
Asher Davidson,
Frank S. Tsung,
Warren B. Mori,
Luís O. Silva
Abstract:
A new generation of laser wakefield accelerators, supported by the extreme accelerating fields generated in the interaction of PW-Class lasers and underdense targets, promises the production of high quality electron beams in short distances for multiple applications. Achieving this goal will rely heavily on numerical modeling for further understanding of the underlying physics and identification o…
▽ More
A new generation of laser wakefield accelerators, supported by the extreme accelerating fields generated in the interaction of PW-Class lasers and underdense targets, promises the production of high quality electron beams in short distances for multiple applications. Achieving this goal will rely heavily on numerical modeling for further understanding of the underlying physics and identification of optimal regimes, but large scale modeling of these scenarios is computationally heavy and requires efficient use of state-of-the-art Petascale supercomputing systems. We discuss the main difficulties involved in running these simulations and the new developments implemented in the OSIRIS framework to address these issues, ranging from multi-dimensional dynamic load balancing and hybrid distributed / shared memory parallelism to the vectorization of the PIC algorithm. We present the results of the OASCR Joule Metric program on the issue of large scale modeling of LWFA, demonstrating speedups of over 1 order of magnitude on the same hardware. Finally, scalability to over $\sim 10^6$ cores, and sustained performance over $\sim 2$ PFlops is demonstrated, opening the way for large scale modeling of laser wakefield accelerator scenarios.
△ Less
Submitted 3 October, 2013;
originally announced October 2013.
-
Numerical instability due to relativistic plasma drift in EM-PIC simulations
Authors:
Xinlu Xu,
Peicheng Yu,
Samual F. Martins,
Frank S. Tsung,
Viktor K. Decyk,
Ricardo A. Fonseca,
Wei Lu,
Luis O. Silva,
Warren B. Mori
Abstract:
The numerical instability observed in the Electromagnetic-Particle-in-cell (EM-PIC) simulations with a plasma drifting with relativistic velocities is studied using both theory and computer simulations. We derive the numerical dispersion relation for a cold plasma drifting with a relativistic velocity and find an instability attributed to the coupling between the beam modes of the drifting plasma…
▽ More
The numerical instability observed in the Electromagnetic-Particle-in-cell (EM-PIC) simulations with a plasma drifting with relativistic velocities is studied using both theory and computer simulations. We derive the numerical dispersion relation for a cold plasma drifting with a relativistic velocity and find an instability attributed to the coupling between the beam modes of the drifting plasma and the electromagnetic modes in the system. The characteristic pattern of the instability in Fourier space for various simulation setups and Maxwell Equation solvers are explored by solving the corresponding numerical dispersion relations. Furthermore, based upon these characteristic patterns we derive an asymptotic expression for the instability growth rate. The asymptotic expression greatly speeds up the calculation of instability growth rate and makes the parameter scan for minimal growth rate feasible even for full three dimensions. The results are compared against simulation results and good agreement is found. These results can be used as a guide to develop possible approaches to mitigate the instability. We examine the use of a spectral solver and show that such a solver when combined with a low pass filter with a cutoff value of $|\vec{k}|$ essentially eliminates the instability while not modifying modes of physical interest. The use of spectral solver also provides minimal errors to electromagnetic modes in the lowest Brillouin zones.
△ Less
Submitted 29 November, 2012; v1 submitted 5 November, 2012;
originally announced November 2012.
-
Anomalously Hot Electrons due to Rescatter of Stimulated Raman Scattering in the Kinetic Regime
Authors:
B. J. Winjum,
J. E. Fahlen,
F. S. Tsung,
W. B. Mori
Abstract:
Using particle-in-cell simulations, we examine hot electron generation from electron plasma waves excited by stimulated Raman scattering and rescattering in the kinetic regime where the wavenumber times the Debye length (kλ_D) is greater than 0.3 for backscatter. We find that for laser and plasma conditions of possible relevance to experiments at the National Ignition Facility (NIF), anomalously e…
▽ More
Using particle-in-cell simulations, we examine hot electron generation from electron plasma waves excited by stimulated Raman scattering and rescattering in the kinetic regime where the wavenumber times the Debye length (kλ_D) is greater than 0.3 for backscatter. We find that for laser and plasma conditions of possible relevance to experiments at the National Ignition Facility (NIF), anomalously energetic electrons can be produced through the interaction of a discrete spectrum of plasma waves generated from SRS (back and forward scatter), rescatter, and the Langmuir decay of the rescatter-generated plasma waves. Electrons are bootstrapped in energy as they propagate into plasma waves with progressively higher phase velocities.
△ Less
Submitted 3 October, 2012;
originally announced October 2012.
-
Convective Raman Amplification of Light Pulses Causing Kinetic Inflation in Inertial Fusion Plasmas
Authors:
Ian N. Ellis,
David J. Strozzi,
Benjamin J. Winjum,
Frank S. Tsung,
Thomas Grismayer,
Warren B. Mori,
Jay E. Fahlen,
Edward A. Williams
Abstract:
We perform 1D particle-in-cell (PIC) simulations using OSIRIS, which model a short-duration (~500/ω0 FWHM) scattered light seed pulse in the presence of a constant counter-propagating pump laser with an intensity far below the absolute instability threshold. The seed undergoes linear convective Raman amplification and dominates over fluctuations due to particle discreteness. Our simulation results…
▽ More
We perform 1D particle-in-cell (PIC) simulations using OSIRIS, which model a short-duration (~500/ω0 FWHM) scattered light seed pulse in the presence of a constant counter-propagating pump laser with an intensity far below the absolute instability threshold. The seed undergoes linear convective Raman amplification and dominates over fluctuations due to particle discreteness. Our simulation results are in good agreement with results from a coupled mode solver when we take into account special relativity and the use of finite size PIC simulation particles. We present linear gain spectra including both effects. Extending the PIC simulations past when the seed exits the simulation domain reveals bursts of large-amplitude scattering in many cases, which does not occur in simulations without the seed pulse. These bursts can have amplitudes several times greater than the amplified seed pulse, and we demonstrate that this large-amplitude scattering is the result of kinetic inflation by examining trapped particle orbits. This large-amplitude scattering is caused by the seed modifying the distribution function earlier in the simulation. We perform some simulations with longer duration seeds, which lead to parts of the seeds undergoing kinetic inflation and reaching amplitudes several times more than the steady-state linear theory results. Simulations with continuous seeds demonstrate that the onset of inflation depends on seed wavelength and incident intensity, and we observe oscillations in the reflectivity at a frequency equal to the difference between the seed frequency and the frequency at which the inflationary SRS grows.
△ Less
Submitted 26 November, 2012; v1 submitted 8 August, 2012;
originally announced August 2012.
-
One-to-one direct modeling of experiments and astrophysical scenarios: pushing the envelope on kinetic plasma simulations
Authors:
R. A. Fonseca,
S. F. Martins,
L. O. Silva,
J. W. Tonge,
F. S. Tsung,
W. B. Mori
Abstract:
There are many astrophysical and laboratory scenarios where kinetic effects play an important role. These range from astrophysical shocks and plasma shell collisions, to high intensity laser-plasma interactions, with applications to fast ignition and particle acceleration. Further understanding of these scenarios requires detailed numerical modelling, but fully relativistic kinetic codes are com…
▽ More
There are many astrophysical and laboratory scenarios where kinetic effects play an important role. These range from astrophysical shocks and plasma shell collisions, to high intensity laser-plasma interactions, with applications to fast ignition and particle acceleration. Further understanding of these scenarios requires detailed numerical modelling, but fully relativistic kinetic codes are computationally intensive, and the goal of one-to-one direct modelling of such scenarios and direct comparison with experimental results is still difficult to achieve. In this paper we discuss the issues involved in performing kinetic plasma simulations of experiments and astrophysical scenarios, focusing on what needs to be achieved for one-to-one direct modeling, and the computational requirements involved. We focus on code efficiency and new algorithms, specifically on parallel scalability issues, namely on dynamic load balancing, and on high-order interpolation and boosted frame simulations to optimize simulation performance. We also discuss the new visualization and data mining tools required for these numerical experiments and recent simulation work illustrating these techniques is also presented.
△ Less
Submitted 14 October, 2008;
originally announced October 2008.
-
Beam loading in the nonlinear regime of plasma-based acceleration
Authors:
M. Tzoufras,
W. Lu,
F. S. Tsung,
C. Huang,
W. B. Mori,
T. Katsouleas,
J. Vieira,
R. A. Fonseca,
L. O. Silva
Abstract:
A theory that describes how to load negative charge into a nonlinear, three-dimensional plasma wakefield is presented. In this regime, a laser or an electron beam blows out the plasma electrons and creates a nearly spherical ion channel, which is modified by the presence of the beam load. Analytical solutions for the fields and the shape of the ion channel are derived. It is shown that very high…
▽ More
A theory that describes how to load negative charge into a nonlinear, three-dimensional plasma wakefield is presented. In this regime, a laser or an electron beam blows out the plasma electrons and creates a nearly spherical ion channel, which is modified by the presence of the beam load. Analytical solutions for the fields and the shape of the ion channel are derived. It is shown that very high beam-loading efficiency can be achieved, while the energy spread of the bunch is conserved. The theoretical results are verified with the Particle-In-Cell code OSIRIS.
△ Less
Submitted 1 September, 2008;
originally announced September 2008.
-
One-to-one full scale simulations of laser wakefield acceleration using QuickPIC
Authors:
J. Vieira,
F. Fiuza,
R. A. Fonseca,
L. O. Silva,
C. Huang,
W. Lu,
M. Tzoufras,
F. S. Tsung,
V. Decyk,
W. B. Mori,
J. Cooley,
T. Antonsen Jr
Abstract:
We use the quasi-static particle-in-cell code QuickPIC to perform full-scale, one-to-one LWFA numerical experiments, with parameters that closely follow current experimental conditions. The propagation of state-of-the-art laser pulses in both preformed and uniform plasma channels is examined. We show that the presence of the channel is important whenever the laser self-modulations do not dominat…
▽ More
We use the quasi-static particle-in-cell code QuickPIC to perform full-scale, one-to-one LWFA numerical experiments, with parameters that closely follow current experimental conditions. The propagation of state-of-the-art laser pulses in both preformed and uniform plasma channels is examined. We show that the presence of the channel is important whenever the laser self-modulations do not dominate the propagation. We examine the acceleration of an externally injected electron beam in the wake generated by 10 J laser pulses, showing that by using ten-centimeter-scale plasma channels it is possible to accelerate electrons to more than 4 GeV. A comparison between QuickPIC and 2D OSIRIS is provided. Good qualitative agreement between the two codes is found, but the 2D full PIC simulations fail to predict the correct laser and wakefield amplitudes.
△ Less
Submitted 23 June, 2008;
originally announced June 2008.
-
Three-dimensional structure of the laser wakefield accelerator in the blowout regime
Authors:
J. F. Vieira,
S. F. Martins,
F. Fiuza,
R. A. Fonseca,
L. O. Silva,
C. Huang,
W. Lu,
M. Tzoufras,
F. S. Tsung,
W. B. Mori,
J. Cooley,
T. Antonsen Jr
Abstract:
Three-dimensional Particle-in-Cell (PIC) simulations with the code QuickPIC are used to illustrate the typical accelerating structures associated with the interaction of an intense laser beam with an underdense plasma in the blowout regime. Our simulations are performed with an externally injected electron beam, positioned in the region of maximum accelerating gradients. As the laser propagates…
▽ More
Three-dimensional Particle-in-Cell (PIC) simulations with the code QuickPIC are used to illustrate the typical accelerating structures associated with the interaction of an intense laser beam with an underdense plasma in the blowout regime. Our simulations are performed with an externally injected electron beam, positioned in the region of maximum accelerating gradients. As the laser propagates in the plasma, almost complete electron cavitation occurs, leading to the generation of accelerating fields in excess of 1 GeV/cm.
△ Less
Submitted 10 March, 2008;
originally announced March 2008.
-
Generating multi-GeV electron bunches using single stage laser wakefield acceleration in a 3D nonlinear regime
Authors:
W. Lu,
M. Tzoufras,
C. Joshi,
F. S. Tsung,
W. B. Mori,
J. Vieira,
R. A. Fonseca,
L. O. Silva
Abstract:
The extraordinary ability of space-charge waves in plasmas to accelerate charged particles at gradients that are orders of magnitude greater than in current accelerators has been well documented. We develop a phenomenological framework for Laser WakeField Acceleration (LWFA) in the 3D nonlinear regime, in which the plasma electrons are expelled by the radiation pressure of a short pulse laser, l…
▽ More
The extraordinary ability of space-charge waves in plasmas to accelerate charged particles at gradients that are orders of magnitude greater than in current accelerators has been well documented. We develop a phenomenological framework for Laser WakeField Acceleration (LWFA) in the 3D nonlinear regime, in which the plasma electrons are expelled by the radiation pressure of a short pulse laser, leading to nearly complete blowout. Our theory provides a recipe for designing a LWFA for given laser and plasma parameters and estimates the number and the energy of the accelerated electrons whether self-injected or externally injected. These formulas apply for self-guided as well as externally guided pulses (e.g. by plasma channels). We demonstrate our results by presenting a sample Particle-In-Cell (PIC) simulation of a 30f sec, 200T W laser interacting with a 0.75cm long plasma with density 1.5*10^18 cm^-3 to produce an ultra-short (10f s) mono-energetic bunch of self-injected electrons at 1.5 GeV with 0.3nC of charge. For future higher-energy accelerator applications we propose a parameter space, that is distinct from that described by Gordienko and Pukhov [Physics of Plasmas 12, 043109 (2005)] in that it involves lower densities and wider spot sizes while keeping the intensity relatively constant. We find that this helps increase the output electron beam energy while keeping the efficiency high.
△ Less
Submitted 22 December, 2006;
originally announced December 2006.
-
A global simulation for laser driven MeV electrons in $50μm$-diameter fast ignition targets
Authors:
C. Ren,
M. Tzoufras,
J. Tonge,
W. B. Mori,
F. S. Tsung,
M. Fiore,
R. A. Fonseca,
L. O. Silva,
J. -C. Adam,
A. Heron
Abstract:
The results from 2.5-dimensional Particle-in-Cell simulations for the interaction of a picosecond-long ignition laser pulse with a plasma pellet of 50-$μm$ diameter and 40 critical density are presented. The high density pellet is surrounded by an underdense corona and is isolated by a vacuum region from the simulation box boundary. The laser pulse is shown to filament and create density channel…
▽ More
The results from 2.5-dimensional Particle-in-Cell simulations for the interaction of a picosecond-long ignition laser pulse with a plasma pellet of 50-$μm$ diameter and 40 critical density are presented. The high density pellet is surrounded by an underdense corona and is isolated by a vacuum region from the simulation box boundary. The laser pulse is shown to filament and create density channels on the laser-plasma interface. The density channels increase the laser absorption efficiency and help generate an energetic electron distribution with a large angular spread. The combined distribution of the forward-going energetic electrons and the induced return electrons is marginally unstable to the current filament instability. The ions play an important role in neutralizing the space charges induced by the the temperature disparity between different electron groups. No global coalescing of the current filaments resulted from the instability is observed, consistent with the observed large angular spread of the energetic electrons.
△ Less
Submitted 25 January, 2006;
originally announced January 2006.