Airy disk

In optics, the Airy disk (or Airy disc) and Airy pattern are descriptions of the best-focused spot of light that a perfect lens with a circular aperture can make, limited by the diffraction of light. The Airy disk is of importance in physics, optics, and astronomy.
The diffraction pattern resulting from a uniformly illuminated, circular aperture has a bright central region, known as the Airy disk, which together with the series of concentric rings around is called the Airy pattern. Both are named after George Biddell Airy. The disk and rings phenomenon had been known prior to Airy; John Herschel described the appearance of a bright star seen through a telescope under high magnification for an 1828 article on light for the Encyclopedia Metropolitana:
Airy wrote the first full theoretical treatment explaining the phenomenon (his 1835 "On the Diffraction of an Object-glass with Circular Aperture").
Mathematically, the diffraction pattern is characterized by the wavelength of light illuminating the circular aperture, and the aperture's size. The appearance of the diffraction pattern is additionally characterized by the sensitivity of the eye or other detector used to observe the pattern.
The most important application of this concept is in cameras, microscopes and telescopes. Due to diffraction, the smallest point to which a lens or mirror can focus a beam of light is the size of the Airy disk. Even if one were able to make a perfect lens, there is still a limit to the resolution of an image created by such a lens. An optical system in which the resolution is no longer limited by imperfections in the lenses but only by diffraction is said to be diffraction limited.
Far from the aperture, the angle at which the first minimum occurs, measured from the direction of incoming light, is given by the approximate formula:
sin ⁡ θ ≈ 1.22 λ d {\displaystyle \sin \theta \approx 1.22{\frac {\lambda }{d}}}
or, for small angles, simply
θ ≈ 1.22 λ d , {\displaystyle \theta \approx 1.22{\frac {\lambda }{d}},}
where θ {\displaystyle \theta } is in radians, λ {\displaystyle \lambda } is the wavelength of the light in meters, and d {\displaystyle {d}} is the diameter of the aperture in meters. The full width at half maximum is given by θ F W H M = 1.029 λ d . {\displaystyle \theta _{\mathrm {FWHM} }=1.029{\frac {\lambda }{d}}.}
Airy wrote this relation as
s = 2.76 a , {\displaystyle s={\frac {2.76}{a}},}
where s {\displaystyle {s}} was the angle of first minimum in seconds of arc, a {\displaystyle {a}} was the radius of the aperture in inches, and the wavelength of light was assumed to be 0.000022 inches (560 nm; the mean of visible wavelengths). This is equal to the angular resolution of a circular aperture. The Rayleigh criterion for barely resolving two objects that are point sources of light, such as stars seen through a telescope, is that the center of the Airy disk for the first object occurs at the first minimum of the Airy disk of the second. This means that the angular resolution of a diffraction-limited system is given by the same formulae.
However, while the angle at which the first minimum occurs (which is sometimes described as the radius of the Airy disk) depends only on wavelength and aperture size, the appearance of the diffraction pattern will vary with the intensity (brightness) of the light source. Because any detector (eye, film, digital) used to observe the diffraction pattern can have an intensity threshold for detection, the full diffraction pattern may not be apparent. In astronomy, the outer rings are frequently not apparent even in a highly magnified image of a star. It may be that none of the rings are apparent, in which case the star image appears as a disk (central maximum only) rather than as a full diffraction pattern. Furthermore, fainter stars will appear as smaller disks than brighter stars, because less of their central maximum reaches the threshold of detection. While in theory all stars or other "point sources" of a given wavelength and seen through a given aperture have the same Airy disk radius characterized by the above equation (and the same diffraction pattern size), differing only in intensity, the appearance is that fainter sources appear as smaller disks, and brighter sources appear as larger disks. This was described by Airy in his original work:
Despite this feature of Airy's work, the radius of the Airy disk is often given as being simply the angle of first minimum, even in standard textbooks. In reality, the angle of first minimum is a limiting value for the size of the Airy disk, and not a definite radius.
If two objects imaged by a camera are separated by an angle small enough that their Airy disks on the camera detector start overlapping, the objects cannot be clearly separated any more in the image, and they start blurring together. Two objects are said to be just resolved when the maximum of the first Airy pattern falls on top of the first minimum of the second Airy pattern (the Rayleigh criterion).
Therefore, the smallest angular separation two objects can have before they significantly blur together is given as stated above by sin ⁡ θ = 1.22 λ d . {\displaystyle \sin \theta =1.22\,{\frac {\lambda }{d}}.}
Thus, the ability of the system to resolve detail is limited by the ratio of λ/d. The larger the aperture for a given wavelength, the finer the detail that can be distinguished in the image.
This can also be expressed as x f = 1.22 λ d , {\displaystyle {\frac {x}{f}}=1.22\,{\frac {\lambda }{d}},} where x {\displaystyle x} is the separation of the images of the two objects on the film, and f {\displaystyle f} is the distance from the lens to the film. If we take the distance from the lens to the film to be approximately equal to the focal length of the lens, we find x = 1.22 λ f d , {\displaystyle x=1.22\,{\frac {\lambda \,f}{d}},} but f d {\displaystyle {\frac {f}{d}}} is the f-number of a lens. A typical setting for use on an overcast day would be f/8 (see Sunny 16 rule). For violet, the shortest wavelength visible light, the wavelength λ is about 420 nanometers (see cone cells for sensitivity of S cone cells). This gives a value for x {\displaystyle x} of about 4 μm. In a digital camera, making the pixels of the image sensor smaller than half this value (one pixel for each object, one for each space between) would not significantly increase the captured image resolution. However, it may improve the final image by over-sampling, allowing noise reduction.
The fastest f-number for the human eye is about 2.1, corresponding to a diffraction-limited point spread function with approximately 1 μm diameter. However, at this f-number, spherical aberration limits visual acuity, while a 3 mm pupil diameter (f/5.7) approximates the resolution achieved by the human eye. The maximum density of cones in the human fovea is approximately 170,000 per square millimeter, which implies that the cone spacing in the human eye is about 2.5 μm, approximately the diameter of the point spread function at f/5.
A circular laser beam with uniform intensity across the circle (a flat-top beam) focused by a lens will form an Airy disk pattern at the focus. The size of the Airy disk determines the laser intensity at the focus.
Some weapon aiming sights (e.g. FN FNC) require the user to align a peep sight (rear, nearby sight, i.e. which will be out of focus) with a tip (which should be focused and overlaid on the target) at the end of the barrel. When looking through the peep sight, the user will notice an Airy disk that will help center the sight over the pin.
Light from a uniformly illuminated circular aperture (or from a uniform, flattop beam) will exhibit an Airy diffraction pattern far away from the aperture due to Fraunhofer diffraction (far-field diffraction).
The conditions for being in the far field and exhibiting an Airy pattern are: the incoming light illuminating the aperture is a plane wave (no phase variation across the aperture), the intensity is constant over the area of the aperture, and the distance R {\displaystyle R} from the aperture where the diffracted light is observed (the screen distance) is large compared to the aperture size, and the radius a {\displaystyle a} of the aperture is not too much larger than the wavelength λ {\displaystyle \lambda } of the light. The last two conditions can be formally written as R > a 2 / λ . {\displaystyle R>a^{2}/\lambda .}
In practice, the conditions for uniform illumination can be met by placing the source of the illumination far from the aperture. If the conditions for far field are not met (for example if the aperture is large), the far-field Airy diffraction pattern can also be obtained on a screen much closer to the aperture by using a lens right after the aperture (or the lens itself can form the aperture). The Airy pattern will then be formed at the focus of the lens rather than at infinity.
Hence, the focal spot of a uniform circular laser beam (a flattop beam) focused by a lens will also be an Airy pattern.
In a camera or imaging system an object far away gets imaged onto the film or detector plane by the objective lens, and the far field diffraction pattern is observed at the detector. The resulting image is a convolution of the ideal image with the Airy diffraction pattern due to diffraction from the iris aperture or due to the finite size of the lens. This leads to the finite resolution of a lens system described above.
The intensity of the Airy pattern follows the Fraunhofer diffraction pattern of a circular aperture, given by the squared modulus of the Fourier transform of the circular aperture:
I ( θ ) = I 0 [ 2 J 1 ( k a sin ⁡ θ ) k a sin ⁡ θ ] 2 = I 0 [ 2 J 1 ( x ) x ] 2 {\displaystyle I(\theta )=I_{0}\left[{\frac {2J_{1}(k\,a\sin \theta )}{k\,a\sin \theta }}\right]^{2}=I_{0}\left[{\frac {2J_{1}(x)}{x}}\right]^{2}}
where I 0 {\displaystyle I_{0}} is the maximum intensity of the pattern at the Airy disc center, J 1 {\displaystyle J_{1}} is the Bessel function of the first kind of order one, k = 2 π / λ {\displaystyle k={2\pi }/{\lambda }} is the wavenumber, a {\displaystyle a} is the radius of the aperture, and θ {\displaystyle \theta } is the angle of observation, i.e. the angle between the axis of the circular aperture and the line between aperture center and observation point. x = k a sin ⁡ θ = 2 π a λ q R , {\displaystyle x=ka\sin \theta ={\frac {2\pi a}{\lambda }}{\frac {q}{R}},} where q is the radial distance from the observation point to the optical axis and R is its distance to the aperture. Note that the Airy disk as given by the above expression is only valid for large R, where Fraunhofer diffraction applies; calculation of the shadow in the near-field must rather be handled using Fresnel diffraction.
However the exact Airy pattern does appear at a finite distance if a lens is placed at the aperture. Then the Airy pattern will be perfectly focussed at the distance given by the lens's focal length (assuming collimated light incident on the aperture) given by the above equations.
The zeros of J 1 ( x ) {\displaystyle J_{1}(x)} are at x = k a sin ⁡ θ ≈ 3.8317 , 7.0156 , 10.1735 , 13.3237 , 16.4706 … . {\displaystyle x=ka\sin \theta \approx 3.8317,7.0156,10.1735,13.3237,16.4706\dots .} From this, it follows that the first dark ring in the diffraction pattern occurs where k a sin ⁡ θ = 3.8317 … , {\displaystyle ka\sin {\theta }=3.8317\dots ,} or
sin ⁡ θ ≈ 3.83 k a = 3.83 λ 2 π a = 1.22 λ 2 a = 1.22 λ d . {\displaystyle \sin \theta \approx {\frac {3.83}{ka}}={\frac {3.83\lambda }{2\pi a}}=1.22{\frac {\lambda }{2a}}=1.22{\frac {\lambda }{d}}.}
If a lens is used to focus the Airy pattern at a finite distance, then the radius q 1 {\displaystyle q_{1}} of the first dark ring on the focal plane is solely given by the numerical aperture A (closely related to the f-number) by q 1 = R sin ⁡ θ 1 ≈ 1.22 R λ d = 1.22 λ 2 A {\displaystyle q_{1}=R\sin \theta _{1}\approx 1.22{R}{\frac {\lambda }{d}}=1.22{\frac {\lambda }{2A}}}
where the numerical aperture A is equal to the aperture's radius d/2 divided by R', the distance from the center of the Airy pattern to the edge of the aperture. Viewing the aperture of radius d/2 and lens as a camera (see diagram above) projecting an image onto a focal plane at distance f, the numerical aperture A is related to the commonly-cited f-number N= f/d (ratio of the focal length to the lens diameter) according to
A = r R ′ = r f 2 + r 2 = 1 4 N 2 + 1 ; {\displaystyle A={\frac {r}{R'}}={\frac {r}{\sqrt {f^{2}+r^{2}}}}={\frac {1}{\sqrt {4N^{2}+1}}};}
for N≫1 it is simply approximated as A ≈ 1 / 2 N . {\textstyle A\approx 1/2N.} This shows that the best possible image resolution of a camera is limited by the numerical aperture (and thus f-number) of its lens due to diffraction.
The half maximum of the central Airy disk (where 2 J 1 ( x ) / x = 1 / 2 {\displaystyle 2J_{1}(x)/x=1/{\sqrt {2}}} ) occurs at x = 1.61633995 … ; {\displaystyle x=1.61633995\dots ;} the 1/e point (where 2 J 1 ( x ) / x = 1 / e {\displaystyle 2J_{1}(x)/x=1/{e}} ) occurs at x = 2.58383899 … , {\displaystyle x=2.58383899\dots ,} and the maximum of the first ring occurs at x = 5.13562230 … . {\displaystyle x=5.13562230\dots .}
The intensity I 0 {\displaystyle I_{0}} at the center of the diffraction pattern is related to the total power P 0 {\displaystyle P_{0}} incident on the aperture by
I 0 = E A 2 A 2 2 R 2 = P 0 A λ 2 R 2 {\displaystyle I_{0}={\frac {\mathrm {E} _{A}^{2}A^{2}}{2R^{2}}}={\frac {P_{0}A}{\lambda ^{2}R^{2}}}}
where E {\displaystyle \mathrm {E} } is the source strength per unit area at the aperture, A is the area of the aperture ( A = π a 2 {\displaystyle A=\pi a^{2}} ) and R is the distance from the aperture. At the focal plane of a lens, I 0 = ( P 0 A ) / ( λ 2 f 2 ) . {\displaystyle I_{0}=(P_{0}A)/(\lambda ^{2}f^{2}).} The intensity at the maximum of the first ring is about 1.75% of the intensity at the center of the Airy disk.
The expression for I ( θ ) {\displaystyle I(\theta )} above can be integrated to give the total power contained in the diffraction pattern within a circle of given size:
P ( θ ) = P 0 [ 1 − J 0 2 ( k a sin ⁡ θ ) − J 1 2 ( k a sin ⁡ θ ) ] {\displaystyle P(\theta )=P_{0}[1-J_{0}^{2}(ka\sin \theta )-J_{1}^{2}(ka\sin \theta )]}
where J 0 {\displaystyle J_{0}} and J 1 {\displaystyle J_{1}} are Bessel functions. Hence the fractions of the total power contained within the first, second, and third dark rings (where J 1 ( k a sin ⁡ θ ) = 0 {\displaystyle J_{1}(ka\sin \theta )=0} ) are 83.8%, 91.0%, and 93.8% respectively.
Classical treatments of the Airy disk and diffraction pattern assume that the incident light is a plane wave that consists of coherent (in phase) photons of the same wavelength that interfere with each other. The famous double slit experiment showed that diffraction patterns could arise even when the coherent photons were so spread out in time that they could not interfere with each other. This led to the quantum mechanical picture that each photon effectively takes all possible paths from a source to a detector. Richard Feynman explained that each path has a complex amplitude that can be thought of as a unit vector that is perpendicular to the path and makes one complete rotation for each wavelength of advance. The detection probability is the square of the modulus of the sum of the complex amplitudes at the detector. Diffraction patterns arise because the paths sum differently at different detector positions. According to these principles the Airy disk and diffraction pattern can be computed numerically by using Feynman photon path integrals to determine the detection probability at different points in the focal plane of a parabolic mirror.
The Airy pattern falls rather slowly to zero with increasing distance from the center, with the outer rings containing a significant portion of the integrated intensity of the pattern. As a result, the root mean square (RMS) spotsize is undefined (i.e. infinite). An alternative measure of the spot size is to ignore the relatively small outer rings of the Airy pattern and to approximate the central lobe with a Gaussian profile, such that
I ( q ) ≈ I 0 ′ exp ⁡ ( − 2 q 2 ω 0 2 ) , {\displaystyle I(q)\approx I'_{0}\exp \left({\frac {-2q^{2}}{\omega _{0}^{2}}}\right)\ ,}
where I 0 ′ {\displaystyle I'_{0}} is the irradiance at the center of the pattern, q {\displaystyle q} represents the radial distance from the center of the pattern, and ω 0 {\textstyle \omega _{0}} is the Gaussian RMS width (in one dimension). If we equate the peak amplitude of the Airy pattern and Gaussian profile, that is, I 0 ′ = I 0 , {\displaystyle I'_{0}=I_{0},} and find the value of ω 0 {\textstyle \omega _{0}} giving the optimal approximation to the pattern, we obtain
where N is the f-number. If, on the other hand, we wish to enforce that the Gaussian profile has the same volume as does the Airy pattern, then this becomes
In optical aberration theory, it is common to describe an imaging system as diffraction-limited if the Airy disk radius is larger than the RMS spotsize determined from geometric ray tracing (see Optical lens design). The Gaussian profile approximation provides an alternative means of comparison: using the approximation above shows that the Gaussian waist ω 0 {\textstyle \omega _{0}} of the Gaussian approximation to the Airy disk is about two-third the Airy disk radius, i.e. 0.84 λ N {\displaystyle 0.84\lambda N} as opposed to 1.22 λ N . {\displaystyle 1.22\lambda N.}
Similar equations can also be derived for the obscured Airy diffraction pattern which is the diffraction pattern from an annular aperture or beam, i.e. a uniform circular aperture (beam) obscured by a circular block at the center. This situation is relevant to many common reflector telescope designs that incorporate a secondary mirror, including Newtonian telescopes and Schmidt–Cassegrain telescopes.
I ( R ) = I 0 ( 1 − ϵ 2 ) 2 ( 2 J 1 ( x ) x − 2 ϵ J 1 ( ϵ x ) x ) 2 {\displaystyle I(R)={\frac {I_{0}}{(1-\epsilon ^{2})^{2}}}\left({\frac {2J_{1}(x)}{x}}-{\frac {2\epsilon J_{1}(\epsilon x)}{x}}\right)^{2}}
where ϵ {\displaystyle \epsilon } is the annular aperture obscuration ratio, or the ratio of the diameter of the obscuring disk and the diameter of the aperture (beam). ( 0 ≤ ϵ < 1 ) , {\displaystyle \left(0\leq \epsilon <1\right),} and x is defined as above: x = k a sin ⁡ ( θ ) ≈ π R λ N {\displaystyle x=ka\sin(\theta )\approx {\frac {\pi R}{\lambda N}}} where R {\displaystyle R} is the radial distance in the focal plane from the optical axis, λ {\displaystyle \lambda } is the wavelength and N {\displaystyle N} is the f-number of the system. The fractional encircled energy (the fraction of the total energy contained within a circle of radius R {\displaystyle R} centered at the optical axis in the focal plane) is then given by:
E ( R ) = 1 ( 1 − ϵ 2 ) ( 1 − J 0 2 ( x ) − J 1 2 ( x ) + ϵ 2 [ 1 − J 0 2 ( ϵ x ) − J 1 2 ( ϵ x ) ] − 4 ϵ ∫ 0 x J 1 ( t ) J 1 ( ϵ t ) t d t ) {\displaystyle E(R)={\frac {1}{(1-\epsilon ^{2})}}\left(1-J_{0}^{2}(x)-J_{1}^{2}(x)+\epsilon ^{2}\left[1-J_{0}^{2}(\epsilon x)-J_{1}^{2}(\epsilon x)\right]-4\epsilon \int _{0}^{x}{\frac {J_{1}(t)J_{1}(\epsilon t)}{t}}\,dt\right)}
For ϵ → 0 {\displaystyle \epsilon \rightarrow 0} the formulas reduce to the unobscured versions above.
The practical effect of having a central obstruction in a telescope is that the central disc becomes slightly smaller, and the first bright ring becomes brighter at the expense of the central disc. This becomes more problematic with short focal length telescopes which require larger secondary mirrors.
A circular laser beam with uniform intensity profile, focused by a lens, will form an Airy pattern at the focal plane of the lens. The intensity at the center of the focus will be I 0 , A i r y = ( P 0 A ) / ( λ 2 f 2 ) {\displaystyle I_{0,Airy}=(P_{0}A)/(\lambda ^{2}f^{2})} where P 0 {\displaystyle P_{0}} is the total power of the beam, A = π D 2 / 4 {\displaystyle A=\pi D^{2}/4} is the area of the beam ( D {\displaystyle D} is the beam diameter), λ {\displaystyle \lambda } is the wavelength, and f {\displaystyle f} is the focal length of the lens.
A Gaussian beam transmitted through a hard aperture will be clipped. Energy is lost and edge diffraction occurs, effectively increasing the divergence. Because of these effects there is a Gaussian beam diameter which maximizes the intensity in the far field. This occurs when the 1 / e 2 {\displaystyle 1/e^{2}} diameter of the Gaussian is 89% of the aperture diameter, and the on axis intensity in the far field will be 81% of that produced by a uniform intensity profile.
The Fourier integral of the circular cross section of radius a {\displaystyle a} is
This is the special case of the Fourier integral of the elliptical cross section with half axes a {\displaystyle a} and b {\displaystyle b} :
where