U Arp1. 1. National Institute of Standards and Technology, Gaithersburg, MD 20899-8410.
Abstract
A method to include the influence of the vertical electron beam emittance onto the calculability of synchrotron radiation is introduced. It accounts for the finite vertical size and angular spread of the electron beam through a convolution procedure. The resulting angular spread of synchrotron radiation can differ significantly from the ideal Schwinger result, depending on the conditions. For the Synchrotron Ultraviolet Radiation Facility detailed results on the influence of the electron emittance for total power and polarization calculations are presented.
A method to include the influence of the vertical electron beam emittance onto the calculability of synchrotron radiation is introduced. It accounts for the finite vertical size and angular spread of the electron beam through a convolution procedure. The resulting angular spread of synchrotron radiation can differ significantly from the ideal Schwinger result, depending on the conditions. For the Synchrotron Ultraviolet Radiation Facility detailed results on the influence of the electron emittance for total power and polarization calculations are presented.
The calibration of light sources using storage rings is based on the calculability of synchrotron radiation [1-4]. The ultimate uncertainty of these calibrations is determined by the accuracy with which all necessary parameters to calculate the radiation output can be determined. Fig. 1 illustrates the problem. The question is, how much radiation passes through the aperture A, which is positioned at distance d from the source point S. If the emittance, i.e., the phase space area occupied by the electron beam [5], of the electron beam and diffraction effects are neglected, this radiation output is a function of electron energy E, bending radius , wavelength , and bandwidth Δ, electron beam current IB, distance d, width ΔX and height ΔY of the aperture, and the angle between the vertical center of the aperture and the orbital plane of the electrons .
Fig. 1
Geometry of the problem: An aperture A is positioned at distance d from the source point S. The aperture is assumed to be rectangular with horizontal (vertical) size ΔX(ΔY).
At the Synchrotron Ultraviolet Radiation Facility SURF III [6] the electron energy and orbital radius are determined by the magnetic flux density B and the radio-frequency RF. The Lorentz-force equation for the ideal orbit [7] gives the electron energy
and the synchronicity condition [7] delivers the orbital radius
e is the elementary charge, me the electron rest mass, and c the speed of light in vacuum. The electron beam current is determined using a well established electron counting procedure [8,9].The determination of the geometrical factors is often quite challenging as well. One has to determine the distance d between the defining aperture A and source point S with high accuracy.
2. Vertical Electron Beam Parameters in Equilibrium
For a weak-focusing storage ring like SURF, most equilibrium parameters can be calculated analytically. The vertical emittance is given by [7]
where Cq is Sands’ quantum excitation constant [10], the vertical betatron function, J the vertical damping partition number, and the orbital radius. For SURF, the vertical betatron function is determined by the orbital radius and the magnetic field index
[11] through
m. The equilibrium vertical Gaussian beam size is given by [7]
and the vertical Gaussian angular spread of the electron beam by [7]
Assuming E = 380 MeV, = 837.224 mm, and n = 0.594 the theoretical numbers for SURF are
In reality the vertical emittance is larger than this theoretical value, because of coupling between the horizontal and vertical motion. From the measured vertical beam size [12] the coupling is estimated to be of order 1 %.
3. Excitation of the Vertical Betatron Motion
The bunch volume at SURF is very small in its natural state, causing unacceptably short lifetimes of the electron beam. To extend the lifetime the vertical betatron oscillation is excited [13,14], causing the vertical angular spread and beam size to increase. If the vertical beam size is measured, the new emittance can be calculated, because the vertical betatron function depends only on the magnetic lattice [5] and will remain unchanged as the beam size changes. The following calculation depends on the assumption that the beam profiles remain Gaussian in shape and the beam is stable [14,15]. The new emittance is
and the new vertical angular spread of the electron beam
For example, if the vertical beam size is enlarged to 1 = 0.425 mm, which corresponds to a vertical full-width-at-half-maximum (FWHM) of 1 mm, the values are 1 = 1.660 × 105
μm μrad and
μrad.
4. Vertical Angular Distribution of Synchrotron Radiation
The vertical angular distribution of the emitted synchrotron radiation can be calculated using Schwinger’s equation [16]. We can separate the contributions into radiation with electrical vector parallel P║() and perpendicular P⊥() to the orbital plane of the electrons. is the wavelength of the emitted radiation and is the angle relative to the orbital plane of the electrons (see Fig. 1). The total power is Ptot() = P║() + P⊥().This formulation is only applicable to SURF and its unique geometry. e is the elementary charge, Δ the horizontal acceptance angle, the orbital radius, relativistic = E/(me
c2) and
, electron beam current IB, and bandwidth Δ. K2/3[ζ()] and K1/3[ζ()] are modified Bessel-functions of fractional order and
with the characteristic wavelength c = 4π /(33).
5. Influence of the Emittance on the Vertical Angular Spread of Synchrotron Radiation
Schwinger’s equation [16] is useful to calculate the synchrotron radiation emission of one electron in a perfect orbit, but does not take into account the electron emittance. This calculated vertical angular distribution of the synchrotron radiation has to be convolved with the vertical angular spread of the electron beam and also the angular spread caused by the finite vertical beam size, which depends on the distance from the tangent point.To simplify the procedure one can convolve the vertical angular spread of the electron beam with the contribution from the beam size first. Both are assumed to be Gaussian and the convolution of two Gaussians with widths 1 and 2 results in a Gaussian with total width
. For the total angular spread caused by the emittance of the electron beam this leads to
where d is the distance from the point of observation to the tangent point.Next, the total vertical angular spread of the emitted synchrotron radiation is the convolution of Ptot() with a Gaussian of width
. The convolution integral is [17]
This convolution integral has to be solved numerically. By applying the convolution theorem one can solve Eq. (12) easily. If X = FT(x) denotes the Fourier-transform of function x and x = IFT(X) its inverse Fourier-transform, the convolution can be written as
In Fig. 2 results are shown for SURF for d = 2500 mm, 1 = 1 mm, and = 100 nm, for both E = 380 MeV and E = 183 MeV.
Fig. 2
Ideal vertical distribution of the synchrotron radiation at = 100 nm, IB = 100 mA, Δ = 0.01 , Ptot() (●), P║() (■), and P⊥() (◆), as well as the calculated distribution taking into account emittance effects
,
, and
. The calculation was performed for a vertical Gaussian beam size y1 = 1 mm at distance d = 2500 mm. Top: E = 183 MeV, = 837.2217 mm, Bottom: E = 380 MeV, = 837.2242 mm. The vertical lines denote the integration limits for the later total flux and polarization calculations.
6. Optical Power Passing Through an Aperture
To calculate the power passing through the aperture A of vertical size ΔY and horizontal size ΔX, positioned at distance d from the source point S, the result of the convolution in Eq. (12) has to be integrated over the vertical angle. The integration over the horizontal acceptance angle is trivial, since the horizontal distribution is flat, and can be replaced by a multiplication [factor Δ = ΔX/d in Eqs. (9) and (10)]. To keep things simple we assume = 0 (vertical center of the aperture is in plane with the electron orbit).
Since the result of Eq. (13) was produced numerically, the integration in Eq. (14) has to be performed numerically as well. In Fig. 3 results are shown for different electron energies and vertical beam sizes. In this example for a 1 mm vertical beam size at 183 MeV and = 100 nm, the change relative to the ideal Schwinger value of the optical power is about 9 %.
Fig. 3
Ideal optical power
divided by the optical power including emittance effects
for d = 2500 mm, ΔX = ΔY = 10 mm. E = 380 MeV, 1 = 1 mm (○), E = 380 MeV, 1 = 0.425 mm (□), E = 183 MeV, 1 = 1 mm (∇), and E = 183 MeV, 1 = 0.425 mm (Δ).
7. Polarization of Radiation Passing Through an Aperture
The degree of linear polarization for synchrotron radiation is defined as [18]
If the degree of linear polarization of all the radiation passing through the aperture A is searched, Eq. (15) has to be integrated
To account for the emittance the two polarization contributions have to be convoluted the same way as the total, and then numerically integrated over the vertical acceptance angle
Again Fig. 4 shows clearly large deviations from the ideal Schwinger values for the polarization.
Fig. 4
Ideal degree of linear polarization of the radiation passing through the aperture A
divided by the degree of linear polarization including the vertical emittance
for d = 2500 mm, ΔX = ΔY = 10 mm. E = 380 MeV, 1 = 1 mm (●), E = 380 MeV, 1 = 0.425 mm (■), E = 183 MeV, 1 = 1 mm (▼), and E = 183 MeV, 1 = 0.425 mm (▲).
8. Conclusions
The influence of the vertical electron beam emittance on the vertical distribution of synchrotron radiation has been analyzed for SURF III.For the total power passing through an aperture and the polarization of the radiation it was found that deviations from the ideal Schwinger values can be of order several percent, depending on the actual conditions. The deviations are most prominent if the vertical acceptance angle Δ is of the same order as the FWHM of the vertical angular distribution, as illustrated in Fig. 2. If all radiation is accepted vertically there is no difference to expect. If the vertical distribution is much wider than the vertical acceptance, the differences are expected to be small.However, it is important to point out that this model only works well when the shape of the beam is Gaussian and no instabilities are present.