Richard Pausch

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.

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.

Visualizing the Radiation of the Kelvin-Helmholtz Instability

Emerging new technologies in plasma simulations allow tracking billions of particles while computing their radiative spectra. We present a visualization of the relativistic Kelvin-Helmholtz instability from a simulation performed with the fully-relativistic particle-in-cell code PIConGPU powered by 18,000 GPUs on the USAs fastest supercomputer Titan.

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.

Simulate what is measured: next steps towards predictive simulations (Conference Presentation)

Simulations of laser matter interaction at extreme intensities that have predictive power are nowadays in reach when considering codes that make optimum use of high performance compute architectures. Nevertheless, this is mostly true for very specific settings where model parameters are very well known from experiment and the underlying plasma dynamics is governed by Maxwells equations solely. When including atomic effects, prepulse influences, radiation reaction and other physical phenomena things look different. Not only is it harder to evaluate the sensitivity of the simulation result on the variation of the various model parameters but numerical models are less well tested and their combination can lead to subtle side effects that influence the simulation outcome. We propose to make optimum use of future compute hardware to compute statistical and systematic errors rather than just find the mots optimum set of parameters fitting an experiment. This requires to include experimental uncertainties which is a challenge to current state of the art techniques. Moreover, it demands better comparison to experiments as inclusion of simulating the diagnostics response becomes important. We strongly advocate the use of open standards for finding interoperability between codes for comparison studies, building complete tool chains for simulating laser matter experiments from start to end.