Codex @codex  0

The Airy disk is the central bright region, up to the first dark ring, in the scalar diffraction pattern of a uniformly illuminated, unobscured circular aperture. The whole pattern includes the weaker surrounding rings; its radius is distinct from its full width at half maximum.

For a circular optical pupil of diameter $D$ under normal, uniform monochromatic plane wave illumination, Fraunhofer diffraction givesHere $I_0$ is the on-axis light intensity, $\theta$ is the angle from the optical axis, $\lambda$ is the wavelength in the propagation medium, and $J_1$ is the order-one Bessel function of the first kind. The paraxial approximation replaces $\sin \theta$ by $\theta$, measured in radians. The first positive zero $j_{1,1}=3.831705970\dots$ givesThis is the angular radius of the Airy disk; the angular diameter is twice as large. The Airy disk derivation explains the origin of the Bessel function of the first kind and of the numerical coefficient.^[ref]

The image and both radial curves show the same calculated light intensity, divided by its peak. The image's logarithmic scale makes the weak rings visible; its colour is a display choice. The linear curve shows how small the ring peaks are compared with the centre. A dark ring marks a zero of wave amplitude, produced by cancellation across the optical pupil; it is not the edge of all transmitted light. The first Airy disk contains about $83.78\%$ of the total radiant flux, as calculated in Airy disk normalization and encircled energy.

For an original design example, take an ideal optical lens in air with focal length $f=100\mathrm{mm}$, circular optical pupil diameter $D=5.00\mathrm{mm}$ and wavelength $\lambda =550\mathrm{nm}$. In the focal plane, the radial coordinate $r\simeq f\theta$ givesThe first-zero diameter is $26.8327\mu \mathrm{m}$. Solving $P(u_{1/2})=1/2$ gives $u_{1/2}=1.616339948\dots$, so the diameter measured by full width at half maximum is onlyThus a quoted “spot diameter” is incomplete without its definition. Doubling $D$ halves both widths at fixed wavelength and focal length; increasing magnification alone changes the image scale without improving object-space discrimination.

The Airy disk supplies a calculable reference for a point spread function, allowing diffraction to be distinguished from defocus and optical aberrations. Its beauty is that one circular boundary determines an entire sequence of bright and dark rings through a Fourier transform. The ideal formula assumes an unobscured pupil, uniform illumination, ideal optical focus, a single wavelength, and a scalar description with small output angles. A blocked centre, nonuniform pupil transmission or a significant pupil wave phase error changes the response. High numerical aperture may require polarization of an electromagnetic wave and a vector field description. The circular-aperture diffraction laboratory experiment tests the far-field pattern directly; the Rayleigh criterion uses its first zero to define a two-source convention.

Write $a=D/2$ and wavenumber $k=2\pi /\lambda$. The pupil Fourier transform, normalized by pupil area, gives the Fraunhofer diffraction wave amplitude. In pupil polar coordinates,The angular average is $J_0$ because $J_0(z)=(2\pi {)}^{-1}{\int }_0^{2\pi }{e}^{-iz\cos \varphi }d\varphi$. The Bessel function of the first kind identity $d\{z{J}_1(z)\}/dz=z{J}_0(z)$ then reduces the radial integral to $A(u)=2J_1(u)/u$. Taking the squared magnitude yields the Airy disk light intensity profile.^[ref]

The apparent division by zero is removable. The power seriesgives $A(u)=1-u^2/8+u^4/192+O(u^6)$ and $P(u)=1-u^2/4+5u^4/192+O(u^6)$. At the centre, every pupil contribution has the same wave phase, so their normalized sum is one. The first dark ring instead corresponds to $J_1(u)=0$ with $u>0$, rather than to the removable zero of $J_1$ at the origin.

An actual teaching setup is described in Laboratory Exercise No. 4: Diffraction, MIT 6.237, Fall 2024, for MIT's 6.237 Modern Optics Project Laboratory, Fall 2024. Preparation includes predicting the aperture-size dependence and the distinction between Fraunhofer diffraction and Fresnel diffraction. A helium–neon laser supplies $632.8\mathrm{nm}$ light; the far-field setup uses a planar incident wavefront, a diffracting object and a distant screen, with beam expansion when needed. Figure 1 also marks an optional CCD camera. Exercise 4.7 asks students to sketch the circular-aperture pattern, measure the Airy disk diameter between correctly identified boundaries, and infer the aperture size. No measured values are supplied. As a regime comparison, exercise 4.9 varies the circular opening's position near a microscope objective optical lens, observes near-field patterns and compares them with the far-field case. This tests whether a change in propagation geometry changes the pattern. The coherent illumination measures one aperture's response, rather than the discrimination of two mutually incoherent optical point sources.^[ref]

For a quantitative interpretation, let $L$ be the aperture-to-screen distance and $d_{\mathrm{spot}}$ the first-zero diameter on the screen. In the paraxial approximation, the independent calculation isFor an explicitly hypothetical measurement, $L=1.80\mathrm{m}$ and $d_{\mathrm{spot}}=5.60\mathrm{mm}$ at $632.8\mathrm{nm}$ would imply $D\simeq 0.496\mathrm{mm}$. The scale $D^2/\lambda \simeq 0.389\mathrm{m}$ is smaller than $L$, but only by a factor of about $4.6$; this alone gives limited confidence in a high-precision Fraunhofer diffraction approximation. It is not a reported MIT measurement.

Additional recommended controls are to keep wavelength fixed, record $L$, centre and uniformly illuminate the aperture, and repeat with another $D$. Compare $d_{\mathrm{spot}}/L$, which should vary inversely with $D$. A camera readout should include blocked-beam background and avoid saturation. A scaled repetition at larger $L$ tests residual Fresnel diffraction; nonuniform illumination, optical aberrations, detector sampling and uncertainty in locating a dark ring limit the inferred $D$. The screen pattern's apparent bright edge depends on display brightness and is not a substitute for the first-zero boundary.

A diffraction grating is a periodic optical structure whose repeated slits or grooves direct light into diffraction orders through wave interference. Its beauty is that many individually broad contributions cooperate to select narrow directions: the spacing determines where the principal peaks occur, while the number of illuminated elements determines how sharply they are selected. This makes a diffraction grating useful in a spectrometer, where separating wavelengths spatially lets us examine an optical spectrum. ^[ref].

Consider a transmission diffraction grating in air, illuminated at normal incidence by a plane wave. Let $d$ be the centre-to-centre slit spacing, $\lambda$ the wavelength in the surrounding medium, and $\theta$ the signed outgoing angle from the forward normal. In Fraunhofer diffraction, where observation directions specify the far-field pattern, adjacent contributions have path difference $d\sin \theta$ and phase difference $\delta =2\pi d\sin \theta /\lambda$. Their wave amplitudes reinforce whenever $\delta =2\pi m$, giving $d\sin \theta =m\lambda$. This normal-incidence geometry should be kept distinct from the incident and reflected angle convention in the grating equation below. ^[ref].

Angular colour separation follows immediately. At fixed nonzero diffraction order, increasing the wavelength moves the peak farther from the normal. With $d=2\mu \mathrm{m}$, first-order light at $450\mathrm{nm}$ leaves at $13.003^{\circ }$, whereas $650\mathrm{nm}$ leaves at $18.966^{\circ }$: a separation of $5.963^{\circ }$. The zero diffraction order stays at $0^{\circ }$ for both wavelengths. Thus grating dispersion comes from an orderly dependence of wave interference on wavelength, rather than an arbitrary assignment of colours to directions.

For $N$ identical slits with uniform coherent illumination, meaning stable relative wave phases across the array, the normalized wave interference factor isThe apparent $0/0$ at a principal peak has limiting value one. A short explanation uses the geometric series: equal wave phases add $N$ equal complex amplitudes; a small directional change distributes those contributions around the complex plane and allows cancellation. With each slit's contribution held fixed, peak scalar wave intensity scales as $N^2$. The normalized plots remove that height change to make the narrowing visible. Stable relative wave phases, or optical coherence, are essential; a sum of independent intensities would lose this array structure. ^[ref].

For finite slit width $a$, the ideal rectangular-slit model multiplies $F_N$ by a broad single-slit diffraction envelope:The sinc function factor, extended continuously at $u=0$, controls the brightness available to the narrow peaks. For example, $a=0.5\mu \mathrm{m}$ gives envelope values about $0.811$, $0.405$, and $0.090$ at diffraction orders $1$, $2$, and $3$ in the $600\mathrm{nm}$ example. If an envelope zero coincides with an allowed diffraction order, that peak disappears. A permitted direction therefore need not produce a visible bright spot. ^[ref].

Historically, Joseph von Fraunhofer used fine parallel wires stretched between threaded rods to make diffraction gratings. Kenyon's historical apparatus account describes both this construction and a later reconstruction with visible light. The example makes the governing idea tangible: regular spacing can turn a simple array of obstacles into an instrument for distinguishing wavelengths. ^[ref].

Diffraction-limited resolution is a scale for distinguishing details with a specified ideal diffraction response and a specified criterion. For two equally bright, mutually incoherent optical point sources imaged through a uniformly transmitting, unobscured circular optical pupil, the Rayleigh criterion places the centre of one Airy disk at the first zero of the other. At small angles in air this separation isHere $D$ is pupil diameter, $\lambda$ is wavelength in air, and $j_{1,1}$ is the first positive zero of the order-one Bessel function of the first kind. This is the radius to the first zero, not the diameter of the central spot. Fraunhofer diffraction gives the radial peak-normalized light intensityThe removable value at zero follows from $J_1(q)\sim q/2$. The pupil Fourier transform reduces to a radial integral by circular symmetry, which produces the Bessel function of the first kind; its first zero supplies the factor $1.22$.^[ref]

In the usual scalar circular-pupil approximation for a fluorescence microscope, the corresponding object-plane Rayleigh criterion isHere ${\lambda }_0$ is the vacuum emission wavelength, and the object-side numerical aperture already includes the medium's refractive index. Do not divide ${\lambda }_0$ by $n$ again. For small cone angles in air, $\mathrm{NA}\simeq D/(2f)$ and image-plane radius $f{\theta }_{\mathrm{R}}$ agrees with this form. The assumptions include a single wavelength, ideal optical focus, and mutually incoherent emitters. High numerical aperture can require a vector description with polarization of an electromagnetic wave and emitter orientation; confocal microscopy has its own excitation-and-detection response. The same $0.61$ formula therefore does not describe every microscope or illumination method.^[ref]

In this original model, ${\lambda }_0=550\mathrm{nm}$ and numerical apertures $0.25$ and $0.50$ give first-zero radii $1.342\mu \mathrm{m}$ and $0.671\mu \mathrm{m}$. Each point spread function is divided by its own peak to compare shape. The pupil axes are transverse wave vector components multiplied by ${\lambda }_0/(2\pi )$, so pupil radius equals numerical aperture. A wider pupil accepts a wider range of transverse wave vectors, producing a narrower image. Increasing magnification alone enlarges the same spot and does not reduce its object-plane radius.

The Rayleigh criterion is a useful visual convention, not an absolute information barrier. For equal incoherent sources with total detected radiant flux $F$ and separation vector $\mathbf{s}$, the ideal image isBelow the Rayleigh criterion this still differs from $Fp(\mathbf{x})$; for small $\mathbf{s}$, symmetry cancels the linear change and the first change is quadratic in separation. Estimating two unknown source positions is then increasingly sensitive to noise, background, relative brightness and uncertainty in the point spread function. Localizing one isolated source and recovering an arbitrary extended object are different inference tasks. Measurements that separate optical modes can also extract separation information differently from direct image light intensity measurement.^[ref]

There remains a real limit for arbitrary images in this model: the optical transfer function is zero above the transmitted spatial frequency band. Its cutoff for scalar incoherent imaging is ${\nu }_c=2\mathrm{NA}/{\lambda }_0$ in object-plane cycles per unit length. The Nyquist–Shannon sampling theorem therefore suggests object-referred detector pixel pitch $p_{\mathrm{pix}}\le {\lambda }_0/(4\mathrm{NA})$, with finer sampling providing margin. Finite detector pixels also integrate over their area, suppressing contrast; coarse sampling can cause aliasing. Finer sampling cannot alone recover spatial frequencies removed by the optics, and Fourier deconvolution amplifies noise where the optical transfer function is small. Photon shot noise, read noise, optical aberrations and defocus can all worsen practical discrimination. A useful resolution claim therefore reports the object model, criterion, wavelength, numerical aperture, sampling and signal-to-noise ratio in photon counting.

The focal ratio, or f-number, of a converging optical instrument focused at infinity is the dimensionless quantitywhere $f>0$ is effective focal length and $D$ is the on-axis entrance pupil diameter. The notation f/4 means $N=4$, not $f=4$ in a particular unit. A $100\mathrm{mm}$ focal length with $D=25\mathrm{mm}$ gives f/4. Closing the aperture stop until $D=12.5\mathrm{mm}$ gives f/8 while retaining the same focal length. The mechanical opening need not equal $D$. ^[ref].

For a thin lens in air with a circular aperture stop at the thin lens, the converging bundle from a distant on-axis point has half-angle $\alpha$ satisfyingThe last step uses the small-angle approximation. For a compound optical lens, the cone is set by the image-side optical pupil and its distance from the receiving plane; the simple triangle is not an exact high-angle rule. At finite object distance a thin lens instead gives $v=f(1+|m|)$, by its imaging equation and transverse magnification, so its working cone corresponds to $v/D=N(1+|m|)$. Nominal focal ratio therefore does not alone specify close-range image irradiance. ^[ref].

Two useful scalings follow under stated conditions. For a fixed uniform incident beam and unchanged transmission, collected radiant flux scales as $D^2$. For an extended source with fixed radiance, the central image irradiance scales approximately as $N^{-2}$ in the paraxial approximation: the solid angle of the cone scales as $(D/f{)}^2$. For an ideal uniformly illuminated unobscured circular entrance pupil, the first dark ring of the point spread function has angular radius and focal plane radiusHere $\lambda$ is wavelength. These are first-zero radii, not diameters or a universal rule for resolving arbitrary objects. A proof outline evaluates the circular aperture's Fraunhofer diffraction using its Fourier transform, then uses the first zero of the resulting Bessel function. ^[ref].

At $\lambda =0.55\mu \mathrm{m}$, the example's f/4 and f/8 first-ring radii are $2.684$ and $5.368\mu \mathrm{m}$. Closing the aperture stop halves $D$, quarters collected radiant flux and doubles this diffraction radius. It can also reduce optical aberrations by excluding outer light rays; actual sharpness reflects both effects. This explains why a smaller aperture need not produce a sharper image. ^[ref].

A proposed optical bench check puts a fixed collimator before a converging thin lens and a photodetector large enough to collect the entire image behind it. Compare two centred aperture stops of known diameters while fixing source output, detector gain and integration time. The integrated signal tests the $D^2$ collection law; a spatially resolved optical camera at the focal plane can separately measure the point spread function. Repeat a dark reading and the first aperture stop, and check that neither the incoming beam nor detector clips either bundle. A resolved pinhole source, imperfect collimation, optical aberrations and detector pixel size can dominate the apparent spot width; a single peak reading does not measure total collected radiant flux.

Consider an ideal coherent monochromatic wave from an unresolved source on the optical axis, a uniform circular entrance pupil of diameter $D$, and ideal relay optics without wavefront errors. In the scalar Fraunhofer diffraction approximation, a Fourier transform carries the pupil's complex amplitude into the image plane. Choose dimensionless pupil coordinates $\mathbf{x}$ in units of $D$ and angular coordinates $\mathbf{u}=D\mathbf{\theta }/\lambda$, where $\lambda$ is the wavelength. DefineThe coronagraphic focal plane mask has transmission $M=1-W$, with $W=1$ inside the occulting disk and $W=0$ outside. Thus the field immediately after the mask is $MA$. Outside the disk $M=1$, so the local wave intensity is still $|A{|}^2$: hiding the core has left those rings intact. Reimaging an unlimited pupil giveswhere $\ast$ is convolution, $L$ is the Lyot stop transmission, and the proportionality constant in the irradiance is omitted. Coordinate reversals from physical relay optics are absorbed into the coordinate choice. These equations use wave amplitudes throughout propagation; applying a Fourier transform to wave intensity would give a different and incorrect model. This is the pupil/image calculation described by Sivaramakrishnan et al.; the circular disk and numerical choices below are our illustrative specialization. ^[ref].

The subtraction explains the beauty of the arrangement: in much of the interior, the original pupil field and the convolution term nearly cancel. Near a boundary, that cancellation is incomplete. The Lyot stop can therefore block a region containing disproportionately much unwanted light. A hard-edged coronagraphic focal plane mask has a transform with extended tails, so the interior need not be perfectly dark. With an unlimited relay, $L=1$ gives $A_{\mathrm{final}}=MA$ exactly: no further ring suppression occurs. A stop covering the full original optical pupil instead uses $L=P$ and already clips field outside the optical pupil; an undersized Lyot stop clips more of the remaining rim.

The toy uses a coronagraphic focal plane mask of angular radius $3\lambda /D$ and a Lyot stop of diameter $0.8D$. The discrete Fourier transform grid has $1024\times 1024$ samples, with 128 across $D$ and angular sampling $0.125\lambda /D$; it approximates the continuous equations on a finite periodic domain. The plot averages wave intensity in annuli of width $0.25\lambda /D$. Over the annulus $4\le |\mathbf{u}|\le 10$, the mean wave intensity divided by the direct stellar peak is approximately $1.40\times 10^{-4}$ without a coronagraph, the same with only the focal mask and an unlimited relay, $4.44\times 10^{-5}$ with the mask and full pupil, and $1.66\times 10^{-6}$ with the mask and $0.8D$ Lyot stop. Thus the undersized opening reduces this particular annular mean about 27-fold relative to the full-pupil comparison. It does not suppress every direction by that factor, and these values do not predict a real instrument's detection limit.

The ideal toy excludes a central obstruction, support vanes, scattering of light, wavefront errors, atmospheric effects, a spread of wavelengths, pointing drift and detector noise. It also contains only the star: no simulated planet or solar corona. For a far off-axis source whose image misses the coronagraphic focal plane mask, the continuous geometric coronagraph throughput of the $0.8D$ opening tends to $0.8^2=0.64$ relative to the full disk. Rejecting stellar light therefore costs some desired light as well. Enlarging the coronagraphic focal plane mask narrows the characteristic scale of its pupil-plane transform by the scaling property of the Fourier transform, but also hides a larger neighborhood of the source.

Fraunhofer diffraction is the far-field approximation for diffraction: contributions from an aperture reach a selected observation direction with approximately linear path-dependent wave phases. Adding their complex amplitudes gives a Fourier transform of the field across the aperture. It also describes the appropriate rear focal plane of an ideal optical lens under coherent plane-wave illumination. The diffraction grating examples use this regime. ^[ref]; ^[ref].

For a slit of width $a$, wavelength $\lambda$ and screen distance $L$, sufficient conditions near the forward direction areThe last two statements are equivalent: the important comparison is with $a^2/\lambda$, rather than just with the slit width. For a finite rectangular aperture, check the corresponding condition for both transverse dimensions. The single-slit diffraction calculation treats the transverse cut associated with width $a$. ^[ref].

A short explanation is that the distance from slit coordinate $x$ to an observation point at distance $R$ and angle $\theta$ isThe linear term supplies the phase difference; the quadratic term supplies a wave phase error of order $a^2/(\lambda R)$, since the wavenumber is $2\pi /\lambda$. Near the forward direction $R\approx L$, so the stated condition makes this error small. This is why a distant screen permits the simple single-slit diffraction formula. ^[ref].

Fresnel diffraction describes diffraction when the quadratic propagation wave phase across an aperture must be retained. Its Fresnel propagator becomes the Fraunhofer diffraction description when that curvature is negligible.

Single-slit diffraction is the spreading of light transmitted by one slit. Consider an ideal rectangular slit of width $a$, illuminated normally by a uniform plane wave with a single wavelength $\lambda$ in the surrounding medium. This monochromatic wave has uniform wave amplitude and equal incident wave phases across the slit, providing the required optical coherence. In the scalar Fraunhofer diffraction model, the normalized angular wave intensity isHere $\theta$ is measured from the forward normal in the plane across the slit width. The sinc function uses the unnormalized convention $\mathrm{sinc}u=\sin u/u$. The limit makes the central value one. This ideal envelope neglects additional direction-dependent propagation and detector factors; near the forward direction it describes the normalized screen irradiance. ^[ref].

The noncentral zeros satisfyThe restriction follows from the range of the sine of a real angle. The first minima have $\sin \theta =\pm \lambda /a$ when $\lambda /a\le 1$. If $a<\lambda$, this model has no noncentral zero in the forward hemisphere; equality puts the first zeros at grazing directions, outside the small-angle approximation and off any finite perpendicular screen. ^[ref].

For $0<\lambda /a<1$, define the central width as the separation between the two first minima. Its angular width and its width on a perpendicular screen areThe angular approximation uses radians and requires $\lambda /a\ll 1$; the screen approximation uses $y\approx L\theta$. The tangent relation is exact straight-line geometry, even though the underlying Fraunhofer diffraction field is approximate. This width is distinct from the full width at half maximum.

For $\lambda =500\mathrm{nm}$, $a=0.1\mathrm{mm}$ and $L=1\mathrm{m}$, the positive first minimum isThe dimensionless far-field parameter is $a^2/(\lambda L)=0.02$. Doubling the slit width to $a=0.2\mathrm{mm}$ gives $W\approx 5.00002\mathrm{mm}$ and parameter $0.08$; both curves below use the Fraunhofer diffraction approximation. Choosing $L=2\mathrm{m}$ halves those parameters and doubles both screen widths.

Each curve is normalized to its own $I_0$; the arrows measure between first minima. The calculation uses $\theta =\arctan (y/L)$ before evaluating the sinc function, rather than replacing $\sin \theta$ by $y/L$. The lower diagram shows the exact screen geometry with the vertical scale enlarged. Normalization removes the change in peak brightness when the slit width changes.

The wave interference explanation is especially simple: divide the slit into equal narrow strips and add their complex amplitudes, rather than their wave intensities. With wavenumber $k=2\pi /\lambda$, the continuous slit givesIntegrating the complex exponential gives the endpoint difference; its continuous value at zero is one. Taking the squared complex modulus gives the stated wave intensity. This is the Fourier transform of an interval indicator expressed in angular coordinates. Feynman's continuous-source argument obtains the same pattern by letting the number of sources tend to infinity. ^[ref].

At the positive first minimum, the wave phases span exactly $2\pi$ across the slit. Every strip in one half has a partner $a/2$ away, with equal wave amplitude and phase difference $\pi$, so the pairs cancel. At $\theta =0$, all the wave phases align instead. ^[ref].

The drawing uses $N=8$ equally weighted midpoint samples $x_j/a=-1/2+(j+1/2)/8$, with $j=0,\dots ,7$. At the first minimum their complex amplitudes are a common rotation of the eight roots of unity, so their head-to-tail polygon closes exactly. This finite sum illustrates cancellation; the continuous integral, rather than eight point sources, defines the slit pattern at general angles. The circle shows the continuous range of wave phases; coloured arrows show the eight samples.

A narrower aperture requires a larger angle to accumulate the cancelling phase difference, so narrowing the slit widens the beam. This is a useful limit on geometrical optics: confining light more tightly at an opening increases its subsequent angular spread. At fixed incident irradiance and slit height, reducing $a$ also reduces transmitted power, exposing a practical balance between diffraction-limited resolution and light throughput. In a diffraction grating, this single-slit diffraction envelope modulates the wave interference between different slits.