skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 10:00 PM ET on Friday, December 8 until 2:00 AM ET on Saturday, December 9 due to maintenance. We apologize for the inconvenience.

Title: Automation and control of laser wakefield accelerators using Bayesian optimization
Abstract Laser wakefield accelerators promise to revolutionize many areas of accelerator science. However, one of the greatest challenges to their widespread adoption is the difficulty in control and optimization of the accelerator outputs due to coupling between input parameters and the dynamic evolution of the accelerating structure. Here, we use machine learning techniques to automate a 100 MeV-scale accelerator, which optimized its outputs by simultaneously varying up to six parameters including the spectral and spatial phase of the laser and the plasma density and length. Most notably, the model built by the algorithm enabled optimization of the laser evolution that might otherwise have been missed in single-variable scans. Subtle tuning of the laser pulse shape caused an 80% increase in electron beam charge, despite the pulse length changing by just 1%.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
Nature Communications
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Ultrafast high-brightness X-ray pulses have proven invaluable for a broad range of research. Such pulses are typically generated via synchrotron emission from relativistic electron bunches using large-scale facilities. Recently, significantly more compact X-ray sources based on laser-wakefield accelerated (LWFA) electron beams have been demonstrated. In particular, laser-driven sources, where the radiation is generated by transverse oscillations of electrons within the plasma accelerator structure (so-called betatron oscillations) can generate highly-brilliant ultrashort X-ray pulses using a comparably simple setup. Here, we experimentally demonstrate a method to markedly enhance the parameters of LWFA-driven betatron X-ray emission in a proof-of-principle experiment. We show a significant increase in the number of generated photons by specifically manipulating the amplitude of the betatron oscillations by using our novel Transverse Oscillating Bubble Enhanced Betatron Radiation scheme. We realize this through an orchestrated evolution of the temporal laser pulse shape and the accelerating plasma structure. This leads to controlled off-axis injection of electrons that perform large-amplitude collective transverse betatron oscillations, resulting in increased radiation emission. Our concept holds the promise for a method to optimize the X-ray parameters for specific applications, such as time-resolved investigations with spatial and temporal atomic resolution or advanced high-resolution imaging modalities, and the generation of X-ray beams with even higher peak and average brightness.

    more » « less
  2. Abstract Temperature control is essential for regulating material properties in laser-based manufacturing. Motion and power of the scanning laser affect local temperature evolution, which in turn determines the a posteriori microstructure. This paper addresses the problem of adjusting the laser speed and power to achieve the desired values of key process parameters: cooling rate and melt pool size. The dynamics of a scanning laser system is modeled by a one-dimensional (1D) heat conduction equation, with laser power as the heat input and heat dissipation to the ambient. Since the model is 1D, length and size are essentially the same. We pose the problem as a regulation problem in the (moving) laser frame. The first step is to obtain the steady-state temperature distribution and the corresponding input based on the desired cooling rate and melt pool size. The controller adjusts the input around the steady-state feedforward based on the deviation of the measured temperature field from the steady-state distribution. We show that with suitably defined outputs, the system is strictly passive from the laser motion and power. To avoid over-reliance on the model, the steady-state laser speed and power are adaptively updated, resulting in an integral-like update law for the feedforward. Moreover, the heat transfer coefficient to the ambient may be uncertain, and can also be adaptively updated. The final form of the control law combines passive error temperature field feedback with adaptive feedforward and parameter estimation. The closed-loop asymptotical stability is shown using the Lyapunov arguments, and the controller performance is demonstrated in a simulation. 
    more » « less
  3. We present results from two new techniques for the generation of meter-scale, low density (∼1017 cm−3 on axis) plasma waveguides, the “two-Bessel” technique, and the “self-waveguiding” technique. Plasma waveguides of this density and length range are needed for demonstration of a ∼10 GeV laser wakefield accelerator module, key for future staging for a ∼TeV lepton collider. Both techniques require the use of high quality ultrashort pulse Bessel beams to efficiently and uniformly ionize hydrogen gas in meter-scale supersonic gas jets via optical field ionization. We review these two techniques, describe our meter-scale gas jets, and present a new method for correction of optical aberrations in Bessel beams. Finally, we briefly present results from recent experiments employing one of our techniques, demonstrating quasi-monoenergetic acceleration of ∼5 GeV electron bunches in 20 cm long, low density plasma waveguides. 
    more » « less
  4. Abstract Due to the highly nonlinear nature of the beam-loading, it is currently not possible to analytically determine the beam parameters needed in a two-bunch plasma wakefield accelerator for maintaining a low energy spread. Therefore in this paper, by using the Broyden–Fletcher–Goldfarb–Shanno algorithm for the parameter scanning with the code QuickPIC and the polynomial regression together with k -fold cross-validation method, we obtain two fitting formulas for calculating the parameters of tri-Gaussian electron beams when minimizing the energy spread based on the beam-loading effect in a nonlinear plasma wakefield accelerator. One formula allows the optimization of the normalized charge per unit length of a trailing beam to achieve the minimal energy spread, i.e. the optimal beam-loading. The other one directly gives the transformer ratio when the trailing beam achieves the optimal beam-loading. A simple scaling law for charges of drive beams and trailing beams is obtained from the fitting formula, which indicates that the optimal beam-loading is always achieved for a given charge ratio of the two beams when the length and separation of two beams and the plasma density are fixed. The formulas can also help obtain the optimal plasma densities for the maximum accelerated charge and the maximum acceleration efficiency under the optimal beam-loading respectively. These two fitting formulas will significantly enhance the efficiency for designing and optimizing a two-bunch plasma wakefield acceleration stage. 
    more » « less
  5. In this work, evolution of parameters of nanosecond repetitively pulsed (NRP) discharges in pin-to-pin configuration in air was studied during the transient stage of initial 20 discharge pulses. Gas and plasma parameters in the discharge gap were measured using coherent microwave scattering, optical emission spectroscopy, and laser Rayleigh scattering for NRP discharges at repetition frequencies of 1, 10, and 100 kHz. Memory effects (when perturbations induced by the previous discharge pulse would not decay fully until the subsequent pulse) were detected for the repetition frequencies of 10 and 100 kHz. For 10 kHz NRP discharge, the discharge parameters experienced significant change after the first pulse and continued to substantially fluctuate between subsequent pulses due to rapid evolution of gas density and temperature during the 100  μs inter-pulse time caused by intense redistribution of the flow field in the gap on that time scale. For 100 kHz NRP discharge, the discharge pulse parameters reached a new steady-state at about five pulses after initiation. This new steady-state was associated with well-reproducible parameters between the discharge pulses and substantial reduction in breakdown voltage, discharge pulse energy, and electron number density in comparison to the first discharge pulse. For repetition frequencies 1–100 kHz considered in this work, the memory effects can be likely attributed to the reduction in gas number density and increase in the gas temperature that cannot fully recover to ambient conditions before subsequent discharge pulses. 
    more » « less