Search

Chapter 1 Automated detection and analysis of particle beams in laser-plasma accelerator simulations Daniela M. Ushizima 1 , Cameron G. Geddes 1 , Estelle Cormier-Michel 1 , E. Wes Bethel 1 , Janet Jacobsen 1 , Prabhat 1 , Oliver Ruebel 1,3 , Gunther Weber 1 , Peter Messmer 2 , Bernd Hamann 1 , Hans Haggen 4 1 Lawrence Berkeley National Laboratory, Berkeley, CA, USA 2 Tech-X Corporation, Boulder, CO, USA 3 University of Kaiserslautern, Germany Introduction Numerical simulations [1] of laser-plasma wakefield (particle) accelerators [2] model the acceler- ation of electrons trapped in plasma oscillations (wakes) left behind when an intense laser pulse propagates through the plasma. The goal of these simulations is to better understand the process involved in plasma wake generation and how electrons are trapped and accelerated by the wake [3, 4]. Understanding of such accelerators, and their development, offer high accelerating gradi- ents, potentially reducing size and cost of new accelerators. One operating regime of interest is where a trapped subset of electrons loads the wake and forms an isolated group of accelerated particles with low spread in momentum and position [3], desirable characteristics for many applications. The electrons trapped in the wake may be accel- erated to high energies, the plasma gradient in the wake reaching up to a gigaelectronvolt per cen- timeter [2]. High-energy electron accelerators power intense X-ray radiation to terahertz sources, and are used in many applications including medical radiotherapy and imaging [1]. To extract information from the simulation about the quality of the beam, a typical approach is to examine plots of the entire dataset, visually determining the adequate parameters necessary

Results are presented on the generation of quasimonoenergetic electron beams with energy up to 1GeV using a 40TW laser and a 3.3 cm-long hydrogen-filled capillary discharge waveguide [1, 2]. Electron beams were not observed without a plasma channel, indicating that self-focusing alone could not be relied upon for effective guiding of the laser pulse. Results are presented of the electron beam spectra, and the dependence of the reliability of producing electron beams as a function of laser and plasma parameters. ©2007 IEEE.

Simulation of laser-plasma accelerator (LPA) experiments is computationally intensive due to the disparate length scales involved. Current experiments extend hundreds of laser wavelengths transversely and many thousands in the propagation direction, making explicit PIC simulations enormously expensive and requiring massively parallel execution in 3D. Simulating the next generation of LPA experiments is expected to increase the computational requirements yet further, by a factor of 1000. We can substantially improve the performance of LPA simulations by modeling the envelope evolution of the laser field rather than the field itself. This allows for much coarser grids, since we need only resolve the plasma wavelength and not the laser wavelength, and therefore larger timesteps can be used. Thus an envelope model can result in savings of several orders of magnitude in computational resources. By propagating the laser envelope in a Galilean frame moving at the speed of light, dispersive errors can be avoided and simulations over long distances become possible. The primary limitation to this envelope model is when the laser pulse develops large frequency shifts, and thus the slowly-varying envelope assumption is no longer valid. Here we describe the model and its implementation, and show rigorous benchmarks for the algorithm, establishing second-order convergence and correct laser group velocity. We also demonstrate simulations of LPA phenomena such as self-focusing and meter-scale acceleration stages using the model.

A summary of progress at Lawrence Berkeley National Laboratory is given on: (1) experiments on down-ramp injection; (2) experiments on acceleration in capillary discharge plasma channels; and (3) simulations of a staged laser wakefield accelerator (LWFA). Control of trapping in a LWFA using plasma density down-ramps produced electron bunches with absolute longitudinal and transverse momentum spreads more than ten times lower than in previous experiments (0.17 and 0.02 MeV Ic FWHM, respectively) and with central momenta of 0.76 +- 0.02 MeV Ic, stable over a week of operation. Experiments were also carried out using a 40 TW laser interacting with a hydrogen-filled capillary discharge waveguide. For a 15 mm long, 200 mu m diameter capillary, quasi-monoenergetic bunches up to 300 MeV were observed. By detuning discharge delay from optimum guiding performance, self-trapping was found to be stabilized. For a 33 mm long, 300 mu m capillary, a parameter regime with high energy bunches, up to 1 Ge V, was found. In this regime, peak electron energy was correlated with the amount of trapped charge. Simulations show that bunches produced on a down-ramn and iniected into a channel-guided LWFA can produce stable beams with 0.2 MeV Ic-class momentum spread at high energies.

Plasma-based particle accelerators can produce and sustain thousands of times stronger acceleration fields than conventional particle accelerators, providing a potential solution to the problem of the growing size and cost of conventional particle accelerators. To facilitate scientific knowledge discovery from the ever growing collections of accelerator simulation data generated by accelerator physicists to investigate next-generation plasma-based particle accelerator designs, we describe a novel approach for automatic detection and classification of particle beams and beam substructures due to temporal differences in the acceleration process, here called acceleration features. The automatic feature detection in combination with a novel visualization tool for fast, intuitive, query-based exploration of acceleration features enables an effective top-down data exploration process, starting from a high-level, feature-based view down to the level of individual particles. We describe the application of our analysis in practice to analyze simulations of single pulse and dual and triple colliding pulse accelerator designs, and to study the formation and evolution of particle beams, to compare substructures of a beam, and to investigate transverse particle loss.

Knowledge discovery from large and complex scientific data is a challenging task. With the ability to measure and simulate more processes at increasingly finer spatial and temporal scales, the growing number of data dimensions and data objects presents tremendous challenges for effective data analysis and data exploration methods and tools. The combination and close integration of methods from scientific visualization, information visualization, automated data analysis, and other enabling technologies—such as efficient data management—supports knowledge discovery from multi-dimensional scientific data. This paper surveys two distinct applications in developmental biology and accelerator physics, illustrating the effectiveness of the described approach.

Because of their ultra-high accelerating gradient, laser plasma based accelerators (LPA) are contemplated for the next generation of high energy colliders and light sources. The upcoming BELLA project will explore acceleration of electron bunches to 10 GeV in a meter long plasma, where a wakefield is driven by a PW-class laser. Particle-in-cell (PIC) simulations are used to design the upcoming experiments. Simulations are challenging because of the disparity of length scale between the laser wavelength (~1 micron) that needs to be resolved and the simulation length (~ 1 m). We report on recent developments of the Laser Envelope Model, a reduced model for laser-plasma interactions that has previously demonstrated orders of magnitude speedup. In particular, we present the implementation of the model in cylindrical coordinates, allowing for quite rapid prototyping of laser acceleration stages. We discuss the performance benefits as well as the limitations and trade-offs of this model. In parallel, high frequency noise in PIC simulations makes it difficult to accurately represent beam energy spread and emittance. We show that calculating the beam self-fields using a static Poisson solve in the beam frame dramatically reduces particle noise, allowing for more accurate simulation of the beam evolution.

This final report describes the work that has been accomplished over the past 5 years under the Community Petascale Project for Accelerator and Simulations (ComPASS) at Tech-X Corporation. Tech-X had been involved in the full range of ComPASS activities with simulation of laser plasma accelerator concepts, mainly in collaboration with LOASIS program at LBNL, simulation of coherent electron cooling in collaboration with BNL, modeling of electron clouds in high intensity accelerators, in collaboration with researchers at Fermilab and accurate modeling of superconducting RF cavity in collaboration with Fermilab, JLab and Cockcroft Institute in the UK.

We describe a modification to the finite-difference time-domain algorithm for electromagnetics on a Cartesian grid which eliminates numerical dispersion error in vacuum for waves propagating along a grid axis. We provide details of the algorithm, which generalizes previous work by allowing 3D operation with a wide choice of aspect ratio, and give conditions to eliminate dispersive errors along one or more of the coordinate axes. We discuss the algorithm in the context of laser-plasma acceleration simulation, showing significant reduction---up to a factor of 280, at a plasma density of ${10}^{23}\text{ }\text{ }{\mathrm{m}}^{\ensuremath{-}3}$---of the dispersion error of a linear laser pulse in a plasma channel. We then compare the new algorithm with the standard electromagnetic update for laser-plasma accelerator stage simulations, demonstrating that by controlling numerical dispersion, the new algorithm allows more accurate simulation than is otherwise obtained. We also show that the algorithm can be used to overcome the critical but difficult challenge of consistent initialization of a relativistic particle beam and its fields in an accelerator simulation.

Because of their ultra-high accelerating gradient, laser plasma based accelerators (LPA) are contemplated for the next generation of high-energy colliders and light sources. The upcoming BELLA project will explore acceleration of electron bunches to 10 GeV in a 1 meter long plasma, where a wakefield is driven by a PW-class laser. Particle-in-cell (PIC) simulations are used to design LPA stages relevant to the upcoming experiments. Simulations in a Lorentz boosted frame are used to gain significant speed up and to simulate the 10 GeV stages at full scale parameters, otherwise impractical. As criteria on energy spread and beam emittance become more stringent, PIC simulations become more challenging as high frequency noise artificially increases those quantities. To reduce numerical noise, we consider using a Poisson solve to calculate the beam self-fields. This method allows correct cancelation of the beam transverse self-forces and prevents artificial emittance growth for a matched beam in a linear focusin...

Higher-order laser modes are analyzed as a method to control focusing forces and improve the electron bunch quality in laser-plasma accelerators. In the linear wake regime, the focusing force is proportional to the transverse gradient of the laser intensity, which can be shaped by a superposition of modes. In particular, the transverse wakefield can be arbitrarily small in a region about the axis by adjusting the laser modes. Plasma channel effects, which prohibit the formation of the controlled-focusing region, can be mitigating by introducing a delay between the modes. Modes with parallel polarization produce a beat interference in the laser intensity, which lead to deflecting forces. This can be avoided by using modes with orthogonal polarization, different frequencies, or short pulses that do not overlap. Particle-in-cell simulations are performed of a laser-plasma accelerator in the quasilinear regime driven by high-order modes. Simulations show that, by including the first-order mode, the matched radius of the electron bunch is substantially increased, which for fixed bunch density and emittance implies an increase in the beam charge.

A number of recent developments have led to simulations of an inverse cyclotron for cooling intense muon beams for neutrino factories and muon colliders. Such a device could potentially act as a novel beam cooling mechanism for muons, and it would be significantly smaller and cheaper than other cooling channel designs. Realistic designs are still being explored, but the first simulations of particle tracking in the inverse cyclotron, with accumulation in the cyclotron core, have been done with electrostatic simulations in the particle-in-cell code VORPAL. We present an overview of the muon inverse cyclotron concept and recent simulation results., Comment: 6 pages, 3 figures, submitted to the Proceedings for the 2010 Advanced Accelerator Concepts workshop, Annapolis, MD

Improved computational methods are essential to the diverse and rapidly developing field of advanced accelerator concepts. We present an overview of some computational algorithms for laser‐plasma concepts and high‐brightness photocathode electron sources. In particular, we discuss algorithms for reduced laser‐plasma models that can be orders of magnitude faster than their higher‐fidelity counterparts, as well as important on‐going efforts to include relevant additional physics that has been previously neglected. As an example of the former, we present 2D laser wakefield accelerator simulations in an optimal Lorentz frame, demonstrating >10 GeV energy gain of externally injected electrons over a 2 m interaction length, showing good agreement with predictions from scaled simulations and theory, with a speedup factor of ∼2,000 as compared to standard particle‐in‐cell.

Laser plasma accelerators are able to produce high quality electron beams from 1 MeV to 1 GeV. The next generation of plasma accelerator experiments will likely use a multi‐stage approach where a high quality electron bunch is first produced and then injected into an accelerating structure. In this paper we present scaled particle‐in‐cell simulations of a 10 GeV stage in the quasi‐linear regime. We show that physical parameters can be scaled to be able to perform these simulations at reasonable computational cost. Beam loading properties and electron bunch energy gain are calculated. A range of parameter regimes are studied to optimize the quality of the electron bunch at the output of the stage.

Compact high-energy linacs are important to applications including monochromatic gamma sources for nuclear material security applications. Recent laser wakefield accelerator experiments at LBNL demonstrated narrow energy spread beams, now with energies of up to 1 GeV in 3 cm using a plasma channel at low density. This demonstrates the production of GeV beams from devices much smaller than conventional linacs, and confirms the anticipated scaling of laser driven accelerators to GeV energies. Stable performance at 0.5 GeV was demonstrated. Experiments and simulations are in progress to control injection of particles into the wake and hence to improve beam quality and stability. Using plasma density gradients to control injection, stable beams at 1 MeV over days of operation, and with an order of magnitude lower absolute momentum spread than previously observed, have been demonstrated. New experiments are post-accelerating the beams from controlled injection experiments to increase beam quality and stability. Thomson scattering from such beams is being developed to provide collimated multi-MeV monoenergetic gamma sources for security applications from compact devices. Such sources can reduce dose to target and increase accuracy for applications including photofission and nuclear resonance fluorescence.

International audience; In laser wakefield accelerators, electron motion is driven by intense forces that depend on the plasma density. Transverse oscillations in the accelerated electron orbits produce betatron radiation. The electron motion and the resulting betatron radiation spectrum can therefore be controlled by shaping the plasma density along the orbit of the electrons. Here, a method based on the use of a plasma with a longitudinal density variation (density depression or step) is proposed to increase the transverse oscillation amplitude and the energy of the electrons accelerated in a wakefield cavity. For fixed laser parameters, by appropriately tailoring the plasma profile, the betatron radiation emitted by these electrons is significantly increased in both flux and energy.

Laser wakefield particle accelerators have shown the potential to generate electric fields thousands of times higher than those of conventional accelerators. The resulting extremely short particle acceleration distance could yield a potential new compact source of energetic electrons and radiation, with wide applications from medicine to physics. Physicists investigate laser-plasma internal dynamics by running particle-in-cell simulations; however, this generates a large dataset that requires time-consuming, manual inspection by experts in order to detect key features such as beam formation. This paper describes a framework to automate the data analysis and classification of simulation data. First, we propose a new method to identify locations with high density of particles in the space-time domain, based on maximum extremum point detection on the particle distribution. We analyze high density electron regions using a lifetime diagram by organizing and pruning the maximum extrema as nodes in a minimum spanning tree. Second, we partition the multivariate data using fuzzy clustering to detect time steps in a experiment that may contain a high quality electron beam. Finally, we combine results from fuzzy clustering and bunch lifetime analysis to estimate spatially confined beams. We demonstrate our algorithms successfully on four different simulation datasets.

Plasma density gradients in a gas jet were used to control the wake phase velocity and trapping threshold in a laser wakefield accelerator, producing stable electron bunches with longitudinal and transverse momentum spreads more than ten times lower than in previous experiments (0.17 and 0.02 MeV/c FWHM, respectively) and with central momenta of 0.76 +- 0.02 MeV/c. Transition radiation measurements combined with simulations indicated that the bunches can be used as a wakefield accelerator injector to produce stable beams with 0.2 MeV/c-class momentum spread at high energies.

Laser driven wakefield accelerators produce accelerating fields thousands of times those achievable in conventional radio-frequency accelerators, offering compactness and ultrafast bunches to potentially extend the frontiers of high energy physics and enable laboratory scale ultrafast radiation sources. Realization of this potential requires understanding of accelerator physics to advance beam performance and stability, and particle simulations model the highly nonlinear, kinetic physics required. One-to-one simulations of experiments provide new insight for optimization and development of 100 MeV to GeV and beyond laser accelerator stages, and on production of reproducible and controllable low energy spread beams with improved emittance (focusability) and energy through control of injection. © 2007 IOP Publishing Ltd.

Laser wakefield accelerators produce accelerating gradients up to hundreds of GeV/m, and recently demonstrated 1-10 MeV energy spread at energies up to 1 GeV using electrons self-trapped from the plasma. Controlled injection and staging may further improve beam quality by circumventing tradeoffs between energy, stability, and energy spread/emittance. We present experiments demonstrating production of a stable electron beam near 1 MeV with hundred-keV level energy spread and central energy stability by using the plasma density profile to control self- injection, and supporting simulations. Simulations indicate that such beams can be post accelerated to high energies, potentially reducing momentum spread in laser accelerators by 100-fold or more.

Numerical simulations of laser wakefield particle accelerators play a key role in the understanding of the complex acceleration process and in the design of expensive experimental facilities. As the size and complexity of simulation output grows, an increasingly acute challenge is the practical need for computational techniques that aid in scientific knowledge discovery. To that end, we present a set of data-understanding algorithms that work in concert in a pipeline fashion to automatically locate and analyze high-energy particle bunches undergoing acceleration in very large simulation datasets. These techniques work cooperatively by first identifying features of interest in individual timesteps, then integrating features across timesteps, and based on the information-derived perform analysis of temporally dynamic features. This combination of techniques supports accurate detection of particle beams enabling a deeper level of scientific understanding of physical phenomena than has been possible before. By combining efficient data analysis algorithms and state-of-the-art data management we enable high-performance analysis of extremely large particle datasets in 3D. We demonstrate the usefulness of our methods for a variety of 2D and 3D datasets and discuss the performance of our analysis pipeline.

It was shown recently that it may be computationally advantageous to perform computer simulations in a Lorentz boosted frame for a certain class of systems. However, even if the computer model relies on a covariant set of equations, it was pointed out that algorithmic difficulties related to discretization errors may have to be overcome in order to take full advantage of the potential speedup. In this paper, we summarize the findings, the difficulties and their solutions, and review the applications of the technique that have been performed to date.

The concept and designs of plasma-based advanced accelerators for high energy physics and photon science are modelled in the SciDAC COMPASS project with a suite of Particle-In-Cell codes and simulation techniques including the full electromagnetic model, the envelope model, the boosted frame approach and the quasi-static model. In this paper, we report the progress of the development of these models and techniques and present recent results achieved with large-scale parallel PIC simulations. The simulation needs for modelling the plasma-based advanced accelerator at the energy frontier is discussed and a path towards this goal is outlined.

One of the central challenges facing visualization research is how to effectively enable knowledge discovery. An effective approach will likely combine application architectures that are capable of running on today's largest platforms to address the challenges posed by large data with visual data analysis techniques that help find, represent, and effectively convey scientifically interesting features and phenomena.

As scientific instruments and computer simulations produce more and more data, the task of locating the essential information to gain insight becomes increasingly difficult. FastBit is an efficient software tool to address this challenge. In this article, we present a summary of the key underlying technologies, namely bitmap compression, encoding, and binning. Together these techniques enable FastBit to answer structured (SQL) queries orders of magnitude faster than popular database systems. To illustrate how FastBit is used in applications, we present three examples involving a high-energy physics experiment, a combustion simulation, and an accelerator simulation. In each case, FastBit significantly reduces the response time and enables interactive exploration on terabytes of data.

Laser- and particle beam-driven plasma wakefield accelerators produce accelerating fields thousands of times higher than radio-frequency accelerators, offering compactness and ultrafast bunches to extend the frontiers of high energy physics and to enable laboratory-scale radiation sources. Large-scale kinetic simulations provide essential understanding of accelerator physics to advance beam performance and stability and show and predict the physics behind recent demonstration of narrow energy spread bunches. Benchmarking between codes is establishing validity of the models used and, by testing new reduced models, is extending the reach of simulations to cover upcoming meter-scale multi-GeV experiments. This includes new models that exploit Lorentz boosted simulation frames to speed calculations. Simulations of experiments showed that recently demonstrated plasma gradient injection of electrons can be used as an injector to increase beam quality by orders of magnitude. Simulations are now also modeling accelerator stages of tens of GeV, staging of modules, and new positron sources to design next-generation experiments and to use in applications in high energy physics and light sources.

Finite plasma temperature can modify the structure of the wake field, reduce the wave-breaking field, and lead to self-trapped electrons, which can degrade the electron bunch quality in a plasma-based accelerator. A relativistic warm fluid theory is used to describe the plasma temperature evolution and alterations to the structure of a nonlinear periodic wave exited in a warm plasma. The trapping threshold for a plasma electron and the fraction of electrons trapped from a thermal distribution are examined using a single-particle model. Numerical artifacts in particle-in-cell models that can mimic the physics associated with finite momentum spread are discussed.

Modeling of laser-plasma wakefield accelerators in an optimal frame of reference [J.-L. Vay, Phys. Rev. Lett. 98, 130405 (2007)] allows direct and efficient full-scale modeling of deeply depleted and beam loaded laser-plasma stages of 10 GeV-1 TeV (parameters not computationally accessible otherwise). This verifies the scaling of plasma accelerators to very high energies and accurately models the laser evolution and the accelerated electron beam transverse dynamics and energy spread. Over 4, 5, and 6 orders of magnitude speedup is achieved for the modeling of 10 GeV, 100 GeV, and 1 TeV class stages, respectively. Agreement at the percentage level is demonstrated between simulations using different frames of reference for a 0.1 GeV class stage. Obtaining these speedups and levels of accuracy was permitted by solutions for handling data input (in particular, particle and laser beams injection) and output in a relativistically boosted frame of reference, as well as mitigation of a high-frequency instability that otherwise limits effectiveness.