#jsDisabledContent { display:none; } My Account |  Register |  Help

Article Id: WHEBN0000170808
Reproduction Date:

 Title: Synchrotron radiation Author: World Heritage Encyclopedia Language: English Subject: Collection: Publisher: World Heritage Encyclopedia Publication Date:

The electromagnetic radiation emitted when charged particles are accelerated radially (\mathbf{a}\perp \mathbf{v}) is called synchrotron radiation. It is produced, for example, in synchrotrons using bending magnets, undulators and/or wigglers. It is similar to cyclotron radiation except that synchrotron radiation is generated by the acceleration of ultrarelativistic charged particles through magnetic fields. Synchrotron radiation may be achieved artificially in synchrotrons or storage rings, or naturally by fast electrons moving through magnetic fields. The radiation produced in this way has a characteristic polarization and the frequencies generated can range over the entire electromagnetic spectrum.

Synchrotron radiation from a bending magnet

## Contents

• History 1
• Properties of synchrotron radiation 2
• Emission mechanism 3
• Synchrotron radiation from accelerators 4
• Synchrotron radiation in astronomy 5
• History of detection 5.1
• Pulsar wind nebulae 5.2
• Formulation 6
• Liénard–Wiechert Field 6.1
• Velocity ⊥ acceleration: synchrotron radiation 6.2
• Example 1: bending magnet 6.4
• Integrating 6.4.1
• Frequency distribution of radiated energy 6.4.2
• Synchrotron radiation emission as a function of the beam energy 6.4.3
• Polarization of synchrotron radiation 6.4.4
• Example 2: undulator 6.5
• Solution of equation of motion and undulator equation 6.5.1
• Radiation from the undulator 6.5.2
• Notes 8
• References 9

## History

General Electric synchrotron accelerator built in 1946, the origin of the discovery of synchrotron radiation. The arrow indicates the evidence of radiation.

Synchrotron radiation was named after its discovery in a General Electric synchrotron accelerator built in 1946 and announced in May 1947 by Frank Elder, Anatole Gurewitsch, Robert Langmuir, and Herb Pollock in a letter entitled "Radiation from Electrons in a Synchrotron".[1] Pollock recounts:

"On April 24, Langmuir and I were running the machine and as usual were trying to push the electron gun and its associated pulse transformer to the limit. Some intermittent sparking had occurred and we asked the technician to observe with a mirror around the protective concrete wall. He immediately signaled to turn off the synchrotron as "he saw an arc in the tube." The vacuum was still excellent, so Langmuir and I came to the end of the wall and observed. At first we thought it might be due to Cherenkov radiation, but it soon became clearer that we were seeing Ivanenko and Pomeranchuk radiation."[2]

1. Broad Spectrum (which covers from microwaves to hard X-rays): the users can select the wavelength required for their experiment.
2. High Flux: high intensity photon beam allows rapid experiments or use of weakly scattering crystals.
3. High Brilliance: highly collimated photon beam generated by a small divergence and small size source (spatial coherence)
4. High Stability: submicron source stability
5. Polarization: both linear and circular
6. Pulsed Time Structure: pulsed length down to tens of picoseconds allows the resolution of process on the same time scale.

## Emission mechanism

When high-energy particles are in rapid motion, including electrons forced to travel in a curved path by a magnetic field, synchrotron radiation is produced. This is similar to a radio antenna, but with the difference that, in theory, the relativistic speed will change the observed frequency due to the Doppler effect by the Lorentz factor, \gamma. Relativistic length contraction then bumps the frequency observed in the lab by another factor of \gamma, thus multiplying the GHz frequency of the resonant cavity that accelerates the electrons into the X-ray range. The radiated power is given by the relativistic Larmor formula while the force on the emitting electron is given by the Abraham–Lorentz–Dirac force. The radiation pattern can be distorted from an isotropic dipole pattern into an extremely forward-pointing cone of radiation. Synchrotron radiation is the brightest artificial source of X-rays. The planar acceleration geometry appears to make the radiation linearly polarized when observed in the orbital plane, and circularly polarized when observed at a small angle to that plane. Amplitude and frequency are however focussed to the polar ecliptic.

Synchrotron radiation may occur in accelerators either as a nuisance, causing undesired energy loss in particle physics contexts, or as a deliberately produced radiation source for numerous laboratory applications. Electrons are accelerated to high speeds in several stages to achieve a final energy that is typically in the GeV range. In the LHC proton bunches also produce the radiation at increasing amplitude and frequency as they accelerate with respect to the vacuum field, propagating photoelectrons, which in turn propagate secondary electrons from the pipe walls with increasing frequency and density up to 7x1010. Each proton may lose 6.7keV per turn due to this phenomenon.[3]

Messier 87's Energetic Jet, HST image. The blue light from the jet emerging from the bright AGN core, towards the lower right, is due to synchrotron radiation.

Synchrotron radiation is also generated by astronomical objects, typically where relativistic electrons spiral (and hence change velocity) through magnetic fields. Two of its characteristics include non-thermal power-law spectra, and polarization.[4]

### History of detection

It was first detected in a jet emitted by Messier 87 in 1956 by Geoffrey R. Burbidge,[5] who saw it as confirmation of a prediction by Iosif S. Shklovsky in 1953, but it had been predicted earlier by Hannes Alfvén and Nicolai Herlofson[6] in 1950.

T. K. Breus noted that questions of priority on the history of astrophysical synchrotron radiation is complicated, writing:

"In particular, the Russian physicist V.L. Ginzburg broke his relationships with I.S. Shklovsky and did not speak with him for 18 years. In the West, Thomas Gold and Sir Fred Hoyle were in dispute with H. Alfven and N. Herlofson, while K.O. Kiepenheuer and G. Hutchinson were ignored by them."[7]
Crab Nebula. The bluish glow from the central region of the nebula is due to synchrotron radiation.

Supermassive black holes have been suggested for producing synchrotron radiation, by ejection of jets produced by gravitationally accelerating ions through the super contorted 'tubular' polar areas of magnetic fields. Such jets, the nearest being in Messier 87, have been confirmed by the Hubble telescope as apparently superluminal, travelling at 6×c (six times the speed of light) from our planetary frame. This phenomenon is caused because the jets are travelling very near the speed of light and at a very small angle towards the observer. Because at every point of their path the high-velocity jets are emitting light, the light they emit does not approach the observer much more quickly than the jet itself. Light emitted over hundreds of years of travel thus arrives at the observer over a much smaller time period (ten or twenty years) giving the illusion of faster than light travel. There is no violation of special relativity.[8]

### Pulsar wind nebulae

A class of astronomical sources where synchrotron emission is important is the pulsar wind nebulae, a.k.a. plerions, of which the Crab nebula and its associated pulsar are archetypal. Pulsed emission gamma-ray radiation from the Crab has recently been observed up to ≥25 GeV,[9] probably due to synchrotron emission by electrons trapped in the strong magnetic field around the pulsar. Polarization in the Crab[10] at energies from 0.1 to 1.0 MeV illustrates a typical synchrotron radiation.

## Formulation

### Liénard–Wiechert Field

where

\mathbf R(t')=\mathbf r-\mathbf r_0(t'),
R(t')=|\mathbf R|,
\hat{\mathbf n}(t')=\mathbf R/R,

which is the unit vector between the observation point and the position of the charge at the retarded time, and t' is the retarded time.

In equation (1), and (2), the first terms fall off as the inverse square of the distance from the particle, and this first term is called the generalized Coulomb field or velocity field. And the second terms fall off as the inverse first power of the distance from the source, and it is called the radiation field or acceleration field. If we ignore the velocity field, the radial component of Poynting's Vector resulted from the Liénard–Wiechert field can be calculated to be

Note that

• The spatial relationship between \vec{\beta} and \dot{\vec{\beta}} determines the detailed angular power distribution.
• The relativistic effect of transforming from the rest frame of the particle to the observer's frame manifests itself by the presence of the factors (1-\vec{\beta}\mathbf{\cdot}\vec{\mathbf{n}}) in the denominator of Eq. (3).
• For ultrarelativistic particles the latter effect dominates the whole angular distribution.

The energy radiated into per solid angle during a finite period of acceleration from t'=T_1 to t'=T_2 is

\frac{\mathrm{d}P}{\mathrm{d}\mathit{\Omega}} = R(t')^2\,[\mathbf{S}(t')\mathbf{\cdot}\hat{\mathbf{n}}(t')]\,\frac{\mathrm{d}t}{\mathrm{d}t'} = R(t')^2\,\mathbf{S}(t')\mathbf{\cdot}\hat{\mathbf{n}}(t')\,[1-\vec{\beta}(t')\mathbf{\cdot}\hat{\mathbf{n}}(t')]

Integrating Eq. (4) over the all solid angles, we get relativistic generalization of Larmor's formula

P=\frac{e^2}{6\pi \varepsilon _0 c}\gamma ^6 \left [ \left | \dot{\vec{\beta }} \right |^2 -\left | \vec{\beta}\times \dot{\vec{\beta }}\right |^2 \right ]\qquad (5)

However, this also can be derived by relativistic transformation of the 4-acceleration in Larmor's formula.

### Velocity ⊥ acceleration: synchrotron radiation

When the electron velocity approaches the speed of light, the emission pattern is sharply collimated forward.

When the charge is in instantaneous circular motion, its acceleration \dot{\vec{\beta}} is perpendicular to its velocity \vec{\beta}. Choosing a coordinate system such that instantaneously \vec{\beta} is in the z direction and \dot{\vec{\beta}} is in the x direction, with the polar and azimuth angles \theta and \phi defining the direction of observation, the general formula Eq. (4) reduces to

\frac{\mathrm{d}P}{\mathrm{d}\mathit{\Omega}} = \frac{q^2}{16\pi^2\epsilon_0 c}\frac{|\dot{\vec{\beta }}|^2}{(1-\beta\cos\theta)^3}\left[1-\frac{\sin^2\theta\cos^2\phi}{\gamma^2(1-\beta\cos\theta)^2}\right]. \qquad (6)

In the relativistic limit (\gamma\gg 1), the angular distribution can be written approximately as

The factors (1-\beta\cos\theta) in the denominators tip the angular distribution forward into a narrow cone like the beam of a headlight pointing ahead of the particle. A plot of the angular distribution (dP/dΩ vs. γθ) shows a sharp peak around θ=0.

Integration over the whole solid angle yields the total power radiated by one electron

P=\frac{e^2}{6\pi\epsilon_0c}\left | \dot{\vec{\beta }} \right |^2\gamma ^4=\frac{e^2c}{6\pi\epsilon_0}\frac{\beta^4\gamma ^4}{\rho ^2}=\frac{e^4}{6\pi\epsilon_0m^4c^5}E^2B^2,\qquad (8)

where E is the electron energy, B is the magnetic field, and ρ is the radius of curvature of the track in the field. Note that the radiated power is proportional to 1/m^4, 1/\rho^2, and B^2. In some cases the surfaces of vacuum chambers hit by synchrotron radiation have to be cooled because of the high power of the radiation.

Using B=\frac{E\beta }{erSin(\alpha)} where \alpha is the angle between the velocity and the magnetic field and r is the radium of the circular acceleration. The power emitted is: P=\frac{e^2}{6\pi\epsilon_0m^4c^5r^2Sin^2(\alpha)}E^4\beta^2 Thus the power emitted increases with the energy to the fourth and decrease with the radius to the square and the mass to the fourth. This radiation is what limit the energy of an electron-positron circular collider. For a proton-antiproton collider, the energy is limited by the maximum magnetic field, this is why the LHC has an energy only 70 times higher than the LEP even though the mass of the proton being 2000 times higher than the mass of the electron, the energy could be 2000 times higher to have the same synchroton radiation.

The energy received by an observer (per unit solid angle at the source) is

\frac{d^2W}{d\Omega }=\int_{-\infty }^{\infty }\frac{d^2P}{d\Omega }dt=c\varepsilon _0\int_{-\infty }^{\infty }\left | R\vec{E}(t) \right |^2dt

Using the Fourier Transformation we move to the frequency space \frac{d^2W}{d\Omega }=2c\varepsilon _0\int_{0 }^{\infty }\left | R\vec{E}(\omega) \right |^2d\omega

Angular and frequency distribution of the energy received by an observer (consider only the radiation field)

\frac{d^3W}{d\Omega d\omega }=2c\varepsilon _0R^2\left | \vec{E}(\omega) \right |^2=\frac{e^2}{4\pi\varepsilon_0 4\pi^2 c}\left | \int_{-\infty}^{\infty}\frac{\hat{n}\times\left [ \left ( \hat{n}-\vec{\beta } \right )\times\dot{\vec{\beta }} \right ]}{\left ( 1-\hat{n}\cdot \vec{\beta } \right )^2}e^{i\omega(t-\hat{n}\cdot\vec{r}(t)/c)}dt\right |^2 \qquad (9)

Therefore, if we know the particle's motion, cross products term, and phase factor, we could calculate the radiation integral. However, calculations are generally quite lengthy (even for simple cases as for the radiation emitted by an electron in a bending magnet, they require Airy function or the modified Bessel functions).

### Example 1: bending magnet

#### Integrating

Trajectory of the arc of circumference

Trajectory of the arc of circumference is

\vec{r}(t)=\left ( \rho \sin\frac{\beta c}{\rho}t, \rho\left ( 1-\cos\frac{\beta c}{\rho}t \right ), 0 \right)

In the limit of small angles we compute

\hat{n}\times\left ( \hat{n}\times\vec{\beta} \right ) =\beta\left [ -\vec{\varepsilon }_\parallel \sin\left ( \frac{\beta c t}{\rho} \right )+\vec{\varepsilon}_\perp \cos\left ( \frac{\beta c t}{\rho}\right )\sin\theta \right ]

\omega\left ( t-\frac{\hat{n}\cdot \vec{r}(t)}{c} \right ) = \omega\left [ t-\frac{\rho}{c}\sin\left ( \frac{\beta c t}{\rho} \right )\cos\theta \right ]

Substituting into the radiation integral and introducing \xi =\frac{\rho \omega}{3c\gamma^3}\left ( 1+\gamma^2 \theta^2 \right )^{3/2}

\frac{d^3 W}{d\Omega d\omega}=\frac{e^2}{16\pi^3\varepsilon_0 c} \left ( \frac{2\omega\rho}{3c\gamma^2} \right )^2 \left ( 1+\gamma^2 \theta^2 \right )^2 \left [ K_{2/3}^2(\xi ) + \frac{\gamma^2 \theta^2}{1+\gamma^2 \theta^2}K_{1/3}^2(\xi)\right ]\qquad (10)

, where the function K is a modified Bessel function of the second kind.

#### Frequency distribution of radiated energy

From Eq.(10), we observe that the radiation intensity is negligible for \xi\gg 1. Critical frequency is defined as the frequency when \xi=\frac{1}{2} and \theta=0. So,

\omega_\text{c}=\frac{3}{2}\frac{c}{\rho}\gamma^3,

and critical angle is defined as the angle for which \xi(\theta_\text{c})\simeq\xi(0)+1 and is approximately

\theta_\text{c}\simeq\frac{1}{\gamma}\left ( \frac{2\omega_\text{c}}{\omega} \right )^{1/3} [11]

For frequencies much larger than the critical frequency and angles much larger than the critical angle, the synchrotron radiation emission is negligible.

Integrating on all angles, we get the frequency distribution of the energy radiated.

\frac{dW}{d\omega}=\oint \frac{d^3 W}{d\omega d\Omega }d\Omega =\frac{\sqrt{3}e^2}{4\pi\varepsilon_0 c}\gamma\frac{\omega}{\omega_\text{c}}\int_{\omega/\omega_\text{c}}^{\infty}K_{5/3}(x)dx

If we define

S(y)\equiv \frac{9\sqrt{3}}{8\pi}y\int_{y}^{\infty}K_{5/3}(x)dx \int_{0}^{\infty}S(y)dy=1

, where y=\frac{\omega}{\omega_\text{c}} . Then,

Note that \frac{dW}{d\omega}\sim \frac{e^2}{4\pi \varepsilon_0c}\left ( \frac{\omega \rho}{c} \right )^{1/3}, if \omega\ll \omega_\text{c}, and \frac{dW}{d\omega}\approx \sqrt{\frac{3\pi}{2}}\frac{e^2}{4\pi\varepsilon_0 c}\gamma\left ( \frac{\omega}{\omega_\text{c}} \right )^{0.5} e^{-\omega/\omega_\text{c}}, if \omega\gg \omega_\text{c}

The formula for spectral distribution of synchrotron radiation, given above, can be expressed in terms of a rapidly converging integral with no special functions involved [12] (see also modified Bessel functions ) by means of the relation:

\int_{\xi}^\infty K_{5/3} (x) dx = \frac{1}{ \sqrt{3}} \, \int_0^\infty \, \frac{9+36x^2+16x^4}{(3+4x^2) \sqrt{1+x^2/3}} \exp \left[- \xi \left(1+\frac{4x^2}{3}\right) \sqrt{1+\frac{x^2}{3}} \right] \ dx

#### Synchrotron radiation emission as a function of the beam energy

Relationship between power radiated and the photon energy

First, define the critical photon energy as \varepsilon_c=\hbar \omega_\text{c}=\frac{3}{2}\frac{\hbar c}{\rho}\gamma^3

Then, the relationship between radiated power and photon energy is shown in the graph on the right side. The higher the critical energy, the more photons with high energies are generated. Note that, there is no dependence on the energy at longer wavelength.

In Eq.(10), the first term K_{2/3}^2(\xi) is the radiation power with polarization in the orbit plane, and the second term \frac{\gamma^2 \theta^2}{1+\gamma^2 \theta^2}K_{1/3}^2(\xi) is the polarization orthogonal to the orbit plane. In the orbit plane \theta=0, the polarization is purely horizontal. Integrating on all frequencies, we get the angular distribution of the energy radiated

\frac{d^2 W}{d\Omega }=\int_{0}^{\infty}\frac{d^3W}{d\omega d\Omega }d\omega =\frac{7e^2 \gamma^5}{64\pi\varepsilon_0\rho}\frac{1}{(1+\gamma^2\theta^2)^{5/2}}\left [1+\frac{5}{7}\frac{\gamma^2\theta^2}{1+\gamma^2\theta^2} \right ] \qquad (12)

Integrating on all the angles, we find that seven times as much energy is radiated with parallel polarization as with perpendicular polarization. The radiation from a relativistically moving charge is very strongly, but not completely, polarized in the plane of motion.

### Example 2: undulator

#### Solution of equation of motion and undulator equation

An undulator consists of a periodic array of magnets, so that they provide a sinusoidal magnetic field.

\vec{B}=\left ( 0, B_0 \sin(k_\text{u} z),0 \right )
undulator

Solution of equation of motion is

\vec{r}(t)=\frac{\lambda_\text{u} K}{2\pi\gamma}\sin \omega_\text{u}t\cdot \hat{x} +\left ( \bar{\beta_z}ct+\frac{\lambda_\text{u}K^2}{16\pi\gamma^2}\cos(2\omega_\text{u}t) \right )\cdot \hat{z}

where, K=\frac{eB_0\lambda_\text{u}}{2\pi mc} , and \bar{\beta_z}=1-\frac{1}{2\gamma^2}\left ( 1+\frac{K^2}{2} \right )

, and the parameter K is called the undulator parameter.

Constructive interference of the beam in the undulator

Condition for the constructive interference of radiation emitted at different poles is

d=\frac{\lambda_\text{u}}{\bar{\beta}}-\lambda_\text{u}\cos\theta=n\lambda

Expanding \cos\theta \approx 1- \frac{\theta^2}{2} and neglecting the terms O(\theta^2) in the resulting equation, one obtains

\lambda_n = \frac{\lambda_\text{u}}{2\gamma^2n}\left ( \frac{1+\frac{K^2}{2}+\gamma^2\theta^2}{\bar{\beta}} \right )

For \bar{\beta} \rightarrow 1 one finally gets

\lambda_n = \frac{\lambda_\text{u}}{2\gamma^2n}\left ( 1+\frac{K^2}{2}+\gamma^2\theta^2 \right )\qquad (13)

This equation is called the undulator equation.

\frac{d^3W}{d\Omega d\omega}=\frac{e^2}{4\pi\varepsilon_0 4\pi^2 c}\left|\int_{-\infty}^{\infty}\frac{\hat{n}\times\left[\left(\hat{n}-\vec{\beta}\right)\times\dot{\vec{\beta}}\right]}{\left(1-\hat{n}\cdot\vec{\beta}\right)^2}e^{i\omega(t-\hat{n}\cdot\vec{r}(t)/c)}dt\right|^2

Using the periodicity of the trajectory, we can split the radiation integral into a sum over N_u terms.

\frac{d^3W}{d\Omega d\omega }=\frac{e^2\omega^2}{4\pi\varepsilon_0 4\pi^2 c}\left | \int_{-\lambda_u/2\bar{\beta}c}^{\lambda_u/2\bar{\beta}c}\hat{n}\times\left ( \hat{n}\times\vec{\beta} \right )e^{i\omega(t-\hat{n}\cdot\vec{r}(t)/c)}dt\right|^2 \left|1+e^{i\delta}+e^{2i\delta}+\cdots +e^{i(N_u-1)\delta} \right |^2 \qquad (14)

, where 　\bar{\beta}=\beta\left ( 1-\frac{K^2}{4\gamma^2} \right )

Peak frequencies become sharp as the number N increases

, and 　　 \delta=\frac{2\pi \omega}{\omega_\text{res} (\theta)},　　　\omega_\text{res}(\theta)=\frac{2\pi c}{\lambda_\text{res}(\theta)},　　and　　　\lambda_\text{res}(\theta)=\frac{\lambda_u}{2\gamma^2}\left ( 1+\frac{K^2}{2}+\gamma^2\theta^2 \right )

Only odd harmonics are radiated on-axis

The radiation integral in an undulator can be written as

\frac{d^3 W}{d\Omega d\omega}=\frac{e^2\gamma^2N^2}{4\pi\varepsilon_0 c} L\left ( N\frac{\Delta \omega_n}{\omega_\text{res}(\theta)} \right ) F_n (K, \theta, \phi) \qquad (15)

The sum of \delta generates a series of sharp peaks in the frequency spectrum harmonics of fundamental wavelength

L\left ( N\frac{\Delta \omega_k}{\omega_\text{res}(\theta)} \right ) =\frac{\sin^2\left ( N\pi\Delta \omega_k / \omega_\text{res}(\theta) \right )}{N^2 \sin^2\left ( \pi\Delta \omega_k/\omega_\text{res}(\theta) \right )}

, and F_n depends on the angles of observations and K

F_n(K,\theta,\phi)\propto \left | \int_{-\lambda_u/2\bar{\beta} c}^{\lambda_u/2\bar{\beta} c}\hat{n}\times\left ( \hat{n}\times\vec{\beta} \right )e^{i\omega(t-\hat{n}\cdot\vec{r}(t)/c)}dt\right|^2

On the axis(\theta=0, \phi=0), the radiation integral becomes

\frac{d^3 W}{d\Omega d\omega}=\frac{e^2\gamma^2N^2}{4\pi\varepsilon_0 c}L\left ( N\frac{\Delta \omega_n}{\omega_\text{res}(0)} \right )F_n(K,0,0)

and,

F_n(K,0,0)=\frac{n^2K^2}{1+K^2/2} \left[J_{\frac{n+1}{2}}(Z)-J_{\frac{n-1}{2}}(Z) \right ]^2

, where Z=\frac{nK^2}{4(1+K^2/2)}

Note that only odd harmonics are radiated on-axis, and as K increases higher harmonic becomes stronger.

## Notes

1. ^ Elder, F. R.; Gurewitsch, A. M.; Langmuir, R. V.; Pollock, H. C., "Radiation from Electrons in a Synchrotron" (1947) Physical Review, vol. 71, Issue 11, pp. 829-830
2. ^ Iwanenko D., Pomeranchuk I., "On the maximal energy attainable in betatron", (1944) Physical Review vol.65, p. 343
3. ^ [2] Synchrotron Radiation Damping in the LHC 2005 Joachim Tuckmantel.
4. ^ Vladimir A. Bordovitsyn, "Synchrotron Radiation in Astrophysics" (1999) Synchrotron Radiation Theory and Its Development, ISBN 981-02-3156-3
5. ^ Burbidge, G. R. "On Synchrotron Radiation from Messier 87. Astrophysical Journal, vol. 124, p. 416"
6. ^ Alfvén, H.; Herlofson, N. "Cosmic Radiation and Radio Stars" Physical Review (1950), vol. 78, Issue 5, pp. 616–616
7. ^ Breus, T. K., "Istoriya prioritetov sinkhrotronnoj kontseptsii v astronomii %t (Historical problems of the priority questions of the synchrotron concept in astrophysics)" (2001) in Istoriko-Astronomicheskie Issledovaniya, Vyp. 26, p. 88 - 97, 262 (2001)
8. ^ Chase, Scott I. "Apparent Superluminal Velocity of Galaxies". Retrieved 22 August 2012.
9. ^ "Observation of Pulsed {gamma}-Rays Above 25 GeV from the Crab Pulsar with MAGIC", Science 21 November 2008: Vol. 322. no. 5905, pp. 1221–1224"
10. ^ Dean et al.,"Polarized Gamma-Ray Emission from the Crab", Science 29 August 2008: Vol. 321. no. 5893, pp. 1183–1185
11. ^ Jackson, John David (1999). Classical Electrodynamics (3rd ed.). Chichester: Wiley. p. 680.
12. ^ M.Kh.Khokonov. Cascade Processes of Energy Loss by Emission of Hard Photons // JETP, V.99, No.4, pp. 690-707 \ (2004).

## References

• Brau, Charles A. Modern Problems in Classical Electrodynamics. Oxford University Press, 2004. ISBN 0-19-514665-4.
• Jackson, John David. Classical Electrodynamics. John Wiley & Sons, 1999. ISBN 0-471-30932-X