Mark E Dieckmann

Multidimensional electron beam-plasma instabilities in the relativistic regime

The interest in relativistic beam-plasma instabilities has been greatly rejuvenated over the past two decades by novel concepts in laboratory and space plasmas. Recent advances in this long-standing field are here reviewed from both theoretical and numerical points of view. The primary focus is on the two-dimensional spectrum of unstable electromagnetic waves growing within relativistic, unmagnetized, and uniform electron beam-plasma systems. Although the goal is to provide a unified picture of all instability classes at play, emphasis is put on the potentially dominant waves propagating obliquely to the beam direction, which have received little attention over the years. First, the basic derivation of the general dielectric function of a kinetic relativistic plasma is recalled. Next, an overview of two-dimensional unstable spectra associated with various beam-plasma distribution functions is given. Both cold-fluid and kinetic linear theory results are reported, the latter being based on waterbag and Maxwell–Juttner model distributions. The main properties of the competing modes (developing parallel, transverse, and oblique to the beam) are given, and their respective region of dominance in the system parameter space is explained. Later sections address particle-in-cell numerical simulations and the nonlinear evolution of multidimensional beam-plasma systems. The elementary structures generated by the various instability classes are first discussed in the case of reduced-geometry systems. Validation of linear theory is then illustrated in detail for large-scale systems, as is the multistaged character of the nonlinear phase. Finally, a collection of closely related beam-plasma problems involving additional physical effects is presented, and worthwhile directions of future research are outlined.

Table-Top Laser-Based Source of Femtosecond, Collimated, Ultrarelativistic Positron Beams

The generation of ultrarelativistic positron beams with short duration (τ(e+) ≃ 30  fs), small divergence (θ(e+) ≃ 3  mrad), and high density (n(e+) ≃ 10(14)-10(15)  cm(-3)) from a fully optical setup is reported. The detected positron beam propagates with a high-density electron beam and γ rays of similar spectral shape and peak energy, thus closely resembling the structure of an astrophysical leptonic jet. It is envisaged that this experimental evidence, besides the intrinsic relevance to laser-driven particle acceleration, may open the pathway for the small-scale study of astrophysical leptonic jets in the laboratory.

Oblique electromagnetic instabilities for a hot relativistic beam interacting with a hot and magnetized plasma

The temperature-dependent fluid model from Phys. Plasmas 13, 042106 (2006) is expanded in order to explore the oblique electromagnetic instabilities, which are driven by a hot relativistic electron beam that is interpenetrating a hot and magnetized plasma. The beam velocity vector is parallel to the magnetic-field direction. The results are restricted to nonrelativistic temperatures. The growth rates of all instabilities but the two-stream instability can be reduced by a strong magnetic field so that the distribution of unstable waves becomes almost one dimensional. For high beam densities, highly unstable oblique modes dominate the spectrum of unstable waves as long as omega(c)/omega(p)less than or similar to 2 gamma(3/2)(b), where omega(c) is the electron gyrofrequency, omega(p) is the electron plasma frequency, and gamma(b) is the relativistic factor of the beam. A uniform stabilization over the entire k space cannot be achieved.

Dynamics of Self-Generated, Large Amplitude Magnetic Fields Following High-Intensity Laser Matter Interaction

The dynamics of magnetic fields with an amplitude of several tens of megagauss, generated at both sides of a solid target irradiated with a high-intensity (~10(19) W/cm(2)) picosecond laser pulse, has been spatially and temporally resolved using a proton imaging technique. The amplitude of the magnetic fields is sufficiently large to have a constraining effect on the radial expansion of the plasma sheath at the target surfaces. These results, supported by numerical simulations and simple analytical modeling, may have implications for ion acceleration driven by the plasma sheath at the rear side of the target as well as for the laboratory study of self-collimated high-energy plasma jets.

Observation and characterization of laser-driven Phase Space Electron Holes

The direct observation and full characterization of a phase space electron hole (EH) generated during laser-matter interaction is presented. This structure, propagating in a tenuous, nonmagnetized plasma, has been detected via proton radiography during the irradiation with a ns laser pulse (Iλ2≈1014 W/cm2) of a gold hohlraum. This technique has allowed the simultaneous detection of propagation velocity, potential, and electron density spatial profile across the EH with fine spatial and temporal resolution allowing a detailed comparison with theoretical and numerical models.

Large-scale numerical simulations of ion beam instabilities in unmagnetized astrophysical plasmas

Collisionless quasiperpendicular shocks with magnetoacoustic Mach numbers exceeding a certain threshold are known to reflect a fraction of the upstream ion population. These reflected ions drive instabilities which, in a magnetized plasma, can give rise to electron acceleration. In the case of shocks associated with supernova remnants (SNRs), electrons energized in this way may provide a seed population for subsequent acceleration to highly relativistic energies. If the plasma is weakly magnetized, in the sense that the electron cyclotron frequency is much smaller than the electron plasma frequency ωp, a Buneman instability occurs at ωp. The nonlinear evolution of this instability is examined using particle-in-cell simulations, with initial parameters which are representative of SNR shocks. For simplicity, the magnetic field is taken to be strictly zero. It is shown that the instability saturates as a result of electrons being trapped by the wave potential. Subsequent evolution of the waves depends on the temperature of the background protons Ti and the size of the simulation box L. If Ti is comparable to the initial electron temperature Te, and L is equal to one Buneman wavelength λ0, the wave partially collapses into low frequency waves and backscattered waves at around ωp. If, on the other hand, Ti≫Te and L=λ0, two high frequency waves remain in the plasma. One of these waves, excited at a frequency slightly lower than ωp, may be a Bernstein–Greene–Kruskal mode. The other wave, excited at a frequency well above ωp, is driven by the relative streaming of trapped and untrapped electrons. In a simulation with L=4λ0, the Buneman wave collapses on a time scale consistent with the excitation of sideband instabilities. Highly energetic electrons were not observed in any of these simulations, suggesting that the Buneman instability can only produce strong electron acceleration in a magnetized plasma.Collisionless quasiperpendicular shocks with magnetoacoustic Mach numbers exceeding a certain threshold are known to reflect a fraction of the upstream ion population. These reflected ions drive instabilities which, in a magnetized plasma, can give rise to electron acceleration. In the case of shocks associated with supernova remnants (SNRs), electrons energized in this way may provide a seed population for subsequent acceleration to highly relativistic energies. If the plasma is weakly magnetized, in the sense that the electron cyclotron frequency is much smaller than the electron plasma frequency ωp, a Buneman instability occurs at ωp. The nonlinear evolution of this instability is examined using particle-in-cell simulations, with initial parameters which are representative of SNR shocks. For simplicity, the magnetic field is taken to be strictly zero. It is shown that the instability saturates as a result of electrons being trapped by the wave potential. Subsequent evolution of the waves depends on the...

Evolution of the fastest-growing relativistic mixed mode instability driven by a tenuous plasma beam in one and two dimensions

Particle-in-cell simulations confirm here that a mixed plasma mode is the fastest growing when a highly relativistic tenuous electron-proton beam interacts with an unmagnetized plasma. The mixed modes grow faster than the filamentation and two-stream modes in simulations with beam Lorentz factors Gamma of 4, 16, and 256, and are responsible for thermalizing the electrons. The mixed modes are followed to their saturation for the case of Gamma=4 and electron phase space holes are shown to form in the bulk plasma, while the electron beam becomes filamentary. The initial saturation is electrostatic in nature in the considered one- and two-dimensional geometries. Simulations performed with two different particle-in-cell simulation codes evidence that a finite grid instability couples energy into high-frequency electromagnetic waves, imposing simulation constraints.

Time-Resolved Characterization of the Formation of a Collisionless Shock

We report on the temporally and spatially resolved detection of the precursory stages that lead to the formation of an unmagnetized, supercritical collisionless shock in a laser-driven laboratory experiment. The measured evolution of the electrostatic potential associated with the shock unveils the transition from a current free double layer into a symmetric shock structure, stabilized by ion reflection at the shock front. Supported by a matching particle-in-cell simulation and theoretical considerations, we suggest that this process is analogous to ion reflection at supercritical collisionless shocks in supernova remnants.

The Formation of a Relativistic Partially Electromagnetic Planar Plasma Shock

Relativistically colliding plasma is modeled by particle-in-cell simulations in one and two spatial dimensions, with an ion-to-electron mass ratio of 400 and a temperature of 100 keV. The energy of an initial quasi-parallel magnetic field is 1% of the plasma kinetic energy. Energy dissipation by a growing wave pulse of mixed polarity, probably an oblique whistler wave, and different densities of the colliding plasma slabs result in the formation of an energetic electromagnetic structure within milliseconds. The structure, which develops for an initial collision speed of 0.9c, accelerates electrons to Lorentz factors of several hundred. A downstream region forms, separating the forward and reverse shocks. In this region, the plasma approaches an energy equipartition between electrons, ions, and the magnetic field. The electron energy spectrum -->N(E) resembles a power law at high energies, with an exponent close to –2.7, or -->N(E) E−2.7. The magnetic field reflects upstream ions, which form a beam and drag the electrons along to preserve the plasma quasineutrality. The forward and reverse shocks are asymmetric due to the unequal slab densities. The forward shock may be representative for the internal shocks of gamma-ray bursts.

The plasma filamentation instability in one dimension: nonlinear evolution

The plasma filamentation instability or beam-Weibel instability generates magnetic fields and accelerates particles in collisionless astrophysical plasma. This instability has been examined with multi-dimensional particle-in-cell (PIC) simulations, demonstrating the formation of current flux tubes. Such simulations could not model a statistically significant number of filaments. Here, we model with a PIC simulation the filamentation instability that is driven by nonrelativistic, cool electron beams in one spatial dimension at an unprecedented resolution. We show unambiguously that the gradient of the magnetic pressure which develops during the quasi-linear evolution of the filamentation instability, gives rise to an electrostatic field component. The interplay of the magnetic and electrostatic fields results in a wavenumber spectrum of the magnetic field that is a power-law, which has been reported previously for multi-dimensional PIC simulations. The magnetic field power spectrum decreases with the exponent -5.7 and that of the electrostatic field with -3.8, yielding a ratio of 3:2. The electromagnetic fields thermalize the electrons. The electrons develop a velocity distribution in the simulation direction that decreases exponentially at low speeds and faster at high speeds. The filamentation instability can thus not efficiently accelerate electrons to high energies. The filaments develop into a stationary final state. The probability distribution of the filament sizes is a Gumbel distribution. In astrophysical settings, this implies that the long exponential tail of this distribution may lead with a reasonable probability to large current filaments, if the filamentation instability develops in a large enough volume. The coherent magnetic fields of large filaments are required to explain the synchrotron emissions of gamma ray bursts.