1. INTRODUCTION
The fast ignition scenario (FIS) for inertial confinement fusion (ICF) has prompted in recent years many theoretical (Tabak et al., Reference Tabak, Clark, Hatchett, Key, Lasinski, Snavely, Wilks, Town, Stephens, Campbell, Kodama, Mima, Tanaka, Atzeni and Freeman2005), numerical and experimental works on relativistic beam plasma instabilities (Kodama et al., Reference Kodama, Norreys, Mima, Dangor, Evans, Fujita, Kitagawa, Krushelnick, Miyakoshi, Miyanaga, Norimatsu, Rose, Shozaki, Shigemori, Sunahara, Tampo, Tanaka, Toyama, Yamanaka and Zepf2001; Roth et al., Reference Roth, Brambrink, Audebert, Blazevic, Clarke, Cobble, Cowan, Fernandez, Fuchs, Geissel, Habs, Hegelich, Karsch, Ledingham, Neely, Ruhl, Schlegel and Schreiber2005; Honrubia et al., Reference Honrubia, Antonicci and Moreno2004; Johzaki et al., Reference Johzaki, Sakagami, Nagatomo and Mima2007). This approach to ICF consists in decoupling the pellet compression from its heating, relaxing the symmetry requirement, and the overall energy of the compression laser. The pellet is thus first compressed, and ignition then left to a relativistic electron beam (REB) generated by a petawatt laser shot on the pre-compressed target. If properly tailored, the REB deposits its energy at the pellet center, igniting the fuel. On its way from the border to the pellet core, the REB is compensated by an electronic return current from the background plasma, generating a two-stream configuration well-known for its instability. Due to the density gradient of the target (the center is about 104 denser than the border), the beam-plasma interaction is collisionless near the REB emitting region, and collisional near the center. Instabilities de-focus the beam, lowering the number of hot electrons eventually reaching the core. By exciting plasma waves, instabilities also transfer energy from the beam to the plasma. While such process is deleterious at the beginning of the beam path, it is welcome at the end where the beam is supposed to stop.
For the collisionless part, it has been established on the one hand that modes propagating perpendicularly (or obliquely) to the beam are the fastest growing ones for typical FIS parameters (Deutsch, Reference Deutsch2004; Bret et al., Reference Bret, Gremillet, Bénisti and Lefebvre2008). The beam is therefore broken-up into finite length filaments, which transverse typical size to the background plasma skin-depth λp = c/ωp. On the other hand, the unstable transport in the dense collisional region reveals a qualitatively different picture: the beam is still filamented, but the typical size of the filaments is now the beam skin-depth λp = c/ωp (Gremillet et al., Reference Gremillet, Bonnaud and Amiranoff2002; Honrubia et al., Reference Honrubia, Antonicci and Moreno2004). Within the FIS context, this means filaments is about 100 times larger than in the collisionless region. The question arises to pinpoint precisely the transition from one regime to another. Assessing this gap is important from the conceptual point of view, and mandatory to describe the beam propagation in the intermediate region.
Some recent work investigated this question for the FIS (Bret et al., Reference Bret, Marín Fernández and Anfray2009). In particular, the influence of partial electronic plasma degeneracy near the pellet core was discarded, at least with respect to the unstable spectrum. Given the number of effect accounted for, this investigation was restricted to a single set of typical FIS parameters. As a consequence, the transition between the two regimes was not documented in details. The goal of the present work is to fill this gap, through a simpler theoretical model highlighting the transition threshold in terms of the main variables. The first section is devoted to the formalism implemented. A detailed study of the transition, between the two regimes then follows, preceding the discussion and conclusion.
2. FORMALISM
We thus consider a relativistic beam of density n b, velocity vb, and Lorentz factor γb = (1 − v b2/c 2)−1/2 passing through a plasma of electronic density n p. The plasma electrons are drifting with velocity vp such as n bvb= n pvp and the plasma ionic density n i is such that n i = n b + n p. The return current velocity v p = (n b/n p)v p can be considered non-relativistic since we are not interested in the fully collisionless region where n b ~ n p. Collision-wise, the electrons from the beam are supposed to be collisionless due to their large velocity (Gremillet et al., Reference Gremillet, Bonnaud and Amiranoff2002). The terms collisional/collisionless rather refers to the background electrons. Their collisionality is here characterized by the plasma electron/ion collision frequency νei.
Although the beam progresses in a density gradient, a Wentzel-Kramers-Brillouin like approximation where homogenous plasma calculations are applied locally, has been found valid for the FIS (Bret et al., Reference Bret, Firpo and Deutsch2006).
The partial degeneracy of the core electrons is neglected, since it has been found that its role on the unstable spectrum is negligible. Finally, the orientation of the perturbation wave vector k needs to be arbitrary. While this is a source of significant analytical difficulties, such framework is necessary if one wishes to capture the most unstable mode (Bret et al., Reference Bret, Firpo and Deutsch2005, Reference Bret, Firpo and Deutsch2007). As will be checked, the fastest growing modes in each regime are usually oblique. An investigation focusing on the filamentation instability with k⊥vb, would thus render improperly the beam response by bypassing the most relevant modes in this respect.
After the background plasma ions, which are assumed at rest, electrons from the beam and the plasma share the same continuity equation,
![{\partial n_j \over \partial t} + \nabla \cdot \lpar n_j {\bf v}_j \rpar =0\comma \;](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn1.gif?pub-status=live)
where the subscript j = b or p for the beam or the plasma. The Euler equation reads for the beam electrons,
![{\partial {\bf p}_b \over \partial t}+ \lpar {\bf v}_b \cdot \nabla \rpar {\bf p}_b = -q \left({\bf E} + {{\bf v}_b\times {\bf B} \over c}\right)-{\nabla P_b \over n_b} \comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn2.gif?pub-status=live)
and for the plasma ones,
![{\partial {\bf v}_p\over \partial t}+ \lpar {\bf v}_p\cdot\nabla \rpar {\bf v}_p= -{q\over m}\left({\bf E} + {{\bf v}_p\times{\bf B}\over c}\right)-\nu {\bf v}_p -{\nabla P_p\over n_p}.](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn3.gif?pub-status=live)
The beam equation is thus collisionless and relativistic, while the plasma is non-relativistic and collisional. The pressure terms are expressed in terms of the temperatures through ∇P j = 3k BT j ∇ n j, where k B is the Boltzmann constant. Such an adiabatic treatment demands sub-relativistic temperatures, a requirement stronger for the beam than for the plasma (Siambis, Reference Siambis1979; Pegoraro & Porcelli, Reference Pegoraro and Porcelli1984). Though lengthy, the derivation of the dispersion equation is quite standard (Bret et al., Reference Bret, Firpo and Deutsch2005) and eventually reads,
![\det {\cal T} = 0 \comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn4.gif?pub-status=live)
where the tensor is reported in Section 6, and expressed in terms of the dimensionless variables,
![{\rm \alpha}={n_{b}\over n_{\,p}} \comma \; {\bf Z}={{\bf k}v_b\over {\rm \omega}_p} \comma \; {\rm \beta}={v_b\over c} \comma \; {\rm \tau}={\nu\over {\rm \omega}_p} \comma \; \rho_{\,j}=\sqrt{{3 k_BT_j\over m v_b^2}}.](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn5.gif?pub-status=live)
Calculations have been conducted aligning the beam velocity vb with the z axis, and considering k = (k x, 0, k z). Components k z and Z z are therefore the parallel ones, while k x and Z x are the perpendicular ones.
3. THE COLLISIONAL/COLLISIONLESS TRANSITION
Figure 1 shows a typical growth-rate map arising from the numerical resolution of the dispersion equation. Modes localized around Z x ~ Z z ~ 1 are collisionless ones and produce filaments of transverse size ~c/ωp. Note their oblique location, impossible to capture if restricting the exploration to the main axis. Unstable modes at small Z are collisional ones, as one can check they vanish when setting ν = 0. The full spectrum is here clearly governed by these collisional modes. The fastest growing mode is found for Z z = 0.014 and Z x = 0.11, producing much larger filaments than the collisional modes. The simple relation between the beam and plasma skin-depths,
![{\rm \lambda}_b={{\rm \lambda}_p\over \sqrt{\rm \alpha}} \comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn6.gif?pub-status=live)
shows that their size fits here the beam skin-depth, as expected when dealing with resistive filamentation (Gremillet et al., Reference Gremillet, Bonnaud and Amiranoff2002).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629041610-47361-mediumThumb-S0263034610000480_fig1g.jpg?pub-status=live)
Fig. 1. (Color online) Typical growth-rate map in terms of (Z x, Z z). Parameters feature α = 10−2, ρp = 4.2 × 10−2, ρp = 0.42, τ = 0.3, and γb = 4.
The unstable spectrum is thus clearly divided into two parts: the “lower” collisional spectrum, and the “upper” collisionless one. Our goal at this juncture is twofold: on the one hand, studying the evolution of the fastest growing mode (and its growth rate) of each part and on the other hand, documenting the transition between the two regimes. In view of the vast numbers of free parameters, we focus on the (τ, α) mapping, choosing for the other variables some FIS relevant values. We thus explore the parameters space,
![\eqalignno{&\rm \alpha\in \lsqb 0 \comma 10^{-1}\rsqb \comma \; \rm \tau\in \lsqb 0 \comma 0.5\rsqb \comma \; \rm \gamma_b=4 \comma \cr & \quad \rm \rho_p=4.2\times 10^{-2} \lpar T_p=1\, {\rm keV} \rpar \;\rm \rho_b=0.42 \lpar T_b=100\, {\rm keV} \rpar. &}](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn7.gif?pub-status=live)
The fastest growing collisional mode Zm and its growth rate have been calculated numerically over the parameters window specified above. Regardless of (τ, α), we find Zm ~ (1,1) so that this part of the spectrum definitely generate skin-depth size structures. Figure 2 pictures the growth rate of these modes in terms of τ for α = 10−5…10−1. Regarding the alpha dependance, a ακ with κ ~ 1 is retrieved, reminiscent of the exact kinetic scaling (Bret et al., Reference Bret, Gremillet and Bénisti2010).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629041828-47119-mediumThumb-S0263034610000480_fig2g.jpg?pub-status=live)
Fig. 2. Growth-rate of the collisionless modes (upper spectrum) in terms of τ for α = 10−5…10−1.
Figure 3 now sketches the largest collisional growth rate in terms of (α, τ). We observe some quasi-constant τ and α scalings, which are numerically found very close to τ1/3 and α2/3.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629041947-70816-mediumThumb-S0263034610000480_fig3g.jpg?pub-status=live)
Fig. 3. (Color online) Growth-rate of the collisional modes (lower spectrum) in terms of τ and α.
The τ trends identified makes it clear that beyond a given collisionality threshold, collisional modes must surpass the collisionless ones. The resulting partition of the (τ,α) domain is pictured in Figure 4, where the beam trajectory from the pellet border to the core is superimposed. Instability wise, the beam clearly starts from the collisionless region to end up in the collisional one. The upper-spectrum is thus relevant at the beginning while the lower one is more important at the end.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629042309-78535-mediumThumb-S0263034610000480_fig4g.jpg?pub-status=live)
Fig. 4. (Color online) Portions of the (τ,α) plane governed by each kind of modes. The red path pictures the beam trajectory from the pellet border to the core. The blue path accounts for the relation between α and τ according to Eqs. (8) and (9). A more detailed analysis (Ren et al., Reference Ren, Tzoufras, Tonge, Mori, Tsung, Fiore, Fonseca, Silva, Adam and Heron2006) gives the green trajectory. In the upper-right corner domain, limited by the gray line, the fastest collisional modes have k ∥ = 0.
Note that the collision parameter τ and the density ratio parameter α are not independent quantities. If restricting the former to electron-ion collisions, the collision frequency νei reads (Huba, Reference Huba2004),
![\nu_{ei}=2.91\times 10^{-6}{n_i \lpar {\rm cm}^{-3} \rpar \over T_p^{3/2} \lpar {\rm eV} \rpar }\ln \Lambda \comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn8.gif?pub-status=live)
with
![\ln \Lambda=23.12 - \ln\left[{n_p^{1/2} \lpar {\rm cm}^{-3} \rpar \over T_p \lpar {\rm eV} \rpar }\right].](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn9.gif?pub-status=live)
The beam path in the (τ, α) plane, as calculated according to the equations above and with T p = 1 keV being pictured in blue in Figure 4. When accounting for a plasma temperature map extracted from a simulation of a pre-compressed target (Ren et al., Reference Ren, Tzoufras, Tonge, Mori, Tsung, Fiore, Fonseca, Silva, Adam and Heron2006), the path recovered is indicated in green. At any rate, the collisional regime is relevant near the core.
4. CHARACTERIZATION OF COLLISIONAL MODES
We finally turn to the most unstable wave-vector analysis. Our goal is mainly to check the size of the structures generated. To this extent, Figure 5 pictures the parallel (Z z) and perpendicular (Z x) components of the most unstable wave-vector. To start with, we find here another kind of transition: for most of the (τ, α) parameters, Zm is oblique, that is, Z z ≠ 0. However, Z z vanishes for τ > 0.5 and α > 10−2 (the shape of the domain is more involved, as evidenced in Fig. 4), a regime that can be denoted as “resistive filamentation” since filamentation modes with k⊥vb are the fastest ones here. This region is represented in the upper-right corner of Figure 4.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160629042609-79565-mediumThumb-S0263034610000480_fig5g.jpg?pub-status=live)
Fig. 5. (Color online) Perpendicular (Z x, left) and parallel (Z z, right) components of the most unstable wave-vector.
The resistive filamentation regime is equally singled out on the perpendicular (Z x) component with a specific scaling. Further numerical analysis shows that Z x scales like α1/2 when Z z = 0, and ασ with σ ~ 1/3 otherwise. Interestingly, the α1/2 scaling is exactly what would be expected according to the beam skin-depth instead of the plasma one. In the oblique collisional regime, we still witness an increase of the filaments size when α decreases, but the scaling is too slow to keep up with the beam skin-depth.
5. CONCLUSION
A model of the REB transport in FIS context has been build to study the bridge between collisionless and collisional instabilities. Considering a relativistic, non-collisional diluted beam, passing through a non-relativistic and collisional return current, two kinds of unstable modes spontaneously arise from the resolution of the dispersion equation. For wave vectors k such as k ~ λp, collisionless unstable modes are triggered and generate plasma skin-depth size structures. At much lower k, collisional modes excited at finite collision frequency νei generate beam skin-depth size structures. When νei increases, collisionless instabilities are mitigated while collisional ones are boosted. As a result, a transition from one regime to the other occurs for a critical collision frequency. This threshold has been determined in terms of the beam-to-plasma density ratio, evidencing the collisional nature of the beam/target-core interaction.
Though the dielectric tensor and the dispersion equation are analytically known, it has not been possible so far to derive analytical expressions of the various quantities involved. Even when focusing on the filamentation modes with k⊥vb, some approximations can yield a simpler dispersion equation, yet still analytically not extractable. The analytical derivation of the various scalings mentioned in this paper remains postponed for the future. Finally, a kinetic theory of the process would be interesting in order to confirm the collisionless/collisional transition and refine its description.
6. TENSOR ELEMENTS
The tensor is symmetric and reads,
![{ \cal T}=\left(\matrix{{\cal T}_{11} & 0 & {\cal T}_{13} \cr 0 & {\cal T}_{22} & 0 \cr {\cal T}_{31} & 0 & {\cal T}_{33}}\right)\comma \;](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn10.gif?pub-status=live)
where,
![\eqalignno{&{\cal T}_{11}=x^2- {Z_z^2\over {\rm \beta} ^2} + {- \lpar x-Z_z \rpar ^2 {\rm \alpha}{\rm \gamma}_b^3+Z_z^2 {\rm \alpha}{\rm \rho}_b^2\over \lpar x-Z_z \rpar ^2 {\rm \gamma}_b^4-{\rm \gamma}_b \lpar Z_z^2+Z_x^2 {\rm \gamma}_b^2 \rpar {\rm \rho}_b^2} \cr & \quad -{Z_z^2 \lpar x+Z_z {\rm \alpha} \rpar \over \lpar Z_x^2+Z_z^2 \rpar \lpar x+Z_z {\rm \alpha} +i {\rm \tau} \rpar } \cr &\quad +{Z_x^2 \lpar x+Z_z {\rm \alpha} \rpar ^2\over \lpar Z_x^2+Z_z^2 \rpar \left(\matrix{x^2+2 x Z_z {\rm \alpha} +Z_z^2 {\rm \alpha}^2 -Z_x^2 {\rm \rho}_p^2-Z_z^2 {\rm \rho}_p^2+i \lpar x+Z_z {\rm \alpha} \rpar {\rm \tau}} \right)} \comma \cr &&}](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn11.gif?pub-status=live)
![{\cal T}_{22}=x^2-{Z_x^2+Z_z^2\over {\rm \beta} ^2}-{{\rm \alpha} \over {\rm \gamma}_b}-{x+Z_z {\rm \alpha} \over x+Z_z {\rm \alpha} +i {\rm \tau}}\comma](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn12.gif?pub-status=live)
![\eqalignno{&{\cal T}_{33}= {P\over {\rm \beta} ^2 {\rm \gamma}_b \lpar \lpar x-Z_z \rpar ^2 {\rm \gamma}_b^3 - \lpar Z_z^2+Z_x^2 {\rm \gamma}_b^2 \rpar {\rm \rho}_b^2 \rpar } -{Z_x^2 \lpar x+Z_z {\rm \alpha} \rpar \over \lpar Z_x^2+Z_z^2 \rpar \lpar x+Z_z {\rm \alpha} +i {\rm \tau} \rpar } \cr &\quad -{ \lpar x Z_z-Z_x^2 {\rm \alpha} \rpar ^2 \over \matrix{\lpar Z_x^2+Z_z^2 \rpar \lpar \lpar x+Z_z {\rm \alpha} \rpar ^2 - \lpar Z_x^2+Z_z^2 \rpar {\rm \rho}_p^2+i \lpar x+Z_z {\rm \alpha} \rpar {\rm \tau} \rpar }} \comma &}](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn13.gif?pub-status=live)
where,
![\eqalignno{P & = 2 x Z_x^2 Z_z {\rm \gamma}_b^4+x^4 {\rm \beta} ^2 {\rm \gamma}_b^4-2 x^3 Z_z {\rm \beta} ^2 {\rm \gamma}_b^4 \cr & \quad +Z_x^2 \lpar -{\rm \gamma}_b^3 \lpar {\rm \alpha}{\rm \beta} ^2+Z_z^2 {\rm \gamma}_b \rpar + \lpar {\rm \alpha}{\rm \beta} ^2+Z_z^2 {\rm \gamma}_b+Z_x^2 {\rm \gamma}_b^3 \rpar {\rm \rho}_b^2 \rpar \cr & \quad -x^2 {\rm \gamma}_b \lpar Z_x^2 {\rm \gamma}_b^3+{\rm \beta} ^2 \lpar {\rm \alpha} +Z_x^2 {\rm \gamma}_b^2 {\rm \rho}_b^2+Z_z^2 \lpar -{\rm \gamma}_b^3+{\rm \rho}_b^2 \rpar \rpar \rpar \comma &}](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn14.gif?pub-status=live)
and,
![\eqalign{&{\cal T}_{13}={\cal T}_{31}\cr&=Z_x {\matrix{\lpar x-Z_z\rpar {\rm \gamma}_b^3 \lpar -{\rm \alpha}{\rm \beta} ^2+\lpar x-Z_z\rpar Z_z {\rm \gamma}_b\rpar -Z_z \lpar {\rm \alpha}{\rm \beta} ^2 +Z_z^2 {\rm \gamma}_b +Z_x^2 {\rm \gamma}_b^3\rpar {\rm \rho}_b^2} \over {\rm \beta} ^2 \lpar \lpar x-Z_z\rpar ^2 {\rm \gamma}_b^4-{\rm \gamma}_b \lpar Z_z^2+Z_x^2 {\rm \gamma}_b^2\rpar {\rm \rho}_b^2\rpar } \cr &\quad + {Z_x Z_z \lpar x+Z_z {\rm \alpha} \rpar \over \lpar Z_x^2+Z_z^2\rpar \lpar x+Z_z {\rm \alpha} +i {\rm \tau} \rpar } \cr &\quad + {Z_x \lpar -x Z_z+Z_x^2 {\rm \alpha} \rpar \lpar x+Z_z {\rm \alpha} \rpar \over \lpar Z_x^2+Z_z^2\rpar \left(\matrix{x^2+2 x Z_z {\rm \alpha} +Z_z^2 {\rm \alpha}^2 -Z_x^2 {\rm \rho}_p^2 - Z_z^2 {\rm \rho}_p^2 +i \lpar x+Z_z {\rm \alpha} \rpar {\rm \tau}} \right)}.}](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20151021111338924-0040:S0263034610000480_eqn15.gif?pub-status=live)
ACKNOWLEDGMENTS
This work has been achieved under projects ENE2009-09276 of the Spanish Ministerio de Educación y Ciencia and PAI08-0182-3162 of the Consejería de Educación y Ciencia de la Junta de Comunidades de Castilla-La Mancha.