Physics of Superpulses in Storage Ring Free-Electron Lasers

Contradictory to the widespread perception, storage ring free-electron lasers with substantial net gain can generate peak lasing power reaching GW levels in the so-called superpulse mode. This power level is sufficient for studies of nonlinear processes and efficient intracavity harmonic generation. This letter describes the physics of the superpulses and a phase-space refreshment in the electron beam responsible for this phenomenon.

All presently existing storage ring free-electron lasers (SR FELs) are oscillators with a two-mirror optical cavity and an optical klystron (OK) [1] comprising two wigglers and a buncher between them.SR FELs can operate in pulsed mode [2] with significant peak optical power.Early theoretical models [1] predicted peak intracavity optical power (discussed in this letter), P opt , to be limited by where E o = γ o mc 2 is the central energy and I e is the peak electron beam current, σ γ ind is RMS spread of γ = E/mc 2 induced by FEL interactions, and e is the electron charge.For a typical SR FEL with E o ∼1 GeV, I e ∼10 A, and σ γ ind /γ o ∼10 −3 , (1) limits peak FEL power to about ten MW.Studies of selfconsistent SR FEL dynamics using the 3D uvfel code [3][4][5] revealed a mechanism of superpulses with peak power two-to-three orders of magnitude exceeding the estimate (1).
Figure 1 shows the typical evolution of a superpulse simulated by the uvfel code [3,4], which is known for reliably predicting the electron and the optical beam evolution in SR FELs [2,[7][8][9].
In this letter, we discuss phenomena responsible for superpulses in SR FEL.We present a self-consistent set of dimensionless nonlinear differential equation describing superpulse and the peak power dependencies on key parameters.We also discuss the physics of superpulses using results of computer simulations for illustration (∼1000 runs of uvfel code with ∼10 6 macroparticles and total number of macroparticle passes through the FEL ∼10 11 are used).An ultrarelativistic electron is described in 6D phase space by its energy E e = E o (1+ε), longitudinal coordinate ξ = v e t − z, transverse coordinates and angles X T ≡ [x, x , y, y ], where x ≡ dx/dz and z is the coordinate along the FEL axis.Without lasing, the electron beam obtains a natural distribution in 6D phase space, which is close to a Gaussian and is the product of the normalized transverse (t) and longitudinal (s) distribution functions f 6D ∼ = N e • f t (X)• f s , where N e is the number of electrons in the e-bunch [10].This separation of variables is valid for an SR FEL located in the dispersionfree straight section and processes much shorter than that of radiation damping time, which we discussed in this letter.In the absence of FEL interactions, longitudinal motion is an oscillation [10] δ where Ω s = 2πQ s , Q s is the synchrotron tune, n is the turn number around the ring, and {δ = ε/σ εo , ζ = ξ/σ zo } are dimensionless coordinates scaled by the natural e-beam energy spread σ εo = mc 2 σ γo and the bunch length σ zo .Trajectories in {δ, ζ} phase-space are simple clock-wise circular rotations (see Figure 3).The FEL wigglers (with period λ w and magnetic field B w ) provide a resonant interaction between electrons and the TEM optical wave at wavelength λ o given by (3)  Passing through an OK, an electron radiates a wave packet with length equal to its total slippage σ zo , where N w is the number of periods in one wiggler and N D is the dimensionless slippage in the buncher [11].When FEL interactions are present, the intracavity optical power builds up from spontaneous radiation, which is accumulated inside the optical cavity and amplified during consequent passes by the circulating e-bunch where 1 is the loss per turn in an optical cavity with mirror reflectivity R 1,2 and P SR is the power of spontaneous radiation into the optical mode.At the beginning of the process (at n = 0 FEL interaction is turned on [2]), the longitudinal FEL gain profile is Gaussian, like the e-beam, with peak gain The complete set of dimensionless selfconsistent equations for a superpulse can be derived for a rather general SR FEL case, using Einstein's relations between spontaneous and induced radiation [12][13][14].The optical field in an FEL can be described as a wave packet with slowly varying complex amplitude and a normalized transverse mode function Variation of electron energy in an FEL is the result of local interaction with a Ξ-long segment of the optical field and is given by an integral along the electron's trajectory For a short wavelength FEL, the correlations between electrons are washed away every turn around the ring.Therefore, optical phases φ = k o (ct i − z i ) of electrons entering the FEL are random, resulting in the famous ratio between the energy loss and the energy diffusion [12,13] Using ( 7) and ( 8), one can write a set of dimensionless equations with ν = n•G th as the independent variable, and optical power is normalized as p o (ζ) = P opt (ζ)/P n by P n = √ 2π P e-beam σ γo /γ o which is of the same order as limit in (1).Equation ( 9) contains dimensionless parameters g = G o /G th (gain parameter), r = Ω s /G th (refreshment parameter), and a dimensionless well-defined function normalized as e −δ 2 /2 (dq(δ)/dδ)dδ = −1 to satisfy (5).Function ξ 1 (X,δ,z) above is the electron's slippage with respect to the wave packet Vladimir N. Litvinenko For the optical beam to be amplified, the derivative of the diffusion coefficient q(δ) should be negative at and around δ = 0.It means that q(δ) reaches its maximum at negative δ and its minimum at positive δ (see Figure 3).N D in the range provides for attainment of maximal peak power with q(δ) ∼1, dq(δ)/dδ ∼1 (the reduction of the peak power caused by finite transverse emittances of the e-beam is well described by F x,y factors in [8]).Figure 2 shows dependencies of peak power on the gain and the refreshment parameters.The power is low close to the lasing threshold (g − 1 1), grows very fast in 1.5 < g < 10 range, and saturates above g = 10.Dependence on the refreshment parameter is less dramatic, but peak power increases about sixfold from r 1, when the synchrotron motion is essentially frozen, to r 1. Detailed studies show that optimal ratio between r and g is near π, that is, when electrons go through half of a syn-chrotron oscillation during one e-growth time ν e = 1/g of the laser power.Enhancement comes from synchrotron oscillations moving electrons through the center region and allowing the short optical pulse to be amplified by entire electron bunch.Further increase of the refreshment parameter is not beneficial and can cause oscillations of optical power.
Figure 3 shows details of the e-beam dynamics.During the superpulse, period of accumulation of spontaneous radiation is followed by a period of exponential growth.The correlation length of the optical field grows with the amplification Γ = G FEL (n)dn as ξ corr ∼ = Ξ √ 1 + Γ, while the envelope of the optical pulse shrinks σ opt ∝ 1/ √ 1 + Γ [6] to about 30% of the initial e-bunch length.The initial power growth by 5 to 6 orders of magnitude (see Figures 1 and 3, turns 1-100) does not affect significantly the e-beam distribution and the gain profile.The amplified optical pulse grown from the "white" noise consists of very short (Δζ = 1 − 2•10 −3 ) spikes with peak power in the spikes exceeding the local average by approximately an order of magnitude, similar to that in SASE FELs [15,16].
Further growth of the optical power affects the e-beam distribution and complicates the dynamics.Energy loss and diffusion in the FEL make the phase-space distribution of electrons asymmetric and time-dependent.The growth of the local energy spread reduces the FEL gain while the asymmetric gain profile and energy diffusion cause a drift of the optical pulse and the electron beam from the center [7,18]see changes in δ , ζ in Figure 3.The optical power reaches maximum when net gain crosses zero at G FEL = G th .After that, the optical power decays, and the electron beam distribution becomes more symmetric.
As shown in "Turn 200" plot in Figure 3, lower energy electrons (−3 < δ < 0) have maximum diffusion, while electrons with higher energy have a low diffusion energy corridor clearly noticeable at 1 < δ < 3.This diffusion asymmetry, which is fundamental to the amplification process in FEL [12,13], produces the top-bottom asymmetry in the Poincaré plots.Electrons at large amplitude of synchrotron oscillations (a s ≥ 1) interact with the intense optical pulse briefly (n int ∝ σ opt /a s Ω s ) while passing near the ζ = 0.The diffusion imprints the level of the FEL power at that moment onto the distribution.Thus, the energy spread seen at angle ψ in Figure 3 shows that the FEL power level happened about Δn ∼ = ψ/Ω s turns in the past.Thus, the segment 0 < ψ < π in the turn 200 graph shows the diffusion imprinted at negative δ between turns 147 and 200.Similarly, the segment π < ψ < 2π shows the diffusion imprinted at positive δ during the same turns.Thus, the snail-like dependence of the amplitude spread is a record of the evolution of the optical pulse power (see clip in Figure 1) as well as a clear indication of the diffusion asymmetry.Such snail-like distribution of electron beam was also observed experimentally [19].
The asymmetry in diffusion creates multipole moments in longitudinal phase space and residual quadrupole oscillations with frequency 2Ω s can be seen in the oscillations of RMS energy spread in Figure 1.
The dimensionless nature of (8) provides for the scalability of the results presented in this letter: SR FELs with the same parameters {r, g, q(δ)} would have the same average values and variations of the dimensionless peak power p o .The peak superpulse power can be then simply obtained through multiplication by P n .Thus, the simulation results presented in this letter can be used for predicting parameters in many SR FELs.Overall, the phenomena of superpulses in SR FELs, which are presently verified both experimentally [18,19] and confirmed by other FEL codes [20], open a range of new opportunities for generating coherent femtosecond Xray pulses in SR FELS and their scientific applications.The OK-4 FEL at Duke produced superpulses at wavelengths of 240-270 nm for efficient coherent harmonic generation from 130 nm to 37 nm [9,18].These experiments were followed by successful harmonic generation in Elettra SR FEL [17].
These experiments also support predictions for superpulses presented in this letter [20].

Figure 2 :
Figure 2: Dependencies of the normalized peak power p o on: (a) the refreshment parameter for g = 8.14, (b) the gain parameter g for r = 10.Other parameters are the same as in Figure 1.Each point in the graphs shows the average value and standard error for eight to sixteen statistically independent runs of uvfel code.

Figure 3 :
Figure 3: Four Poincaré plots (electrons are red dots) in (δ, ζ) phase space during the superpulse (parameters are as in Figure 1).The clips show the ζ-profiles of the optical power (blue) with power scales 10 W for turn 1, 10 MW for turn 100, and 1 GW for turns 200 and 250.The Turn 200 plot has additional clip on the right with the FEL energy diffusion function.The RMS length of optical pulse is σ opt = 0.35 at turn 200.The circle and arrow show the trajectory of synchrotron motion for a particle with a s = 10.