Waves in composites and metamaterials/Airy theory

The content of these notes is based on the lectures by Prof. Graeme W. Milton (University of Utah) given in a course on metamaterials in Spring 2007.

Introduction edit

The Descartes theory of rainbows cannot explain a number of observed phenomena, one of which is the existence of supernumerary bands of alternating green and purple.

The Airy theory of rainbows attempts to address that issue by considering the wave nature of light. Figure 1 shows a figure of a water droplet with the incident planar wavefront getting converted into a curved wavefront due to the spherical shape of the raindrop. Because the curved wavefront turns back on itself, there is interference between preceding and succeeding wavefronts. The regions of constructive and destructive interference give rise to supernumerary bands.

The outgoing wave front in the figure has been computed at the locus of points on the outgoing rays that are at a constant distance from the incident wave front. Airy determined that the outgoing wave appears to come from a cubic wavefront as shown in the figure.

 
Figure 1. Curvature of the wavefront due to the shape of a raindrop.

An image of the interference patterns observed in a rainbow can be seen in Figure 2.

 
Figure 2. Interference and the formation of supernumerary bands (from Adam, 2002).

Airy's Theory edit

In Airy's theory, we start with the cubic wavefront and choose the point of inflection on the cubic wavefront at the origin of coordinates. The wavefront is oriented tangent to the  -axis as shown in Figure 3.

 
Figure 3. Coordinate system for Airy's theory of the rainbow.

The equation of the cubic wavefront is (see Humphreys (1964) for details)

 

where   is the drop radius,   is the refractive index of water, and   is the number of times the ray traverses the raindrop.

Airy made the following approximations:

  1. The virtual wavefront   can be treated as a cubic.
  2. The amplitude along the wavefront is constant.
  3. The virtual wavefront is cylindrical.
  4. The wave can be treated as a scalar wave.

Note that coherence of natural light is assumed when we think of a planar wavefront with all the rays in phase.

To calculate the emergent radiation at large distances, Airy employed Huygen's principle, i.e., each point on the wavefront is treated as a source of radiation. The change in phase due to the different distances traveled by the waves is then used to compute the intensity of the rainbow.

To determine the change in phase we start of with the diagram shown in Figure 4. We wish to know the difference in phase between the waves starting at point   and those starting at point   when they reach a distant observer at point  , i.e., we want ( ) at  . The point   is far enough away that the rays   and   may be assumed to be parallel.

 
Figure 4. Geometry of phase change in Airy's theory.

From the geometry in Figure 4, we see that  . Therefore,

 

where   is the wavelength. Now,   and  . Therefore,

 

Substituting equation (1) into the above equation and rearranging, we get

 

If   is the time it takes for the wave starting from   to reach point  , and if   is the period of the wave, then the phase of the wave at   is

 

Therefore,

 

Let us assume that the waveform is a sine wave (we could alternatively use a cosine waveform). Then the amplitude at point   due to the wave starting at point   is

 

Integrating contributions over the wavefront, we find that the amplitude at   is

 

where   is an differential element along the wavefront. Assuming that   is constant and that   along the effective portion of the front, we can write

 

Define

 

Then,

 

Using the identity  , we have

 

or,

 

where

 

Because the wavelength of light is approximately 550 nm and the diameter of a raindrop is approximately 1 mm, the value of   changes rapidly with   beyond a very small distance around  . So the net contribution from most of the wavefront (except from a small region around  ) is essentially zero.

Hence, only a small error is incurred by extending the limits of integration to  . Therefore,

 

Recall that the series expansion

 

Since   is an odd function of  , the series expansion of   involves odd powers of  . Hence   is an odd function of   and its integral from   to   is zero, i.e.,   = 0.

To evaluate the integral for   we can write

 

or

 

Apply a change of variables

 

to get

 

We can write the above equation as

 

where

 

The function   is called the Airy function. Hints on how these integrals can be evaluated in terms of Bessel functions can be found in Adams (2002a).

Next, plugging back into equation (2), the amplitude at point   is

 

Since a cylindrical waveform is assumed, we can compute the intensity at point P using the relation

 

A plot of the intensity versus the scattering angle is shown in Figure 5 Khare74.

 
Figure 5. Intensity vs Scattering Angle from Airy Theory for a rainbow (from Khare, 1974).

The dark band occurs between the primary and the secondary rainbows because of the interference predicted by Airy's theory (see Figure 6). The figure shows the intensity for a single frequency of light. If multiple frequencies are plotted on the same graph, the peaks and troughs for each frequency occur at different locations and with different wavelengths (see for example Humphreys, 1964). However, when the water droplets are very small, the peaks for all the frequencies overlap (as do the troughs) and alternating bands of white and dark are observed. This accounts for the behavior of fogbows. Similar behavior is also observed for particles scattered by atoms.

 
Figure 6. Airy's theory explanation of Alexander's dark band.

For light that is polarized in certain directions, Airy's theory does not predict the correct behavior (see Khare, 1974 for examples). This is problematic because light from a rainbow is almost completely polarized. Luckily, for the dominant polarization which constitutes most of natural light, Airy's theory works remarkably well. For the opposite polarization Airy's theory predicts minima where there should be maxima. This has to do with Brewster's angle.

When an unpolarized wave of light is incident on the interface between two materials with different refractive indices, some of the energy is reflected at the boundary and some is transmitted. However, if the incident wave is polarized such that the electric field is parallel to the plane containing a ray of light and the normal to the interface, an interesting thing occurs. For a certain angle of incidence, no light is reflected from the interface. This angle is called Brewster's angle ( ) (see Figure 7).

 
Figure 7. Brewster's Angle and the polarization of light.

If   is the angle made by the transmitted wave to the normal to the interface, then at the Brewster angle

 

By Snell's law

 

Therefore, at Brewster's angle,

 

For an air-water interface with   and  , we get a value of   = 53 . This angle is very close to the angle that the rainbow ray makes with the normal to the surface of a droplet at the point where it is reflected inside the droplet. Hence only certain polarizations get reflected at this point and other polarizations are transmitted. At angles greater than the Brewster angle, the radiation is 90  out of phase with light reflected at less than the Brewster angle.

Airy's assumption that the amplitude is constant along the cubic wavefront breaks down due to this polarization.

References edit

Adam (2002)    
J. A. Adam. Like a bridge over colored water: A mathematical review of The Rainbow Bridge: Rainbows in Art, Myth, and Science. AMS Notices, 49(11):1360--1371, 2002.

Adam (2002a)    
J. A. Adam. The mathematical physics of rainbows and glories. Physics Reports, 356:229--365, 2002.

Humphreys (1964)    
W. J. Humphreys. Physics of the Air. Dover, New York, 1964.

Khare (1974)    
V. Khare and H. M. Nussenzveig. Theory of the rainbow. Physical Review Letters, 33(16):976--980, 1974.