INTRODUCTION
One notices a rather high level of current interest for the interaction physics of intense and relativistic electron beams (REB) with superdense deuterium + tritium (DT) plasmas. Corresponding target electron n p density ranges between 1023 electrons/cm−3 for the outer layer of compressed DT core and 1026 electrons/cm−3 at the center. Such a concern is basically motived by the fast ignition scenario (FIS) for inertial confinement fusion (ICF), driven by intense and laser produced REB with number density n b~1022 electrons/cm−3 (Tabak et al., Reference Tabak, Hammer, Glinsky, Kruer, Wilks, Woodworth, Campbell, Perry and Mason1994; Deutsch et al., Reference Deutsch, Furukawa, Mima and Nishihara1996, Reference Deutsch, Furukawa, Mima and Nishihara2000; Eliezer et al., Reference Eliezer, Murakami and Martinez-Val2007; Johzaki et al., Reference Johzaki, Sakagami, Nagatomo and Mima2007; Flippo et al., Reference Flippo, Hegelich, Albright, Yin, Gautier, Letzring, Schollmeier, Schreiber, Schulze and Fernandez2007; Hora, Reference Hora2007; Yu et al., Reference Yu, Yu, Xu, Tian, Chen and Wong2007). Such a scheme outlines to the extreme the separation of compression driven by lasers, heavy ion beams or Z pinches, from the controlled ignition through REB, which thus relieves to a large extent the energy request on initial precompression. The successful completion of such an approach to ICF demands a careful control of any electromagnetic instability susceptible to evolve very swiftly within the REB propagation, and nearly instantaneously, a significant amount of the incoming REB translation energy. That process is not only quantitative; it is also highly likely to occur through a filamentation with a substantial degradation of beam focusing capabilities (Bret et al., Reference Bret, Firpo and Deutsch2007; Malik, Reference Malik2007).
Recently, our group has devoted a very thorough scrutiny to the interplay of filamentation (Weibel, Reference Weibel1959) electromagnetic instabilities with the more mundane two-stream (TS) and longitudinal instability (Deutsch, Reference Deutsch2004; Firpo et al., Reference Firpo, Lifschitz, Lefebvre and Deutsch2006). As a result, a previously unnoticed oblique ridge has thus been unraveled in the space delimited by the growth rate and the two corresponding wave number coordinates. However, those studies are often restricted to a collisionless situation. A complete exploration, including collisions with TS altogether, is also worked out.
Here, we intend to focus attention on the very large densities implied in a realistic beam-target interaction of FIS interest (Tabak et al., Reference Tabak, Hammer, Glinsky, Kruer, Wilks, Woodworth, Campbell, Perry and Mason1994; Kodama et al., Reference Kodama, Norreys, Mima, Dangor, Evans, Fujita, Kitagawa, Krushelnik, Miyakoshi, Norimatsu, Rose, Shozaki, Shigemon, Sunahara, Tampo, Tanaka, Toyama, Yamanaka and Zepf2001). Despite that beam particle density n b always fulfill a weak beam approximation n b << n p, it should be taken around n b~1022 electrons/cm−3 to secure on a 10 ns time scale several tens of kJ to the off-center hot spot requested for DT ignition. Such a situation implies that the usual collisionless assumption used in working out electromagnetic instabilities through the Vlasov equation be removed. At this juncture, a few earlier endeavors dedicated to the interplay of electromagnetic instabilities with electron-electron collisions should be recalled (Taguchi et al., Reference Taguchi, Antonsen, Liu and Mima2001; Kono & Ichikawa, Reference Kono and Ichikawa1973; Mason, Reference Mason2006).
LINEAR AND QUASI-LINEAR FILAMENTATION GROWTH
Let us now consider a current neutral beam-plasma system. The relativistic electron beam propagates with the velocity vdb and the plasma return current flows with vdp. It is reasonable to assume that an electromagnetic mode has k normal to vdb, perturbed electric field E parallel to vdb, and perturbed magnetic field B normal to both vdb and E. So, the total asymmetric f 0 consists of non-relativistic background electrons and relativistic beam electrons:
![\eqalign{\hbox{f}_0 \lpar \hbox{p}\rpar & = {\hbox{n}_{\rm p} \over 2 \pi \hbox{m} \lpar \theta_{\rm x}^{\rm p} \theta_{\rm y}^{\rm p}\rpar ^{1/2}} \exp \left(-{\lpar \hbox{p}_{\rm x} + \hbox{p}_{\rm d}^{\rm p}\rpar ^2 \over 2 \hbox{m} \theta_{\rm x}^{\rm p}} - {\hbox{p}_{\rm y}^2 \over 2\hbox{m}\theta_{\rm y}^{\rm p}}\right)\cr & \quad + {\hbox{n}_{\rm b} \over 2 \pi \hbox{m} \gamma_{\rm b} \lpar \theta_{\rm x}^{\rm b} \theta_{\rm y}^{\rm b}\rpar ^{1/2}} \exp \left(-{\lpar \hbox{p}_{\rm x} + \hbox{p}_{\rm d}^{\rm b}\rpar ^2 \over 2\hbox{m}\gamma_{\rm b} \theta_{\rm x}^{\rm b}} - {{\rm p}_{\rm y}^2 \over 2\hbox{m} \gamma_{\rm b} \theta_{\rm y}^{\rm b}}\right).}](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021082330820-0288:S0263034608000189_eqn1.gif?pub-status=live)
with γb = (1 + p2/mc)2)1/2 and m is the electron rest mass. Here θx, θy are the temperature components parallel to the x and y directions, P d is the drift momentum, and superscripts p an b represent the plasma electron and the beam electron, respectively.
The filamentation growth rate is then qualified by the asymmetry parameters
![\hbox{A} = (\theta_{\rm x}^{\rm p} + \theta_{\rm d}^{\rm p^2}/\hbox{m}\rpar /\theta_{\rm y}^{\rm p} \quad \hbox{B} = \lpar \theta_{\rm x}^{\rm b} + \theta_{\rm d}^{\rm b^2}/\hbox{m}\gamma_{\rm b}\rpar /\theta_{\rm y}^{\rm b} \eqno \lpar 2\rpar](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021082330820-0288:S0263034608000189_eqn2.gif?pub-status=live)
denoting, respectively, to the target plasma (p) and the REB (b). These linear parameters may by refined through mode-mode coupling and quasi-linear theory (Deutsch et al., Reference Deutsch, Bret, Firpo and Fromy2005) by including the shifted temperatures
![\eqalign{& \theta_{\rm y}^{\rm p}/{\rm m} \rightarrow \theta_{\rm y}^{\rm p} \; {\rm m}+ {\omega_{\rm b}^2 \over \omega_{\rm p}^2} {\hbox{p}_{\rm d}^{\rm b^2} \over \hbox{m}^2 \gamma_{\rm b}^2} \cr & \theta_{\rm y}^{\rm b}/\hbox{m}\gamma_{\rm b} \rightarrow \theta_{\rm y}^{\rm b}/\hbox{m}\gamma_{\rm b} + {\omega_{\rm b}^{2} \over \omega_{\rm p}^{2}} {\hbox{p}_{\rm d}^{\rm b^2} \over \hbox{m}^{2}\gamma_{\rm b}^{2}}, \cr} \eqno \lpar 3\rpar](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021082330820-0288:S0263034608000189_eqn3.gif?pub-status=live)
where
![\hbox{p}_{\rm d}^{\rm b} = \hbox{mV}_{\rm b}.](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021082330820-0288:S0263034608000189_eqnU1.gif?pub-status=live)
Using these parameters in the relevant linear dispersion relations, we obtain the filamentation growth rates depicted in Figures 1 to 4 as linear (L), and quasi-linear (QL), respectively. We then emphasize relatively hot T p = (1 keV) outer DT layers and 1-MeV REB with a thermal T b = 1 keV. Maximum GR δmax occurs close to the skin depth wavelength c/ωp. δmax is also shown steadily decaying with increasing Np. Also, quasi-linear (“turbulent”) GR profiles stand beneath linear ones.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160709194743-30318-mediumThumb-S0263034608000189_fig1g.jpg?pub-status=live)
Fig. 1. WEI GR profiles with T p = T b = 1 keV. E b = 1 MeV in terms of n and n1. (a) L, (b) QL. N p/N b = 10.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160709194743-34280-mediumThumb-S0263034608000189_fig2g.jpg?pub-status=live)
Fig. 2. Same caption as Figure 1 with N p/N b = 100.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160709194743-08928-mediumThumb-S0263034608000189_fig3g.jpg?pub-status=live)
Fig. 3. Same caption as Figure 1 with N p/N b = 1000.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160709194743-15324-mediumThumb-S0263034608000189_fig4g.jpg?pub-status=live)
Fig. 4. Same caption as Figure 1 with N p/N b = 10000.
COLLISIONS CONTRIBUTION
In view of the very high densities reached by the precompressed pellet, it seems rather mandatory to pay attention to collisions between e-beam and target ions. Very often, these contributions are neglected in many standard treatments of the Weibel instability. A large consensus even asserts that collisions remain helpless against transverse electromagnetic growth rates. Nonetheless, no one has given ground to the validity of such a claim in the present extreme situation we are now considering. To clear up this point, we included in the present treatment, a Krook collision term with the right-hand side of the corresponding Vlasov equation.
We use a split Krook term υ = υp + υb. Both terms essentially document the electron-electron collision frequency
![\nu_{\rm p\comma b} = 2.91 \times 10^{-6} {\hbox{n}_{\rm p\comma b} \lpar \hbox{cm}^{-3}\rpar \over \hbox{T}_{\rm p\comma b} \lpar \hbox{eV}\rpar ^{3/2}} \ell {\rm n}\Lambda_{\rm p\comma b}\comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021082330820-0288:S0263034608000189_eqn4.gif?pub-status=live)
in terms of usual Coulomb logarithm. n p should also include a small electron-ion term υp/2.5, so one finally gets
![\nu_{\rm p} = 3.31 \times 10^{-6} {\hbox{n}_{\rm p} \lpar \hbox{cm}^{-3}\rpar \over \hbox{T}_{\rm p} \lpar \hbox{eV}\rpar ^{3/2}} \ell \hbox{n}\Lambda_{\rm p}\comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021082330820-0288:S0263034608000189_eqn5.gif?pub-status=live)
![\nu_{\rm b} = 2.91 \times 10^{-6} {\hbox{n}_{\rm b} \lpar \hbox{cm}^{-3}\rpar \over \hbox{T}_{\rm b} \lpar \hbox{eV}\rpar ^{3/2}} \ell {\rm n}\Lambda_{\rm b}.](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021082330820-0288:S0263034608000189_eqn6.gif?pub-status=live)
![\ell \hbox{n} \Lambda_{\rm p\comma b} = \ell \hbox{n}\lsqb 9 \hbox{N}_{\rm D}\rsqb \comma \;](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021082330820-0288:S0263034608000189_eqnU2.gif?pub-status=live)
with
![\hbox{N}_{\rm D} = {1\comma 72 \times 10^{9} \hbox{T}_{\rm p\comma b}^{3/2} \lpar \hbox{eV}\rpar \over {\rm n}_{\rm p\comma b}^{1/2} \lpar \hbox{cm}^{-3}\rpar }.](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021082330820-0288:S0263034608000189_eqnU3.gif?pub-status=live)
Target plasma is taken in thermal equilibrium (T e = T i) with for a DT mixture. Ion-ion collision frequency,
![\nu_{\rm p} \cong {2.9 \times 10^{-6} \over 60} {\hbox{n}_{\rm p} \lpar \hbox{cm}^{-3}\rpar \over \hbox{T}_{\rm p} \lpar \hbox{eV}\rpar ^{3/2}} \ell \hbox{n}\Lambda_{\rm p}\comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021082330820-0288:S0263034608000189_eqnU4.gif?pub-status=live)
is usually negligible. In Figures 1–4, n ≡ υp while n1 ≡ υb.
The isotropic target temperature is taken to be systematically larger than the corresponding Fermi T F accounting for nonnegligible quantum effects. We observe that as target electron density increases from 1023 up to 1025 cm−3, the discrepancy between linear L and quasi-linear QL growth rates progressively disappear. Also, collision-corrected growth rates δ have a tendency to rise above the collisionless ones for 10−2 ≤ kc/ωp ≤ 1. However, only the collision-corrected ones abruptly decay to negative (stabilization) values for 1 ≤ kc/ωp ≤ 10. So even in this rather extreme case of transverse electromagnetic instability, the beam-target collision asserts their usual damping role, as a priori expected on empirical grounds.
A WIDER PERSPECTIVE
Up to now we arbitrarily restricted attention to the filamentation instability in the given beam-plasma system. Corresponding E and k may also be permuted to produce the TS and the Weibel instabilities. More generally, it is also possible to include any nonnormal relative orientation between E and k (Bret et al., Reference Bret, Firpo and Deutsch2004, Reference Bret, Firpo and Deutsch2006). One of our most conspicuous finding is the appearance of fastest growing modes at oblique propagation (Bret et al., Reference Bret, Firpo and Deutsch2004, Reference Bret, Firpo and Deutsch2006), denoted as TS-filamentation instability (TSF).
As far as the REB-target plasma is considered, it is recalled that the target remains non-relativistic with a temperature below 10 keV, while a fiducial REB, let say a 4 MeV one (γb = 5) can feature a 500 keV relativistic energy spread. In the linear stage, filamentation is generally studied under some simplifying ab-initio transverse approximation of the dielectric tensor, so that filamentation instability is attributed to the exponential growth of unstable electromagnetic purely transverse modes (k · E = 0) with wave vector k normal to the beam (Bret et al., Reference Bret, Firpo and Deutsch2004; Deutsch et al., Reference Deutsch, Bret, Firpo and Fromy2005). It is also common to refer to this instability as Weibel instability (Bret et al., Reference Bret, Firpo and Deutsch2004, Reference Bret, Firpo and Deutsch2005a, Reference Bret, Firpo and Deutsch2005b), through the original mode (Weibel, Reference Weibel1959) has its wave vector along the beam. Figure 5 sketches the original definitions of various modes. As long as the beam is not relativistic, the largest instability it undergoes is the TS one, where the second “stream” is the return current it generates in the plasma. But in the relativistic regime, the “filamentation” growth rate eventually exceeds the TS one and is supposed to induce beam filamentation.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160709194743-16646-mediumThumb-S0263034608000189_fig5g.jpg?pub-status=live)
Fig. 5. Weibel, two-steam and filamentation modes.
In reality, the beam suffers many more instabilities at the same time. Indeed, filamentation, Weibel or TS instabilities pertain to various orientations of the wave vector and various kinds of waves (transverse or longitudinal), but in the real world, the beam-plasma system triggers every possible modes allowed by Maxwell equations with a wide range of wave vector orientation. Among all the triggered modes, the unstable ones shall start growing exponentially, while the most unstable ones shall mostly shape the beam. When it comes to knowing how the beam is eventually affected when entering the plasma, one needs therefore to answer two questions: (1) which is the most unstable mode all over, the k space for the system investigated? and (2) how shall this mode shape the beam? Following the guideline built by these two questions, we can assert that the so-called “filamentation” instability is not the fastest growing instability, even in the relativistic regime, so that it is not the answer to the first question. As for the second question, we shall see that this instability would not produce the observed effects anyway, even if it were the stronger one. For clarity, we shall keep labeling the most unstable transverse mode with wave vector normal to the beam as the “filamentation” mode, though our point is precisely that it does not filament.
We now write the growth rate of the three instabilities with the first order temperature corrections for a collisionless plasma and non-relativistic beam and plasma temperatures. As previously mentioned, there is a continuum of unstable modes so that there is not one TS mode, nor is there one filamentation mode. But when speaking about the “TS growth rate” for example, we actually mean “the maximum growth rate found on the TS profile.” With this in mind, growth rates for the TS, the filamentation, and the TS/filamentations instabilities are (Bret et al., Reference Bret and Deutsch2005c; Bret & Deutsch, Reference Bret and Deutsch2006) (α = N b/N p),
![\eqalign{\delta_{\rm TS} & = \omega_{\rm p} {\sqrt{3} \over 2^{4/3}} {\alpha^{1/3} \over \gamma_{\rm b}}\comma \; \cr \delta_{\rm TSF} & = \omega_{\rm p} {\sqrt{3} \over 2^{4/3}} \left({\alpha \over \gamma_{\rm b}}\right)^{1/3}\comma \cr \delta_{\rm F} & = \omega_{\rm p} \beta {\sqrt{\alpha} \over \gamma_{\rm b}} \left(1 - {\rho_{\rm b} \over \sqrt{\alpha \gamma_{\rm b}}}\right).}](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021082330820-0288:S0263034608000189_eqn7.gif?pub-status=live)
It is then readily seen than the largest growth rate for a relativistic beam is the TSF one in view of its 1/γb1/3 scaling. Turning to the z component of the phase velocity for these three modes, we know it is equal to the beam velocity for the TS and TSF modes (Bret et al., Reference Bret, Firpo and Deutsch2005c; Bret & Deutsch, Reference Bret and Deutsch2006) whereas, as previously mentioned, it vanishes for the Filamentation instability.
Figure 6 displays the growth rates on the TSF branch in the (k⊥, k//) plane for some zero or finite plasma thermal velocities Vtp = Y tp// = V tp⊥ and some zero or finite beam transverse thermal velocities V tb⊥. These curves clearly show that when temperatures are accounted for they act to control the instability domain, damping the small wavelengths perturbations along the filamentation direction (k// = 0), and deforming the growth rates surface so that a maximum growth rate appears for a finite oblique wave vector k M.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160709194743-34722-mediumThumb-S0263034608000189_fig6g.jpg?pub-status=live)
Fig. 6. Growth rates on the TSF branch in terms of Z = kVb/ωp with . (a) Cold beam – cold plasma. (b) Hot beam – cold plasma ad (c) Hot beam – hot plasma. Parameters are α = 0.05 and γb = 4 for (abc). Vtb⊥ = Vb/10 for (b,c) and Vtp = Vb/10 for (c).
EARLY STAGE EVOLUTION
We shall consider the out-of-equilibrium initial value dynamical problem taking place when a radically inhomogeneous electron forward current is launched into plasma and is still not current compensated (Firpo et al., Reference Firpo, Lifschitz, Lefebvre and Deutsch2006). We shall focus on this early stage where collisions may be neglected. Ions will be assumed to form a fixed neutralizing background. In order to simplify both the analysis and the numerical PIC computations, we shall consider the system to be infinite along the beam direction z. We remove any z dependence by assuming also that plasma density n pe is uniform and constant. At time t = 0, a REB is switched on in the plasma. Maxwell equations are linear and can thus be solved for all time to yield the electromagnetic fields as functions of the sources, namely beam and plasma current densities.
We consider initial beam density and velocity of the form and
with
. Let us introduce here the beam radius r b, the electron skin depth λs ≡ c/ωp, their ratio η ≡ rb/λs, and let us redefine α = nb0/npe, β0 = υb0z/c and the initial relativistic Lorentz factor
. Let us consider some initially mono-kinetic beam having density functions of the form
. This enables the study of the influence of the beam edge gradients as
.
Our analysis shows that depending on its initial radial density and velocity distribution, the shape of an electron beam propagating in plasma may evolve into a transient ring structure. This results from the natural evolution of the system and not from the usually invoked filamentation instability. If its transverse temperature is low enough, filamentation instability eventually proceeds. The observation of the ring structure is favored by sharp beam edges and not too large beam radius (compared to the electron skin depth). It is not generic which may explain the variety of experimental observations.
Figure 7 presents the longer term evolution of the transverse components of beam velocity for a mono-kinetic beam with a larger value of α (α = 0.15) and smaller γ0 (γ0 = 3). As previously, there is an initial phase, between t = 0 and 20 ωp−1, during which radial velocity grows fast and poloidal velocity remains small. For ωpet ≅ 20, we can see in the inset of Figure 7, that beam density presents a clear ring structure at its edge (Taguchi et al., Reference Taguchi, Antonsen, Liu and Mima2001). When the beam current is partially neutralized, filamentation instability starts, breaking the initial azimuthal system symmetry and producing the exponential growth of poloidal beam velocity.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160709194743-42225-mediumThumb-S0263034608000189_fig7g.jpg?pub-status=live)
Fig. 7. Early time evolution of the modulus of the radial average of the transversal beam velocity (βr = υr/c and βθ = νθ/c) in lin-log scale for α = 0.15, γo = 3, η = 10 and κ = 20. Beam and plasma transverse temperatures are 5 keV.
CONCLUDING REMARKS
Indulging now into a more global perspective, one could wonder about further refinements of growth rate calculations, beyond the present quasi-linear formalism. A preliminary step in this direction could involve a three-wave interaction mechanism. However, the numerous and very intense PIC simulations devoted to the fully nonlinear regime obviously demand a more ambitious theoretical framework. Despite that such a framework in not yet available, one could speculate that present quasi-linear improvements have only the status of higher order corrections to an initial asymptotic approximation to the complete nonlinear growth rate. In such an occurrence, it appears problematic to bridge straightforwardly the wide gap between present quasi-linear treatment and a fully nonlinear one.
Presently, it appears much more profitable to enlarge the linear approximation by superimposing three electromagnetic instabilities: Weibel, filamentation, and TS, while allowing the wave-number vector to take any orientation. Such an approach demonstrates that maximum growth rates appear on an oblique ridge between transverse and parallel directions with respect to the initial beam. It seems to us that quasi-linear corrections to this combined approach of electromagnetic instabilities might be worth while to pursue.
Moreover, it seems to us that more and still unraveled nonlinear phenomena lie in store in the very early stage of the beam-plasma interaction when counter streaming ionic flow has not yet reached a full current neutralization regime.