Diff
Diff
In this chapter we extend the ideas of the previous chapter to the case of
waves in more than one dimension. The extension of the sine wave to higher
dimensions is the plane wave. Wave packets in two and three dimensions
arise when plane waves moving in different directions are superimposed.
Diffraction results from the disruption of a wave which is impingent upon
an object. Those parts of the wave front hitting the object are scattered,
modified, or destroyed. The resulting diffraction pattern comes from the
subsequent interference of the various pieces of the modified wave. A knowl-
edge of diffraction is necessary to understand the behavior and limitations of
optical instruments such as telescopes.
Diffraction and interference in two and three dimensions can be manipu-
lated to produce useful devices such as the diffraction grating.
1
CHAPTER 1. WAVES IN TWO AND THREE DIMENSIONS 2
y
Paul
By
B
C
Cy George
Ay A
Mary x
Ax Bx
Cx
Ax = A cos θ
A
Ay = A sinθ
θ
x
Figure 1.2: Definition sketch for the angle θ representing the orientation of
a two dimensional vector.
CHAPTER 1. WAVES IN TWO AND THREE DIMENSIONS 3
the idea of vector addition. The tail of vector B is collocated with the head
of vector A, and the vector which stretches from the tail of A to the head of
B is the sum of A and B, called C in figure 1.1.
The quantities Ax , Ay , etc., represent the Cartesian components of the
vectors in figure 1.1. A vector can be represented either by its Cartesian
components, which are just the projections of the vector onto the Cartesian
coordinate axes, or by its direction and magnitude. The direction of a vector
in two dimensions is generally represented by the counterclockwise angle of
the vector relative to the x axis, as shown in figure 1.2. Conversion from one
form to the other is given by the equations
Cx = Ax + Bx Cy = Ay + By . (1.3)
Ax = Cx − Bx Ay = Cy − By . (1.4)
x
y
B
A
θ
Ax
where θ is the angle between the two vectors. In two dimensions an alternate
expression for the dot product exists in terms of the Cartesian components
of the vectors:
A · B = Ax Bx + Ay By . (1.6)
It is easy to show that this is equivalent to the cosine form of the dot product
when the x axis lies along one of the vectors, as in figure 1.3. Notice in
particular that Ax = |A| cos θ, while Bx = |B| and By = 0. Thus, A · B =
|A| cos θ|B| in this case, which is identical to the form given in equation (1.5).
All that remains to be proven for equation (1.6) to hold in general is to
show that it yields the same answer regardless of how the Cartesian coor-
dinate system is oriented relative to the vectors. To do this, we must show
that Ax Bx + Ay By = A′x Bx′ + A′y By′ , where the primes indicate components
in a coordinate system rotated from the original coordinate system.
CHAPTER 1. WAVES IN TWO AND THREE DIMENSIONS 5
y
y’
R
Y’ x’
Y
X’ θ
b
a x
X
Figure 1.4: Definition figure for rotated coordinate system. The vector R has
components X and Y in the unprimed coordinate system and components
X ′ and Y ′ in the primed coordinate system.
Figure 1.4 shows the vector R resolved in two coordinate systems rotated
with respect to each other. From this figure it is clear that X ′ = a + b.
Focusing on the shaded triangles, we see that a = X cos θ and b = Y sin θ.
Thus, we find X ′ = X cos θ + Y sin θ. Similar reasoning shows that Y ′ =
−X sin θ + Y cos θ. Substituting these and using the trigonometric identity
cos2 θ + sin2 θ = 1 results in
A′x Bx′ + A′y By′ = (Ax cos θ + Ay sin θ)(Bx cos θ + By sin θ)
+ (−Ax sin θ + Ay cos θ)(−Bx sin θ + By cos θ)
= Ax Bx + Ay By (1.7)
thus proving the complete equivalence of the two forms of the dot product
as given by equations (1.5) and (1.6). Multiply out the above expression to
verify this.
A numerical quantity that doesn’t depend on which coordinate system
is being used is called a scalar . The dot product of two vectors is a scalar.
However, the components of a vector, taken individually, are not scalars,
since the components change as the coordinate system changes. Since the
laws of physics cannot depend on the choice of coordinate system being used,
CHAPTER 1. WAVES IN TWO AND THREE DIMENSIONS 6
wave vector
λ
|k|
wave fronts
Figure 1.5: Definition sketch for a plane sine wave in two dimensions. The
wave fronts are constant phase surfaces separated by one wavelength. The
wave vector is normal to the wave fronts and its length is the wavenumber.
we insist that physical laws be expressed in terms of scalars and vectors, but
not in terms of the components of vectors.
In three dimensions the cosine form of the dot product remains the same,
while the component form is
A · B = Ax Bx + Ay By + Az Bz . (1.8)
The equation for the displacement associated with a plane sine wave in three
dimensions at some instant in time is
Since wave fronts are lines or surfaces of constant phase, the equation defining
a wave front is simply k · x = const.
In the two dimensional case we simply set kz = 0. Therefore, a wave front,
or line of constant phase φ in two dimensions is defined by the equation
This can be easily solved for y to obtain the slope and intercept of the wave
front in two dimensions.
As for one dimensional waves, the time evolution of the wave is obtained
by adding a term −ωt to the phase of the wave. In three dimensions the
wave displacement as a function of both space and time is given by
The frequency depends in general on all three components of the wave vector.
The form of this function, ω = ω(kx , ky , kz ), which as in the one dimensional
case is called the dispersion relation, contains information about the physical
behavior of the wave.
Some examples of dispersion relations for waves in two dimensions are as
follows:
kx kx kx
ω ω ω
θ
−π/2 π/2
1/2 1/2
(k x2 + ky2 ) (k x2 + ky2 )
Figure 1.6: Contour plots of the dispersion relations for three kinds of waves
in two dimensions. In the upper panels the curves show lines or contours
along which the frequency ω takes on constant values. Contours are drawn
for equally spaced values of ω. For light and ocean waves the frequency
depends only on the magnitude of the wave vector, whereas for gravity waves
it depends only on the wave vector’s direction, as defined by the angle θ in
the upper right panel. These dependences for each wave type are illustrated
in the lower panels.
Brunt-Väisälä frequency.
Contour plots of these dispersion relations are plotted in the upper pan-
els of figure 1.6. These plots are to be interpreted like topographic maps,
where the lines represent contours of constant elevation. In the case of fig-
ure 1.6, constant values of frequency are represented instead. For simplicity,
the actual values of frequency are not labeled on the contour plots, but are
represented in the graphs in the lower panels. This is possible because fre-
quency depends only on wave vector magnitude (kx2 + ky2 )1/2 for the first two
examples, and only on wave vector direction θ for the third.
−∆k ∆k
k0
k1 k2
λ
w
Figure 1.7: Wave fronts and wave vectors (k1 and k2 ) of two plane waves with
the same wavelength but oriented in different directions. The vertical bands
show regions of constructive interference where wave fronts coincide. The
vertical regions in between have destructive interference, and hence define
the lateral boundaries of the beams produced by the superposition. The
quantities k0 and ∆k are also shown.
Applying the trigonometric identity for the sine of the sum of two angles (as
we have done previously), equation (1.15) can be reduced to
where
k0 = (k1 + k2 )/2 ∆k = (k2 − k1 )/2. (1.17)
This is in the form of a sine wave moving in the k0 direction with phase
speed cphase = ω/|k0 | and wavenumber |k0 |, modulated in the ∆k direction
by a cosine function. The lines of destructive interference are normal to
∆k. The distance w between lines of destructive interference is the distance
between successive zeros of the cosine function in equation (1.16), implying
CHAPTER 1. WAVES IN TWO AND THREE DIMENSIONS 11
35
30
25
20
y
15
10
0
−30 −20 −10 0 10 20 30
x
Figure 1.8: Example of beams produced by two plane waves with the same
wavelength moving in different directions. The wave vectors of the two waves
are k = (±0.1, 1.0). Regions of positive displacement are lighter, while re-
gions of negative displacement are darker.
w = π/|∆k|. (1.18)
−∆k
y
∆k
k0
k1
k2
λ1 λ2
w x
Figure 1.9: Wave fronts and wave vectors (k1 and k2 ) of two plane waves with
different wavelengths oriented in different directions. The slanted bands show
regions of constructive interference where wave fronts coincide. The slanted
regions in between have destructive interference, and as previously, define the
lateral limits of the beams produced by the superposition. The quantities k0
and ∆k are also shown.
35
30
25
20
y
15
10
0
−30 −20 −10 0 10 20 30
x
Figure 1.10: Example of beams produced by two plane waves with wave
vectors differing in both direction and magnitude. The wave vectors of the
two waves are k1 = (−0.1, 1.0) and k2 = (0.1, 0.9). Regions of positive
displacement are lighter, while regions of negative displacement darker.
with the same frequency. In general, the beam orientation (and the lines of
constructive interference) are not perpendicular to the wave fronts. This only
occurs when the wave frequency is independent of wave vector direction.
ky
orientation of
wave fronts
orientation of
beam
−∆k
∆k
k0
k1
k2
kx
const.
frequency
curve
Figure 1.11: Illustration of factors entering the addition of two plane waves
with the same frequency. The wave fronts are perpendicular to the vector av-
erage of the two wave vectors, k0 = (k1 +k2 )/2, while the lines of constructive
interference, which define the beam orientation, are oriented perpendicular
to the difference between these two vectors, ∆k = (k2 − k1 )/2.
ky
α
kx
Figure 1.12: Illustration of wave vectors of plane waves which might be added
together.
CHAPTER 1. WAVES IN TWO AND THREE DIMENSIONS 15
40
k =1.0; max =0.8
35
30
25
y
20
15
10
0
30 20 10 0 10 20 30
x
Figure 1.13: Plot of the displacement field h(x, y) from equation (1.19) for
αmax = 0.8 and k = 1.
X
= hi sin[kx sin(αi ) + ky cos(αi )] (1.19)
i
where we have assumed that kxi = k sin(αi ) and kyi = k cos(αi ). The param-
eter k = |k| is the magnitude of the wave vector and is the same for all the
waves. Let us also assume in this example that the amplitude of each wave
component decreases with increasing |αi |:
40
k =1.0; max =0.2
35
30
25
y
20
15
10
0
30 20 10 0 10 20 30
x
Figure 1.14: Plot of the displacement field h(x, y) from equation (1.19) for
αmax = 0.2 and k = 1.
the subsequent spread of the wave is much smaller than in the case of figure
1.13.
We conclude that a superposition of plane waves with wave vectors spread
narrowly about a central wave vector which points in the y direction (as in
figure 1.14) produces a beam which is initially broad in x but for which the
breadth increases only slightly with increasing y. However, a superposition
of plane waves with wave vectors spread more broadly (as in figure 1.13)
produces a beam which is initially narrow in x but which rapidly increases
in width as y increases.
The relationship between the spreading angle αmax and the initial breadth
of the beam is made more understandable by comparison with the results for
the two-wave superposition discussed at the beginning of this section. As
indicated by equation (1.18), large values of kx , and hence α, are associated
with small wave packet dimensions in the x direction and vice versa. The
superposition of two waves doesn’t capture the subsequent spread of the
beam which occurs when many waves are superimposed, but it does lead to
a rough quantitative relationship between αmax (which is just tan−1 (kx /ky ) in
the two wave case) and the initial breadth of the beam. If we invoke the small
CHAPTER 1. WAVES IN TWO AND THREE DIMENSIONS 17
angle approximation for α = αmax so that αmax = tan−1 (kx /ky ) ≈ kx /ky ≈
kx /k, then kx ≈ kαmax and equation (1.18) can be written w = π/kx ≈
π/(kαmax ) = λ/(2αmax ). Thus, we can find the approximate spreading angle
from the wavelength of the wave λ and the initial breadth of the beam w:
x
m = +1.5
slit A L1
m = +1
L θ
111
000 m = +0.5
000
111 L2
incident 000
111
000
111
wave front d 111
000
000
111
000
111
m=0
000
111
000
111
000
111 m = −0.5
d sin θ
m = −1
slit B
L0
m = −1.5
Figure 1.16: Definition sketch for the double slit. Light passing through slit
B travels an extra distance to the screen equal to d sin θ compared to light
passing through slit A.
the spreading angle of the beams is large and the diffraction pattern from
each slit individually is a cylindrical wave spreading out in all directions, as
illustrated in figure 1.13. The cylindrical waves from the two slits interfere,
resulting in oscillations in wave intensity at the screen on the right side of
figure 1.16.
Constructive interference occurs when the difference in the paths traveled
by the two waves from their originating slits to the screen, L2 − L1 , is an
integer multiple m of the wavelength λ: L2 − L1 = mλ. If L0 ≫ d, the lines
L1 and L2 are nearly parallel, which means that the narrow end of the dark
triangle in figure 1.16 has an opening angle of θ. Thus, the path difference
between the beams from the two slits is L2 −L1 = d sin θ. Substitution of this
into the above equation shows that constructive interference occurs when
d sin θ = mλ, m = 0, ±1, ±2, . . . (two slit interference). (1.22)
Destructive interference occurs when m is an integer plus 1/2. The integer
m is called the interference order and is the number of wavelengths by which
the two paths differ.
2.0
1.5
1.0
0.5
0.010 5 0 5 10
x
8
6
4
2
0 10 5 0 5 10
x
250
200
intensity
150
100
50
0 10 5 0 5 10
x
λ
k
π /4
x
in the grating:
∆λ 1
> . (1.24)
λ 2mn
1.7 Problems
1. Point A is at the origin. Point B is 3 m distant from A at 30◦ coun-
terclockwise from the x axis. Point C is 2 m from point A at 100◦
counterclockwise from the x axis.
2. For the vectors in the previous problem, find D1 · D2 using both the
cosine form of the dot product and the Cartesian form. Check to see
if the two answers are the same.
ky
kb kc
1
ka
-1 1 2 kx
ω=4
-1 ω=8
ω = 12
8. Two gravity waves have the same frequency, but slightly different wave-
lengths.
(a) If one wave has an orientation angle θ = π/4 radians, what is the
orientation angle of the other? (See figure 1.6.)
(b) Determine the orientation of lines of constructive interference be-
tween these two waves.
11. A laser beam from a laser on the earth is bounced back to the earth
by a corner reflector on the moon.
(a) Engineers find that the returned signal is stronger if the laser beam
is initially spread out by the beam expander shown in figure 1.22.
Explain why this is so.
(b) The beam has a diameter of 1 m leaving the earth. How broad is
it when it reaches the moon, which is 4 × 105 km away? Assume
the wavelength of the light to be 5 × 10−7 m.
(c) How broad would the laser beam be at the moon if it weren’t
initially passed through the beam expander? Assume its initial
diameter to be 1 cm.
(a) How thick does the glass have to be to slow down the incoming
wave so that it lags the wave going through the other slit by a
phase difference of π? Take the wavelength of the light to be
λ = 6 × 10−7 m.
(b) For the above situation, describe qualitatively how the diffraction
pattern changes from the case in which there is no glass in front
of one of the slits. Explain your results.
beam spreader
laser
1m
(a) Qualitatively sketch the two slit diffraction pattern from this source.
Sketch the pattern for each wavelength separately.
(b) Qualitatively sketch the 16 slit diffraction pattern from this source,
where the slit spacing is the same as in the two slit case.