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

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

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 ${\displaystyle x}$-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)

${\displaystyle {\text{(1)}}\qquad y={\cfrac {h}{3~a^{2}}}~x^{3}~;~~h={\cfrac {(p^{2}-1)^{2}}{p^{2}~(n^{2}-1)}}~{\sqrt {\cfrac {p^{2}-n^{2}}{n^{2}-1}}}}$

where ${\displaystyle a}$ is the drop radius, ${\displaystyle n}$ is the refractive index of water, and ${\displaystyle p}$ is the number of times the ray traverses the raindrop.

1. The virtual wavefront ${\displaystyle A'B'}$ 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 ${\displaystyle O}$ and those starting at point ${\displaystyle M}$ when they reach a distant observer at point ${\displaystyle P}$, i.e., we want (${\displaystyle \varphi _{M}-\varphi _{O}}$) at ${\displaystyle P}$. The point ${\displaystyle P}$ is far enough away that the rays ${\displaystyle MP}$ and ${\displaystyle OP}$ 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 ${\displaystyle TP\approx RP}$. Therefore,

${\displaystyle \varphi _{M}-\varphi _{O}={\cfrac {2~\pi }{\lambda }}(OR-MT)}$

where ${\displaystyle \lambda }$ is the wavelength. Now, ${\displaystyle OR=x~\sin \theta }$ and ${\displaystyle MT=y~\cos \theta }$. Therefore,

${\displaystyle \varphi _{M}-\varphi _{O}={\cfrac {2~\pi }{\lambda }}(x~\sin \theta -y~\cos \theta )~.}$

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

${\displaystyle \varphi _{M}=\varphi _{O}-{\cfrac {2~\pi }{\lambda }}\left({\cfrac {h~x^{3}}{3~a^{2}}}~\cos \theta -x~\sin \theta \right)~.}$

If ${\displaystyle t}$ is the time it takes for the wave starting from ${\displaystyle O}$ to reach point ${\displaystyle P}$, and if ${\displaystyle T}$ is the period of the wave, then the phase of the wave at ${\displaystyle P}$ is

${\displaystyle \varphi _{O}={\cfrac {2~\pi ~t}{T}}~.}$

Therefore,

${\displaystyle \varphi _{M}=2~\pi ~\left[{\cfrac {t}{T}}-{\cfrac {1}{\lambda }}\left({\cfrac {h~x^{3}}{3~a^{2}}}~\cos \theta -x~\sin \theta \right)\right]~.}$

Let us assume that the waveform is a sine wave (we could alternatively use a cosine waveform). Then the amplitude at point ${\displaystyle P}$ due to the wave starting at point ${\displaystyle M}$ is

${\displaystyle A_{M}=K_{M}~\sin \varphi _{M}=K_{M}~\sin \left\{2~\pi ~\left[{\cfrac {t}{T}}-{\cfrac {1}{\lambda }}\left({\cfrac {h~x^{3}}{3~a^{2}}}~\cos \theta -x~\sin \theta \right)\right]\right\}~.}$

Integrating contributions over the wavefront, we find that the amplitude at ${\displaystyle P}$ is

${\displaystyle V=\int A_{M}~ds=\int K_{M}~\sin \left\{2~\pi ~\left[{\cfrac {t}{T}}-{\cfrac {1}{\lambda }}\left({\cfrac {h~x^{3}}{3~a^{2}}}~\cos \theta -x~\sin \theta \right)\right]\right\}~{\text{d}}s}$

where ${\displaystyle {\text{d}}s}$ is an differential element along the wavefront. Assuming that ${\displaystyle K=K_{M}}$ is constant and that ${\displaystyle y\ll x}$ along the effective portion of the front, we can write

${\displaystyle V\approx K~\int \sin \left\{2~\pi ~\left[{\cfrac {t}{T}}-{\cfrac {1}{\lambda }}\left({\cfrac {h~x^{3}}{3~a^{2}}}~\cos \theta -x~\sin \theta \right)\right]\right\}~{\text{d}}x~.}$

Define

${\displaystyle \delta (x):={\cfrac {2~\pi }{\lambda }}\left({\cfrac {h~x^{3}}{3~a^{2}}}~\cos \theta -x~\sin \theta \right)~.}$

Then,

${\displaystyle V=K~\int \sin \left({\cfrac {2~\pi ~t}{T}}-\delta \right)~{\text{d}}x~.}$

Using the identity ${\displaystyle \sin(a-b)=\sin a~\cos b-\cos a~\sin b}$, we have

${\displaystyle V=K~\int \left(\sin {\cfrac {2~\pi ~t}{T}}~\cos \delta -\cos {\cfrac {2~\pi ~t}{T}}~\sin \delta \right)~{\text{d}}x}$

or,

${\displaystyle {\text{(2)}}\qquad V=A~\sin {\cfrac {2~\pi ~t}{T}}-B~\cos {\cfrac {2~\pi ~t}{T}}}$

where

${\displaystyle A:=K~\int \cos \delta ~{\text{d}}x~;~~B:=K~\int \sin \delta ~{\text{d}}x~.}$

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

Hence, only a small error is incurred by extending the limits of integration to ${\displaystyle [-\infty ,\infty ]}$. Therefore,

${\displaystyle A\approx K~\int _{-\infty }^{\infty }\cos \delta ~{\text{d}}x~;~~B\approx K~\int _{-\infty }^{\infty }\sin \delta ~{\text{d}}x~.}$

Recall that the series expansion

${\displaystyle \sin \delta =\delta -{\cfrac {\delta ^{3}}{3!}}+{\cfrac {\delta ^{5}}{5!}}+\dots ~.}$

Since ${\displaystyle \delta }$ is an odd function of ${\displaystyle x}$, the series expansion of ${\displaystyle \sin \delta }$ involves odd powers of ${\displaystyle x}$. Hence ${\displaystyle \sin \delta }$ is an odd function of ${\displaystyle x}$ and its integral from ${\displaystyle -\infty }$ to ${\displaystyle \infty }$ is zero, i.e., ${\displaystyle B}$ = 0.

To evaluate the integral for ${\displaystyle A}$ we can write

${\displaystyle A\approx 2~K~\int _{0}^{\infty }\cos \delta ~{\text{d}}x}$

or

${\displaystyle A\approx 2~K~\int _{0}^{\infty }\cos \left[{\cfrac {2~\pi }{\lambda }}\left({\cfrac {h~x^{3}}{3~a^{2}}}~\cos \theta -x~\sin \theta \right)\right]~{\text{d}}x~.}$

Apply a change of variables

${\displaystyle z={\cfrac {4~x}{\lambda ~u}}~\sin \theta ~;~~u=\left({\cfrac {4~h}{3~a^{2}~\lambda }}~x^{3}~\cos \theta \right)^{1/3}}$

to get

${\displaystyle A\approx 2~K~\int _{0}^{\infty }\left({\cfrac {3~a^{2}~\lambda }{4~h~\cos \theta }}\right)^{1/3}~\cos \left[{\cfrac {\pi }{2}}~(u^{3}-z~u)\right]~du~.}$

We can write the above equation as

${\displaystyle A\approx M~f(z)}$

where

${\displaystyle M:=2~K\left({\cfrac {3~a^{2}~\lambda }{4~h~\cos \theta }}\right)^{1/3}~;~~f(z):=\int _{0}^{\infty }\cos \left[{\cfrac {\pi }{2}}~(u^{3}-z~u)\right]~du~.}$

The function ${\displaystyle f(z)}$ 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 ${\displaystyle P}$ is

${\displaystyle V=M~f(z)~\sin {\cfrac {2~\pi ~t}{T}}~.}$

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

${\displaystyle I=A^{2}=M^{2}~[f(z)]^{2}~.}$

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 (${\displaystyle \theta _{B}}$) (see Figure 7).

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

If ${\displaystyle \theta _{t}}$ is the angle made by the transmitted wave to the normal to the interface, then at the Brewster angle

${\displaystyle \theta _{B}+\theta _{t}={\cfrac {\pi }{2}}~.}$

By Snell's law

${\displaystyle {\cfrac {\sin \theta _{t}}{\sin \theta _{B}}}={\cfrac {n_{1}}{n_{2}}}~.}$

Therefore, at Brewster's angle,

${\displaystyle \cot \theta _{B}={\cfrac {n_{1}}{n_{2}}}~.}$

For an air-water interface with ${\displaystyle n_{1}=1}$ and ${\displaystyle n_{2}=1.33}$, we get a value of ${\displaystyle \theta _{B}}$ = 53${\displaystyle ^{o}}$. 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${\displaystyle ^{o}}$ 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.