# Diffraction

Diffraction is the bending and spreading of waves when they meet an obstruction. It can occur with any type of wave, including sound waves, water waves, and electromagnetic waves such as light and radio waves. Diffraction also occurs when any group of waves of a finite size is propagating; for example, a narrow beam of light waves from a laser must, because of diffraction of the beam, eventually diverge into a wider beam at a sufficient distance from the laser. As a simple example of diffraction, if you speak into one end of a cardboard tube, the sound waves emerging from the other end spread out in all directions, rather than propagating in a straight line like a stream of water from a garden hose.

## Introduction

Diffraction is one particular type of wave interference, caused by the partial obstruction or lateral restriction of a wave. Not all interference is diffraction; for example, sound waves emitted by two stereo speakers will interfere with each other if they are of the same frequency and have a definite phase relationship, but this is not diffraction. Diffraction will not occur if the wave is not coherent, and diffraction effects become weaker (and ultimately undetectable) as the size of obstruction is made larger and larger compared to the wavelength. In well-defined cases, a diffraction pattern may be observed.

Diffraction is not the same as refraction, although both are phenomena in which a wave does not propagate in a single direction. Refraction is not an interference phenomenon, and, e.g., can occur without coherence.

It is the diffraction of "particles," such as electrons, which stood as one of the powerful arguments in favor of quantum mechanics. It is possible to observe diffraction of particles such as neutrons or electrons and hence we are able to infer the existence of wave-particle duality. Indeed, this diffraction is a useful tool; the wavelengths of these particle-waves are small enough that they are used as probes of the atomic structure of crystals. See electron diffraction and neutron diffraction.

File:Doubleslitdiffraction.png
Double-slit diffraction

File:Laserdiffraction.jpg
Double-slit diffraction
(red laser light)

File:Diffraction2vs5.jpg
2-slit and 5-slit diffraction

The most conceptually simple example of diffraction is double-slit diffraction in which both slits have relatively narrow widths compared to the wavelength of the wave. Suppose, for the sake of visualization, that these are water waves. After passing through the slits, two overlapping patterns of semicircular ripples are formed, as shown in the first figure. Where a crest overlaps with a crest, a double-height crest will be formed; this is constructive interference. Constructive interference also occurs where a trough overlaps another trough. However, when a trough and a crest overlap, they cancel out; the interference is destructive. The second figure shows the result of this process with light waves of a single wavelength originating from a laser. The constructive-interference locations are called maxima, because they have maximum brightness. The destructive-interference locations are the minima. Historically, the first proof that light was a wave phenomenon came from the double-slit experiment of Thomas Young.

Several qualitative observations can be made:

• When the dimensions of the diffracting object are reduced, the angular spacing of the diffraction pattern is increased in inverse proportion. (More precisely, this is true of the sines of the angles.)
• The diffraction angles are invariant under scaling; that is, they depend only on the ratio of the wavelength to a dimension, a, of the diffracting object.
• When the diffracting object is repeated, the effect is to narrow each maximum, concentrating its energy within a narrower range of angles. The third figure, for example, shows a comparison of a double-slit pattern with a pattern formed by five slits, both sets of slits having the same spacing, a, between the center of one slit and the next.

## Mathematical description

It is mathematically easier to consider the case of far-field or Fraunhofer diffraction, where the diffracting obstruction is many wavelengths distant from the point at which the wave is measured. The more general case is known as near-field or Fresnel diffraction, and involves more complex mathematics. As the observation distance is increased the results predicted by the Fresnel theory converge towards those predicted by the simpler Fraunhofer theory. This article considers far-field diffraction, which is commonly observed in nature.

Quantitatively, the angular positions of the minima in multiple-slit diffraction are given by the equation

$\sin \theta = \frac{\lambda}{a} m,$

where m is an integer that labels the order of each minimum. The central maximum is two orders wide, however, so m = 0, θ = 0 is the absolute maximum of the distribution and intensity functions. This is a form of Bragg's law (see below).

File:Diffraction1.png
Graph and image

### Quantitative analysis of single-slit diffraction

As an example, we will now derive an exact equation for the intensity of the diffraction pattern as a function of angle in the case of single-slit diffraction.

Consider a monochromatic complex plane wave $\Psi^\prime$ of wavelength λ incident on a slit of width a.

Let the slit lie in the x′-y′ plane, with its center at the origin.

Then we shall assume that diffraction generates a complex wave ψ, traveling radially in the r direction away from the slit, and this is given by:

$\Psi = \int_{slit} \frac{i}{r\lambda} \Psi^\prime e^{-ikr}\,dslit$

let (x′,y′,0) be a point inside the slit over which we are integrating; and let (x,0,z) be the location at which we are computing the intensity of the diffraction pattern. The slit extends from $x^\prime=-a/2$ to $+a/2\,$, and from $y'=-\infty$ to $\infty$.

The distance r from the slot is:

$r = \sqrt{\left(x - x^\prime\right)^2 + y^{\prime2} + z^2}$
$r = z \left(1 + \frac{\left(x - x^\prime\right)^2 + y^{\prime2}}{z^2}\right)^\frac{1}{2}$

We assume Fraunhofer diffraction, so that $z >> \big|\left(x - x^\prime\right)\big|$. In other words, the distance to the target is much larger than the diffraction width on the target. By the binomial expansion rule, ignoring terms quadratic and higher, we can estimate our quantity on the right to be:

$r \approx z \left( 1 + \frac{1}{2} \frac{\left(x - x^\prime \right)^2 + y^{\prime 2}}{z^2} \right)$
$r \approx z + \frac{\left(x - x^\prime\right)^2 + y^{\prime 2}}{2z}$

We see that our 1/r in front of the equation is non-oscillatory, i.e. its contribution to the magnitude of the intensity is small compared to our exponential factors. Therefore, we will lose little accuracy by approximating it as z.

 $\Psi \,$ $= \frac{i \Psi^\prime}{z \lambda} \int_{-\frac{a}{2}}^{\frac{a}{2}}\int_{-\infty}^{\infty} e^{-ik\left[z+\frac{ \left(x - x^\prime \right)^2 + y^{\prime 2}}{2z}\right]} \,dx^\prime \,dy^\prime$ $= \frac{i \Psi^\prime}{z \lambda} e^{-ikz} \int_{-\frac{a}{2}}^{\frac{a}{2}}e^{-ik\left[\frac{\left(x - x^\prime \right)^2}{2z}\right]} \,dx^\prime \int_{-\infty}^{\infty} e^{-ik\left[\frac{y^{\prime 2}}{2z}\right]} \,dy^\prime$ $=\Psi^\prime \sqrt{\frac{i}{z\lambda}} e^\frac{-ikx^2}{2z} \int_{-\frac{a}{2}}^{\frac{a}{2}}e^\frac{ikxx^\prime}{z} e^\frac{-ikx^{\prime 2}}{2z} \,dx^\prime$

To make things cleaner, we will use a placeholder 'C' to denote constants in our equation. It is important to keep in mind that C can contain imaginary numbers, thus our wave function will be complex, however at the end, we will bracket our ψ, which will eliminate any imaginary components.

Now, in Fraunhoffer diffraction, $kx^{\prime 2}/z$ is small, so $e^\frac{-ikx^{\prime 2}}{2z} \approx 1$. The same approximation holds for $e^\frac{-ikx^2}{2z}$. Thus, taking $C = \Psi^\prime \sqrt{\frac{i}{z\lambda}}$, we have:

 $\Psi\,$ $= C \int_{-\frac{a}{2}}^{\frac{a}{2}}e^\frac{ikxx^\prime}{z} \,dx^\prime$ $=C \frac{\left(e^\frac{ikax}{2z} - e^\frac{-ikax}{2z}\right)}{\frac{2ikax}{2z}}$

Now we note that $\sin x = \frac{e^{ix} - e^{-ix}}{2i}$ and $\sin \theta = \frac{x}{z}$.

$\Psi = C \frac{\sin\frac{ka\sin\theta}{2}}{\frac{ika\sin\theta}{2}}$

Now, substituting in $\frac{2\pi}{\lambda} = k$, the intensity I of the diffracted waves at an angle θ is given by:

 $I(\theta)\,$ $= \frac{c}{8\pi}\big|\Psi(\theta)\big|^2$ $= I_0 \langle \Psi \Big|\Psi \rangle \,$ $= I_0 {\left[ \operatorname{sinc} \left( \frac{\pi a}{\lambda} \sin \theta \right) \right] }^2$

where the sinc function is given by sinc(x) = sin(x)/x.

### Quantitative analysis of n-slit diffraction

$\Psi = \int_{slit} \frac{i}{r\lambda} \Psi^\prime e^{-ikr}\,dslit$

Consider n slits in the prime plane of the equal size (a, $\infty$, 0) and spacing d spread along the x′ axis. As above, the distance r from the slit 1 is:

$r = z \left(1 + \frac{\left(x - x^\prime\right)^2 + y^{\prime2}}{z^2}\right)^\frac{1}{2}$

To generalize this to n slits, we make the observation that while z and y remain constant, x′ shifts by

$x_{j=0 \cdots n-1}^{\prime} = x_0^\prime - j d$

Thus

$r_j = z \left(1 + \frac{\left(x - x^\prime - j d \right)^2 + y^{\prime2}}{z^2}\right)^\frac{1}{2}$

and the sum of all n contributions to the wave function is:

$\Psi = \sum_{j=0}^{N-1} C \int_{-\frac{a}{2}}^{\frac{a}{2}} e^\frac{ikx\left(x^\prime - jd\right)}{z} e^\frac{-ik\left(x^\prime - jd\right)^2}{2z} \,dx^\prime$

Again noting that $\frac{k\left(x^\prime -jd\right)^2}{z}$ is small, so $e^\frac{-ik\left(x^\prime -jd\right)^2}{2z} \approx 1$, we have:

 $\Psi\,$ $= C\sum_{j=0}^{N-1} \int_{-\frac{a}{2}}^{\frac{a}{2}} e^\frac{ikx\left(x^\prime - jd\right)}{z} \,dx^\prime$ $= C \sum_{j=0}^{N-1} \frac{\left(e^{\frac{ikax}{2z} - \frac{ijkxd}{z}} - e^{\frac{-ikax}{2z}-\frac{ijkxd}{z}}\right)}{\frac{2ikax}{2z}}$ $= C \sum_{j=0}^{N-1} e^\frac{ijkxd}{z} \frac{\left(e^\frac{ikax}{2z} - e^\frac{-ikax}{2z}\right)}{\frac{2ikax}{2z}}$ $= C \frac{\sin\frac{ka\sin\theta}{2}}{\frac{ka\sin\theta}{2}} \sum_{j=1}^{N-1} e^{ijkd\sin\theta}$

Now, we can use the following identity

$\sum_{j=0}^{N-1} e^{x j} = \frac{1 - e^{Nx}}{1 - e^x}.$

Substituting into our equation, we find:

 $\Psi\,$ $= C \frac{\sin\frac{ka\sin\theta}{2}}{\frac{ka\sin\theta}{2}}\left(\frac{1 - e^{iNkd\sin\theta}}{1 - e^{ikd\sin\theta}}\right)$ $= C \frac{\sin\frac{ka\sin\theta}{2}}{\frac{ka\sin\theta}{2}}\left(\frac{e^{-iNkd\frac{\sin\theta}{2}}-e^{iNkd\frac{\sin\theta}{2}}}{e^{-ikd\frac{\sin\theta}{2}}-e^{ikd\frac{\sin\theta}{2}}}\right)\left(\frac{e^{iNkd\frac{\sin\theta}{2}}}{e^{ikd\frac{\sin\theta}{2}}}\right)$ $= C \frac{\sin\frac{ka\sin\theta}{2}}{\frac{ka\sin\theta}{2}}\frac{\frac{e^{-iNkd \frac{\sin\theta}{2}} - e^{iNkd\frac{\sin\theta}{2}}}{2i}}{\frac{e^{-ikd\frac{\sin\theta}{2}} - e^{ikd\frac{\sin\theta}{2}}}{2i}} \left(e^{i(N-1)kd\frac{\sin\theta}{2}}\right)$ $= C \frac{\sin\left(\frac{ka\sin\theta}{2}\right)}{\frac{ka\sin\theta}{2}} \frac{\sin\left(\frac{Nkd\sin\theta}{2}\right)} {\sin\left(\frac{kd\sin\theta}{2}\right)}e^{i\left(N-1\right)kd\frac{\sin\theta}{2}}$

We now make our k substitution as before and represent all non-oscillating constants by the I0 variable as in the 1-slit diffraction and bracket the result. Remember that

$\langle e^{ix} \Big| e^{ix}\rangle\ = e^0 = 1$

This allows us to discard the tailing exponent and we have our answer:

$I\left(\theta\right)\,$
$=I_0 \left[ \operatorname{sinc} \left( \frac{\pi a}{\lambda} \sin \theta \right) \right]^2$
$\left[\frac{\sin\left(\frac{N\pi d}{\lambda}\sin\theta\right)}{\sin\left(\frac{\pi d}{\lambda}\sin\theta\right)}\right]^2$

## Other cases

### Bragg diffraction

Diffraction from multiple slits, as described above, is similar to what occurs when waves are scattered from a periodic structure, such as atoms in a crystal or rulings on a diffraction grating. Each scattering center (e.g., each atom) acts as a point source of spherical wavefronts; these wavefronts undergo constructive interference to form a number of diffracted beams. The direction of these beams is described by Bragg's law:

$m \lambda = 2 d \sin \theta,\,$

where λ is the wavelength, d is the distance between scattering centers, θ is the angle of diffraction and m is an integer known as the order of the diffracted beam. Bragg diffraction is used in X-ray crystallography to deduce the structure of a crystal from the angles at which X-rays are diffracted from it. Since the diffraction angle θ is dependent on the wavelength λ, diffaction gratings impart angular dispersion on a beam of light.

The most common demonstration of Bragg diffraction is the spectrum of colors seen reflected from a compact disc: the closely-spaced tracks on the surface of the disc form a diffraction grating, and the individual wavelengths of white light are diffracted at different angles from it, in accordance with Bragg's law.

### Diffraction limit of telescopes

File:Zboo lucky image 1pc.png
The Airy disc around each of the stars from the 2.56m telescope aperture can be seen in this lucky image of the binary star zeta Boötis.

For diffraction through a circular aperture, there is a series of concentric rings surrounding a central Airy disc. The mathematical result is similar to a radially symmetric version of the equation given above in the case of single-slit diffraction.

A wave does not have to pass through an aperture to diffract; for example, a beam of light of a finite size also undergoes diffraction and spreads in diameter. This effect limits the minimum size d of spot of light formed at the focus of a lens, known as the diffraction limit:

$d = 2.44 \lambda \frac{f}{a},\,$

where λ is the wavelength of the light, f is the focal length of the lens, and a is the diameter of the beam of light, or (if the beam is filling the lens) the diameter of the lens. (See Rayleigh criterion).

By use of Huygens' principle, it is possible to compute the diffraction pattern of a wave from any arbitrarily shaped aperture. If the pattern is observed at a sufficient distance from the aperture, it will appear as the two-dimensional Fourier transform of the function representing the aperture.