This tutorial is meant to complement the first lecture of the course High Energy Astrophysics. Here you can practice your knowledge of the topic after completing Lecture 1.

Sometimes you will find some collapsable/expandable paragraphs marked in blue (questions) or orange (quizzes). These are quick tests to check your knowledge as you follow the tutorial. Click on the paragraph to read the answer.

Radiative Transfer Theory Basics

Radiative Transfer Approximation

Electromagnetic radiation is described by the Maxwell equations. Energy is transported by a transverse wave in the electric and magnetic fields. The Maxwell equations need to be solved in order to find how these waves interact and propagate. When photons propagate through macroscopic media, the solutions of the Maxwell equations become too complex to solve for all practical purposes. There are two simplifications that can be done without losing generality:

  • The first is that of geometric optics, with photons approximated as light rays. Here the scale of a system L greatly exceeds the wavelength of radiation \lambda: L \gg \lambda. The wave nature of photons is considered only to the extent that it introduces light bending on surfaces and changing of refraction index.
  • If the index of refraction is n=\rm constant everywhere (or at least it can be safely approximated as n=\rm constant everywhere), then we talk about radiative transfer, which represents a further simplification with respect to geometric optics.

Specific Intensity Definition

Exercise 1

Let’s consider a pinhole camera, which is nothing but a small hole in a box, able to generate an image without using a lens. Let’s assume that the hole is circular with a diameter d. Let’s say that the box has a depth L so that the rays have to cross a distance L \gg d before reaching the image plane. If a source of light with a specific intensity (sometimes also called surface brightness) I_\nu is placed at a certain distance from the pinhole, then we can calculate what specific flux F_\nu is received at the image plane.

Pinhole camera

What is the meaning of “specific” in the specific intensity and flux?

In this exercise, we want to calculate the flux F_\nu as a function of I_\nu(\theta, \phi). Here \theta represents the angle in the x-y plane between the source and the pinhole, whereas \phi is the azimuthal angle. We want to show that:

    \[ F_\nu = \frac{\pi\,cos^4\theta}{4\,f^2}I_\nu(\theta, \phi)\]

where f = L/d is the focal length.

Solution 1

Before reading the full solution, consider the following hints.

Hint 1

Hint 2

Solution

Exercise 2

A thin disk of material is emitting radiation with specific intensity I_\nu. The disk is “thin” in the sense that its height H \ll R, where R is the disk radius. The radiation emitted is collected by a detector placed at a distance r \gg R. The angle between the line of sight and the direction normal to the disk is \beta, whereas the disk is seen at an angle \alpha. Calculate the flux measured by the detector. The solution should be a function of both \alpha and \beta.

The geometry of the problem.

Solution 2

As usual, let’s start with some hints.

Hint 1

Hint 2

Hint 3

Solution

The Constancy of Intensity in Free Space

If the rays propagate in free space, then the (specific) intensity does not change with distance. This might sound counter-intuitive at first. Indeed we know very well that the inverse square law tells us that, in free space, the flux of an object decreases as the square of the distance. So why do we say that the surface intensity remains constant? The keyword here is surface, meaning that we are observing an extended object and not a point source. Remember that the apparent size of an extended object, i.e., its solid angle, decreases with the square of the distance so that you squeeze more and more surface elements into a given solid angle. So the “number” of surface elements emitting light in the field of view increases as the square of the distance, whereas the flux we receive decreases by the square of the distance. The two effects cancel each other out and we have a constant surface intensity.

In the image on the left, the spheres represent an extended object emitting light at distances D=1, 2, and 4. The circle is the detector (e.g., eye, telescope, camera, etc.) with a fixed size. As the distance increases the surface area of the sphere scales as D^{-2}. However, more surface elements fall within the same solid angle, thus preserving the surface intensity.

Surface Intensity and the Moon

It might be worth insisting on this point a bit further with an even simpler example. Imagine walking in the countryside and looking at a row of nearly identical trees in perspective, like in the figure below.

Suppose you are close to the first tree, about 5 meters away, whereas the last tree you see is at about 500 meters. Since the flux scales as D^{-2}, you receive about 10,000 fewer photons from the last tree than from the first tree. However, you use the same exposure time to record the optimal picture, be it with a camera or with your eyes. You don’t see the trees fading away with distance!

This happens because all trees are roughly illuminated in the same way, so they all have the same surface brightness. Their flux decreases as D^{-2}, but so does their solid angle and the two effects cancel out.

Look at the picture in the panel below. The trees here are illuminated by a flashlight. Their intensity clearly decreases with distance. Why?

(Credits: Saturated)

Solution

Exercise 3

The specific intensity (or surface brightness) of an object remains constant in free space. Give now an expression that describes how the specific intensity changes when there is some absorbing material along the line of sight (to be clear: you need to write the radiative transfer equation when only absorption is present and solve the equation by giving an expression for I_\nu(s)). Call the absorption coefficient \alpha_\nu and the initial specific intensity as I_{\nu,0}.
(NOTE: you can also choose to write and solve the transport equation in terms of the optical depth \tau_\nu rather than \alpha_\nu and s).

Do you remember the definition of the absorption coefficient?

Solution 3

Let’s start by writing the equation of transport for the absorption-only case. This means we will set the emission coefficient j_\nu to zero and ignore any scattering.

What are the dimensions of the absorption coefficient?

The differential expression for the equation of transport is therefore d I_\nu = -\alpha_\nu I_\nu ds. Integrating this expression we get I_\nu(s) = I_\nu(s_0) e^{-\int_{s_0}^{s}\alpha_\nu(s')ds'.

Source Function

In the previous exercise, we mentioned that you can use the optical depth in place of the absorption coefficient \alpha_\nu and the length ds. Remember that the infinitesimal optical depth, for the case of absorption only, is defined as d\tau_\nu = \alpha_\nu\cdot ds.

Can you write the transfer equation as a function of the optical depth?

The ratio between emission coefficient and absorption coefficient is called source function S_\nu = \frac{j_\nu}{\alpha_\nu} and has the same dimensions of specific intensity. The source function works very well together with the optical depth, since it tells how the specific intensity varies along the optical depth scale.

The source function is not only a useful quantity to rewrite the equation of radiative transport in a simpler way, but it has also a profound physical meaning. The source function, having the dimensions of a specific intensity, can be thought of as the local input of radiation. For example when the optical depth \tau_\nu \gg 1 then I_\nu \approx S_\nu which means that matter far from the location of interest is giving a negligibly small contribution to the specific intensity. In other words, the source function is the value approached by the specific intensity given sufficient optical depth.

The Meaning of the Source Function

Let’s illustrate the meaning of the source function with an astrophysical example. An astrophysical source has intrinsic specific brightness I_{\nu,0}=10^{-6} \rm\,erg\,s^{-1}\,cm^{-2}\,Hz^{-1}\,ster^{-1} and emits a number of rays that cross a medium with optical depth \tau_\nu. After the usual absorption and emission processes, the modified specific brightness reaches the observer that detects the photons with a receiver. For the medium, let’s choose the sky and approximate it as a uniform slab of gas whose optical thickness is \ll 1, for example, \tau_\nu = 0.02 at \nu = 30 GHz. Let’s assume that the specific intensity of the sky at 30 GHz is 1.7\times 10^{-7}\rm\,erg\,s^{-1}\,cm^{-2}\,Hz^{-1}\,ster^{-1}. From the equation of transport, what would be the specific brightness observed at 30 GHz with respect to the initial specific brightness I_{\nu,0}?

Hint 1

Solution