Interaction of Ultra Relativistic e- e+ Fireball Beam with Plasma

Ab initio simulations of the propagation in a plasma of a soon to be available relativistic electron-positron beam or fireball beam provide an effective mean for the study of microphysics relevant to astrophysical scenarios. We show that the current filamentation instability associated with some of these scenarios reaches saturation after only 10 cm of propagation in a typical laboratory plasma with a density 10^17/cc. The different regimes of the instability, from the purely transverse to the mixed mode filamentation, can be accessed by varying the background plasma density. The instability generates large local plasma gradients, intense transverse magnetic fields, and enhanced emission of radiation. We suggest that these effects may be observed experimentally for the first time.

Ab initio simulations of the propagation in a plasma of a soon to be available relativistic electronpositron beam or fireball beam provide an effective mean for the study of microphysics relevant to astrophysical scenarios. We show that the current filamentation instability associated with some of these scenarios reaches saturation after only 10 cm of propagation in a typical laboratory plasma with a density ∼ 10 17 cm −3 . The different regimes of the instability, from the purely transverse to the mixed mode filamentation, can be accessed by varying the background plasma density. The instability generates large local plasma gradients, intense transverse magnetic fields, and enhanced emission of radiation. We suggest that these effects may be observed experimentally for the first time. Several astrophysical scenarios lead to extreme physical regimes, typically observed on Earth in the form of radiation and cosmic rays. These regimes encompass a set of phenomena such as magnetic field generation, shock formation, energy transfer processes, and non-thermal particle acceleration (for a review, see [1]). In the particular case of the fireball model of gamma ray bursts (GRBs) [2], the kinetic energy of an ultra-relativistic plasma shell, with an arbitrary mixture of electronspositrons-ions (e − e + p + ), is converted into radiation as moving shells collide, but the specific conversion mechanism is still an open question. Relativistic flows are also frequent in shock waves and pulsar wind nebulae, where relativistic shells interact with a background plasma [3]. It is very difficult to reproduce these astrophysical conditions in the laboratory, and the studies of the nonlinear physical phenomena are essentially simplified analytical models and numerical simulations (see [4][5][6] and references therein). Identifying the laboratory conditions that can validate the conclusions reached in previous studies is thus of paramount importance. Progress in laser technology, for instance, already suggests the possibility of exploring experimentally scaled-down astrophysical phenomena in laser-plasma interactions [7]. The occurrence of the current filamentation instability has recently been observed in experiments with a low energy electron beam and a capillary discharge plasma [8].
In this Letter we focus on a scenario similar to that widely believed to be present in, and at the origin of GRBs, by examining the collision of a relativistic e − e + beam or neutral plasma (that we call a fireball beam) mimicking a realistic plasma shell, with a static plasma consisting of e − and p + . The interaction leads to current filamentation instability (CFI), or Weibel instability [9,10], which generates very large magnetic fields as * Electronic address: muggli@mpp.mpg.de † Electronic address: luis.silva@ist.utl.pt the beam plasma interaction evolves. The self-consistent evolution of electric (E) and magnetic (B) fields, and the resulting radiation generation as particles propagate in CFI driven turbulence are observed. Here, we consider conditions that will soon become available in the laboratory, and focus on the direct observation of the CFI of an e − e + neutral beam. We also examine the possible experimental detection of the nonlinear stage of this instability as a function of its key parameters, through imaging of the associated plasma gradients, detection of the beam radiation, and the B-field of the beam filaments. Our results show that the role and the dynamics of the plasma microinstabilities associated with the collision of a relativistic fireball with a plasma may be probed in the laboratory.
Recent theoretical [10,11] and numerical results [12,13] have shown the relevance of the CFI for GRBs scenarios and for the onset of relativistic shocks in unmagnetized plasmas [5,6]. Moreover, the consequences of the excitation of the mixed mode, or tilted filamentation [14], in the long time evolution of the generated E/Bfields remains to be addressed; it was suggested that this preferential mode excitation will lead to significant beam spraying [12], but the impact on the saturated level of the fields was not discussed. Experimental evidence for the radiation signatures from Weibel turbulence [15], the energy transfer rate from the fireball to the B-field, and the long time evolution of the self-generated E and Bfields are critical inputs to the existing models, to perform direct comparisons with astronomical observations, and to assess the relevance of the CFI to relativistic astrophysics.
Relativistic e − beams are available in many laboratories around the world, while e + beams are not. Recent numerical studies of a plasma-based, e − e + accelerator concept based on the plasma wakefield accelerator (PWFA) [16] have revealed that it may be advantageous to accelerate a e + bunch on the wake driven by an e − bunch [17]. Ultra-relativistic e − and e + bunches suitable to test this acceleration scheme are available at the arXiv:1306.4380v1 [physics.plasm-ph] 18 Jun 2013 SLAC National Accelerator Laboratory. For this test, the distance between the e − and the e + bunch must be adjustable and on the order of the plasma wavelength or about 100 µm. A double or sailboat magnetic chicane has been developed [18], that allows for the adjustment of the spacing between the two bunches, and may be used to overlap the two bunches with equal charge, both in space and time to effectively create a relativistic fireball beam. This will make possible the first ever collision between relativistic neutral plasmas in the laboratory: a relativistic e − e + plasma onto an e − p + plasma at rest, separating the effects of the space charge fields associated with a charged beam.
We investigate the propagation of the SLAC fireball beam in a pre-formed plasma with numerical simulations performed with the fully relativistic, fully electromagnetic, and massivelly parallel particle-in-cell (PIC) code OSIRIS [19]. This simulation framework has been extensively used for studies of laser/beam plasma interaction (e.g., [20]), and astrophysical regimes (e.g., [12,21]), among others. The system is studied numerically with a 205 × 205 × 82 µm 3 window moving at the speed of light along the z-direction, and discretized in 400 × 400 × 80 cells with absorbing boundary conditions for the fields and for the particles in the transverse x, y directions. The fireball beam is defined with Gaussian profiles in all directions with rms sizes: σ x = σ y = 2σ z = 2 c/ω pe = 20.4 µm, where ω pe = (n e e 2 / 0 m e ) 1/2 is the e − pulsation of the rest plasma with density n e = 2.7×10 17 cm −3 . The standard beam used in the simulations has 1.8 × 10 10 e − , and the same number of e + , all with an incoming energy of 29 GeV and a normalized emittance of 2 × 10 −5 mrad, corresponding to a peak beam density n b = n e , and a transverse thermal spread v th /c = 1.7 × 10 −5 . A total of ∼ 7 × 10 7 simulation particles (plasma and beam) is pushed for ∼ 10 4 c/ω pe 10 cm of pre-formed plasma (20 cm were also simulated to confirm the saturated state parameters). The time step is 0.033/ω pe . The neutrality of the beam guaranties its propagation at constant radius. A background of fixed p + is assumed for the pre-formed plasma: quantitative variations below 1% were obtained for the standard case when using mobile p + . Note also that, as in astrophysics, these parameters correspond to a collisionless fireball-plasma interaction: ν ei /ω pe = O(10 −18 ), where ν ei is the beam e − , e + -background p + collision frequency. Fig. 1a-c show the structure of the fireball beam after 10 cm propagation in the laboratory plasma, or, equivalently, to the propagation of a fireball with a density 1 cm −3 in > 50 km in the background density of 1 cm −3 . The CFI generates well-defined current (and density) filaments, which size increases as the beam propagates in the plasma, and may grow to a thickness above 5 µm 0.5 c/ω pe . These conditions correspond to a beam with σ r /(c/ω pe ) 2, which explains the few filaments obtained at saturation, reached when the filaments coalescence ceases and the B-field energy remains constant. The large currents associated with the beam filaments generate local B-fields up to 2 MGauss (Fig. 1d). The space charge separation, also associated with the filaments, leads to radial E-fields as high as 5 × 10 8 V/cm. The presence of oblique modes/tilted filamentation [14] is clear on Fig. 1c showing that, as expected, the beam can excite a combination of transverse (filamentation) and longitudinal (two-stream-like) instabilities. Finally, a system of filaments is present in the background plasma, behind the beam, evidencing a 3D structure. As the filaments merge, the space-charge separation leads to the plasma blowout and to the generation of strong E-fields.
This scenario is in stark contrast with that of a pure e − beam interacting with the same plasma, as in the recent PWFA experiments [22]. Simulations and experiments for an e − beam with the same transverse size show that, for these parameters, the beam drives strong plasma wakefields that focus the beam to a narrow radius after one quarter betatron wavelength λ β /4 = (2 3/2 πγ 1/2 c/ω pe )/4 5 mm, and the beam envelope experiences oscillations along the plasma with period λ β /2. No CFI is observed under these circumstances.
In Fig. 2, we present the evolution with propagation distance in the plasma of the total normalized energy in the B-field, B for different beam/plasma parameters, illustrating the exponential growth and saturation within the 10 cm range. The growth rates (Γ std /ω pe 2.0 × 10 −3 , Γ hot /ω pe 1.7 × 10 −3 , Γ high /ω pe 2.3 × 10 −3 ) are within the range predicted for this configuration for the purely transverse CFI ( with β th = v th /c the particle thermal rms spread of velocity [12]). A more detailed analysis reveals, however, that for higher plasma densities (keeping the beam density fixed) the growth rate is higher, but the saturated level of the B-field is lower. The former is an indication of the spatial-temporal character of the instability in this configuration, while the latter is an evidence for the different saturation mechanisms involved when the mixed mode/tilted filamention is dominant [14]. The finite transverse dimension of the beam (not its length) determines: (i) the longest wavenumber that can be excited, (ii) the typical noise source for the instability. Since the beam is quite cold, for wavenumbers such that k ≤ c/ω pe the growth rate is already close to its maximum value. The finite length of the beam impacts the two-stream mode (or in the more general form of the filamentation instability) the oblique mode. However, there is no theory for the excitation of these modes for finitelength finite-width modes, and thus this work motivates further theoretical developments of a spatio-temporal theory for the fireball beam since it does not exist [23]. The spatial-temporal theory for the two-stream instability [24] predicts that an E-field perturbation excited at the vacuum/plasma interface (x = 0 and t = 0) will grow with ∝ exp[iω pe ψ] exp[(3 √ 3/4)(ψx 2 n e /n) 1/3 ω pe /γ] where ψ = t − z/v b is the distance to the head of the beam. In the beam region, for the same distance in the laboratory and assuming the same initial perturbation in the longitudinal E-field, the amplified field in the stan- dard case is approximately twice the amplified field in the high plasma density case. This is consistent with what we observe in the simulations, namely the fact that the mixed mode has clearly developed more strongly in the high density scenario. The coupling of the excited longitudinal field with the transverse field leads to the excitation of the mixed mode [14], as clearly seen in the high density case (Fig. 2): the filaments are tilted, which indicates that the particles can detrap more easily, leading to a lower current and thus to a lower saturated B-field.
The interaction of the relativistic e − and e + with the x + B 2 y + B 2 z ) normalized to the kinetic energy of the particles p = (γ0 − 1)V b , (V b the volume of the beam) for different beam and plasma parameters. Values are normalized to B0, where B0 is the field when the growth becomes exponential (after ∼ 0.1 mm, or 10/ωpe). Standard case (solid line): fireball beam with 2 × 10 −5 m-rad emittance, in a plasma with ne = 2.7 × 10 17 cm −3 (which also defines the baseline density for the normalization). The dotted line illustrates the linear growth rate. Slices of the density in the middle of the beam after ≈ 1.5 cm of plasma (plotted in blue) illustrate the difference in the instability structure. The inset includes the trajectories of two fireball electrons for the standard case.
B-fields confining the current filaments leads to the emission of synchrotron radiation. The oscillatory motion of the charges in the transverse directions due to the radial E-fields associated with the filaments of opposite charges leads to the emission of betatron radiation. In both cases, the radiation is incoherent with the wiggler strength parameter K ≡ a β = γk β r 1, for k β the betatron wavenumber and r the orbit amplitude. The spectra have a photon critical energy E syn = ω syn = where η 1 describes the typical radius of the filaments (in units of c/ω pe ). Simulation results indicate that these two radiation processes might not be distinguishable [25], at least for the initial stage, since E ⊥ and B ⊥ grow together. After a significant field growth in the CFI driven turbulence, however, the field structure may lead to different spectral signatures, as previously hinted in [26,27]. We briefly describe some of the particular aspects of the fireball beam diagnostic implementation. For PWFA applications the separation between the e − and the e + bunches must be of the order of a plasma wavelength (∼ 100 µm). Such a small spacing between the bunches can be achieved with two interleaved magnetic chicanes with a coarse path length difference of the order of the bunch separation in the accelerator (∼ 5 cm), and with fine magnetic adjustments [18]. These adjustments may also be used to overlap the two bunches in time and create the relativistic, neutral e − e + fireball beam. The optimal temporal overlap is achieved by minimizing the coherent transition radiation the bunches emit when traversing a thin metallic foil located after the double chicane. The transverse overlap is obtained by imaging the incoherent optical transition radiation the bunches emit when traversing two thin foils located before and after the plasma. For PWFA experiments the beam ionizes a lithium vapor and creates the plasma by fieldionization [22], while the neutral fireball beam requires a pre-ionized plasma. Pre-ionization can be achieved by photo-ionization of a lithium vapor with an ultra-violet laser pulse [28].
The filamentation of the beam is the most obvious indication of the CFI occurrence. The filaments, however, have a relatively small transverse size of ∼ 5 µm and, because of their emittance, diverge and overlap rapidly when exiting the plasma. To detect them inside the plasma, the strong plasma e − density gradients associated with the beam filamentation can be visualized with Schlieren shadowgraphy [29] using a laser pulse propagating perpendicularly to the fireball beam path. The laser light is weakly deflected by the index of refraction variations corresponding to the e − density modulation.
The filamentation of the beam results in the generation of large B-fields in the plane perpendicular to the filaments themselves (see Fig.1d), which can be visualized by analyzing the polarization of a linearly polarized probe laser pulse traveling perpendicularly to the fireball beam, the same that is used for the Schlieren shadowgraphy. The laser light experiences Faraday rotation caused by the component of the filaments B-fields parallel to the laser propagation direction (see Fig.1e). Even though the B-field pattern is related to the structure of the random filaments, the effect computed from the simulation results for B y (y) produces an image similar to Fig. 1e and is clearly visible. Faraday rotation has been used to sample the B-fields generated in a laser wakefield experiment with similar parameters [30]. The plasma density and B-field structure may be sampled along the e − beam path by moving the intersection point between the probe laser pulse and the plasma with a time resolution equal to the laser pulse length (fs) and a longitudinal resolution of the order of the probe beam size (mm), thereby giving access to the growth of the instability. The excess radiation associated with the oscillation of the e − and e + in the B-field and in the filaments can be directly observed using standard x-ray detection methods, similar to those that were used to detect synchrotron or betatron radiation in PWFA experiments [31]. Finally, the B-field growth occurs at the expense of beam energy. In our simulations the beam looses 6-11% energy, i.e., 2-3 GeV, in the standard and high temperature cases, respectively. These beam energy changes can be measured using an imaging magnetic spectrometer, as in previous PWFA experiments [32].
In conclusion, we have shown that the e − e + or fireball beam and plasma system that will be developed for PWFA experiments may also be used to produce in the laboratory a scenario relevant to test the very important microphysics issues of relativistic astrophysical phenomena. As a result of the CFI, the incoming fireball beam filamentation occurs over a plasma length of only a few cm. The current filaments generate large B-fields that lead to the enhanced emission of synchrotron and betatron radiation. Initial considerations indicate that the beam filamentation, the B-field generation, and the associated beam energy loss, as well as the excess radiation can in principle be observed in a single experiment. We