Correlated Energy Spread Compensation in Multi-Stage Plasma-Based Accelerators

The extreme electromagnetic fields sustained by plasma-based accelerators allow for energy gain rates above 100 GeV/m but are also an inherent source of correlated energy spread. This severely limits the usability of these devices. Here we propose a novel compact concept which compensates the induced energy correlation by combining plasma accelerating stages with a magnetic chicane. Particle-in-cell and tracking simulations of a particular 1.5 m-long setup with two plasma stages show that 5.5 GeV bunches with a final relative energy spread of $1.2\times10^{-3}$ (total) and $5.5\times10^{-4}$ (slice) could be achieved while preserving sub-micron emittance. This at least one order of magnitude below current state-of-the-art and paves the way towards applications such as Free-Electron Lasers.

The extreme electromagnetic fields sustained by plasma-based accelerators allow for energy gain rates above 100 GeV/m but are also an inherent source of correlated energy spread.This severely limits the usability of these devices.Here we propose a novel compact concept which compensates the induced energy correlation by combining plasma accelerating stages with a magnetic chicane.Particle-in-cell and tracking simulations of a particular 1.5 m-long setup with two plasma stages show that 5.5 GeV bunches with a final relative energy spread of 1.2 × 10 −3 (total) and 5.5 × 10 −4 (slice) could be achieved while preserving sub-micron emittance.This at least one order of magnitude below current state-of-the-art and paves the way towards applications such as Free-Electron Lasers.
Plasma-based accelerators (PBAs), driven either by charged particle beams (plasma wakefield accelerator, PWFA [1]) or intense laser pulses (laser wakefield accelerator, LWFA [2]), are able to sustain accelerating gradients in excess of 100 GeV/m [3].These extreme gradients are orders of magnitude higher than those achievable with radiofrequency technology and offer a path towards miniaturized particle accelerators with ground-breaking applications in science, industry and medicine [4].
Steady progress over the past decades has led to the successful demonstration of electron bunches with multi-GeV energy [5][6][7][8], micron-level emittance [9,10] and kiloampere current [11,12].However, the high amplitude and short wavelength (∼ 100 µm) of the wakefields naturally imprint a longitudinal energy correlation (or chirp) along the accelerated (witness) bunch, leading to a large relative energy spread typically on the 1-10 % range [13].This is a long-standing issue for PBAs which critically impacts the beam quality [14], particularly for applications such as Free-Electron Lasers (FELs) [15] where a relative energy spread 0.1% is required [16].
Solving this issue is therefore key for demonstrating the usability of these devices.A well known concept for mitigating the correlated energy spread is that of beam loading [17][18][19], in which the witness bunch itself is used to flatten the slope of the accelerating fields.This, however, relies on a very precise shaping of the current profile and has yet to be demonstrated with the desired performance.Furthermore, since the optimal profile depends on the wakefield structure, a certain energy spread will always develop in LWFAs, where the wakefield experienced by the bunch will change due to the laser evolution [20] as well as dephasing [21].Alternative ideas have also been proposed in order to achieve, in average, a flat accelerating gradient.These include modulating [22] or tailoring [23] the plasma density profile as well as injecting a secondary bunch [24], but they show limited success or remain to be experimentally realized.A different approach contemplates stretching the bunch in order to minimize the slice energy spread [25].
In this Letter we propose a novel concept for compen-sating the correlated energy spread by taking advantage of the naturally occurring energy chirp.The scheme, illustrated in Fig. 1, consists mainly on two identical plasma accelerating stages joined by an intermediate magnetic chicane in which the longitudinal energy correlation of the bunch is inverted.Thus, the energy chirp generated in the first stage is compensated in the second.Numerical simulations with the Particle-in-Cell (PIC) code FBPIC [26] as well as the tracking codes AS-TRA [27] and CSRtrack [28] show that multi-GeV beams with unprecedented energy spread could be obtained with this method.Although LWFA stages are used here, the core idea behind the scheme would be equally valid in the case of a particle driver.
In order to introduce this concept, the blowout regime [29,30] of plasma acceleration will be considered.In this case the laser or beam driver is able to expel all background plasma electrons, leaving behind an ion cavity with uniform focusing gradient, K = (m/2ec)ω 2 p , and an approximately constant longitudinal electric field slope, p , along most of the accelerating phase.Here ω p = n p e 2 /m 0 is the plasma frequency, e and m the electron charge and mass, 0 the vacuum permittivity and n p the unperturbed plasma density.In order to describe the position and energy of the particles along the accelerator it is also useful to introduce the speed-of-light coordinate, ξ = z − ct, as well as the relativistic Lorentz factor, γ = 1/ 1 − (v/c) 2 , where t is the time and v and z are, respectively, the particle velocity and longitudinal position in the laboratory frame.Within the generated cavity, electrons perform transverse oscillations (known as betatron motion) with a frequency ω β (t) = eK/mγ(t), while their energy evolves as γ(t) = γ 0 − (e/mc)E z t.
For a particle bunch with average energy γ(t) = γ(t) centered at ξ, the longitudinal chirp can be expressed as χ(t) = ∆ξ∆γ(t) / ∆ξ 2 γ(t), where ∆γ(t) = γ(t)−γ(t) and ∆ξ = ξ − ξ.A simple expression for the chirp evolution within a plasma stage can be obtained if a constant E z is assumed.This yields which tends asymptomatically to χ = E z /E z and where χ 0 and γ0 are the initial bunch chirp and energy.If the bunch length is σ z = ∆ξ 2 , this induces a correlated energy spread σ corr γ (t)/γ(t) = χ(t)σ z .In a 2-stage accelerator as in Fig. 1, the accumulated chirp after a first stage of length L p,1 for an initially unchirped bunch will be χ 1 = −(e/mc 2 )E z,1 L p,1 /γ 1 .Therefore, if the longitudinal phase space of the bunch is inverted at this point such that χ1 = −(σ z,1 /σ z,2 )χ 1 is obtained, the correlated energy spread could be compensated in a following stage fulfilling E z,2 L p,2 = (mc 2 /e) χ1 γ1 .For a symmetric inversion (σ z,1 = σ z,2 ), using two identical plasma stages (same E z and L p ) would be the simplest setup.
This longitudinal phase space inversion can be performed with a conventional chicane.As illustrated in Fig. 2, this device is composed by 4 dipole magnets in which particles undergo an energy-dependent trajectory bend.With respect to a hypothetical reference particle with γ = γ ref , those with γ > γ ref experience less bending and therefore a shorter path length, while the opposite occurs for those with γ < γ ref .Defining δ = (γ − γ ref )/γ ref , the path length differences after the chicane, ∆ξ ch , can be expressed with respect to the reference particle as where T 566 −3/2R 56 [31].To first order, the R 56 coefficient can be simply determined as R 56 = ∆ξ ch /δ = ∆ξ ch /χ∆ξ (assuming γ ref = γ).These path length differences allow for a certain control of the longitudinal phase space.For example, full bunch compression can be achieved if ∆ξ ch exactly compensates the initial offsets with respect to the bunch center, ∆ξ.Similarly, inverting the chirp, i.e inverting the bunch along ξ, can be achieved if ∆ξ ch = 2∆ξ.This implies that a chicane with R 56 = 2/χ is required.
However, successfully performing this process requires χ, and therefore E z , to be uniform along the bunch.Thus, this scheme is ideally suited for weakly beamloaded wakefields, where E z is not perturbed by the bunch, or where the beam-loading effect linearly modifies the accelerating fields [19].Additionally, in order to maintain a purely linear chirp after the chicane, the higher order terms in Eq. ( 2) need to be minimized.This is required because the non-linear contributions cannot be compensated in the second stage and therefore will lead to increased energy spread.In particular, keeping the higher order effects below 1%, i.e |T 566 δ/R 56 | 10 −2 , requires σ z 0.015c/ω p if typical blowout fields are assumed (E z = (m/2e) ω 2 p and E z = (mc/e)ω p ).These contributions can also be mitigated by sextupole magnets [32] or by imprinting non-linearities in the bunch to compensate those from the chicane, e.g., by optimizing a non-uniform E z along the bunch.Still, for a non-uniform χ, the scheme could be designed to mitigate the chirp at the bunch core in order to minimize the slice energy spread for applications such as FELs.
Once the chicane R 56 is determined, the magnet length L m and bending angle θ experienced by the reference particle can be directly determined from R 56 = −2θ 2 (L d + 2L m /3) [31], where L d , as defined in Fig. 2, is the length of the drift space between the first and second as well as the third and fourth dipoles.The magnetic field strength can then be obtained as B = (mc/e)θγ ref /L m .Assuming a bunch with χ = E z /E z and typical blowout fields, it can be obtained that R 56 4c/ω p 1 mm (for n p 10 16 cm −3 ) which, considering L m ∼ L d ∼ 0.1 m leads to θ 0.1 rad.The high χ characteristic of PBAs therefore allows for a very compact chicane design (∼ 1 m) while requiring a very small bending angle.This greatly minimizes the impact of Coherent Synchrotron Radiation (CSR) [33] on the beam parameters.
A possible implementation of this scheme is shown in Fig. 1.Two LWFA stages accelerating an externally injected bunch are joined by a magnetic chicane, including active plasma lenses (APLs) [34] to transport the beam and plasma mirrors [35] to couple the laser pulses in and out.The two required laser drivers could originate from splitting a single original pulse and therefore be intrinsically synchronized.APLs consist on circular, gas-filled capillaries with typically sub-mm radius, R c , on which a multi-kV discharge is applied by electrodes at both ends, causing a breakdown of the gas.A current is then driven though the ionized plasma, generating radially symmetric focusing fields with kT/m gradients which are linear up to a radius r R c /2 [34].The strong focusing fields make these devices ideal for transporting the strongly divergent and high energy spread bunches coming out of the first stage, as they can be placed very close to it, thus mitigating emittance growth in the drift [14], and they can focus on both planes, which significantly reduces their chromaticity with respect to other focusing systems [34].
The plasma mirrors are necessary in order to maintain the compact footprint the design, as they can be placed close to the laser focal spot.In comparison, conventional mirrors would have to be placed several meters away in order to avoid damage from the multi-TW lasers required by LWFAs.As an example, tape-based plasma mirrors have been successfully used to remove and couple in laser pulses [36].However, due to their thickness, they could also perturb the witness bunch due to scattering as well as the plasma created by the incident laser.A promising alternative are liquid crystal films [37], which can have a thickness down to the few-nm scale.A specific implementation is not considered here.
This scheme, apart from offering an energy spread compensation, could also reduce the sensitivity to critical issues such as the timing jitter between laser driver and witness bunch [38] or the hosing instability [39].The timing jitter is one of the main challenges of external injection, as it translates into a large energy jitter at the LWFA exit due to the large E z .However, thanks to the chicane in this scheme, an injection offset with respect to the ideal phase in the first LWFA would translate, to first order in Eq. ( 2), into the opposite offset at the second LWFA, thus providing a stable average accelerating field and energy output.Furthermore, since the bunches are accelerated with a large energy chirp, the hosing instability is also mitigated [40,41].
In order to test the performance of this scheme, startto-end simulations for a particular set of parameters have been performed.The LWFA stages and APLs have been simulated using the spectral, quasi-3D PIC code FBPIC, while the tracking code ASTRA has been used for the remaining beamline elements taking into account 3D spacecharge effects.Additionally, CSRtrack has also been used to account for CSR effects in the chicane.
Motivated by the parameters from the EuPRAXIA de-sign study [42], the simulated setup aims at providing 5 GeV electron beams suitable for FEL applications, i.e., a peak current in the kA range, sub-micron emittance and, specially, an energy spread 0.1%.For this purpose, we consider an externally injected Gaussian electron bunch with an initial energy of 250 MeV, a 0.5% energy spread with no chirp, a normalized transverse emittance n,x = 0.5 µm rad, 10 pC of charge, a FWHM duration τ = 5 fs and a peak current I peak 2 kA.The bunch transverse size, σ x , is matched [43,44] to the plasma focusing fields in order to prevent emittance growth.This requires the beam beta function, β x = γσ 2 x / n,x , to satisfy β x = c/ω β at the LWFA entrance.The normalized emittance is defined as n,x = ( x 2 p 2 x − xp x 2 ) 1/2 /mc, where p x is the transverse particle momentum.Electron bunches within this range of parameters can be produced with conventional accelerators [45].The two identical LWFAs have a length L p = 8 cm and a parabolic transverse density profile n p = n p,0 + r 2 /πr e w 4 0 , where n p,0 = 10 17 cm −3 is the on-axis plasma density, r the radial coordinate, r e the classical electron radius and w 0 the spot size of the laser driver.Plasma cells in this range of parameters have been recently demonstrated [46] and provide ideal guiding properties for the driving laser pulse.For simplicity, a longitudinal flat-top plasma density profile has been considered, although the presence of smooth plasma-to-vacuum transitions would be beneficial for electron beam matching [47] and emittance growth minimization [48].This choice of electron beam and plasma parameters also helps in reducing below the 10 −4 level the energy spread generated in the LWFAs due to slice mixing from betatron motion [49].Each LWFA is driven by a 40 J, 0.75 PW laser pulse with a peak normalized vector potential a 0 = 3, a spot size w 0 = 50 µm and a FWHM duration τ 0 = 50 fs.This laser can be successfully guided throughout the LWFAs, which provide an energy gain of ∼ 2.6 GeV each.The APLs are placed 3 cm away from the LWFAs, leaving enough space for the plasma mirrors.They provide a focusing gradient of 3 kT/m, as measured experimentally [34], and have a length of 6.6 cm, optimized to achieve a beam waist at the chicane center.A plasma density n APL p = 10 15 cm −3 has been considered in order to minimize the impact of beamdriven wakefields [50].The chicane has a total length of 1.2 m, with L d = 12.5 cm and dipoles with L m = 20 cm and B = 0.54 T for a bending angle of θ = 0.011 rad.
An overview of the simulation results can be seen in Fig. 3.The electron bunch leaves the first LWFA with preserved emittance, an energy of ∼ 2.9 GeV and a chirp χ −0.031 µm −1 , which induces a total relative energy spread ∼ 2%.As a consequence, the projected emittance grows after the accelerating stage until the beam divergence is controlled by the APL.The maximum beam size within the lens is σ r 15 µm R c /2 and therefore experiences linear focusing.ASTRA and CSRtrack simulations show that the influence of space-charge and CSR on the beam parameters is negligible thanks to its GeV energy and the small bending angle.The beam is then focused in the following APL and injected into the second LWFA, where it gains an additional ∼ 2.6 GeV for a final energy of ∼ 5.5 GeV while reducing its relative energy spread down to 0.12%.As seen in Fig. 4, the relative energy spread is reduced by a factor ∼ 20 with respect to a case with no chirp compensation.However, it is only a factor ∼ 5 smaller than the value at injection, while the energy has increased more than 20 times.This is due the non-linear contributions in the chicane, as seen in Fig. 4b, which increase the final projected energy spread.
The emittance evolution in the drifts and APLs is well controlled, achieving a final value of 0.69 µm rad (∼ 37% increase).This growth arises from the large energy chirp along the bunch, which causes individual slices to diverge (or converge) at different rates in the drifts and to have a different betatron frequency in the APLs, leading to increased projected emittance.However, it should be noted that this growth is not the same in both APLs, but it is more moderate in the first one.This is due to the beam-induced focusing wakefields which, for a short bunch, grow linearly towards the back [50] and, for a case with χ < 0, can mitigate the projected emittance growth by equalizing ω β along the bunch.This suggests that wakefields in APLs, which are typically regarded as a key limitation of these devices [50], can also be useful and could be optimized for emittance preservation in bunches with a large negative energy chirp.Another consequence of this slice decoherence is that β x will evolve differently along the bunch.Therefore, not all slices will be matched to the focusing fields in the second LWFA, causing the oscillations in β x seen in Fig. 3a.
Regarding the slice parameters, which are key for FEL applications, this scheme achieves a final average relative energy spread and emittance of 5.5 × 10 −4 and 0.58 µm rad, respectively, assuming 1 µm slices.This energy spread is at least one order of magnitude lower than in state-of-art LWFAs and would satisfy the requirements for an X-ray FEL.The growth of the slice parameters in the chicane arises from the reduction in bunch length, which reaches its minimum at the chicane center and becomes shorter than a single slice.Thus, the slice parameters approach those of the whole bunch.The presented scheme therefore offers a path towards ultra-low energy spread (∼ 10 −4 ) beams while reducing the sensitivity to issues such as the timing jitter for external injection and the hosing instability.The concept is ideally suited for PBAs with weakly or linearly beamloaded wakefields.Electron beams obtained with this method could enable ground breaking applications, such as compact FELs, and could be accelerated to higher energies thanks to its modular design, allowing other applications such as plasma-based colliders.This could be achieved either by introducing additional plasma stages while keeping a single central chicane, or by repeating multiple sections of one chicane every two plasma accelerating modules.Furthermore, this scheme could also be modified as an extreme bunch compressor by halving the chicane R 56 .CSRtrack simulations with the same bunch produced in the first LWFA show that ultra-short bunches with a duration down to ∼ 70 attoseconds and a ∼ 50 kA peak current could be obtained, with possible applications in ultrafast science [51].This project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No. 653782.We acknowledge the use of the High-Performance Cluster (Maxwell) at DESY for the numerical simulations.We would like to thank T. Mehrling, S.K. Barber, S. Steinke and M. Dohlus for useful discussions.

FIG. 1 .
FIG. 1. Overview of the proposed accelerator concept (not to scale).

FIG. 2 .
FIG.2.Working principle of a magnetic chicane with R56 = 2/χ.The bunch longitudinal phase space is shown at the chicane entrance (a), middle (b) and exit (c).Darker color implies higher energy.

FIG. 3 .
FIG. 3. Evolution of the beta function (a), normalized emittance (b), energy spread (c) and energy (d) along the beamline.The emittance and beta function in the chicane are corrected for dispersion up to second order.

FIG. 4 .
FIG.4.Comparison of the final bunch properties with respect to a case with acceleration in a single LWFA (same driver and plasma profile but ∼ 20 cm long).The bunch current (a), longitudinal phase space (b) and energy profile (c) normalized to the peak number of counts, N , are shown.The total energy spread evolution during acceleration can be seen in (d).