Light-Matter Interaction near the Schwinger Limit Using Tightly Focused Doppler-Boosted Lasers.

Strong-field quantum electrodynamics (SF QED) is a burgeoning research topic dealing with electromagnetic fields comparable to the Schwinger field (≈1.32×10^{18} V/m). While most past and proposed experiments rely on reaching this field in the rest frame of relativistic particles, the Schwinger limit could also be approached in the laboratory frame by focusing to its diffraction limit the light reflected by a plasma mirror irradiated by a multipetawatt laser. We explore the interaction between such intense light and matter with particle-in-cell simulations. We find that the collision with a relativistic electron beam would enable the study of the nonperturbative regime of SF QED, while the interaction with a solid target leads to a profusion of SF QED effects that retroact on the interaction. In both cases, relativistic attosecond pair jets with high densities are formed.

Broadband spectral combining of three pulse-shaped fiber amplifiers with 42fs compressed pulse duration.

We demonstrate ultra-broadband spectral combining of ultrashort pulses from Yb-doped fiber amplifiers, with coherently spectrally synthesized pulse shaping, to achieve tens-of-fs pulses. This method can fully compensate for gain narrowing and high order dispersion over broad bandwidth. We produce 42fs pulses by spectrally synthesizing three chirped-pulse fiber amplifiers and two programmable pulse shapers across an 80nm overall bandwidth. To the best of our knowledge, this is the shortest pulse duration achieved from a spectrally combined fiber system at one-micron wavelength. This work provides a path toward high-energy, tens-of-fs fiber chirped-pulse amplification systems.

Exascale applications: skin in the game.

As noted in Wikipedia, skin in the game refers to having 'incurred risk by being involved in achieving a goal', where 'skin is a synecdoche for the person involved, and game is the metaphor for actions on the field of play under discussion'. For exascale applications under development in the US Department of Energy Exascale Computing Project, nothing could be more apt, with the skin being exascale applications and the game being delivering comprehensive science-based computational applications that effectively exploit exascale high-performance computing technologies to provide breakthrough modelling and simulation and data science solutions. These solutions will yield high-confidence insights and answers to the most critical problems and challenges for the USA in scientific discovery, national security, energy assurance, economic competitiveness and advanced healthcare. This article is part of a discussion meeting issue 'Numerical algorithms for high-performance computational science'.

Absorption of charged particles in perfectly matched layers by optimal damping of the deposited current.

Perfectly matched layers (PMLs) are widely used in particle-in-cell simulations, in order to absorb electromagnetic waves that propagate out of the simulation domain. However, when charged particles cross the interface between the simulation domain and the PMLs, a number of numerical artifacts can arise. In order to mitigate these artifacts, we introduce a PML algorithm whereby the current deposited by the macroparticles in the PML is damped by an analytically derived optimal coefficient. The benefits of this algorithm are illustrated in practical simulations. In particular, it is shown that this algorithm is well suited for particles exiting the box in near-normal incidence, in the sense that the fields behave as if the exiting particle is propagating in an infinite vacuum.

Overcoming timestep limitations in boosted-frame particle-in-cell simulations of plasma-based acceleration.

Explicit electromagnetic particle-in-cell (PIC) codes are typically limited by the Courant-Friedrichs-Lewy (CFL) condition, which implies that the timestep multiplied by the speed of light must be smaller than the smallest cell size. In the case of boosted-frame PIC simulations of plasma-based acceleration, this limitation can be a major hindrance, as the cells are often very elongated along the longitudinal direction and the timestep is thus limited by the small, transverse cell size. This entails many small-timestep PIC iterations and can limit the potential speed-up of the boosted-frame technique. Here, by using a CFL-free analytical spectral solver, and by mitigating additional numerical instabilities that arise at large timestep, we show that it is possible to overcome traditional limitations on the timestep and thereby realize the full potential of the boosted-frame technique over a much wider range of parameters.

Scalable spectral solver in Galilean coordinates for eliminating the numerical Cherenkov instability in particle-in-cell simulations of streaming plasmas.

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.

Laser and electron deflection from transverse asymmetries in laser-plasma accelerators.

We report on the deflection of laser pulses and accelerated electrons in a laser-plasma accelerator (LPA) by the effects of laser pulse front tilt and transverse density gradients. Asymmetry in the plasma index of refraction leads to laser steering, which can be due to a density gradient or spatiotemporal coupling of the laser pulse. The transverse forces from the skewed plasma wave can also lead to electron deflection relative to the laser. Quantitative models are proposed for both the laser and electron steering, which are confirmed by particle-in-cell simulations. Experiments with the BELLA Petawatt Laser are presented which show controllable 0.1-1 mrad laser and electron beam deflection from laser pulse front tilt. This has potential applications for electron beam pointing control, which is of paramount importance for LPA applications.

Elimination of numerical Cherenkov instability in flowing-plasma particle-in-cell simulations by using Galilean coordinates.

Particle-in-cell (PIC) simulations of relativistic flowing plasmas are of key interest to several fields of physics (including, e.g., laser-wakefield acceleration, when viewed in a Lorentz-boosted frame) but remain sometimes infeasible due to the well-known numerical Cherenkov instability (NCI). In this article, we show that, for a plasma drifting at a uniform relativistic velocity, the NCI can be eliminated by simply integrating the PIC equations in Galilean coordinates that follow the plasma (also sometimes known as comoving coordinates) within a spectral analytical framework. The elimination of the NCI is verified empirically and confirmed by a theoretical analysis of the instability. Moreover, it is shown that this method is applicable both to Cartesian geometry and to cylindrical geometry with azimuthal Fourier decomposition.

WarpIV: In Situ Visualization and Analysis of Ion Accelerator Simulations.

The generation of short pulses of ion beams through the interaction of an intense laser with a plasma sheath offers the possibility of compact and cheaper ion sources for many applications--from fast ignition and radiography of dense targets to hadron therapy and injection into conventional accelerators. To enable the efficient analysis of large-scale, high-fidelity particle accelerator simulations using the Warp simulation suite, the authors introduce the Warp In situ Visualization Toolkit (WarpIV). WarpIV integrates state-of-the-art in situ visualization and analysis using VisIt with Warp, supports management and control of complex in situ visualization and analysis workflows, and implements integrated analytics to facilitate query- and feature-based data analytics and efficient large-scale data analysis. WarpIV enables for the first time distributed parallel, in situ visualization of the full simulation data using high-performance compute resources as the data is being generated by Warp. The authors describe the application of WarpIV to study and compare large 2D and 3D ion accelerator simulations, demonstrating significant differences in the acceleration process in 2D and 3D simulations. WarpIV is available to the public via https://bitbucket.org/berkeleylab/warpiv. The Warp In situ Visualization Toolkit (WarpIV) supports large-scale, parallel, in situ visualization and analysis and facilitates query- and feature-based analytics, enabling for the first time high-performance analysis of large-scale, high-fidelity particle accelerator simulations while the data is being generated by the Warp simulation suite. This supplemental material https://extras.computer.org/extra/mcg2016030022s1.pdf provides more details regarding the memory profiling and optimization and the Yee grid recentering optimization results discussed in the main article.

Absolute measurement of electron-cloud density in a positively charged particle beam.

Clouds of stray electrons are ubiquitous in particle accelerators and frequently limit the performance of storage rings. Earlier measurements of electron energy distribution and flux to the walls provided only a relative electron-cloud density. We have measured electron accumulation using ions expelled by the beam. The ion energy distribution maps the depressed beam potential and gives the dynamic cloud density. Clearing electrode current reveals the static background cloud density, allowing the first absolute measurement of the time-dependent electron-cloud density during the beam pulse.