Issue |
J. Eur. Opt. Society-Rapid Publ.
Volume 18, Number 1, 2022
|
|
---|---|---|
Article Number | 3 | |
Number of page(s) | 4 | |
DOI | https://doi.org/10.1051/jeos/2022004 | |
Published online | 06 July 2022 |
Short Communication
It is a sufficient condition only, not a necessary and sufficient condition, for decomposing wavefront aberrations
Department of Mechanical Engineering, National Cheng Kung University, Tainan 70101, Taiwan
* Corresponding author: pdlin@mail.ncku.edu.tw
Received:
14
March
2022
Accepted:
22
June
2022
The classic equation for decomposing the wavefront aberrations of axis-symmetrical optical systems has the form,where 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.
Key words: Aberrations / Imaging systems / Wavefront
© The Author(s), published by EDP Sciences, 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The wavefront and ray aberrations of axis-symmetrical systems have attracted significant attention in the literature [1–12]. 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]),(1)where 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,
(2)
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,(3)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]):(4)
(5)
(6)
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 , ρ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.
![]() |
Fig. 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,(7)
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,(8)
Mathematically, the power of ρ should be greater than or equal to the power of cosϕ, i.e., 2p + m ≥ m, which yields,(9)
One then has the following inequality from the sum of q = j + p + m (see Eq. (2)) and equation (8):(10)
The intersection of equations (9) and (10) defines the possible range of p. That is,(11)
Equation (11) shows that the integer index p starts at p = 0 and ends at,(12)where 〈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:
(13)
Furthermore, from equation (11), the possible upper limit of m which yields 0 ≤ q − m/2 is,(14)
The intersection of equations (7) and (14) then shows the possible domain of integer m for a given q. That is,(15)
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:(16)
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 , that component generated from equation (16) does not exist.
Consider 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 Table 1 with the existing literature, it is found that three components in Table 1 (i.e., C244
ρ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 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.
Components of secondary aberrations with q = 3 in axis-symmetrical system.
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, C264
ρ6(cosϕ)4, C355
ρ5(cosϕ)5, C175h0ρ7(cosϕ)5, C266
ρ6(cosϕ)6 and C177h0ρ7(cosϕ)7) are missing from equation (1).
Components of tertiary aberrations with q = 4 in axis-symmetrical system.
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.
Conflict of interest
The author declares no conflicts of interest.
Funding
Ministry of Science and Technology, Taiwan (MOST) (109-2221-E-006-045).
References
- Buchdahl H.A. (1968) Optical aberration coefficients, Dover. [Google Scholar]
- Welford W.T. (1972) A new total aberration formula, J. Mod. Opt. 19, 719–727. [NASA ADS] [Google Scholar]
- Lin P.D. (2022) Determination of secondary wavefront aberrations in axis-symmetrical optical systems, Appl. Comput. Math. 11, 3, 60–68. [Google Scholar]
- Sasián J. (2013) Introduction to aberrations in optical imaging systems, Cambridge University. [Google Scholar]
- Mahajan V.N. (1998) Optical imaging and aberrations. Part I ray geometrical optics, SPIE Press. [CrossRef] [Google Scholar]
- Welford W.T. (1986) Aberrations of optical systems, Adam Hilger. [Google Scholar]
- Oleszko M., Hambach R., Gros H. (2017) Decomposition of the total wave aberration in generalized optical systems, J. Opt. Soc. Am. 34, 10, 1856–1864. [NASA ADS] [CrossRef] [Google Scholar]
- Johnson R.B. (1993) Balancing the astigmatic fields when all other aberrations are absent, Appl. Optics 32, 19, 3494–3496. [NASA ADS] [CrossRef] [Google Scholar]
- Kingslake R., Johnson R.B. (2010) Lens design fundamentals, 2nd edn., Academic. [Google Scholar]
- Lin P.D., Johnson R.B. (2019) Seidel aberration coefficients: An alternative computation method, Opt. Express 27, 14, 19712–19725. [NASA ADS] [CrossRef] [Google Scholar]
- Lin P.D. (2020) Seidel primary ray aberration coefficients for objects placed at finite and infinite distances, Opt. Express 28, 9, 19740–19754. [NASA ADS] [CrossRef] [Google Scholar]
- Lin P.D. (2021) Taylor series expansion method for determination of primary wavefront aberrations of axis-symmetrical optical systems, OPTIK, Int. J. Light Elect. Opt. 248, Article 168134. [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
![]() |
Fig. 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. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.