• Journal of the European Optical Society-Rapid Publications
  • Vol. 18, Issue 1, 2022004 (2022)
Psang Lin*
Author Affiliations
  • Department of Mechanical Engineering, Taiwan Cheng Kung University, Tainan 70101, China
  • show less
    DOI: 10.1051/jeos/2022004 Cite this Article
    Psang Lin. It is a sufficient condition only, not a necessary and sufficient condition, for decomposing wavefront aberrations[J]. Journal of the European Optical Society-Rapid Publications, 2022, 18(1): 2022004 Copy Citation Text show less

    Abstract

    The classic equation for decomposing the wavefront aberrations of axis-symmetrical optical systems has the form,W(h0,ρ,ϕ)=j=0p=0m=0C2j+m2p+mm(h0)2j+m(ρ)2p+m(cosϕ)mwhere j, p and m are non-negative integers, ρ and ϕ are the polar coordinates of the pupil, and h0 is the object height. However, one non-zero component of the aberrations (i.e., C133h0ρ3cos3ϕ) is missing from this equation when the image plane is not the Gaussian image plane. This implies that the equation is a sufficient condition only, rather than a necessary and sufficient condition, since it cannot guarantee that all of the components of the aberrations can be found. Accordingly, this paper presents a new method for determining all the components of aberrations of any order. The results show that three and six components of the secondary and tertiary aberrations, respectively, are missing in the existing literature.The classic equation for decomposing the wavefront aberrations of axis-symmetrical optical systems has the form,W(h0,ρ,ϕ)=j=0p=0m=0C2j+m2p+mm(h0)2j+m(ρ)2p+m(cosϕ)mwhere j, p and m are non-negative integers, ρ and ϕ are the polar coordinates of the pupil, and h0 is the object height. However, one non-zero component of the aberrations (i.e., C133h0ρ3cos3ϕ) is missing from this equation when the image plane is not the Gaussian image plane. This implies that the equation is a sufficient condition only, rather than a necessary and sufficient condition, since it cannot guarantee that all of the components of the aberrations can be found. Accordingly, this paper presents a new method for determining all the components of aberrations of any order. The results show that three and six components of the secondary and tertiary aberrations, respectively, are missing in the existing literature.

    1 Introduction

    The wavefront and ray aberrations of axis-symmetrical systems have attracted significant attention in the literature [112]. The usual equation for decomposing the monochromatic wavefront aberration W(h0, ρ, ϕ) into different orders and components is given as (e.g., Eq. (3.31b) of [5]),W(h0,ρ,ϕ)=j=0p=0m=0C2j+m2p+mm(h0)2j+m(ρ)2p+m(cosϕ)mwhere j, p and m are non-negative integers, ρ and ϕ are the polar coordinates of the pupil, and h0 is the object height. The sum of the powers of h0 and ρ gives the order of the related component. That is,2q=2j+p+m

    For example, if the piston term is included, the primary (i.e., fourth-order W4th) aberrations are obtained from equation (1) with 2q = 4 as,W4th=C400h04+C040ρ4+C131h0ρ3cosϕ+C222h02ρ2(cosϕ)2+C220h02ρ2+C311h03ρcosϕwhere the six components of the equation represent the piston term and the spherical, coma, astigmatism, field curvature, and distortion aberrations, respectively. The composition ability of the primary aberrations from equation (1) is echoed by Buchdahl [1], who computed the Buchdahl aberration coefficients to determine the wavefront and ray aberrations of axis-symmetrical systems. However, the question arises as to whether equation (1) provides all the components of the various order wavefront aberrations in an axis-symmetrical optical system.

    2 Decomposition of wavefront aberrations

    Equation (1) is based on the fact that the aberration function W(h0, ρ, ϕ) must satisfy the following three equations related to the fundamental axis-symmetrical nature of axis-symmetrical systems (e.g., p. 154 of [5]):W0,ρ,ϕ=W0,-ρ,ϕWh0,ρ,ϕ=Wh0,ρ,-ϕWh0,ρ,ϕ=W-h0,ρ,π+ϕ=W-h0,ρ,π-ϕ

    Table Infomation Is Not EnableTable Infomation Is Not Enable

    Equation (4) indicates that the aberration function of an on-axis object must be radially symmetric, and hence implies that the components of W(h0, ρ, ϕ) that do not depend on h0 should vary as ρ2 (or its integer power). Equation (5) states that W(h0, ρ, ϕ) must be a function of cosϕ. Finally, equation (6) shows that W(h0, ρ, ϕ) must equal W(−h0, ρ, π + ϕ) for an object with height h0 above the optical axis and W(−h0, ρ, π − ϕ) for an object with height h0 below the optical axis. Hence, those terms that depend on ϕ should be a function of h0ρcosϕ. Combining this with ϕ-independent terms, it follows that W(h0, ρ, ϕ) must consist of terms containing h02, ρ2 and h0ρcosϕ factors, to have a sufficient condition given by equation (1). Note that a sufficient condition is taken here to mean that any term generated by equation (1) is a component of an aberration.

    The present group recently proposed a method for determining the aberrations of axis-symmetrical optical systems [12]. It is shown in Figure 1 that when the image plane is not the Gaussian image plane, a non-zero component (i.e., C133h0ρ3(cosϕ)3) is missing from equation (1), even though it satisfies equations (4)(6). This implies that equation (1) alone does not guarantee that all of the components of the aberrations can be found. In other words, equation (1) is not a necessary and sufficient condition for determining all the components of the aberrations in an axis-symmetrical optical system.

    The variation of W133 = C133h0ρ3/λ (where h0 = 17 mm, ρ = 21 mm, and λ = 550 µm) versus the separation of image plane Vimage for the optical system of [12]. The Gaussian image plane of this system is located at Vimage = 92.088474 mm when the object is placed at P0z = −200 mm. This figure shows that C133h0ρ3(cosϕ)3 has non-zero value when the image plane is not the Gaussian image plane. It also shows that C133 = 0 when Vimage = 92.088474 mm, indicating C133h0ρ3(cosϕ)3 is the defocus component of primary aberrations.

    Figure 1.The variation of W133 = C133h0ρ3/λ (where h0 = 17 mm, ρ = 21 mm, and λ = 550 µm) versus the separation of image plane Vimage for the optical system of [12]. The Gaussian image plane of this system is located at Vimage = 92.088474 mm when the object is placed at P0z = −200 mm. This figure shows that C133h0ρ3(cosϕ)3 has non-zero value when the image plane is not the Gaussian image plane. It also shows that C133 = 0 when Vimage = 92.088474 mm, indicating C133h0ρ3(cosϕ)3 is the defocus component of primary aberrations.

    Thus, a second question arises as to how all of these components may be found. To address this question, it is first necessary to realize that the power of cosϕ should be non-negative. That is,m0

    Without any loss of generality, the power of h0 can be confined to a non-negative integer value in order to have 2j + m ≥ 0. That is,j-m/2

    Mathematically, the power of ρ should be greater than or equal to the power of cosϕ, i.e., 2p + m ≥ m, which yields,p0

    One then has the following inequality from the sum of q = j + p + m (see Eq. (2)) and equation (8):pq-m/2

    The intersection of equations (9) and (10) defines the possible range of p. That is,0pq-m/2

    Equation (11) shows that the integer index p starts at p = 0 and ends at,pmax=q-m/2where 〈q − m/2〉 is the maximum non-negative integer value of p for a given m and q. In other words, index p belongs to the following set:p0,1,2,...,pmax

    Furthermore, from equation (11), the possible upper limit of m which yields 0 ≤ q − m/2 is,m2q

    The intersection of equations (7) and (14) then shows the possible domain of integer m for a given q. That is,0m2q

    Given the preceding derivations, it is possible to obtain all components of any order (say, the (2q)th order) wavefront aberration for an axis-symmetrical optical system by the following equation when j = q − (p + m) is used:W2qth-order(h0,ρ,ϕ)=p=0p=pmaxm=0m=2qC2q-2p-m2p+mm(h0)2q-2p-m(ρ)2p+m(cosϕ)m

    Note that, as shown in equation (4), components of the aberrations that do not depend on h0 should vary as ρ2, or as its integer power only. In other words, if 2q − 2p – m = 0 and m0, that component generated from equation (16) does not exist.

    Consider {L-End} Table 1 below, which shows all the components of the secondary aberration (q = 3) of an axis-symmetrical optical system for illustration purposes. The entries of the first and second columns are the values of m and p obtained from equations (15) to (13), respectively. Meanwhile, the entries of the third column denote the sequence (2j + m, 2p + m, m) = (2q − 2p − m, 2p + m, m) for each value of m. The fourth column shows the aberration component for each value of m (if it exists). Comparing {L-End} Table 1 with the existing literature, it is found that three components in {L-End} Table 1 (i.e., C244h02ρ4(cosϕ)4, C153h0ρ5(cosϕ)3 and C155h0ρ5(cosϕ)5) are not included among the secondary aberrations given in the literature despite satisfying equations (4)(6). In order to validate {L-End} Table 1, the methodology proposed in [12] was extended to determine the values of all the secondary aberrations listed in the right-hand column of the table [3]. The results confirmed that all of the secondary aberrations possessed non-zero values.

    The method in this study was further applied to determine all the components of the tertiary aberrations (q = 4) (see Table 2). Comparing the results in Table 2 with those in Table 3–3 of [5], it is found that six components (i.e., C173h0ρ7(cosϕ)3, C264h02ρ6(cosϕ)4, C355h03ρ5(cosϕ)5, C175h0ρ7(cosϕ)5, C266h02ρ6(cosϕ)6 and C177h0ρ7(cosϕ)7) are missing from equation (1).

    3 Conclusions

    The wavefront aberrations W(h0, ρ, ϕ) of axis-symmetrical systems are generally decomposed into their various components using equation (1). However, the numerical results presented in [12] show that this equation cannot guarantee that all of the components of the primary aberrations can be found. In other words, the equation is a sufficient condition only, not a necessary and sufficient condition.

    Accordingly, this study has presented a method for determining the possible domains of the non-negative integer indices, m and p, in equation (1) such that all of the components of the aberrations can be found. It has been shown that the index j computed from equation (2) may be negative. Furthermore, three and six new components of the secondary and tertiary aberrations of an axis-symmetrical system have been found, where these components all satisfy the equations describing the fundamental axis-symmetrical nature of axis-symmetrical systems. Overall, the method proposed in this study provides a systematic and robust approach for ensuring that all of the components of any order wavefront aberration in an axis-symmetrical system can be found.

    References

    [1] H.A. Buchdahl. Optical aberration coefficients(1968).

    [2] W.T. Welford. A new total aberration formula. J. Mod. Opt., 19, 719-727(1972).

    [3] P.D. Lin. Determination of secondary wavefront aberrations in axis-symmetrical optical systems. Appl. Comput. Math., 11, 60-68(2022).

    [4] J. Sasián. Introduction to aberrations in optical imaging systems(2013).

    [5] V.N. Mahajan. Optical imaging and aberrations(1998).

    [6] W.T. Welford. Aberrations of optical systems(1986).

    [7] M. Oleszko, R. Hambach, H. Gros. Decomposition of the total wave aberration in generalized optical systems. J. Opt. Soc. Am., 34, 1856-1864(2017).

    [8] R.B. Johnson. Balancing the astigmatic fields when all other aberrations are absent. Appl. Optics, 32, 3494-3496(1993).

    [9] R. Kingslake, R.B. Johnson. Lens design fundamentals(2010).

    [10] P.D. Lin, R.B. Johnson. Seidel aberration coefficients: An alternative computation method. Opt. Express, 27, 19712-19725(2019).

    [11] P.D. Lin. Seidel primary ray aberration coefficients for objects placed at finite and infinite distances. Opt. Express, 28, 19740-19754(2020).

    [12] P.D. Lin. Taylor series expansion method for determination of primary wavefront aberrations of axis-symmetrical optical systems. OPTIK, Int. J. Light Elect. Opt., 248, Article 168134(2021).

    Psang Lin. It is a sufficient condition only, not a necessary and sufficient condition, for decomposing wavefront aberrations[J]. Journal of the European Optical Society-Rapid Publications, 2022, 18(1): 2022004
    Download Citation