Mathematical Description of Single-Slit Diffraction and Circular Aperture Diffraction

When geometric optics proves insufficient to account for phenomena such as light interference, polarization, and diffraction, the Huygens-Fresnel principle offers a microscopic resolution. This fundamental principle postulates that every point on a wavefront acts as a source of secondary spherical wavelets, and the optical disturbance at any subsequent point is the coherent superposition of these wavelets. Applying this principle to specific aperture geometries—specifically single-slit and circular configurations—yields the core mathematical formulations that govern light intensity distribution. This article explores the derivation, key parameters, and physical implications of these two classical diffraction models.
Consider a narrow slit of width $a$, illuminated by a monochromatic plane wave of wavelength $\lambda$ propagating perpendicularly to the slit plane. According to the Huygens-Fresnel principle, the optical intensity at a screen point $P$ (defined by the diffraction angle $\theta$) is determined by the interference of all secondary wavelets originating across the aperture.

To simplify the mathematical treatment under the Fraunhofer diffraction regime, we assume the distance from the aperture to the observation screen is vastly larger than the slit dimensions. This allows us to neglect higher-order phase variations caused by path-length differences. Consequently, the phase difference $\Delta \phi$ between wavelets originating from position $x$ across the slit is directly proportional to the path difference $\delta = x \sin\theta$.

Integrating the contributions of all infinitesimal wavelets across the aperture yields the well-known intensity distribution function:
$$ I(\theta) = I_0 \left( \frac{\sin \beta}{\beta} \right)^2 $$
where $\beta = \frac{\pi a \sin\theta}{\lambda}$ represents the normalized phase factor.

This equation reveals the hallmark signatures of single-slit diffraction: a central maximum that is both the widest and most intense ($I_0$), flanked by symmetric secondary fringes of rapidly diminishing intensity. Minima occur when the phase factor satisfies $\beta = n\pi$ ($n = \pm 1, \pm 2, \dots$), leading to the condition for destructive interference:
$$ a \sin\theta_n = n\lambda $$
This relationship demonstrates that longer wavelengths or smaller apertures amplify the diffraction angle, rendering the wave nature of light more pronounced.

Mathematical Description of Circular Aperture Diffraction: The Airy Disk and Bessel Functions

Transitioning from a linear obstruction to a circular aperture shifts the problem from a one-dimensional integral to a two-dimensional domain. For a circular opening of diameter $D$ (radius $r_0 = D/2$), the resulting diffraction pattern exhibits rotational symmetry, manifesting as concentric bright and dark rings rather than linear fringes. Due to this complex boundary geometry, the trigonometric functions utilized in single-slit analysis are replaced by Bessel functions.

Within the Fraunhofer approximation, the intensity distribution generated by a circular aperture is governed by the Bessel function of the first kind:
$$ I(\theta) = I_0 \left[ \frac{J_1(k r_0 \sin\theta)}{k r_0 \sin\theta / 2} \right]^2 $$
In engineering and applied optics, this is frequently expressed using the dimensionless variable $u = \frac{\pi D \sin\theta}{\lambda}$, simplifying the intensity profile to:
$$ I(\theta) = I_0 \left( \frac{J_1(u)}{u/2} \right)^2 $$

The defining feature of circular aperture diffraction is the central bright disc known as the Airy disk, which establishes the theoretical resolution limit for optical instruments. The first dark ring corresponds to the first root of the Bessel function where $J_1(u) = 0$, occurring approximately at $u = 3.8317$. This yields the primary angular resolution criterion for circular apertures:
$$ \sin\theta_{min} = 1.22 \frac{\lambda}{D} $$

Comparative Physical Significance and Applications

While the mathematical frameworks of single-slit and circular aperture diffraction differ in dimensionality, both originate from the same foundation of Huygens-Fresnel superposition. Their distinct mathematical profiles stem entirely from geometrical symmetry: an infinitely long one-dimensional slit constrains phase accumulation to a single axis (producing a sinc-squared distribution), whereas a two-dimensional circular boundary necessitates radial integration, naturally generating Bessel function solutions.

In practical optical systems, circular diffraction fundamentally constrains imaging resolution. According to Rayleigh's criterion, two point sources are considered just resolvable when the center of the Airy disk of one source coincides with the first minimum of the other. This mathematical constraint directly dictates optical design and aperture selection: increasing the aperture diameter $D$ compresses the diffraction angle $\theta_{min}$, thereby enhancing resolution. Conversely, over-constricting the aperture amplifies diffraction blurring, posing an intrinsic limit independent of aberrations like defocus or spherical aberration.

Furthermore, the single-slit model serves as a vital analytical tool for understanding diffraction gratings and scanning slit systems, whereas the circular aperture model is indispensable for evaluating microscopes, telescopes, and photographic lenses. Mastery of the mathematical derivations and physical principles underlying both models provides an essential foundation for exploring wave optics and photon propagation, enabling engineers and physicists to accurately predict optical performance and optimize experimental design.