Accurate modeling of plasma acceleration with arbitrary order pseudospectral particleincell methods
Abstract
Particle in Cell (PIC) simulations are a widely used tool for the investigation of both laser and beamdriven plasma acceleration. It is a known issue that the beam quality can be artificially degraded by numerical Cherenkov radiation (NCR) resulting primarily from an incorrectly modeled dispersion relation. Pseudospectral solvers featuring infinite order stencils can strongly reduce NCR  or even suppress it  and are therefore well suited to correctly model the beam properties. For efficient parallelization of the PIC algorithm, however, localized solvers are inevitable. Arbitrary order pseudospectral methods provide this needed locality. Yet, these methods can again be prone to NCR. Here in this paper, we show that acceptably low solver orders are sufficient to correctly model the physics of interest, while allowing for parallel computation by domain decomposition.
 Authors:
 Univ. of Hamburg (Germany). Center for FreeElectron Laser Science & Dept. of Physics
 Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
 Publication Date:
 Research Org.:
 Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
 Sponsoring Org.:
 USDOE Office of Science (SC), High Energy Physics (HEP) (SC25)
 OSTI Identifier:
 1393128
 Alternate Identifier(s):
 OSTI ID: 1349384
 Grant/Contract Number:
 AC0205CH11231
 Resource Type:
 Journal Article: Accepted Manuscript
 Journal Name:
 Physics of Plasmas
 Additional Journal Information:
 Journal Volume: 24; Journal Issue: 3; Journal ID: ISSN 1070664X
 Publisher:
 American Institute of Physics (AIP)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Finite difference time domain calculations; Particleincell method; Dispersion relations; Fourier transforms; Maxwell equations
Citation Formats
Jalas, S., Dornmair, I., Lehe, R., Vincenti, H., Vay, J. L., Kirchen, M., and Maier, A. R. Accurate modeling of plasma acceleration with arbitrary order pseudospectral particleincell methods. United States: N. p., 2017.
Web. doi:10.1063/1.4978569.
Jalas, S., Dornmair, I., Lehe, R., Vincenti, H., Vay, J. L., Kirchen, M., & Maier, A. R. Accurate modeling of plasma acceleration with arbitrary order pseudospectral particleincell methods. United States. doi:10.1063/1.4978569.
Jalas, S., Dornmair, I., Lehe, R., Vincenti, H., Vay, J. L., Kirchen, M., and Maier, A. R. Mon .
"Accurate modeling of plasma acceleration with arbitrary order pseudospectral particleincell methods". United States.
doi:10.1063/1.4978569. https://www.osti.gov/servlets/purl/1393128.
@article{osti_1393128,
title = {Accurate modeling of plasma acceleration with arbitrary order pseudospectral particleincell methods},
author = {Jalas, S. and Dornmair, I. and Lehe, R. and Vincenti, H. and Vay, J. L. and Kirchen, M. and Maier, A. R.},
abstractNote = {Particle in Cell (PIC) simulations are a widely used tool for the investigation of both laser and beamdriven plasma acceleration. It is a known issue that the beam quality can be artificially degraded by numerical Cherenkov radiation (NCR) resulting primarily from an incorrectly modeled dispersion relation. Pseudospectral solvers featuring infinite order stencils can strongly reduce NCR  or even suppress it  and are therefore well suited to correctly model the beam properties. For efficient parallelization of the PIC algorithm, however, localized solvers are inevitable. Arbitrary order pseudospectral methods provide this needed locality. Yet, these methods can again be prone to NCR. Here in this paper, we show that acceptably low solver orders are sufficient to correctly model the physics of interest, while allowing for parallel computation by domain decomposition.},
doi = {10.1063/1.4978569},
journal = {Physics of Plasmas},
number = 3,
volume = 24,
place = {United States},
year = {Mon Mar 20 00:00:00 EDT 2017},
month = {Mon Mar 20 00:00:00 EDT 2017}
}

Accurate Modeling of LaserPlasma Accelerators with ParticleInCell Codes
ParticleInCell (PIC) codes are often used to study systems where the details of phasespace are important; for example, self trapping or optical injection in laserplasma accelerators. Here we investigate the numerical heating and macroparticle trajectory errors in 2D PIC simulations of laserplasma accelerators. The effects of grid resolution and laser polarization on the momentum spread and on subsequent spurious trapping in a plasma wave is studied. It is shown that when the laser is polarized in the plane of the simulation, which mimics the 3D behavior, the macroparticles are subject to trajectory errors resulting in a high momentum spread. Themore » 
Modeling singlefrequency laserplasma acceleration using particleincell simulations: The physics of beam breakup
The authors investigate electron acceleration from spacecharge waves driven by singlefrequency lasers using a fully explicit particleincell (PIC) model. The two dimensional (2D) simulations model {approximately}100 fs pulses at densities near n = 4 {times} 10{sup 19} cm{sup {minus}3} for 1{micro}m lasers. The pulses are found to break up due to a combination of longitudinal and transverse bunching of the laser intensity via Raman forward scattering type instabilities. The ponderomotive force of these intensity modulations generates large amplitude plasma waves. Large numbers of selftrapped electrons and multiple Raman forward scattering satellites are observed. The relevance of these simulations to experimentsmore » 
Asynchronous communication in spectralelement and discontinuous Galerkin methods for atmospheric dynamics – a case study using the HighOrder Methods Modeling Environment (HOMMEhomme_dg_branch)
The scalability of computational applications on current and nextgeneration supercomputers is increasingly limited by the cost of interprocess communication. We implement nonblocking asynchronous communication in the HighOrder Methods Modeling Environment for the time integration of the hydrostatic fluid equations using both the spectralelement and discontinuous Galerkin methods. This allows the overlap of computation with communication, effectively hiding some of the costs of communication. A novel detail about our approach is that it provides some data movement to be performed during the asynchronous communication even in the absence of other computations. This method produces significant performance and scalability gains in largescalemore »Cited by 1