Download JUAS_2014_RB_synchrotron_radiation_I - Indico

Synchrotron radiation
R. Bartolini
John Adams Institute for Accelerator Science, University of Oxford
and
Diamond Light Source
JUAS 2014
27-31 January 2014
Contents
Introduction to synchrotron radiation
properties of synchrotron radiation
synchrotron light sources
angular distribution of power radiated by accelerated particles
angular and frequency distribution of energy radiated
radiation from undulators and wigglers
Storage ring light sources
electron beam dynamics in storage rings
radiation damping and radiation excitation
emittance and brilliance
low emittance lattices
diffraction limited storage rings
R. Bartolini, JUAS, 27-31 January 2014
2/71
What is synchrotron radiation
Electromagnetic radiation is emitted by charged particles when accelerated
The electromagnetic radiation emitted when the charged particles are
accelerated radially (v  a) is called synchrotron radiation
It is produced in the synchrotron radiation sources using bending magnets
undulators and wigglers
R. Bartolini, JUAS, 27-31 January 2014
3/71
Synchrotron radiation sources properties (I)
Broad Spectrum which covers from microwaves to hard X-rays:
the user can select the wavelength required for experiment;
either with a monochromator
or adjusting the emission wavelength of insertion devices
synchrotron light
R. Bartolini, JUAS, 27-31 January 2014
4/71
Synchrotron radiation sources properties (II)
High Flux: high intensity photon beam, allows rapid experiments or use of
weakly scattering crystals;
Flux = Photons / ( s  BW)
High Brilliance (Spectral Brightness): highly collimated photon beam
generated by a small divergence and small size source
Brilliance = Photons / ( s 
mm2

mrad2
 BW )
Partial coherence in SRs
Full T coherence in FELs
Polarisation: both linear and circular (with IDs)
10s ps in SRs
Pulsed Time Structure: pulsed length down to
High Stability: submicron source stability in SR
… and it can be computed!
10s fs in FELs
diamond
1.E+20
1.E+18
1.E+16
2
2
Brightness (Photons/sec/mm /mrad /0.1%)
Peak Brilliance
1.E+14
1.E+12
X-rays from Diamond
will be 1012 times
brighter than from
an X-ray tube,
105 times brighter
than the SRS !
1.E+10
1.E+08
1.E+06
X-ray
tube
60W bulb
Candle
1.E+04
1.E+02
R. Bartolini, JUAS, 27-31 January 2014
6/71
X-ray sources
Courtesy K. Wille
R. Bartolini, JUAS, 27-31 January 2014
7/71
Layout of a synchrotron radiation source (I)
Electrons are generated and
accelerated in a linac, further
accelerated to the required energy in
a booster and injected and stored in
the storage ring
The circulating electrons emit an
intense beam of synchrotron radiation
which is sent down the beamline
R. Bartolini, JUAS, 27-31 January 2014
8/71
Layout of a synchrotron radiation source (II)
Evolution of synchrotron radiation sources (I)
• First observation:
1947, General Electric, 70 MeV synchrotron
• First user experiments:
1956, Cornell, 320 MeV synchrotron
• 1st generation light sources: machine built for High
Energy Physics or other purposes used parasitically for
synchrotron radiation
• 2nd generation light sources: purpose built synchrotron
light sources, SRS at Daresbury was the first dedicated
machine (1981 – 2008)
• 3rd generation light sources: optimised for high brilliance
with low emittance and Insertion Devices; ESRF, Diamond,
…
R. Bartolini, JUAS, 27-31 January 2014
10/71
Evolution of synchrotron radiation sources (II)
• 4th generation light sources: photoinjectors LINAC based Free Electron
Laser sources;
FLASH (DESY) 2007
LCLS (SLAC) 2009
SACLA (Japan) 2011
Elettra (Italy) 2012
and in the near(?) future
• 4th generation light sources storage ring based: diffraction limited storage
rings
• …and even a 5th generation with more compact and advanced accelerator
technologies e.g. based on laser plasma wakefield accelerators
R. Bartolini, JUAS, 27-31 January 2014
11/71
3rd generation storage ring light sources
1992
1993
1994
1996
1997
1998
2000
2004
2006:
2008
2009
2011
ESRF, France (EU)
ALS, US
TLS, Taiwan
ELETTRA, Italy
PLS, Korea
MAX II, Sweden
APS, US
LNLS, Brazil
Spring-8, Japan
BESSY II, Germany
ANKA, Germany
SLS, Switzerland
SPEAR3, US
CLS, Canada
SOLEIL, France
DIAMOND, UK
ASP, Australia3 GeV
MAX III, Sweden
Indus-II, India
SSRF, China
PETRA-III, Germany
ALBA, Spain
6 GeV
1.5-1.9 GeV
1.5 GeV
2.4 GeV
2 GeV
1.5 GeV
7 GeV
1.35 GeV
8 GeV
1.9 GeV
2.5 GeV
2.4 GeV
3 GeV
2.9 GeV
2.8 GeV
3 GeV
700 MeV
2.5 GeV
3.4 GeV
6 GeV
3 GeV
ESRF
SSRF
3rd generation storage ring light sources
under construction or planned
> 2013
NLSL-II
NSLS-II, US
3 GeV
MAX-IV, Sweden
SOLARIS, Poland
SESAME, Jordan
1.5-3 GeV
1.5 GeV
2.5 GeV
TPS, Taiwan
CANDLE, Armenia
3 GeV
3 GeV
PEP-X, USA
Spring8-II, Japan
4.5 GeV
6 GeV
Max-IV
R. Bartolini, JUAS, 27-31 January 2014
13/71
Diamond aerial views
June 2003
Oct 2006
Main components of a storage ring
Dipole magnets to bend the electrons
Sextupole magnets to focus off-energy
electrons (mainly)
Quadrupole magnets to focus the electrons
RF cavities to replace energy losses due to
the emission of synchrotron radiation
Main components of a storage ring
Insertion devices (undulators) to
generate high brilliance radiation
R. Bartolini, JUAS, 27-31 January 2014
Insertion devices (wiggler) to
reach high photon energies
16/71
Accelerator physics and technology challenges
Photon energy
Brilliance
Flux
Stability
Polarisation
Time structure
R. Bartolini, JUAS, 27-31 January 2014
Ring energy
Small Emittance
Insertion Devices
High Current; Feedbacks
Vibrations; Orbit Feedbacks; Top-Up
Short bunches; Short pulses
17/71
Brilliance and low emittance
The brilliance of the photon beam is determined (mostly) by the electron beam
emittance that defines the source size and divergence
brilliance 
flux
4 2 x  x ' y  y '
2
 x   x2,e   ph
,e
 x '   x2',e   '2ph,e
 x   x  x  ( Dx  ) 2
 x'   x  x  ( D' x   ) 2
Brilliance with IDs
Thanks to the progress with IDs technology storage ring light sources can cover
a photon range from few tens of eV to tens 10 keV or more with high brilliance
Medium energy storage rings with In-vacuum undulators operated at low gaps (e.g.
5-7 mm) can reach 10 keV with a brilliance of 1020 ph/s/0.1%BW/mm2/mrad2
R. Bartolini, JUAS, 27-31 January 2014
19/71
Many ways to
use x-rays
photo-emission
(electrons)
electronic structure
& imaging
diffraction
crystallography
& imaging
scattering
SAXS
& imaging
absorption
from the synchrotron
to the detector
Spectroscopy
EXAFS
XANES
& imaging
fluorescence EXAFS
XRF imaging
R. Bartolini, JUAS, 27-31 January 2014
20/71
Applications
Medicine, Biology, Chemistry, Material Science, Environmental Science and more
Biology
Reconstruction of the 3D structure of a nucleosome
with a resolution of 0.2 nm
The collection of precise information on the molecular
structure of chromosomes and their components can
improve the knowledge of how the genetic code of
DNA is maintained and reproduced
R. Bartolini, JUAS, 27-31 January 2014
Archeology
A synchrotron X-ray beam at the SSRL facility
illuminated an obscured work erased, written over
and even painted over of the ancient
mathematical genius Archimedes, born 287 B.C.
in Sicily.
X-ray fluorescence imaging revealed the
hidden text by revealing the iron contained in
the ink used by a 10th century scribe. This xray image shows the lower left corner of the
page.
21/71
Life science examples: DNA and myoglobin
Franklin and Gosling used a X-ray tube:
Brilliance was 108 (ph/sec/mm2/mrad2/0.1BW)
Exposure times of 1 day were typical (105 sec)
e.g. Diamond provides a brilliance of 1020
100 ns exposure would be sufficient
Nowadays pump probe experiment in life science are
performed using 100 ps pulses from storage ring light
sources: e.g. ESRF myoglobin in action
Photograph 51
Franklin-Gosling
DNA (form B)
1952
Lienard-Wiechert potentials (II)
The general solutions for the wave equation driven by a time varying
charge and current density read (in the Lorentz gauge)
1
( x, t ) 
4 0

V
( x  x ' , t ret ) 3
d x'
| x  x'|
where ret means retarted
t ret  t 
A(x, t ) 
| x ( t )  x ( t ret ) |
c
0
4

V
J ( x  x ' , t ret ) 3
d x'
| x  x'|
(see next slide)
Using the properties of the Dirac deltas we can integrate and obtain the
Lienard-Wiechert potentials
1
( x, t ) 
4 0


e


 (1    n )R  ret
A( x, t ) 

1 
e


4 0 c  (1    n )R  ret
These are the potentials of the em fields generated by the charged particle
in motion.
The trajectory itself is determined by external electric and magnetic fields
R. Bartolini, JUAS, 27-31 January 2014
23/71
Lienard-Wiechert Potentials (III)
[ ]ret means computed at time t’
t  t '
R (t ' )
c
Potentials and fields at position x at time t are determined by the
characteristic of the electron motion at a time t’
t’ – t is the time it takes for the em radiation to travel the distance R(t’)
i.e. grey is the position of the electron at time t
R. Bartolini, JUAS, 27-31 January 2014
24/71
The Lienard-Wiechert fields
The electric and magnetic fields are computed from the potentials using
B  A
E
A
 
t
and are called Lienard-Wiechert fields
The computation has to be done carefully since the potentials depend on t’,
not t. The factor dt/dt’ represents the Doppler factor. We get

e 
n 
e  n  (n   )   
E (x, t) 




2
3 2
40   (1    n ) R  ret 40 c  (1    n ) 3 R 
ret
velocity field 
1
R
2
acceleration field
If we consider the acceleration field we have
B ( x ,t ) 

1

nE
c
1
R

E  B  nˆ
and the correct dependence as 1/R as for radiation field
R. Bartolini, JUAS, 27-31 January 2014
25/71

rit
Power radiated
Power radiated by a particle on a surface is the flux of the Poynting vector
S
1
0
E B
  (S )(t )   S ( x , t )  n d

Angular distribution of radiated power
d 2P
 ( S  n )(1  n   )R 2
d
radiation emitted by the particle
We will analyse two cases:
acceleration orthogonal to the velocity → synchrotron radiation
acceleration parallel to the velocity → bremmstrahlung
R. Bartolini, JUAS, 27-31 January 2014
26/71
Synchrotron radiation: non relativistic motion (I)
Assuming
  0 and substituting the acceleration field
 n  (n   ) 
Eacc ( x , t ) 


40 c 
R
 ret
e
The angular distribution of the power radiated is
given by
2
d 2P
1
e2
)

REacc 
n

(
n


d  0 c
4 2  0 c
2
Working out the double cross product
We have
R. Bartolini, JUAS, 27-31 January 2014
27/71
Synchrotron radiation: non relativistic motion (II)
Since
where  is the angle between the acceleration and the observation
direction, we finally get
The angular distribution of power reads
2
d 2P
e2

2


sin

2
d (4 )  0c
R. Bartolini, JUAS, 27-31 January 2014
28/71
Synchrotron radiation: non relativistic motion (III)
Integrating over the angles gives the total radiated power
This integral gives
P
e2
60 c

2
Larmor’s formula
It shows that radiation is emitted when the particle is accelerated.
Using

p
mc
we have (to be used later for the generalisation to the relativistic case)
e2
dp
P
6 0 m 2 c 3 dt
2
R. Bartolini, JUAS, 27-31 January 2014
29/71
Synchrotron radiation: relativistic motion (I)
In the relativistic case the total radiated power is computed in the same way.
Using only the acceleration field (large R)

e  n  (n   )   


E ( x, t ) ~
3
4 0 c  (1    n ) R 

 ret
The angular distribution of the power emitted is (use the retarted time!)
d 2P
1

RE acc
d  o c

n  ( n   )   
2


2
e
(1    n ) 
(1  n   ) 5
42  0 c
2
The emission is peaked in the direction of the velocity
The pattern depends on the details of velocity and acceleration but it is
dominated by the denominator
R. Bartolini, JUAS, 27-31 January 2014
30/71
velocity  acceleration: synchrotron radiation (I)
Assuming
  
d 2P
e2

d 4 2  0 c 2
and substituting the acceleration field

n  (n   )  
(1  n   )5

2
e 
2
2
1

4 2  0c (1   cos  )3

sin 2  cos 2  
1  2
2
  (1   cos  ) 
cone aperture
 1/
When the electron velocity approaches the speed of light, the emission
pattern is sharply collimated forward
R. Bartolini, JUAS, 27-31 January 2014
31/71
velocity  acceleration: synchrotron radiation (II)
Courtesy K. Wille
R. Bartolini, JUAS, 27-31 January 2014
32/71
Total radiated power via synchrotron radiation
Integrating over the whole solid angle we obtain the total instantaneous
power radiated by one electron
2
e2  2 E
e2
dp 2 e 2c  4
e4

4
2 2
P
  





E
B
4
2 3
2
4 5
60c
60c
Eo
60 m c dt
60 
60 m c
e2
4
2
• Strong dependence 1/m4 on the rest mass
• proportional to 1/2 ( is the bending radius)
• proportional to B2 (B is the magnetic field of the bending dipole)
The radiation power emitted by an electron beam in a storage ring is very high.
The surface of the vacuum chamber hit by synchrotron radiation must be
cooled.
R. Bartolini, JUAS, 27-31 January 2014
33/71
velocity || acceleration: bremsstrahlung
Assuming
 | | 
dP
e2

d 4 2  0 c
and substituting the acceleration field

n  ( n   )  
( 1  n   )5

2
e 
2
2
sin 2 

4 2  0 c ( 1   cos  )5
 1
 3
 max  arccos 
 1  15  1  21
2
Integrating over the angles as before gives the total radiated power
e2
dp
P
6 0 m 2 c 3 dt
2
R. Bartolini, JUAS, 27-31 January 2014
34/71
Comparison of radiation from linear and
circular trajectories
Back to the general expression for the acceleration field, integrating over
the angles gives the total radiated power
P
e2
60c

 6 (  ) 2  (    ) 2

Relativistic generalization of
Larmor’s formula
The total radiated power can also be computed by relativistic transformation
of the 4-acceleration in Larmor’s formula
e2
dp
P
6 0 m 2 c 3 dt
2
with
t
and
dp
dp

dt
dt
We build the relativistic invariant
 dp
P

2 3
6 0 m c  d
e2
2
1  dE 
 2

c  d 
R. Bartolini, JUAS, 27-31 January 2014
2


35/71
Comparison of radiation from linear and
circular trajectories
The particle energy is
E 2  (mc2 ) 2  p 2c2
Therefore
dE
dp
E
 c2p
d
d
and
dE
dp
v
d
d
Inserting in the formula we get
For a linear trajectory
e2
dp
P
6 0 m 2 c 3 dt
2
R. Bartolini, JUAS, 27-31 January 2014
36/71
Comparison of radiation from linear and
circular trajectories
Repeating for a circular trajectory
 dp
P

2 3
6 0 m c  d
e2
2
1  dE 
 2

c  d 
2


The particle energy is now constant
E 2  (mc2 ) 2  p 2c2
Therefore
dE
0
d
Inserting in the formula we get the total radiated power in a circular
trajectory
2
This is 2 larger than
2
e
2 dp
P

the linear case
2 3
dt
6 0 m c
P(v || a)  1/2 P(v  a)
R. Bartolini, JUAS, 27-31 January 2014
37/71
Energy loss via synchrotron radiation emission
in a storage ring
In the time Tb spent in the bendings
the particle loses the energy U0
2 e 2  4
U 0   Pdt  PTb  P

c
3 0 
i.e. Energy Loss per turn (per electron)
e 2 4
E (GeV ) 4
U 0 (keV ) 
 88.46
3 0 
 (m)
Power radiated by a beam of average
current Ib: this power loss has to be
compensated by the RF system
I b  Trev
e
4
e
E (GeV ) 4 I ( A)
P(kW ) 
I b  88.46
3 0 
 (m)
N tot 
Power radiated by a beam of average
e 4
L(m) I ( A) E (GeV ) 4
P(kW ) 
LI b  14.08
current Ib in a dipole of length L
60  2
 (m) 2
(energy loss per second)
R. Bartolini, JUAS, 27-31 January 2014
38/71
The radiation integral (I)
The energy received by an observer (per unit solid angle at the source) is


d 2W
d 2P
 
dt  c 0  | RE (t ) |2 dt
d  d

Using the Fourier Transform we move to the frequency space

d 2W
 2c 0  | RE ( ) |2 d
d
0
Angular and frequency distribution of the energy received by an observer
2
d 3W
2 ˆ
 2 0cR E ( )
dd
Neglecting the velocity fields and assuming the observer in the far field:
n constant, R constant
2


3
2
dW
e
n  ( n   )   i ( t  n  r ( t ) / c )
Radiation Integral

e
dt
2
2

dd 40 4 c  (1  n   )

R. Bartolini, JUAS, 27-31 January 2014

39/71
The radiation integral (II)
The radiation integral can be simplified to [see Jackson]
d 3W
e 2 2

dd 40 4 2c

2
i ( t  n r ( t ) / c )
n

(
n


)
e
dt


How to solve it?
 determine the particle motion
r (t );  (t );  (t )
 compute the cross products and the phase factor
 integrate each component and take the vector square modulus
Calculations are generally quite lengthy: even for simple cases as for the
radiation emitted by an electron in a bending magnet they require Airy
integrals or the modified Bessel functions (available in MATLAB)
R. Bartolini, JUAS, 27-31 January 2014
40/71
Radiation integral for synchrotron radiation
Trajectory of the arc of circumference [see Jackson]
 

c 
c
r (t )    1  cos t , sin
t , 0 
 

 

In the limit of small angles we compute


 ct 
 ct 
    cos
 sin  
n  (n   )     || sin 
  
  


   ct 

 n  r (t ) 
 cos  
 t 
   t  sin 
c 

 c   

Substituting into the radiation integral and introducing
2


2 2 3/ 2

1


 
3
3c

d 3W
e 2  2 
 2 2
2 2 2
2
2



1



K
(

)

K
(

)
1/ 3
 2/3

2 2
d d 16 3 0 c  3c 2 
1






R. Bartolini, JUAS, 27-31 January 2014

41/71
Critical frequency and critical angle
2

d 3W
e 2  2 
 2 2
2 2 2
2
2



1



K
(

)

K
(

)
1/ 3
 2/3

2 2
d d 16 3 0 c  3c 2 
1







Using the properties of the modified Bessel function we observe that the
radiation intensity is negligible for  >> 1


2 2 3/ 2

1


   1
3c 3
Higher frequencies
have smaller critical
angle
3c 3

2
Critical frequency
c 
Critical angle
1  
c   c 

1/ 3
For frequencies much larger than the critical frequency and angles much
larger than the critical angle the synchrotron radiation emission is negligible
R. Bartolini, JUAS, 27-31 January 2014
42/71
Frequency distribution of radiated energy
Integrating on all angles we get the frequency distribution of the energy
radiated
dW
d 3I
3e 2 

d 

d 
d

d

4

c
C
0
4
dW
e 2   



d 40c  c 

K

 
/
5/ 3
( x)dx
C
1/ 3
 << c
1/ 2

dW
3 e

   e  / c
d
2 40 c  c 
2
 >> c
often expressed in terms of the
function S() with  = /c

9 3
S ( ) 
  K 5 / 3 ( x)dx
8 

 S ( )d  1
0

dW
3e 2 
2 e 2

K 5 / 3 ( x)dx 
S ( )

d 40c c  / c
9 0c
R. Bartolini, JUAS, 27-31 January 2014
43/71
Frequency distribution of radiated energy
It is possible to verify that the integral over the frequencies agrees with
the previous expression for the total power radiated [Hubner]



U 0 1 dW
1 2e 2
e 2c  4
P
 
d 
c   d  K 5 / 3 ( x)dx 
Tb Tb 0 d
Tb 9 0c 0
6 0c  2

The frequency integral extended up to the critical frequency contains half
of the total energy radiated, the peak occurs approximately at 0.3c
It is also convenient to define the critical
photon energy as
 c  c 
3 c 3

2 
For electrons, the critical energy in practical
units reads
50% 50%
E [GeV ]3
 c [keV ]  2.218
 0.665  E[GeV ]2  B [T ]
[m]
R. Bartolini, JUAS, 27-31 January 2014
44/71
Heuristic derivation of critical frequency
Synchrotron radiation is emitted in an arc of circumference with radius , Angle of
emission of radiation is 1/ (relativistic argument), therefore
T 

c
transit time in the arc of dipole
During this time the electron travels a distance

s  cT 

The time duration of the radiation pulse seen by the observer is the difference
between the time of emission of the photons and the time travelled by the electron
in the arc
s   1   1   2

  T 
   1 

c c    c  (1   ) 2c 3
The width of the Fourier
transform of the pulse is
 
R. Bartolini, JUAS, 27-31 January 2014
2c 3

45/71
Polarisation of synchrotron radiation
2

d 3W
e 2  2 
 2 2
2 2 2
2
2



1



K
(

)

K
(

)
1/ 3
 2/3

2 2
d d 16 3 0 c  3c 2 
1







Polarisation in
the orbit plane
Polarisation orthogonal
to the orbit plane
In the orbit plane  = 0, the polarisation is purely horizontal
Integrating on all frequencies we get the angular distribution of the energy
radiated

d 2W
d 3I
7e 2 5
1

d 
d 0 d d
640  (1   2 2 )5 / 2
 5  2 2 
1 
2 2
 7 1   
Integrating on all the angles we get a polarization on the plan of the orbit 7
times larger than on the plan perpendicular to the orbit
R. Bartolini, JUAS, 27-31 January 2014
46/71
Synchrotron radiation emission as a
function of beam the energy
Dependence of the frequency distribution of the energy radiated via synchrotron
emission on the electron beam energy
Critical frequency
3c 3
c 

2
Critical angle
No dependence on
the energy at longer
wavelengths
1  
c   c 

1/ 3
Critical energy
 c  c 
R. Bartolini, JUAS, 27-31 January 2014
3 c 3

2 
47/71
Brilliance with IDs (medium energy light sources)
Brilliance dependence
with current
with energy
with emittance
Medium energy storage rings with in-vacuum undulators operated at low gaps
(e.g. 5-7 mm) can reach 10 keV with a brilliance of 1020 ph/s/0.1%BW/mm2/mrad2
R. Bartolini, JUAS, 27-31 January 2014
48/71
Radiation from undulators and wigglers
Radiation emitted by undulators and wigglers
Types of undulators and wigglers
R. Bartolini, JUAS, 27-31 January 2014
49/71
Undulators and wigglers
Periodic array of magnetic
poles providing a sinusoidal
magnetic field on axis:
B  (0, B0 sin( ku z ), 0,)
Solution of equation of motions:


u K
u K 2
 ˆ
r (t )  
sin u t  xˆ    z ct 
cos(
2

t
)
u  z
2
2
16


Undulator
eB0 u
parameter
2mc
1  K2 

 z  1  2 1 
2 
2 
K
Constructive interference of radiation emitted at different poles
u
 u cos   n


u  K 2
n  2 1 
  2 2 
2 n 
2

d
R. Bartolini, JUAS, 27-31 January 2014
50/71
From the lecture on synchrotron radiation
Continuous spectrum characterized by c
= critical energy
bending magnet - a “sweeping searchlight”
c(keV) = 0.665 B(T)E2(GeV)
eg: for B = 1.4T E = 3GeV c = 8.4 keV
(bending magnet fields are usually lower
~ 1 – 1.5T)
wiggler - incoherent superposition K > 1
Max. angle of trajectory > 1/
Quasi-monochromatic spectrum with
peaks at lower energy than a wiggler
u  K 2  u
n 
1
 2
2 
2n 
2  n
undulator - coherent interference K < 1
Max. angle of trajectory < 1/
R. Bartolini, JUAS, 27-31 January 2014
nE[GeV ]2
 n (eV )  9.496
 K2 
u [m] 1 

2


51/71
Radiation integral for a linear undulator (I)
The angular and frequency distribution of the energy emitted by a wiggler is
computed again with the radiation integral:
dW
e

d  d 40 4 2c
3
2
2
2

i ( t  nˆ r / c )
ˆ
ˆ
n

(
n


)
e
dt


Using the periodicity of the trajectory we
can split the radiation integral into a sum
over Nu terms
dW
e

d d 40 4 2 c
3
2
2
2
u / 2  c
i ( N u 1)
i ( t  nˆ r / c )
i
i 2
ˆ
ˆ
n

(
n


)
e
dt
1

e

e

...

e

2
u / 2  c
where

2
res ( )
2c
res ( ) 
res ( )
R. Bartolini, JUAS, 27-31 January 2014

u  K 2
res ( )  2 1 
  2 2 
2 
2

52/71
Radiation integral for a linear undulator (II)
The radiation integral in an undulator or a wiggler can be written as
d 3W
e 2 2 N 2

d d
40c

 

 Fn ( K , ,  )
L N
 res ( ) 
 =  - n res()
The sum on  generates a series of sharp
peaks in the frequency spectrum harmonics
of the fundamental wavelength

  sin 2 ( N / res ( ))
  2 2
L N
 res ( )  N sin ( / res ( ))
The integral over one undulator period generates
a modulation term Fn which depends on the
angles of observations and K
Fn ( K , ,  ) 
0 / 2  c

nˆ  (nˆ   )e



0
2
i ( t  nˆ r / c )
dt
/2 c
R. Bartolini, JUAS, 27-31 January 2014
53/71
Radiation integral for a linear undulator (II)
e.g. on axis ( = 0,  = 0):
d 3W
e 2 2 N 2

d d
40c

n2 K 2 
Fn ( K ,0,0) 
J
(
Z
)

J
(
Z
)
 n1

n1
(1  K 2 / 2)  2
2

2

 
 Fn ( K ,0,0)
L N
 res (0) 
nK 2
Z
4(1  K 2 / 2)
Off-axis radiation
contains many
harmonics
Only odd harmonic are radiated on-axis;
as K increases the higher harmonic becomes stronger
R. Bartolini, JUAS, 27-31 January 2014
54/71
Angular patterns of the radiation emitted on
harmonics
Angular spectral flux as a function of frequency for a linear undulator; linear
polarisation solid, vertical polarisation dashed (K = 2)
u  K 2 

1  2 1 
2 
2 
Fundamental
wavelength emitted
by the undulator
2 
1
2
2nd harmonic, not
emitted on-axis !
R. Bartolini, JUAS, 27-31 January 2014
2 
1
3
3rd harmonic, emitted
on-axis !
55/71
Synchrotron radiation emission from a
bending magnet
Dependence of the frequency distribution of the energy radiated via synchrotron
emission on the electron beam energy
Critical frequency
3c 3
c 

2
Critical angle
No dependence on
the energy at longer
wavelengths
1  
c   c 

Critical energy
 c  c 
R. Bartolini, JUAS, 27-31 January 2014
1/ 3
3 c 3

2 
56/71
Undulators and wigglers (large K)
Radiated intensity emitted vs K
For large K the wiggler spectrum
becomes similar to the bending magnet
spectrum, 2Nu times larger.
Fixed B0, to reach the bending magnet
critical wavelength we need:
K
1
2
10
20
n
1
5 383 3015
R. Bartolini, JUAS, 27-31 January 2014
57/71
Undulator tuning curve (with K)
high K
K
eB0 u
2mc
Undulator
parameter
low K
K decreases by
opening the gap
of the undulator
(reducing B)
Brightness of a 5 m undulator 42 mm period with maximum K = 2.42 (ESRF)
Varying K one varies the wavelength emitted at various harmonics (not all
wavelengths of this graph are emitted at a single time)
R. Bartolini, JUAS, 27-31 January 2014
58/71
Spectral brightness of undulators of wiggler
Comparison of undulators for a 1.5 GeV ring for three harmonics (solid dashed
and dotted) compared with a wiggler and a bending magnet (ALS)
R. Bartolini, JUAS, 27-31 January 2014
59/71
Diamond undulators and wiggler
Spectral brightness for undulators and wigglers in state-of-the-art 3rd
generation light sources
R. Bartolini, JUAS, 27-31 January 2014
60/71
Summary of radiation characteristics of
undulators or wiggler
Undulators have weaker field or shorter periods (K< 1)
Produce narrow band radiation and harmonics
Intensity is proportional to Nu2
Wigglers have higher magnetic field (K >1)
Produce a broadband radiation
Intensity is proportional to Nu
R. Bartolini, JUAS, 27-31 January 2014
61/71
Type of undulators and wigglers
Electromagnetic undulators: the field is generated by current carrying coils;
they may have iron poles;
Permanent magnet undulators: the field is generated by permanent magnets
Samarium Cobalt (SmCo; 1T) and Neodymium Iron Boron (NdFeB; 1.4T); they
may have iron poles (hybrid undulators);
APPLE-II: permanent magnets arrays which can slide allowing the polarisation
of the magnetic field to be changed from linear to circular
In-vacuum: permanent magnets arrays which are located in-vacuum and
whose gap can be closed to very small values (< 5 mm gap!)
Superconducting wigglers: the field is generated by superconducting coils and
can reach very high peak fields (several T, 3.5 T at Diamond)
R. Bartolini, JUAS, 27-31 January 2014
62/71
Electromagnetic undulators (I)
HU64 at SOLEIL:
Period 64 mm
14 periods
variable polarisation electromagnetic
undulator
R. Bartolini, JUAS, 27-31 January 2014
Min gap 19 mm
Photon energy < 40 eV (1 keV with EM undulators)
63/71
Electromagnetic undulators (II)
Depending on the way the coil power supplies are powered it can generate
linear H, linear V or circular polarisations
R. Bartolini, JUAS, 27-31 January 2014
64/71
Permanent magnet undulators
Halback
configuration
hybrid
configuration with
steel poles
R. Bartolini, JUAS, 27-31 January 2014
65/71
In-vacuum undulators
U27 at Diamond
27 mm, 73 periods 7 mm gap, B
= 0.79 T; K = 2
R. Bartolini, JUAS, 27-31 January 2014
66/71
Apple-II type undulators (I)
HU64 at Diamond; 33 period of 64 mm; B = 0.96 T;
gap 15 mm; Kmax = 5.3
R. Bartolini, JUAS, 27-31 January 2014
67/71
Apple-II type undulators (II)
Four independent arrays of permanent magnets
Diagonally opposite arrays move longitudinal, all arrays move vertically
Sliding the arrays of magnetic pole it is possible to control the polarisation of the
radiation emitted
R. Bartolini, JUAS, 27-31 January 2014
68/71
Superconducting Wigglers
Superconducting wigglers
are used when a high
magnetic field is required
3 - 10 T
They need a cryogenic
system to keep the coil
superconductive
Nb3Sn and NbTi wires
SCMPW60 at Diamond
3.5 T coils cooled at 4 K
24 period of 64 mm
gap 10 mm
Undulator K = 21
R. Bartolini, JUAS, 27-31 January 2014
69/71
Summary and bibliography
Accelerated charged particles emit electromagnetic radiation
Synchrotron radiation is stronger for light particles and is emitted by
bending magnets in a narrow cone within a critical frequency
Undulators and wigglers enhance the synchrotron radiation emission
J. D. Jackson, Classical Electrodynamics, John Wiley & sons.
E. Wilson, An Introduction to Particle Accelerators, OUP, (2001)
M. Sands, SLAC-121, (1970)
R. P. Walker, CAS CERN 94-01 and CAS CERN 98-04
K. Hubner, CAS CERN 90-03
J. Schwinger, Phys. Rev. 75, pg. 1912, (1949)
B. M. Kincaid, Jour. Appl. Phys., 48, pp. 2684, (1977).
R.P. Walker: CAS 98-04, pg. 129
A. Ropert: CAS 98-04, pg. 91
P. Elleaume in Undulators, Wigglers and their applications, pg. 69-107
R. Bartolini, JUAS, 27-31 January 2014
70/71