Geometrical Irradiance Changes in a Symmetric Optical System
Total Page:16
File Type:pdf, Size:1020Kb
Geometrical irradiance changes in a symmetric optical system Dmitry Reshidko Jose Sasian Dmitry Reshidko, Jose Sasian, “Geometrical irradiance changes in a symmetric optical system,” Opt. Eng. 56(1), 015104 (2017), doi: 10.1117/1.OE.56.1.015104. Downloaded From: https://www.spiedigitallibrary.org/journals/Optical-Engineering on 11/28/2017 Terms of Use: https://www.spiedigitallibrary.org/terms-of-use Optical Engineering 56(1), 015104 (January 2017) Geometrical irradiance changes in a symmetric optical system Dmitry Reshidko* and Jose Sasian University of Arizona, College of Optical Sciences, 1630 East University Boulevard, Tucson 85721, United States Abstract. The concept of the aberration function is extended to define two functions that describe the light irra- diance distribution at the exit pupil plane and at the image plane of an axially symmetric optical system. Similar to the wavefront aberration function, the irradiance function is expanded as a polynomial, where individual terms represent basic irradiance distribution patterns. Conservation of flux in optical imaging systems is used to derive the specific relation between the irradiance coefficients and wavefront aberration coefficients. It is shown that the coefficients of the irradiance functions can be expressed in terms of wavefront aberration coefficients and first- order system quantities. The theoretical results—these are irradiance coefficient formulas—are in agreement with real ray tracing. © 2017 Society of Photo-Optical Instrumentation Engineers (SPIE) [DOI: 10.1117/1.OE.56.1.015104] Keywords: aberration theory; wavefront aberrations; irradiance. Paper 161674 received Oct. 26, 2016; accepted for publication Dec. 20, 2016; published online Jan. 23, 2017. 1 Introduction the normalized field H~ and aperture ρ~ vectors. The field vec- In the development of the theory of aberrations in optical tor and the aperture vector may be defined in either the object imaging systems, emphasis has been given to the study of or the image spaces. Two vectors uniquely specify a ray image aberrations, which are described as wave, angular, propagating in the lens system. The wavefront aberration transverse, or longitudinal quantities.1 The light irradiance function is expanded into a polynomial of the rotational variation, specifically at the exit pupil plane and at the invariants as dot products of the field and aperture vectors, H~ H~ H~ ρ~ ρ~ ρ~ image plane of an optical system, is a radiometric aspect specifically ð · Þ, ð · Þ, and ð · Þ, and to the fourth H~ ρ~ of a system that is also of interest. order of approximation on and is written as follows: EQ-TARGET;temp:intralink-;e001;326;425 X The geometrical irradiance calculation is useful for deter- ~ ~ ~ j ~ m n WðH; ρ~Þ¼ Wk;l;mðH · HÞ ðH · ρ~Þ ðρ~ · ρ~Þ mining the apodization of the wavefront at the exit pupil. The j;m;l point spread function (PSF) describes the response of an im- ~ ~ ~ aging system to a point object. An accurate diffraction cal- ¼ W000 þ W200ðH · HÞþW111ðH · ρ~Þ ’ culation of a system s PSF requires knowledge not only of W ρ~ ρ~ W ρ~ ρ~ 2 the wavefront phase but also of its amplitude. These, the þ 020ð · Þþ 040ð · Þ ~ ~ 2 phase and amplitude, are usually calculated geometrically þ W131ðH · ρ~Þðρ~ · ρ~ÞþW222ðH · ρ~Þ at the exit pupil of the system. While the wavefront phase ~ ~ is evaluated by calculating the optical path length for þ W220ðH · HÞðρ~ · ρ~Þ rays, the wavefront amplitude is estimated by the square W H~ H~ H~ ρ~ W H~ H~ 2; (1) root of the geometrical irradiance. þ 311ð · Þð · Þþ 400ð · Þ The image plane illumination or relative illumination is where each aberration coefficient Wk;l;m represents the ampli- another characteristic that may have a strong impact on tude of basic wavefront deformation forms and subindices j, m, the performance of a lens system. and n represent integers k ¼ 2j þ m and l ¼ 2n þ m.3 The light irradiance distribution at the exit pupil plane Similarly, a pupil aberration function WðH;~ ρ~Þ can be defined and/or at the image plane of an optical system can be derived to describe the aberration between the pupils. This function is from basic radiometric principles, such as conservation of constructed by interchanging the role of the field and aperture flux.2 In this paper, we provide a study of the relationship vectors, and to fourth order, it is as follows: between irradiance at these two planes and the system’s EQ-TARGET;temp:intralink-;e002;326;217 X wavefront aberration coefficients. Our study is based on geo- ~ ~ ~ j ~ m n WðH; ρ~Þ¼ Wk;l;mðH · HÞ ðH · ρ~Þ ðρ~ · ρ~Þ metrical optics. Edge diffraction effects are not considered in j;m;l this study. An extended object that has a uniform or ~ Lambertian radiance is assumed. The geometrical irradiance ¼ W000 þ W200ðρ~ · ρ~ÞþW111ðH · ρ~Þ that we discuss describes the illumination distribution on the W H~ H~ W H~ H~ 2 image plane or on the exit pupil plane of an optical system. þ 020ð · Þþ 040ð · Þ The concept of the wavefront aberration function is well ~ ~ ~ ~ 2 þ W131ðH · HÞðH · ρ~ÞþW222ðH · ρ~Þ established. The wavefront aberration function WðH;~ ρ~Þ of ~ ~ an axially symmetric system gives the geometrical wavefront þ W220ðH · HÞðρ~ · ρ~Þ deformation (from a sphere) at the exit pupil as a function of ~ 2 þ W311ðH · ρ~Þðρ~ · ρ~ÞþW400ðρ~ · ρ~Þ ; (2) *Address all correspondence to: Dmitry Reshidko, E-mail: dmitry@optics. arizona.edu 0091-3286/2017/$25.00 © 2017 SPIE Optical Engineering 015104-1 January 2017 • Vol. 56(1) Downloaded From: https://www.spiedigitallibrary.org/journals/Optical-Engineering on 11/28/2017 Terms of Use: https://www.spiedigitallibrary.org/terms-of-use Reshidko and Sasian: Geometrical irradiance changes in a symmetric optical system where the pupil aberration coefficients are barred to distin- assumed that the irradiance is greater than zero for any guish them from the image aberration coefficients and point defined by the normalized field H~ and aperture ~ρ subindices j, m,andn represent integers k ¼ 2j þ m vectors. and l ¼ 2n þ m. The question we pose and answer is: what is the relation- In analogy to the wavefront aberration function, we define ship between the wavefront aberration functions WðH;~ ρ~Þ the irradiance function IðH;~ ρ~Þ that gives the irradiance at the and WðH;~ ρ~Þ and the irradiance functions IðH;~ ρ~Þ and image plane of an optical system. To the fourth order of I¯ðH;~ ~ρÞ? We use conservation of flux in an optical system approximation, this irradiance function is expressed as to determine the relationship between wavefront and irradi- follows: ance coefficients. We also discuss in some detail the particu- lar case of relative illumination at the image plane. Our EQ-TARGET;temp:intralink-;e003;63;653 X ~ ~ ~ j ~ m n IðH; ρ~Þ¼ Il;k;mðH · HÞ · ðH · ρ~Þ · ðρ~ · ρ~Þ theoretical results are in agreement with real ray tracing. j;m;n Overall, the paper provides new insights into the irradiance I I H~ H~ I H~ ρ~ I ρ~ ρ~ changes in an optical system and furthers the theory of ¼ 000 þ 200ð · Þþ 111ð · Þþ 020ð · Þ aberrations. 2 ~ ~ 2 þ I040ðρ~ · ρ~Þ þ I131ðH · ρ~Þðρ~ · ρ~ÞþI222ðH · ρ~Þ ~ ~ ~ ~ ~ þ I220ðH · HÞðρ~ · ρ~ÞþI311ðH · HÞðH · ρ~Þ 2 Radiative Transfer in an Optical System ~ ~ 2 In this section, we review radiometry and show how aberra- þ I400ðH · HÞ ; (3) tions relate to conservation of optical flux Φ in an optical system. Figure 2 illustrates the basic elements of an axially where the irradiance coefficients Ik;l;m represent basic illumi- symmetric optical system and the geometry defining the nation distribution patterns at the image plane, as shown in transfer of radiant energy. Rays from a differential area Fig. 1, and subindices j, m, and n represent integers dA in the object plane pass through the optical system k ¼ 2j þ m and l ¼ 2n þ m. We also define another irradi- and converge at the conjugate area dA 0 in the image ance function I¯ðH;~ ρ~Þ that gives the irradiance of the beam at plane. The differential cross sections of the beam dS in the exit pupil plane of an optical system: the entrance pupil plane and dS 0 in the exit pupil plane EQ-TARGET;temp:intralink-;e004;63;469 X I¯ H;~ ρ~ I¯ H~ H~ j H~ ρ~ m ρ~ ρ~ n are optically conjugated to the first order. ð Þ¼ l;k;mð · Þ ð · Þ ð · Þ In a lossless and passive optical system, the element of j;m;l radiant flux dΦ is conserved through all cross sections of ¯ ¯ ¯ ~ ¯ ~ ~ ¼ I000 þ I200ðρ~ · ρ~ÞþI111ðH · ρ~ÞþI020ðH · HÞ the beam. The equation for the conservation of flux is written ¯ ~ ~ 2 ¯ ~ ~ ~ as follows: þ I040ðH · HÞ þ I131ðH · HÞðH · ρ~Þ ¯ ~ 2 ¯ ~ ~ dAdS cos4ðθÞ dA 0dS 0 cos4ðθ 0Þ þ I222ðH · ρ~Þ þ I220ðH · HÞðρ~ · ρ~Þ dΦ L L dΦ 0; (5) EQ-TARGET;temp:intralink-;e005;326;421 ¼ o ¼ o ¼ t2 t 02 ¯ ~ ¯ 2 þ I311ðH · ρ~Þðρ~ · ρ~ÞþI400ðρ~ · ρ~Þ ; (4) where Lo is the source radiance. We assume that the source is ¯ where each irradiance coefficient Il;k;m represents basic Lambertian, and, consequently, Lo is constant. The object apodization distributions at the exit pupil, as shown in space angle θ is between the ray connecting dA and dS Fig. 1, and subindices j, m, and n represent integers and the optical axis of the lens system. Similarly, θ 0 is k ¼ 2j þ m and l ¼ 2n þ m.