Consider a point source emitting radiation isotropically in all directions. By placing a spherical shells and at radii and , we can write, by conserving the energy:

If we assume one of the shells, say , is fixed, then we can write

which is the inverse square law. Now represents a quantity which is representative of all the rays of radiation passing through the surface. A single ray is a mathematical idealization and by definition has zero cross-section area to carry some of that energy. However, we may consider some infinitesimal area through which a bundle of rays is passing perpendicularly.

Now this bundle can carry energy as it has some cross-sectional area, which we can write as

We define as the specific intensity and is the angle between the vector normal to and direction of . Here it is useful to introduce a unit vector which points in the direction of the propagation of radiation. The solid angle element represents an infinitesimal range of directions centered around . Thus, is the angle between, the propagation direction and the normal to the surface element .

The way to understand this expression is quite simple. We want to know how much energy

  1. is going to cross the infinitesimally small area ,
  2. towards direction given by ,
  3. in time ,
  4. in the frequency range .
    The factor accounts for the fact that the rays do not necessarily cross the surface perpendicularly. The effective area presented to the radiation is the projected area . In particular, when the radiation is normally incident, , and the projected area is simply .

Because it gives a more detailed description of the radiation field, it is much more instructive to derive other quantities (like flux) from the specific intensity. One way to derive other quantities is to calculate the moments of .

Moments of Specific Intensity

Before we derive any moments, let us write equation for a bit more generally in terms of the unit normal . Here is the unit vector normal to the surface element . It describes the orientation of the surface, whereas describes the direction in which the radiation is traveling. These are therefore two different vectors.

The geometrical factor can be written as

so the expression above is the same as the previous expression. More precisely, we will generally write the specific intensity as , since its angular dependence is on the propagation direction , while describes the projection of the surface onto the direction of propagation.

This essentially means all the rays are aligned with . This statement refers to the special case of normally incident radiation, where . In the general case, and need not be aligned.

Next is to ask about the nature of the independent variable for which we will be calculating the moments for. For specific intensity, it is more constructive to think in terms of direction. From the definition, it is clear that is not a spatial density, but rather a function of direction1. Therefore, it makes sense to derive angular moments.

Zeroth moment (Radiation Energy Density)

By definition, zeroth moment of a quantity is something which is does not depend of our independent variable (direction in our case). So we want to derive the quantity which is describes how much radiation is present locally, regardless of the direction. This essentially defines energy density at that location.

To calculate this energy density, we can once again use the isotropic scenario, and imagine a point which is absorbing incoming radiation from all directions. If in time , they cross an area while traveling at speed , then, we can write energy density for the photons in the frequency bin as

Here we have used . The distance travelled by the radiation in time is , so the volume swept out by the radiation through the area is . Now, this energy density is only the contribution from the rays aligned with , so we can integrate over all directions to get the total energy density. More precisely, this is the contribution from radiation within the particular solid-angle element . To obtain the total energy density, we sum the contributions from all possible propagation directions.

The angular integral here is over the entire sphere of possible propagation directions, so

A closely associated quantity to is the mean intensity, which is defined as:

Consequently, we can write:

Finally, we can also calculate total radiation density

First moment (Flux)

The zeroth moment tells us how much radiation energy is present locally, but it does not tell us whether the radiation is preferentially travelling in some direction. To describe the net transport of radiation energy, we need to retain one power of the propagation direction .

Recall that is the unit normal to the surface, while is the direction of propagation of the radiation. Thus,

From

the energy crossing the surface per unit area, per unit time, and per unit frequency is

Integrating over all directions gives us the net flux through the surface:

or, equivalently,

The word net is important here. Radiation travelling in the direction of contributes positively, while radiation traveling in the opposite direction contributes negatively. For an isotropic field, is independent of direction. In this case, radiation traveling in opposite directions contributes equally but with opposite signs, and therefore

Thus, an isotropic radiation field can have a non-zero energy density but zero net flux. We can also integrate over the frequencies to get the integrated flux:

Second moment (Momentum flux)

Radiation carries both energy and momentum. For a photon with energy , its momentum has magnitude

Therefore, we can similarly ask for the rate at which radiation momentum crosses a surface.

For radiation travelling at an angle to the surface normal, there is one factor of from the projected area and another factor of because only the component of momentum normal to the surface contributes.

Thus, the momentum flux is

This is called the second moment because the angular dependence now contains two powers of the direction cosine, .

For isotropic radiation,

Using

we obtain,

so,

This is the familiar relation between the energy density and pressure of an isotropic radiation field. Like flux, we can again integrate over all frequencies to get:

Specific intensity along a ray

Specific intensity remains constant along a ray if there is no absorption and emission (in other words, through free space) by simply the virtue of conservation of energy. This does not mean flux is constant along a ray. Flux can decrease because the solid angle subtended by the source decreases, even though the intensity of each individual ray remains unchanged.

Footnotes

  1. Note that we are trying to define as a density function. So although, it does depend on , it is not a density function in . ↩