Import standard modules:


In [2]:
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
from IPython.display import HTML 
HTML('../style/course.css') #apply general CSS


Out[2]:

Import section specific modules:


In [3]:
from IPython.display import Image

In [4]:
HTML('../style/code_toggle.html')


Out[4]:
The raw code for this notebook is by default hidden for easier reading. To toggle on/off the raw code, click here.

1.6.1 Synchrotron Emission:

Sychrotron emission is one of the most commonly encountered forms of radiation found from astronomical radio sources. This type of radiation originates from relativistic particles get accelerated in a magnetic field.

An understanding of generation of Synchrotron emission depends fundamentally on special relativistic effects. Here we wil not go into the details; but try to give an (somewhat handwaving) idea of the underlying physics. An accelerating charge emitts radiation (as we have seen in section 1.2.1) - here the acceleration is done through the ambient magnetic fields. The non-relativistic Larmor formula for the radiated power is:

$$P= \frac{2}{3}\frac{q^{2}a^{2}}{c^{3}}$$

.

Since the acceleration is a result of the magnetic fields, we get:

$$P=\frac{2}{3}\frac{q^{2}}{c^{3}}\frac{v_{\perp}^{2}B^{2}q^{2}}{m^{2}} \qquad$$

where $v_{\perp}$ is the component of velocity of the particle perpendicular to the magnetic field, $m$ is the mass of the charged particle, $q$ is it's charge and $a$ is the acceleration it is undergoing. This is essentially the cyclotron radiation. The relativistic effects modify this into :

$$P = \gamma^{2} \frac{2}{3}\frac{q^{2}}{c^{3}}\frac{v_{\perp}^{2}B^{2}q^{2}}{m^{2}c^{2}} = \gamma^{2} \frac{2}{3}\frac{q^{4}}{c^{3}}\frac{v_{\perp}^{2}B^{2}}{m^{2}c^{2}} \qquad ,$$

where $$\gamma = \frac{1}{\sqrt{1+v^{2}/c^{2}}} = \frac{E}{mc^{2}} \qquad ,$$

is a measure of the energy of the particle. For particles to be considered relativistic and ultrarelativistic, $\gamma >> 1$, typically with a value of $100-1000$. Since $v_{\perp}=vsin\alpha$, with $\alpha$ being the angle between the magnetic field and the velocity of the particle, the radiated power can be written as:

$$P=\gamma^{2} \frac{2}{3}\frac{q^{4}}{c^{3}}\frac{v^{2}B^{2}sin\alpha^{2}}{m^{2}c^{2}} \qquad .$$

From this equation it can be seen that the total power radiated by the particle depends on the strength of the magnetic field and that the higher the energy of the particle, the more the power it radiates away.

In analogy with the non-relativistic case, there is a frequency of gyration. This refers to the path the charged particle follows while being accelerated in a magnetic field. The figure below illustrates the idea.


In [5]:
Image(filename='figures/drawing.png', width=300)


Out[5]:

The frequency of gyration in the non-relativistic case is simply

$$\omega = \frac{qB}{mc} \qquad .$$

For synchrotron radiation, this gets modified to

$$\omega_{G}= \frac{qB}{\gamma mc} \qquad ,$$

since it is the relativistic mass, i.e. $\gamma mc$ which we should consider in this case.

In the non-relativistic case (i.e. cyclotron radiation), the frequency of gyration is also the frequency of radiation. If this was also the case for the synchrotron radiation, for magnetic fields of galactic strengths (a few microGauss or so), the resultant frequency would be less than a Hertz ! Luckily for us, the relativistic beaming and doppler effects come into play; the actual frequency of the observed radiation gets bumped up by a factor of ~ $\gamma^{3}$, which brings it in the radio regime. This frequency, known also as the 'critical frequency' is at most of the emission takes place. It is given by

$$\nu_{c} \propto \gamma^{3}\nu_{G} \propto E^{2}$$

.

So far we have discussed a single particle emitting synchrotron radiation. However, what we want to know is the case of an ensemble of raditing particles. Since the Synchrotron emission depends only on the magnetic field and the energy of the particle, given a (more or less) uniform magnetic field, what we need is a distribution of the particles with energy. Denoting the distribution N(E) (number of particles per unit volume per solid angle) of the particles, the resultant spectrum from the ensemble of the particles is: $$ \epsilon(E) dE = N(E) P(E) dE $$.

The usual assumption made about the distribution N(E) (based also on the observed cosmic ray distribution) is that of a power law, i.e.

$$N(E)dE=E^{-\alpha}dE \qquad .$$

Plugging in this and remembering that $P(E) \propto \gamma^{2} \propto E^{2}$, we get

$$ \epsilon(E) dE \propto E^{2-\alpha} dE \qquad .$$

Shifting to frequncy domain, $$\epsilon(\nu) \propto \nu^{(1-\alpha)/2} \qquad .$$

The usual value for $\alpha$ is 5/2 and since flux $S_{\nu} \propto \epsilon_{\nu}$, $$S_{\nu} \propto \nu^{-0.75} \qquad .$$

This shows that the synchrotron flux is also a power law, if the underlying distribution of particles is a power law.

This is approximately valid for 'fresh' collection of radiating particles. However, as mentioned above, the higher energy particles lose energy through radiation much faster than the lower energy particles. This means that the distribution of particles over time gets steeper at higher frequncies (which is where the contribution from the high energy particles comes in). This steepening of the spectral index is a typical feature of older plasma in astrophysical scenarios.

1.2.5.1 Sources of Synchrotron Emission:

So where do we actually see synchrotron emission? As mentioned above, the prerequisites are magnetic fields and relativistic particles. These conditions are satisfied in a variety of situations. A prime example is that of lobes of radio galaxies. The lobes contain relativistic plasma in magnetic fields of strength ~ $\mu$G. The origin of the plasma and the magnetic field is thought to be in ultimately in the activity taking place in the central part of the radio galaxies, which hosts a supermassive black hole. The figure below shows a radio image of the radio galaxy nearest to us, Cygnus A. The jets, which carry relativistic plasma/particles originate from the centre of the host galaxy (emission from there is marked as 'core'), collide with the surrounding medium at the places which show bright spots of emission, labelled as "hotspots". The plasma forming the radio lobe streams back from the hotspots, which means that the 'freshest' plasma resides in the hotspots and the surrounding regions of the lobes; by contrast, the plasma in lobes closest to the centre is the oldest.


In [6]:
Image(filename='figures/cygnusA.png')


Out[6]:

The non-thermal nature of the emission can be seen easily by measuring the spectrum of the radio emission from the lobes. The plots below show the spectrum of the lobes of Cygnus A, from a frequency of 150 MHz to 14.65 GHz. It can be seen that these can be explained better by a synchrotron spectrum rather than a spectrum of thermal origins. Another observation can be made from the plot that the spectral index of the lobe emission seems steeper than the flux from the hotspot. This is consistent with fresh plasma being the one from the hotspot and the older plasma residing in the lobes.


In [7]:
# Data taken from Steenbrugge et al.,2010, MNRAS
freq=(151.0,327.5,1345.0,4525.0,8514.9,14650.0)
flux_L=(4746,2752.7,749.8,189.4,83.4,40.5)
flux_H=(115.7,176.4,69.3,45.2,20.8,13.4)
fig,ax = plt.subplots()
ax.loglog(freq,flux_L,'bo--',label='Lobe Flux')
ax.loglog(freq,flux_H,'g*-',label='Hotspot Flux')
ax.legend()
ax.set_xlabel("Frequency (MHz)")
ax.set_ylabel("Flux (Jy)")


Out[7]:
<matplotlib.text.Text at 0x7f6b485ba690>

As mentioned before, Synchrotron emission is an example of non-thermal continuum emission - emission occuring over a wide range of frequencies. In the next section, we'll see an example of spectral line emission - emission which is confined to a narrow range in frequencies. The particular example which we'll study is the famous 21 cm line.

1.6.2 Inverse Compton:

Another important mechanism of radiation for radio astronomical objects is the inverse Compton process. The Compton process is one through which a high energy photon loses energy through scattering. The inverse Compton, as the name implies, is the process through which a low energy photon (in this context, a radio photon) gets upscattered to higher energies (and consequently, higher frequencies). This process in important in several astronomical scenarios.

We'll again resort to handwaving arguments to give a flavour of the physics and arrive at the final result. We start with a situation analoguos to the derivation of Bremsstrahlung radiation. This is a situation where photon-electron interaction is occuring, instead of the scattering of an electron from an ion.

The 'acceleration' of the electron in this case is: $$\mathbf{a} = \frac{e}{m_{e}}\mathbf{E'} $$, with $\mathbf{E'}$ being the electric field of the incident photon (Note that this is in the rest frame of the electron). Then the Larmor formula gives the total radiated power as: $$P = \frac{2\,e^{4}}{3m_{e}^{2}c^{3}} \mathbf{E'}^{2}$$. In the lab frame, this turns out to be: $$P=2\frac{8\pi}{3}\frac{e^{4}}{m_{e0}^{2}c^{4}}\gamma^{2}(\mathbf{E}+\frac{\mathbf{v}}{\mathbf{c}}\times\,\mathbf{B})^{2}\frac{c}{8\pi}$$, where we apply the Lorentz transform for the electric field. The expression for the radiated power can be put in a more elegant form: $$P = \sigma_{T}\gamma^{2}c\,u_{rad}$$, where $u_{rad}$ is the energy density of the photon and $\sigma_{T}$ is the Thomson cross section of the electron. This is the energy lost by the electron and it is also the energy gained by the photon in the scattering process, so: $$\frac{\Delta E}{E} = (\gamma^{2}-1)\sim \frac{v^{2}}{c^{2}}$$

This process takes place with a low energy photon and high energy plasma. One of the main examples is photons from the cosmic microwave background interacting with hot intracluster medium. As we learned from $\S$ 1.5, the CMB has a temperature $T = 2.725 K $. Hot intracluster medium, on the other hand, has a temperature of $\sim 10^7 - 10^8$ K. The photons from the CMB interact with this medium, which results in the photons getting shifted to higher energy. This effectively marks a change in the CMB temperature in the direction of the cluster.

The change in the temperature can be written as: $$\frac{\Delta T}{T} = \frac{2kT}{m_{e}c^{2}}\sigma_{T}\,N_{e}L$$.