Fraunhofer Diffraction Effects on Total Power for a Planckian Source

An algorithm for computing diffraction effects on total power in the case of Fraunhofer diffraction by a circular lens or aperture is derived. The result for Fraunhofer diffraction of monochromatic radiation is well known, and this work reports the result for radiation from a Planckian source. The result obtained is valid at all temperatures.


Introduction
Fraunhofer diffraction by a circular lens or aperture is a ubiquitous phenomenon in optics in general and radiometry in particular. Figure 1 illustrates two practical situations in which Fraunhofer diffraction occurs. In the first example, diffraction limits the ability of a lens or other focusing optic to focus light. According to geometrical optics, it is possible to focus rays incident on a lens to converge at a focal point. In practice, even with the focal point at the center of a finite-sized circular detector, some of the light incident on the lens fails to be collected, giving rise to a "diffraction loss." The second example can occur when using a source, such as a blackbody cavity source viewed through a small "pinhole" aperture, to calibrate optical systems intended for observing celestial objects or for use in remote-sensing applications. In this case, the angle subtended by the cavity opening at the aperture is larger than the angle subtended by the detector pupil on the dark side of the aperture [1]. Because the former angle is larger than the latter one, geometrical optics suggests that the total power detected depends only on the blackbody temperature and a geometrical factor related to the pinhole aperture, detector optics, and relative separa-tion, because the detector pupil is overfilled. However, diffraction leads to losses in the total power reaching the detector.
All of the above diffraction losses have been a subject of considerable interest, and they have been considered by Blevin [2], Boivin [3], Steel, De, and Bell [4], and Shirley [5]. The formula for the relative diffraction loss in spectral power is already well known, and it is noted below. However, a formula for the diffraction loss in the total power for the case of a Planckian source such as a star or blackbody appears to have been given only in the high-temperature limit. This article reports a formula for the diffraction loss in the power for such a source at all temperatures. The formula is useful for predicting the loss in cases such as the examples discussed, and it can be used as an independent check for numerically calculated diffraction losses.
For the radiation present at a given wavelength , the diffraction loss can be described in terms of a unit-less parameter, v = 2⌿R / . R is the radius of the lens or aperture. The angle ⌿ is defined in either of two ways: either 2⌿ is the full angle subtended by the detection pupil at the focusing optic, or 2⌿ is the full angle Two optical systems are shown, for which the (spectral) power detected is subject to Fraunhofer diffraction. In the first system (a), a lens concentrates radiation on the detector. In the second system (b), radiation from a blackbody cavity may pass through the blackbody cavity opening, pass through the pinhole aperture, and reach the detector pupil. subtended by the blackbody cavity opening at the pinhole aperture. For a source in thermal equilibrium at temperature T , the effects of temperature on the distribution of spectral power enter relevant equations through the ratio, c 2 /(T ) ≡ Av . Here, c 2 = hc /k = 1.438 7752(25) ϫ 10 Ϫ2 m K is the second radiation constant, where the number in parenthesis is the one-standard-deviation uncertainty in the last two digits. This implicitly defines another unit-less parameter, A = c 2 /(2⌿RT ). The relative effects of diffraction on the total power only depend on this one parameter, A .
As a final note, if the detection pupil subtends a considerable angle, the diffraction losses can be approximately given as the weighted average of losses arising, in the case of an infinitesimal detector, for values of ⌿ varying between ⌿ minus the half-angle subtended by the detector and ⌿ plus the same half-angle. Analogous weighting can be used, in the context of the first example, to account for the finite angular diameter of a source (smaller than 2⌿ ). Such averaging is discussed elsewhere [5] but is not discussed further in this work, which henceforth assumes a single effective value of ⌿ .

Derivation of Formula
In the Fraunhofer case, the ratio of the true spectral power to the spectral power predicted by geometrical optics is Here, J m (v ) is a cylindrical Bessel function. It is necessary to incorporate this wavelength-dependent factor when evaluating the ratio of the total power detected to the total power predicted by geometrical optics. For a source whose spectral output obeys the Planck distribution, this unit-less ratio is given by The value of the denominator is familiar, being given by Here, (n ) is the Riemann zeta function. When evaluating the numerator, two techniques have been found to be helpful: one for the "low-temperature" case (defined as A > 4) and one for the "high-temperature" case (defined as A < 4). By such definitions, "high" and "low" temperature cases arise depending on the average diffraction loss for a source at temperature T . This depends on T and the geometry through A . In either case, one first makes use of the relation, from which follows,

High-Temperature Case
To evaluate the numerator in the high-temperature case, the relation in Eq. (2) is rewritten to render with Use of Eq. (5) and the relation, one has Making the abbreviation, w = 4/(A 2 + 4), one has where E (w ) and K (w ) are respectively complete elliptic integrals of the first and second kind, defined according to the convention, 7 (In a different convention, w is replaced by w 2 in the integrand but nowhere else.) From the properties of elliptic integrals [8], one may obtain and The latter result is easily obtained, because Q 1 (A ) has already been written as the third derivative of an expression with respect to A . The closed-form result in Eq. (14) may be directly substituted in Eq. (8), after which one may sum over n . Because the summand can be integrated according to Eq. (15), however, the Euler-Maclaurin formula may be used to obtain a more easily evaluated, asymptotic expression for Q (A ): where B i is again a Bernoulli number. Here a new function, has been introduced as the summand, so that the largest part of Q 1 (A ) is extracted from the sum prior to summation. The parameters m and N are both positive integers. R m (N , A ), the remainder (i.e., error) term, usually decreases initially with increasing m , but then rapidly diverges. It is therefore good to choose a suitable combination of N and m , for a given value of A , to have an efficient yet accurate result. Increasing N should usually improve results, and one obtains an error in F (A ) estimated to be as small as 10 Ϫ14 for A < 4 if one uses N = 60 and m = 5.
Explicitly, Q (A ) is given by One should include only the first m Ϫ 1 terms in the parenthetic expression preceding the second sum. Two abbreviations, ␣ = NA and = 4/(4 + ␣ 2 ), have been introduced. In the second sum, the first several constants are given by , Likewise, the first several e -and k -functions are given by

Evaluation of Formula
In the high-temperature limit, this gives a result of the form and evaluating the successive terms in Eq. (21) may be both difficult and unfruitful, because the series appears to be both very complicated and slow to converge. Note that the leading terms in Eq. (21) are consistent with other works cited, such as the work of Blevin. In the low-temperature limit, one has F (A ) = 10 2 21A 2 Ϫ 4 4A 4 + иии.
(22) Figure 2 shows F (A ) for a range of values of A bridging the low-and high-temperature regions. The solid curve indicates "exact" results obtained using Eq. (7) for A > 4 and Eq. (18) for A < 4. The dashed curves indicate approximate results for low-and high-temperature limits, obtained using the respective formulas in Eqs. (21-22). Likewise, Table 1 shows sample values of F (A ) over a similar range. For 0 < A < 0.2 and A > 8, inclusion of terms shown in Eqs. (21-22) yields an error in F (A ) smaller than 0.001. While this immediate discussion helps provide a sense of the behavior of F (A ), accurate values should be found using Eq. (7) and Eq. (18).