Fraunhofer diffraction

A geometrical diagram used to derive Fraunhofer condition at which Fraunhofer diffraction is valid.

The derivation of Fraunhofer condition here is based on the geometry described in the right box.[8] The diffracted wave path r₂ can be expressed in terms of another diffracted wave path r₁ and the distance b between two diffracting points by using the law of cosines;

$${\displaystyle {r_{2}}={\left(r_{1}^{2}+b^{2}-2b{r_{1}}\cos \left({\frac {\pi }{2}}+\theta \right)\right)}^{\frac {1}{2}}={r_{1}}{\left(1+{\frac {b^{2}}{r_{1}^{2}}}+2{\frac {b}{r_{1}}}\sin \theta \right)}^{\frac {1}{2}}.}$$

This can be expanded by using the binomial series for ${\displaystyle {\left(1+x\right)}^{\frac {1}{2}}}$ by letting ${\displaystyle x={\frac {b^{2}}{r_{1}^{2}}}+2{\frac {b}{r_{1}}}\sin \theta }$,

$${\displaystyle {r_{2}}={r_{1}}\left(1+{\frac {b}{r_{1}}}\sin \theta +{\frac {b^{2}}{2r_{1}^{2}}}\cos ^{2}\theta +\cdots \right)={r_{1}}+b\sin \theta +{\frac {b^{2}}{2r_{1}}}\cos ^{2}\theta +\cdots ~.}$$

The phase difference between waves propagating along the paths r₂ and r₁ are, with the wavenumber where λ is the light wavelength,

$${\displaystyle k{r_{2}}-k{r_{1}}=kb\sin \theta +k{\frac {b^{2}}{2r_{1}}}\cos ^{2}\theta +\cdots .}$$

If ${\displaystyle k{\frac {b^{2}}{2{r_{1}}}}\cos ^{2}\theta =\pi {\frac {b^{2}}{\lambda r_{1}}}\cos ^{2}\theta \ll \pi }$ so ${\displaystyle {\frac {b^{2}}{\lambda r_{1}}}\cos ^{2}\theta \ll 1}$, then the phase difference is ${\displaystyle kr_{2}-kr_{1}\approx kb\sin \theta }$. The geometrical implication from this expression is that the paths r₂ and r₁ are approximately parallel with each other. Since there can be a diffraction plane - observation plane diffracted wave path which angle with respect to a straight line parallel to the optical axis is close to 0, this approximation condition can be further simplified as ${\displaystyle {\frac {b^{2}}{\lambda }}\ll L}$ where L is the distance between two planes along the optical axis. Due to the fact that an incident wave on a diffracting plane is effectively a plane wave if ${\displaystyle {\frac {b^{2}}{\lambda }}\ll L}$ where L is the distance between the diffracting plane and the point wave source is satisfied, Fraunhofer condition is ${\displaystyle {\frac {b^{2}}{\lambda }}\ll L}$ where L is the smaller of the two distances, one is between the diffracting plane and the plane of observation and the other is between the diffracting plane and the point wave source.

Plane wave focused by a lens.

In the far field, propagation paths for wavelets from every point on an aperture to a point of observation are approximately parallel, and a positive lens (focusing lens) focuses parallel rays toward the lens to a point on the focal plane (the focus point position on the focal plane depends on the angle of the parallel rays with respect to the optical axis). So, if a positive lens with a sufficiently long focal length (so that differences between electric field orientations for wavelets can be ignored at the focus) is placed after an aperture, then the lens practically makes the Fraunhofer diffraction pattern of the aperture on its focal plane as the parallel rays meet each other at the focus.[9]

Graph and image of single-slit diffraction

The width of the slit is W. The Fraunhofer diffraction pattern is shown in the image together with a plot of the intensity vs. angle θ.[10] The pattern has maximum intensity at θ = 0, and a series of peaks of decreasing intensity. Most of the diffracted light falls between the first minima. The angle, α, subtended by these two minima is given by:[11]

$${\displaystyle \alpha \approx {\frac {2\lambda }{W}}}$$

Thus, the smaller the aperture, the larger the angle α subtended by the diffraction bands. The size of the central band at a distance z is given by

$${\displaystyle d_{f}={\frac {2\lambda z}{W}}}$$

Diffraction glass with 300 lines per millimeter

For example, when a slit of width 0.5 mm is illuminated by light of wavelength 0.6 μm, and viewed at a distance of 1000 mm, the width of the central band in the diffraction pattern is 2.4 mm.

The fringes extend to infinity in the y direction since the slit and illumination also extend to infinity.

If W < λ, the intensity of the diffracted light does not fall to zero, and if D << λ, the diffracted wave is cylindrical.

Geometry of single-slit diffraction

We can find the angle at which a first minimum is obtained in the diffracted light by the following reasoning. Consider the light diffracted at an angle θ where the distance CD is equal to the wavelength of the illuminating light. The width of the slit is the distance AC. The component of the wavelet emitted from the point A which is travelling in the θ direction is in anti-phase with the wave from the point B at middle of the slit, so that the net contribution at the angle θ from these two waves is zero. The same applies to the points just below A and B, and so on. Therefore, the amplitude of the total wave travelling in the direction θ is zero. We have:

$${\displaystyle \theta _{\text{min}}\approx {\frac {CD}{AC}}={\frac {\lambda }{W}}.}$$

The angle subtended by the first minima on either side of the centre is then, as above:

$${\displaystyle \alpha =2\theta _{\text{min}}={\frac {2\lambda }{W}}.}$$

There is no such simple argument to enable us to find the maxima of the diffraction pattern.

Continuous broadside array of point sources of length a.

We can develop an expression for the far field of a continuous array of point sources of uniform amplitude and of the same phase. Let the array of length a be parallel to the y axis with its center at the origin as indicated in the figure to the right. Then the differential field is:[12]

$${\displaystyle dE={\frac {A}{r_{1}}}e^{i\omega [t-(r_{1}/c)]}dy={\frac {A}{r_{1}}}e^{i(\omega t-\beta r_{1})}dy}$$

where ${\displaystyle \beta =\omega /c=2\pi /\lambda }$. However ${\displaystyle r_{1}=r-y\sin \theta }$ and integrating from ${\displaystyle -a/2}$ to ${\displaystyle a/2}$,

$${\displaystyle E\simeq A'\int _{-a/2}^{a/2}e^{i\beta y\sin \theta }dy}$$

where ${\displaystyle A'={\frac {Ae^{i(\omega t-\beta r)}}{r}}}$.

Integrating we then get

$${\displaystyle E={\frac {2A'}{\beta \sin \theta }}\sin \left({\frac {\beta a}{2}}\sin \theta \right)}$$

Letting ${\displaystyle \psi ^{'}=\beta a\sin \theta =\alpha _{r}\sin \theta }$ where the array length in radians is ${\displaystyle a_{r}=\beta a=2\pi a/\lambda }$, then,[12]

$${\displaystyle E=A'a{\frac {\sin(\psi ^{'}/2)}{\psi ^{'}/2}}}$$

Computer simulation of Fraunhofer diffraction by a rectangular aperture

The form of the diffraction pattern given by a rectangular aperture is shown in the figure on the right (or above, in tablet format).[13] There is a central semi-rectangular peak, with a series of horizontal and vertical fringes. The dimensions of the central band are related to the dimensions of the slit by the same relationship as for a single slit so that the larger dimension in the diffracted image corresponds to the smaller dimension in the slit. The spacing of the fringes is also inversely proportional to the slit dimension.

If the illuminating beam does not illuminate the whole vertical length of the slit, the spacing of the vertical fringes is determined by the dimensions of the illuminating beam. Close examination of the double-slit diffraction pattern below shows that there are very fine horizontal diffraction fringes above and below the main spot, as well as the more obvious horizontal fringes.

Computer simulation of the Airy diffraction pattern

The diffraction pattern given by a circular aperture is shown in the figure on the right.[14] This is known as the Airy diffraction pattern. It can be seen that most of the light is in the central disk. The angle subtended by this disk, known as the Airy disk, is

$${\displaystyle \alpha \approx {\frac {1.22\lambda }{W}}}$$

where W is the diameter of the aperture.

The Airy disk can be an important parameter in limiting the ability of an imaging system to resolve closely located objects.

Intensity of a plane wave diffracted through an aperture with a Gaussian profile

The diffraction pattern obtained given by an aperture with a Gaussian profile, for example, a photographic slide whose transmissivity has a Gaussian variation is also a Gaussian function. The form of the function is plotted on the right (above, for a tablet), and it can be seen that, unlike the diffraction patterns produced by rectangular or circular apertures, it has no secondary rings.[15] This technique can be used in a process called apodization—the aperture is covered by a Gaussian filter, giving a diffraction pattern with no secondary rings.

The output profile of a single mode laser beam may have a Gaussian intensity profile and the diffraction equation can be used to show that it maintains that profile however far away it propagates from the source.[16]

Double-slit fringes with sodium light illumination

In the double-slit experiment, the two slits are illuminated by a single light beam. If the width of the slits is small enough (less than the wavelength of the light), the slits diffract the light into cylindrical waves. These two cylindrical wavefronts are superimposed, and the amplitude, and therefore the intensity, at any point in the combined wavefronts depends on both the magnitude and the phase of the two wavefronts.[17] These fringes are often known as Young's fringes.

The angular spacing of the fringes is given by

$${\displaystyle \theta _{\text{f}}=\lambda /d.}$$

The spacing of the fringes at a distance z from the slits is given by[18]

$${\displaystyle w_{\text{f}}=z\theta _{f}=z\lambda /d,}$$

where d is the separation of the slits.

The fringes in the picture were obtained using the yellow light from a sodium light (wavelength = 589 nm), with slits separated by 0.25 mm, and projected directly onto the image plane of a digital camera.

Double-slit interference fringes can be observed by cutting two slits in a piece of card, illuminating with a laser pointer, and observing the diffracted light at a distance of 1 m. If the slit separation is 0.5 mm, and the wavelength of the laser is 600 nm, then the spacing of the fringes viewed at a distance of 1 m would be 1.2 mm.

Geometry for far-field fringes

The difference in phase between the two waves is determined by the difference in the distance travelled by the two waves.

If the viewing distance is large compared with the separation of the slits (the far field), the phase difference can be found using the geometry shown in the figure. The path difference between two waves travelling at an angle θ is given by

$${\displaystyle d\sin \theta \approx d\theta .}$$

When the two waves are in phase, i.e. the path difference is equal to an integral number of wavelengths, the summed amplitude, and therefore the summed intensity is maximal, and when they are in anti-phase, i.e. the path difference is equal to half a wavelength, one and a half wavelengths, etc., then the two waves cancel, and the summed intensity is zero. This effect is known as interference.

The interference fringe maxima occur at angles

$${\displaystyle d\theta _{n}=n\lambda ,\quad n=0,\pm 1,\pm 2,\ldots }$$

where λ is the wavelength of the light. The angular spacing of the fringes is given by

$${\displaystyle \theta _{\text{f}}\approx \lambda /d.}$$

When the distance between the slits and the viewing plane is z, the spacing of the fringes is equal to zθ and is the same as above:

$${\displaystyle w=z\lambda /d.}$$

Diffraction of a laser beam by a grating

A grating is defined in Born and Wolf as "any arrangement which imposes on an incident wave a periodic variation of amplitude or phase, or both".

A grating whose elements are separated by S diffracts a normally incident beam of light into a set of beams, at angles θ_n given by:[19]

$${\displaystyle ~\sin \theta _{n}={\frac {n\lambda }{S}},\quad n=0,\pm 1,\pm 2,\ldots }$$

This is known as the grating equation. The finer the grating spacing, the greater the angular separation of the diffracted beams.

If the light is incident at an angle θ₀, the grating equation is:

$${\displaystyle \sin \theta _{n}={\frac {n\lambda }{S}}+\sin \theta _{0},\quad n=0,\pm 1,\pm 2,\ldots }$$

The detailed structure of the repeating pattern determines the form of the individual diffracted beams, as well as their relative intensity while the grating spacing always determines the angles of the diffracted beams.

The image on the right shows a laser beam diffracted by a grating into n = 0, and ±1 beams. The angles of the first order beams are about 20°; if we assume the wavelength of the laser beam is 600 nm, we can infer that the grating spacing is about 1.8 μm.

A simple grating consists of a series of slits in a screen. If the light travelling at an angle θ from each slit has a path difference of one wavelength with respect to the adjacent slit, all these waves will add together, so that the maximum intensity of the diffracted light is obtained when:

$${\displaystyle W\sin \theta =n\lambda ,\quad n=0,\pm 1,\pm 2,\ldots }$$

This is the same relationship that is given above.