Fraunhofer diffraction: Difference between revisions

	Fraunhofer diffraction occurs when:
	– aperture or slit size, – wavelength, – distance from the aperture

Browse history interactively

← Previous edit Next edit →

Content deleted Content added

VisualWikitext

Inline

Revision as of 13:52, 3 January 2021

In optics, the Fraunhofer diffraction equation is used to model the diffraction of waves when the diffraction pattern is viewed at a long distance from the diffracting object (in the far-field region), and also when it is viewed at the focal plane of an imaging lens.^[1]^[2] In contrast, the diffraction pattern created near the object (in the near field region) is given by the Fresnel diffraction equation.

The equation was named in honor of Joseph von Fraunhofer although he was not actually involved in the development of the theory.^[3]

This article explains where the Fraunhofer equation can be applied, and shows the form of the Fraunhofer diffraction pattern for various apertures. A detailed mathematical treatment of Fraunhofer diffraction is given in Fraunhofer diffraction equation.

Equation

When a beam of light is partly blocked by an obstacle, some of the light is scattered around the object, light and dark bands are often seen at the edge of the shadow – this effect is known as diffraction.^[4] These effects can be modelled using the Huygens–Fresnel principle. Huygens postulated that every point on a primary wavefront acts as a source of spherical secondary wavelets and the sum of these secondary wavelets determines the form of the proceeding wave at any subsequent time. Fresnel developed an equation using the Huygens wavelets together with the principle of superposition of waves, which models these diffraction effects quite well.

It is not a straightforward matter to calculate the displacement (amplitude) given by the sum of the secondary wavelets, each of which has its own amplitude and phase, since this involves addition of many waves of varying phase and amplitude. When two waves are added together, the total displacement depends on both the amplitude and the phase of the individual waves: two waves of equal amplitude which are in phase give a displacement whose amplitude is double the individual wave amplitudes, while two waves which are in opposite phases give a zero displacement. Generally, a two-dimensional integral over complex variables has to be solved and in many cases, an analytic solution is not available.^[5]

The Fraunhofer diffraction equation is a simplified version of the Kirchhoff's diffraction formula and it can be used to model the light diffracted when both a light source and a viewing plane (the plane of observation) are effectively at infinity with respect to a diffracting aperture.^[6] With the sufficiently distant light source from the aperture, the incident light to the aperture is a plane wave so that the phase of the light at each point on the aperture is the same. The phase of the contributions of the individual wavelets in the aperture varies linearly with position in the aperture, making the calculation of the sum of the contributions relatively straightforward in many cases.

With a distant light source from the aperture, the Fraunhofer approximation can be used to model the diffracted pattern on a distant plane of observation from the aperture (far field). Practically it can be applied to the focal plane of a positive lens.

Far field

When the distance between the aperture and the plane of observation (on which the diffracted pattern is observed) is large enough so that the optical path lengths from edges of the aperture to a point of observation differ much less than the wavelength of the light, then propagation paths for individual wavelets from every point on the aperture to the point of observation can be treated as parallel. This is often known as the far field and is defined as being located at a distance which is significantly greater than $W 2 /λ$ , where $λ$ is the wavelength and $W$ is the largest dimension in the aperture. The Fraunhofer equation can be used to model the diffraction in this case.^[7]

For example, if a 0.5 mm diameter circular hole is illuminated by a laser with 0.6 μm wavelength, the Fraunhofer diffraction equation can be employed if the viewing distance is greater than 1000 mm.

Focal plane of a positive lens as the far field plane

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.^[8]

Examples of Fraunhofer diffraction

In each of these examples, the aperture is illuminated by a monochromatic plane wave at normal incidence.

Diffraction by a slit of infinite depth

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 $θ$ .^[9] 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:^[10]

\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

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

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.

Semi-quantitative analysis 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:

\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:

\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.

Single-slit diffraction of Electric Field using Huygens' Principle

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:^[11]

$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 $\beta =\omega /c=2\pi /\lambda$ . However $r_{1}=r-y\sin \theta$ and integrating from $-a/2$ to $a/2$ ,

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

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

Integrating we then get

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

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

$E={\frac {\sin(\psi ^{'}/2)}{\psi ^{'}/2}}$ ^[11]

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).^[12] 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.

Diffraction by a circular aperture

Computer simulation of the Airy diffraction pattern

The diffraction pattern given by a circular aperture is shown in the figure on the right.^[13] 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

\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.

Diffraction by 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.^[14] 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.^[15]

Diffraction by a double slit

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.^[16] These fringes are often known as Young's fringes.

The angular spacing of the fringes is given by

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

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

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.

Semi-quantitative explanation of double-slit 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

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

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

\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:

w=z\lambda /d.

Diffraction 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:^[18]

~\sin \theta _{n}=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 $θ 0$ , the grating equation is:

\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.

Semi-quantitative explanation

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:

W\sin \theta =n\lambda ,\quad n=0,\pm 1,\pm 2,\ldots

This is the same relationship that is given above.

References

^ Born & Wolf, 1999, p. 427.
^ Jenkins & White, 1957, p288
^ http://scienceworld.wolfram.com/biography/Fraunhofer.html
^ Heavens and Ditchburn, 1996, p. 62
^ Born & Wolf, 1999, p. 425
^ Jenkins & White, 1957, Section 15.1, p. 288
^ Lipson, Lipson and Lipson, 2011, p. 203
^ Hecht, 2002, p. 448
^ Hecht, 2002, Figures 10.6(b) and 10.7(e)
^ Jenkins & White, 1957, p. 297
^ ^a ^b Kraus, John Daniel; Marhefka, Ronald J. (2002). Antennas for all applications. McGraw-Hill. ISBN 9780072321036.
^ Born & Wolf, 1999, Figure 8.10
^ Born & Wolf, 1999, Figure 8.12
^ Hecht, 2002, Figure 11.33
^ Hecht, 2002, Figure 13.14
^ Born & Wolf, 1999, Figure 7.4
^ Hecht, 2002, eq. (9.30).
^ Longhurst, 1957, eq.(12.1)

^[1]

Sources

Born M & Wolf E, Principles of Optics, 1999, 7th Edition, Cambridge University Press, ISBN 978-0-521-64222-4
Heavens OS and Ditchburn W, Insight into Optics, 1991, Longman and Sons, Chichester ISBN 978-0-471-92769-3
Hecht Eugene, Optics, 2002, Addison Wesley, ISBN 0-321-18878-0
Jenkins FA & White HE, Fundamentals of Optics, 1957, 3rd Edition, McGraw Hill, New York
Lipson A., Lipson SG, Lipson H, Optical Physics, 4th ed., 2011, Cambridge University Press, ISBN 978-0-521-49345-1
Longhurst RS, Geometrical and Physical Optics,1967, 2nd Edition, Longmans, London

External links

^ Goodman, Joseph W. (1996). Introduction to Fourier Optics (second ed.). Singapore: The McGraw-HillCompanies, Inc. p. 73. ISBN 0-07-024254-2.

[1] Born & Wolf, 1999, p. 427.

[2] Jenkins & White, 1957, p288

[3] ttp://scienceworld.wolfram.com/biography/Fraunhofer.html

[4] Heavens and Ditchburn, 1996, p. 62

[5] Born & Wolf, 1999, p. 425

[6] Jenkins & White, 1957, Section 15.1, p. 288

[7] Lipson, Lipson and Lipson, 2011, p. 203

[8] Hecht, 2002, p. 448

[9] Hecht, 2002, Figures 10.6(b) and 10.7(e)

[10] Jenkins & White, 1957, p. 297

[:0-11] Kraus, John Daniel; Marhefka, Ronald J. (2002). Antennas for all applications. McGraw-Hill. ISBN 9780072321036.

[12] Born & Wolf, 1999, Figure 8.10

[13] Born & Wolf, 1999, Figure 8.12

[14] Hecht, 2002, Figure 11.33

[15] Hecht, 2002, Figure 13.14

[16] Born & Wolf, 1999, Figure 7.4

[17] Hecht, 2002, eq. (9.30).

[18] Longhurst, 1957, eq.(12.1)

[19] Goodman, Joseph W. (1996). Introduction to Fourier Optics (second ed.). Singapore: The McGraw-HillCompanies, Inc. p. 73. ISBN 0-07-024254-2.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

[16]

[17]

[18]

[1]

@@ Line 74: / Line 74: @@
 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 [[Electric field|field]] is:<ref name=":0">{{Cite book|url=https://books.google.com/books?id=NRxTAAAAMAAJ|title=Antennas for all applications|last1=Kraus|first1=John Daniel|last2=Marhefka|first2=Ronald J.|date=2002|publisher=McGraw-Hill|isbn=9780072321036|language=en}}</ref>
-<math>dE=\frac{A}{r_1}e^{j \omega [t-(r_1/c)]}dy=\frac{A}{r_1}e^{j(\omega t-\beta r_1)}dy</math>
+<math>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</math>
 where <math>\beta=\omega/c=2\pi /\lambda</math>. However <math>r_1=r-y\sin\theta</math>and integrating from <math>-a/2</math> to <math>a/2</math>,
-<math>E=A^'\int\limits_{-a/2}^{a/2} e^{j\beta y \sin\theta} dy</math>
+<math>E=A^'\int\limits_{-a/2}^{a/2} e^{i\beta y \sin\theta} dy</math>
-where <math>A^'=\frac{Ae^{j(\omega t-\beta r)}}{r_1}</math>.
+where <math>A^'=\frac{Ae^{i(\omega t-\beta r)}}{r_1}</math>.
 Integrating we then get