Next Article in Journal
The Semiclassical Limit of the Gailitis Formula Applied to Electron Impact Broadening of Spectral Lines of Ionized Atoms
Next Article in Special Issue
Self-Organization in Cold Atoms Mediated by Diffractive Coupling
Previous Article in Journal
Photoionization Cross-Sections of Carbon-Like N+ Near the K-Edge (390–440 eV)
Previous Article in Special Issue
Few-Cycle Infrared Pulse Evolving in FEL Oscillators and Its Application to High-Harmonic Generation for Attosecond Ultraviolet and X-ray Pulses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

A Review of High-Gain Free-Electron Laser Theory

1
Dipartimento di Fisica “Aldo Pontremoli”, Università degli Studi di Milano, Via Celoria 16, 20133 Milano, Italy
2
Centro de Laseres Pulsados (CLPU), Parque Cientifico, 37185 Salamanca, Spain
*
Author to whom correspondence should be addressed.
Submission received: 6 March 2021 / Revised: 7 May 2021 / Accepted: 10 May 2021 / Published: 12 May 2021
(This article belongs to the Special Issue Collective Atomic and Free-Electron Lasing)

Abstract

:
High-gain free-electron lasers, conceived in the 1980s, are nowadays the only bright sources of coherent X-ray radiation available. In this article, we review the theory developed by R. Bonifacio and coworkers, who have been some of the first scientists envisaging its operation as a single-pass amplifier starting from incoherent undulator radiation, in the so called self-amplified spontaneous emission (SASE) regime. We review the FEL theory, discussing how the FEL parameters emerge from it, which are fundamental for describing, designing and understanding all FEL experiments in the high-gain, single-pass operation.

1. Basic Concepts

The free-electron laser is essentially a device that transforms the kinetic energy of a relativistic electron beam (e-beam) into e.m. radiation [1,2,3,4]. The e-beam passing through a transverse periodic magnetic field oscillates in a direction perpendicular to the magnetic field and the propagation axis, and emits radiation confined in a narrow cone along the propagation direction. The periodic magnetic field is provided by the so-called undulator, an insertion device usually realized with two arrays of permanent magnets with alternating polarities (see Figure 1) or with two helical coils with current circulating in opposite directions. The wavelength of the emitted radiation depends on the undulator period, on the strength of the magnetic field and on the electron energy. This means that FELs can be continuously tuned in wavelength, ranging from microwaves (cm) to X-rays (Å); this is one of the main advantages of FELs with respect to atomic and molecular lasers, where the wavelength of the radiation field is fixed by the quantum transition between two atomic or molecular states, and it has, in general, a small tunability. Moreover, another advantage of the FEL with respect to atomic and molecular lasers is that its main processes happen in vacuum, with no thermal dispersion or breakdown effects in the active medium. The FEL is a powerful source of tunable, coherent e.m. radiation [5]. In general, the FEL radiation can be of two different kinds: basic spontaneous emission, i.e., synchrotron radiation, which comes from the direct interaction of the e-beam with the undulator magnetic field with no injected field, while stimulated emission occurs when also a seed radiation field co-propagates with the electron beam. We use the term "spontaneous emission" in analogy with the atomic laser physics (i.e., emission from the spontaneous decay of the atomic excited state) although its origin here is different: it corresponds to the classical incoherent undulator emission from accelerated charges.

1.1. FEL Spontaneous Emission

It is well-known that a free charge cannot radiate because of the energy–momentum conservation; the word “free” in FEL describes the state of the electrons used as active medium: the electrons are not bound as in atomic or molecular lasers. In the FEL, the electrons are not really “free” because they are accelerated by the periodic magnetic field. This accelerated motion generates a radiation field which is particularly intense and directed in a small cone around the instantaneous velocity vector. The FEL radiation is really a synchrotron radiation, i.e., the radiation emitted by an electric charge moving at relativistic speed when a transverse force is applied to it. This radiation power is greater than that due to a parallel force by a factor γ 2 , where
γ = 1 1 β 2
is the Lorentz factor, γ m c 2 the electron energy, m the electron mass unit and β = v / c .
The electrons inside the undulator are periodically deflected by the Lorenz force
F = e v × B w
exerted by the magnetic field B w on electrons with charge e , traveling at speed v . Since the magnetic field has a periodic alternated polarization, the electrons will “wiggle”, i.e., oscillate transversely along the undulator axis. Its main features are the following:
(i)
The intensity is proportional to the electrons’ current, i.e., the radiation is incoherent ( I N e where N e is the number of electrons).
(ii)
The emitted radiation is confined in a narrow cone along the direction of electrons motion (that will be identify with the z-axis) within an angle of order of 1 / γ .
(iii)
It is a narrow-band radiation, with on-axis spectral distribution
d 2 I d Ω d ω = 4 N w 2 γ 2 a w 2 sinc 2 π N w ω ω s ω s
where sinc ( x ) = sin ( x ) / x ,
a w = e λ w B w 2 π m c
is the undulator parameter, λ w is the undulator period and N w is the number of undulator periods. The spectrum is peaked around a spontaneous frequency ω s = 2 π c / λ s where, for on-axis radiation,
λ s = 1 β β λ w
where β = v / c and v is the average longitudinal velocity. The resonant condition (5) can be derived requiring that during the time necessary for an electron to travel an undulator period λ w the electromagnetic wave slips over it by a radiation wavelength λ s ,
λ w v = λ w + λ s c .
The radiation line-width, from Equation (3), is
Δ ω ω s 1 N w
The above result can be easily understood in the (average) longitudinal electron rest frame: here each electron “sees” an N w -periods undulator magnetic field as an N w -periods counter-propagating pseudo-radiation field (known as “Weizsacker–Williams Approximation” [6,7]), with Lorentz contracted wavelength λ w = λ w / γ . Hence, it oscillates N w times, emitting a sinusoidal wave train of length N w λ w at a wavelength λ s = λ w . In other terms, it acts as a “relativistic mirror” where the radiation is back-reflected. From this picture, we obtain the same result of Equations (3)–(7). In fact, by Lorentz-transforming the incident and reflecting wavelengths λ w and λ s back to the laboratory frame, we obtain the relation (5). Moreover, it is well-known that the Fourier transform of a plane-wave truncated after N w oscillations is a sinc-function with line-width Δ ω / ω s = 1 / N w .
From Equation (6), we can write ( k w + k s ) d z / d t ω s = 0 , (where k w = 2 π / λ w and v = d z / d t ) or equivalently d ( k w + k r ) z ω r t / d t = 0 , showing that the resonant relation can be obtained also by imposing the relative phase
θ = ( k w + k s ) z ω s t
of the electron in the undulator and e.m. fields to be constant. As we will show later, the module of electron transverse velocity v = c β is approximately:
β a w γ .
Thus, from 1 / γ 2 = 1 β 2 β 2 , it follows
1 γ 2 = 1 + a w 2 γ 2 .
where γ 2 = 1 / ( 1 β 2 ) . Finally, using the resonant condition (5) in the ultrarelativistic limit γ 1 , we obtain:
λ s = 1 β β λ w λ w 2 γ 2 λ w 1 + a w 2 2 γ 2 .
This relation shows the high tunability of the FEL; in fact, the wavelength λ s can be tuned by varying either the electron energy γ or the undulator magnetic field B w .

1.2. FEL Stimulated Emission

Stimulated emission takes place when a radiation field with wavelength λ λ s co-propagates with the electron beam inside the undulator. From the resonant relation (11) one can define the resonant electron energy
γ r = λ w ( 1 + a w 2 ) 2 λ
If the electron energy and the radiation wavelength preserve the relation (12), then the relative phase (8) between the transverse oscillations of the electron and the radiation remains constant. Depending on the value of this relative phase, one of these processes occurs for each electron:
(a)
the electron gives energy to the field and decelerates, i.e., stimulated emission which provides “gain”,
(b)
the electron takes energy from the field and accelerates, i.e., absorption.
If the first of these two processes dominates, then the injected radiation field is amplified, as in the first amplifier experiment [8]; moreover, if the undulator is long enough or if the process happens in an optical cavity, then the spontaneous emission is amplified, as in the first FEL oscillator experiment [9]. This oversimplified picture of the FEL gain process is at the single-particle level and leads to the Madey’s small-gain regime [10]: let us consider a “long” electron pulse so that the slippage of the radiation over the electrons can be neglected (steady-state regime). Then, the initial electron phases θ are randomly distributed over a radiation wavelength. As a consequence, for a nearly mono-energetic and resonant electron beam, on average half of the electrons will decelerate and half of the electrons will accelerate, with the result that no net gain will occur. The Madey’s small-gain regime occurs when we inject, in a “short” undulator, an electron beam with average energy slightly above resonance, γ 0 > γ r , such that gain (slightly) prevails over absorption [11].

1.3. High-Gain Regime and SASE

More generally, electrons will communicate with each other via the common radiation field. In fact, if the undulator is long enough and the electron current is high enough, then the electrons will start to bunch within an optical wavelength: electrons faster than γ r will decelerate, slower electrons will accelerate, so that the electron energy will be driven toward resonance. This energy modulation, after a transient time called “lethargy”, will become a space modulation, i.e., the electrons will start to longitudinally bunch in pancakes on the scale of the radiation wavelength (micro-bunching), around a phase that produces gain. Since most of electrons will have, at this stage, nearly the same phase, they will emit coherently. The quantity measuring how strongly bunched the electrons are is the bunching parameter [12]:
b 1 N e j = 1 N e e i θ j e i θ .
Given the physical meaning and mathematical definition of the phase θ , it should be clear how b is the measure of the longitudinal modulation of the electron beam on the scale of the radiation wavelength. A bunching equal to zero represents a completely uniform/random distribution of phases, while an ideal bunching of | b | = 1 can only be possible with all electrons perfectly in phase. We shall see the existence of a collective instability for the system, which leads to electron self-bunching (up to a value | b | 0.8 ) and to exponential growth of radiation until saturation, which set a limit on the conversion of kinetic electron energy into radiation energy; this is high-gain steady-state regime [12,13]. In this regime, the peak of the radiated power is proportional to N e 4 / 3 (see Equation (99)), whereas in the incoherent radiation case the power is only proportional to N e .
New effects appear when the propagation effect, due to the different velocity of the electrons and the radiation beam, is taken into account. As a consequence of this effect (known as slippage), the radiation pulse moves ahead with respect to the electron pulse by a radiation wavelength λ r every undulator period λ w , i.e., by l s = N w λ r at the end of a N w period undulator. From analytical and numerical studies, it has been shown [14,15,16,17] that when the slippage length l s is larger than a properly defined cooperation length l c (defined as the slippage in a gain length, see later), an FEL can operate in the superradiant regime, in which the peak of the radiation power scales as N e 2 . Notice that radiation intensities scaling as N e 2 may also arise from coherent synchrotron radiation emitted by electrons which have been pre-bunched by an external source. The radiation field emitted by these pre-bunched electrons sum up coherently to give a N e 2 scaling; this superradiant emission is typical from coherently prepared systems. In FEL the electrons enter the undulator in an unprepared state, so that emission is initially proportional to N e (incoherent emission). Then, the electrons begin to bunch on interacting with the spontaneous radiation and undulator field, which, with the additional slippage effect, leads to a radiated intensity proportional to N e 2 . Roughly speaking, the slippage introduces a dissipative mechanism for the radiation which, as soon it is emitted, it quickly escapes from the interaction volume inhibiting the re-absorption process. This feature is fundamental for the occurring of superradiance [18]. Superradiance in a single-pass high-gain FEL has been observed experimentally since 2007 [19].
The presence of the collective instability (i.e., an exponential growth with a rate depending on N e ) allows an FEL to operate also in the absence of a seed signal (i.e., self-emission), starting from the fluctuations in the initial particle phases; this mode operation is called self-amplified spontaneous emission (SASE). In the SASE-mode, the radiation pulse contains a random superposition of narrow superradiant spikes [20]. In the time domain, the noisy pattern of the bunching along the electron beam leads to a superradiant spiking in the radiation pulse. The initial irregular spiking, seeded by non-uniformities on the scale of the radiation wavelength, cleans up and tends to a more regular pattern, with one spike every cooperation length. The position of the spikes, however, is random and depends strongly on the initial noise pattern of the bunching.
The SASE-FEL mode operation is important at wavelengths not accessible to coherent radiation sources, such as in the XUV and X region. The realization of a single-pass high-gain FEL in the SASE mode has been realized in the last decades and it is presently one (if not the only one) bright and powerful X-ray source in the world. Such sources emit X-ray pulses with a broad spectrum composed by many random superradiant spikes [20] (see the example shown in Figure 2). This drawback can be alleviated by different techniques, for instance by a seed signal injected at the entrance of the undulator, forcing the electrons to emit in the narrow-band spectrum imposed by the seed [21,22]. However, a more fundamental way to "clean" the SASE spectrum has been proposed operating an FEL in a quantum regime [23], where the SASE spectrum reduces to a single narrow line, strongly increasing the longitudinal coherence of a SASE-FEL-based X-ray source [24]. However, the price to pay is more stringent experimental requirements of the electron beam quality, in terms of emittance and energy spread, at the borderline of the present electron beam technology.

1.4. Quantum FEL

Most of the properties of FEL have been analyzed and found remaining in a strictly classical framework. Additionally, if the original proposal of Madey was based on a quantum description of FEL [10], the derived expression for the gain (known as “Madey gain”) is independent on the Planck constant h, i.e., is classical; however, in attempt to extend the range of the FEL operation into the X-ray region, the quantum effects must be taken into account. In fact, at these wavelengths, the electron recoil due to the emission of a photon becomes comparable to the emission line-width, and a quantum mechanical treatment of the electron–photon interaction becomes compulsory. It is well known that the borderline between the classical and quantum world is defined by the Heisenberg Uncertainly Principle (HUP), which can be written in terms of the position and momentum uncertainty, Δ p Δ x , or in terms of the energy and time uncertainty, Δ E Δ t (where = h / 2 π ). When the Heisenberg inequality is strong (i.e., ≫) quantum effects can be neglected and the system behaves classically, conversely when the inequality tends to similarity (i.e., ∼) quantum effects become relevant and a quantum theory is necessary. Assuming that for a relativistic electron with energy E e = m c 2 γ , the interaction time is proportional to the inverse of the emitted photon frequency Δ t ω r 1 and the energy–time uncertainly relation can be written as follows
δ γ λ c 2 π λ r
where λ r = 2 π c / ω r is the emitted photon wavelength, and λ c = h / m c = 0.024 Å is the Compton wavelength. This relation shows that when reducing the emission wavelength, the strong inequality is not satisfied anymore. Therefore, for a complete understanding of the basic FEL process, we need a quantum theory which describes the interaction between the beam electrons and the emitted photons in terms of discrete momentum exchange. Since a quantum theory should tend to a classical description when the inequality in the HUP is strong (i.e., ≫), we should look for a parameter ruling the transition between the classical and the quantum regime. This parameter can be found starting from Equation (14) which can be written as follows
m c δ γ k r 1
where k r = 2 π / λ r . The dimensionless left side of Equation (15) can be chosen as our transition parameter, because it represents the ratio between the maximum classical electron momentum spread m c δ γ and the photon recoil momentum k r . The classical regime occurs when the electron momentum recoil is greater than the photons recoil, i.e., ( m c δ γ ) k r , whereas the quantum effects become important in the opposite. We observe that the energy spread δ γ has a two-fold relevance: (a) if it refers to the initial energy spread, i.e., the energy spread of the electron beam injected in the undulator, the condition set a restriction on the quality of the electron beam; (b) if it refers to the electron energy spread induced by the momentum recoil due to the photon emission, the condition set a restriction on the strength of the coupling between electron and light. Since in the high-gain regime the maximum induced energy spread is δ γ / γ ρ (where ρ is the FEL parameter, see later), then the quantum FEL regime occurs for [24,26]
ρ ¯ m c γ ρ k r 1 .
The consequence that this condition implies for an experimental realization of the quantum FEL regime will be discussed in Section 4.3. However, we outline that this review presents mainly the well consolidated classical FEL theory, without discussing in details the relatively young theory of quantum FEL, which, due to some yet unsolved aspects, deserves a more appropriate presentation in an evolving research context.

2. Classical Model of Equations

From the discussion of the previous section, it results that a proper classical theory describing a high-gain FEL amplifier must be a many-particle theory. Moreover, the electron dynamics should be self-consistently related to the evolution of the e.m. field dynamics. One can set a self-consistent scheme for the FEL dynamics (similarly to that of the Maxwell–Bloch equations in laser physics [27]) which couple the Maxwell equations with the Newton–Lorentz equations for charged particles moving at relativistic speed in an e.m. field; alternatively one can use a Hamiltonian approach, as done in [28,29,30]. We will obtain the electron dynamics equations starting from a full 3D relativistic Hamiltonian of a single electron interacting with a circular polarized e.m. field and counter-propagating laser beam (e.m. undulator) (instead of the usual static undulator). This choice is motivated by recent increase of interest toward the laser undulator. However, it will be easy to adapt the obtained results to the usual static undulator. The evolution equations of the e.m. field are derived from the Maxwell equations in the slowly varying envelope approximation (SVEA, see later). The longitudinal electrostatic interaction between electrons (i.e., space-charge effects) will be also taken into account, but only in the 1D approximation.

3. 3D FEL Model

The circularly polarized laser undulator and the radiation fields are characterized by the transverse vector potential A which can be written in terms of the polarization vector e ^ = ( x ^ + i y ^ ) / 2 as:
a = e m c A = a L + a r .
a L = e ^ 2 a L e i k L ( z + c t ) + c . c
a r = i e ^ 2 a r e i k r ( z c t ) + c . c
a L , r = e m c 2 k L , r E L , r ( x , z )
where E L , r and k L , r = 2 π / λ L , r are the electric field and the wave number of the laser undulator and of the radiation field with frequency ω L , r = c k L , r , respectively.

3.1. 3D Hamiltonian

The time-dependent Hamiltonian equations can be derived from the modified Hamiltonian principle [28]
δ t 1 t 2 ( p x d x / d t + p y d y / d t + p z d z / d t H ) d t = 0 .
Since we are interested in the systems evolution along the z axis, we change the independent variable from t to z, and using H = E (where E is the total energy), we obtain:
δ z 1 z 2 ( p x d x / d z + p y d y / d z E d t / d z + p z ) d z = 0 .
In Equation (21), ( x , p x ) , ( y , p y ) , ( t , E ) appear as the canonical variables with respect to a new Hamiltonian H 1 = p z . Hence, we can write
d x d z = p z p x , d y d z = p z p y , d p x d z = p z x , d p y d z = p z y ,
d t d z = p z E , d E d z = p z t .
Let H be the relativistic Hamiltonian for an electron interacting with e.m. field
H = c p + e A 2 + m 2 c 2 e V = m c 2 γ e V = E .
where E z = d V / d z . The second equation in (23) yields:
d γ d z = 1 m c 2 p z t e z ,
e z = e m c 2 E z = d Φ d z , Φ = e m c 2 V
and the new Hamiltonian is
H ( x , P x , y , P y , c t , γ ; z ) = γ 2 1 P x 2 P y 2 | a | 2 2 P · a + Φ
where P = γ ( d x / d z ) + a and H = P z are respectively the transverse and longitudinal momenta, in m c units. The above Hamiltonian can be simplified with the following approximations:
  • the fast oscillating term 2 P · a can be neglected,
  • the ultra-relativistic limit, γ 1 + | P | 2 + | a | 2 is assumed,
  • the small term | a r | 2 | a L | 2 is neglected.
With 1–3, the Hamiltonian (27) is approximated by
H ( x , P , c t , γ ; z ) = γ + 1 2 γ 1 + | P | 2 + | a L | 2 i a L * · a r e i θ c . c . + Φ
where θ = k r + k L z c k r k L t is electron phase in the laser-undulator and radiation potential. From the Hamiltonian (28), we obtain the equations of motion:
d x d z = H P x = P x γ
d y d z = H P y = P y γ
d P x d z = H x = 1 2 γ x | a L | 2 i a L * a r e i θ c . c .
d P y d z = H y = 1 2 γ y | a L | 2 i a L * a r e i θ c . c .
c d t d z = H γ = 1 + 1 2 γ 2 1 + P x 2 + P y 2 + | a L | 2 i a L * a r e i θ c . c .
d γ d z = H ( c t ) + e z = k r 2 γ a l * a r e i θ + c . c . + e z .
Notice that, using the definition of the total momenta of an electron in e.m. potential normalized to m c , P β γ = ( d x / d z ) γ a , we obtain the transverse velocity of the electron
β = a γ + d x d z
In the paraxial approximation, d x / d z 1 and Equation (35) reduces to the simple expression β = a / γ .

3.2. Maxwell Evolution Equations

The evolution of the e.m. transverse potential a r and of the longitudinal electrostatic field e z = Φ / z can be determined from the rescaled Maxwell equation, expressed in the Coulomb gauge ( · a r = 0 ):
2 + 2 z 2 1 c 2 2 t 2 a r = μ 0 e m c J
2 + 2 z 2 Φ = e ϵ 0 m c ϱ
where J and ϱ are the transverse density current and the charge density, respectively, for a beam of N e electrons:
J = e c j = 1 N β δ ( 3 ) x x j ( t ) , ϱ = e j = 1 N δ ( 3 ) x x j ( t ) .

3.2.1. Vector Potential

The evolution of the transverse e.m. field can be obtained in the following way: substituting the electron current (38) in (36), we get
D + D + 2 a r = e 2 ϵ 0 m c 2 j = 1 N β δ 2 x x j ( t ) δ z z j ( t )
where we have defined D ± = / z ± ( 1 / c ) / t . The radiation beam which propagates in the opposite direction with respect to the electron motion can be neglected because it does not interact resonantly with the electrons. The complex amplitude a r ( x , z ) is assumed to be a slowly varying function of z and t; this means that the fast variation on the radiation wavelength scale does not affect significantly the field enveloped amplitude. This is known as slowly varying envelope approximation (SVEA):
a r ( x , z , t ) z k r | a r ( x , z , t ) |
a r ( x , z , t ) t ω r | a r ( x , z , t ) |
Notice that this approximation is based on the presence, in the FELs dynamics, of two different scales; one of the order of the radiation wavelength, and the other on the order of the interaction length, normally three or four orders larger than the radiation wavelength. More precisely (as we will show in Section 3.7), the second scale length coincides with the “cooperation” length l c ; moreover, the possibility to separate the two different scales can be done more rigorously by using a multiple scaling approach [31]. As a result of (SVEA), the second-order derivatives in longitudinal and temporal coordinates in Equation (39) can be neglected and
D + D + 2 a r 2 e ^ k r e i k r z c t D + + 2 2 i k r a r + c . c .
Substituting Equation (35) for the electron transverse velocity β , and projecting the Equation (39) on the e ^ direction, we obtain
D + + 2 i 2 k r a r = e 2 ϵ 0 m c 2 k r j = 1 N e a L ( x , z ) e i θ i a r ( x , z ) γ j + ( d x / d z ) e i k ( z c t ) × δ x x j ( t ) δ z z j ( t )
We average Equation (42) on the fast scale θ , over a longitudinal dimension l b several wavelengths long l b = s λ r (where s is an integer number). In particular, we integrate both sides of Equation (42) by 1 / l b R χ s λ r 2 , s λ r 2 d z or similarly by 1 / ( k r l b ) R χ s π , s π d θ where χ a , b is equal to one within [ a , b ] and zero otherwise; then the term on the right side of Equation (42) becomes:
1 k r l b R d z χ s π , s π δ ( z z j ( t ) ) e i θ = 1 l b R d θ χ s π , s π δ ( θ θ j ( t ) ) e i θ = 1 l b χ s π , s π e i θ j ( t ) .
Then, defining the average over N e electrons of a generic function as:
f ( θ , x ) = Σ N e j = 1 N e f ( θ j , x ) δ ( x x j ( t ) )
where Σ is the transverse section of the electron beam and the total longitudinal electron density n = χ s π , s π N e / l b , we obtain:
D + + 2 2 i k r a r = k 2 ω p ω r 2 a L e i θ γ i a r γ ( d x / d z ) e i k ( z c t )
where ω p e 2 n e / ϵ 0 m e is the plasma frequency and n e = n / Σ is the total electron density. The first average in Equation (44) shows that the e.m. field evolution is ruled by the sum of all electron phases. As previously discussed, this average term is zero if the electron phases are homogeneously distributed. Conversely, if most of the electrons have the same phase, then the average is different from zero, and it may drive the radiation dynamics. The last two terms on the right side of Equation (44) are usually neglected; the second term is proportional to a r which is usually much smaller compared to a L ; the third term is fast oscillating and can be neglected too. Then, the evolution equation for the self-consistent radiation field can be written as:
z + 1 c t + 2 2 i k r a r = k 2 ω p ω r 2 a L e i θ γ

3.2.2. Space Charge Effects

Here we investigate the longitudinal micro-bunch space charge effects, i.e., the repulsion between neighbor electrons at the scale of the radiation wavelength. The macroscopic charge effects are not considerate because they become important only for very small e-beam energies. A full investigation of the space charge effect would be useful since the local repulsion between electrons could inhibit the establishment of the micro-bunching in the FEL process. However, a complete solution of the Equation (37) is difficult; rather, here we investigate the paraxial approximation of Equation (37), in which only the longitudinal micro-bunching space charge effects are taken into account, neglecting the transverse operator . Similarly to the evolution equation of the e.m. vector potential a r , we assume the SVEA approximation in Equation (37), such that the evolution of the longitudinal electrostatic field in the paraxial approximation can be written as follows:
e z z = e 2 m ϵ 0 c 2 n e ( θ ) N e j = 1 N δ θ θ j
where n e ( θ ) is the local electron density. Expanding n e ( θ ) in a Fourier series with respect to θ , we obtain
e z z = ω p 2 c 2 n = 1 e i n θ e i n θ + c . c .
where we omitted the term n = 0 in the sum, since we neglect the macroscopic space-charge effect. Finally, the longitudinal electrostatic field is:
e z ( θ ) = k r ω p ω r 2 n = 1 e i n θ e i n θ i n + c . c .
Notice that the first harmonic term (i.e., n = 1 ) is proportional to the bunching factor b = e i θ . Clearly, the longitudinal space-charge force is appreciable only if the electron density is large enough. Similarly to the transverse e.m. field, the e z -field is proportional to the sum over the electron phases. This fact shows that the space-charge effects are negligible also if the electron phases are homogeneously distributed. We observe also that if we sum the space-charge force over the electrons, the result is zero; this means that the space charge force is an internal force with respect to the electron system once that the macroscopic effects have been neglected. Then, the main contribution to the space-charge force comes from its first harmonic term, and we write [32]:
e z ( θ ) = 2 k r ω p ω r 2 1 N e j = 1 N e sin ( θ θ j )
Equation (45) and Equations (29)–(34) together with Equation (49) form a self-consistent system of equations for the 3D dynamics of N e electrons interacting with radiation and laser undulator fields. The basic physics of FEL can be understood more clearly in the 1D scheme approximation. For this reason, before discussing the full 3D model, we review the one-dimensional theory and its results.

3.3. 1D FEL Model

The 1D approach is exhaustive for almost everything of the FELs physics. Moreover, at the same time it offers a very clear simple picture of the basic mechanism of the FEL process. In the previous section, a closed set of 3D equations representing the evolution of the whole system has been obtained. The one-dimensional approximation consists of neglecting any dependence on transverse spatial coordinates (Equations (29)–(32)), so that the dynamical variables, i.e., the electron phases θ j , the electron energy γ j and the dimensionless radiation field amplitude a r , depend only on z and t. The electron dynamic Equations (33) and (34) in the 1D approximation become:
c d t j d z = 1 + 1 2 γ j 2 1 + | a L | 2 i a L a r e i θ j c . c .
d γ j d z = a L k r 2 γ j a r e i θ j + c . c . + e z ( θ j ) .
with j = 1 , , N e . The 1D field evolution is obtained neglecting the transverse dependence in the radiation field a r ( z , x , t ) a r ( z , t ) and assuming a uniform magnetic undulator a L ( z , x ) = a w ; then, by integrating both side of Equation (45) by 1 / Σ Σ d x , where Σ is the transverse section of the beam, we obtain
z + 1 c t a r = k 2 ω p ω 2 a w e i θ γ
where the plasma frequency is ω p e 2 n e / ϵ 0 m e . The evolution equation for each electron phase can be obtained from Equation (50)
d θ j d z = ( k r + k w ) k r c d t j d z = k w γ j 2 γ r 2 γ j 2
where γ r = k r ( 1 + a w 2 ) / 2 k w is the resonant energy introduced in the previous section, (Equation (12)). We assume that during the interaction with the e.m. field, the energy of each electron remains close to the resonant energy, γ γ r (this assumption is known as the “Compton limit”), i.e.,
γ j 2 γ r 2 γ j 2 2 γ j γ r γ r 1
With this assumption, we define a new energy-variable η j = γ j γ r γ r . The self-consistent equations for the electron phases and energies and the radiation field are:
d θ j d z = 2 k w η j
d η j d z = a w k 2 γ r 2 a r e i θ j + c . c . i k r γ r ω p ω 2 e i θ e i θ j c . c .
z + 1 c t a r = k a w 2 γ r ω p ω 2 e i θ .

3.4. Universal Scaling

One of the most useful features of the 1D FEL model is the possibility to introduce in Equations (55)–(57) dimensionless variables so that any physical parameters will not appear explicitly in the equations. This allows a general analysis of the FEL scaling laws and a simple interpretation of the basic physics of the process. We start defining the fundamental FEL-parameter [12]
ρ = 1 γ r a w ω p 4 c k w 2 / 3 .
In terms of ρ , the system of coupled evolution Equations (55)–(57) can be set in a dimensionless form by introducing the following variables and parameters.
p ¯ j = η j ρ A r = ω r ω p ρ γ r a r
z ¯ = z l g , t ¯ = c t l g , l g = 1 2 k w ρ = λ w 4 π ρ
where l g is the gain length. Note that in this universal scaling
ρ | A r | 2 = ϵ 0 | E 0 | 2 2 m c 2 γ r n e = P r a d P b e a m
is the ratio between the e.m. and the electron power, i.e., is the FEL efficiency. Using this “universal scaling”, we obtain:
d θ j d z ¯ = p ¯ j
d p ¯ j d z ¯ = A r e i θ j + c . c . + σ e i θ e i θ j c . c .
z ¯ + 1 c t ¯ A r = e i θ
where σ = 4 ρ ( 1 + a w 2 ) / a w 2 is the space-charge parameter. Note that σ is proportional to ρ , so it can be neglected in the Compton limit ρ 1 (this is because η j = ρ p ¯ j 1 ). As we have anticipated before, the 1D FEL equations in the Compton limit assume a dimensionless form, with the advantage of being solvable without specifying any operating parameters. Once solved, the scaling can be reversed to find the real physical quantities needed for a particular experimental set-up. The 1D FEL Equations (62)–(64) form a Maxwell-pendulum model: in fact, writing the complex field as A r = | A r | exp ( i ϕ r ) and neglecting the space-charge term, we obtain
d 2 θ j d z ¯ 2 = 2 | A r | cos ( θ j + ϕ r )
Of course, this equation greatly differs from an ordinary pendulum equation since amplitude and phase of the field are not constant, but their evolution is determined by Equation (64).

3.5. Steady State Regime

The set of Equations (62)–(64) can be written in a simpler way by a transformation of coordinates
z ¯ = z ¯ , z ¯ 1 = z ¯ v t ¯ 1 β
The differential operators of the FEL equations change as follows:
d d z ¯ z ¯ + 1 v t ¯ z ¯
z ¯ + 1 c t ¯ z ¯ + z ¯ 1
where in the first step of the Equation (67) we assumed z ¯ v t ¯ with ( v the longitudinal electron average velocity). Relation (68) shows that the time derivative can be neglected if the difference between the electrons velocity and the speed of light is unimportant ( 1 β ) 0 v c ) or alternatively if the interaction time is small compared with the slippage time. This limit can be easily understood from Equation (68) by normalizing z ¯ to the undulator length L w and t ¯ to the electron pulse duration τ b = l b / v ; the ratio between the time and space derivative coefficients is:
L w ( 1 β ) l b β = N w λ r l b = l s l b
where l s = N w λ r is the slippage length and we have used the resonant condition λ r = λ w ( 1 β ) / β . If l b l s , the time derivative can be neglected (steady-state regime):
d θ j d z ¯ = p ¯ j
d p ¯ j d z ¯ = A r e i θ j + c . c .
d A r d z ¯ = e i θ .

Constants of Motion

The set of Equations (70)–(72) admits two constants of motion. The first is obtained by differentiating the average momentum
d d z ¯ 1 N e j = 1 N e p ¯ j = d d z ¯ | A r | 2 p ¯ + | A r | 2 = C
where C is a constant. The meaning of this conservation law becomes clear if we return to physical quantities, using the relation (59) and (61):
p ¯ + | A r | 2 = 1 N e j = 1 N e Δ γ j ρ γ r + ϵ 0 | E 0 | 2 2 m c 2 γ r n e ρ m c 2 j = 1 N e γ j + ϵ 0 2 V | E 0 | 2 = C
where n e = N e / V . The above relation describes the energy conservation law (the energy of the undulator has been supposed constant) and shows that the energy of the emitted radiation field is provided by the kinetic energy of the electron beam.
The second constant of motion is the total Hamiltonian of the system with 2 N e + 1 variables:
H ( θ j , p ¯ j , Re ( A r ) , Im ( A r ) ) j = 1 , N e = j = 1 N e p ¯ j 2 2 i j = 1 N e A r e i θ j + c . c .
This Hamiltonian can be written in a more compact form, in which the real and imaginary part of the radiation field are included as canonical variables ( θ 0 , p ¯ 0 ). In fact, defining A r = θ 0 + i p ¯ 0 / 2 N e , we obtain:
H ( θ j , p ¯ j ) j = 0 , N e = j = 1 N e p ¯ j 2 2 2 N e j = 1 N e θ 0 sin ( θ j ) + p ¯ 0 cos ( θ j )
from which the equations of motion (70)–(72) follow as
d θ j d z ¯ = H p ¯ j , d p ¯ j d z ¯ = H θ j , for j = 0 , N e

3.6. Linear Analysis

In FELs at the undulator entrance ( z ¯ = 0 ), the electron beam exhibits a narrow longitudinal energy distribution. For simplicity, we assume that each electron has the same initial energy γ j ( z ¯ = 0 ) = γ 0 , so that the dimensionless momentum p ¯ at time z ¯ = 0 is given by
p ¯ ( 0 ) = γ 0 γ r ρ γ r δ
We define this value as the detuning parameter  δ . It is particularly useful to redefine our variables so that the initial condition for p ¯ is zero:
p ¯ j = p ¯ j δ θ j = θ j δ z ¯ A r = A r e i δ z ¯
In this way, the detuning parameter appears explicitly in the equations which (dropping the primes) read:
d θ j d z ¯ = p ¯ j
d p ¯ j d z ¯ = A r e i θ j + c . c .
d A r d z ¯ = e i θ + i δ A r
This set of equations can be linearized in terms of three collective variables [12]
A = A r , field amplitude
B = e i θ , bunching
P = p ¯ e i θ , momentum bunching
Neglecting the second-order terms e i 2 θ and p ¯ 2 e i θ , we obtain a closed set of linear equations
d B d z ¯ = i P
d P d z ¯ = A
d A d z ¯ = B + i δ A
which can be reduced to
d 3 A d z ¯ 3 i δ d 2 A d z ¯ 2 i A = 0 .
Assuming a solution of the form A e i λ z ¯ , one obtains the following dispersion relation which rules the stability of the system:
λ 3 δ λ 2 + 1 = 0 .
When Equation (90) has three real roots, then the system is stable, but if it has one real root and two complex-conjugate ones, then one of the latter will cause an exponential growth of the field until non-linear effects come into play. If instead of choosing an ideal cold beam, we assume an initial energy distribution f ( p 0 ) with a finite width, then (90) generalize into
λ δ + + f ( p 0 ) ( λ + p 0 ) 2 d p 0 = 0 .
The above integral can be analytically solved, for instance in the case of a rectangular or Lorenz distribution [16]. For the rectangular case with half-width δ γ , we obtain:
( λ δ ) ( λ 2 μ 2 ) + 1 = 0
where μ δ γ / ρ γ r is the energy spread parameter.
As the exponential behavior of A r ( z ¯ ) is determined by the imaginary part of the complex root of Equation (92), it is relevant to plot it as a function of the detuning parameter δ for different values of the energy spread μ . Figure 3 suggests some immediate considerations:
  • given a spread μ , the optimal gain occurs for the specific detuning shift;
  • energy spread ( μ > 0 ) lowers the growth rate, and shift the resonance to δ = μ ;
  • the width of the gain curve shrinks as 1 / μ .
Coming back to physical variables, this means that in order to preserve the exponential gain, for μ = 0 , the e-beam must satisfy:
δ γ γ r ρ
The solution for the cold beam case ( μ = 0 ) at resonance ( δ = 0 ) is:
λ 3 = 1 λ 1 = 1 , λ 2 = 1 + 3 i 2 , λ 3 = 1 3 i 2
and the scaled field is
A ( z ¯ ) = 1 3 k = 1 3 A ( 0 ) i B ( 0 ) λ k + i P ( 0 ) λ k 2 e i λ k z ¯
For z ¯ 1 (exponential growth regime) the growing mode λ 3 dominates over the oscillatory λ 1 and decaying λ 2 modes, so that
A ( z ¯ ) 1 3 A ( 0 ) i B ( 0 ) λ 3 + i P ( 0 ) λ 3 2 e i λ 3 z ¯
The first term in the bracket corresponds to the amplification of an external input signal, while the second and the third terms correspond to some initial bunching. As it results from the linear analysis, an exponential instability of the emitted radiation field takes place in the FEL, until non-linear effects saturate this growth. The linear solution of the exponential growing mode fits well the numerical solution of the full non-linear system (70)–(72).
In Figure 4, the radiated power is plotted, as obtained from the numerical integration of the 1D non-linear equations. The figure shows that after an initial lethargy, the field power grows exponentially, reaching a saturation after a certain number of gain lengths (depending on the initial conditions). After saturation, the amplification process is replaced by an oscillatory energy exchange between the electrons and the radiation field. In the case of a seed signal intensity | A 0 | 2 , the asymptotic linear solution is given by the
| A | 2 1 9 | A 0 | 2 exp 3 z l g ,
so that the l g as defined in Equation (60) corresponds effectively to the FEL gain length. Since saturation occurs at | A | 2 1 , from the Equation (97) we can estimate the saturation length as
z s a t l g 3 ln 9 | A 0 | 2 λ w ρ
Furthermore, the result | A | 2 1 means that the saturation value is independent of the initial conditions. Since we know that | A | 2 is proportional to | E 0 | 2 / ρ n e , and ρ n 1 / 3 then
| A | 2 | E 0 | 2 n e 4 / 3 | E 0 | 2 n e 4 / 3
i.e., the intensity is proportional to n e 4 / 3 instead of n e : this implies the existence of a collective behavior in the electron beam. The independence of the saturated field amplitude on its initial value gives us another important information, related again to the importance of the ρ parameter. Since in Equation (61) we have defined the efficiency η of the FEL as η = ρ | A | 2 , then the measure of FEL efficiency is given by ρ . In this section we have shown that in the context of the 1D theory, ρ is the only parameter which determines the constraints on the FEL:
Δ γ / γ ρ Energy spread η = P r a d / P b e a m ρ Efficiency l g 1 / ρ Gain length

3.7. Superradiant Regime

The steady-state regime that we just described is based on the assumption that the slippage is negligible ( l b l s , see Equation (69)): the undulator is not long enough to appreciate the difference in velocity between the electrons and the radiation, so that all sections of the electron beam evolve almost identically and the peak power of radiation scale as n e 4 / 3 (see Equation (99)). When slippage is taken into account, the FEL can operate in a different regime of cooperative emission, the superradiant regime [14,15,16,33,34], where the peak power scale as n e 2 . The slippage modifies substantially the interaction process between the radiation and the electrons: in fact, the radiation propagates with respect to the electrons, interacting with different slices of the electron beam, such that there is a region, near the trailing edge of the electron pulse and of length l s , where the electrons emit radiation without being affected by the radiation produced by the other electrons behind them. Let us introduce a new characteristic length which is useful for the analysis of the propagation effect, the cooperation length,
l c = λ r 4 π ρ .
Its meaning can be better understood using the resonant relation (5), which allows to write it in terms of the relativistic parameter β :
l c = l g 1 β β
so that it can be interpreted as the slippage in a gain length l g . We define the electron pulse to be long or short with respect to the cooperation length. In general, if the electron pulse is long enough ( l b l c ), the superradiant and the steady state instability will be present together: superradiance occurs in the region near the trailing edge of the electron pulse (slippage region) while the steady-state emission occurs in the remaining part of the beam: this is named “strong superradiance” because the peak power is greater than that in the steady-state regime.
Figure 5 shows the results of a numerical simulation integrating the 1D FEL equations with slippage:
θ j z ¯ = p ¯ j
p ¯ j z ¯ = A r e i θ j + c . c .
z ¯ + z 1 A r = e i θ + i δ A r
where z ¯ = z / l g and z 1 = ( z c β t ) / l c . In Figure 5, the electron beam length is l b = 30 l c , with a flat current profile in the interval 0 < z 1 < 30 and the detuning δ = 0 . We observe a steady-state intensity in z ¯ < z 1 < 30 and a superradiant peak growing near the trailing edge ( z 1 = 0 ), in the slippage region 0 < z 1 < z ¯ (see Figure 5b). The intensity observed in the region 30 < z 1 < 30 + z ¯ is the trace of the radiation propagating forward, in front to the electron leading edge z 1 = 30 . The peak of the superradiant pulse in the slippage region strongly exceeds the maximum value of the steady-state intensity, | A r | 1.4 .
In the short-bunch regime ( l b l c ) the radiation emitted by electrons escapes from the electron pulse in a length shorter than the gain length, so that the steady-state saturation never occurs: this regime is called “weak superradiance” because the peak power is lower than that in the steady-state regime.
Figure 6 shows the results of a simulation with l b = 0.1 l c , δ = 0 and different values of z ¯ . In this case, almost all the radiation intensity is outside of the electron pulse, propagated from the trailing edge z 1 = 0.1 , and no steady-state emission occurs. Weak superradiance in high-gain single-pass FELs has been observed experimentally in the past years by several groups [19,35].
It is possible also to observe pure superradiance by tuning the system out of resonance: in fact, while the steady-state regime needs resonance (i.e., δ = 0 ) to produce exponential gain, superradiance weakly depends on the detuning, since it has a much larger gain bandwidth. Therefore, when the system is detuned in such a way to prevent steady-state radiation, the superradiant instability travels forward over unperturbed electrons, extracting energy from them with an even greater efficiency than in the steady-state regime.
This case is shown in Figure 6, with the same parameters as Figure 5 but with a detuning δ = 2 , such that the steady-state emission is inhibited. In this case the superradiant pulse, born near the trailing edge of the electron beam, travels unperturbed over the electrons toward the leading edge, with a peak intensity growing as z ¯ 2 . It has been demonstrated [15,34,36] that the superradiant pulse, as clearly observed in Figure 6 and Figure 7, is described by a self-similar solution of Equations (103)–(105). It is obtained by setting θ j ( z ¯ , z 1 ) = θ 1 j ( y ) , p ¯ j ( z ¯ , z 1 ) = p 1 j ( y ) / z 1 and A r ( z ¯ , z 1 ) = z 1 A 1 ( y ) , where y = z 1 ( z ¯ z 1 ) and θ 1 j ( y ) , p 1 j ( y ) and A 1 ( y ) are the solutions of the following ordinary differential equations:
d θ 1 j d y = p 1 j
d p 1 j d y = A 1 e i θ 1 j + c . c .
y 2 d A 1 d y + A 1 = e i θ 1
Figure 8 shows | A 1 | 2 as it results from the numerical solution of Equations (106)–(108).
From the definition of the dimensionless radiation amplitude, the power of the pulse is P r = ρ P b e a m z 1 2 | A 1 ( y ) | 2 ρ 3 P b e a m N e 2 , i.e., superradiant, and its width decreases as 1 / N e . In the weak superradiance, z 1 = l b / l c and y l b / l c z ¯ , so that | A r | 2 ( l b / l c ) 2 | A 1 ( l b / l c z ¯ ) | 2 (see Figure 6). The peak power is proportional to l b 2 and the pulse width is proportional to 1 / l b . In the strong superradiance, the peak pulse grows as z ¯ 2 and shrinks as 1 / z ¯ when it propagates over the electron beam (see Figure 7).

3.8. SASE Operation

The self-amplified spontaneous emission (SASE) operation for an FEL is made up of three basic ingredients [20]:
  • high gain instability;
  • propagation effects, i.e., “slippage”;
  • start-up from noise.
The first experimental observation of the high-gain regime, also starting from noise, was carried out in the microwave range using a wave guide in the Livermore experiment [37]. Presently, short wavelength FEL, which amplify incoherent shot-noise via SASE are of greater interest worldwide as a source of ultra-bright coherent X-ray radiation. SASE FEL in the X-ray region has been demonstrated firstly in LCLS project facility [38] at Standford, CA, USA, in the European X-ray FEL [39] at Hamburg, Germany and in Riken, Japan [40]. Actually, several X-ray FELs operating in the SASE regime have been and will be realized over the world [5]. The most important features of the SASE FEL operation are determined by the following characteristic lengths [20]:
l b bunch length l c cooperation length l b gain length
An important parameter determining the evolution of the system is given by the ratio between l b and l c (the number of cooperation lengths in a bunch length). If the ratio l b / l c > 2 π (long bunch case), then the radiation pulse contains many spikes, each one having a maximum duration corresponding to about 2 π l c with a width of order of 1 / l b and large intensity fluctuations (see Figure 9).
On the contrary, if the ratio l b / l c 2 π (short bunch case), only a single radiation pulse is present in this case, with no inner spikes. In both cases, superradiance occurs for sufficiently long undulators. The final result is an almost chaotic temporal pulse structure, with a broad spectral width and with a number of spikes of the order of the number of cooperation lengths in the electron bunch (see Figure 10 for the power spectrum of the intensity profile shown in Figure 9 vs. ω ¯ = ( ω ω s ) / 2 ρ ω s , where ω s = c k w γ 2 / ( 1 + a w 2 ) is the spontaneous frequency).
Finally, in order to model correctly this spiking behavior, we must take into account not only the slippage between radiation electron pulse, but also the finite bunch length; the electron bunch behaves as a l b / l c statistical independent zones, each of them giving rise to a superradiant spike which grows and narrows extracting energy from electrons within a cooperation length [20,41].
In the “so-called” quantum SASE regime [24], a completely different behavior occurs, the “classical” random spiking behavior almost disappears and a strong narrowing of the spectrum occurs. This new phenomenon is called “quantum purification”, also if its experimental evidence is yet to come.

4. From 1D to 3D

The description of an N e electron beam can be made introducing, in the transverse phase space ( x , η ) (where η = d x d z and z is the longitudinal coordinate along the beam direction), the second-order moments, as for instance the rms beam size σ e 2 = | x | 2 , the rms beam angular divergence σ e 2 = | η | 2 and the mixed terms x · η . The rms beam size (for example in x direction) evolves in the free space as:
σ e ( z ) = ϵ x β x * + z 2 β x *
where β x * = x 2 / ϵ x is one of the Twiss parameters and ϵ x = x 2 η x 2 x η x 2 is the rms emittance along x, which is conserved in the free space evolution and for a linear transport system; z = 0 corresponds to the point where x η x = 0 and the beam size is minimum (“beam waist”). The radiation beam can be described in a similar way: in fact, using the free space diffraction formula, it is possible to show that the size of a Gaussian beam evolves as:
σ r ( z ) = λ r 4 π Z r + z 2 Z r
where Z r = 4 π σ r 2 / λ r is the Rayleigh range. The evolution equations of the electron and radiation beam sizes are formally equivalent if we identify ε r a d = λ r / 4 π as the rms emittance of the Gaussian radiation beam and β * as the Rayleigh range of the electron beam. Furthermore, the free space evolution of the counter-propagating laser beam size σ L ( z ) follows the same relation of Equation (111) where instead of Z r we substitute Z L = 4 π σ L 2 / λ L where λ L is the laser beam wavelength.

4.1. Transverse Effects

From this picture, we can set some reasonable condition that must be satisfied in a 3D geometry:
  • The matching between electron and radiation beam requires that the beam waist and the Rayleigh range of each other must be comparable:
    σ e σ r
    β * Z r k r ϵ r 1 2 σ e σ r 2
  • The electron beam should be contained in the laser beam and the electron beam should not diverge appreciably in a Rayleigh range Z L
    σ e σ L
    β * Z L k L ϵ r 1 2 σ e σ L 2
The condition (113) is known as the “Pellegrini criterium” [42] and it can be relaxed in particular conditions, depending on the ratio between the electron and the radiation beam size (112). Conditions (114) and (115) are more stringent and concern the matching between the e-beam and the laser undulator profile. If conditions (114) and (115) are violated, inhibition of the gain process occurs. If we use a magnetic undulator instead of the laser undulator, then the homogeneous condition (115) can be written as k w ϵ r < 1 , which is usually satisfied. In fact, for an undulator period of the order of λ w = 1   cm and emittance ε r 10 7 m - rad , k w ϵ r 10 4 . An FEL experiment operating with a laser undulator does not require a focusing system for the electron beam, since the interaction between electrons and laser beam occurs in free space near the beam waist zone, i.e., the interaction length L i n t is of the order of few laser Rayleigh ranges ( L i n t Z L ). Hence, the gain process occurs if the condition (115) is satisfied. Furthermore, if a magneto-static undulator instead of a laser undulator is used, then the interaction length is usually more than one meter, and the electron beam must be conveniently focused to keep a nearly constant beam size. A linear magneto-static undulator provides a “natural” focusing effect in one transverse direction and needs external focusing (by quadrupole magnets) on the other transverse direction. This focusing force induces an oscillation of the electron motion with a wavelength much longer than the undulator period λ w . For an electron beam matched in two planes, the beam size (averaged over the undulator period λ w ) remains constant along the undulator, while individual electrons perform a periodic transverse motion, called “betatron motion”.
A potential deleterious effect for an FEL is the apparent energy spread induced by the betatron motion, i.e., by the electron transverse velocities. The relation between energy spread and transverse velocity can be made more transparent taking into account the transverse dimension in the resonant condition; it straightforward to demonstrate that an observer far away in the η direction from the axis sees a slightly different period, i.e.,
λ r = λ w 2 γ r 2 1 + a w 2 + η 2 γ 2
where λ w is replaced by λ L / 2 for a laser undulator. The above relation agrees with the resonance condition (11) when the observer angle tends to zero. From the transverse resonant relation (116), using the relation Δ γ / γ Δ λ / 2 λ , we can argue that the energy spread normalized to the FEL parameter ρ has the following contributions:
  • different longitudinal momentum distribution, (see Equation (93))
    Δ γ ρ γ r 1 D 1
  • off-axis variation of the undulator parameter
    Δ γ ρ γ r a w 1 2 ρ Δ a w 2 1 + a w 2
  • angular divergence of the beam
    Δ γ ρ γ r 1 2 ρ η 2 γ r 2 1 + a w 2
These effects are called “non homogeneous effects” and give rise to a broadening of the resonant condition.

4.2. Full 3D Model

A most useful scheme for modeling a beam of N e electrons interacting with e.m. field is the Maxwell–Vlasov scheme. Following this scheme, we introduce the transverse phase space variables ( x , η ) and the longitudinal phase space variables ( θ , p ¯ ). The transverse phase-space distribution, for a given phase and longitudinal momentum, can be described by a distribution function F ( x , η ) with a non-negative value equal to the number of electrons per unit area at the transverse phase-space point ( x , η ) . A Gaussian distribution function is normally adopted for its simplicity. Hence, the number of electrons within a d x d η at the beam waist (such that x · η = 0 ) is
F ( x , η ) d x d η = N e 2 π ϵ r exp x 2 2 σ 2 η 2 2 σ 2 d x d η
The evolution equation of the electron beam distribution function can be obtained using the Liouville theorem f / z ¯ = H , f where H , f is the Poisson bracket. In the previous sections, we have obtained a three dimensional Hamiltonian with its associated equations of motion, Equations (29)–(34), and an evolution equation for the radiated e.m. field (45). Now, we introduce the “3D Universal Scaling”:
θ = k r + k L z c ( k r k L ) t , z ¯ = z l g , x ¯ = x σ , p ¯ = γ γ r ρ γ r , p ¯ = σ ϵ r η , η = d x d z , a L = a 0 g ( z ¯ , x )
where a 0 = a w , g ( z ¯ , x ¯ ) is the transverse laser profile, σ is the rms electron beam radius at the beam waist and ϵ r = σ σ is the rms beam emittance. With this universal scaling the equations of motion (29)–(34) can be written:
d x ¯ d z ¯ = b ( 1 ρ p ¯ ) p ¯ d p ¯ d z ¯ = 2 ρ X ( 1 ρ p ¯ ) x ξ 2 ρ 2 | g | 2 i g * A e i θ c . c . d θ d z ¯ = p ¯ + ξ 2 ρ ( 1 | g | 2 ) b X 4 p ¯ 2 d p ¯ d z ¯ = ( 1 ρ p ¯ ) g * A e i θ + c . c . + E .
Whereas the 1D theory is ruled only by the FEL parameter ρ , instead, in the 3D theory, we must introduce new parameters which characterize the transversal dynamics:
X = 2 k r ϵ r , b = l g β * = l g ϵ r σ 2 ξ = a 0 2 1 + a 0 2 .
In the Compton limit, the terms proportional to ρ (as for instance in ( 1 ρ p ¯ ) 1 ) can be neglected, since the value of ρ ranges from about 10 5 to 10 3 . With this simplification, Equation (122) becomes:
d x ¯ d z ¯ = b p ¯ ¯ d p ¯ ¯ d z ¯ = ξ ρ X x | g | 2 d θ d z ¯ = p ¯ + ξ 2 ρ ( 1 | g | 2 ) b X 4 p ¯ 2 d p ¯ d z ¯ = g * A r e i θ + c . c . + E .
They can be obtained by the following Hamiltonian
H ¯ = p ¯ 2 2 + b 2 p ¯ 2 + p ¯ ξ 2 ρ ( 1 | g | 2 ) b X 4 p ¯ 2 + ξ ρ X | g | 2 i g * A r e i θ c . c . + Φ ¯ .
where Φ ¯ = γ r / ( ρ 2 ( 1 + a 0 2 ) ) Φ is such that E = Φ ¯ / θ . Now, using the Liouville theorem, the Hamiltonian (125), the Maxwell Equation (45) and 3D universal scaling (121), we obtain the following Maxwell–Vlasov system of equations for an e-beam interacting with a laser undulator and a radiation field:
f z ¯ + b p ¯ · x ¯ f + p ¯ + ξ 2 ρ ( 1 | g | 2 ) b X 4 p ¯ 2 f θ + E g * A e i θ + c . c . f p ¯ ξ ρ X x ¯ | g | 2 p ¯ f = 0 z ¯ + z ¯ 1 i a x ¯ 2 A r = g π + π d θ R d p ¯ R 2 d 2 p ¯ e i θ f + i δ A r
Here a = l g / Z r = b / X and Z r = 4 π σ 2 / λ r are respectively the diffraction parameter and the Rayleigh range of the emitted radiation, with a transverse radius equal to the electron beam radius. The distribution function f and the radiation field A r have the following dependence f = f ( θ , p ¯ , x ¯ , p ¯ , z ¯ , z ¯ 1 ) , A r = A r ( θ , x ¯ , z ¯ , z ¯ 1 ) , and g = g ( x ¯ , z ¯ ) is the laser undulator profile. The term b p ¯ x ¯ corresponds, with the unscaled variables, to η x and describes the transverse drift of the beam, responsible for instance of the beam size increase away from the beam waist; the comparison of the two diffraction terms a b (which corresponds to Equation (113)) rules the matching condition between the electron and radiation beam. In fact, assuming the same initial spot size, if b < a the electron beam is contained into the radiation beam, which corresponds to the “Pellegrini condition” X < 1 [42], otherwise the radiation is confined into the electron beam size. The phase term of Equation (125) contains three terms which, with the unscaled variables, can be written as:
p ¯ = γ γ r ρ γ r = Δ γ ρ γ r 1 D 1 ξ 2 ρ ( 1 | g | 2 ) 1 2 ρ Δ a w 2 1 + a w 2 Δ γ ρ γ r a w b X 4 p ¯ 2 1 2 ρ ϑ 2 γ r 2 1 + a w 2 Δ γ ρ γ r
and are responsible of the “non homogeneous effects” (117)–(119). Finally, the last term d p / d z = ξ / ( 2 ρ X ) x | g | 2 corresponds with unscaled variables to a focusing force due to the laser undulator profile η ¯ ˙ = a w 2 / ( 2 γ 0 ) 2 x | g | 2 . Notice that p ¯ · x ¯ = p ¯ x / x ¯ + p ¯ y / y ¯ and the average in the field evolution equation is replaced by the average over an ensemble, O = O f ( ξ ) d ξ · with f ( ξ ) d ξ = 1 . This closed set of equations can be solved numerically choosing a Gaussian initial condition at f ( z ¯ = 0 ) = f 0 :
f 0 ( θ , p ¯ , x ¯ , p ¯ ) exp x ¯ + b z ¯ 0 p ¯ 2 2 p ¯ 2 2
where z ¯ 0 is the waist position. Typically, a Gaussian phase-space profile corresponds to a thermal distribution.

Optical Guiding

One of the most peculiar results of the high-gain FEL theory is the presence of optical guiding, i.e., the existence of exponentially growing modes which have a profile independent of the longitudinal coordinate [43,44]. This is particularly important in the case of long undulators, where the radiation tends to affected by diffraction. Neglecting emittance, energy spread, betatron oscillation and slippage, the linear FEL equations taking into account diffraction in the paraxial approximations are obtained by combining Equations (86)–(88) and (125) to give
2 B z ¯ 2 = i A z ¯ i a r ¯ 2 A = u ( r ¯ ) B + i δ A
where A and B depend on the transverse position r ¯ and u ( r ¯ ) is the normalized transverse electron distribution. Looking for solutions which depend exponentially on z ¯ as exp ( i λ z ¯ ) and are azimuthally symmetric, Equation (127) becomes
( λ δ ) A a 1 r ¯ d d r ¯ r ¯ d A d r ¯ = u ( r ¯ ) λ 2 A
Outside the electron distribution ( u = 0 ), the solution (normalized to unity) is
A = H 0 ( ϕ r ¯ )
where H 0 is the Hankel function of zero order of the first kind, and ϕ 2 = ( δ λ ) / a . For large r ¯ , the asymptotic limit is
A 2 π ϕ 1 / 2 exp ( i ϕ r ¯ i π / 4 ) .
In order that A 0 for r ¯ , we require that Im ϕ > 0 . Since we want an exponentially growing mode, then I m λ < 0 and Im ϕ 2 > 0 , so that it must be also Re ϕ > 0 . Considering an uniform sharp-edge distribution, with u = 1 / π for r ¯ < 1 and u = 0 for r ¯ > 0 , the mode amplitude satisfying Equation (128) inside the electron beam and normalized to unity is
A = J 0 ( χ r ¯ )
where J 0 is the Bessel function of zero order of the first kind and
χ 2 = 1 a δ λ 1 λ 2 = ϕ 2 1 a 3 ( ϕ 2 δ ^ ) 2
where δ ^ = δ / a . When a is small (i.e., Z r l g ), diffraction is unimportant and the gain is given by the 1D theory, where λ is solution of the cubic ( λ δ ) λ 2 + 1 = 0 . When a is large, the gain is less than that expected from 1D theory. The values of λ are determined imposing the continuity of the A and its derivative ar the boundary r ¯ = 1 , finding that this implies
ϕ H 1 ( ϕ ) J 0 ( χ ) = χ J 1 ( χ ) H 0 ( ϕ ) ,
looking for values of ϕ and χ which are simultaneous solutions of Equations (132) and (133). A detailed numerical study has been done in [43], showing that the condition a = 1 (i.e., Z r = l g ) measures the distance over which diffraction establishes transverse coherence in the high-gain FEL. If Z r l g , one-dimensional theory gives the correct gain, but there is not full transverse coherence. If Z r < l g , the fundamental transverse mode has gain substantially larger than that of the other modes, the laser mode area is much larger than the electron beam area and a full transverse coherence is established.

4.3. Quantum Regime of FEL

We conclude with some remarks about the experimental constraints of the quantum FEL regime, as described in [24,45]. It has been proposed [46,47] that one can use, instead of a magneto-static undulator, a high-power laser beam as undulator for the X-ray emission region. In fact, since λ r λ w / 2 γ 2 , the wavelength of the emitted radiation can be reduced or increasing the electron energy or, alternatively, reducing the undulator period. An infrared laser, coming from high-power N d or C O 2 lasers, with a wavelength λ L of 1 or 10 m could yield X-ray FEL radiation with electrons energy about a hundred MeV.
λ r λ L 4 γ 2 0.1   nm
Now, with a simple argument, we show why a laser undulator is advisable for X-ray FEL emission in the classical scheme, while its is compulsory in the quantum regime. As we have discussed, the transition between quantum and classical regime is ruled by the quantum FEL (QFEL) parameter [24]
ρ ¯ = ρ m c γ k r = ρ γ λ r λ c
where λ c = h / m c is the Compton wavelength and h is the Planck constant. The QFEL-parameter is proportional to the ratio between the electron energy and the single-photon energy and the quantum regime occurs when ρ ¯ 1 . Using the resonant relation and the definition of ρ ¯ in Equation (135), the condition for the quantum regime can be written in terms of ρ :
ρ 2 λ c λ r λ w ( 1 + a w 2 )
We have seen that to reach the high-gain regime, a number of period of the order of N w 1 ρ is required, so that the undulator length must be:
L w = N w λ w λ w ρ λ r λ w 3 ( 1 + a w 2 ) λ c
For an X-ray FEL experiment with λ r 1 Å, using a magnetic undulator with a period λ w 1   cm and E = 3.5 Gev, we obtain from Equation (137) L w 3 km, which is possible but impracticable. If we substitute the magneto-static undulator with a counter-propagating laser beam, we must substitute in the resonant relation and otherwise λ w λ L / 2 , and for a typical laser with wavelength of 1 μ m and an electron energy E = 2.5 Mev, we obtain L w 2 mm. However, the use of a laser undulator would require an exceptional stability in intensity and frequency in the interaction region, which makes it still very challenging [48].

5. Conclusions

We reviewed the theory of high-gain FEL, following the approach developed by R. Bonifacio and coworkers in the 1980s, which underpins the major advances in understanding that led to the next generation of X-ray FELs. The review started with a brief introduction to the basic concepts and quickly developing the equations for “spontaneous” FEL radiation, followed by a description of stimulated FEL emission, which gives an intuitive understanding of the role of bunching, which is also necessary to understand the high-gain regime. Finally, the introductory section introduces the basic concepts of the quantum FEL where electron recoil resulting of photon emission dominates the resonant interaction. The theoretical model, based on the many-particle evolution of the electrons moving in the undulator and self-consistently coupled to the emitted radiation field, has the advantage of showing the emergence of few fundamental parameters describing the main features of the FEL devices. In particular, the 1D dynamics of a Compton (i.e., high-energy) classical FEL is all described by the fundamental FEL parameter ρ , measuring efficiency, gain length l g = λ w / 4 π ρ (where λ w is the undulator period), energy spread and gain bandwidth in the high-gain regime.
A 3D model of the FEL is developed based on a 3D Hamiltonian and Maxwell’s equations. This is expanded to include the role of space-charge effects, which become important for FELs operating in the Raman regime. To gain an intuitive understanding of FEL theory, we discussed in detail the 1D model, using suitably scaled parameters to develop a “universal” set of equations. Then we examined the steady-state regime, where slippage is neglected, and we reviewed the results of the linear analysis and a general dispersion relationship, showing the exponential collective instability. Including the slippage effect, we discussed the FEL superradiance, which is characterized by the cooperation length, l c = ( λ r / λ w ) l g . Superradiance in FELs can be classified as weak superradiance for short electron bunches ( l b l c ) and strong superradiance for long electron bunches ( l b l c ). At the core of the superradiant emission is the existence of a self-similar solution, i.e., a solution which maintains the same profile growing and shrinking when it propagates along the electron bunch. From the self-similar solution, the main scaling laws with the electron beam length, current and undulator length follow naturally. Superradiance is also the the core of the self-amplified spontaneous emission (SASE) operation, very important for the development of the X-ray sources. We described the main features of SASE in the 1D limit, showing some typical temporal profiles and spectra of the radiation intensity. Next, we described three-dimensional effects, as optical guiding and transverse electron motion due to emittance. These effects are characterized by two other dimensionless parameters, a = l g / Z r and b = l b / β * for radiation and electron diffraction, where Z r is the Rayleigh range and β * the betatron length. Finally, we presented the main concept of the quantum regime of FEL, which is characterized by the quantum FEL parameter ρ ¯ = ρ ( m c γ / k r ) , that can be interpreted as the average number of photons emitted by each electron. The quantum regime of FEL occurs when ρ ¯ < 1 . The major constraints for the experimental realization of a quantum FEL have been discussed.
In this review, we preferred to present only the classical theory of FEL, since it is well established and confirmed by many experiments and is at the base of several numerical codes developed for assisting the experiments [49,50,51,52]. About the quantum FEL, it has not been experimentally proved yet and its studies are still under progress [53,54], leaving the subject still in an advanced research context, with the hope that in the near future, due to the fast progresses in particle-beam and laser technology, it may become a reality. However, since the quantum FEL stems from an original idea of Rodolfo Bonifacio, we have presented the main basic concepts of the quantum FEL, but without entering into a detailed description of its quantum theory.
Since FEL has been studied since more than fifty years, many other authors have contributed to the advance of the FEL physics, both theoretically and experimentally, and it would be a hard task to make reference of all of them. About that, several excellent reviews and books on FELs exist [1,2,3,4]. The main aim of this paper was to review the main contributions made to the FEL physics by Rodolfo Bonifacio, to which this work is dedicated.

Author Contributions

All authors conceived the study. N.P. directed the research. All authors contributed to writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Freund, H.P.; Antonsen, T.M. Principles of Free-Electron Lasers; Springer: Berlin, Germany, 1992. [Google Scholar]
  2. Saldin, E.; Schneidmiller, E.; Yurkov, M.V. The Physics of Free Electron Lasers; Springer: Berlin, Germany, 1999. [Google Scholar]
  3. Schmüser, P.; Dohlus, M.; Rossbach, J.; Behrens, C. Free-electron lasers in the ultraviolet and X-ray regime. In Springer Tracts in Modern Physics; Springer: Berlin, Germany, 2014; Volume 258. [Google Scholar]
  4. Kim, K.J.; Huang, Z.; Lindberg, R. Synchrotron Radiation and Free-Electron Lasers; Cambridge University Press: Cambridge, UK, 2017. [Google Scholar]
  5. Ueda, K. X-Ray Free-Electron Laser. Appl. Sci. 2018, 8, 879. [Google Scholar] [CrossRef] [Green Version]
  6. Weizsäcker, C.V. Ausstrahlung bei Stößen sehr schneller Elektronen. Z. Phys. 1934, 88, 612–625. [Google Scholar] [CrossRef]
  7. Williams, E.J. Correlation of Certain Collision Problems with Radiation Theory; Levin & Munksgaard: Aboulevard, Denmark, 1935. [Google Scholar]
  8. Elias, L.; Fairbank, W.; Madey, J.M.J.; Schwettman, H.A.; Smith, T.I. Observation of stimulated emission of radiation by relativistic electrons in a spatially periodic transverse magnetic field. Phys. Rev. Lett. 1976, 36, 717. [Google Scholar] [CrossRef] [Green Version]
  9. Deacon, D.A.; Elias, L.; Madey, J.M.; Ramian, G.; Schwettman, H.; Smith, T.I. First operation of a free-electron laser. Phys. Rev. Lett. 1977, 38, 892. [Google Scholar] [CrossRef]
  10. Madey, J.M. Stimulated emission of bremsstrahlung in a periodic magnetic field. J. Appl. Phys. 1971, 42, 1906–1913. [Google Scholar] [CrossRef]
  11. Colson, W.B. Free Electron Laser Theory; Technical Report; Berkeley Research Associates, Inc.: Richmond, CA, USA, 1986. [Google Scholar]
  12. Bonifacio, R.; Pellegrini, C.; Narducci, L.M. Collective instabilities and high-gain regime in a free electron laser. Opt. Commun. 1984, 50, 373–378. [Google Scholar] [CrossRef] [Green Version]
  13. Bonifacio, R.; Casagrande, F. Classical and quantum treatment of amplifier and superradiant free-electron laser dynamics. J. Opt. Soc. Am. B 1985, 2, 250–258. [Google Scholar] [CrossRef]
  14. Bonifacio, R.; McNeil, B.; Pierini, P. Superradiance in the high-gain free-electron laser. Phys. Rev. A 1989, 40, 4467. [Google Scholar] [CrossRef]
  15. Bonifacio, R.; Piovella, N.; McNeil, B. Superradiant evolution of radiation pulses in a free-electron laser. Phys. Rev. A 1991, 44, R3441. [Google Scholar] [CrossRef]
  16. Bonifacio, R.; Casagrande, F.; Cerchioni, G.; de Salvo Souza, L.; Pierini, P.; Piovella, N. Physics of the high-gain FEL and superradiance. Riv. Nuovo Cimento 1990, 13, 1–69. [Google Scholar] [CrossRef]
  17. Bonifacio, R.; Corsini, R.; De Salvo, L.; Pierini, P.; Piovella, N. New effects in the physics of high-gain free-electron lasers; a proposed experiment and possible applications. Riv. Nuovo Cimento 1992, 15, 1. [Google Scholar] [CrossRef]
  18. Dicke, R.H. Coherence in Spontaneous Radiation Processes. Phys. Rev. 1954, 93, 99–100. [Google Scholar] [CrossRef] [Green Version]
  19. Watanabe, T.; Wang, X.; Murphy, J.; Rose, J.; Shen, Y.; Tsang, T.; Giannessi, L.; Musumeci, P.; Reiche, S. Experimental characterization of superradiance in a single-pass high-gain laser-seeded free-electron laser amplifier. Phys. Rev. Lett. 2007, 98, 034802. [Google Scholar] [CrossRef]
  20. Bonifacio, R.; De Salvo, L.; Pierini, P.; Piovella, N.; Pellegrini, C. Spectrum, temporal structure, and fluctuations in a high-gain free-electron laser starting from noise. Phys. Rev. Lett. 1994, 73, 70. [Google Scholar] [CrossRef]
  21. Yu, L.H. Generation of intense uv radiation by subharmonically seeded single-pass free-electron lasers. Phys. Rev. A 1991, 44, 5178. [Google Scholar] [CrossRef]
  22. Allaria, E.; Appio, R.; Badano, L.; Barletta, W.; Bassanese, S.; Biedron, S.; Borga, A.; Busetto, E.; Castronovo, D.; Cinquegrana, P.; et al. Highly coherent and stable pulses from the FERMI seeded free-electron laser in the extreme ultraviolet. Nat. Photonics 2012, 6, 699–704. [Google Scholar] [CrossRef]
  23. Bonifacio, R.; Piovella, N.; Robb, G.; Schiavi, A. Quantum regime of free electron lasers starting from noise. Phys. Rev. Spec. Top. Accel. Beams 2006, 9, 090701. [Google Scholar] [CrossRef] [Green Version]
  24. Bonifacio, R.; Piovella, N.; Cola, M.; Volpe, L.; Schiavi, A.; Robb, G. The quantum free-electron laser. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 2008, 593, 69–74. [Google Scholar] [CrossRef] [Green Version]
  25. Zhao, Z.; Wang, D.; Gu, Q.; Yin, L.; Gu, M.; Leng, Y.; Liu, B. Status of the SXFEL Facility. Appl. Sci. 2017, 7, 607. [Google Scholar] [CrossRef] [Green Version]
  26. Kling, P.; Giese, E.; Endrich, R.; Preiss, P.; Sauerbrey, R.; Schleich, W.P. What defines the quantum regime of the free-electron laser? New J. Phys. 2015, 17, 123019. [Google Scholar] [CrossRef]
  27. Scully, M.O.; Zubairy, M.S. Quantum Optics; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  28. Bonifacio, R.; Casagrande, F.; Pellegrini, C. Hamiltonian model of a free electron laser. Opt. Commun. 1987, 61, 55–60. [Google Scholar] [CrossRef]
  29. Tran, T.M.; Wurtele, J.S. Free-electron laser simulation techniques. Phys. Rep. 1990, 195, 1–21. [Google Scholar] [CrossRef]
  30. Volpe, L. 3D Quantum Theory of Free Electron Lasers. Ph.D. Thesis, Università degli Studi di Milano, Milan, Italy, 2007. [Google Scholar]
  31. Moore, G.T.; Scully, M.O. Coherent dynamics of a free-electron laser with arbitrary magnet geometry. I. General formalism. Phys. Rev. A 1980, 21, 2000. [Google Scholar] [CrossRef]
  32. Murphy, J.; Pellegrini, C.; Bonifacio, R. Collective instability of a free electron laser including space charge and harmonics. Opt. Commun. 1985, 53, 197–202. [Google Scholar] [CrossRef]
  33. Bonifacio, R.; Maroli, C.; Piovella, N. Slippage and superradiance in the high-gain FEL: Linear theory. Opt. Commun. 1988, 68, 369–374. [Google Scholar] [CrossRef]
  34. Bonifacio, R.; Souza, L.D.S.; Pierini, P.; Piovella, N. The superradiant regime of a FEL: Analytical and numerical results. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 1990, 296, 358–367. [Google Scholar] [CrossRef]
  35. Giannessi, L.; Bellaveglia, M.; Chiadroni, E.; Cianchi, A.; Couprie, M.; Del Franco, M.; Di Pirro, G.; Ferrario, M.; Gatti, G.; Labat, M.; et al. Superradiant cascade in a seeded free-electron laser. Phys. Rev. Lett. 2013, 110, 044801. [Google Scholar] [CrossRef] [Green Version]
  36. Piovella, N. A hyperbolic secant solution for the superradiance in free electron lasers. Opt. Commun. 1991, 83, 92–96. [Google Scholar] [CrossRef]
  37. Orzechowski, T.; Anderson, B.; Fawley, W.; Prosnitz, D.; Scharlemann, E.; Yarema, S.; Hopkins, D.; Paul, A.; Sessler, A.; Wurtele, J. Microwave radiation from a high-gain free-electron laser amplifier. Phys. Rev. Lett. 1985, 54, 889. [Google Scholar] [CrossRef]
  38. Emma, P.; Akre, R.; Arthur, J.; Bionta, R.; Bostedt, C.; Bozek, J.; Brachmann, A.; Bucksbaum, P.; Coffee, R.; Decker, F.J.; et al. First lasing and operation of an ångstrom-wavelength free-electron laser. Nat. Photonics 2010, 4, 641–647. [Google Scholar] [CrossRef]
  39. Altarelli, M.; Brinkmann, R.; Chergui, M. The European X-ray Free-Electron Laser; Technical Design Report; DESY XFEL Project Group, European XFEL Project Team, Deutsches Elektronen-Synchrotron: Hamburg, Germany, 2007. [Google Scholar]
  40. Ishikawa, T.; Aoyagi, H.; Asaka, T.; Asano, Y.; Azumi, N.; Bizen, T.; Ego, H.; Fukami, K.; Fukui, T.; Furukawa, Y.; et al. A compact X-ray free-electron laser emitting in the sub-ångström region. Nat. Photonics 2012, 6, 540–544. [Google Scholar] [CrossRef]
  41. Saldin, E.L.; Schneidmiller, E.A.; Yurkov, M.V. Statistical and coherence properties of radiation from X-ray free-electron lasers. New J. Phys. 2010, 12, 035010. [Google Scholar] [CrossRef]
  42. Pellegrini, C. Progress toward a soft X-ray FEL. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 1988, 272, 364–367. [Google Scholar] [CrossRef]
  43. Moore, G.T. The high-gain regime of the free electron laser. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 1985, 239, 19–28. [Google Scholar] [CrossRef]
  44. Scharlemann, E.; Sessler, A.M.; Wurtele, J. Optical guiding in a free electron laser. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 1985, 239, 29–35. [Google Scholar] [CrossRef] [Green Version]
  45. Bonifacio, R.; Piovella, N.; Cola, M.; Volpe, L. Experimental requirements for X-ray compact free electron lasers with a laser wiggler. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 2007, 577, 745–750. [Google Scholar] [CrossRef]
  46. Gea-Banacloche, J.; Moore, G.; Schlicher, R.; Scully, M.; Walther, H. Soft X-ray free-electron laser with a laser undulator. IEEE J. Quantum Electron. 1987, 23, 1558–1570. [Google Scholar] [CrossRef]
  47. Dobiasch, P.; Meystre, P.; Scully, M. Optical wiggler free-electron X-ray laser in the 5 Å region. IEEE J. Quantum Electron. 1983, 19, 1812–1820. [Google Scholar] [CrossRef]
  48. Steiniger, K.; Bussmann, M.; Pausch, R.; Cowan, T.; Irman, A.; Jochmann, A.; Sauerbrey, R.; Schramm, U.; Debus, A. Optical free-electron lasers with traveling-wave Thomson-scattering. J. Phys. B At. Mol. Opt. Phys. 2014, 47, 234011. [Google Scholar] [CrossRef]
  49. Reiche, S. GENESIS 1.3: A fully 3D time-dependent FEL simulation code. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 1999, 429, 243–248. [Google Scholar] [CrossRef]
  50. Fawley, W.M. A User Manual for GINGER and Its Post-Processor XPLOTGIN; Technical Report; Lawrence Berkeley National Lab (LBNL): Berkeley, CA, USA, 2002.
  51. Saldin, E.; Schneidmiller, E.; Yurkov, M. FAST: A three-dimensional time-dependent FEL simulation code. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 1999, 429, 233–237. [Google Scholar] [CrossRef]
  52. Campbell, L.; McNeil, B. Puffin: A three dimensional, unaveraged free electron laser simulation code. Phys. Plasmas 2012, 19, 093119. [Google Scholar] [CrossRef] [Green Version]
  53. Piovella, N.; Cola, M.; Volpe, L.; Schiavi, A.; Bonifacio, R. Three-dimensional Wigner-function description of the quantum free-electron laser. Phys. Rev. Lett. 2008, 100, 044801. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Debus, A.; Steiniger, K.; Kling, P.; Carmesin, C.M.; Sauerbrey, R. Realizing quantum free-electron lasers: A critical analysis of experimental challenges and theoretical limits. Phys. Scr. 2019, 94, 074001. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Scheme of free-electron laser. Here the undulator is built by permanent magnets, where the period corresponds to the distance between two red arrows. The electrons wiggle along the undulator axis and emit a radiation burst in the forward direction.
Figure 1. Scheme of free-electron laser. Here the undulator is built by permanent magnets, where the period corresponds to the distance between two red arrows. The electrons wiggle along the undulator axis and emit a radiation burst in the forward direction.
Atoms 09 00028 g001
Figure 2. Typical SASE spectrum, in this case from the simulations of an experiment at λ = 2 nm (from [25]).
Figure 2. Typical SASE spectrum, in this case from the simulations of an experiment at λ = 2 nm (from [25]).
Atoms 09 00028 g002
Figure 3. | Im λ | vs. δ and different values of μ : (a) 0, (b) 0.5 , (c) 3, (d) 5, (e) 7 and (f) 10.
Figure 3. | Im λ | vs. δ and different values of μ : (a) 0, (b) 0.5 , (c) 3, (d) 5, (e) 7 and (f) 10.
Atoms 09 00028 g003
Figure 4. Dimensionless radiation intensity | A | 2 vs. z / l g for δ = 0 and A 0 = 10 4 .
Figure 4. Dimensionless radiation intensity | A | 2 vs. z / l g for δ = 0 and A 0 = 10 4 .
Atoms 09 00028 g004
Figure 5. Strong superradiance, with l b = 30 l c , δ = 0 and different values of z ¯ : (a) z ¯ = 8 ; (b) z ¯ = 17 ; (c) z ¯ = 34 .
Figure 5. Strong superradiance, with l b = 30 l c , δ = 0 and different values of z ¯ : (a) z ¯ = 8 ; (b) z ¯ = 17 ; (c) z ¯ = 34 .
Atoms 09 00028 g005
Figure 6. Weak superradiance, with l b = 0.1 l c , δ = 0 and different values of z ¯ : (a) z ¯ = 50 ; (b) z ¯ = 75 ; (c) z ¯ = 125 .
Figure 6. Weak superradiance, with l b = 0.1 l c , δ = 0 and different values of z ¯ : (a) z ¯ = 50 ; (b) z ¯ = 75 ; (c) z ¯ = 125 .
Atoms 09 00028 g006
Figure 7. Strong pure superradiance, with l b = 30 l c , δ = 2 and different values of z ¯ : (a) z ¯ = 11 ; (b) z ¯ = 21 ; (c) z ¯ = 34 .
Figure 7. Strong pure superradiance, with l b = 30 l c , δ = 2 and different values of z ¯ : (a) z ¯ = 11 ; (b) z ¯ = 21 ; (c) z ¯ = 34 .
Atoms 09 00028 g007
Figure 8. Superradiant self-similar solution | A 1 | 2 vs. y = z 1 ( z ¯ z 1 ) .
Figure 8. Superradiant self-similar solution | A 1 | 2 vs. y = z 1 ( z ¯ z 1 ) .
Atoms 09 00028 g008
Figure 9. SASE mode operation: typical temporal pattern of the intensity | A r | 2 vs. z 1 for a long beam with l b = 50 l c , at different undulator lengths z ¯ : (a) z ¯ = 0 ; (b) z ¯ = 50 ; (c) z ¯ = 100 .
Figure 9. SASE mode operation: typical temporal pattern of the intensity | A r | 2 vs. z 1 for a long beam with l b = 50 l c , at different undulator lengths z ¯ : (a) z ¯ = 0 ; (b) z ¯ = 50 ; (c) z ¯ = 100 .
Atoms 09 00028 g009
Figure 10. SASE mode operation: power spectrum of the intensity shown in Figure 9 vs. ω ¯ = ( ω ω s ) / 2 ρ ω s ; (a) z ¯ = 0 ; (b) z ¯ = 50 ; (c) z ¯ = 100 .
Figure 10. SASE mode operation: power spectrum of the intensity shown in Figure 9 vs. ω ¯ = ( ω ω s ) / 2 ρ ω s ; (a) z ¯ = 0 ; (b) z ¯ = 50 ; (c) z ¯ = 100 .
Atoms 09 00028 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Piovella, N.; Volpe, L. A Review of High-Gain Free-Electron Laser Theory. Atoms 2021, 9, 28. https://doi.org/10.3390/atoms9020028

AMA Style

Piovella N, Volpe L. A Review of High-Gain Free-Electron Laser Theory. Atoms. 2021; 9(2):28. https://doi.org/10.3390/atoms9020028

Chicago/Turabian Style

Piovella, Nicola, and Luca Volpe. 2021. "A Review of High-Gain Free-Electron Laser Theory" Atoms 9, no. 2: 28. https://doi.org/10.3390/atoms9020028

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop