Download Free Electron Lasers Introduction Undulator Radiation Low

Free Electron Lasers
Peter Schmüser
University of Hamburg
Introduction
Undulator Radiation
Low-Gain Free Electron Laser
High-Gain FEL, SASE Principle
Experimental Results
– Typeset by FoilTEX –
1
– Typeset by FoilTEX –
2
– Typeset by FoilTEX –
3
– Typeset by FoilTEX –
4
Big advantage of FEL: the wavelength is tunable by changing the electron energy.
– Typeset by FoilTEX –
5
Reflectivity at normal incidence (%)
An optical cavity is no longer possible for wavelengths below 100 nm.
100
Al
80
Ag
Au
60
40
SiC
20
0
10
2
3
4
5 6 7 8
2
3
100
4
5 6 7 8
1000
ncidence (%)
Photon wavelength (nm)
100
Al
80
– Typeset by FoilTEX –
Ag
6
Principle of a Self Amplified Spontaneous Emission (SASE) Free Electron Laser
– Typeset by FoilTEX –
7
Undulator Radiation
Electron motion in undulator
Schematic view of electron motion in an undulator magnet
Call W = Ekin + mec2 the total relativistic energy of the electron.
p
2
Lorentz factor, normalized velocity: γ = W/(mec ) , β = 1 − 1/γ 2
The average velocity in z direction is less than βc owing to the sinodoidal trajectory
(proof in FEL Course)
v̄z ≡ β̄ c =
1
eB0λu
2
1 − 2 (1 + K /2) c with K =
2γ
2πmec
K is called the undulator parameter
– Typeset by FoilTEX –
8
Synchrotron radiation is emitted inside a cone with opening angle 1/γ (see lecture by
Rivkin)
Undulator: K ≤ 1, the electron trajectory is inside the radiation cone. Therefore, the
photons emitted by a given electron at various positions along the undulator interfere with
each other. This has the important consequence that the radiation is monochromatic,
in contrast to synchrotron radiation in bending magnets.
Physical Processes
in aindependently
Free Electron Laserand do not interfere.
Note: different2. electrons
radiate
B
A
B
e–
1/γ
Figure 2.1.: Emission of radiation in an undulator.
In the TTF undulator, the deviation from the straight orbit is only 10 µm. Synchrotron radiation is emitted by relativistic electrons in a cone with opening angle
1/γ . In an undulator, the the maximum angle of the particle velocity with respect to
the undulator axis α = arctan(vx/vz ) is always smaller than the opening angle of the
radiation, therefore the radiation field may add coherently. In a wiggler, αmax > 1/γ,
– Typeset by FoilTEX –
9
Lorentz transformation to moving coordinate system
Consider a coordinate system (x∗, y ∗, z ∗) moving with the average velocity v̄z = β̄c
of the electron. The undulator period appears shortened due to the relativistic length
contraction
λ∗u = λu/γ
In the moving system the electron carries out a harmonic oscillation in x direction with
the frequency
ω ∗ = γ c 2π/λu
(Superimposed is a small longitudinal oscillation, which will be ignored here, it leads to
higher harmonics in the radiation).
In the moving system the electron emits dipole radiation with a frequency ω ∗ = γ ωu
(with ωu = c 2π/λu) and a wavelength λ∗ = λu/γ
Remember: λu is the undulator period, i.e. the distance between two equal poles.
Take typical values: λu = 25 mm, γ = 1000
⇒ λ∗ = 25 µm.
– Typeset by FoilTEX –
10
Transformation of radiation into laboratory system
We are interested in the wavelength of the light emitted in forward direction. The
Lorentz transformation of the photon energy (this can also be considered as the relativistic
Doppler shift) reads
~ω ∗ = γ~ω`
After a little algebra we get for the wavelength of the undulator light (see FEL Course)
λ` =
– Typeset by FoilTEX –
λu
2
(1
+
K
/2)
2
2γ
11
Properties of undulator radiation
An electron passing an undulator with Nu periods produces a wavetrain with Nu
oscillations.
Finite wave train
(here with 10 periods)
The time duration of the wave train is T = Nuλ`/c. Due to the finite duration the
1
radiation is not monochromatic but containsSpectral
a frequency
spectrum which is obtained by
intensity for a wave train
Fourier transformation (see FEL Course). with Nu = 100 periods
The spectral
intensity is
I( ω ) 0.5
0
2
sin
ξ
∆ω T
ω − ω`
0.96
0.98
1
1.02
1.04
I(ω) ∝
with ξ =
= π Nu
ξ
2
ω`
ω
ω0
It has a maximum at ω = ω` and a width proportional to 1/Nu.
– Typeset by FoilTEX –
12
1
Spectral intensity for a wave train
with Nu = 100 periods
I( ω ) 0.5
0
0.96 0.98
1
1.02 1.04
ω
ω0
In the figure, the normalized intensity is plotted as a function of ω/ω`. The total radiation
power (integrated over frequency) is the same as in a bending magnet (see Rivkin):
Prad
2α~c2γ 4
=
3ρ2
Main differences to synchrotron radiation in bending magnets: (a) undulator radiation
is confined to a narrow spectral line, (b) the radiation is well collimated. Note, however,
that different electrons radiate indepedently both in bending magnets and in undulators,
hence the intensity depends linearly on the number N of electrons per bunch:
PN = N · P1
– Typeset by FoilTEX –
13
Low-Gain FEL
Energy transfer from electron to light wave
Consider “seeding” by an external light source with wavelength λ`
The light wave is co-propagating with the relativistic electron beam. It is described by a
plane electromagnetic wave
Ex(z, t) = E0 cos(k`z − ω`t) with k` = ω`/c
– Typeset by FoilTEX –
14
Question: can there be a continuous energy transfer from electron beam to light wave?
The electron energy W = γmec2 changes in time dt by
dW = ~v · F~ = −evx(t)Ex(t)dt
The average electron speed in z direction is v̄z = c 1 − 2γ1 2 (1 + K 2/2) < c
To determine the condition for resonant energy transfer we compute the electron and
light travel times for a half period of the undulator:
tel = λu/(2v̄z ), tlight = λu/(2c)
Continuous energy transfer happens if ω`(tel − tlight) = π
(Remark: also 3π, 5π . . . are possible, leading to higher harmonics of the radiation)
vx
electron trajectory
Ex
Ex
– Typeset by FoilTEX –
z
vx
light wave
15
From this condition we compute the light wavelength (see FEL Course)
λ` =
λu
2γ 2
2
1+
K
2
This wavelength is identical with the undulator radition wavelength (in forward direction).
The quantitative treatment of the energy transfer from the electron to the light wave
is presented in the FEL Course. Here I quote the results. Introducing so so-called
ponderomotive phase:
ψ ≡ (k` + ku)z − ω`t
one can show that the time-variation of the electron γ factor and of the phase are
dγ
eE0K
=−
sin ψ
2
dt
2mecγr
dψ
γ − γr
= 2kuc
dt
γr
where the “resonant” gamma-factor is defined by the condition
2
λu
K
λ` = 2 1 +
2γr
2
– Typeset by FoilTEX –
16
Note: if the electron energy is equal to Er = γr mec2 then the undulator radiation
produced by the electron beam has exactly the wavelength of the seed laser. In an FEL,
however, one has to run the e-beam at a slightly higher energy E = γmec2 > Er in
order to amplify the light wave.
The combination of the two first order equations yields the “Pendulum Equation” of the
low-gain FEL
eE0Kku
2
2
ψ̈ + Ω sin ψ = 0
with Ω =
meγr2
– Typeset by FoilTEX –
17
Phase space representation
The two equations
dγ
eE0K
=−
sin ψ
2
dt
2mecγr
dψ
γ − γr
= 2kuc
dt
γr
can be used to plot the trajectories in the (ψ, γ) phase space. There is a close analogy
with the motion of a mathematical pendulum. At small amplitude we get a harmonic
oscillation. With increasing angular momentum the motion becomes unharmonic. At
very large angular momentum one gets a rotation (unbounded motion).
ψ ∼ (γ − γr)
Rotation
ψ
0
ψ
–π
– Typeset by FoilTEX –
0
ψ
π
Oscillation
18
– Typeset by FoilTEX –
19
The..mathematical
10 , 9.97
8
treatment of the energy transfer between electron and light wave is
quite involved,
see the FEL Course. The essential results is that the FEL-gain is given
2
sin( ξ )
I ( ξby
)
2
sin(
sin(
ξ
)
ξ
)
2
2
2
3
2
ξ
G
2 .u λu n.ecos(dξ ) sin
2. ξ
π( eξ ) K N
ω − ω`
2
3
G(ξ) = −
·
with ξ = π Nu
4ε0mec2γξ r3
dξ
ξ2 ξ
ω`
Madey Theorem
The FEL gain curve is obtained by taking the negative derivative of the line-shape curve
of undulator radiation.
spectral line of undulator
gain of FEL
1
1
0.5
G( ξ )
I( ξ ) 0.5
0
0
0.5
0
10
0
ξ
– Typeset by FoilTEX –
10
1
10
0
10
ξ
20
High-Gain FEL
The essential feature of the high-gain FEL is that a large number of electrons radiate
coherently. In that case, the intensity of the radiation field grows quadratically with the
number of particles: IN = N 2 I1.
– Typeset by FoilTEX –
21
Big problem: the particle bunches are much longer than the FEL wavelength, it appears
impossible to produce intense electron bunches with a length λ`.
The way out of this dilemma is given by the process of microbunching:
Electrons which lose energy to the light wave travel a longer path in the undulator,
electrons which gain energy from the light wave travel a shorter path. The result is a
modulation of the longitudinal velocity. This velocity modulation leads eventually to a
concentration of the electrons in slices which are much shorter than λ` .
a)
0.2
0
–0.1
–0.2
0.2
0.1
x / mm
x / mm
0.1
b)
0
–0.1
–2π
0
ψ
2π
4π
–0.2
c)
0.1
x / mm
0.2
0
–0.1
–2π
0
ψ
2π
4π
–0.2
–2π
0
ψ
2π
4π
The particles within a micro-bunch radiate coherently. The resulting strong radiation
field enhances the micro-bunching even further.
Result: “collective instability”, exponential growth of radiation power
– Typeset by FoilTEX –
22
The ultimate power is P ∝ Nc2 where Nc is the number of particles in a coherence
region
typical value Nc ≈ 106 ⇒ PF EL = 106Pundulator
Coherent action is what counts:
– Typeset by FoilTEX –
23
An approximate analytic description of the high-gain FEL requires the self-consistent
solution of the coupled pendulum equations and the inhomogeneous wave equation for
the electromagnetic field of the light wave.
In the 1D-FEL theory the dependencies on the transverse coordinates x, y are
disregarded.
The wave equation for the radiation field Ex reads
∂ 2 Ex
1 ∂ 2 Ex
∂jx
−
=
µ
0
∂z 2
c2 ∂t2
∂t
where the current density ~j is generated by the electron bunch moving on its cosine-like
trajectory.
In addition, one has to consider the longitudinal space charge field Ez which is generated
by the gradually evolving periodic charge density modulation. After a lot of tedious
mathematical steps and several simplifying assumptions one arrives at a third-order
differential equation for the “slowly varying amplitude” of the electric field of the light
wave:
d3Ẽx
d2Ẽx
2 2 dẼx
3
−
4ik
η
−
4k
η
−
iΓ
Ẽx(z) = 0
u
u
dz 3
dz 2
dz
– Typeset by FoilTEX –
24
Here we have introduced the gain parameter Γ and the relative energy deviation η
Γ=
2 2
µ0K e kune
4γ 3me
13
γ − γr
η=
γr
and assumed that the electron beam has negligible energy spread.
This third-order differential equation can be solved analytically. For the case γ = γr one
obtains
!
!
√
√
i+ 3
i− 3
Γz + A3 exp
Γz
Ẽx(z) = A1 exp (−iΓz) + A2 exp
2
2
The second term exhibits exponential growth as a function
√ of the position z in the
3
undulator.
The
electric
field
grows
exponentially
as
exp(
2 Γz), the power grows as
√
exp( 3Γz).
The gain parameter Γ is related to two parameters which are in widespread use: the
Pierce parameter and the power gain length
ρpierce =
– Typeset by FoilTEX –
λu Γ
4π
1
Lg = √
3Γ
25
The above calculations, which have been sketched only very briefly, indicate that there
is an onset of an “instability”, leading to a progressing microbunching and an exponential
increase in radiation power along the undulator. A quantitave treatment requires
elaborate numerical simulations.
– Typeset by FoilTEX –
26
Simulation of microbunching
a)
0.2
0
–0.1
–0.2
0.2
0.1
x / mm
x / mm
0.1
b)
0
–0.1
–2π
0
ψ
2π
4π
–0.2
c)
0.1
x / mm
0.2
0
–0.1
–2π
0
ψ
2π
4π
–0.2
–2π
0
ψ
2π
4π
Experimental observation of microbunching at the 60 µm FEL Firefly, Stanford
– Typeset by FoilTEX –
27
Exponential growth of radiation power and progressing microbunching in a long undulator
– Typeset by FoilTEX –
28
First lasing at the TESLA Test Facility (TTF) Free Electron Laser
Bjorn Wiik Price 2000 for Evgeni Saldin, Evgeni Schneidmiller, Mikhail Yurkov
– Typeset by FoilTEX –
29
The FEL as a wavelength tunable laser:
λu
λ = 2 (1 + K 2/2)
2γ
– Typeset by FoilTEX –
30
SASE FELs:
State of the art
All observations agree with
theor. expectations/
computer models
VISA
ATF/BNL
840 nm
September 2000
TTF-FEL
DESY
80-120 nm
March 2001
– Typeset by FoilTEX –
LEUTL
APS/ANL
385 nm
September 2001
31
3.7.
Life Sciences
V-157
Figure 3.7.5.: Coulomb explosion of a T4 lysozyme molecule (H: white; C: grey; N: blue; O: red;
S: yellow) induced by the radiation damage caused by a 3 × 10 12 photon per (0.1 µm) 2 pulse of
12.4 keV energy. The FWHM of the pulse was 50 fs. The molecule is shown at the beginning, in
the middle and after the pulse. Even after half of the pulse has passed, the distortions are small.
After the pulse the Coulomb explosion is under way [ 7]. The distortion of the molecule during the
time of the pulse is considerably smaller for lower flux densities during the pulse (see also Fig.
3.7.4).
– Typeset by FoilTEX –
32
– Typeset by FoilTEX –
33