4 Search Results

Discretizing Maxwell's equations in Galilean (comoving) coordinates allows the derivation of a pseudospectral solver that eliminates the numerical Cherenkov instability for electromagnetic particle-in-cell simulations of relativistic plasmas flowing at a uniform velocity. Here we generalize this solver by incorporating spatial derivatives of arbitrary order, thereby enabling efficient parallelization by domain decomposition. This allows scaling of the algorithm to many distributed compute units. We derive the numerical dispersion relation of the algorithm and present a comprehensive theoretical stability analysis. The method is applied to simulations of plasma acceleration in a Lorentz-boosted frame of reference.

Intense X-ray free-electron lasers (XFELs) can rapidly excite matter, leaving it in inherently unstable states that decay on femtosecond timescales. The relaxation occurs primarily via Auger emission, so excited-state observations are constrained by Auger decay. In situ measurement of this process is therefore crucial, yet it has thus far remained elusive in XFELs owing to inherent timing and phase jitter, which can be orders of magnitude larger than the timescale of Auger decay. Here we develop an approach termed ‘self-referenced attosecond streaking’ that provides subfemtosecond resolution in spite of jitter, enabling time-domain measurement of the delay between photoemission and Augermore » emission in atomic neon excited by intense, femtosecond pulses from an XFEL. FUrthermore, using a fully quantum-mechanical description that treats the ionization, core-hole formation and Auger emission as a single process, the observed delay yields an Auger decay lifetime of \(2.2_{ - 0.3}^{ + 0.2}\) fs for the KLL decay channel.« less

Representing the electrodynamics of relativistically drifting particle ensembles in discrete, co-propagating Galilean coordinates enables the derivation of a Particle-In-Cell algorithm that is intrinsically free of the numerical Cherenkov instability for plasmas flowing at a uniform velocity. Application of the method is shown by modeling plasma accelerators in a Lorentz-transformed optimal frame of reference.

Particle in Cell (PIC) simulations are a widely used tool for the investigation of both laser- and beam-driven 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. Pseudo-spectral 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 pseudo-spectral methods provide this needed locality. Yet, these methods can again be pronemore » 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.« less