Search Results (10 total)

Unphysical heating and macroparticle trapping that arise in the numerical modeling of laser wakefield accelerators using particle-in-cell codes are investigated. A dark current free laser wakefield accelerator stage, in which no trapping of background plasma electrons into the plasma wave should occur, and a highly nonlinear cavitated wake with self-trapping, are modeled. Numerical errors can lead to errors in the macroparticle orbits in both phase and momentum. These errors grow as a function of distance behind the drive laser and can be large enough to result in unphysical trapping in the plasma wake. The resulting numerical heating in intense short-pulse laser-plasma interactions grows much faster and to a higher level than the known numerical grid heating of an initially warm plasma in an undriven system. The amount of heating, at least in the region immediately behind the laser pulse, can, in general, be decreased by decreasing the grid size, increasing the number of particles per cell, or using smoother interpolation methods. The effect of numerical heating on macroparticle trapping is less severe in a highly nonlinear cavitated wake, since trapping occurs in the first plasma wave period immediately behind the laser pulse.

The effects of radiation reaction on electron beam dynamics are studied in the context of plasma-based accelerators. Electrons accelerated in a plasma channel undergo transverse betatron oscillations due to strong focusing forces. These oscillations lead to emission by the electrons of synchrotron radiation, with a corresponding energy loss that affects the beam properties. An analytical model for the single particle orbits and beam moments including the classical radiation reaction force is derived and compared to the results of a particle transport code. Since the betatron amplitude depends on the initial transverse position of the electron, the resulting radiation can increase the relative energy spread of the beam to significant levels (e.g., several percent). This effect can be diminished by matching the beam into the channel, which could require micron sized beam radii for typical values of the beam emittance and plasma density.

A warm, relativistic fluid theory of a nonequilibrium, collisionless plasma is developed to analyze nonlinear plasma waves excited by intense drive beams. The maximum amplitude and wavelength are calculated for nonrelativistic plasma temperatures and arbitrary plasma wave phase velocities. The maximum amplitude is shown to increase in the presence of a laser field. These results set a limit to the achievable gradient in plasma-based accelerators.

We develop a low-temperature fluidlike plasma model without recourse to a collisional closure. The equations are closed by treating the momentum spread asymptotically. This model inherits the Hamiltonian structure, including Casimir invariants of the Vlasov–Maxwell theory. We study temperature evolution in the wake of an intense laser pulse propagating in a plasma. We show that the momentum spread is intrinsically anisotropic and that, for conditions corresponding to recent experiments, modest heating occurs.

The effect of asymmetric laser pulses on electron yield from a laser wakefield accelerator has been experimentally studied using >10¹⁹   cm^-3 plasmas and a 10 TW, >45  fs, Ti∶Al₂O₃ laser. The laser pulse shape was controlled through nonlinear chirp with a grating pair compressor. Pulses (76 fs FWHM) with a steep rise and positive chirp were found to significantly enhance the electron yield compared to pulses with a gentle rise and negative chirp. Theory and simulation show that fast rising pulses can generate larger amplitude wakes that seed the growth of the self-modulation instability, and that frequency chirp is of minimal importance for the experimental parameters.

Spontaneous radiation emitted from relativistic electrons undergoing betatron motion in a plasma-focusing channel is analyzed, and applications to plasma wake-field accelerator experiments and to the ion-channel laser (ICL) are discussed. Important similarities and differences between a free electron laser (FEL) and an ICL are delineated. It is shown that the frequency of spontaneous radiation is a strong function of the betatron strength parameter a_β, which plays a role similar to that of the wiggler strength parameter in a conventional FEL. For a_β≳1, radiation is emitted in numerous harmonics. Furthermore, a_β is proportional to the amplitude of the betatron orbit, which varies for every electron in the beam. The radiation spectrum emitted from an electron beam is calculated by averaging the single-electron spectrum over the electron distribution. This leads to a frequency broadening of the radiation spectrum, which places serious limits on the possibility of realizing an ICL.

We present 2D simulations of both beam-driven and laser-driven plasma wakefield accelerators, using the object-oriented particle-in-cell code XOOPIC, which is time explicit, fully electromagnetic, and capable of running on massively parallel supercomputers. Simulations of laser-driven wakefields with low \(∼10¹⁶ W/cm²\) and high \(∼10¹⁸ W/cm²\) peak intensity laser pulses are conducted in slab geometry, showing agreement with theory and fluid simulations. Simulations of the E-157 beam wakefield experiment at the Stanford Linear Accelerator Center, in which a 30 GeV electron beam passes through 1 m of preionized lithium plasma, are conducted in cylindrical geometry, obtaining good agreement with previous work. We briefly describe some of the more significant modifications to XOOPIC required by this work, and summarize the issues relevant to modeling relativistic electron-neutral collisions in a particle-in-cell code.

Nonparaxial propagation of ultrashort, high-power laser pulses in plasma channels is examined. In the adiabatic limit, pulse energy conservation, nonlinear group velocity, damped betatron oscillations, self-steepening, self-phase modulation, and shock formation are analyzed. In the nonadiabatic limit, the coupling of forward Raman scattering (FRS) and the self-modulation instability (SMI) is analyzed and growth rates are derived, including regimes of reduced growth. The SMI is found to dominate FRS in most regimes of interest.

A method is described for predicting statistical properties of turbulence. Collections of Fourier amplitudes are represented by nonuniformly spaced modes with enhanced coupling coefficients. The statistics of the full dynamics can be recovered from the time-averaged predictions of the reduced model. A Liouville theorem leads to inviscid equipartition solutions. Excellent agreement is obtained with two-dimensional forced-dissipative pseudospectral simulations. For the two-dimensional enstrophy cascade, logarithmic corrections to the high-order structure functions are observed.

The quantum Liouville equation for an n-level atomic system driven by external fields has a nontrivial kinematic structure; the quantities tr ρ^j, j  =  1,2,…,n remain constant in time, independent of the Hamiltonian. These invariants are physically significant; the qualitative character of the solution depends on their existence. A generic numerical method will not, in general, preserve these invariants. We present a numerical technique which evolves the density matrix via unitary transformations thus exactly preserving these invariants to all orders in the time step.