Search

Particle-in-cell simulations are among the most essential tools for the modeling and optimization of laser-plasma accelerators, since they reproduce the physics from first principles. However, the high computational cost associated with them can severely limit the scope of parameter and design optimization studies. Here, we show that a multitask Bayesian optimization algorithm can be used to mitigate the need for such high-fidelity simulations by incorporating information from inexpensive evaluations of reduced physical models. In a proof-of-principle study, where a high-fidelity optimization with fbpic is assisted by reduced-model simulations with wake-t, the algorithm demonstrates an order-of-magnitude speedup. This opens a path for the cost-effective optimization of laser-plasma accelerators in large parameter spaces, an important step toward fulfilling the high beam quality requirements of future applications.

Current high-power laser amplifiers use chirped-pulse amplification to prevent damage to their solid-state components caused by intense electromagnetic fields. To increase laser power further requires ever larger and more expensive devices. The Raman backscatter instability in plasma facilitates an alternative amplification strategy without the limitations imposed by material damage thresholds. Plasma-based amplification has been experimentally demonstrated, but only with relatively low efficiency. Further progress requires extensive use of numerical simulations, which usually need significant computational resources. Here we present particle-in-cell (PIC) simulation techniques for accurately simulating Raman amplification using a moving window with suitable boundary conditions, reducing computational cost. We show that an analytical model for matched pump propagation in a parabolic plasma channel slightly overestimates amplification as pump laser intensity is increased. However, a method for loading data saved from separate pump-only simulations demonstrates excellent agreement with full PIC simulation. The reduction in required resources will enable parameter scans to be performed to optimize amplification, and stimulate efforts toward developing viable plasma-based laser amplifiers. The methods may also be extended to investigate Brillouin scattering, and for the development of laser wakefield accelerators. Efficient, compact, low-cost amplifiers would have widespread applications in academia and industry.

Similarly to laser or x-ray beams, the interaction of sufficiently intense particle beams with neutral gases will result in the creation of plasma. In contrast to photon-based ionization, the strong unipolar field of a particle beam can generate a plasma where the electron population receives a large initial momentum kick and escapes, leaving behind unshielded ions. Measuring the properties of the ensuing Coulomb exploding ions—such as their kinetic energy distribution, yield, and spatial distribution—can provide information about the peak electric fields that are achieved in the electron beams. Particle-in-cell simulations and analytical models are presented for high-brightness electron beams of a few femtoseconds or even hundreds of attoseconds, and transverse beam sizes on the micron scale, as generated by today’s free electron lasers. Different density regimes for the utilization as a potential diagnostics are explored, and the fundamental differences in plasma dynamical behavior for e-beam or photon-based ionization are highlighted. By measuring the dynamics of field-induced ions for different gas and beam densities, a lower bound on the beam charge density can be obtained in a single shot and in a noninvasive way. The exponential dependency of the ionization yield on the beam properties can provide unprecedented spatial and temporal resolution, at the submicrometer and subfemtosecond scales, respectively, offering a practical and powerful approach to characterizing beams from accelerators at the frontiers of performance.

Simulations of ionization-induced injection in a laser driven plasma wakefield show that high-quality electron injectors in the 50–200 MeV range can be achieved in a gas cell with a tailored density profile. Using the PIC code Warp with parameters close to existing experimental conditions, we show that the concentration of N_{2} in a hydrogen plasma with a tailored density profile is an efficient parameter to tune electron beam properties through the control of the interplay between beam loading effects and varying accelerating field in the density profile. For a given laser plasma configuration, with moderate normalized laser amplitude, a_{0}=1.6 and maximum electron plasma density, n_{e0}=4×10^{18} cm^{-3}, the optimum concentration results in a robust configuration to generate electrons at 150 MeV with a rms energy spread of 4% and a spectral charge density of 1.8 pC/MeV.

Control of the properties of laser-plasma-accelerated electron beams that were injected along a shock-induced density downramp through precision tailoring of the density profile was demonstrated using a 1.8 J, 45 fs laser interacting with a mm-scale gas jet. The effects on the beam spatial profile, steering, and absolute energy spread of the density region before the shock and tilt of the shock were investigated experimentally and with particle-in-cell simulations. By adjusting these density parameters, the electron beam quality was controlled and improved while the energy (30–180 MeV) and energy spread (2–11 MeV) were independently tuned. Simple models that are in good agreement with the experimental results are proposed to explain these relationships, advancing the understanding of downramp injection. This technique allows for high-quality electron beams with percent-level energy spread to be tailored based on the application.

The dynamics of electron acceleration driven by laser wakefield is studied in detail using the particle-in-cell code WARP with the objective to generate high-quality electron bunches with narrow energy spread and small emittance, relevant for the electron injector of a multistage accelerator. Simulation results, using experimentally achievable parameters, show that electron bunches with an energy spread of ∼11% can be obtained by using an ionization-induced injection mechanism in a mm-scale length plasma. By controlling the focusing of a moderate laser power and tailoring the longitudinal plasma density profile, the electron injection beginning and end positions can be adjusted, while the electron energy can be finely tuned in the last acceleration section.

Thanks to their compactness and unique properties, laser-wakefield accelerators are currently considered for several innovative applications. However, many of these applications—and especially those that require beam transport—are hindered by the large divergence of laser-accelerated beams. Here we propose a collimating concept that relies on the strong radial electric field of the laser-wakefield to reduce this divergence. This concept utilizes an additional gas jet, placed after the laser-wakefield accelerator. When the laser pulse propagates through this additional gas jet, it drives a wakefield which can refocus the trailing electron bunch. Particle-in-cell simulations demonstrate that this approach can reduce the divergence by at least a factor of 3 for realistic electron bunches.

Transverse emittance is a crucial feature of laser-wakefield accelerators, yet accurately reproducing its value in numerical simulations remains challenging. It is shown here that, when the charge of the bunch exceeds a few tens of picocoulombs, particle-in-cell (PIC) simulations erroneously overestimate the emittance. This is mostly due the interaction of spurious Cherenkov radiation with the bunch, which leads to a steady growth of emittance during the simulation. A new computational scheme is proposed, which is free of spurious Cherenkov radiation. It can be easily implemented in existing PIC codes and leads to a substantial reduction of the emittance growth.

Laser plasma accelerators (LPA) offer an exciting possibility to deliver high energy, high brightness electrons beams in drastically smaller distance scales than is typical for conventional accelerators. As such, LPAs draw considerable attention as potential drivers for next generation light sources and for a compact linear collider. In order to asses the viability of an LPA source for a particular application, the brightness of the source should be properly characterized. In this paper, we present charge dependent transverse emittance measurements of LPA sources using both ionization injection and shock induced density down ramp injection, with the latter delivering smaller transverse emittances by a factor of two when controlling for charge density. The single shot emittance method is described in detail with a discussion on limitations related to second order transport effects. The direct role of space charge is explored through a series of simulations and found to be consistent with experimental observations. [ABSTRACT FROM AUTHOR]

Propagating a relativistically intense, negatively chirped laser pulse (the bandwidth  >150 nm) in a plasma channel makes it possible to generate background-free, comb-like electron beams—sequences of synchronized bunches with a low phase-space volume and controlled energy spacing. The tail of the pulse, confined in the accelerator cavity (an electron density ‘bubble’), experiences periodic focusing, while the head, which is the most intense portion of the pulse, steadily self-guides. Oscillations of the cavity size cause periodic injection of electrons from the ambient plasma, creating an electron energy comb with the number of components, their mean energy, and energy spacing dependent on the channel radius and pulse length. These customizable electron beams enable the design of a tunable, all-optical source of pulsed, polychromatic γ-rays using the mechanism of inverse Thomson scattering, with up to  ∼10−5 conversion efficiency from the drive pulse in the electron accelerator to the γ-ray beam. Such a source may radiate  ∼107 quasi-monochromatic photons per shot into a microsteradian-scale cone. The photon energy is distributed among several distinct bands, each having sub-30% energy spread, with a highest energy of 12.5 MeV. [ABSTRACT FROM AUTHOR]

The recent advances in developing compact laser plasma accelerators that deliver high quality electron beams in a more reliable way offer the possibility to consider their use in designing a compact free electron laser (FEL). Because of the particularity of these beams (especially concerning the divergence and the energy spread), specific electron beam handling is proposed in order to achieve FEL amplification. [ABSTRACT FROM AUTHOR]

This paper introduces a formulation of the particle-in-cell (PIC) method for the modeling of relativistic plasmas, that leverages the ability of the pseudospectral analytical time-domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle (applied to second-order polynomial dependencies here). The formulation is applied to a modified set of Maxwell's equations that was proposed earlier in the context of divergence cleaning, and to recently proposed extensions of the PSATD-PIC algorithm. Detailed analysis and testings revealed that, under some condition, the formulation can expand the range of numerical parameters under which PIC simulations are stable and accurate when modeling relativistic plasmas such as, e.g., plasma-based particle accelerators.

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.

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.

Generating high-quality laser-plasma accelerated electron beams requires carefully balancing a plethora of physical effects and is therefore challenging-both conceptually and in experiments. Here, we use Bayesian optimization of key laser and plasma parameters to flatten the longitudinal phase space of an ionization-injected electron bunch via optimal beam loading. We first study the concept with particle-in-cell simulations and then demonstrate it in experiments. Starting from an arbitrary set point, the plasma accelerator autonomously tunes the beam energy spread to the subpercent level at 254 MeV and 4.7  pC/MeV spectral density. Finally, we study a robust regime, which improves the stability of the laser-plasma accelerator and delivers sub-five-percent rms energy spread beams for 90% of all shots.

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.

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.

The beam hosing instability is analyzed theoretically for a witness beam in the quasilinear regime of plasma accelerators. In this regime, the hosing instability saturates, even for a monoenergetic bunch, at a level much less than standard scalings predict. Analytic expressions are derived for the saturation distance and amplitude and are in agreement with numerical results. Saturation is due to the natural head-to-tail variations in the focusing force, including the self-consistent transverse beam loading.

Single-shot, charge-dependent emittance measurements of electron beams generated by a laser plasma accelerator (LPA) reveal that shock-induced density down-ramp injection produces beams with normalized emittances a factor of 2 smaller than beams produced via ionization injection. Such a comparison is made possible by the tunable LPA setup, which allows electron beams with nearly identical central energy and peak spectral charge density to be produced using the two distinct injection mechanisms. Parametric measurements of this type are essential for the development of LPA-based applications which ultimately require high charge density and low emittance.

With the advent of petawatt class lasers, the very large laser intensities attainable on target should enable the production of intense high-order Doppler harmonics from relativistic laser-plasma mirror interactions. At present, the modeling of these harmonics with particle-in-cell (PIC) codes is extremely challenging as it implies an accurate description of tens to hundreds of harmonic orders on a broad range of angles. In particular, we show here that due to the numerical dispersion of waves they induce in vacuum, standard finite difference time domain (FDTD) Maxwell solvers employed in most PIC codes can induce a spurious angular deviation of harmonic beams potentially degrading simulation results. This effect was extensively studied and a simple toy model based on the Snell-Descartes law was developed that allows us to finely predict the angular deviation of harmonics depending on the spatiotemporal resolution and the Maxwell solver used in the simulations. Our model demonstrates that the mitigation of this numerical artifact with FDTD solvers mandates very high spatiotemporal resolution preventing realistic three-dimensional (3D) simulations even on the largest computers available at the time of writing. We finally show that nondispersive pseudospectral analytical time domain solvers can considerably reduce the spatiotemporal resolution required to mitigate this spurious deviation and should enable in the near future 3D accurate modeling on supercomputers in a realistic time to solution.

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.

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.

Laser-plasma technology promises a drastic reduction of the size of high-energy electron accelerators. It could make free-electron lasers available to a broad scientific community and push further the limits of electron accelerators for high-energy physics. Furthermore, the unique femtosecond nature of the source makes it a promising tool for the study of ultrafast phenomena. However, applications are hindered by the lack of suitable lens to transport this kind of high-current electron beams mainly due to their divergence. Here we show that this issue can be solved by using a laser-plasma lens in which the field gradients are five order of magnitude larger than in conventional optics. We demonstrate a reduction of the divergence by nearly a factor of three, which should allow for an efficient coupling of the beam with a conventional beam transport line.

The capability of plasmas to sustain ultrahigh electric fields has attracted considerable interest over the last decades and has given rise to laser-plasma engineering. Today, plasmas are commonly used for accelerating and collimating relativistic electrons, or to manipulate intense laser pulses. Here we propose an ultracompact plasma undulator that combines plasma technology and nanoengineering. When coupled with a laser-plasma accelerator, this undulator constitutes a millimetre-sized synchrotron radiation source of X-rays. The undulator consists of an array of nanowires, which are ionized by the laser pulse exiting from the accelerator. The strong charge-separation field, arising around the wires, efficiently wiggles the laser-accelerated electrons. We demonstrate that this system can produce bright, collimated and tunable beams of photons with 10-100 keV energies. This concept opens a path towards a new generation of compact synchrotron sources based on nanostructured plasmas.

The transverse properties of an electron beam are characterized by two quantities, the emittance which indicates the electron beam extent in the phase space and the angular momentum which allows for nonplanar electron trajectories. Whereas the emittance of electron beams produced in a laser-plasma accelerator has been measured in several experiments, their angular momentum has been scarcely studied. It was demonstrated that electrons in a laser-plasma accelerator carry some angular momentum, but its origin was not established. Here we identify one source of angular-momentum growth and we present experimental results showing that the angular-momentum content evolves during the acceleration.

Laser-wakefield acceleration constitutes a promising technology for future electron accelerators. A crucial step in such an accelerator is the injection of electrons into the wakefield, which will largely determine the properties of the extracted beam. We present here a new paradigm of colliding-pulse injection, which allows us to generate high-quality electron bunches having both a very low emittance (0.17 mm·mrad) and a low energy spread (2%), while retaining a high charge (~100 pC) and a short duration (3 fs). In this paradigm, the pulse collision provokes a transient expansion of the accelerating bubble, which then leads to transverse electron injection. This mechanism contrasts with previously observed optical injection mechanisms, which were essentially longitudinal. We also specify the range of parameters in which this new type of injection occurs and show that it is within reach of existing high-intensity laser facilities.

Laser-plasma accelerators can produce high-quality electron beams, up to giga electronvolts in energy, from a centimetre scale device. The properties of the electron beams and the accelerator stability are largely determined by the injection stage of electrons into the accelerator. The simplest mechanism of injection is self-injection, in which the wakefield is strong enough to trap cold plasma electrons into the laser wake. The main drawback of this method is its lack of shot-to-shot stability. Here we present experimental and numerical results that demonstrate the existence of two different self-injection mechanisms. Transverse self-injection is shown to lead to low stability and poor-quality electron beams, because of a strong dependence on the intensity profile of the laser pulse. In contrast, longitudinal injection, which is unambiguously observed for the first time, is shown to lead to much more stable acceleration and higher-quality electron beams.

Many theoretical and experimental studies suggest that range expansions can have severe consequences for the gene pool of the expanding population. Due to strongly enhanced genetic drift at the advancing frontier, neutral and weakly deleterious mutations can reach large frequencies in the newly colonized regions, as if they were surfing the front of the range expansion. These findings raise the question of how frequently beneficial mutations successfully surf at shifting range margins, thereby promoting adaptation towards a range-expansion phenotype. Here, we use individual-based simulations to study the surfing statistics of recurrent beneficial mutations on wave-like range expansions in linear habitats. We show that the rate of surfing depends on two strongly antagonistic factors, the probability of surfing given the spatial location of a novel mutation and the rate of occurrence of mutations at that location. The surfing probability strongly increases towards the tip of the wave. Novel mutations are unlikely to surf unless they enjoy a spatial head start compared to the bulk of the population. The needed head start is shown to be proportional to the inverse fitness of the mutant type, and only weakly dependent on the carrying capacity. The precise location dependence of surfing probabilities is derived from the non-extinction probability of a branching process within a moving field of growth rates. The second factor is the mutation occurrence which strongly decreases towards the tip of the wave. Thus, most successful mutations arise at an intermediate position in the front of the wave. We present an analytic theory for the tradeoff between these factors that allows to predict how frequently substitutions by beneficial mutations occur at invasion fronts. We find that small amounts of genetic drift increase the fixation rate of beneficial mutations at the advancing front, and thus could be important for adaptation during species invasions.