To predict the evolution of electron clouds and their effect on the beam, the high energy physics community has relied so far on the complementary use of 'buildup' and 'single/multi-bunch instability' reduced descriptions. The former describes the evolution of electron clouds at a given location in the ring, or 'station', under the influence of prescribed beams and external fields [1], while the latter (sometimes also referred as the 'quasi-static' approximation [2]) follows the interaction between the beams and the electron clouds around the accelerator with prescribed initial distributions of electrons, assumed to be concentrated at a number of discrete 'stations' around the ring. Examples of single bunch instability codes include HEADTAIL [3], QuickPIC [4, 5], and PEHTS [6]. By contrast, a fully self-consistent approach, in which both the electron cloud and beam distributions evolve simultaneously under their mutual influence without any restriction on their relative motion, is required for modeling the interaction of high-intensity beams with electron clouds for heavy-ion beam-driven fusion and warm-dense matter science. This community has relied on the use of Particle-In-Cell (PIC) methods through the development and use of the WARP-POSINST code suite [1, 7, 8]. The development of novel numerical techniques (including adaptive mesh refinement, and a new 'drift-Lorentz' particle mover for tracking charged particles in magnetic fields using large time steps) has enabled the first application of WARP-POSINST to the fully self-consistent modeling of beams and electron clouds in high energy accelerators [9], albeit for only a few betatron oscillations. It was recently observed [10] that there exists a preferred frame of reference which minimizes the number of computer operations needed to simulate the interaction of relativistic objects. This opens the possibility of reducing the cost of fully self-consistent simulations for the interaction of ultrarelativistic beams with electron cloud by orders of magnitude. The computational cost of the fully self-consistent mode is then predicted to be comparable to that of the quasi-static mode, assuming that several stations per betatron period are needed. During the workshop, there was some debate about the number of stations per betatron period that are needed when using the quasi-static mode. The argument was made that if there is less than one station per betatron period, then artificial resonances can be triggered and the resulting emittance growth provides an upper bound. The emittance growth thus obtained will fall either above or below the operational requirements of the machine. In the latter case, one can conclude that the electron effect that has been simulated is of no concern. However, if the emittance growth that was obtained is above the threshold, then the results become inconclusive, and simulations which resolve the betatron motion are then needed. In this case, according to [10], the fully self-consistent approach becomes an option. The aim of this paper is to investigate whether this option is indeed practical.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2012

Non-paraxial Pulse Broadening in a Solenoid Focusing Element

Abstract not provided

43 PARTICLE ACCELERATORS↗

DOE OSTI2011

Modeling of 10 GeV-1 TeV laser-plasma accelerators using Lorentz boosted simulations

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.

43 PARTICLE ACCELERATORS↗

DOE OSTI2011

Numerical methods for instability mitigation in the modeling of laser wakefield accelerators in a Lorentz-boosted frame

Modeling of laser-plasma wakefield accelerators in an optimal frame of reference has been shown to produce orders of magnitude speed-up of calculations from first principles. Obtaining these speedups required mitigation of a high-frequency instability that otherwise limits effectiveness. In this paper, methods are presented which mitigated the observed instability, including an electromagnetic solver with tunable coefficients, its extension to accommodate Perfectly Matched Layers and Friedman's damping algorithms, as well as an efficient large bandwidth digital filter. It is observed that choosing the frame of the wake as the frame of reference allows for higher levels of filtering or damping than is possible in other frames for the same accuracy. Detailed testing also revealed the existence of a singular time step at which the instability level is minimized, independently of numerical dispersion. A combination of the techniques presented in this paper prove to be very efficient at controlling the instability, allowing for efficient direct modeling of 10 GeV class laser plasma accelerator stages. The methods developed in this paper may have broader application, to other Lorentz-boosted simulations and Particle-In-Cell simulations in general.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2011

NDCX-II, an Induction Linac for HEDP and IFE Research

The Heavy Ion Fusion Science Virtual National Laboratory in the USA is constructing a new Neutralized Drift Compression eXperiment (NDCX-II) at LBNL. This facility is being developed for high energy density physics and inertial fusion energy research. The 12 m long induction linac in NDCX-II will produce a Li{sup +} beam pulse, at energies of 1.2-3 MeV, to heat target material to the warm dense matter regime ({approx} 1 eV). By making use of special acceleration voltage waveforms, 2.5T solenoid focusing, and neutralized drift compression, 20 - 50 nC of beam charge from the ion source will be compressed longitudinally and radially to achieve a subnanosecond pulse length and mm-scale target spot size. The original Neutralized Drift Compression Experiment (NDCX-I) has successfully demonstrated simultaneous radial and longitudinal compression by imparting a velocity ramp to the ion beam, which then drifts in a neutralizing plasma to and through the final focussing solenoid and onto the target. At higher kinetic energy and current, NDCX-II will offer more than 100 times the peak energy fluence on target of NDCX-I. NDCX-II makes use of many parts from the decommissioned Advanced Test Accelerator (ATA) at LLNL. It includes 27 lattice periods between the injector and the neutralized drift compression section (Figure 1). There are 12 energized induction cells, 9 inactive cells which provide drift space, and 6 diagnostic cells which provide beam diagnostics and pumping. Custom pulsed power systems generate ramped waveforms for the first 7 induction cells, so as to quickly compress the beam from 600 ns at the injector down to 70 ns. After this compression, the high voltages of the ATA Blumleins are then used to rapidly add energy to the beam. The Blumleins were designed to match the ferrite core volt-seconds with pulses up to 250 kV and a fixed FWHM of 70 ns. The machine is limited to a pulse repetition rate of once every 20 seconds due to cooling requirements. The NDCX-II beam is highly space-charge dominated. The 1-D ASP code was used to synthesize high voltage waveform for acceleration, while the 3-D Warp particle-in-cell code was used for detailed design of the lattice. The Li{sup +} ion was chosen because its Bragg Peak energy (at {approx} 2 MeV) coincides with the NDCX-II beam energy. The 130 keV injector will have a 10.9 cm diameter ion source. Testing of small (0.64 cm diameter) lithium doped alumino-silicate ion sources has demonstrated the current density ({approx} 1 mA/cm{sup 2}) used in the design, with acceptable lifetime. A 7.6 cm diameter source has been successfully produced to verify that the coating method can be applied to such a large emitting area. The ion source will operate at {approx} 1275 C; thus a significant effort was made in the design to manage the 4 kW heating power and the associated cooling requirements. In modifying the ATA induction cells for NDCX-II, the low-field DC solenoids were replaced with 2.5 T pulsed solenoids. The beam pipe diameter was decreased in order to reduce the axial extent of the solenoid fringe fields and to make room for water cooling. In addition, an outer copper cylinder (water-cooled) was used to exclude the solenoid magnetic flux from the ferrite cores. Precise alignment is essential because the beam has a large energy spread due to the rapid pulse compression, such that misalignments lead to corkscrew deformation of the beam and reduced intensity at focus. A novel pulsed-wire measurement method is used to align the pulsed solenoid magnets. Alignment accuracy has been demonstrated to within 100 {micro}m of the induction cell axis. The neutralized drift compression region after the last induction cell is approximately 1.2 m long and includes ferroelectric plasma sources (FEPS) fabricated by PPPL similar to those successfully operating in NDCX-I. The 8-T final focus pulsed solenoid, filtered cathodic arc plasma sources (FCAPS), and target chamber from NDCX-I are to be relocated to NDCX-II. The NDCX-II project started in July 2009 and is expected to complete in fall of 2011. As future funds become available, additional induction cells and pulsed power systems will be added to increase the beam energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2011

INERTIAL FUSION DRIVEN BY INTENSE HEAVY-ION BEAMS

Intense heavy-ion beams have long been considered a promising driver option for inertial-fusion energy production. This paper briefly compares inertial confinement fusion (ICF) to the more-familiar magnetic-confinement approach and presents some advantages of using beams of heavy ions to drive ICF instead of lasers. Key design choices in heavy-ion fusion (HIF) facilities are discussed, particularly the type of accelerator. We then review experiments carried out at Lawrence Berkeley National Laboratory (LBNL) over the past thirty years to understand various aspects of HIF driver physics. A brief review follows of present HIF research in the US and abroad, focusing on a new facility, NDCX-II, being built at LBNL to study the physics of warm dense matter heated by ions, as well as aspects of HIF target physics. Future research directions are briefly summarized.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2011

Effects of hyperbolic rotation in Minkowski space on the modeling of plasma accelerators in a Lorentz boosted frame

The effects of hyperbolic rotation in Minkowski space resulting from the use of Lorentz boosted frames of calculation on laser propagation in plasmas are analyzed. Selection of a boost frame at the laser group velocity is shown to alter the laser spectrum, allowing the use of higher boost velocities. The technique is applied to simulations of laser driven plasma wakefield accelerators, which promise much smaller machines and whose development requires detailed simulations that challenge or exceed current capabilities. Speedups approaching the theoretical optima are demonstrated, producing the first direct simulations of stages up to 1 TeV. This is made possible by a million times speedup thanks to a frame boost with a relativistic factor {gamma}{sub b} as high as 1300, taking advantage of the rotation to mitigate an instability that limited previous work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2010

Electron-Driven Fast-Ignition Modeling with Realistic Electron-Source

Abstract not provided

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2010

Effects of Hyperbolic Rotation in Minkowski Space on the Modeling of Plasma Accelerators in a Lorentz Boosted Frame

Laser driven plasma accelerators promise much shorter particle accelerators but their development requires detailed simulations that challenge or exceed current capabilities. We report the first direct simulations of stages up to 1 TeV from simulations using a Lorentz boosted calculation frame resulting in a million times speedup, thanks to a frame boost as high as gamma = 1300. Effects of the hyperbolic rotation in Minkowski space resulting from the frame boost on the laser propagation in the plasma is shown to be key in the mitigation of a numerical instability that was limiting previous attempts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2010

Modeling laser wakefield accelerators in a Lorentz boosted frame

Modeling of laser-plasma wakefield accelerators in an optimal frame of reference [1] is shown to produce orders of magnitude speed-up of calculations from first principles. Obtaining these speedups requires mitigation of a high frequency instability that otherwise limits effectiveness in addition to solutions for handling data input and output in a relativistically boosted frame of reference. The observed high-frequency instability is mitigated using methods including an electromagnetic solver with tunable coefficients, its extension to accomodate Perfectly Matched Layers and Friedman's damping algorithms, as well as an efficient large bandwidth digital filter. It is shown that choosing theframe of the wake as the frame of reference allows for higher levels of filtering and damping than is possible in other frames for the same accuracy. Detailed testing also revealed serendipitously the existence of a singular time step at which the instability level is minimized, independently of numerical dispersion, thus indicating that the observed instability may not be due primarily to Numerical Cerenkov as has been conjectured. The techniques developed for Cerenkov mitigation prove nonetheless to be very efficient at controlling the instability. Using these techniques, agreement at the percentage level is demonstrated between simulations using different frames of reference, with speedups reaching two orders of magnitude for a 0.1 GeV class stages. The method then allows direct and efficient full-scale modeling of deeply depleted laser-plasma stages of 10 GeV-1 TeV for the first time, verifying the scaling of plasma accelerators to very high energies. 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗