David L. Bruhwiler

High-quality electron beams from a laser wakefield accelerator using plasma-channel guiding

Laser-driven accelerators, in which particles are accelerated by the electric field of a plasma wave (the wakefield) driven by an intense laser, have demonstrated accelerating electric fields of hundreds of GV m-1 (refs 1–3). These fields are thousands of times greater than those achievable in conventional radio-frequency accelerators, spurring interest in laser accelerators as compact next-generation sources of energetic electrons and radiation. To date, however, acceleration distances have been severely limited by the lack of a controllable method for extending the propagation distance of the focused laser pulse. The ensuing short acceleration distance results in low-energy beams with 100 per cent electron energy spread, which limits potential applications. Here we demonstrate a laser accelerator that produces electron beams with an energy spread of a few per cent, low emittance and increased energy (more than 109 electrons above 80 MeV). Our technique involves the use of a preformed plasma density channel to guide a relativistically intense laser, resulting in a longer propagation distance. The results open the way for compact and tunable high-brightness sources of electrons and radiation.

Production of high-quality electron bunches by dephasing and beam loading in channeled and unchanneled laser plasma accelerators

High-quality electron beams, with a few 109 electrons within a few percent of the same energy above 80 MeV, were produced in a laser wakefield accelerator by matching the acceleration length to the length over which electrons were accelerated and outran (dephased from) the wake. A plasma channel guided the drive laser over long distances, resulting in production of the high-energy, high-quality beams. Unchanneled experiments varying the length of the target plasma indicated that the high-quality bunches are produced near the dephasing length and demonstrated that channel guiding was more stable and efficient than relativistic self-guiding. Consistent with these data, particle-in-cell simulations indicate production of high-quality electron beams when trapping of an initial bunch of electrons suppresses further injection by loading the wake. The injected electron bunch is then compressed in energy by dephasing, when the front of the bunch begins to decelerate while the tail is still accelerated.

Particle-in-cell simulations of tunneling ionization effects in plasma-based accelerators

Plasma-based accelerators can sustain accelerating gradients on the order of 100 GV/m. If the plasma is not fully ionized, fields of this magnitude will ionize neutral atoms via electron tunneling, which can completely change the dynamics of the plasma wake. Particle-in-cell simulations of a high-field plasma wakefield accelerator, using the OOPIC code [D. L. Bruhwiler et al., Phys. Rev. ST Accel. Beams 4, 101302 (2001)], which includes field-induced tunneling ionization of neutral Li gas, show that the presence of even moderate neutral gas density significantly degrades the quality of the wakefield. The tunneling ionization model in OOPIC has been validated via a detailed comparison with experimental data from the l’OASIS laboratory [W.P. Leemans et al., Phys. Rev. Lett. 89, 174802 (2002)]. The properties of a wake generated directly from a neutral gas are studied, showing that one can recover the peak fields of the fully ionized plasma simulations, if the density of the electron drive bunch is increased...

Diffusion of particles in a slowly modulated wave

Abstract Particle motion in the trapping/detrapping regime of a slowly modulated wave is studied. It is shown that separatrix crossing causes diffusion of the adiabatic invariant J . As a result, the trajectories are ergodic throughout separatrix-swept phase space. A lowest-order diffusion coefficient, D 0 ( J ), is calculated by neglecting correlations between separatrix crossings. This theory predicts that the characteristic diffusion time, τ, should scale as the inverse cube of the modulation frequency. Numerical simulations confirm this scaling of τ, but indicate that correlations significantly affect the rate of diffusion.

Energy spectrum of particles accelerated near a magnetic x line

We study the acceleration of test particles near a static magnetic x line with a uniform electric field and a strong component of the magnetic field, B ∥, parallel to the x line. The energy spectrum of the accelerated particles is found analytically in the nonrelativistic limit, showing good agreement with numerical simulations. At high energies, the spectrum decays exponentially with a characteristic energy very different from that found in studies assuming B ∥=0

Numerical modeling of laser tunneling ionization in explicit particle-in-cell codes

Methods for the calculation of laser tunneling ionization in explicit particle-in-cell codes used for modeling laser-plasma interactions are compared and validated against theoretical predictions. Improved accuracy is obtained by using the direct current form for the ionization rate. Multi level ionization in a single time step and energy conservation have been considered during the ionization process. The effects of grid resolution and number of macro-particles per cell are examined. Implementation of the ionization algorithm in two different particle-in-cell codes is compared for the case of ionization-based electron injection in a laser-plasma accelerator.

Computationally efficient methods for modelling laser wakefield acceleration in the blowout regime

Electron self-injection and acceleration until dephasing in the blowout regime is studied for a set of initial conditions typical of recent experiments with 100-terawatt-class lasers. Two different approaches to computationally efficient, fully explicit, 3D particle-in-cell modelling are examined. First, the Cartesian code vorpal (Nieter, C. and Cary, J. R. 2004 VORPAL: a versatile plasma simulation code. J. Comput. Phys. 196, 538) using a perfect-dispersion electromagnetic solver precisely describes the laser pulse and bubble dynamics, taking advantage of coarser resolution in the propagation direction, with a proportionally larger time step. Using third-order splines for macroparticles helps suppress the sampling noise while keeping the usage of computational resources modest. The second way to reduce the simulation load is using reduced-geometry codes. In our case, the quasi-cylindrical code calder-circ (Lifschitz, A. F. et al. 2009 Particle-in-cell modelling of laser-plasma interaction using Fourier decomposition. J. Comput. Phys. 228(5), 1803-1814) uses decomposition of fields and currents into a set of poloidal modes, while the macroparticles move in the Cartesian 3D space. Cylindrical symmetry of the interaction allows using just two modes, reducing the computational load to roughly that of a planar Cartesian simulation while preserving the 3D nature of the interaction. This significant economy of resources allows using fine resolution in the direction of propagation and a small time step, making numerical dispersion vanishingly small, together with a large number of particles per cell, enabling good particle statistics. Quantitative agreement of two simulations indicates that these are free of numerical artefacts. Both approaches thus retrieve the physically correct evolution of the plasma bubble, recovering the intrinsic connection of electron self-injection to the nonlinear optical evolution of the driver.

Electron Cooling of RHIC

We report progress on the R&D program for electron-cooling of the Relativistic Heavy Ion Collider (RHIC). This electron cooler is designed to cool 100 GeV/nucleon at storage energy using 54 MeV electrons. The electron source will be a superconducting RF photocathode gun. The accelerator will be a superconducting energy recovery linac. The frequency of the accelerator is set at 703.75 MHz. The maximum electron bunch frequency is 9.38 MHz, with bunch charge of 20 nC. The R&D program has the following components: The photoinjector and its photocathode, the superconducting linac cavity, start-to-end beam dynamics with magnetized electrons, electron cooling calculations including benchmarking experiments and development of a large superconducting solenoid. The photoinjector and linac cavity are being incorporated into an energy recovery linac aimed at demonstrating ampere class current at about 20 MeV.

Computational studies and optimization of wakefield accelerators

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.

Characteristics of an envelope model for laser-plasma accelerator simulation

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.