Synchrotron Radiation

I was recently at CERN for one of their accelerator schools, learning about new and potentially disruptive plasma-based particle accelerator technology (disclaimer alert – the subject of my PhD). In honour of such a famous institution, I though I’d write about the worst enemy of CERN (and the friend of plasma accelerators) – synchrotron radiation.

Starting at the beginning then, a charged particle at rest is surrounded by an electric field. If the particle moves, the electric field must be dragged along with the particle. The electric field is therefore changing with time, and so generates a magnetic field (which is also generated by the motion of the charge – a current). This magnetic field is also changing with time, creating an electric field etc.

However, these changes in the fields aren’t communicated instantaneously but must move at the speed of light. A disturbance in one of the fields thus leads to ripples in both fields, streaking away from the particle at the speed of light. A charged particle in motion thus emits electromagnetic radiation. Delving into the mathematics more closely will show you that only acceleration of a particle will lead to electromagnetic radiation which radiates away to ‘infinity’.

Unfortunately this mathematics gets a little involved, so I’ll just state the result for the electric field of a particle moving along some arbitrary path. This is known as the Lienard-Wiechert Field, and is an absolute bugger to derive

\mathbf{E} = \frac{q}{4\pi\epsilon_0}\left(\frac{\mathbf{n} - \mathbf{\beta}}{\gamma^2R^2(1 - \mathbf{\beta}\cdot\mathbf{n})^3} + \frac{\mathbf{n} \times( (\mathbf{n} - \mathbf{\beta}) \times \dot{\mathbf{\beta}}))}{cR(1 - \mathbf{\beta}\cdot\mathbf{n})^3}\right)

The magnetic field has a similarly involved expression, so rather than writing it out I’ll point you to the Wikipedia page.

In this expression, \mathbf{\beta} is the velocity of the particle divided by the speed of light, \gamma is the relativistic factor of the particle, R is the distance from the particle to the observer and \mathbf{n} is the unit vector from the particle to the observer.

This is fairly complicated to look at, but there are a couple of interesting things we can pick out. The intensity of the radiation will scale like E^2, so the first term will drop off like R^{-4}. The total area of a sphere receiving this radiation at R only has area 4\pi R^2, so the total energy received goes like R^{-2}. This is then a non-radiative term – the particle is producing an electric field but this field isn’t propagating out to the observer. In fact, this term simply acts to move the field along with the particle.

The second term is more interesting, the intensity scales like R^{-2}, so the energy received by a sphere at radius R is constant – this is just electromagnetic radiation. From the \dot{\mathbf{\beta}} appearing we know it is nonzero only when the particle is accelerating.

We can integrate this expression over all observer angles and find the total power emitted by an electron travelling in a circle – this is the situation at CERN where charged particles travel in a giant ring of radius \rho:

P = \frac{q^2c\beta^4\gamma^4}{6\pi\epsilon_0\rho^2}

The killer for CERN is that \gamma^4 factor – if you double the energy of a particle, it radiates 16 times more energy. In the case of electrons, this very quickly becomes such a high power that the accelerator can’t accelerate any more – the electrons radiate away their energy quicker than it can be pumped in. For protons of the same energy, their \gamma^4 factor is 10^{13} times smaller, so radiate much less. This is why the LHC is so big, and uses protons rather than electrons. There are plans to build an even bigger LHC, unimaginatively called the Very Large Hadron Collider, or to accelerate electrons in a straight line so that they radiate much less energy at the proposed International Linear Collider. Both of these propositions are very expensive though, and very large, so there is much interest in smaller accelerators like plasma accelerators.

Anyway, that’s a lot of words and not much to look at, so how about some soothing animations? One thing which I glossed over above is the fact that these expressions for the electromagnetic fields need to be evaluated at the retarded time. In the words of my old undergraduate lecturer, lets do this in a retarded kind of way.

Obtain the retarded time by tracing back along the lightcone from the observation point.


In one dimension, assume we are observing a particle at time t_0 and observer position x_0. The particle emits radiation at the speed of light, so when we see radiation it has been emitted by the particle some time ago. We should take this into account when evaluating the appropriate electromagnetic fields, and evaluate at the retarded time. This is defined implicitly by the equation

t_r = t_0 - \frac{|x_0 - x(t_r)|}{c}

and encodes the fact that electromagnetic disturbances travel at the speed of light. In the general case this equation must be solved numerically if we want to use the Lienard-Wiechart potentials to make some pretty pictures. In effect we are travelling back from the observation point to the emission point along a lightcone. Once we intersect the particle path, we know the retarded time and the position/velocity/acceleration of the particle at the retarded time.

This sounds a bit tricky, but luckily for you I’ve had plenty of wet Saturday afternoons recently to do some simulations. First off, a particle going round in a circle with \gamma = 1.2:

The top left represents the Poynting flux S, or the radiation which would be observed. The retarded time is plotted beneath, and the two components of the electric field are plotted on the right. You can see how the retarded time is trailing the particle, and spreads out in a circle at the speed of light.

We can turn up the energy a bit, now \gamma = 5:

Here the electron is much closer to the speed of light, and is consequently chasing it’s own lightcone – notice how there is a much larger jump in retarded time (from yellow to green). This causes the Poynting flux to be squeezed into a smaller time window, it’s much more compressed, intense, and energetic. The temporal shortening of the radiation pulse corresponds to a broad bandwidth, and it is precisely this pulse of radiation that corresponds to broadband synchrotron radiation.

What about something less extreme, an electron travelling at constant velocity in a straight line, again at \gamma = 1.2:

We now see the action of the first term in the expression for the electric field – the field is being dragged along with the particle, but otherwise not radiating. Because the electron is travelling significantly slower than the speed of light, the field looks relatively undistorted. If you squint you’ll notice that the transverse field E_y is stronger than the longitudinal field, which is a manifestation of special relativity (Lorentz transforming a uniform field back into the lab frame where the charge is moving).

Back to something more interesting, a simple dipole motion at \gamma = 1.2:

We see what we hoped for, oscillating radiation perpendicular to the dipole, and zero far-field radiation parallel to the dipole. In the near-field (close to the dipole), there are electric field components everywhere, but it is clear that these must correspond to non-radiative terms and they die off rapidly.

Finally, let’s turn to a use for this radiation: undulators/wigglers. These are machines which purposefully oscillate an electron bunch to force synchrotron radiation emission, where the electron travels with high velocity in one direction and oscillates in the other:

You’ll see that unlike the dipole, the radiation field is asymmetric – the radiation in the direction of the electron has much higher frequency. In fact, the frequency is 4\gamma^2 times the electron oscillation frequency (relativistic Doppler shift), which means for energetic electrons we quickly reach the hard x-ray range. This is great for imaging purposes, we can oscillate electrons and produce hard, collimated, bright x-ray beams. This is another area where plasma accelerators help to shrink down the required technology and beamline length.

Bumping up to a \gamma of 5, my animations now don’t resolve the radiation wavelength, and the high-energy radiation is compressed into a short pulse ahead of the electron.


10 thoughts on “Synchrotron Radiation

  1. I just finished taking an undergraduate course in Electrodynamics (Griffiths!) and this is incredibly enlightening after staring at the equations in the textbook for so long, thanks! It would also be interesting to see the electric field plotted as a vector field on one plot, rather than componentwise.


    1. Thanks, glad you like the animations! I was hoping that solving for the retarded time explicitly might turn up some interesting very-near-field effects, but everything looks pretty smooth. I’ve been experimenting with plotting the field lines, I might come back to this in the future because they look pretty good.


    1. Sure, I used Matlab to calculate the fields over a discrete 2D grid from the instantaneous position/velocity/acceleration of a single charge moving along a prescribed path. Each frame the fields were plotted using the Matlab command ‘imagesc’ and colormap ‘parula’, and saved as a .png image. I then concatenated the images into some very large .gifs using software called ImageJ, though there are many ways to make animated GIFs. If you’re not a university student, Matlab is very expensive, so I’d suggest using e.g. Python and numpy if you wanted to try for yourself. And thanks!


  2. Hello!

    I’m trying to reproduce this cool plots in Mathematica, but i’m unable to do that.

    I’m using a trayectory of the form:

    Cos(t Sqrt[v] 2 pi)
    Sin( t Sqrt[v] 2 pi)

    Where v is a parameter for the velocity field that makes that the velocity remains between 0 and c. I’m solving everything following this cool blog post:

    where the author give a code for the Vector Field and the retarded time solver. But when i’m plotting the components of the vector i cannot get you nice plots.

    I’m using ArcTan as the scaling function for the plot, but when i plot S= E^2 (Poynting) i only get a blob, not those nice lines that appear in the plots.

    Can you provide some help? Or maybe the matlab code to contrast.

    Thank you very much!


    1. Hi, thanks for the link, and good work!. It seems you’re having trouble with getting the colouring of the plot correct. To get the colour scales above, what I actually plot is log(abs(E)). If you want to keep the distinction between positive and negative values, you could also colour the plot by sign(E) * log(abs(E)).


  3. Hi! The work you did is very very very very nice!!! I would like to use the illustrations above for a free book that i’m writing (for sure i will cite the source!). Do you authorize me to use the images (or provide me a sequence of illustration that you think would fit better in an A4 book paper size). It could also be nice to provide for free to the reader the MATLAB source code (only if not longer than a full A4 page).


Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s