1. Introduction
Plasma accelerators that rely on high-amplitude plasma wakefields are promising several orders of magnitude higher acceleration gradients than in conventional radio-frequency accelerators (Malka et al. Reference Malka, Faure, Gauduel, Lefebvre, Rousse and Phuoc2008; Esarey, Schroeder & Leemans Reference Esarey, Schroeder and Leemans2009). The two main types of plasma accelerators are laser–wakefield accelerators (LWFAs) based on driving a plasma wake with a short intense laser pulse (Tajima & Dawson Reference Tajima and Dawson1979) and plasma–wakefield accelerators (PWFAs) based on using a short high-current particle bunch as a driver (Chen et al. Reference Chen, Dawson, Huff and Katsouleas1985; Rosenzweig et al. Reference Rosenzweig, Breizman, Katsouleas and Su1991). LWFAs have demonstrated rapid growth both in terms of the energy of accelerated electrons (reaching energies of 8 GeV at a distance of 20 cm Gonsalves et al. Reference Gonsalves, Nakamura, Daniels, Benedetti, Pieronek, de Raadt, Steinke, Bin, Bulanov and van Tilborg2019) and stability of beam properties (Faure et al. Reference Faure, Rechatin, Norlin, Lifschitz, Glinec and Malka2006). For PWFA, experimental studies have demonstrated high-gradient acceleration (Blumenfeld et al. Reference Blumenfeld, Clayton, Decker, Hogan, Huang, Ischebeck, Iverson, Joshi, Katsouleas and Kirby2007) and efficient energy transfer from the driver to the witness (Litos et al. Reference Litos, Adli, An, Clarke, Clayton, Corde, Delahaye, England, Fisher and Frederico2014; Lindstrøm et al. Reference Lindstrøm, Garland, Schröder, Boulton, Boyle, Chappell, D'Arcy, Gonzalez, Knetsch and Libov2021).
Despite some advantages such as the long dephasing length and a more stable wakefield structure, PWFAs saw comparatively less development than LWFAs because the sources of short (comparable to plasma wavelength) electron bunches with a high enough current to excite a nonlinear wake were not widely available. Recently, the concept of hybrid LWFA–PWFA multistaged plasma accelerators (Hidding et al. Reference Hidding, Königstein, Osterholz, Karsch, Willi and Pretzler2010; Kurz et al. Reference Kurz, Heinemann, Gilljohann, Chang, Couperus Cabadağ, Debus, Kononenko, Pausch, Schöbel and Assmann2021; Foerster et al. Reference Foerster, Döpp, Haberstroh, Grafenstein, Campbell, Chang, Corde, Couperus Cabadağ, Debus and Gilljohann2022) based on the idea of using electron bunches from the first LWFA stage to drive a second PWFA stage has started gaining popularity, broadening the possibilities for experimental PWFA studies. LWFA-produced electron bunches are very short and high current, so they naturally have excellent parameters to drive a highly nonlinear (blowout) plasma wake. The density-downramp injection technique also proved to be an effective way of injecting electrons in the PWFA stage (Grebenyuk et al. Reference Grebenyuk, Martinez de la Ossa, Mehrling and Osterhoff2014; Martinez de la Ossa et al. Reference Martinez de la Ossa, Hu, Streeter, Mehrling, Kononenko, Sheeran and Osterhoff2017; Xu et al. Reference Xu, Li, An, Dalichaouch, Yu, Lu, Joshi and Mori2017; Zhang et al. Reference Zhang, Huang, Marsh, Xu, Li, Hogan, Yakimenko, Corde, Mori and Joshi2019; Couperus Cabadağ et al. Reference Couperus Cabadağ, Pausch, Schöbel, Bussmann, Chang, Corde, Debus, Ding, Döpp and Foerster2021). Even though the electron-driven second stage cannot significantly surpass the performance of the first LWFA stage in terms of the total energy of the accelerated electrons, due to its stable nature it can serve as a ‘quality booster’ by generating bunches with improved properties (Martinez de la Ossa et al. Reference Martinez de la Ossa, Hu, Streeter, Mehrling, Kononenko, Sheeran and Osterhoff2017; Foerster et al. Reference Foerster, Döpp, Haberstroh, Grafenstein, Campbell, Chang, Corde, Couperus Cabadağ, Debus and Gilljohann2022). However, although relatively high energy transfer efficiency is achieved and high-quality beams with hundreds of megaelectronvolts are predicted, the energy stability for such accelerator in previous research is still comparable to the single-stage LWFAs (Foerster et al. Reference Foerster, Döpp, Haberstroh, Grafenstein, Campbell, Chang, Corde, Couperus Cabadağ, Debus and Gilljohann2022), and more studies are required to understand the parameter dependence and physics behind this stage.
For PWFA to be efficient, one of the important steps is to optimise the driver-to-witness energy transfer efficiency. The beam currents of both the driver and the witness beams play decisive roles in optimising the energy transfer efficiency. Few studies so far have systematically focused on controlling the beam currents produced in plasma-based accelerators. The production in LWFA of the current profile finely tuned for the use in the second stage was studied in Hue et al. (Reference Hue, Wan, Levine and Malka2023), but studies for PWFAs are lacking.
In this paper, we focus on the performance of a PWFA based on density-downramp injection of electrons. One important step in calculating and optimising the energy transfer efficiency is to understand the beam dynamics. Previous research described such important phenomena as the hosing instability experienced at the tail of the bunch (Huang et al. Reference Huang, Lu, Zhou, Clayton, Joshi, Mori, Muggli, Deng, Oz and Katsouleas2007), beam head erosion (Zhou et al. Reference Zhou, Clayton, Huang, Joshi, Lu, Marsh, Mori, Katsouleas, Muggli and Oz2007; Li et al. Reference Li, Adli, England, Frederico, Gessner, Hogan, Litos, Walz, Muggli and An2012; An et al. Reference An, Zhou, Vafaei-Najafabadi, Marsh, Clayton, Joshi, Mori, Lu, Adli and Corde2013), energy depletion (Muggli et al. Reference Muggli, Blumenfeld, Clayton, Decker, Hogan, Huang, Ischebeck, Iverson, Joshi and Katsouleas2010) and the transformer ratio (Blumenfeld et al. Reference Blumenfeld, Clayton, Decker, Hogan, Huang, Ischebeck, Iverson, Joshi, Katsouleas and Kirby2010). Yet, driver parameters used in these studies corresponded mostly to electron beams produced by conventional linear accelerators, which are very different from beams produced in LWFA stage for hybrid accelerators. In this article, the dynamics of flattop–current electron drivers with a total beam charge of 100–500 pC, 100–300 MeV energy typical for the LWFA-produced bunches (Hue et al. Reference Hue, Wan, Levine and Malka2023) is studied. In § 2, the stability of the wakefield and driver-to-plasma energy transfer efficiency are discussed. In § 3, we numerically study density-downramp injection in the PWFA stage and the dependence of the current of the injected electron bunch on the parameters of the driver, the plasma and the downramp. Despite having such a multidimensional parameter space, we demonstrate that the injected current is determined mostly by one parameter, the effective current $J_{\textrm {eff}} = J_b (k_p \xi _b)^{2/3}$
, which combines both the properties of the driver (its current $J_b$
and length $\xi _b$
) and the plasma (wavenumber $k_p$
). The dependence of the witness current on this parameter is linear. We also show the limitation of this scaling for longer driver bunches which cannot efficiently excite a nonlinear wake.
2. Driver stability and driver-to-plasma efficiency
In this section, we study the evolution and propagation of an electron beam in PWFA assuming fully preionised plasma. For sufficiently short, tightly focused and high-current beams, the wakefield is excited in the strongly nonlinear (‘bubble’ or ‘blowout’) regime (Rosenzweig et al. Reference Rosenzweig, Breizman, Katsouleas and Su1991; Pukhov & Meyer-ter-Vehn Reference Pukhov and Meyer-ter-Vehn2002; Lotov Reference Lotov2004), when a cavity (a bubble) devoid of plasma electrons is formed behind the driver. The excitation of the bubble by an electron beam can be self-consistently described by a model by Golovanov et al. (Reference Golovanov, Kostyukov, Reichwein, Thomas and Pukhov2021) which can be used to calculate the shape of the bubble as well as the distribution of all the fields in it based solely on the charge density distribution of the driver. As this model is based on the relativistic limit of the theory by Lu et al. (Reference Lu, Huang, Zhou, Mori and Katsouleas2006), it is strictly valid only in the case when the transverse size of the bubble $R_{\textrm {bub}}$
is large in terms of plasma units, $k_p R_{\textrm {bub}} \gg 1$
, which is not necessarily the case for the parameters considered in the paper. A more accurate model for comparatively small bubbles ($k_p R_{\textrm {bub}} \sim 1$
) was recently proposed in Golovanov et al. (Reference Golovanov, Kostyukov, Pukhov and Malka2023), but it lacks simple analytical scalings. Theoretical models also cannot fully describe the driver evolution and self-injection, and thus cannot completely replace numerical simulations. Still, we can use the predictions of the models to compare with the simulation results.
To study the evolution of the driver in PWFA, we perform numerical simulations of the beam–plasma interaction using the 3D quasistatic particle-in-cell (PIC) code QuickPIC (Huang et al. Reference Huang, Decyk, Ren, Zhou, Lu, Mori, Cooley, Antonsen and Katsouleas2006). The beam has the mean energy of 250 MeV and the charge of 137 pC (corresponding to the total energy of 17 mJ) with a flattop longitudinal current profile with the length $\xi _b = 13.4\,\mathrm {\mu } {\rm m}$
(or the duration 45 fs corresponding to the current of 3.1 kA) and a Gaussian transverse profile $\exp (-r^2/2\sigma _r^2)$
with radius $\sigma _r = 0.52\,\mathrm {\mu } {\rm m}$
and normalised emittance of $\epsilon _{x,y} = 0.42\,{\rm \pi} {\rm mm}\,{\rm mrad}$
defined as $\epsilon _{x} = (mc)^{-1}\sqrt {\left \langle x^2 \right \rangle \left \langle p_x^2 \right \rangle - \left \langle x p_x \right \rangle ^2}$
. The peak number density of the beam $n_b = 3.75\times 10^{19}\,{\rm cm}^{-3}$
. The chosen parameters correspond to typical electron beams generated by the density-downramp injection in the first LWFA stage and are taken from simulations in Hue et al. (Reference Hue, Wan, Levine and Malka2023). In fact, as will be briefly explained in the next section, the tunability of the LWFA-produced driver parameters is fairly limited.
The plasma density $n_0$
is chosen to be much lower than the density of the beam $n_b$
, so that the driver excites a highly nonlinear plasma wakefield or a bubble (see figure 1). In the simulation, we use a box with the size of $6 k_p^{-1} \times 6 k_p^{-1} \times 7 k_p^{-1}$
and $512 \times 512 \times 256$
cells (the beam propagates along $z$
), where $k_p = \omega _p / c$
is the plasma wavenumber corresponding to the plasma density $n_0$
, $\omega _p = (n_0 e^2/m \varepsilon _0)^{1/2}$
and $\varepsilon _0$
is the vacuum permittivity. The number of particles per cell for the background plasma is $4$
, and for the beam the total number of macroparticles is $2^{21}$
. The simulation timestep is equal to $5 \omega _p^{-1}$
.
Figure 1. Electron density $n_{e}$
distributions and the longitudinal electric fields $E_z$
on the axis in a wake driven by an electron driver at different propagation distances. The dashed lines show the initial distribution of the accelerating field at 0 mm. For the driver, the colour shows the energy of particles. The dotted lines at 0 mm show the analytical solution according to the model in Golovanov et al. (Reference Golovanov, Kostyukov, Pukhov and Malka2023). Vertical dotted lines show the position of the maximum decelerating field $E_z$
. Plasma density $n_0 = 3.125\times 10^{17}\,{\rm cm}^{-3}$
.
The transverse emittance of the beam is chosen to be close to matched to the plasma density to prevent significant changes of the transverse size (the matched plasma density for the chosen beam parameters is $2.74\times 10^{17}\,{\rm cm}^{-3}$
). In this wakefield, the driver generally experiences two forces: the focusing force from the ion column (linear in the distance $r$
from the axis inside the bubble and leading to betatron oscillations of the electrons of the driver) and the decelerating longitudinal force. We observe that the accelerating structure remains highly stable before rapidly collapsing when electrons of the driver start dephasing due to deceleration to very low energies. This inherent stability owes to the fact that the changing transverse distribution of the almost matched driver does not influence the structure of the bubble, whereas the longitudinal velocity of the ultrarelativistic particles of the bunch stays almost equal to the speed of light until some of the electrons are decelerated to subrelativistic energies, corresponding to the moment of collapse (compare the first three distances in figure 1 before the collapse to the fourth after the collapse began).
The propagation length corresponding to this collapse is much shorter than the distance of particle loss due to head erosion or beam hosing instability, so it effectively determines the PWFA stage length. Head erosion usually happens due to the initial emittance of the driver in the region at the head of the bunch where the gas is not yet ionised by its fields and, thus, the head propagates effectively in vacuum and experiences no focusing from the plasma (Zhou et al. Reference Zhou, Clayton, Huang, Joshi, Lu, Marsh, Mori, Katsouleas, Muggli and Oz2007; Li et al. Reference Li, Adli, England, Frederico, Gessner, Hogan, Litos, Walz, Muggli and An2012; An et al. Reference An, Zhou, Vafaei-Najafabadi, Marsh, Clayton, Joshi, Mori, Lu, Adli and Corde2013). The Coulomb self-force is proportional to $\gamma ^{-2}$
and can be usually ignored for ultrarelativistic bunches. When the plasma is preionised, head erosion can still happen because the wakefield providing the focusing force is not yet strong enough at the head of the bunch, but its effect is much smaller than in the initially neutral gas (An et al. Reference An, Zhou, Vafaei-Najafabadi, Marsh, Clayton, Joshi, Mori, Lu, Adli and Corde2013). As the head erodes, it still continues to drive a lower-amplitude wake behind it, which makes the spread of this erosion to other parts of the bunch very slow. As can be seen from figure 1, it leads only to a small shift in the phase of the nonlinear wake over the propagation distance and cannot significantly affect the acceleration process.
Another effect which is believed to limit the length of PWFAs is the hosing instability which should lead to the exponential growth in oscillations of the beam centroid and the corresponding growth of oscillations of the bubble (Huang et al. Reference Huang, Lu, Zhou, Clayton, Joshi, Mori, Muggli, Deng, Oz and Katsouleas2007). However, recent research suggests that energy depletion of the driver as well as its energy spread or energy chirp can efficiently suppress and saturate this instability due to the dephasing of betatron oscillations (Mehrling et al. Reference Mehrling, Fonseca, de la Ossa and Vieira2017, Reference Mehrling, Fonseca, Martinez de la Ossa and Vieira2019). In our simulations, hosing is not taken into account as the initial charge density distribution is ideally symmetric.
Thus, energy depletion of the driver is the main mechanism which determines the acceleration length. It leads to electrons reaching subrelativistic energies, after which significant portion of the electrons of the driver are quickly lost and the accelerating structure collapses (see figure 1 at 14.4 mm). For a bunch with no energy chirp, this process happens at the point of the peak decelerating electric field the value of which can be estimated from the solution based on the model in Golovanov et al. (Reference Golovanov, Kostyukov, Pukhov and Malka2023) (shown with dotted lines in figure 1 at 0 mm). The length at which the accelerating structure collapse happens can, thus, be estimated as the length at which an electron with the kinetic energy $K \approx \gamma m c^2$
is decelerated to zero energy in the field $E_{\textrm {max}}$
,
where $\gamma$
is the Lorentz factor, $m$
is the electron mass, $e > 0$
is the elementary charge and $c$
is the speed of light. The comparison to values observed in the simulations (see table 1) provides a fairly good estimate for the collapse length.
Table 1. The collapse length of the electron driver $L_{\textrm {col}}$
observed in quasistatic PIC simulations and estimated from (2.1) as well as the efficiency $\eta$
calculated according to (2.2) and estimated based on the model from Golovanov et al. (Reference Golovanov, Kostyukov, Pukhov and Malka2023) for different plasma densities $n_0$
.
When the collapse of the accelerating structure begins, the driver still has part of its energy left (as decelerating field is non-uniform, and different parts of the driver experience different field values), which limits driver-to-plasma energy transfer efficiency defined as the percentage of the bunch initial energy spent on generating the wakefield by the moment the collapse begins,
where averaging is performed over the bunch particles. Assuming that the bunch is monoenergetic, it can also be estimated as $\eta \approx \left \langle E_z \right \rangle /E_{{\textrm {max}}}$
, so the efficiency mostly reflects how uniform the distribution of the decelerating electric field inside the driver is. The comparison between the actual efficiency observed in simulations to the estimated efficiency based on the field distribution $E_z$
according to the model by Golovanov et al. (Reference Golovanov, Kostyukov, Pukhov and Malka2023) is given in table 1. The efficiency changes depending on the plasma density, so it is required to carefully choose the plasma density for the second PWFA stage in order to increase the energy transfer efficiency. For a very high plasma density, the driving bunch can become too long compared with the plasma wavelength (determined by the dimensionless value $k_p \xi _b$
in table 1), leading to its tail being accelerated in the created bubble and a significant drop in the driver-to-plasma energy transfer efficiency due to non-uniformity of the field distribution. The efficiency can, thus, be increased by lowering the plasma density. In very low-density plasmas, the efficiency also starts to slightly drop due to the non-uniformity of the field at the front of a short driver. In addition, lower densities can be less desirable due to lower acceleration gradients, which increase the size of the accelerator. Therefore, there is an optimal range of plasma densities at which the beam-to-plasma efficiency is high. As results of numerical simulations show (table 1), for the considered electron driver, the optimal plasma density of the second PWFA stage is around $10^{18}\,{\rm cm}^{-3}$
, corresponding to the length of the driver in plasma units $k_p \xi _b \sim 3$
.
3. Downramp injection
In this section, we study the dependence of the current profile of the witness electron beam produced by density-downramp injection on the driver and plasma parameters. A recent study on downramp injection for LWFA (Hue et al. Reference Hue, Wan, Levine and Malka2023) reports that the injected beam current is strongly influenced by the laser intensity and, consecutively, the wakefield strength. Unlike for the laser driven counterpart, the electron driver can be considered non-evolving during the beam injection process in the PWFA stage since the beam evolution scale is the period of betatron oscillations $\sqrt {2 \gamma } \lambda _p$
, usually of millimetre to centimetre scale, much longer than the downramp length, usually in the range of hundreds of micrometres or less. As a result, the location of the downramp does not play such a crucial role as in the LWFA case.
3.1. The influence of downramp steepness
Now we consider the influence of downramp steepness on the parameters of the injected beam. As quasistatic codes such as QuickPIC used in § 2 cannot self-consistently describe the injection of particles into the wakefield, we perform PIC simulations with FBPIC (Lehe et al. Reference Lehe, Kirchen, Andriyash, Godfrey and Vay2016) which is a spectral quasi-3-D PIC code with angular mode decomposition. In the simulations, we use a simulation box of the $80\,\mathrm {\mu }$
m size with 2048 cells in the longitudinal $z$
direction and the $50\,\mathrm {\mu }$
m size with 640 cells in the radial $r$
direction. The number of the angular modes is equal to 2, and the total number of particles per cell is 45 (5 in the angular direction and $3 \times 3$
in the $zr$
plane). The timestep is equal to the longitudinal grid cell size, ${\rm \Delta} t = {\rm \Delta} z/c$
, and simulations are performed in a moving window with the velocity of $c$
.
The downramp is modelled as a linear change of plasma density from $n_0$
to $n_0/2$
of length $L$
. For a given downramp steepness, particles are injected only when their energy is high enough, the required energy grows with the decrease of the downramp steepness (Xu et al. Reference Xu, Li, An, Dalichaouch, Yu, Lu, Joshi and Mori2017). The energy of the particles depends only on the nonlinearity of the created bubble (which will be quantified later) irrespective of other driver properties.
Four group of simulations are presented in figure 2 that explain the role that the nonlinearity of the bubble (which depends on the density ratio $n_b/n_0$
between the driver and the plasma for a fixed driver shape and size) and $L$
play on the witness injection. The injected beam currents are plotted and the corresponding wakefield structures are illustrated in the insets of figure 2(c,g). As figure 2(b) demonstrates, for a fixed downramp length $L$
, the injection occurs only when the bubble nonlinearity is strong enough. In addition, once the condition for injection is reached, the witness beam current is barely influenced by the downramp steepness for the same bubble strength (the same driver), as shown in figure 2(c). Similar to the results in Ekerfelt et al. (Reference Ekerfelt, Hansson, González, Davoine and Lundh2017), Silva et al. (Reference Silva, Helm, Vieira, Fonseca and Silva2019) and Hue et al. (Reference Hue, Wan, Levine and Malka2023), a current peak is observed at the head of the witness bunch. This happens due to the nonlinear phase mixing caused by a sharp plasma density transition when the downramp begins, and the peak can be mitigated with a smooth density transition at the beginning of the downramp, usually found in experimental cases. Behind the head of the bunch with respect to the comoving coordinate $\xi = z - ct$
, the injected current stabilises at a constant value. Therefore, in the following considerations, we define the value of the injected witness current $J_w$
as the average current (independent of coordinate $\xi$
) of the constant part of the profile. As figure 2(c,g) and corresponding insets show, the length and thus the charge of the injected beam get slightly smaller with the increase in the transition length $L$
(in particular, in figure 2(c), the charges of the injected bunches are 209, 175 and 148 pC, whereas in figure 2(g) the charges are 39, 28 and 21 pC). However, it is still mostly determined by the ratio between the initial and the final plasma density (fixed at 2 in our case), as this ratio determines the increase in the length of the bubble during the transition. The value of the injected beam current $J_w$
is even less dependent on the transition length (except the green line in figure 2(g) where the downramp length is already close to the injection threshold and injection does not result in a flat current profile).
Figure 2. Current profiles of injected beam $J_w$
corresponding to different lengths of plasma density downramp $L$
and different wakefield nonlinearity $n_b/n_0$
. The insets are wakefield structures for the plotted simulations. Plasma density before the downramp $n_0 = 5\times 10^{18}\,{\rm cm}^{-3}$
; $n_0/2$
after the downramp. The driver has the Gaussian shape with $\sigma _z = 2.5\,\mathrm {\mu } {\rm m}$
and $\sigma _r = 0.5\,\mathrm {\mu } {\rm m}$
.
As shown by figure 2(a,b), the higher value of ratio $n_b/n_0$
(and, thus, the wakefield nonlinearity) leads to the higher current $J_w$
. From the comparison of figures 2(a) and 2(b), one can see that shorter and steeper downramp enables injection at lower driver charge and lower nonlinearity. No injection at all is observed for the case where $n_b/n_0=10$
shown in (b), and the higher the wakefield nonlinearity $n_b/n_0$
is, the stronger the injected witness current becomes.
3.2. Injected beam current scaling
As shown in the previous section, as long as the density downramp provides stable injection, its properties do not significantly affect the injected current $J_w$
(defined now as the average value of the constant part of the injected bunch current profile), so it should mostly be determined by the parameters of the driver and their relation to the plasma density. It is suggested by figure 2 that the higher values of $n_b/n_0$
which should correspond to stronger nonlinearity lead to a stronger injected beam current. As we want to investigate the dependence of $J_w$
on the nonlinearity of the bubble, we need to introduce a measure of it first. In the previous section, we used the ratio $n_b/n_0$
as this measure, which is only suitable for a fixed shape and size of the driver. A more general measure requires introducing the energy properties of the bubble.
Following Lotov (Reference Lotov2004) and Golovanov et al. (Reference Golovanov, Kostyukov, Reichwein, Thomas and Pukhov2021, Reference Golovanov, Kostyukov, Pukhov and Malka2023), we introduce the quantity $\varPsi (\xi ) = \int (c W - S_z)\,{\rm d}^2r_\perp$
which depends on the comoving coordinate $\xi = z - ct$
and contains the integral over the transverse plane of the energy density $W$
and the longitudinal energy flux $S_z$
of both the electromagnetic field and the plasma particles. The value of $\varPsi$
is also equal to the total energy flux in the comoving window (Lotov Reference Lotov2004). As it has the dimension of power, we will refer to $\varPsi$
as ‘the power of the bubble’. In the absence of energy exchange with bunches, $\varPsi$
is a conserved property in any wake; drivers lead to the increase of $\varPsi$
, whereas accelerated witnesses decrease it. As a property describing the energetic properties of the bubble, it is equal to the total power of deceleration felt by the driving bunch and is also equal to the maximum achievable power of acceleration for the witness bunch. For the blowout regime of plasma wakefield, $\varPsi$
is fully determined by the size of the bubble $R_{\textrm {bub}}$
and grows with it; in the limit of a large bubble size ($k_p R_{\textrm {bub}} \gg 1$
), $\varPsi \propto (k_p R_{\textrm {bub}})^4$
(see Tzoufras et al. Reference Tzoufras, Lu, Tsung, Huang, Mori, Katsouleas, Vieira, Fonseca and Silva2008; Golovanov et al. Reference Golovanov, Kostyukov, Reichwein, Thomas and Pukhov2021). The power of the bubble $\varPsi$
serves as quantification of nonlinearity of the wake and it is the most important property of the bubble. Regardless of the shape the driver which creates a bubble with a certain value of $\varPsi$
, the effect of the bubbles with the same power on acceleration of particles and on downramp injection will be mostly the same. Therefore, we can expect that injection only depends on the value of $\varPsi$
.
According to Golovanov et al. (Reference Golovanov, Kostyukov, Reichwein, Thomas and Pukhov2021), for a blowout wakefield excited by a sufficiently tightly focused ($k_p r_b \ll 1$
) driver with the charge $Q$
and a flattop current profile of length $\xi _b$
, the power of the bubble in the large-bubble limit ($k_p R_{\textrm {bub}} \gg 1$
) is calculated using the following formula:
where $J_A = 4 {\rm \pi}\varepsilon _0\;m c^3 / e \approx 17\,{\rm kA}$
is the Alfvén current. It can be rewritten as
where we introduce the quantity $J_{\textrm {eff}}$
which we call ‘the effective current’ of the driver:
For high-current ($J_{\textrm {eff}} \sim J_A$
) sufficiently short electron bunches, the last factor in (3.2) is close to $1$
, and $\varPsi \propto J_{\textrm {eff}}^{3/2}$
. In general, this factor describes the weakening of the bubble when the bunch becomes long enough compared with the plasma wavelength. In the limiting case when the back of the bunch is already located in the accelerating phase, the bubble effectively transfers the energy from one part of the driver to another, limiting the possible efficiency of accelerating the witness, which corresponds to this factor tending to 0. When it becomes negative, this solution for $\varPsi$
is formally incorrect, which corresponds to the situation when the driver cannot fit inside the bubble.
Since $J_{\textrm {eff}}$
fully determines $\varPsi$
and, thus, the properties of the bubble for sufficiently short bunches, it should be the only combination of the driver and plasma parameters affecting the injection process. Although the model in Golovanov et al. (Reference Golovanov, Kostyukov, Reichwein, Thomas and Pukhov2021) is not necessarily strictly valid in the case of smaller-size bubbles, because this dependence holds in the important limit of large bubbles, we make a conjecture that $J_{\textrm {eff}}$
should be the most important combination of the driver's parameters.
To explore the influence of $J_{\textrm {eff}}$
on the injected witness current $J_w$
, we perform numerical PIC simulations with FBPIC (Lehe et al. Reference Lehe, Kirchen, Andriyash, Godfrey and Vay2016) for two types of drivers: 3-D Gaussian beams which are usually used to model the driver beam produced by a conventional accelerator (Joshi et al. Reference Joshi, Adli, An, Clayton, Corde, Gessner, Hogan, Litos, Lu and Marsh2018; D'Arcy et al. Reference D'Arcy, Aschikhin, Bohlen, Boyle, Brümmer, Chappell, Diederichs, Foster, Garland and Goldberg2019), and flattop longitudinal bunch which models the LWFA produced driver beam (Couperus Cabadağ et al. Reference Couperus Cabadağ, Pausch, Schöbel, Bussmann, Chang, Corde, Debus, Ding, Döpp and Foerster2021; Foerster et al. Reference Foerster, Döpp, Haberstroh, Grafenstein, Campbell, Chang, Corde, Couperus Cabadağ, Debus and Gilljohann2022). To prevent any effects due to the evolution of the driver beam, we artificially freeze it by setting its initial particle energy to 50 GeV. The ratio between the plasma density before and after the downramp is fixed to 2 (going from $n_0$
to $n_0/2$
) and the downramp length is between $20\,\mathrm {\mu }$
m and $160\,\mathrm {\mu }$
m. For Gaussian drivers with the density profile defined as $n_b \exp [-r^2 / 2 \sigma _r^2 - (\xi -\xi _0)^2 / 2 \sigma _z^2]$
, the beam transverse size $\sigma _r = 0.5\,\mathrm {\mu }$
m. Flattop drivers correspond to the cylinder of constant density $n_b$
, length $\xi _b$
and fixed radius $r_b = 0.6\,\mathrm {\mu }$
m.
As the effective current $J_{\textrm {eff}}$
defined by (3.3) depends on the current of the driver $J_b$
and its length $\xi _b$
, we study the relationship $J_w$
and $J_{\textrm {eff}}$
by varying these two parameters. Of course, the value of $J_{\textrm {eff}}$
also depends on the plasma density, so we use the plasma density $n_0$
before the downramp to calculate the plasma wavenumber $k_p$
used in (3.3). Because definition (3.3) is written for flattop profiles, for Gaussian beams we use $\xi _b = \sqrt {2} \sigma _z$
and the peak value of $J_b$
. This combination was shown to provide a good approximation for $\varPsi$
in simulations (not presented in the paper). The results for varying the driver beam current $J_b$
at fixed length are shown in figure 3, and the results for varying the driver bunch length while the current is fixed are presented in figure 4 for different plasma densities.
Figure 3. The relationship between $J_w$
and $J_{\textrm {eff}}$
for different types of beams. Each line corresponds to the same value of plasma density (in ${\rm cm}^{-3}$
in the legend), and the driver beam current $J_b$
is varied while the beam longitudinal size ($\sigma _z$
or $\xi _b$
) remains the same. The dashed black line corresponds to the linear dependence on $J_w$
.
Figure 4. The relationship between $J_w$
and $J_{\textrm {eff}}$
for different types of beams and different beam currents. Each dashed line corresponds to the same value of plasma density (in ${\rm cm}^{-3}$
in the legend), and the driver beam longitudinal size $\xi _b$
is varied while the driver current $J_b$
remains the same. The dashed black line corresponds to the linear dependence on $J_w$
. The dotted line shows the prediction according to (3.4) with a correction $c_2 = 0.8$
for the Gaussian shape and $0.4$
for the rectangular shape. The vertical lines correspond to the $J_{\textrm {eff}}$
value above which the length of the driver becomes comparable to the plasma wavelength ($\sigma _z = \lambda _p/2\sqrt {2}$
for the Gaussian shape and $\xi _b = \lambda _p$
for the rectangular shape).
For the fixed beam length (figure 3), we observe a linear dependence of the injected current on the effective current, which seems to indicate that $J_w \propto \varPsi ^{2/3} \propto J_{\textrm {eff}}$
. The slope of the linear dependence depends only on the shape of the driver and remains the same for different plasma densities and lengths of the driver. We observe that the witness current becomes slightly lower for denser plasmas (thus, larger values of $k_p$
), which corresponds to the lowering efficiency due to the driver length in plasma units, as predicted by the length factor in (3.2). At very low $J_{\textrm {eff}}$
, the injected current goes to 0, which corresponds to the injection threshold discussed in § 3.1.
To study the dependence of the witness current $J_w$
on the driver length $\xi _b$
in more details, we also plot the dependence of $J_w$
on the effective current $J_{\textrm {eff}}$
for a fixed current of the driver, which means that the increase in $J_{\textrm {eff}}$
corresponds to the increase in its length (figure 4). For smaller $J_{\textrm {eff}}$
and, thus, shorter driver lengths, the already found linear trend with the same slope is again observed for both types of drivers. But for higher $J_{\textrm {eff}}$
and longer driver length, we see the deviation of $J_w$
from the linear trend and the decline in it. This behaviour is qualitatively consistent with (3.2) which predicts that the elongation of the driver lowers the efficiency of exciting a bubble, so we can expect that the witness current $J_w \propto \varPsi ^{2/3}$
behaves as
By choosing the value of $c_2$
, the qualitative behaviour of $J_{\textrm {eff}}$
can be recreated using this formula (see the dotted lines in figure 4).
However, this formula cannot show the full picture, because downramp injection for longer drivers is more complex than for short drivers. According to (3.2), the power of the bubble changes as the driver propagates through the density downramp which changes the plasma wavenumber $k_p$
. For short drivers, $\varPsi \propto J_{\textrm {eff}}^{3/2} \propto n_p^{1/2}$
scales exactly the same for all drivers with the same $J_{\textrm {eff}}$
and, thus, the injected current is still fully determined by $J_{\textrm {eff}}$
. However, longer drivers which have a lowered efficiency in dense plasma begin exciting the bubble more effectively when transitioning into lower-density plasma at the downramp, as their length in plasma units becomes small. An example of such a driver is shown in figure 5: in the initial higher density it is too long to fit inside the bubble, so during the density downramp it passes through the point of having almost zero efficiency of exciting a bubble ($\varPsi \approx 0$
), and then $\varPsi$
starts growing again as the driver becomes shorter than the excited bubble. In addition to the complex behaviour of the nonlinearity of the bubble, the corresponding velocity of the back of the bubble which determines the injection threshold will also significantly depend on the length of the driver. The dynamics of downramp injection for such cases cannot be reduced to a simple formula such as linear dependence on $J_{\textrm {eff}}$
or even corrected (3.4) and requires the full description of the evolution of the bubble in the downramp. However, as explained in § 2, these cases are suboptimal for efficient utilisation of the driver's energy and should be avoided. In more optimal cases, when the driver is comparatively short, the linear dependence of $J_w$
and $J_{\textrm {eff}}$
holds.
Figure 5. Electron density distribution in a wake excited by a long rectangular-shape driver before and after the downramp. The downramp provides transition from $n_0$
to $n_0/2$
where $n_0 = 2.5\times 10^{18}\,{\rm cm}^{-3}$
. The parameters of the driver: $n_b = 40 n_0$
, $\xi _b = 42.5\,\mathrm {\mu }$
m and $r_b = 0.6\,\mathrm {\mu }$
m.
4. Conclusion
We have studied a PWFA driven by a short high-current electron bunch. The accelerating structure in this case remains stable until the moment when some of the electrons of the driver lose all their energy, which is the condition which determines the acceleration length optimising the driver-to-plasma energy transfer efficiency. The dependence of the witness bunch generated by density-downramp injection on the parameters of the driver has been studied. We have shown that a steeper downramp enables injection for a weaker driver, but the current of the injected bunch does not significantly depend on the downramp length as long as the criterion for the injection is met. The witness current of the injected bunch mostly depends on the effective current of the driver $J_{\textrm {eff}}$
which combines both the driver's current and its length in plasma units. The dependence of the injected witness current on the effective current is linear for reasonably short drivers.
Acknowledgements
Editor Luís O. Silva thanks the referees for their advice in evaluating this article.
Funding
This work was supported by the Fondation Jacques Toledano, the Schwartz Reisman Center for Intense Laser Physics, and by EIC ebeam4therapy grant.
Declaration of interests
The authors report no conflict of interest.
Data availability statement
The simulation data are available on request from the corresponding author, A.G.
Open access
