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.

Peak normalization $P(0)=1$ does not make the integral of a point spread function equal to one. In the paraxial approximation to the focal plane, define $b=\pi D/(\lambda f)$ and $u=br$. The unit-integral point spread function normalization isHere $p$ has units of inverse area; multiplying by the detected radiant flux gives image-plane irradiance. Its encircled energy, meaning the fraction inside radius $R=U/b$, isThis is an integral over area, so its radial weight is $u$. An unweighted integral of a line through the centre is a different quantity.

A short proof uses the Bessel function of the first kind identities $J_0^{\prime }=-J_1$ and $J_1^{\prime }=J_0-J_1/u$. They giveBoth expressions for $E$ vanish at zero. Since these Bessel functions of the first kind tend to zero as $U\to \infty$, $E(\infty )=1$, establishing the normalization. At the first dark ring,A circular detector region ending at the first zero therefore misses about $16.22\%$ of the ideal radiant flux, even though the surrounding rings have low peak light intensity. This matters when converting measured counts into total source brightness or choosing an image aperture for a photodetector.

Beam divergence describes the angular spread of a beam away from its optical axis; specify the radius convention and whether the quoted angle is a half-angle or full angle. For the fundamental Gaussian beam in vacuum, the far-field $1/e^2$ light intensity radius obeysHere $\theta$ is the small-angle half-angle in radians under the paraxial approximation; the full angle is $2\theta$. The geometric asymptote has angle $\arctan (w_0/z_R)$, approximated by $w_0/z_R$ in this regime. Other beam profiles need their own beam divergence definitions and propagation laws.

The beam waist is the narrowest cross-section of a freely propagating Gaussian beam. Its radius $w_0$ uses the wave amplitude $1/e$, or light intensity $1/e^2$, convention; its diameter is $2w_0$. Specifying a “spot size” without saying radius or diameter can introduce a factor-of-two error. Within the scalar paraxial approximation, the wavefront at the beam waist is planar. The minimum radius sets both the Rayleigh range and the far-field beam divergence.

In coherent imaging, complex wave amplitudes add with their relative wave phases preserved. The measured light intensity includes wave interference between the contributions.

For an isolated fluorescent molecule at displacement $\mathbf{r}$ from the scanning optical focus, write $E(\mathbf{r})$ for relative excitation light intensity and $D_a(\mathbf{r})$ for the fraction of its collected emission passed by a centred detection pinhole of radius $a$. The expected response, apart from a common brightness and collection factor, isHere $\mathbf{u}$ is position in the detection-pinhole plane and $p_{\mathrm{em}}$ is the emitted light intensity distribution there conditional on collection, with units of inverse area. This point spread function normalization retains throughput: $D_a$ is a fraction, not a curve rescaled to have peak one. The product expresses excitation followed by detection; it is an intensity/probability product, not multiplication of coherent wave amplitudes from different molecules. ^[ref]

Assume linear, unsaturated fluorescence, independent emitters whose light intensities add, a stationary specimen, fixed collection efficiency and quantum efficiency, and negligible changes from scattering, absorption or bleaching during a scan. Different excitation and emission wavelengths generally give different excitation and detection responses. For a sufficiently small pinhole, $D_a$ is approximately its area times $p_{\mathrm{em}}$ at the centre; this is the familiar excitation-times-detection point spread function. It is not universally the square of one point spread function. A spatially invariant response further allows the image to be expressed as a convolution with the emitter distribution; depth-dependent optical aberrations can invalidate that approximation.

For an explicitly synthetic axial example, put $t=z/z_0$ and choose an on-axis excitation $E(t)=1/(1+t^2)$, inspired by a Gaussian beam. Stipulate a circular Gaussian function for the detection-plane light intensity, with standard deviation $s(t)=s_0\sqrt{1+t^2}$. In dimensionless coordinates $\mathbf{v}=\mathbf{u}/s_0$, put $\alpha =a/s_0$. ThenThe expression for $D_{\alpha }$ follows by integrating $p$ over a disk using polar coordinates. The assumptions specify a toy model with matched axial scales; its parameters are not a calibration for a real objective optical lens. It omits Airy disk rings, vector diffraction, and unequal excitation and emission wavelengths.

For a vanishingly small pinhole, the peak-normalized synthetic response tends to $(1+t^2{)}^{-2}$, whose full width at half maximum is $2\sqrt{\sqrt{2}-1}{z}_0\simeq 1.2872z_0$. With all collected emission accepted it tends to $(1+t^2{)}^{-1}$, of width $2z_0$. Meanwhile $H_{\alpha }(0)=1-e^{-{\alpha }^2/2}\to 0$ as the pinhole closes. The narrower normalized curve therefore accompanies a vanishing detected signal; peak normalization must not conceal that loss.

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.

A diffraction-limited system has an optical response close to the ideal diffraction response for its specified optical pupil, wavelength and optical focus. Improving the optical aberrations further makes little difference to that response. This describes the optics: a poorly sampled or noisy detector can still deliver an image with much worse usable diffraction-limited resolution.

The optical pupil determines which wave amplitudes and wave phases contribute to an image point. In a scalar, paraxial approximation at the focal plane of an optical lens in air, write the pupil transmission as $T(u,v)e^{i\varphi (u,v)}$, after removing the ideal focusing phase. The amplitude spread function is proportional to its scaled Fourier transform:Here $(u,v)$ are pupil coordinates, $(x,y)$ are image coordinates, $f$ is focal length, and $\lambda$ is the wavelength in air. For a uniformly transmitting, unobscured circular optical pupil at ideal optical focus, the result is an Airy disk with surrounding rings. A blocked centre, support structures, or nonuniform transmission changes the ideal point spread function; “diffraction limited” does not specify one universal spot shape. The existing Fourier model of a Lyot coronagraph illustrates how pupil and focal masks deliberately alter this response.^[ref]

Optical aberrations add pupil wave phase errors. Defocus contributes a quadratic pupil wave phase in the paraxial approximation, redistributing light intensity from the best-focused spot. The Strehl ratio compares the peak with the ideal response for the same pupil and total radiant flux. A broad optical spectrum produces a weighted sum of wavelength-dependent point spread functions; chromatic aberration can additionally change optical focus with wavelength. Thus the operating conditions belong with any claim that a system is diffraction limited.

For an elliptically polarized wave, the tip of its transverse vector traces an ellipse at a fixed point during each wave period. For polarization of an electromagnetic wave, this is the electric field vector. With the real components in polarization of waves, assume $A_x,A_y>0$ and $\sin \delta \ne 0$. Eliminating time givesA short derivation uses $\varphi =\omega t$ and the cosine addition formula to obtain $E_y/A_y=(E_x/A_x)\cos \delta -\sin \varphi \sin \delta$; squaring and using the sine and cosine identity ${\sin }^2\varphi =1-{\cos }^2\varphi$ gives the displayed equation. Outside phase quadrature, the mixed term gives a tilted ellipse; $A_x,A_y$ are component wave amplitudes and need not be its semi-major and semi-minor axes.

In phase quadrature, $\delta =\pm \pi /2$, the mixed term vanishes. For $(A_x,A_y)=(2,1)V/m$, the ellipse is $E_x^2/4+E_y^2=1$, its semi-major and semi-minor axes have lengths $2$ and $1V/m$, and its electric field magnitude varies between $1$ and $2V/m$. At $\varphi =\pi /4$ with $\delta =\pi /2$, the electric field is $(\sqrt{2},-1/\sqrt{2})V/m$ and has magnitude $\sqrt{5/2}V/m$. Equal wave amplitudes instead give circular polarization, a special case of elliptical polarization with constant magnitude. As $\delta$ approaches $0$ or $\pi$, the ellipse collapses to a line segment and gives linear polarization; this limit is degenerate rather than a nondegenerate ellipse.

Changing $\delta$ from $\pi /2$ to $-\pi /2$ reverses traversal of the same ellipse. In the displayed $x,y$ plot, $\delta =\pi /2$ moves from $(A_x,0)$ toward negative $y$ as time increases. A label such as right- or left-handed circular polarization additionally requires a stated viewing direction and naming convention: viewing along propagation and viewing against propagation reverse the apparent rotation. ^[ref].

The ellipse records both the relative wave amplitudes and the phase difference of two transverse oscillations. This separates a change of polarization of waves from a change of overall wave amplitude: multiplying both components by the same positive factor enlarges the ellipse while preserving its shape and orientation.

A rotating ideal linear polarizer can probe this geometry. For the $2:1$ ellipse in phase quadrature, orthogonal projection onto a pass axis at angle $\theta$ gives $E_{\mathrm{out}}(t)=2\cos \theta \cos \omega t-\sin \theta \sin \omega t$, in $V/m$. Averaging its square over a wave period gives $(4{\cos }^2\theta +{\sin }^2\theta )/2$, while the incoming squared electric field magnitude averages to $5/2$. Thus the transmitted fraction of light intensity is $(4{\cos }^2\theta +{\sin }^2\theta )/5$: it is $4/5$ at $0^{\circ }$, $1/2$ at $45^{\circ }$ and $1/5$ at $90^{\circ }$. There is no single incident line of linear polarization to insert into Malus's law; projecting the electric field and averaging its square is the appropriate generalization. This measurement cannot distinguish $\delta =\pi /2$ from $\delta =-\pi /2$, because reversing the traversal leaves all these squared averages unchanged.

Consider one selected polarization of waves, represented by complex fields $E_1(t)$ and $E_2(t)$ at two sampling positions. Suppose their second moments are finite and stationary: $J_j=⟨|E_j(t){|}^2⟩$ and the following uncentered wave-field correlation do not depend on the absolute time $t$. DefineThe star denotes the complex conjugate; brackets denote an ensemble average. A long detector average estimates this quantity only if it adequately samples the fluctuations. The superscript $1$ identifies a field correlation: it does not mean the autocorrelation of the measured light intensity. The one-position version is $g^{(1)}(\tau )=⟨E^{\ast }(t)E(t+\tau )⟩/⟨|E{|}^2⟩$, so $g^{(1)}(0)=1$ for every nonzero field covered by this definition; this alone says nothing about its correlation at nonzero delay. ^[ref].

The Cauchy-Schwarz inequality gives $0\le |g_{12}^{(1)}(\tau )|\le 1$. At a specified delay, magnitude one means $E_2(t+\tau )=a{E}_1(t)$ in mean square for a fixed nonzero complex $a$: the samples have a reproducible relation in both wave amplitude and wave phase. Magnitude zero means the averaged field cross term vanishes. It need not mean statistical independence. A deterministic phase difference changes the argument of $g^{(1)}$, not its complex modulus; a distorted but fixed wavefront can therefore remain fully coherent. ^[ref].

To see the observable, recombine the samples in the same detected spatial mode with matched wave polarizations. With a common positive conversion factor $C$, let $I_j=C{J}_j$ be the single-arm light intensities reaching that detector. Adding $E_1(t)+e^{i\varphi }{E}_2(t+\tau )$ givesThis follows simply by expanding the squared complex modulus and averaging, then taking the maximum and minimum as $\varphi$ spans a full cycle. The extra wave phase scan must leave $I_1,I_2$ and the magnitude of the correlation effectively unchanged. Thus interferometric visibility equals $|g^{(1)}|$ only for equal nonzero light intensities and ideal mode matching. For $I_2/I_1=1/4$, perfectly coherent beams have interferometric visibility $0.8$; a measured value $0.4$ would imply $|g^{(1)}|=0.5$ under these assumptions. Orthogonal linear polarizations, background light, unresolved interference fringes or vibrations can reduce measured contrast further. One absent arm gives zero contrast and makes the normalized two-sample correlation undefined.

A fundamental Gaussian beam is an ideal beam of light whose transverse wave amplitude has a Gaussian function profile. Here consider a circularly symmetric beam in vacuum, with propagation along $z$, two transverse coordinates $x,y$, and radius $r=\sqrt{x^2+y^2}$ from the optical axis. With the beam waist at $z=0$, total power $P>0$, and vacuum wavelength $\lambda$, its time-averaged light intensity, meaning irradiance in $\mathrm{W}{\mathrm{m}}^{-2}$, isHere $w_0$ is the minimum beam radius, and $z_R$ is the Rayleigh range. Far from the beam waist, $w(z)\sim \theta |z|$ with beam divergence half-angle $\theta =\lambda /(\pi w_0)$. These are the fundamental circular Gaussian beam laws within the scalar paraxial approximation, appropriate when $w_0\gg \lambda$ and the beam divergence is small. The underlying envelope solves the parabolic wave equation exactly; this does not make it an exact vector solution of Maxwell equations. ^[ref].

The radius convention matters. Writing $U$ for the complex wave amplitude envelope, the same $w(z)$ is the wave amplitude $1/e$ radius and the light intensity $1/e^2$ radius:At $r=w(z)$ these ratios are $e^{-1}\approx 0.367879$ and $e^{-2}\approx 0.135335$, respectively. The light intensity $1/e$ radius is instead $w(z)/\sqrt{2}$. None of these radii is a sharp edge: the ideal Gaussian function has nonzero tails at every finite radius. Kogelnik and Li explicitly use the wave amplitude $1/e$ radius in their fundamental-mode treatment. ^[ref].

For a numerical example, take vacuum wavelength $\lambda =632.8\mathrm{nm}$ and beam waist radius $w_0=0.500\mathrm{mm}$. Substitution gives Rayleigh range $z_R\approx 1.24114754\mathrm{m}$ and far-field beam divergence half-angle $\theta \approx 0.402852992\mathrm{mrad}$. At $z=z_R$, the beam radius is $w=\sqrt{2}{w}_0\approx 0.707107\mathrm{mm}$ and the on-axis irradiance is half its waist value. At the same wavelength, halving the beam waist radius to $0.250\mathrm{mm}$ quarters the Rayleigh range to $0.310286885\mathrm{m}$ and doubles the beam divergence half-angle to $0.805705984\mathrm{mrad}$. At fixed power, that smaller beam waist also quadruples the peak irradiance.

The contours in the left panel are $r=\pm w(z)$ in a plane through the optical axis; they mark an light intensity level, not an outer boundary. The right panel compares the exact model radii with the far-field asymptotes $w=\theta z$. The horizontal and vertical axes use different length units, so the drawn slopes are not physical angles.

The normalization expresses conserved power during free propagation. In transverse polar coordinates, an annulus has area $2\pi rdr$, soA proof outline is the substitution $s=2r^2/w(z{)}^2$ in this integral. Spreading reduces peak irradiance as $w^{-2}$ while increasing transverse area as $w^2$, leaving power unchanged. In particular, the disk of radius $w(z)$ contains $1-e^{-2}\approx 0.8646647168$, or $86.4665\%$, of the total power; the local light intensity at its rim is only $13.5335\%$ of the on-axis light intensity. These are different quantities.

The left panel compares normalized wave amplitude magnitude and irradiance; the right panel integrates the transverse area to show enclosed power. For an ideal absorbing circular aperture of radius $a$ centered on the optical axis, the same fraction gives the power passing the aperture at that plane. An aperture of radius $w$ clips about $13.5335\%$ of the incident power; radius $2w$ clips only $e^{-8}\approx 0.000335463$, or $0.0335463\%$. After clipping, extra diffraction means the transmitted field no longer obeys the unmodified ideal Gaussian beam law. Every finite aperture cuts some ideal tail. Kogelnik and Li distinguish their large-aperture beam analysis from their finite-aperture diffraction analysis. ^[ref].

The useful tradeoff is that a small bright spot cannot remain narrow indefinitely. An optical lens can create a smaller beam waist, but diffraction then makes the beam spread faster. The particularly simple relation $w_0\theta =\lambda /\pi$ ties localization directly to beam divergence. A proof outline for the propagation law starts with $U(x,y,0)=U_0{e}^{-(x^2+y^2)/w_0^2}$. The parabolic wave equation separates into two transverse factors, and a Gaussian integral in each coordinate givesHere $q$ is dimensionless and $k$ is the wavenumber; the carrier convention is $e^{i(kz-\omega t)}$, with angular frequency $\omega$. Taking the squared complex modulus gives the stated $w(z)$ and its $w^{-2}$ light intensity scaling. Each transverse factor contributes $q^{-1/2}$, explaining the $q^{-1}$ prefactor here and the different prefactor in the preserved Gaussian beam with one transverse coordinate.

Historically, Gaussian beam theory developed alongside early laser resonator theory and lens-guided transmission. In their 1966 review, Kogelnik and Li credit Boyd and Gordon, and Goubau and Schwering, with early beam treatments; they summarize earlier work rather than claim the whole construction as a new discovery. Their paper brings beam propagation, focusing and finite-aperture diffraction into a common practical account. ^[ref].