Alexander Debus

PIConGPU: A Fully Relativistic Particle-in-Cell Code for a GPU Cluster

The particle-in-cell (PIC) algorithm is one of the most widely used algorithms in computational plasma physics. With the advent of graphical processing units (GPUs), large-scale plasma simulations on inexpensive GPU clusters are in reach. We present an implementation of a fully relativistic plasma PIC algorithm for GPUs based on the NVIDIA CUDA library. It supports a hybrid architecture consisting of single computation nodes interconnected in a standard cluster topology, with each node carrying one or more GPUs. The internode communication is realized using the message-passing interface. The simulation code PIConGPU presented in this paper is, to our knowledge, the first scalable GPU cluster implementation of the PIC algorithm in plasma physics.

Radiative signatures of the relativistic Kelvin-Helmholtz instability

We present a particle-in-cell simulation of the relativistic Kelvin-Helmholtz Instability (KHI) that for the first time delivers angularly resolved radiation spectra of the particle dynamics during the formation of the KHI. This enables studying the formation of the KHI with unprecedented spatial, angular and spectral resolution. Our results are of great importance for understanding astrophysical jet formation and comparable plasma phenomena by relating the particle motion observed in the KHI to its radiation signature. The innovative methods presented here on the implementation of the particle-in-cell algorithm on graphic processing units can be directly adapted to any many-core parallelization of the particle-mesh method. With these methods we see a peak performance of 7.176 PFLOP/s (double-precision) plus 1.449 PFLOP/s (single-precision), an efficiency of 96% when weakly scaling from 1 to 18432 nodes, an efficiency of 68.92% and a speed up of 794 (ideal: 1152) when strongly scaling from 16 to 18432 nodes.

Linear and non-linear Thomson-scattering x-ray sources driven by conventionally and laser plasma accelerated electrons

Compact tuneable sources of ultrashort hard x-ray pulses can be realized by Thomson scattering, taking advantage of the comparatively short wavelength of a scattered laser pulse with respect to the period length of conventional undulators. Here, we present a detailed analysis and optimization of the efficiency of linear and non-linear Thomson scattering when the process is driven with relativistic laser pulses and when the conventional accelerator is replaced by a laser-plasma wakefield accelerator.

Optical free-electron lasers with Traveling-Wave Thomson-Scattering

We present a fully analytic model of an all-optical free electron laser (OFEL) undulator based on the Traveling-Wave Thomson-Scattering (TWTS) scheme. The TWTS undulator provides for sub-mm undulator wavelengths, does not require any material or plasma to generate or contain the undulator field and allows for sub-meter saturation lengths. Starting from a fully analytic description of the three-dimensional TWTS field we derive the OFEL pendulum equation for electrons in the TWTS field and discuss the constraints on laser and electron pulse parameters that have to be fulfilled for OFEL operation. We conclude in applying the TWTS OFEL to the realization of compact free electron laser sources at 13.5 nm and 0.2 nm using laser and electron sources in reach of present day technologies.

First results with the novel petawatt laser acceleration facility in Dresden

We report on first commissioning results of the DRACO Petawatt ultra-short pulse laser system implemented at the ELBE center for high power radiation sources of Helmholtz-Zentrum Dresden-Rossendorf. Key parameters of the laser system essential for efficient and reproducible performance of plasma accelerators are presented and discussed with the demonstration of 40 MeV proton acceleration under TNSA conditions as well as peaked electron spectra with unprecedented bunch charge in the 0.5 nC range.

Demonstration of a beam loaded nanocoulomb-class laser wakefield accelerator

Laser-plasma wakefield accelerators have seen tremendous progress, now capable of producing quasi-monoenergetic electron beams in the GeV energy range with few-femtoseconds bunch duration. Scaling these accelerators to the nanocoulomb range would yield hundreds of kiloamperes peak current and stimulate the next generation of radiation sources covering high-field THz, high-brightness X-ray and γ-ray sources, compact free-electron lasers and laboratory-size beam-driven plasma accelerators. However, accelerators generating such currents operate in the beam loading regime where the accelerating field is strongly modified by the self-fields of the injected bunch, potentially deteriorating key beam parameters. Here we demonstrate that, if appropriately controlled, the beam loading effect can be employed to improve the accelerator’s performance. Self-truncated ionization injection enables loading of unprecedented charges of ∼0.5 nC within a mono-energetic peak. As the energy balance is reached, we show that the accelerator operates at the theoretically predicted optimal loading condition and the final energy spread is minimized.Higher beam quality and stability are desired in laser-plasma accelerators for their applications in compact light sources. Here the authors demonstrate in laser plasma wakefield electron acceleration that the beam loading effect can be employed to improve beam quality by controlling the beam charge.

Identifying the linear phase of the relativistic Kelvin-Helmholtz instability and measuring its growth rate via radiation

For the relativistic Kelvin-Helmholtz instability (KHI), which occurs at shear interfaces between two plasma streams, we report results on the polarized radiation over all observation directions and frequencies emitted by the plasma electrons from ab initio kinetic simulations. We find the polarization of the radiation to provide a clear signature for distinguishing the linear phase of the KHI from its other phases. During the linear phase, we predict the growth rate of the KHI radiation power to match the growth rate of the KHI to a high degree. Our predictions are based on a model of the vortex dynamics, which describes the electron motion in the vicinity of the shear interface between the two streams. Albeit the complex and turbulent dynamics happening in the shear region, we find excellent agreement between our model and large-scale particle-in-cell simulations. Our findings pave the way for identifying the KHI linear regime and for measuring its growth rate in astrophysical jets observable on earth as well as in laboratory plasmas.

Quantitatively consistent computation of coherent and incoherent radiation in particle-in-cell codes—A general form factor formalism for macro-particles

Abstract Quantitative predictions from synthetic radiation diagnostics often have to consider all accelerated particles. For particle-in-cell (PIC) codes, this not only means including all macro-particles but also taking into account the discrete electron distribution associated with them. This paper presents a general form factor formalism that allows to determine the radiation from this discrete electron distribution in order to compute the coherent and incoherent radiation self-consistently. Furthermore, we discuss a memory-efficient implementation that allows PIC simulations with billions of macro-particles. The impact on the radiation spectra is demonstrated on a large scale LWFA simulation.

Scaling EUV and X-ray Thomson sources to optical free-electron laser operation with traveling-wave Thomson scattering (Conference Presentation)

Traveling-Wave Thomson-Scattering (TWTS) allows for the realization of optical free-electron lasers (OFELs) from the interaction of short, high-power laser pulses with brilliant relativistic electron bunches. The laser field provides the optical undulator which is traversed by the electrons. In order to achieve coherent amplification of radiation through electron microbunching the interaction between electrons and laser must be maintained over hundreds to thousands of undulator periods. Traveling-Wave Thomson-Scattering is the only scattering geometry so far allowing for the realization of optical undulators of this length which is at the same time scalable from extreme ultraviolet to X-ray photon energies. TWTS is also applicable for the realization of incoherent high peak brightness hard X-ray to gamma-ray sources which can provide orders of magnitude higher photon output than classic head-on Thomson sources. In contrast to head-on Thomson sources TWTS employs a side-scattering geometry where laser and electron propagation direction of motion enclose an angle. Tilting the laser pulse front with respect to the wave front by half of this interaction angle optimizes electron and laser pulse overlap. In the side-scattering geometry the tilt of the pulse-front compensates the spatial offset between electrons and laser pulse-front which would be present otherwise for an electron bunch far from the interaction point where it overlaps with the laser pulse center. Thus the laser pulse-front tilt ensures continuous overlap between laser pulse and electrons while these traverse the laser pulse cross-sectional area. This allows to control the interaction distance in TWTS by the laser pulse width rather than laser pulse duration as is the case for head-on Thomson scattering. Utilizing petawatt class laser pulses with millimeter to centimeter scale width allows for the realization of compact optical undulators with thousands of periods. When laser pulses for TWTS are prepared, care has to be taken of laser dispersion. Especially for scenarios featuring interaction angles of several ten to over one hundred degree the angular dispersion originating from laser pulse-front tilt can significantly prolong the pulse duration during the interaction which leads to a decrease in optical undulator amplitude and eventually terminates the interaction long before the target interaction distance is reached. In the talk it is shown how a pair of two gratings can be used to first generate the pulse-front tilt and second control and compensate dispersion during the interaction by utilizing the plane of optimum compression. Furthermore an experimental setup strategy is presented allowing for an interaction outside the laser pulse focus. This is a necessity for TWTS OFELs requiring focusing to reach optical undulator strengths on the order of unity since the centimeter scale laser pulse width at the interaction point result in turn in Rayleigh lengths on the order of one hundred meter and thus in laser focusing distances of several hundred meter. The talk shows how an out-of-focus interaction geometry utilizing strong focusing of the incident laser pulse needs to be designed in order to regain compactness by reducing the focusing distance by one to two orders of magnitude.

Investigation of electron dynamics in an ionization-injection laser-wakefield accelerator via betatron radiation (Conference Presentation)

The injection process of electrons into the plasma cavity in laser-wakefield accelerators is a nonlinear process that strongly influences the property of the accelerated electrons. During the acceleration electrons perform transverse (betatron) oscillations around the axis. This results in the emission of hard x-ray radiation (betatron radiation) whose characteristics depend directly on the dynamic of the accelerated electrons. Thus, betatron radiation can be utilized as a powerful diagnostic tool to investigate the acceleration process inside the wakefield. Here we describe our recent LWFA experiments deploying ionization induced injection technique carried out with the Draco Ti:Sapphire laser. We focused 30 fs short pulses down to a FWHM spot size of 19 μm resulting in a normalized vacuum laser intensity a0 = 3.3 on a gas target. The target, which was a supersonic gas jet, provided a flat plasma profile of 3mm length. By varying the plasma density from 2x10^18 cm^-3 to 5x10^18 cm^-3 and the laser pulse energy from 1.6 J to 3.4 J we were able to tune the electron bunch and betatron parameters. Electron spectra were obtained by acquiring an energy resolved and charge calibrated electron profile after detection from the beam axis by a permanent magnetic dipole. Simultaneously, a back-illuminated and deep-depleted CCD placed on axis recorded the emitted x-ray photons with energies up to 20keV. Equipped with an 2D spectroscopy technique based on single pixel absorption events, we reconstructed the corresponding energy resolved x-ray spectrum for every shot and deduced the betatron source size at the plasma exit. Combining the data of the electron and betatron spectrum, we compare the characteristics of the betatron spectra for different electron bunches. In our experiments we recorded a total number of 25x10^4 photons per shot within a divergence angle of 1 mrad and betatron radii in the order of 1 μm. Finally, we compare our results with simulated spectra from the parallel classical radiation calculator Clara2 that is based on the Liénard-Wiechert potentials.