Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

10.3: Synchrotron Radiation

( \newcommand{\kernel}{\mathrm{null}\,}\)

Now let us consider a charged particle being accelerated in the direction perpendicular to its velocity u (for example by the magnetic component of the Lorentz force), so that its speed  u, and hence the magnitude  p of its momentum, do not change. In this case, the second term in the square brackets of Eq. (39) vanishes, and it yields

 P=Z0q26πm2c2(dpdτ)2=Z0q26πm2c2(dpdtret)2γ2.

Comparing this expression with Eq. (40), we see that for the same acceleration magnitude, the electromagnetic radiation is a factor of  γ2 larger. For modern accelerators, with  γ104105, such a factor creates an enormous difference. For example, if a particle is on a cyclotron orbit in a constant magnetic field (as was analyzed in Sec. 9.6), both u and \ p=γmu obey Eq. (9.150), so that

 |dpdtret|=ωcp=uRp=β2γmc2R,

(where  R is the orbit’s radius), so that for the power of this synchrotron radiation, Eq. (43) yields

 P=Z0q26πβ4γ4c2R2.Synchrotron radiation: total power

According to Eq. (9.153), at a fixed magnetic field (in particle accelerators, limited to a-few-tesla fields of the beam-bending magnets), the synchrotron orbit radius  R scales as  γ, so that according to Eq. (45),  P scales as  γ2, i.e. grows as the square of the particle’s energy  Eγ. For example, for typical parameters of the first electron cyclotrons (such as the General Electric’s machine in that the synchrotron radiation was first noticed in 1947),  R1 m,E0.3GeV(γ600), Eq. (45) gives a very modest electron energy loss per one revolution:  PTP(2πR/u)2πPR/c1keV. However, already by the mid-1970s, electron accelerators, with  R100 m, could give each particle energy  E10GeV, and the energy loss per revolution grew to  10MeV, becoming the major energy loss mechanism. For proton accelerators, such energy loss is much less of a problem, because  γ of an ultra-relativistic particle (at fixed  E) is proportional to  1/m, so that the estimates, at the same  R, should be scaled back by  (mp/me)41013. Nevertheless, in the giant modern accelerators such as the already mentioned LHC (with  R4.3 km and  E up to 7 TeV), the synchrotron radiation loss per revolution is rather noticeable  (PT6keV), leading not as much to particle deceleration as to a substantial photoelectron emission from the beam tube’s walls, creating harmful defocusing effects.

However, what is bad for particle accelerators and storage rings is good for the so-called synchrotron light sources – the electron accelerators designed especially for the generation of intensive synchrotron radiation – with the spectrum extending well beyond the visible light range. Let us analyze the angular and spectral distributions of such radiation. To calculate the angular distribution, let us select the coordinate axes as shown in Fig. 5, with the origin at the current location of the orbiting particle, the z-axis directed along its instant velocity (i.e. the vector  β), and the x-axis, toward the orbit’s center.

Screen Shot 2022-03-03 at 7.25.57 PM.pngFig. 10.5. The synchrotron radiation problem’s geometry.

In the general case, when the unit vector n toward the radiation’s observer is not within any of the coordinate planes, it has to be described by two angles – the polar angle  θ, and the azimuthal angle  φ between the x-axis and the projection 0P of the vector n to the plane  [x,y]. Since the length of the segment 0P is  sinθ, the Cartesian components of the relevant vectors are as follows:

 n={sinθcosφ,sinθsinφ,cosθ},β={0,0,β}, and ˙β={˙β,0,0}.

Plugging these expressions into the general Eq. (30), we get

 Synchrotron radiation: angular distributiondPdΩ=2Z0q2π2|˙β|2γ6f(θ,φ), where f(θ,φ)18γ6(1βcosθ)3[1sin2θcos2φγ2(1βcosθ)2],

According to this result, just as at the linear acceleration, in the ultra-relativistic limit, most radiation goes into a narrow cone (of a width  Δθγ1<<1) around the vector  β, i.e. around the instant direction of particle’s propagation. For such small angles, and  γ>>1,

 f(θ,φ)1(1+γ2θ2)3[14γ2θ2cos2φ(1+γ2θ2)2].

The left panel of Fig. 6 shows a color-coded contour map of this angular distribution  f(θ,φ), as observed on a distant plane normal to the particle’s instant velocity (in Fig. 5, parallel to the plane  [x,y]), while its right panel shows the factor  f as a function of  θ in two perpendicular directions: within the particle’s
rotation plane (in the direction parallel to the x-axis, i.e. at  φ=0) and perpendicular to this plane (along the y-axis, i.e. at  φ=±π/2). The result shows, first of all, that, in contrast to the case of linear acceleration, the narrow radiation cone is now not hollow: the intensity maximum is reached at  θ=0, i.e. exactly in the direction of the particle’s motion direction. Second, the radiation cone is not axially symmetric: within the particle rotation plane, the intensity drops faster (and even has nodes at  θ=±1/γ).

Screen Shot 2022-03-03 at 7.36.01 PM.png
Fig. 10.6. The angular distribution of the synchrotron radiation at  γ>>1.

As Fig. 5 shows, the angular distribution (47) of the synchrotron radiation was calculated for the (inertial) reference frame whose origin coincides with the particle’s position at this particular instant, i.e. its radiation pattern is time-independent in the frame moving with the particle. This pattern enables a semi-quantitative description of the radiation by an ultra-relativistic particle from the point of view of a stationary observer: if the observation point is on (or very close to) the rotation plane,11 it is being “struck” by the narrow radiation cone once each rotation period  T2πR/c, each “strike” giving a field pulse of a short duration  Δtret<<1/ωc - see Fig. 7.12

Screen Shot 2022-03-03 at 7.40.10 PM.png
Fig. 10.7. (a) The synchrotron radiation cones (at  γ>>1) for two close values of  tret, and (b) the in-plane component of the electric field, observed in the rotation plane, as a function of time  t – schematically.

The evaluation of the time duration  Δt of each pulse requires some care: its estimate  Δtret1/γωc is correct for the duration of the retarded time interval during that its cone is aimed at the observer. However, due to the time compression effect discussed in detail in Sec. 1 and described by Eq. (16), the pulse duration as seen by the observer is a factor of  1/(1β) shorter, so that

 Δt=(1β)Δtret1βγωc1γ3ωcγ3T, for γ>>1.

From the Fourier theorem, we can expect the frequency spectrum of such radiation to consist of numerous  (Nγ3>>1) harmonics of the particle rotation frequency  ωc, with comparable amplitudes. However, if the orbital frequency fluctuates even slightly  (δωc/ωc>1/N1/γ3), as it happens in most practical systems, the radiation pulses are not coherent, so that the average radiation power spectrum may be calculated as that of one pulse, multiplied by the number of pulses per second. In this case, the spectrum is continuous, extending from low frequencies all the way to approximately

 ωmax1/Δtγ3ωc.

In order to verify and quantify this result, let us calculate the spectrum of radiation, due to a single pulse. For that, we should first make the general notion of the radiation spectrum quantitative. Let us represent an arbitrary electric field (say that of the synchrotron radiation we are studying now), observed at a fixed point r, as a function of the observation time  t, as a Fourier integral:13

 E(t)=+Eωeiωtdt.

This expression may be plugged into the formula for the total energy of the radiation pulse (i.e. of the loss of particle’s energy  E) per unit solid angle:14

 dEdΩ+Sn(t)R2dt=R2Z0+|E(t)|2dt.

This substitution, followed by a natural change of the integration order, yield

 dEdΩ=R2Z0+ωωdω+ωωdωEωEω+dtei(ω+ω)t.

But the inner integral (over  t) is just  2πδ(ω+ω).15 This delta function kills one of the frequency integrals (say, one over  ω), and Eq. (53) gives us a result that may be recast as

 dEdΩ=+0I(ω)dω, with I(ω)4πR2Z0EωEω4πR2Z0EωEω,

where the evident frequency symmetry of the scalar product  EωEω has been utilized to fold the integral of  I(ω) to positive frequencies only. The first of Eqs. (54) makes the physical sense of the function  I(ω) very clear: this is the so-called spectral density of the electromagnetic radiation (per unit solid angle).16

To calculate the spectral density, we can express the function  Eω via  E(t) using the Fourier transform reciprocal to Eq. (51):

 Eω=12π+E(t)eiωtdt.

In the particular case of radiation by a single point charge, we may use here the second (radiative) term of Eq. (19):

 Eω=12πq4πε01cR+[n×{(nβ)×˙β}(1βn)3]reteiωtdt.

Since the vectors n and  β are more natural functions of the radiation (retarded) time  tret, let us use Eqs. (5) and (16) to exclude the observation time  t from this integral:

 Eω=q4πε012π1cR+[n×{(nβ)×˙β}(1βn)2]retexp{iω(tret+Rretc)}dtret.

Assuming that the observer is sufficiently far from the particle,17 we may treat the unit vector n as a constant and also use the approximation (8.19) to reduce Eq. (57) to

 Eω=12πq4πε01cRexp{iωrc}+[n×{(nβ)×˙β}(1βn)2exp{iω(tnrc)}]retdtret.

Plugging this expression into Eq. (54), and then using the definitions  c1/(ε0μ0)1/2 and  Z0(μ0/ε0)1/2, we get18

 I(ω)=Z0q216π3|+[n×{(nβ)×˙β}(1βn)2exp{iω(tnrc)}]retdtret|2.

This result may be further simplified by noticing that the fraction before the exponent may be represented as a full derivative over  tret,

 [n×{(nβ)×˙β}(1βn)2]ret[n×{(nβ)×dβ/dt}(1βn)2]retddt[n×(n×β)1βn]ret,

and working out the resulting integral by parts. At this operation, the time differentiation of the parentheses in the exponent gives  d[tretnr(tret)/c]/dtret=(1nu/c)ret(1βn)ret, leading to the cancellation of the remaining factor in the denominator and hence to a very simple general result:19

 Relativistic radiation: spectral densityI(ω)=Z0q2ω216π3|+[n×(n×β)exp{iω(tnrc)}]retdtret|2.

Now returning to the particular case of the synchrotron radiation, it is beneficial to choose the origin of time  tret so that at  tret=0, the angle  θ between the vectors n and  β takes its smallest value  θ0, i.e., in terms of Fig. 5, the vector n is within the plane  [y,z]. Fixing this direction of the axes, so that they do not move in further times, we can redraw that figure as shown in Fig. 8.

Screen Shot 2022-03-03 at 8.54.46 PM.png

Fig. 10.8. Deriving the synchrotron radiation’s spectral density. The vector n is static within the plane  [y,z], while the vectors  r(tret) and  βret  rotate, within the plane  [x,z], with the angular velocity  ωc of the particle.

In this “lab” reference frame, the vector n does not depend on time, while the vectors  r(tret) and  βret do depend on it via the angle  αωctret:,

 n={0,sinθ0,cosθ0},r(tret)={R(1cosα),0,Rsinα},βret{βsinα,0,βcosα}.

Now an easy multiplication yields

 [n×(n×β)]ret=β{sinα,sinθ0cosθ0cosα,sin2θ0sinα},

 [exp{iω(tnrc)}]ret=exp{iω(tretRccosθ0sinα)}.

As we already know, in the (most interesting) ultra-relativistic limit  γ>>1, most radiation is confined to short pulses, so that only small angles  αωcΔtretγ1 may contribute to the integral in Eq. (61). Moreover, since most radiation goes to small angles  θθ0γ1, it makes sense to consider only such small angles. Expanding both trigonometric functions of these small angles, participating in parentheses of Eq. (64), into the Taylor series, and keeping only the leading terms, we get

 tretRccosθ0sinαtretRcωctret+Rcθ202ωctret+Rcω3c6t3ret.

Since  (R/c)ωc=u/c=β1, in the two last terms we may approximate this parameter by 1. However, it is crucial to distinguish the difference of the two first terms, proportional to  (1β)tret, from zero; as we have done before, we may approximate it with  tret/2γ2. On the right-hand side of Eq. (63), which does not have such a critical difference, we may be bolder, taking20

 β{sinα,sinθ0cosθ0cosα,sin2θ0sinα}{α,θ0,0}{ωctret,θ0,0}.

As a result, Eq. (61) is reduced to

 I(ω)=Z0q216π3|axnx+ayny|2Z0q216π3(|ax|2+|ay|2),

where  ax and  ay are the following dimensionless factors:

Synchrotron radiation: spectral density

 axω+ωctretexp{iω2((θ20+γ2)tret+ω2c3t3ret)}dtret,ayω+θ0exp{iω2((θ20+γ2)tret+ω2c3t3ret)}dtret,

that describe the frequency spectra of two components of the synchrotron radiation, with mutually perpendicular directions of polarization. Defining the following dimensionless parameter

 νω3ωc(θ20+γ2)3/2,

which is proportional to the observation frequency, and changing the integration variable to  ξωctret/(θ20+γ2)1/2, the integrals (68) may be reduced to the modified Bessel functions of the second kind, but with fractional indices:

 ax=ωωc(θ20+γ2)+ξexp{32iν(ξ+ξ33)}dξ=23i(θ20+γ2)1/2νK2/3(ν)ay=ωωcθ0(θ20+γ2)1/2+exp{32iν(ξ+ξ33)}dξ=23θ0θ20+γ2νK1/3(ν)

Figure 9a shows the dependence of the Bessel factors, defining the amplitudes ax and ay, on the normalized observation frequency  ν. It shows that the radiation intensity changes with frequency relatively slowly (note the log-log scale of the plot!) until the normalized frequency, defined by Eq. (69), is increased beyond ~1. For most important observation angles  θ0γ this means that our estimate (50) is indeed correct, though formally the frequency spectrum extends to infinity.21

Screen Shot 2022-03-04 at 12.23.18 AM.png
Fig. 10.9. The frequency spectra of: (a) two components of the synchrotron radiation, at a fixed angle  θ0, and (b) its total (polarization- and angle-averaged) intensity.

Naturally, the spectral density integrated over the full solid angle exhibits a similar frequency behavior. Without performing the integration,22 let me just give the result (also valid for  γ>>1 only) for the reader’s reference:

 4πI(ω)dΩ=34πq2γζζK5/3(ξ)dξ, where ζ23ωωcγ3.

Figure 9b shows the dependence of this integral on the normalized frequency  ζ. (This plot is sometimes called the “universal flux curve”.) In accordance with the estimate (50), it reaches the maximum at

 ζmax0.3, i.e. ωmaxωc2γ3.

For example, in the National Synchrotron Light Source (NSLS-II) in the Brookhaven National Laboratory, near the SBU campus, with the ring circumference of 792 m, the electron revolution period  T is 2.64 μs. Calculating  ωc as  2π/T2.4×106 s1, for the achieved  γ6×103(E3GeV), we get  ωmax3×1017s1 (the photon energy  ωmax200eV), corresponding to soft X-rays. In the light of this estimate, the reader may be surprised by Fig. 10, which shows the calculated spectra of the radiation that this facility was designed to produce, with the intensity maxima at photon energies up to a few keV.

Screen Shot 2022-03-04 at 12.41.55 AM.png
Fig. 10.10. Design brightness of various synchrotron radiation sources of the NSLS-II facility. For bend magnets and wigglers, the “brightness” may be obtained by multiplication of the spectral density  I(ω) from one electron pulse, calculated above, by the number of electrons passing the source per second. (Note the non-SI units used by the synchrotron radiation community.) However, for undulators, there is an additional factor due to the partial coherence of radiation – see below. (Adapted from NSLS-II Source Properties and Floor Layout, was available online at https://www.bnl.gov/ps/docs/pdf/SourceProperties.pdf in 2011-2010.)

The reason for this discrepancy is that in the NLLS-II, and in all modern synchrotron light sources, most radiation is produced not by the circular orbit itself (which is, by the way, not exactly circular, but consists of a series of straight and bend-magnet sections), but by such bend sections, and the devices called wigglers and undulators: strings of several strong magnets with alternating field direction (Fig. 11), that induce periodic bending (wiggling”) of the electron’s trajectory, with the synchrotron radiation emitted at each bend.

Screen Shot 2022-03-04 at 12.43.52 AM.pngFig. 10.11. The generic structure of the wigglers, undulators, and free-electron lasers. (Adapted from http://www.xfel.eu/overview/how_does_it_work/. )

The difference between the wigglers and the undulators is more quantitative than qualitative: the former devices have a larger spatial period  λu (the distance between the adjacent magnets of the same polarity, see Fig. 11), giving enough space for the electron beam to bend by an angle larger than  γ1, i.e.
larger than the radiation cone’s width. As a result, the radiation arrives to an in-plane observer as a periodic sequence of individual pulses – see Fig. 12a.

Screen Shot 2022-03-04 at 12.46.53 AM.pngFig. 10.12. Waveforms of the radiation emitted by (a) a wiggler and (b) an undulator – schematically.

The shape of each pulse, and hence its frequency spectrum, are essentially similar to those discussed above,23 but with much higher local values of  ωc and hence  ωmax – see Fig. 10. Another difference is a much higher frequency of the pulses. Indeed, the fundamental Eq. (16) allows us to calculate the time distance between them, for the observer, as

 Δtttret Δtret (1β)λuu12γ2λuc<<λuc,

where the first two relations are valid at  λu<<R (the relation typically satisfied very well, see the numbers in Fig. 10), and the last two relations assume the ultra-relativistic limit. As a result, the radiation intensity, which is proportional to the number of the poles, is much higher than that from the bend magnets – see Fig. 10 again.

The situation is different in undulators – similar structures with a smaller spatial period  λu, in which the electron’s velocity vector oscillates with an angular amplitude smaller than  γ1. As a result, the radiation pulses overlap (Fig. 12b) and the radiation waveform is closer to the sinusoidal one. As a result, the radiation spectrum narrows to the central frequency24

 ω0=2πΔt2γ22πcλu.

For example, for the LSNL-II undulators with  λu=2 cm, this formula predicts a radiation peak at phonon energy  ω04keV, in a reasonable agreement with the results of quantitative calculations, shown in Fig. 10.25 Due to the spectrum narrowing, the intensity of undulator’s radiation is higher than that of wigglers using the same electron beam.

This spectrum-narrowing trend is brought to its logical conclusion in the so-called free-electron lasers26 whose basic structure is the same as that of wigglers and undulators (Fig. 11), but the radiation at each beam bend is so intense and narrow-focused that it affects the electron motion downstream the radiation cone. As a result, the radiation spectrum narrows around the central frequency (74), and its power grows as a square of the number N of electrons in the structure (rather than proportionately to N in wigglers and undulators).

Finally, note that wigglers, undulators, and free-electron lasers may be also used at the end of a linear electron accelerator (such as SLAC) which, as was noted above, may provide extremely high values of  γ, and hence radiation frequencies, due to the smallness of radiation energy losses at the electron acceleration stage. Very unfortunately, I do not have time/space to discuss (very interesting) physics of these devices in more detail.27


Reference

11 It is easy (and hence is left for the reader’s exercise) to show that if the observation point is much off-plane (say, is located on the particle orbit’s axis), the radiation is virtually monochromatic, with frequency  ωc. (As we know from Sec. 8.2, in the non-relativistic limit  u<<c this is true for any observation point.)

12 The fact that the in-plane component of each electric field’s pulse  E(t) is antisymmetric with respect to its central point, and hence vanishes at that point (as Fig. 7b shows), readily follows from Eq. (19).

13 In contrast to the single-frequency case (i.e. a monochromatic wave), we may avoid taking the real part of the complex function  (Eωeiωt) by requiring that in Eq. (51),  Eω=Eω. However, it is important to remember the factor 1⁄2 required for the transition to a monochromatic wave of frequency  ω0 and with real amplitude  E0:  Eω=E0[δ(ωω0)+δ(ω+ω0)]/2

14 Note that the expression under this integral differs from  dP/dΩ defined by Eq. (29) by the absence of the term  (1βn)=tret/t – see Eq. (16). This is natural because now we are calculating the wave energy arriving at the observation point r during the time interval  dt rather than  dtret.

15 See, e.g. MA Eq. (14.4).

16 The notion of spectral density may be readily generalized to random processes – see, e.g., SM Sec. 5.4.

17 According to the estimate (49), for a synchrotron radiation’s pulse, this restriction requires the observer to be much farther than  ΔrcΔtR/γ3 from the particle. With the values  R104 and  γ105 mentioned above,  Δr1011m, so this requirement is satisfied for any realistic radiation detector.

18 Note that for our current purposes of calculation of the spectral density of radiation by a single particle, the factor  exp{iωr/c} has got canceled. However, as we have seen in Chapter 8, this factor plays the central role at the interference of radiation from several (many) sources. Such interference is important, in particular, in undulators and free-electron lasers – the devices to be (qualitatively) discussed below.

19 Actually, this simplification is not occasional. According to Eq. (10b), the expression under the derivative in the last form of Eq. (60) is just the transverse component of the vector potential  A (give or take a constant factor), and from the discussion in Sec. 8.2 we know that this component determines the electric dipole radiation of the particle, which dominates the radiation field in our current case of a particle with a non-zero electric charge.

20 This expression confirms that the in-plane (x) component of the electric field is an odd function of  tret and hence of  tt0 (see its sketch in Fig. 7b), while the normal  (y) component is an even function of this difference. Also, note that for an observer exactly in the rotation plane  (θ0=0) the latter component equals zero for all times – the fact which could be predicted from the very beginning because of the evident mirror symmetry of the problem with respect to the particle rotation plane.

21 The law of the spectral density decrease at large  ν may be readily obtained from the second of Eqs. (2.158), which is valid even for any (even non-integer) Bessel function index  n:  axayν1/2exp{ν}. Here the exponential factor is certainly the most important one.

22 For that, and many other details, the interested reader may be referred, for example, to the fundamental review collection by E. Koch et al. (eds.) Handbook on Synchrotron Radiation (in 5 vols.), North-Holland, 1983-1991, or to a more concise monograph by A. Hofmann, The Physics of Synchrotron Radiation, Cambridge U. Press, 2007.

23 Indeed, the period  λu is typically a few centimeters (see the numbers in Fig. 10), i.e. is much larger than the interval  ΔrR/γ3 estimated above. Hence the synchrotron radiation results may be applied locally, to each electron beam’s bend. (In this context, a simple problem for the reader: use Eqs. (19) and (63) to explain the difference between shapes of the in-plane electric field pulses emitted at opposite magnetic poles of the wiggler, which is schematically shown in Fig. 12a.)

24 This important formula may be also derived in the following way. Due to the relativistic length contraction (9.20), the undulator structure period as perceived by beam electrons is  λ=λu/γ, so that the central frequency of the radiation in the reference frame moving with the electrons is  ω0=2πc/λ=2πcγ/λu. For the lab-frame observer, this frequency is Doppler-upshifted in accordance with Eq. (9.44):  ω0=ω0[(1+β)/(1β)]1/22γω0, giving the same result as Eq. (74).

25 Some of the difference is due to the fact that those plots show the spectral density of the number of photons  n=E/ω per second, which peaks at a frequency below that of the density of power, i.e. of the energy  E per second.

26 This name is somewhat misleading, because in contrast to the usual (“quantum”) lasers, a free-electron laser is essentially a classical device, and the dynamics of electrons in it is very similar to that in vacuum-tube microwave generators, such as the magnetrons briefly discussed in Sec. 9.6.

27 The interested reader may be referred, for example, to either P. Luchini and H. Motz, Undulators and Free-electron Lasers, Oxford U. Press, 1990; or E. Salin et al., The Physics of Free Electron Lasers, Springer, 2000.


This page titled 10.3: Synchrotron Radiation is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Konstantin K. Likharev via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?