Next Article in Journal
A Configurable 64-Channel ASIC for Cherenkov Radiation Detection from Space
Next Article in Special Issue
Asymmetric Dual-Grating Dielectric Laser Accelerator Optimization
Previous Article in Journal / Special Issue
Opportunities for Bright Beam Generation at the Argonne Wakefield Accelerator (AWA)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Scaling Laws for Radiofrequency-Based Pulse Compression in Ultrafast Electron Diffraction Beamlines

Department of Physics and Astronomy, University of California, Los Angeles, CA 90095, USA
*
Author to whom correspondence should be addressed.
Instruments 2023, 7(4), 49; https://doi.org/10.3390/instruments7040049
Submission received: 10 October 2023 / Revised: 19 November 2023 / Accepted: 20 November 2023 / Published: 29 November 2023

Abstract

:
We present an envelope equation-based approach to obtain analytical scaling laws for the shortest pulse length achievable using radiofrequency (RF)-based bunch compression. The derived formulas elucidate the dependencies on the electron beam energy and beam charge and reveal how relativistic energies are strongly desirable to obtain bunches containing 1 million electrons with single-digit femtosecond pulse lengths. However, the non-linearities associated with the RF curvature and the beam propagation in drift spaces significantly limit the attainability of extreme compression ratios. Therefore, an additional higher frequency RF cavity is implemented, which linearizes the bunch compression, enabling the generation of ultrashort beams in the sub-femtosecond regime.

1. Introduction

Ultrafast electron scattering requires the generation of very short electron bunches to capture the fastest physical processes [1,2]. Due to the repulsive effect of space-charge forces, one critical challenge in this field is related to packing as many electrons as possible in a short bunch [3]. In ultrafast electron diffraction (UED), pushing the electron energy to relativistic levels has helped in minimizing the space-charge effects, concurrently bringing other advantages such as longer penetration depths, reduced group velocity mismatch, and suppressed inelastic scattering background [4,5,6,7]. Over recent years, UED beamlines have seen continuous improvement in the achievable temporal resolution due to the introduction of techniques borrowed from accelerator physics based on the use of a time-dependent radiofrequency (RF) electric field to compress the electron bunch during its propagation in the beamline [8]. RF compression using 3 GHz resonant cavities has been applied to both non-relativistic and relativistic electron beamlines for UED [9,10,11], yielding bunch lengths down to the single-digit fs in the latter case.
While the discussion in this paper focuses on the electron bunch length, it is important to recognize that there are many additional factors, other than the temporal duration of the probe pulse, that contribute to the actual temporal resolution limit in a specific UED setup such as temporal jitter, group velocity mismatch, and laser pulse length. For example, to counteract the additional temporal jitter introduced by RF-based compression, naturally synchronized laser-generated higher-frequency waves have been used to impart an energy chirp on the beam in more complex coupling structures and drive compression dynamics [12,13].
In any case, though, to push the boundary of the UED technique, it is critical to understand the limits in beam compression and how the various beamline parameters such as charge, energy, cavity voltage, and frequency affect the shortest bunch duration achievable. The minimum bunch length at the sample results from a complex interplay between the details of the bunching dynamics and the longitudinal space-charge forces in the beam so that typically UED practitioners resolve to particle tracking simulation codes to design the beamline and predict the beam dynamics. Generally, there is agreement with experimental results [14]. Still, particle simulations only deal with specific beamline setups, typically lack generality, and might not offer an immediate answer to how to improve compression in a given configuration.
It would be beneficial to have a unified formalism describing beam dynamics in RF-compression UED beamlines, covering both relativistic and non-relativistic cases while including the space-charge effects. To this end, we employ the longitudinal envelope equation formalism to highlight the interplay between longitudinal emittance and space-charge forces on pulse evolution. The single-particle dynamics presentation builds on previous works of Floettman and Zeitler [15,16] that pointed out the role of the non-linearities in the beam compression process. The collective effects are then considered in the approximation that the beam aspect ratio remains constant along the beamline, thus decoupling the longitudinal dynamics from the transverse beam size evolution. While this is a somewhat restrictive assumption, it is experimentally relevant (the beams in UED are usually focused transversely and longitudinally at the sample) because this approach yields an upper bound estimate for the minimum bunch length. In this case, space-charge forces are over-estimated for the situation in which the transverse spot size is kept large during the compression. Using the constant aspect ratio approximation, we can extend on the previous work and obtain analytical formulas for the minimum bunch length at the longitudinal waist that are valid in the presence of space charge. The expressions presented can then be used to guide the system optimization, compare parameter choices at different facilities, and evaluate mechanisms for further improving the bunch length.
The simple cartoon in Figure 1 illustrates the dynamics under study. Essentially, a finite electron beam propagates through an RF buncher cavity where electromagnetic fields oscillate with angular frequency ω = k c at the zero-crossing phase, where k is the cavity wave number and c is the speed of light. Ideally, the input bunch length satisfies k σ z < < 1 , where σ z is the electron beam rms bunch length, so only a small phase window of the wave is sampled by the beam and the chirp imparted on the beam is predominantly linear. However, in our discussion, we keep the higher-order terms in the energy modulation expansion to elucidate their role in the final bunch length. In the propagation region after the buncher, due to the strong energy chirp, the tail of the beam begins to catch up, while the head of the beam slows down. Finally, at some location downstream of the buncher, ideally arranged to be the sample plane or the interaction point of the UED experiment, the minimum bunch length occurs when the phase space distribution is vertically aligned.
We strive to keep all the formulas in the paper as general as possible. For example, by not assuming the ratio of beam velocity to the speed of light, β = 1 , so that formulas can be applied to different RF compression setups (non-relativistic, MeV UED beamlines, as well as higher-frequency compression schemes) once the parameters are scaled accordingly. For this reason, we use two different example cases, loosely based on the UED beamlines at the UCLA Pegasus laboratory [10], to benchmark the agreement between the analytical framework and particle tracking simulations. The reference parameters used for this study are reported in Table 1.
The particle tracking software utilized to benchmark the analytical framework presented in this paper is General Particle Tracer (GPT) [14]. GPT is a full 3D particle tracking simulation software that enables a more complete study of 3D and non-linear effects of charged particle beam dynamics in electromagnetic fields and hosts a variety of visualization tools for windows-based computers. GPT also offers a selection of default fieldmaps of common beamline elements such as RF cavities and charged particle optics, which can be easily implemented to provide preliminary beamline models. Typically, as more realistic fieldmaps are generated, they can be implemented as custom elements in the simulation.
The paper is structured as follows. We first introduce the envelope equation formalism to describe the evolution of the bunch duration through the system [17]. We approximate the buncher as a thin lens, compute the main contribution to emittance growth and analyze the ballistic dynamics in the drift while neglecting space charge, thus essentially re-obtaining the results first presented in Floettmann [15], then in Zeitler et al. [16]. We then include the space-charge repulsion term in the envelope equation for a Gaussian beam using the approximation of a constant aspect ratio and arrive at an expression for the minimum achievable bunch length. We build on the formalism using it to describe different current profiles. Finally, in the light of the findings in this paper, we review the use of an additional higher-frequency RF cavity for compensation of non-linearities to reach sub-fs bunch lengths as originally proposed by Floettmann [15].

2. Longitudinal Envelope Equation

2.1. Envelope Equation

Among many possible choices for defining bunch duration (e.g., full width half maximum FWHM, or full width containing 50% of the charge FW50), in our discussion, we select the second-order moment of the longitudinal profile, or root mean square (RMS) bunch length, as the primary quantity to follow the evolution of in the RF compression beamline. The RMS is defined by
σ z = z 2 ,
where z is the longitudinal particle coordinate relative to the center of the bunch. represents an expectation value over the beam distribution function. All the other definitions (i.e., FWHM and FW50) are simply proportional to the RMS bunch length using a proper order-of-unity pre-factor which depends on the beam distribution shape. The RMS beam size evolves along the beamline axial coordinate s according to equation [18]
σ z = z z σ z + ϵ z , z 2 σ z 3 ,
where quantity ϵ z , z 2 is the square of the RMS longitudinal trace space emittance, explicitly written in terms of the distribution moments as
ϵ z , z 2 = z 2 z 2 z z 2 ,
and z is the relative velocity deviation from the average beam velocity. If the particle dynamics are linear in ( z , z ) coordinates, then the trace space emittance is a conserved quantity [19].
Ideally, the longitudinal equation of motion does not depend on transverse coordinates. On the other hand, the force acting on the particles might have a transverse dependence (i.e., z = F ( z , r , s ) ) and in general the first term on the right side of Equation (2) becomes proportional to z F ( z , r , s ) . Heuristically, we separate the contributions to the longitudinal force into a term associated with the external RF fields and a term proportional to the beam current, which encodes the effect of the space-charge force, i.e., F = F R F + F s c . After averaging over the transverse beam distribution, the latter depends on the horizontal and vertical RMS transverse sizes σ x = σ y = σ , effectively coupling the longitudinal and transverse envelope equations.
In most practical cases, we can model the RF buncher as a thin lens, and F R F ( z ) provides an impulse force proportional to δ ( s ) . In this case, the only effect of the RF buncher is to modify the initial conditions for the envelope equation through the application of a negative energy chirp to the phase space, which causes the beam to begin the bunching process (i.e., σ z 0 < 0 ). If the applied energy chirp is non-linear (which is the most common case), then the change in longitudinal emittance and the initial conditions must be evaluated at the exit of the buncher.
Meanwhile, the space-charge term z F s c ( s , σ z , σ ) acts over the whole time of flight in the drift following the buncher. This term is related to the first-order Taylor expansion of the space-charge-induced longitudinal electric field near the center of the bunch. We discuss this in detail later in the paper, but we can anticipate some of the results derived below to orient this initial discussion. Typically, this term can be evaluated analytically for simple longitudinal distributions and written as
z F s c ( s , σ z , σ ) σ z = K L σ z 2 ,
where K L g N r c / β 2 γ 5 is the longitudinal perveance proportional to the classical electron radius, r c , the number of electrons in the bunch N, and g is a geometry factor which depends on the beam aspect ratio.
In what follows, we assume a constant aspect ratio for the beam, a transverse RMS spot size that decreases along the beamline proportionally to the longitudinal RMS bunch length. This assumption simplifies the s and σ dependencies in K l and solves the longitudinal envelope equation independently from the transverse dynamics. The approximation works well even if the longitudinal and transverse spot sizes are not exactly proportional to each other along the beamline, but the aspect ratio remains within a factor of two of the initial value. Of course, this type of transverse focusing may not occur in all setups. In those cases, the transverse beam size remains large. Thus, the aspect ratio can increase by orders of magnitude during compression, and the approximation fails, but it does yield an overestimate of the waist size in the presence of space-charge effects. In addition, it is essential to note that for general beam distributions, the non-linearities in the bunch self-fields cause emittance growth during the propagation, violating the premises of our approach, which uses the envelope equation to treat the problem and assumes a constant longitudinal trace space emittance. Instead, one must self-consistently track all particles in the field generated by the charge distribution. Nevertheless, if one desires ultrashort bunch lengths, this situation should be avoided because space-charge effects should be relatively small.

2.2. Solution in a Drift

Let us start the discussion from the most straightforward and physically relevant case in which the propagation occurs in a drift and space-charge forces can be neglected, i.e., F = 0 . The shortest bunch length along the line can then be found by recasting Equation (2) as
1 2 d d s σ z 2 = ϵ z , z 2 σ z 3 σ z ,
which can be integrated exactly
σ z f 2 σ z 0 2 = ϵ z , z 2 1 σ z 0 2 1 σ z f 2 .
The waist position is a local minimum for σ z ; thus, we can set σ z f = 0 . Right after the buncher, we can write σ z 0 = z 0 z 0 σ z 0 σ z 0 f where the first term accounts for any incoming correlations in the phase space and the second term accounts for the linear chirp imparted by the buncher cavity. After substituting into Equation (6), the waist size can be written as
σ z f = 1 1 ϵ z , z 2 z 0 z 0 σ z 0 σ z 0 f 2 + 1 σ z 0 2 f ϵ z , z σ z 0 1 f z 0 z 0 σ z 0 2 1 ,
where we neglect the term 1 σ z 0 2 which for large compression factors is always much smaller than 1 σ z 2 .
By inspecting the factor accounting for initial correlations, we can see that an initial negative chirp (i.e., z 0 z 0 < 0 ) effectively shortens the focal length and final bunch length in the system. On the other hand, in most cases, an RF buncher is added to the system to achieve strong compression and f z 0 z 0 σ z 0 2 < < 1 , i.e., incoming correlations are small compared to the linear correlation imparted by the RF fields. In this case, the final bunch length at the waist is given by f ϵ z , z σ z 0 and is simply proportional to the focal length times the energy spread, assuming that thermal contributions dominate the longitudinal emittance. However, we see that the non-linear correlations imparted by the buncher significantly distort the trace space and dominate the emittance in the final drift. Ultimately, beams with smaller longitudinal emittance enable reaching shorter bunch lengths. It also follows that one can achieve proportionally shorter final bunch durations by decreasing the focal length f of the RF buncher.

2.3. Single Particle Dynamics and Non-Linear Phase-Space Correlations in the RF Buncher

Let us now look more closely at the details of the energy chirp imparted by the RF buncher on the beam distribution. The main assumption here is that the cavity fields act on the electrons by adding an energy kick with sinusoidal dependence on the initial longitudinal particle position z 0 as
Δ γ = α sin ( k z 0 ) ,
where α = e V 0 / m c 2 and e V 0 is the cavity voltage or the maximum energy gain seen by an ideally phased particle, and k = k 0 / β is the RF angular wave number divided by the normalized longitudinal velocity. The phase of the cavity is tuned so that the center of the bunch experiences no net energy gain and particles at the tail gain energy, while particles at the head of the bunch lose energy. There are two distinct sources of non-linearities in the trace space dynamics resulting from the applied energy change to the particles. First, for finite duration input bunches, the curvature of the RF wave causes significant non-linear effects in the trace space. In addition, the relativistic relation between normalized velocity and beam energy β = 1 1 / γ 2 adds an important degree of non-linearity to the transport as pointed out in Zeitler et al. [16]. Following the discussion therein, Taylor-expanding the relative velocity deviation Δ β β in terms of the energy deviation, we can write
Δ β β = m η m Δ γ m ,
where η m are proportional to the mth derivatives d m d γ m β and in particular
η 1 = 1 β 2 γ 3 ,
η 2 = 2 3 γ 2 2 γ 6 β 4 ,
η 3 = 2 5 γ 2 + 4 γ 4 2 γ 9 β 6 ,
where γ and β are the mean values of the normalized energy and velocity distributions, respectively. Coefficients η m scale as γ ( m + 2 ) so that at high relativistic energies, the higher-order non-linear terms in the transport can be neglected. Considering the lowest-order dynamics, we can simply replace sin ( k z 0 ) with k z 0 ( k z 0 ) 3 6 , and truncate the series to obtain
Δ β β η 1 α ( k z 0 ) + η 2 α 2 ( k z 0 ) 2 + η 1 α 6 η 3 α 3 ( k z 0 ) 3 .
We verify this expression at high energy (4.6 MeV) and low energy (150 keV) by considering a particle tracking simulation of the buncher configuration listed in Table 1 with an initial bunch length of 195 μm and 1.87 mm, respectively. The buncher is modeled by a 2.856 GHz standing wave cylindrically symmetric TM010 cavity, supporting the first transverse magnetic standing wave mode with an axial electric field component for particle acceleration, with an amplitude adjusted to reach a longitudinal focus 1.88 m and 1 m downstream, respectively, for the high- and low-energy cases. Since the RF cavity length is 0.05 m, it is reasonable to approximate it as a thin lens. The longitudinal phase spaces from GPT at the exit of the buncher are shown in Figure 2a,c for the high- and low-energy cases, respectively, with subtracted linear correlations. The quality of the agreement between GPT and our analytical framework can be assessed by comparing the distributions with the lines corresponding to Equation (11) which are also shown. The parameters chosen for these examples highlight the different possibilities for the dominant non-linearity in the system. In the high-energy case, the relativistic effects are responsible for the parabolic shape seen in the simulation, while in the low-energy case, the injected bunch length is longer, and the third-order non-linearity associated with the sinusoidal RF fields is the main effect in the beam distribution shape.
The convenience of working in the trace space is the linearity of the dynamics in the drift which fully preserves the trace space area. Explicitly, in the drift after the buncher, the longitudinal particle position can be written as
z = z 0 + s Δ β β = z 0 + s n = 1 η n Δ γ n .
The initial coordinate z 0 is expressed in terms of the induced energy modulation Δ γ by inverting Equation (8). Then, Δ γ is Taylor expanded in terms of Δ β / β . Substituting into Equation (12), keeping only terms up to the third order, we can write
z s 1 η 1 α k Δ β β + η 2 η 1 3 α k Δ β 2 β 2 η 1 α 6 ( η 3 2 η 2 2 / η 1 ) α 3 η 1 4 α 4 k Δ β 3 β 3 .
The longitudinal waist occurs where the linear chirp is cancelled at distance s = 1 η 1 α k along the beamline, allowing us definition of the buncher longitudinal focal length
f = 1 η 1 α k = m 0 c 2 γ 3 β 2 e V 0 k ,
which indicates that very high voltage cavities are needed to obtain short focal lengths for relativistic electrons. It is also useful to note the k-dependence of this expression which favors the use of very high frequencies for this application.
At the focal plane, the residual correlation is quadratic or cubic in Δ β depending on the relative importance of the non-linearity in the drift propagation concerning the RF curvature. As discussed above, lower beam energies and longer input bunches tend to show higher third-order non-linearities, while relativistic energies typically have dominant second-order contributions. Predictions from Equation (13) can be again verified by comparison to the phase spaces at the temporal waist plane from the same GPT simulation, as shown in Figure 2b,d.

2.4. Emittance Growth Mechanisms and the Relationship between Different Longitudinal Phase Space Definitions

2.4.1. ( z , z ) Trace Space Emittance

Since the drift dynamics in the trace space are entirely linear, and the emittance growth is all accrued in the buncher, the envelope equation formalism is a convenient choice to follow the RMS bunch length evolution. To evaluate the RMS emittance growth induced by the RF compressor, we start from an initial longitudinal phase space with RMS emittance ϵ z 0 , z 0 . After the energy chirp is applied, the single-particle velocity variation maps to z 0 z 0 + Δ β / β , where Δ β / β represents the velocity variation imparted by the buncher, which is correlated with particle position. In the thin lens approximation, the particles do not change position as the beam proceeds through the cavity.
The moments of the new distribution can be calculated and the relation between initial and final emittance after the buncher written as
ϵ z z 2 = ϵ z 0 z 0 2 + ϵ R F 2 = ϵ z 0 z 0 2 + z 0 2 Δ β β 2 z 0 Δ β β 2 .
We note that each respective expectation value needed to calculate the emittance, with the exception of z 0 2 = σ z 0 2 , requires integrating sin m ( k z ) or z sin m ( k z ) over the beam distribution. Assuming an initial Gaussian current profile, these integrals have closed-form expressions up to arbitrary order of m, but keeping the leading contributions to the emittance growth, we obtain
ϵ R F 2 σ z 0 2 2 η 2 2 α 4 k 4 σ z 0 4 + 1 6 η 1 α 6 η 3 α 3 2 k 6 σ z 0 6 .
In most cases, this expression is much larger than the initial longitudinal emittance because the buncher is fundamentally inducing a large velocity spread (with correspondingly significant non-linear contributions) to achieve strong compression. As long as the space-charge forces are negligible in the drift, ϵ R F is also equal to the final emittance and can be used to calculate the shortest bunch length achievable at the waist using Equation (7). The emittance growth in the buncher and its preservation in the drift are shown in Figure 3, where the trace space RMS emittance evolution calculated from GPT is plotted along the beamline. The black dotted line shows Equation (15), which provides a good approximation for the final emittance after a thin lens buncher cavity located at the origin.

2.4.2. ( z , δ ) Phase Space

It is important to note at this point that if instead we utilized more common choices of defining the longitudinal trace space in terms of the relative energy spread δ = Δ γ γ or momentum spread Δ p z / p z , the drift dynamics would become highly non-linear, especially for mildly relativistic particles.
Initially, the trace space emittance ϵ z , z is related to the ( z , δ ) emittance by the following relationship:
ϵ z 0 , δ = γ 2 β 2 ϵ z 0 , z 0 .
This relationship holds because z η 1 Δ γ and explains the order of magnitude difference in absolute emittance values in Figure 3. It is then fairly common to see in the literature the envelope equation written in terms of the ( z , δ ) emittance:
σ z = z z σ z + ϵ z , δ 2 β 2 γ 2 σ z 3 .
Nevertheless, in a drift, the particle positions evolve according to Equation (12) where the transport is inherently non-linear, especially for mildly relativistic particles, causing emittance growth and limiting the usefulness of the envelope equation approach. Ultimately, for large enough initial energy spreads, the higher-order terms proportional to η m lead to significant ( z , δ ) emittance growth, which can be estimated using the same techniques as in the previous subsection:
ϵ z δ 2 = ϵ z 0 , δ 0 2 + 2 s 2 η 2 2 α 6 k 6 σ z 6 ,
where ϵ z 0 , δ 0 is the emittance at the beginning of the drift. This expression predicts a nearly linear growth with propagation distance for a small initial emittance. This is shown in Figure 3 where the ( z , δ ) phase space evolution from GPT is compared with Equation (18) with the inclusion of the initial emittance as well.
The seemingly counterintuitive behavior of the longitudinal emittance (growing linearly in the drift) is the main reason we adopt the ( z , z ) trace space emittance in calculating the final bunch length when using the envelope equation formalism. Finally, for completeness, we observe that if the un-normalized momentum ( z , Δ p z / p z ) was used as a trace space variable, all expressions could be simply modified substituting z = 1 γ 0 2 Δ p z / p z . Nevertheless, due to the relativistic non-linear relation between momentum and velocity, even in this case, we would have significant emittance growth in drift propagation.

2.5. Bunch Length Limit in Absence of Space-Charge Effects

This formalism clarifies how the minimum achievable bunch duration depends on the main beamline parameters (when space-charge effects can be neglected). Inserting the trace space longitudinal emittance estimate (Equation (18)) into the envelope-equation solution (Equation (7)), we obtain an expression written as the quadrature sum of the different contributions to the final emittance between (i) the initial uncorrelated relative energy spread σ δ , (ii) the non-linearities introduced by the relativistic correction to the transport or (iii) the RF-induced emittance.
In order to facilitate comparison between beamlines of different energies, we can rewrite the terms as a function of the buncher’s focal length, which is a practically helpful parameter related to standard requirements on dimensions of the sample chamber, pumping geometry, and transverse optics. Assuming that one term in the sum is much larger than the others, we can synthesize this result as
σ z f max f σ z 0 = f β 2 γ 2 σ δ 2 | η 2 | η 1 α k σ z 0 2 3 2 γ 2 2 σ z 0 2 f 1 6 k 2 σ z 0 3 .
We note that the total quadrature sum of these three expressions should be used in those cases where two or more contributing terms have a similar magnitude.
In the first case, non-linearities in the transport can be neglected, and the final bunch length is simply proportional to the initial relative velocity spread σ z 0 . For the same uncorrelated relative energy spread σ δ , a relativistic energy system has a clear advantage to achieve ultrashort bunch lengths due to the inverse square γ dependence in the formula. In reality, as the energy increases, the non-linearities due to the relativistic dynamics in the drift would likely become the dominant contribution to the final emittance and bunch length. In this regime, after approximating η 2 / η 1 2 3 / 2 γ 2 and using the definition of the focal length (Equation (14)), the final bunch duration scales as the square of the beam energy. For longer initial bunch lengths, the curvature of the RF dominates the shape of the final phase space and sets the limit for the shortest bunch duration achievable. This contribution is essential for both relativistic and non-relativistic energies, and it strongly favors the use of lower RF frequencies. In all THz compression experiments carried out so far, this term has been the principal limit to the final bunch length [12,20].
The analytical formulas summarized in Equation (19) were found in excellent agreement with GPT in both the relativistic and non-relativistic regimes. The results are shown in Figure 4, where (a) is corresponds to the 4.6 MeV and (b) the 150 keV cases, respectively. In the low-energy case, the cubic non-linearity from the RF curvature dominates due to the longer initial bunch length.
In the absence of space charge, the shortest bunch durations are obtained by minimizing the emittance growth in the temporal lens, which can be achieved by using shorter input beams. At the same time, as we see further, decreasing the initial bunch length also increases the initial peak current so that the space-charge effects in the longitudinal envelope equation can no longer be neglected. As a result, space-charge effects begin to take over at a certain point and prevent further bunch compression. Thus, there must be an optimum initial bunch length to inject, which exactly balances space charge and RF emittance growth.

3. Space-Charge Limits to Compression

3.1. An Example of Geometry Factor Calculation: Gaussian Distribution Case

In order to add the effect of space charge on the bunch length evolution into the envelope equation formalism, we need to compute a reasonable representation for the z F s c ( z , σ z , σ ) term.
For simplicity, we derive the self-field of an azimuthally symmetric 3D Gaussian beam assuming that the distribution function remains a Gaussian profile throughout evolution. Admittedly, this approximation is especially poor at the waist where the linear chirp is removed, and what remains is a non-linear distribution in trace space with a characteristic current spike and a temporal profile strongly asymmetric and far from a regular Gaussian. In this situation, the emittance growth caused by the higher-order terms in the electric field profile also makes the envelope equation approach less helpful in describing the bunch evolution, and we would need to revert back to self-consistent particle tracking simulations.
Nonetheless, the simple expression for the force obtained below can at the very least be used to estimate the evolution of the bunch length and the point in the system where space-charge effects become dominant in the dynamics. To validate all of the estimates made in what follows, we compare the results with simulations of the bunching process utilizing GPT’s spacecharge3Dmesh algorithm for the nominal high-energy beam parameters shown in Table 1 and a beam charge of 16 fC (10 5 electrons).
The charge density in the beam rest-frame can be written as
ρ = Q exp r 2 2 σ r 2 z 2 2 σ z 2 ( 2 π ) 3 / 2 σ r 2 σ z .
In order to obtain the electric field components, we Fourier transform the charge density and write the potential as
ϕ ( x , y , z ) = ρ ˜ ( k x , k y , k z ) ϵ 0 ( k x 2 + k y 2 + k z 2 ) exp ( i k · r ) d 3 k ( 2 π ) 3 ,
where ρ ˜ is the Fourier transform of the charge density which also has a Gaussian shape.
We can then expand the complex exponential as a Taylor series in k · r . After dropping the constant term, which is immaterial for the field profile, we also note that all the odd terms vanish by the symmetry of ρ ˜ . The second-order term yields the uncorrelated linear electric field components. They can be written as
E ( 1 ) = Q / 2 σ z ( 2 π ) 3 / 2 ϵ 0 σ 2 u ( A ) r r ^ + Q ( 2 π ) 3 / 2 ϵ 0 σ z 3 v ( A ) z z ^ ,
where A = σ / σ z is the beam aspect ratio (in its rest frame). The predicted field gradients are found in excellent agreement with core field gradients extracted from GPT simulation as shown in the example in Figure 5a. The geometry factors (plotted for reference in Figure 5b are given by
u ( A ) = ξ ( A ) ( 1 ξ ( A ) 2 ) coth 1 1 ξ ( A ) ξ ( A ) 3 ,
v ( A ) = coth 1 1 ξ ( A ) ξ ( A ) ξ ( A ) 3 ,
where ξ ( A ) = 1 A 2 .
The generalized force term in the envelope equation can be calculated using z = γ / γ 3 β 2 = e E z / γ 3 β 2 m c 2 and then taking the average over the beam distribution of z E z . Before doing so, we express the longitudinal electric field just calculated in the beam rest frame in terms of laboratory frame quantities. This can be achieved by rescaling the longitudinal coordinate and RMS moments by γ , i.e, E z ( z , σ z , σ ) E z ( γ z , γ σ z , σ ) . Finally, the space-charge term in the envelope equation is given by
z F s c ( z , σ z , σ ) σ z = g ( A ) N r c β 2 γ 5 σ z 2 = K L σ z 2 ,
where g ( A ) = 1 2 π v ( A ) is an order of unity factor which takes into account the geometry effects, r c is the classical electron radius and N the number of electrons in the bunch; thus, we define the longitudinal perveance K L as anticipated in the first section. In Figure 5a, we compare the result of averaging the longitudinal electric field over the beam distribution with the on-axis linear component, showing how averaging reduces the slope of the z , E z correlation by the constant term 2 / 4 .
In general, the dependence of the geometry factor g ( A ) on the beam aspect ratio requires self-consistently solving the transverse and longitudinal envelope equations as a coupled system, which typically can be achieved through numerical integration.
An attractive simplification occurs if we use the approximation that the aspect ratio remains constant (or nearly constant) throughout the beamline, which corresponds to the time when the beam is simultaneously transversely focusing and longitudinally compressing. In this case, the system is decoupled, and we can then utilize g ( A ) to integrate the longitudinal envelope equation analytically. This approximation turns out to be quite acceptable, as in many UED setups it is required to have a small transverse spot at the sample in the same plane of the temporal waist.
In Figure 6, we plot the evolution of the bunch length and transverse spot along the beamline in the case where optics are arranged so that the aspect ratio during compression is kept close to unity. For this case, we also use GPT simulations to numerically compute the exact expression for the space charge term in the envelope equation at any given location along the beamline and find good agreement with the analytical expression in Equation (25) with A = 1 .
Importantly, when transverse focusing is not applied (or weak) as the beam reaches its minimum bunch length, the aspect ratio increases to values larger than unity along the beamline. Indeed, Figure 5b suggests that our approximation would lead to overestimating the space-charge force.

3.2. Effect of the Longitudinal Space-Charge Force on the Minimum Bunch Length

Using these results, the longitudinal envelope equation in presence of the space charge can be written as
σ z = K L σ z 2 + ϵ z , z 2 σ z 3 ,
and following the same steps that led to Equation (7), we can directly integrate to calculate the bunch length at the waist, obtaining
σ z f = ϵ z , z 2 σ z 0 2 ϵ z , z 2 f 2 + K L 2 K L ,
where we once again assume that the initial bunch length is much larger than the final one.
In the limit that K L becomes negligible, Equation (27) yields back the zero space-charge solution discussed in the previous section. Conversely, if space charge dominates the bunch length evolution, one can expand the formula for large K L to obtain
σ z f = 2 f 2 K L σ z 0 2 ,
which is linear in the bunch charge but most importantly decreases as σ z 0 2 , driving the initial bunch length towards larger values. We note that if by some clever scheme (for example, by pre-compensating using an X-band RF cavity as discussed in the last section of the paper) the RF emittance growth in the buncher is eliminated, a longer initial bunch length could be used, reducing space-charge effects and potentially allowing for reaching very short bunch lengths.
The optimal initial conditions are to compromise between the tendency to minimize the buncher’s RF emittance growth and lowering the initial peak current. An estimate for the ideal initial bunch length can be obtained by equating the asymptotic dependencies of the waist size in the space-charge dominated regime and the emittance dominated regime, respectively.
For example, in the case where the RF emittance is dominated by the quadratic non linearity, we can set Equation (28) equal to the second expression in Equation (19), and the optimum initial bunch size is given by
σ z 0 4 2 K L η 1 | η 2 | α 3 k 3 .
This result can be substituted into Equation (27), yielding for the minimum bunch duration at the temporal waist:
σ z f 2.72 K L | η 2 | η 1 3 α k .
In the high energy limit ( γ 0 1 ) ,
σ z f 3 g N f r c γ 0 3 .
In Figure 7a, we show the final waist as a function of the injected bunch length obtained from GPT, and overlay Equation (27) for the case when the total charge is Q = 10 5 e = 16 fC. The dashed red line is the estimate obtained from Equation (30), which predicts an RMS size of 0.8 μm, underestimating by 20% the GPT prediction at 1 μm. For large enough initial bunch length, space-charge effects are negligible, and the emittance growth in the buncher dominates the waist size. When the initial bunch length is smaller than the optimum, i.e., in the space-charge dominated regime, the minimum bunch length increases as predicted by Equation (28). In Figure 7b, the phase space at the waist for the optimal injection case is shown. The quadratic correlation in the phase space is still visible, but an imprint of the space-charge field also appears.
The slight discrepancy between the simulation and the analytical prediction is due to the unaccounted space-charge-induced emittance growth. As we shorten the initial bunch length, the non-linearities in the space charge field are responsible for significant emittance growth, thus disrupting the applicability of the analytical result, which relies on the assumption of constant emittance in the drift. This is elucidated in Figure 7c,d, where we show the trace space emittance evolution for a few different compression cases and the final emittance as a function of input bunch length, respectively. In (c), it can be seen that the longitudinal emittance after the buncher is changing in the drift due to non-linear space-charge forces. The increase is modest when the input size is sufficiently large and more pronounced for smaller initial bunch lengths. In (d), the longitudinal emittances calculated from Equation (15) are compared with the final emittances at the waist from GPT. The emittance growth induced by the space charge is proxied by subtracting the RF growth from the emittance at the focus position in quadrature. Therein, it can be seen that the optimum input bunch length (72 μm, in this case) occurs at the onset of the space-charge-induced emittance growth. Eventually, when the initial bunch length becomes too short, the emittance growth is dominated by the space charge, so the GPT simulation results are not identical to the analytical predictions.
In Figure 8a, the comparison of Equation (30) with GPT simulations performed while setting the initial bunch length to Equation (29) as the charge is varied from 10 5 e 10 6 e shows good agreement when using a scaled longitudinal perveance to take into account the emittance growth from non-linear space-charge forces. In Figure 8b, we utilize Equation (30) to visualize the dependence on energy and charge when the focal length of the bunching system is set to 1.88 m. Due to the energy dependence intrinsic in the focal length, as we increase the beam energy, keeping the focal length constant becomes a technological feat involving considerations of the breakdown limit in RF cavities and available power sources at higher frequencies.
We note that the optimal initial bunch length identified above is the RMS spot size at the entrance of the buncher. Depending on the gun technology employed to generate the electron beam, there might be different ways to tune this quantity (i.e., changing the laser pulse length, operating at different phase or gradient in an RF gun, controlling the transverse dynamics), which might have an additional effect on the actual beam emittance or phase space correlations.
To move beyond these limits, we would need to reduce the space-charge effects (working to either minimize the geometry factor g in the longitudinal perveance or utilize a charge distribution with more linear self-fields) or minimize the emittance growth in the buncher, for example, pre-compensating for the non-linearities in the input phase space. In principle, both options are feasible and are the subjects of the last two sections of this paper.
We should also observe that the presented solution of the envelope equation by direct integration does not yield the position of the longitudinal waist along the beamline, which is an essential experimental parameter in the case where space-charge forces significantly contribute to the dynamics. The space charge pressure pushes the waist position downstream of the zero charge case, with a shift increasing for more significant beam charges, but that remains relatively small for the optimal injection case. Looking at the simulation for the case where the input bunch length is 72 μm, the longitudinal waist is found at 1.911 m from the buncher, 3.1 cm downstream of the zero charge location.

4. Bunch Compression Limits for Different Charge Distributions

So far, we assumed a Gaussian temporal profile for the input electron bunch. It is worth investigating the prospect of an ideal phase space such as the uniformly filled ellipsoid with transverse and longitudinal dimensions a and z m , respectively. The line charge density for this distribution is
λ ( z ) = 3 Q 4 z m 1 z 2 z m 2 ,
where Q is the bunch charge, z m is related to the RMS bunch length of the beam as z m = 5 σ z , and likewise a = 5 σ . Repeating the calculation from Section 2.4.1 of the RF emittance at the exit of the buncher with this line charge density yields
ϵ R f 2 σ z 0 2 8 7 η 2 2 α 4 k 4 σ z 0 4 + 50 1323 α η 1 6 α 3 η 3 2 k 6 σ z 0 6 .
This expression is similar to what we previously calculated in Equation (15), but with smaller coefficients (respectively, by 55% and 22%) since, for the same RMS bunch length, the parabolic current profile extends over a smaller interval of RF phases than the long-tailed Gaussian distribution, thus reducing the associated RMS emittance growth.
Reiser offers the perveance and geometry parameters of the longitudinal field evolution [18]. The longitudinal field is given by
E z = v ( A ) 2 π ϵ 0 λ z .
Then, the generalized longitudinal force derived from this field yields the perveance
K L , u = g u ( A ) N r c β 2 γ 5 ,
where now, g u ( A ) = 3 v ( A ) / 5 5 . The prefactors of the RF emittance growth and the perveance are smaller than those obtained for Gaussian, so this distribution yields substantial improvements to the optimum bunch length. To this end, we can use the perveance factor and again use the same steps discussed in the previous section to estimate the optimum initial bunch length and the shortest final duration of the beam.
We repeat, for this case, the benchmarking simulation study by initializing in GPT a uniformly filled ellipsoidal distribution with optimal input bunch length and compare the results with the analytical predictions as a function of beam charge as shown in Figure 8a. The agreement is excellent, and it is worthwhile to note that in this case, there is no need to rescale the longitudinal perveance as the space-charge forces are more linear than the Gaussian, so there is less induced emittance growth. Compared to the Gaussian current profile, this distribution improves the optimum bunch length by nearly two.

5. X-Band Cavity Compensation

5.1. Analytical Estimates

An alternative option to obtain even shorter bunch lengths is using an additional higher-frequency cavity to compensate for the non-linearities imparted on the longitudinal trace space by the buncher RF fields. The use of a higher frequency linearizing cavity has been proposed originally in the context of RF compression by Floettmann [15], and is also very common in FEL beamlines [21] to compensate for the non-linearities in magnetic compressors.
In this section, we consider a setup closely mimicking the current configuration of the UCLA Pegasus beamline. The setup consists of a 1.6 cell S-band (2.856 GHz) gun, an X-band (9.6 GHz) cavity situated at 1.1 m from the cathode plane, and an S-band (2.856 GHz) 11 cell linac. The Linac has a phase chosen to velocity bunch the beam. The linac entrance is 1.4 m downstream of the cathode plane. In Figure 9, we display an illustration of the beamline.
An analytical expression for the compensation condition can be found simply imposing the cancellation of the second-order coefficients in the dependence of the relative velocity from the initial longitudinal position yielding
α x = 2 | η 2 | η 1 k k x 2 α 2 ,
where we assume the linearizing cavity to operate at 180 degrees from the crest (i.e., Δ γ x = α x cos ( k x z 0 ) ) and α x and k x are the normalized voltage and wavenumber of the X-band cavity, respectively. In order to minimize the amount of beam deceleration and the power requirements for the linearizer, it helps to maximize the ratio between the linearizer and buncher cavity frequencies. In addition, for a fixed focal length of the bunching system, α x γ 5 ; thus, lowering beam energy reduces the required voltage. By adding the phase of the linearizer as a degree of freedom, we could simultaneously cancel second- and third-order distortions.
Assuming that the x-band cavity fully compensates the second-order contribution to the emittance and that the cubic RF non-linearity is retained in Equation (27), we can obtain analytical estimates for the minimum achievable bunch length following the same process presented in the previous section. In Figure 10, we plot the optimal solutions for beams of charge 160 fC with Gaussian and uniformly fill ellipsoidal distributions along with the analytical prediction for the uncompensated cases (dashed). The addition of the linearizer cavity allows for achieving bunch lengths significantly shorter (by up to a factor of five) than the single cavity case. The uncorrelated energy spread sets the thermal limit on the final bunch length even when the space charge is on. In the calculations, we set the uncorrelated energy spread, σ δ 10 5 , and for convenience, we show as a dashed line the corresponding thermal limit (0.13 μm).
These results are validated using GPT simulations where we start from an initially uncorrelated coasting beam injected into a beamline with the linearizer cavity and the buncher. A solenoid keeps the aspect ratio between 0.5 and 0.8 as the beam is compressing. The uniformly filled ellipsoidal density analytical prediction matches the simulation results well, while the analytical curve for the Gaussian distribution again underestimates the simulation results due to the non-linear space-charge-induced emittance growth.
In the region along these curves where the space-charge effects are more relevant (i.e., left of the minimum), the final bunch length should increase linearly with charge according to the longitudinal perveance expression from Equation (28). We verify this by simulating an initial bunch length of 180 μm for multiple beam charge between 16 fC and 160 fC. The results are shown in Figure 11 and match well with the analytical predictions. The intercept of the line is due to the finite initial emittance due to thermal and transverse effects [22].

5.2. Start-to-End Simulations

So far, our discussion has assumed an initially flat and uncorrelated longitudinal phase space. However, in a realistic system, the beam out of the gun would typically present already some correlations (often of non-linear nature) between energy and time along the bunch. To analyze the effects of these, we present start-to-end simulations of the entire beamline, including the RF photo-injector section with nominal parameters listed in Table 2.
The results for minimum waist size as a function of the laser pulse length illuminating the cathode are shown as red squares in Figure 10 and establish a good agreement with the analytical predictions. However, the start-to-end simulations depend on many parameters, such as gun solenoid setting, initial laser spot size, cathode thermal emittance, etc. Consequently, the transverse phase space at the entrance of the RF bunching section was not perfectly matched to the idealized simulations. As a result, the aspect ratio at the focus was 2.8 instead of 1, so the square red dots fall slightly below the analytical estimates, and a minimum bunch length of 0.28 μm or 940 as can be reached.
In Figure 9, we display the simulated phase spaces at various positions along the beamline after each cavity. Figure 12 elucidates the dependence upon the injected bunch length. In the optimal case (center), the initial thermal emittance and the emittance growth imparted by the space charge fields limit the final bunch length. The third-order distortions limit the bunch length if the initial bunch length is too long (right). On the other hand, if the initial injection becomes too short (left), then the final phase space is diluted further by transverse dependence on the longitudinal field, and the focus position significantly moves from the waist plane of the zero space charge case.
Finally, we note that while this scheme provides a direct path to sub-fs bunch lengths, the synchronization of the drive signals for two different resonant frequency RF cavities is a significant technological challenge. However, addressing this challenge is the first step before achieving high-quality compensation in a user-facility UED beamline. Therefore, the case is presented here mainly as an example of the insights offered by the analytical framework that we developed, which suggests that properly shaping the drive laser and compensating the emittance growth does provide a path towards much shorter bunch lengths.

6. Conclusions

In conclusion, in this paper, we developed an analytical framework that yields a simple estimate for the minimum bunch length achievable in an RF compression beamline. Besides UED, the formulas in this paper might be helpful in the optimization of RF compression for other ultrashort electron beam applications, including radiation generation and injection into very high-frequency THz-driven [23] and laser-driven accelerators [24].
The envelope formalism allowed us accurate estimation of the final bunch length; our analysis showed that we can integrate the competing non-linear effects into the envelope evolution accurately with relative ease. Specifically, we evaluated longitudinal emittance growth caused by RF curvature and relativistic beam transport, then developed an excellent approximation to the space-charge force in the envelope equation. The results indicate that an optimum initial bunch length condition (which can be satisfied by adjusting the laser pulse length on the cathode, for example) compromises the RF-induced emittance growth and the effects of the longitudinal self-fields. The simplicity of the reported expressions mainly stems from the fact that we approximated the coupling of the transverse and longitudinal space charge dynamics with a simple constant order-of-unity geometry factor in the longitudinal perveance. We also limited the expansion yielding the non-linear terms in the emittance growth to second and third orders. Still, in principle, we could have evaluated all the higher-order terms (and compensated if enough independent knobs were added/available on the beamline).
Although realistic beams produced by photoinjectors typically present more complex phase-space distributions, the initial conditions assumed in the derivations are an initially unchirped longitudinal phase space. Nonetheless, the results obtained here still provide valuable estimates of compression limits in a given configuration, which prove helpful as a starting point for numerical optimizations. In addition, the scaling laws proved capable of guiding parameter choices in the design of new setups. Most importantly, these results highlighted the main contributions to the final bunch length and suggested possible paths to improve the compression further and achieve sub-fs bunch lengths.

Author Contributions

Conceptualization and Methodology: P.D. and P.M. Validation: P.D. validated expressions with particle tracking code. Writing—Original Draft Preparation: Both P.D. and P.M. contributed to writing the main text. Visualization: P.D. generated all figures. Supervision: P.M. supervised the development of the paper. Funding Acquisition: P.M. acquired funding. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially funded by the National Science Foundation (NSF) (PHY-1734215); National Science Foundation under the STROBE Science and Technology Center Grant No. DMR-1548924 and the STROBE Consortium MRI grant DMR-1828705.

Acknowledgments

The authors want to acknowledge D. Filippetto, A. Kogar and J. Maxson for insightful discussions leading to the calculations in this paper over the years.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zewail, A.H. 4D Ultrafast Electron Diffraction, Crystallography, and Microscopy. Annu. Rev. Phys. Chem. 2006, 57, 65–103. [Google Scholar] [CrossRef] [PubMed]
  2. Sciaini, G.; Miller, R.D. Femtosecond electron diffraction: Heralding the era of atomically resolved dynamics. Rep. Prog. Phys. 2011, 74, 096101. [Google Scholar] [CrossRef]
  3. Siwick, B.J.; Dwyer, J.R.; Jordan, R.E.; Miller, R.D. An atomic-level view of melting using femtosecond electron diffraction. Science 2003, 302, 1382–1385. [Google Scholar] [CrossRef] [PubMed]
  4. Weathersby, S.; Brown, G.; Centurion, M.; Chase, T.; Coffee, R.; Corbett, J.; Eichner, J.; Frisch, J.; Fry, A.; Gühr, M.; et al. Mega-electron-volt ultrafast electron diffraction at SLAC National Accelerator Laboratory. Rev. Sci. Instrum. 2015, 86, 073702. [Google Scholar] [CrossRef] [PubMed]
  5. Zhu, P.; Zhu, Y.; Hidaka, Y.; Wu, L.; Cao, J.; Berger, H.; Geck, J.; Kraus, R.; Pjerov, S.; Shen, Y.; et al. Femtosecond time-resolved MeV electron diffraction. New J. Phys. 2015, 17, 063004. [Google Scholar] [CrossRef]
  6. Manz, S.; Casandruc, A.; Zhang, D.; Zhong, Y.; Loch, R.A.; Marx, A.; Hasegawa, T.; Liu, L.C.; Bayesteh, S.; Delsim-Hashemi, H.; et al. Mapping atomic motions with ultrabright electrons: Towards fundamental limits in space-time resolution. Faraday Discuss. 2015, 177, 467–491. [Google Scholar] [CrossRef] [PubMed]
  7. Filippetto, D.; Qian, H. Design of a high-flux instrument for ultrafast electron diffraction and microscopy. J. Phys. At. Mol. Opt. Phys. 2016, 49, 104003. [Google Scholar] [CrossRef]
  8. Van Oudheusden, T.; De Jong, E.; Van der Geer, S.; ’t Root, W.O.; Luiten, O.; Siwick, B. Electron source concept for single-shot sub-100 fs electron diffraction in the 100 keV range. J. Appl. Phys. 2007, 102, 093501. [Google Scholar] [CrossRef]
  9. Van Oudheusden, T.; Pasmans, P.; Van Der Geer, S.; De Loos, M.; Van Der Wiel, M.; Luiten, O. Compression of subrelativistic space-charge-dominated electron bunches for single-shot femtosecond electron diffraction. Phys. Rev. Lett. 2010, 105, 264801. [Google Scholar] [CrossRef] [PubMed]
  10. Maxson, J.; Cesar, D.; Calmasini, G.; Ody, A.; Musumeci, P.; Alesini, D. Direct measurement of sub-10 fs relativistic electron beams with ultralow emittance. Phys. Rev. Lett. 2017, 118, 154802. [Google Scholar] [CrossRef] [PubMed]
  11. Zhao, L.; Wang, Z.; Lu, C.; Wang, R.; Hu, C.; Wang, P.; Qi, J.; Jiang, T.; Liu, S.; Ma, Z.; et al. Terahertz streaking of few-femtosecond relativistic electron beams. Phys. Rev. X 2018, 8, 021061. [Google Scholar] [CrossRef]
  12. Snively, E.; Othman, M.; Kozina, M.; Ofori-Okai, B.; Weathersby, S.; Park, S.; Shen, X.; Wang, X.; Hoffmann, M.; Li, R.; et al. Femtosecond compression dynamics and timing jitter suppression in a THz-driven electron bunch compressor. Phys. Rev. Lett. 2020, 124, 054801. [Google Scholar] [CrossRef] [PubMed]
  13. Kealhofer, C.; Schneider, W.; Ehberger, D.; Ryabov, A.; Krausz, F.; Baum, P. All-optical control and metrology of electron pulses. Science 2016, 352, 429–433. [Google Scholar] [CrossRef]
  14. De Loos, M.; Van der Geer, S. General Particle Tracer: A new 3D code for accelerator and beamline design. In Proceedings of the 5th European Particle Accelerator Conference, Sitges, Barcelona, 10–14 June 1996; p. 1241. [Google Scholar]
  15. Floettmann, K. Generation of sub-fs electron beams at few-MeV energies. Nucl. Instrum. Methods Phys. Res. Sect. Accel. Spectrom. Detect. Assoc. Equip. 2014, 740, 34–38. [Google Scholar] [CrossRef]
  16. Zeitler, B.; Floettmann, K.; Grüner, F. Linearization of the longitudinal phase space without higher harmonic field. Phys. Rev. Spec. Top.-Accel. Beams 2015, 18, 120102. [Google Scholar] [CrossRef]
  17. Bacci, A.; Rossi, A.R. Ultra-short electron bunches by Velocity Bunching as required for plasma wave accelerations. Nucl. Instrum. Methods Phys. Res. Sect. Accel. Spectrom. Detect. Assoc. Equip. 2014, 740, 42–47. [Google Scholar] [CrossRef]
  18. Reiser, M.; O’Shea, P. Theory and Design of Charged Particle Beams; Wiley Online Library: Hoboken, NJ, USA, 1994; Volume 312. [Google Scholar] [CrossRef]
  19. Wiedemann, H. Particle Accelerator Physics; Springer Nature: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  20. Zhao, L.; Tang, H.; Lu, C.; Jiang, T.; Zhu, P.; Hu, L.; Song, W.; Wang, H.; Qiu, J.; Jing, C.; et al. Femtosecond relativistic electron beam with reduced timing jitter from THz driven beam compression. Phys. Rev. Lett. 2020, 124, 054802. [Google Scholar] [CrossRef] [PubMed]
  21. Emma, P. X-Band RF Harmonic Compensation for Linear Bunch Compression in the LCLS; Rep. ILCLS-TN-01-1; SLAC: Stanford, CA, USA, 2001. [Google Scholar] [CrossRef]
  22. Bartnik, A.; Gulliford, C.; Hoffstaetter, G.H.; Maxson, J. Ultimate bunch length and emittance performance of an MeV ultrafast electron diffraction apparatus with a dc gun and a multicavity superconducting rf linac. Phys. Rev. Accel. Beams 2022, 25, 093401. [Google Scholar] [CrossRef]
  23. Nanni, E.A.; Huang, W.R.; Hong, K.H.; Ravi, K.; Fallahi, A.; Moriena, G.; Miller, R.D.; Kärtner, F.X. Terahertz-driven linear electron acceleration. Nat. Commun. 2015, 6, 8486. [Google Scholar] [CrossRef] [PubMed]
  24. Cesar, D.; Maxson, J.; Shen, X.; Wootton, K.; Tan, S.; England, R.; Musumeci, P. Enhanced energy gain in a dielectric laser accelerator using a tilted pulse front laser. Opt. Express 2018, 26, 29216–29224. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Illustration of RF ballistic bunching scheme. A velocity chirp is imparted on an electron using an RF cavity so that the tail of the beam has a higher energy than the head. During the following drift, the particles in the tail catch up with the particles in the head, resulting in strong longitudinal compression.
Figure 1. Illustration of RF ballistic bunching scheme. A velocity chirp is imparted on an electron using an RF cavity so that the tail of the beam has a higher energy than the head. During the following drift, the particles in the tail catch up with the particles in the head, resulting in strong longitudinal compression.
Instruments 07 00049 g001
Figure 2. (Left) Trace spaces of the beam at the exit of the prebuncher after the linear chirp is subtracted from the distribution for high-energy (a) and low-energy (c) cases compared with the analytical predictions from Equation (11) (green curves). (Right) Longitudinal trace spaces at the temporal waist for the high-energy (b) and low-energy (d) cases compared with the predictions from Equation (13) (green curves). The current profiles at the focus are also shown in black.
Figure 2. (Left) Trace spaces of the beam at the exit of the prebuncher after the linear chirp is subtracted from the distribution for high-energy (a) and low-energy (c) cases compared with the analytical predictions from Equation (11) (green curves). (Right) Longitudinal trace spaces at the temporal waist for the high-energy (b) and low-energy (d) cases compared with the predictions from Equation (13) (green curves). The current profiles at the focus are also shown in black.
Instruments 07 00049 g002
Figure 3. Comparison between analytical estimates (dashed lines) and GPT simulations (solid lines) for an initial bunch length of 0.65 ps. The trace space emittance after the buncher is shown in blue, and emittance growth in ( z , δ ) phase space is shown in orange.
Figure 3. Comparison between analytical estimates (dashed lines) and GPT simulations (solid lines) for an initial bunch length of 0.65 ps. The trace space emittance after the buncher is shown in blue, and emittance growth in ( z , δ ) phase space is shown in orange.
Instruments 07 00049 g003
Figure 4. Final bunch length as a function of the initial bunch length for the (a) high-energy and (b) low-energy cases. The analytical curves are also shown and are found to be in very good agreement with the simulations.
Figure 4. Final bunch length as a function of the initial bunch length for the (a) high-energy and (b) low-energy cases. The analytical curves are also shown and are found to be in very good agreement with the simulations.
Instruments 07 00049 g004
Figure 5. (a) Longitudinal field of the bunch compared with the linear field component of the 3D Gaussian field and with the beam-averaged generalized force term that appears in the envelope equation. The total charge in the bunch in this simulation is Q = 10 5 e . (b) Geometry factors for the transverse (blue) and longitudinal (purple) field plotted as a function of the rest frame aspect ratio.
Figure 5. (a) Longitudinal field of the bunch compared with the linear field component of the 3D Gaussian field and with the beam-averaged generalized force term that appears in the envelope equation. The total charge in the bunch in this simulation is Q = 10 5 e . (b) Geometry factors for the transverse (blue) and longitudinal (purple) field plotted as a function of the rest frame aspect ratio.
Instruments 07 00049 g005
Figure 6. The evolution of longitudinal and transverse beam sizes are shown in purple. The space charge envelope equation term is shown in green, and its approximation is shown as a dotted line with aspect ratio fixed, i.e., A = 1.
Figure 6. The evolution of longitudinal and transverse beam sizes are shown in purple. The space charge envelope equation term is shown in green, and its approximation is shown as a dotted line with aspect ratio fixed, i.e., A = 1.
Instruments 07 00049 g006
Figure 7. (a) Final waist size plotted versus initial bunch length compared with GPT. The green and black lines show the asymptotic behavior in the regimes where space charge and RF emittance growth dominate the dynamics. The blue line is the analytical expression that considers all effects discussed in this paper (b) The optimum phase space found from the scan in (a) for an input bunch length of 72 μm. The RMS bunch length is 1 μm. Particles are color-coded as a function of their radial coordinate (blue corresponds to on-axis). (c) Evolution of the emittance along the line. The growth observed in GPT is due to the non-linearities of the space-charge field. (d) Final emittance in GPT as a function of initial bunch length. The relative importance of the space-charge contribution to the emittance can be inferred by subtracting the expected emittance growth due to the buncher dynamics non-linearities. The cross-over point (where space charge becomes the dominant effect) can be used to estimate the optimal injection condition and hence the minimum final bunch length achievable for a given setup.
Figure 7. (a) Final waist size plotted versus initial bunch length compared with GPT. The green and black lines show the asymptotic behavior in the regimes where space charge and RF emittance growth dominate the dynamics. The blue line is the analytical expression that considers all effects discussed in this paper (b) The optimum phase space found from the scan in (a) for an input bunch length of 72 μm. The RMS bunch length is 1 μm. Particles are color-coded as a function of their radial coordinate (blue corresponds to on-axis). (c) Evolution of the emittance along the line. The growth observed in GPT is due to the non-linearities of the space-charge field. (d) Final emittance in GPT as a function of initial bunch length. The relative importance of the space-charge contribution to the emittance can be inferred by subtracting the expected emittance growth due to the buncher dynamics non-linearities. The cross-over point (where space charge becomes the dominant effect) can be used to estimate the optimal injection condition and hence the minimum final bunch length achievable for a given setup.
Instruments 07 00049 g007
Figure 8. (a) Minimum bunch length analytical estimates compared with GPT simulation as the charge in the 4.6 MeV energy beam is varied for Gaussian (purple) and uniformly filled ellipsoidal (blue) beam distributions. (b) Minimum bunch length versus γ and N for the Gaussian distribution assuming a constant focal length f = 1.88 m.
Figure 8. (a) Minimum bunch length analytical estimates compared with GPT simulation as the charge in the 4.6 MeV energy beam is varied for Gaussian (purple) and uniformly filled ellipsoidal (blue) beam distributions. (b) Minimum bunch length versus γ and N for the Gaussian distribution assuming a constant focal length f = 1.88 m.
Instruments 07 00049 g008
Figure 9. (a) Illustration of beamline setup for the RF emittance growth compensation. A short x-band cavity is used to compensate for the curvature in ( z , z ) space imparted by the S-band linac and gun and linearize the output trace space. (b) GPT simulation of phase space at the exit of the S-band gun, the exit of the X-band linearizer, the exit of the buncher, and at the longitudinal focus.
Figure 9. (a) Illustration of beamline setup for the RF emittance growth compensation. A short x-band cavity is used to compensate for the curvature in ( z , z ) space imparted by the S-band linac and gun and linearize the output trace space. (b) GPT simulation of phase space at the exit of the S-band gun, the exit of the X-band linearizer, the exit of the buncher, and at the longitudinal focus.
Instruments 07 00049 g009
Figure 10. Analytical predictions for the RMS bunch length at the longitudinal waist for 160 fC charge and 4.6 MeV energy beams plotted with respect to initial bunch length for compensated (solid) and uncompensated (dashed) cases. Space-charge effects are taken into account assuming a constant aspect ratio equal to 0.6 for Gaussian (purple) and uniformly filled ellipsoidal (blue) distribution. These results are compared with GPT simulations of the linearizer beamline in compensation mode. The red, purple, and blue squares show the results of start-to-end simulations for varying laser pulse lengths following the beam from the cathode located in an S-band RF gun.
Figure 10. Analytical predictions for the RMS bunch length at the longitudinal waist for 160 fC charge and 4.6 MeV energy beams plotted with respect to initial bunch length for compensated (solid) and uncompensated (dashed) cases. Space-charge effects are taken into account assuming a constant aspect ratio equal to 0.6 for Gaussian (purple) and uniformly filled ellipsoidal (blue) distribution. These results are compared with GPT simulations of the linearizer beamline in compensation mode. The red, purple, and blue squares show the results of start-to-end simulations for varying laser pulse lengths following the beam from the cathode located in an S-band RF gun.
Instruments 07 00049 g010
Figure 11. Comparison of predicted bunch length dependence on charge for fixed input bunch length. Application to the Gaussian line charge.
Figure 11. Comparison of predicted bunch length dependence on charge for fixed input bunch length. Application to the Gaussian line charge.
Instruments 07 00049 g011
Figure 12. Final phase spaces for a beam shorter than, equal to and longer than the optimal bunch length, respectively. The phase spaces are color coded with respect to the radial coordinate.
Figure 12. Final phase spaces for a beam shorter than, equal to and longer than the optimal bunch length, respectively. The phase spaces are color coded with respect to the radial coordinate.
Instruments 07 00049 g012
Table 1. Simulation beam parameters.
Table 1. Simulation beam parameters.
ParameterHigh EnergyLow Energy
Focal length1.88 m1 m
Beam kinetic energy4.6 MeV150 keV
Norm. transverse emittance100 nm8.3 nm
RMS transverse beam size100 μm100 μm
Cavity Frequency2.856 GHz2.856 GHz
Relative energy spread 10 5 10 5
Table 2. Parameters of start-to-end simluation.
Table 2. Parameters of start-to-end simluation.
ParameterValue
Charge 10 6 e
Laser Spot Size10 μm
Cathode MTE0.5 eV
Optimal laser pulse length0.95 ps (rms)
Gun Accelerating Gradient94.7 MV/m
Gun Phase35.5°
Linearizer accelerating voltage1.8 MV
Linearizer phase173.5°
Buncher accelerating voltage6.75 MV
Buncher phase101°
Final kinetic energy4.5 MeV
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Denham, P.; Musumeci, P. Analytical Scaling Laws for Radiofrequency-Based Pulse Compression in Ultrafast Electron Diffraction Beamlines. Instruments 2023, 7, 49. https://doi.org/10.3390/instruments7040049

AMA Style

Denham P, Musumeci P. Analytical Scaling Laws for Radiofrequency-Based Pulse Compression in Ultrafast Electron Diffraction Beamlines. Instruments. 2023; 7(4):49. https://doi.org/10.3390/instruments7040049

Chicago/Turabian Style

Denham, Paul, and Pietro Musumeci. 2023. "Analytical Scaling Laws for Radiofrequency-Based Pulse Compression in Ultrafast Electron Diffraction Beamlines" Instruments 7, no. 4: 49. https://doi.org/10.3390/instruments7040049

APA Style

Denham, P., & Musumeci, P. (2023). Analytical Scaling Laws for Radiofrequency-Based Pulse Compression in Ultrafast Electron Diffraction Beamlines. Instruments, 7(4), 49. https://doi.org/10.3390/instruments7040049

Article Metrics

Back to TopTop