Fraunhofer diffraction
Fraunhofer diffraction is a type of wave interference pattern in optics that occurs when coherent light waves pass through an aperture or scatter off an obstacle, with the light source and observation screen positioned at effectively infinite distances from the diffracting element, resulting in plane wavefronts and a far-field approximation.[1][2] This condition leads to a linear phase variation across the aperture and uniform intensity contributions from each point, producing characteristic patterns such as bright and dark fringes that scale uniformly with distance.[3] Named after the German physicist and optician Joseph von Fraunhofer (1787–1826), who pioneered its study through experiments with diffraction gratings in 1821–1822, the phenomenon laid foundational work for spectroscopy by enabling precise wavelength measurements of light.[4][5] In contrast to near-field Fresnel diffraction, which involves curved wavefronts and distance-dependent patterns, Fraunhofer diffraction simplifies analysis by treating incoming and outgoing waves as parallel, often achieved practically using lenses to focus the pattern at a finite distance.[4][3] Mathematically, the intensity distribution for a single slit, for example, follows I(\theta) = I_0 \left( \frac{\sin \alpha}{\alpha} \right)^2, where \alpha relates to the slit width, diffraction angle \theta, and wavelength \lambda, while multi-slit or grating patterns exhibit principal maxima at angles satisfying \sin \theta = m \lambda / d (with m as the order and d as the spacing).[3][4] This far-field regime underpins applications in Fourier optics, where the diffraction pattern represents the Fourier transform of the aperture function, facilitating techniques in holography, beam analysis, and modern imaging systems like laser-based spectrometers.[4] Fraunhofer's innovations, building on earlier wave theories by Huygens and Fresnel, shifted optics from geometric to physical principles, influencing electromagnetic wave understanding advanced by Maxwell in the 1870s.[4][5]Fundamentals
Definition and conditions
Fraunhofer diffraction refers to the far-field regime of wave diffraction, in which a planar incident wavefront illuminates an aperture or diffracting object, and the resulting diffraction pattern is observed in a plane at effectively infinite distance from the object, yielding an intensity distribution that depends solely on the observation angle rather than on specific positions in the observation plane.[2][6] This contrasts with near-field diffraction, known as Fresnel diffraction, which involves curved wavefronts and position-dependent patterns closer to the diffracting object.[2] The Fraunhofer condition requires that the distances from the source to the diffracting object and from the object to the observation plane be much larger than the square of the object's characteristic dimension divided by the wavelength of the light, mathematically expressed as z \gg \frac{a^2}{\lambda}, where z is the propagation distance, a is the aperture size, and \lambda is the wavelength.[6] This ensures that the incoming and outgoing wavefronts can be treated as planar across the extent of the aperture, simplifying the diffraction analysis to a Fourier transform relationship between the aperture transmittance and the far-field amplitude.[2] The theory assumes illumination by a monochromatic, coherent plane wave, which maintains phase coherence across the wavefront and allows for straightforward superposition of contributions from different parts of the aperture.[2] Under these conditions, the diffracted intensity in the far field is given by I(\theta) = |U(\theta)|^2, where U(\theta) is the complex amplitude as a function of the angular coordinate \theta.[6] This diffraction regime is named after Joseph von Fraunhofer (1787–1826), who pioneered its study through experiments with diffraction gratings starting in 1821, enabling precise wavelength measurements including the dark absorption lines observed earlier.[7][2]Comparison with Fresnel diffraction
Fresnel diffraction describes the near-field propagation of light from an aperture, where the observation plane is at a finite distance, leading to curved wavefronts and the inclusion of quadratic phase factors in the wave amplitude. These phase factors, arising from the variation in path lengths, are typically expressed as \exp\left(i \frac{\pi x^2}{\lambda z}\right), where x is the transverse coordinate, \lambda is the wavelength, and z is the propagation distance.[8] This results in diffraction patterns that depend on the specific distance z, producing interference fringes whose size and structure vary with position, such as bright and dark bands near shadow edges.[4] In contrast, Fraunhofer diffraction emerges as a limiting case of Fresnel diffraction when the propagation distance z becomes sufficiently large, satisfying z \gg \frac{a^2}{\lambda} (where a is the aperture dimension), such that the quadratic phase factor approximates to 1, simplifying the wavefront to effectively plane.[9] Under this far-field condition, the phase variation across the aperture is linear, and the diffraction pattern corresponds to the Fourier transform of the aperture function, independent of z in terms of angular distribution.[10] Practically, Fraunhofer patterns exhibit scaling of their angular width proportional to \lambda / a, unaffected by distance, whereas Fresnel patterns show distance-dependent fringe spacing that changes as the observation plane moves.[8] Fraunhofer diffraction is particularly suited for analyzing angular spectra in far-field optics, such as in spectrometers or telescope resolution limits, where plane-wave assumptions hold.[9] Fresnel diffraction, however, applies to near-field scenarios like edge diffraction in shadows or computational holography, where spherical wavefront curvature must be accounted for to model finite-distance effects accurately.[4]| Aspect | Fresnel Diffraction | Fraunhofer Diffraction |
|---|---|---|
| Field Region | Near field (z \lesssim a^2 / \lambda) | Far field (z \gg a^2 / \lambda) |
| Wavefront Assumption | Spherical (curved) | Plane |
| Phase Factors | Includes quadratic terms, e.g., \exp(i \pi x^2 / \lambda z) | Linear phase only; quadratic terms negligible |
| Pattern Dependence | On distance z; fringes scale with z | Angular distribution independent of z; scales with \lambda / a |
| Typical Applications | Shadow edges, near-field holography | Angular spectra, far-field imaging |
Mathematical formulation
Fraunhofer approximation
The Fraunhofer approximation describes the diffracted field in the far-field regime, where the observation distance z greatly exceeds both the wavelength \lambda and the square of the aperture dimensions divided by \lambda, ensuring plane wave-like propagation from the aperture. This approximation simplifies the general diffraction integral by neglecting quadratic phase terms in the propagation kernel, resulting in a direct relationship to the spatial Fourier transform of the aperture function. The Fraunhofer condition is met when the incident wave can be treated as plane, which is valid in the far field or when using a lens to focus the pattern. The diffraction amplitude U(\theta_x, \theta_y) at observation angles \theta_x and \theta_y (small angles in the paraxial regime) is given by U(\theta_x, \theta_y) = \frac{i}{\lambda z} \iint_{\text{aperture}} t(x,y) \exp\left[-i \frac{2\pi}{\lambda} (u x + v y)\right] \, dx \, dy, where t(x,y) is the complex aperture transmission function, and u = x'/z \approx \theta_x, v = y'/z \approx \theta_y represent the spatial coordinates scaled by the propagation distance, effectively serving as angular variables.[12][13] The corresponding intensity pattern is then I(\theta_x, \theta_y) = |U(\theta_x, \theta_y)|^2, which provides the observable diffraction distribution independent of the overall phase factors.[12] This formulation reveals that the Fraunhofer diffraction amplitude is proportional to the two-dimensional Fourier transform of the aperture function t(x,y), evaluated at spatial frequencies f_x = u / \lambda = \theta_x / \lambda and f_y = v / \lambda = \theta_y / \lambda. The exponential term \exp\left[-i 2\pi (f_x x + f_y y)\right] directly embodies this transform, highlighting the duality between spatial coordinates in the aperture plane and frequency coordinates in the observation plane.[12][13] Regarding normalization and units, the prefactor i / (\lambda z) ensures dimensional consistency for the electric field amplitude: if t(x,y) is dimensionless (as for amplitude transmission), the integral over area yields units of length squared, and division by \lambda z (length squared) results in a field strength with units of volts per meter, assuming consistent SI conventions for the incident field. Often, the expression is normalized such that the on-axis amplitude (\theta_x = \theta_y = 0) relates to the total aperture area for energy conservation.[12] The validity of the Fraunhofer approximation relies on the paraxial approximation, limiting it to small diffraction angles where \theta_x, \theta_y \ll 1, and the scalar wave assumption, which treats the field as a scalar quantity neglecting vectorial polarization effects. These conditions hold under monochromatic illumination and when evanescent waves are negligible, typically for observation distances satisfying z \gg D^2 / \lambda, where D is the aperture diameter.[12][13]Derivation from wave optics
The derivation of Fraunhofer diffraction from wave optics relies on the Huygens-Fresnel principle, which posits that every point on a wavefront acts as a source of secondary spherical wavelets, with the resulting field at a distant point determined by superposition.[14] Under Kirchhoff's boundary conditions—where the field and its normal derivative in the aperture match the incident wave, and both vanish behind an opaque screen—the principle leads to the Fresnel-Kirchhoff diffraction formula for the scalar field U(P) at observation point P: U(P) = \frac{1}{i \lambda} \iint_{\Sigma} \frac{U_0(x', y')}{r} \exp(i k r) \cos \chi \, dx' dy', where \lambda is the wavelength, k = 2\pi / \lambda, U_0(x', y') is the incident field over the aperture \Sigma, r is the distance from aperture element (x', y', 0) to P(x, y, z), and \cos \chi is the obliquity factor accounting for the angle \chi between the aperture normal and the line to P.[15] This formula assumes a scalar electromagnetic field and monochromatic plane-wave illumination, treating secondary sources as isotropic within the paraxial approximation.[14] In the Fraunhofer regime, the observation point lies in the far field, where the distance z satisfies z \gg a^2 / \lambda (with a the aperture dimension) to ensure quadratic phase terms remain small. The distance r is then approximated by expanding the exact form r = \sqrt{(x - x')^2 + (y - y')^2 + z^2} and neglecting higher-order terms beyond the linear approximation: r \approx z - (x x' + y y') / z.[15] This binomial expansion is valid because the phase variation k (x'^2 + y'^2)/(2z) \ll \pi across the aperture, linearizing the propagation paths. Similarly, the amplitude factor simplifies to $1/r \approx 1/z, reflecting the geometric dilution of intensity.[14] The phase term follows directly: \exp(i k r) \approx \exp(i k z) \exp\left[ -i k (x x' + y y') / z \right], isolating a common propagating phase \exp(i k z) outside the integral and leaving a bilinear phase within that depends on aperture and observation coordinates. For small observation angles (paraxial condition), the obliquity factor \cos \chi \approx 1, as the direction to P aligns nearly parallel to the aperture normal, eliminating directional bias in secondary wave contributions.[15] Substituting these into the Fresnel-Kirchhoff formula yields the Fraunhofer diffraction integral: U(P) = \frac{\exp(i k z)}{i \lambda z} \iint_{\Sigma} U_0(x', y') \exp\left[ -i k (x x' + y y') / z \right] dx' dy', which represents the Fourier transform of the aperture field U_0, with spatial frequencies f_x = x / (\lambda z) and f_y = y / (\lambda z).[14] In a lens-based setup, where the observation plane is at the focal length f of a converging lens (effectively setting z = f), the prefactor scales as $1/(i \lambda f), and the transform coordinates adjust to f_x = x / (\lambda f), f_y = y / (\lambda f), preserving the far-field pattern in the focal plane.[15] These approximations hold under the scalar field assumption, neglecting vectorial effects like polarization, and treat secondary sources as isotropic point radiators modulated only by the incident field, consistent with the Huygens-Fresnel principle for isotropic media.[14]Experimental observation
Lens-based far-field setup
In a lens-based far-field setup for Fraunhofer diffraction, a coherent light source, typically a laser, is first spatially filtered (e.g., using a pinhole) and collimated—often with a separate collimating lens—to produce a plane wave that uniformly illuminates an aperture or diffracting object. The aperture is positioned in the collimated beam, ideally at the front focal plane of a subsequent positive (converging) lens (or immediately in front of the lens as a practical approximation). The lens collects the diffracted wavefront emanating from the aperture and focuses it onto a screen or detector placed at the lens's focal plane, a distance f behind the lens, where the Fraunhofer diffraction pattern is observed as an intensity distribution.[16] The converging lens performs an optical Fourier transform on the aperture's transmission function, directing parallel rays (corresponding to specific diffraction angles) to distinct points in the focal plane.[17] This focal plane acts as the effective Fraunhofer observation plane, where the angular spectrum of the diffracted light is spatially mapped to linear coordinates via the small-angle approximation \theta \approx x'/f, with \theta representing the diffraction angle and x' the transverse position on the screen.[18] The schematic of the setup consists of the light source and spatial filtering to achieve coherence, followed by a collimating lens producing a plane wave, the aperture at (or near) the front focal plane of the converging lens, the converging lens itself, and the observation screen precisely at the rear focal plane.[19] This arrangement effectively meets the Fraunhofer condition by optically simulating infinite propagation distance. One key advantage is the compactness of the setup, which fits within a standard laboratory bench (often with f on the order of 10–50 cm), avoiding the impracticality of physically extending the observation distance to infinity as in the pure far-field case.[20] In terms of intensity, the Fraunhofer pattern in this lens configuration scales with the focal length f rather than the propagation distance z, such that the on-axis intensity is proportional to $1/f^2, reflecting the lens's focusing action and providing a practical means to control pattern size by varying f.[18][17] This adjustment facilitates quantitative measurements, as the pattern's linear scale directly ties to the lens properties, enabling straightforward comparison with theoretical predictions.Infinite distance approximation
The infinite distance approximation in Fraunhofer diffraction refers to the scenario where the light source is effectively at infinite distance from the diffracting aperture, producing plane wavefronts, and the observation screen is also positioned at a sufficiently large distance to capture the far-field pattern. This setup satisfies the Fraunhofer conditions without optical elements like lenses, relying instead on geometric separation where the distance z to the screen greatly exceeds z \gg D^2 / \lambda, with D as the aperture diameter and \lambda the wavelength.[21] Such conditions are naturally met in astronomical observations, where stellar or solar light acts as a point source at infinity, illuminating the aperture (e.g., a telescope's entrance pupil) with parallel rays, and the diffraction pattern forms in the far field.[22] Historically, Joseph von Fraunhofer utilized this approximation in his early 19th-century experiments with telescopes to analyze diffraction effects in solar spectra, leveraging the sun's remote position to observe far-field patterns from gratings and apertures without needing finite-distance corrections.[23] In laboratory settings, a direct far-field realization involves placing the aperture and screen at large separations, such as across an extended bench or using mirrors to fold the path, ensuring the source-to-aperture and aperture-to-screen distances both satisfy \rho, r \gg D^2 / 4\lambda.[24] To approximate plane wave illumination without an infinite source, collimated laser beams are commonly employed, as their low divergence mimics wavefronts from a distant point source. A helium-neon laser, often expanded to uniformly illuminate the aperture, enables Fraunhofer patterns to be observed at practical distances (e.g., several meters) while maintaining the far-field criterion.[24] This method avoids lenses by directly projecting the beam onto the aperture, though precise alignment of the beam path is essential to prevent distortions from slight divergence.[21] Computational simulations provide a theoretical means to evaluate the infinite distance approximation by numerically computing the Fraunhofer diffraction integral, which represents the far-field limit as a Fourier transform of the aperture function, bypassing physical distance constraints entirely. These simulations, implemented in tools like MATLAB, allow precise prediction of patterns for arbitrary apertures at effectively infinite z, facilitating analysis where experimental setups are impractical. Practical limitations arise in real implementations: in astronomical contexts, atmospheric turbulence introduces phase fluctuations that broaden the diffraction pattern beyond the ideal Airy disk, imposing a "seeing" limit typically larger than the diffraction limit for ground-based observations.[22] In controlled environments, challenges include maintaining alignment over long paths and achieving uniform illumination, often requiring vibration isolation and extended optical tables, though these can be mitigated by adaptive techniques or shorter effective paths via mirrors.[24] This lens-free approach contrasts with focal-plane methods but demands greater spatial resources for fidelity.Diffraction patterns from basic apertures
Rectangular slit
In Fraunhofer diffraction, a rectangular slit is modeled as a one-dimensional aperture with uniform transmission, defined by the aperture function t(x) = \rect\left(\frac{x}{a}\right), where a is the slit width and the rect function equals 1 for |x| < a/2 and 0 otherwise. This setup assumes the slit is infinitely extended in the perpendicular direction, focusing on diffraction in one dimension. The far-field amplitude distribution U(\theta) is proportional to the Fourier transform of the aperture function, yielding U(\theta) \propto a \sinc\left( \frac{\pi a \sin\theta}{\lambda} \right), where \lambda is the wavelength and \theta is the observation angle from the slit normal.[25] The resulting intensity pattern is given by I(\theta) = I_0 \left[ \frac{\sin \beta}{\beta} \right]^2, where I_0 is the central intensity and \beta = \frac{\pi a \sin\theta}{\lambda}.[25] This sinc-squared profile features a bright central maximum at \theta = 0, flanked by secondary maxima of decreasing intensity and minima where the argument \beta = m\pi for integer m = \pm 1, \pm 2, \dots. The first minima occur at \sin\theta = \pm \frac{\lambda}{a}, or approximately \theta \approx \pm \frac{\lambda}{a} for small angles, defining the angular half-width of the central maximum as \frac{\lambda}{a} and the full width as \frac{2\lambda}{a}.[25] Qualitatively, the pattern arises from the Huygens-Fresnel principle, treating each point along the slit as a secondary spherical wave source with equal amplitude under plane-wave illumination. At the center (\theta = 0), path lengths from all points are equal, leading to constructive interference and maximum intensity. Away from the center, path differences accumulate, causing phase variations that result in constructive peaks at secondary maxima and complete destructive interference at the minima, where the slit divides into an integer number of equal segments with opposing phases.[26] In multi-slit configurations, the single-slit diffraction pattern provides the overall intensity envelope, modulating the finer interference fringes from the slit array and determining their visibility.[18]Circular aperture
Fraunhofer diffraction through a circular aperture of diameter D produces a radially symmetric intensity pattern due to the aperture's transmission function t(r) = \circ(r/D), where the circ function \circ is 1 for r < D/2 and 0 otherwise. This setup models the ideal case of a plane wave incident on a hard-edged circular opening, leading to interference in the far field that forms concentric rings. The pattern's symmetry arises from the aperture's rotational invariance, distinguishing it from one-dimensional slit diffraction. The complex amplitude in the observation plane is obtained via the Hankel transform (the radial form of the Fourier transform) of the aperture function. For small angles, this yields a distribution proportional to J_1(\pi D \sin \theta / \lambda) / (\pi D \sin \theta / \lambda), where J_1 is the first-order Bessel function of the first kind, \theta is the angular coordinate from the optical axis, and \lambda is the wavelength. This expression was first derived by George Biddell Airy to describe the image of a point source through a telescope objective. The Airy pattern relates to the Fourier transform property of Fraunhofer diffraction, where the far-field pattern represents the Fourier transform of the aperture transmittance. The intensity distribution, or Airy pattern, is then given by I(\theta) = I_0 \left[ \frac{2 J_1(\alpha)}{\alpha} \right]^2, where \alpha = \pi D \sin \theta / \lambda and I_0 is the central intensity. This squared form arises from the modulus of the amplitude, producing a central bright spot (Airy disk) surrounded by alternating dark and bright concentric rings of diminishing intensity. The first dark ring, marking the boundary of the Airy disk, occurs where J_1(\alpha) = 0 at \alpha \approx 3.832, corresponding to an angular radius \theta \approx 1.22 \lambda / D for small \theta. Lord Rayleigh established this as the criterion for resolving two point sources, where their Airy disks just touch, setting the diffraction limit for optical resolution. The two-dimensional isotropic nature of the pattern manifests as circular rings, with the central Airy disk enclosing approximately 84% of the total energy in the diffraction pattern. This energy concentration highlights the pattern's efficiency in focusing light, though the rings contain the remaining 16%. In applications such as telescopes, the $1.22 \lambda / D limit directly constrains the smallest resolvable angular separation, fundamental to astronomical imaging and the design of large-aperture instruments.Gaussian aperture
The Gaussian aperture features a smooth, radially symmetric transmittance function for the electric field, t(x,y) = \exp\left[ -\frac{x^2 + y^2}{2 \sigma^2} \right], where \sigma characterizes the spatial extent of the aperture, with the intensity profile decaying as \exp\left[ -(x^2 + y^2)/\sigma^2 \right]. This soft-edged form arises naturally in scenarios involving Gaussian beam illumination, avoiding the discontinuities inherent in hard-edged apertures.[27] Under the Fraunhofer approximation, the far-field amplitude distribution is the Fourier transform of the aperture transmittance, yielding a self-similar Gaussian profile: U(\theta_x, \theta_y) \propto \exp\left[ -\frac{2\pi^2 \sigma^2 (\theta_x^2 + \theta_y^2)}{\lambda^2} \right], where \lambda is the wavelength and \theta_x, \theta_y are the angular coordinates. The resulting intensity pattern is I(\theta) \propto \exp\left[ -\frac{4\pi^2 \sigma^2 \theta^2}{\lambda^2} \right], with \theta = \sqrt{\theta_x^2 + \theta_y^2}, producing a rotationally symmetric distribution that narrows with increasing \sigma.[27] This intensity exhibits a smooth exponential decay without any zeros, sidelobes, or concentric rings, in stark contrast to the oscillatory patterns from abrupt-edged apertures. The lack of ripples stems from the analytic continuity of the Gaussian function, ensuring minimal energy scattering into secondary maxima.[27] Such properties render the Gaussian aperture ideal for laser beam propagation, where the self-similar diffraction preserves beam quality over long distances without unwanted artifacts that could degrade focusing or coupling efficiency.[27]Interference in multiple apertures
Double slit
The double-slit configuration in Fraunhofer diffraction involves two narrow, parallel slits of finite width a, separated by a center-to-center distance d, illuminated by a plane wave of monochromatic light with wavelength \lambda.[28] In the far-field approximation, the diffraction pattern arises from the superposition of wavelets emanating from all points across both slits, resulting in an interference pattern modulated by the single-slit diffraction envelope.[28] The path difference between waves from corresponding points in the two slits determines the phase shift, leading to constructive interference when this difference is an integer multiple m\lambda (where m = 0, \pm 1, \pm 2, \ldots), or d \sin \theta = m \lambda, with \theta the angle from the central axis.[28] The resulting amplitude in the far field is proportional to the single-slit amplitude modulated by the two-slit interference factor:U(\theta) \propto \operatorname{sinc}(\beta) \cos\left(\frac{\delta}{2}\right),
where \beta = \frac{\pi a \sin \theta}{\lambda} accounts for the phase variation across each slit width, and \delta = \frac{2\pi d \sin \theta}{\lambda} is the phase difference between the slits.[28] The intensity pattern is then the square of this amplitude:
I(\theta) = I_0 \operatorname{sinc}^2(\beta) \cos^2\left(\frac{\delta}{2}\right),
with I_0 the intensity at the central maximum.[28] The \cos^2(\delta/2) term produces rapid oscillations of interference fringes, while the \operatorname{sinc}^2(\beta) envelope from the individual slits broadens the overall pattern and suppresses higher-order fringes.[28] The angular spacing between adjacent interference maxima is \Delta \theta \approx \lambda / d for small angles, where the approximation \sin \theta \approx \theta holds.[28] This spacing decreases as d increases, leading to finer fringes, but the single-slit envelope limits the visible orders, with missing fringes occurring when an interference maximum coincides with an envelope minimum (e.g., when a \sin \theta = m' \lambda for integer m' \neq 0).[28] In the limiting case where the slit width a \to 0, the single-slit envelope becomes uniform, reducing the pattern to the pure two-source interference of Young's double-slit experiment, with equally spaced fringes extending indefinitely.[28]