Convergence and differentiation of Zernike expansion: application for an analysis of odd-order surfaces

Takao Tanabe; Masato Shibuya; Kazuhisa Maehara

doi:10.1117/1.OE.55.3.035101

1 March 2016 Convergence and differentiation of Zernike expansion: application for an analysis of odd-order surfaces

Takao Tanabe, Masato Shibuya, Kazuhisa Maehara

Author Affiliations +

Optical Engineering, Vol. 55, Issue 3, 035101 (March 2016). https://doi.org/10.1117/1.OE.55.3.035101

Abstract

Odd-order surfaces have begun to be used in optics. In order to investigate the aberration characteristics of such surfaces, Zernike expansion is widely used since it directly and explicitly corresponds to wavefront aberrations. Since the Zernike expansion of an odd-order surface contains an infinite number of terms, the convergence of the expanded sum and the possibility of termwise derivatives are not explicitly guaranteed mathematically. We give a complete proof for these problems. For an application of this result, we analyze the aberration characteristics of odd-order surfaces and present their effectiveness in optical design.

1. Introduction

New types of aspheric surfaces and their mathematical expressions have been proposed in optical design.¹^–⁵ In particular, the effectiveness of odd-order surfaces, which are expressed with terms including an odd-integer power of the radius, has begun to be recognized via the designs of viewfinders,⁶ microscopes,⁷ camera lenses,⁸ projection display optics,⁹ EUV optical systems,¹⁰^–¹² and so on. Since aberration characteristics of odd-order surfaces differ largely from those of conventional aspheric surfaces, it is meaningful and important to analyze the features of odd-order surfaces mathematically.

Shibuya et al.¹¹ pointed out that not only surface shape but also surface slope should be fully approximated to express an optical surface because the refraction of light is determined by the tangent of the surface. From this viewpoint, by considering Taylor expansion, they proved that no odd-order surfaces can be fully represented by even-order terms and thus deduced that the odd-order surfaces have different aberration characteristics from conventional even-order surfaces. They concluded that this peculiarity of odd-order surfaces leads to effectiveness in optical design.¹¹^,¹² However, they ignored the possibility that a finite number of Zernike polynomials fully approximates odd-order surface shapes and slopes.

To analyze aberration characteristics of aspheric surfaces mathematically, Zernike expansion is the most common tool because each Zernike polynomial corresponds to a specific geometrical aberration. Zernike polynomials were originally introduced by Zernike¹³^,¹⁴ as the eigenfunctions of a rotationally symmetrical and self-adjoint differential equation on the unit disk. The orthogonality and the $L^{2}$ -completeness of Zernike polynomials are direct consequences of this definition.¹⁵ Another definition was given by Bhatia and Wolf¹³^,¹⁶ from the viewpoint of symmetry. In this literature, convergence of Zernike expansion is not explicitly discussed, but is regarded as an implicit condition. A mathematician, Szegö,¹⁷ presented an explicit and strong solution for the convergence problem. That is, uniform convergence is fulfilled for a wide class of orthogonal expansion.

Braat and Janssen¹⁸ gave an example of explicit Zernike expansion coefficients by factorial functions or the gamma function. For derivatives of Zernike polynomials, many researchers such as Nijboer¹⁹ obtained the recurring formulae of derivatives to convert wavefront aberrations into geometrical aberrations. This derivation was based on the three term equalities of Gauss’s hypergeometric function. Janssen²⁰ discussed the derivative formulae in connection with a Laplace operator. He presented the effectiveness of his method in the calculation and measurement of geometrical aberrations. He also mentioned that this method is applied to solve the Neumann problem of differential equations. By studying the algebraic structure of differential operators on the complex unit disk, Wünsche²¹ generalized these discussions and derived a generalized form of orthogonal polynomials in the analogy of quantum theory.

As mentioned above, expansion by orthogonal polynomials and its derivatives have been deeply investigated. However, as far as we know, a concrete formulation for odd-order surfaces cannot be found. In particular, both the convergence of Zernike expansion and the convergence of termwise differentiation of Zernike expansion are not discussed in any literature at all.

The goal of this study is to prove that both shapes and slopes of odd-order surfaces are fairly approximated by finite numbers of Zernike polynomials. To show this, we address the convergence of Zernike expansion and their derivatives specifically for odd-order surfaces. Section 2 defines Zernike polynomials and describes their derivation formulae. We present the expansion coefficients of odd-order surfaces and an estimation of the decreasing speed of Zernike coefficients mainly for odd-order surfaces. In Sec. 3, we originally give a proof of the convergence for Zernike expansion in surface shapes and slopes for odd-order surfaces using the result of Sec. 2. In Sec. 4, we practically show the effectiveness of the expansion formula by its numerical estimation. By applying the result to lens design of a Schmidt surface, we demonstrate that both the corrector shape and slope are fully approximated by a finite number of Zernike polynomials.

In addition, the resultant approximation of an odd-order surface is not an example of Taylor expansion, which will be briefly explained in Sec. 2.1. Thus, the result of our method does not contradict the consequences of Ref. 11, which describes the impossibility of Taylor expansion of odd-order surfaces.

To express aspherical surfaces, expansion into power series, Zernike expansion, Qcon polynomials, and so on have been used. The method described in this paper allows us to analyze odd-order surfaces by Zernike polynomials and power series. Expressing odd-order surfaces by Qcon surfaces is our future issue.

2. Definition of Zernike Polynomials and Their Basic Properties

2.1.

Definition of Rotationally Symmetrical Zernike Polynomials

Rotationally symmetrical Zernike polynomials $Z_{2 n}^{0} (r)$ are given as¹³

Eq. (1)

Z_{2 n}^{0} (r) = Q_{n} (t) = \frac{{(- 1)}^{n}}{n!} \frac{d^{n}}{d t^{n}} t^{n} {(1 - t)}^{n},

where

t = r^{2}

,

r

is the normalized radial coordinate, and

n

is a non-negative integer.

This definition of the right-side differential is called the Rodrigues formulae.¹⁵ Explicitly, one obtains the first six polynomials as

Eq. (2)

Q_{0} (t) = 1, Q_{1} (t) = 2 t - 1, Q_{2} (t) = 6 t^{2} - 6 t + 1, Q_{3} (t) = 20 t^{3} - 30 t^{2} + 12 t - 1, Q_{4} (t) = 70 t^{4} - 140 t^{3} + 90 t^{2} - 20 t + 1, Q_{5} (t) = 252 t^{5} - 630 t^{4} + 560 t^{3} - 210 t^{2} + 30 t - 1 .

Their orthogonality is expressed as

Eq. (3)

\int_{0}^{1} Q_{n} (t) Q_{m} (t) d t = \frac{1}{2 n + 1} δ_{m, n},

where

δ_{m, n}

is Kronecker’s delta.

Since ${Q_{n} (t)}$ is a complete orthogonal set on the unit interval of [0, 1], any $L^{2}$ function $f (t)$ is expanded by ${Q_{n} (t)}$ . Namely,

Eq. (4)

f (t) = a_{0} Q_{0} (t) + a_{1} Q_{1} (t) + a_{2} Q_{2} (t) + \dots,

where the coefficients are

Eq. (5)

a_{n} = (2 n + 1) \int_{0}^{1} f (t) Q_{n} (t) d t .

The expansion by orthogonal polynomials such as Eq. (4) is also expressed as a polynomial $f (t) = c_{0} + c_{1} t + c_{2} t^{2} + \dots c_{N} t^{N}$ by taking the terms up to $N$ and rearranging the order of terms by Eq. (2). Although this formal expansion does not give an example of Taylor expansion, this is actually a fair approximation of the original function $f (t)$ . (If the expression $f (t) = c_{0} + c_{1} t + c_{2} t^{2} + \dots$ is a Taylor expansion, the coefficients $c_{k}$ s do not depend on the number of terms. However, in this case, the constant term $c_{0} = \sum_{k = 0}^{N} {(- 1)}^{k} a_{k}$ obviously depends on the number of terms $N$ . Hence this is not an example of Taylor expansion.) This fact will be demonstrated in Sec. 3.3 Eq. (31) using $f (t) = t^{3 / 2}$ as an example.

2.2.

Theoretical Estimation of the Decreasing Speed of Zernike Coefficients

For a generic monomial $f (t) = r^{2 α} = t^{α} (α > 0)$ , the coefficients of Eq. (5) are expressed by the gamma function.¹⁸ That is,

Eq. (6)

a_{n} = \frac{(2 n + 1) Γ {(α + 1)}^{2}}{Γ (α + n + 2) Γ (α - n + 1)} .

Hence, the expansion becomes

Eq. (7)

f (t) = t^{α} = \sum_{n = 0}^{\infty} \frac{(2 n + 1) Γ {(α + 1)}^{2}}{Γ (α + n + 2) Γ (α - n + 1)} Q_{n} (t) .

Braat and Janssen¹⁸ gave an equivalent formulation by analytical extension of an even-order term. We present another derivation of this expression by use of the beta function $B (p, q)$ . Substituting $f (t) = t^{α}$ into Eq. (5) and integrating by part,

Eq. (8)

a_{n} = (2 n + 1) \frac{{(- 1)}^{n}}{n!} \int_{0}^{1} t^{α} \frac{d^{n}}{d t^{n}} t^{n} {(1 - t)}^{n} d t = \frac{2 n + 1}{n!} α (α - 1) \dots (α - n + 1) \int_{0}^{1} t^{α - n} t^{n} {(1 - t)}^{n} d t = \frac{2 n + 1}{n!} \frac{Γ (α + 1)}{Γ (α - n + 1)} B (α + 1, n + 1) = (2 n + 1) \frac{Γ {(α + 1)}^{2}}{Γ (α - n + 1) Γ (α + n + 2)} .

The last equality is derived by

B (p, q) = Γ (p) Γ (q) / Γ (p + q)

.

In order to evaluate the convergence, we estimate the decreasing speed of $a_{n}$ . For representing an odd-order aspherical surface, let $α$ be a noninteger and $α \geq 1 / 2$ . The minimum $α = 1 / 2$ corresponds to the first-order surface or the cone. By Euler’s reflection formula $Γ (z) Γ (1 - z) = \frac{\sin π z}{π}$ ,

Eq. (9)

Γ (α - n + 1) = \frac{\sin π (n - α)}{π} \frac{1}{Γ (n - α)} .

By the definition of the gamma function of

Γ (z + 1) = z Γ (z)

,

Eq. (10)

Γ (α + n + 2) = (α + n + 1) (α + n) Γ (n + α) .

The denominator of Eq. (6) is

Eq. (11)

Γ (α + n + 2) Γ (α - n + 1) = \frac{\sin π (n - α)}{π} (α + n + 1) (α + n) \frac{Γ (n + α)}{Γ (n - α)} .

Since

α \geq 1 / 2

, the inequality

n + α > n + α - 1 \geq n - α

holds. Since

Γ (x) > Γ (y)

for positive numbers,

x > y > 0

. Then,

Γ (n + α)

is evaluated as

Eq. (12)

Γ (n + α) = (n + α - 1) Γ (n + α - 1) > (n + α - 1) Γ (n - α) where n > α .

Therefore,

Eq. (13)

| Γ (α + n + 2) Γ (α - n + 1) | > | \frac{\sin π (n - α)}{π} | (α + n + 1) (α + n) (n + α - 1) .

Thus,

Eq. (14)

| a_{n} | = | \frac{(2 n + 1) Γ {(α + 1)}^{2}}{Γ (α + n + 2) Γ (α - n + 1)} | < M | \frac{(2 n + 1)}{(α + n + 1) (α + n) (α + n - 1)} |,

where

M = | \frac{\sin π (n - α)}{π} |

is a positive number. Hence,

| a_{n} |

decreases faster than or equally to

n^{- 2}

.

2.3.

Another Representation of Coefficients for Odd-Order Surfaces and a More Strict Theoretical Estimation of Their Decreasing Speed

We present a more detailed formulation for the odd-order case $(r^{2 k - 1} = t^{\frac{2 k - 1}{2}})$ . As far as we know, this concrete estimation has not been studied. If $n \geq k$ , substituting the equalities of the gamma function of half-integer¹⁵

Eq. (15)

Γ (\frac{1}{2} + m) = \frac{(2 m)!}{4^{m} m!} \sqrt{π}, Γ (\frac{1}{2} - m) = \frac{{(- 4)}^{m} m!}{(2 m)!} \sqrt{π}

into Eq. (6), we obtain

Eq. (16)

a_{n} = 4^{n + k + 1} \frac{(n + k + 1)!}{{2 (n + k + 1)}!} \frac{{2 (n - k)}!}{{(- 4)}^{n - k} (n - k)!} (2 n + 1) {(\frac{(2 k)!}{4^{k} k!})}^{2} = {(- 1)}^{n - k} 4^{2 k + 1} \frac{(n + k + 1) \dots (n - k + 1)}{(2 n + 2 k + 2) \dots (2 n - 2 k + 1)} (2 n + 1) {(\frac{(2 k)!}{4^{k} k!})}^{2} .

Since Eq. (16) includes

2 k + 1

terms of

n

in the numerator and

4 k + 1

terms of

n

in the denominator,

a_{n}

decreases as fast as

n^{- 2 k}

, which is faster than the estimation described in Eq. (14).

For the third-order surface, substituting $k = 2$ into Eq. (16), we obtain an explicit expression as

Eq. (17)

a_{n} = {(- 1)}^{n} (2 n + 1) \frac{9}{16} \times \frac{2^{10} (n + 3) (n + 2) (n + 1) n (n - 1)}{(2 n + 6) (2 n + 5) (2 n + 4) (2 n + 3) (2 n + 2) (2 n + 1) 2 n (2 n - 1) (2 n - 2) (2 n - 3)} .

Thus, the coefficient

a_{n}

decreases as fast as

n^{- 4}

.

Furthermore for the first-order surface, which describes the cone, substituting $k = 1$ into Eq. (16), we obtain an explicit expression as

Eq. (18)

a_{n} = {(- 1)}^{n - 1} (2 n + 1) \times \frac{(n + 2) (n + 1) n (n - 1)}{(2 n + 4) (2 n + 3) (2 n + 2) (2 n + 1) 2 n (2 n - 1)} .

This decreases as fast as $n^{- 2}$ , which corresponds to the minimum speed of decrease described in Eq. (14).

2.4.

Derivatives of Zernike Polynomials

Derivatives of Zernike polynomials have been already thoroughly discussed.¹⁹^–²¹ In this paper, we present a simple formulation for the rotationally symmetrical case.

According to Appendix A, it is enough to discuss the derivatives by $t = r^{2}$ instead of the radial coordinate $r$ . Since $Q_{n} (t)$ is a polynomial to the $n$ ’th power, its derivative ${Q_{n}}^{'} (t)$ must be described as a linear combination of ${Q_{k} (t)}_{k = 0}^{n - 1}$ . Thus,

Eq. (19)

Q_{n}^{'} (t) = c_{0} Q_{0} (t) + c_{1} Q_{1} (t) + \dots c_{n - 1} Q_{n - 1} (t) .

By the orthogonality shown in Eq. (3), the coefficient $c_{m}$ is

Eq. (20)

c_{m} = (2 m + 1) \int_{0}^{1} {Q_{n}}^{'} (t) Q_{m} (t) d t = (2 m + 1) {[Q_{n} (t) Q_{m} (t)]}_{0}^{1} - (2 m + 1) \int_{0}^{1} Q_{n} (t) Q_{m}^{'} (t) d t .

Since $Q_{n} (0) = {(- 1)}^{n}$ and $Q_{n} (1) = 1$ , the first term of the right side is $(2 m + 1) {{(- 1)}^{n + m} - 1}$ . The second term is zero because the order of ${Q_{m}}^{'} (t)$ is less than $Q_{n} (t)$ . Thus

Eq. (21)

c_{m} = (2 m + 1) {{(- 1)}^{n + m} - 1} .

Therefore, Eq. (19) is represented as

Eq. (22)

Q_{2 k + 1}^{'} (t) = 2 \sum_{j = 0}^{k} (4 j + 1) Q_{2 j} (t), Q_{2 k}^{'} (t) = 2 \sum_{j = 1}^{k} (4 j - 1) Q_{2 j - 1} (t) .

For example, the first six terms of this derivation are definitely shown as

Eq. (23)

Q_{0}^{'} = 0, Q_{1}^{'} = 2 Q_{0}, Q_{2}^{'} = 6 Q_{1}, Q_{3}^{'} = 2 Q_{0} + 10 Q_{2}, Q_{4}^{'} = 6 Q_{1} + 14 Q_{3}, Q_{5}^{'} = 2 Q_{0} + 10 Q_{2} + 18 Q_{4} .

2.5.

Sums of Expansion Coefficients

When we consider $f (t) = t^{α}$ , then $f (0) = 0$ and $f (1) = 1$ . Since $Q_{n} (0) = {(- 1)}^{n}$ and $Q_{n} (1) = 1$ , when substituting $t = 0, 1$ into Eq. (4), the following relations hold:

Eq. (24)

f (0) = a_{0} - a_{1} + a_{2} - \dots = \sum_{n = 1}^{\infty} {(- 1)}^{n} a_{n} = 0, f (1) = a_{0} + a_{1} + a_{2} + \dots = \sum_{n = 1}^{\infty} a_{n} = 1 .

By adding and subtracting the two equations of Eq. (24)

Eq. (25)

\sum_{k = 1}^{\infty} a_{2 k} = \sum_{k = 1}^{\infty} a_{2 k - 1} = \frac{1}{2} .

3. Uniform Convergence of Zernike Expansion in Shape and Slope for Odd-Order Aspherical Surface

In general, let us consider a function $g (t)$ and its expansion $g (t) = \sum h_{n} W_{n} (t)$ , where ${W_{n} (t)}$ is an arbitrary sequence of functions. Even though $g (t) = \sum h_{n} W_{n} (t)$ converges uniformly, its termwise differential $\sum h_{n} {W_{n}}^{'} (t)$ does not necessarily converge to $g^{'} (t)$ .

As shown in Sec. 2, when $f (t)$ represents an odd-order surface shape, its Zernike expansion $f (t) = \sum a_{n} Q_{n} (t)$ contains an infinite number of terms. In optical ray tracing, since the direction of the exit ray is determined by the surface slope, it is expected that $\sum a_{n} {Q_{n}}^{'} (t)$ converges to the surface slope $f^{'} (t)$ . In this section, we prove the convergence of both $f (t) = \sum a_{n} Q_{n} (t)$ and $f^{'} (t) = \sum a_{n} {Q_{n}}^{'} (t)$ .

3.1.

Convergence of Surface Shapes

We present a proof of the uniform convergence of Eq. (7). For an odd-order surface, $| a_{n} |$ decreases faster than or equally to $n^{- 2}$ . Since $| Q_{n} (t) | \leq 1$ , there exists a positive number $C$ and an integer $N$ such that

Eq. (26)

| a_{0} Q_{0} (t) + a_{1} Q_{1} (t) + a_{2} Q_{2} (t) \dots | \leq \sum_{n = 0}^{\infty} | a_{n} | \leq \sum_{n = 0}^{N} | a_{n} | + C \sum_{n = N + 1}^{\infty} \frac{1}{n^{2}} < \sum_{n = 0}^{N} | a_{n} | + \frac{C π^{2}}{6} .

Thus, we conclude that the convergence of the expansion Eq. (7) for any odd-order surface is uniform and absolute by the Weierstrass M-test.¹⁵

3.2.

Convergence of Surface Slopes

To discuss the convergence of the surface slope of an odd-order surface, let us consider the monomial $f (t) = t^{\frac{2 k - 1}{2}} = r^{2 k - 1}$ on the interval [0, 1]. To guarantee the continuity of the derivative, we suppose that the parameter $k \geq 2$ (higher than the third-order surface). Since the case $k = 1$ corresponds to the first-order surface or the cone, the original function $f (t) = \sqrt{t}$ has a singular point at the origin. Although this case is considered as an exceptional case in this paper, similar discussion will be applied to excluding the origin in Sec. 3.4.

Since the derivative $f^{'} (t)$ is also continuous, this function is also expanded as

Eq. (27)

f^{'} (t) = b_{0} Q_{o} (t) + b_{1} Q_{1} (t) + b_{2} Q_{2} (t) + \dots = \sum_{k = 0}^{\infty} b_{n} Q_{n} (t),

By the analogy of Fourier analysis of derivatives, the convergence in slopes is proven if the following two conditions are satisfied:

Condition 1. Convergence of slope: $\sum_{k = 0}^{\infty} b_{n} Q_{n} (t)$ converges to $f^{'} (t)$ .
Condition 2. Consistency of termwise differential: the equality
Eq. (28)
$\sum_{k = 0}^{\infty} a_{n} Q_{n}^{'} (t) = \sum_{k = 0}^{\infty} b_{n} Q_{n} (t)$

holds.

Since $f^{'} (t)$ is also the odd-order surface of order $(2 k - 3) / 2$ , the proof of condition 1 is the direct conclusion of the discussion of Sec. 3.1.

To show condition 2, we give a straightforward calculation of the left side of Eq. (28). According to the Zernike derivatives formulae Eq. (23), by differentiating Eq. (4) formally

Eq. (29)

{(a_{0} Q_{0} (t) + a_{1} Q_{1} (t) + a_{2} Q_{2} (t) + \dots)}^{'} = 2 (a_{1} + a_{3} + a_{5} + \dots) Q_{0} + 6 (a_{2} + a_{4} + a_{6} + \dots) Q_{1} + 10 (a_{3} + a_{5} + a_{7} + \dots) Q_{2} + 14 (a_{4} + a_{6} + a_{8} + \dots) Q_{3} + 18 (a_{5} + a_{7} + a_{9} + \dots) Q_{4} \dots

By the procedure described in Appendix B and Eq. (25), the coefficients

b_{k}

are

Eq. (30)

b_{0} = 1 = 2 (a_{1} + a_{3} + a_{5} + \dots), b_{1} = 6 (\frac{1}{2} - a_{0}) = 6 (a_{2} + a_{4} + a_{6} + \dots), b_{2} = 10 (\frac{1}{2} - a_{1}) = 10 (a_{3} + a_{5} + a_{7} + \dots), b_{3} = 14 (\frac{1}{2} - a_{0} - a_{2}) = 14 (a_{4} + a_{6} + a_{8} + \dots), b_{4} = 18 (\frac{1}{2} - a_{1} - a_{3}) = 18 (a_{5} + a_{7} + a_{9} + \dots), b_{5} = 22 (\frac{1}{2} - a_{0} - a_{2} - a_{4}) = 22 (a_{6} + a_{8} + a_{10} + \dots) .

Thus, the right side of Eq. (28) is identical to the left side of this equation.

3.3.

Numerical Confirmation for $f (t) = t^{3 / 2} = r^{3}$

By Eq. (17), we obtain

Eq. (31)

t^{\frac{3}{2}} ≅ \frac{2}{5} Q_{0} (t) + \frac{18}{35} Q_{1} (t) + \frac{2}{21} Q_{2} (t) - \frac{2}{165} Q_{3} (t) + \frac{18}{5005} Q_{4} (t) - \frac{2}{1365} Q_{5} (t) = \frac{4}{2145} + \frac{28}{143} t + \frac{224}{143} t^{2} - \frac{224}{143} t^{3} + \frac{168}{143} t^{4} - \frac{24}{65} t^{5} .

Also by Eq. (30), the coefficients of derivatives are calculated as follows:

Eq. (32)

b_{0} = 1, b_{1} = 6 \times [\frac{1}{2} - a_{0}] = \frac{3}{5}, b_{2} = 10 \times [\frac{1}{2} - a_{1}] = - \frac{1}{7}, b_{3} = 14 \times [\frac{1}{2} - (a_{0} + a_{2})] = \frac{1}{15}, b_{4} = 18 \times [\frac{1}{2} - (a_{1} + a_{3})] = - \frac{3}{77} .

Thus, the expansion of the derivatives is

Eq. (33)

\frac{d}{d t} (t^{3 / 2}) ≅ Q_{0} (t) + \frac{3}{5} Q_{1} (t) - \frac{1}{7} Q_{2} (t) + \frac{1}{15} Q_{3} (t) - \frac{3}{77} Q_{4} (t) .

The fitting errors of Eqs. (31) and (32) are shown in Fig. 1

Fig. 1

Fitting errors in shape and slope for the third order surface. ( $N = 5$ ).

To estimate the maximum error versus the number of Zernike polynomials, we present the list of coefficients up to $n = 15$ in Table 1. Since the maximum errors occur at the origin, the errors can be estimated by $\sum_{k = 0}^{n} {(- 1)}^{n} a_{n}$ for the shape and $\sum_{k = 0}^{n} {(- 1)}^{n} b_{n}$ for the slope. Both of them converge to zero when $n \to \infty$ .

Table 1

Errors in shape and slope for the third order surface (N≤15).

n	an	∑k=0n(−1)kak	bn	∑k=0n(−1)kbk
0	$4.0000 E - 01$	$4.0000 E - 01$	$1.0000 E + 00$	$1.0000 E + 00$
1	$5.1429 E - 01$	$- 1.1429 E - 01$	$6.0000 E - 01$	$4.0000 E - 01$
2	$9.5238 E - 02$	$- 1.9048 E - 02$	$- 1.4286 E - 01$	$2.5714 E - 01$
3	$- 1.2121 E - 02$	$- 6.9264 E - 03$	$6.6667 E - 02$	$1.9048 E - 01$
4	$3.5964 E - 03$	$- 3.3300 E - 03$	$- 3.8961 E - 02$	$1.5152 E - 01$
5	$- 1.4652 E - 03$	$- 1.8648 E - 03$	$2.5641 E - 02$	$1.2587 E - 01$
6	$7.1301 E - 04$	$- 1.1518 E - 03$	$- 1.8182 E - 02$	$1.0769 E - 01$
7	$- 3.8970 E - 04$	$- 7.6209 E - 04$	$1.3575 E - 02$	$9.4118 E - 02$
8	$2.3135 E - 04$	$- 5.3074 E - 04$	$- 1.0526 E - 02$	$8.3591 E - 02$
9	$- 1.4615 E - 04$	$- 3.8459 E - 04$	$8.4034 E - 03$	$7.5188 E - 02$
10	$9.6917 E - 05$	$- 2.8768 E - 04$	$- 6.8650 E - 03$	$6.8323 E - 02$
11	$- 6.6834 E - 05$	$- 2.2084 E - 04$	$5.7143 E - 03$	$6.2609 E - 02$
12	$4.7595 E - 05$	$- 1.7325 E - 04$	$- 4.8309 E - 03$	$5.7778 E - 02$
13	$- 3.4821 E - 05$	$- 1.3843 E - 04$	$4.1379 E - 03$	$5.3640 E - 02$
14	$2.6067 E - 05$	$- 1.1236 E - 04$	$- 3.5842 E - 03$	$5.0056 E - 02$
15	$- 1.9903 E - 05$	$- 9.2455 E - 05$	$3.1348 E - 03$	$4.6921 E - 02$

Since $f (t) = t^{3 / 2}$ and $f^{'} (t) = (3 / 2) t^{1 / 2}$ , the decreasing speed of the coefficients is as fast as $n^{- 4}$ in shapes by Eq. (17) and $n^{- 2}$ in slopes by Eq. (18), respectively. Thus, Zernike expansions both in shape and in slope converge to the original functions. The behaviors of the maximum errors as functions of the number of Zernike terms are shown in Fig. 2.

Fig. 2

Convergence of the maximum errors for $f (t) = t^{3 / 2}$ .

3.4.

Approximation for the Cone $f (t) = t^{1 / 2} = r$

The cone is not an example of smooth odd-order functions. However, as discussed in Secs. 2.3 and 3.1, the Zernike expansion of the cone converges uniformly to $f (t) = t^{1 / 2}$ . Thus, the cone ure is approximated enough by the Zernike expansion including the origin. Moreover, the similar discussions in Secs. 3.2 and 3.3 can be applied to the cone except at the singular point or the origin. However, the Zernike expansion of $f^{'} (t)$ diverges at the origin because its expansion gives a diverging series.

Actually, by Eq. (6), the expansion coefficients of $f^{'} (t)$ are

Eq. (34)

b_{n} = \frac{1 / 2 (2 n + 1) Γ {(1 / 2)}^{2}}{Γ (3 / 2 + n) Γ (1 / 2 - n)} .

By Eq. (15), the denominator is

Eq. (35)

Γ (\frac{3}{2} + n) Γ (\frac{1}{2} - n) = \frac{(2 n + 2)!}{4^{n + 1} (n + 1)!} \sqrt{π} \times \frac{{(- 4)}^{n} n!}{(2 n)!} \sqrt{π} = \frac{{(- 1)}^{n} (2 n + 2) (2 n + 1) π}{4 (n + 1)} = {(- 1)}^{n} \frac{2 n + 1}{2} π .

Since

Γ (1 / 2) = \sqrt{π}

,

Eq. (36)

b_{n} = {(- 1)}^{n} .

Thus,

Eq. (37)

\sum_{k = 0}^{\infty} b_{n} Q_{n} (0) = \sum_{k = 0}^{\infty} {(- 1)}^{n} {(- 1)}^{n} = \sum_{k = 0}^{\infty} 1 = \infty .

Hence, the sum

| \sum_{k = 0}^{\infty} b_{n} Q_{n} (t) |

does not converge at the origin.

Another explanation is possible as follows. Since the cone is described as $f (t) = t^{1 / 2} = | r |$ on the interval [ $- 1, 1$ ], the differential is

Eq. (38)

\frac{d f}{d r} = {\begin{cases} - 1 (r < 0) \\ 1 (r > 0) \end{cases} .

This conclusion intuitively meets the fact that the cone has a singular point at the vertex.

4. Design Example

4.1.

Schmidt Surface

In this section, we evaluate the approximation accuracy of the Zernike expansion and discuss the aberration properties of odd-order surfaces via the optical design of the Schmidt corrector plate. The corrector plate is incorporated into the Schmidt camera and corrects the spherical aberration of the primary mirror. The optical layout of the Schmidt camera is shown in Fig. 3. Even though corrector surfaces are conventionally designed only with even-order terms containing the power,²²^,²³ the design with odd-order terms has been known as the Schmidt surface²⁴.

Fig. 3

Layout of Schmidt camera.

The specifications of the design example are shown in Table 2 and the basic lens data are shown in Table 3. The surface indicated by (*) in Table 3 is the corrector surface. The aperture stop is also placed there. For simplicity, we employ only the third term and even-order terms for the corrector plate. Thus, the aspherical sag of surface (*) is given as follows:

Eq. (39)

z = A_{2} r^{2} + A_{3} r^{3} + A_{4} r^{4} + A_{6} r^{6} + \dots

Note that all optical design and analysis is achieved with Code-V^TM.

Table 2

Design specification.

Features	Requirements
Wavelength	587.56 nm
Focal length	200 mm
$F$ number	3
Half FOV	7 deg
Residual wavefront aberration	$< 0.0156 λ$ RMS
Residual wavefront aberration	( $Strehl ratio > 0.99$ )

Table 3

Basic lens data.

	R	d	Glass
0	Infinity	Infinity
1	Infinity	10	$n = 1.51633$
2(*)	Infinity	200
3	Infinity	200
4	$- 400$	$- 200.000$	Reflect
5	Variable	Variable

4.2.

Design Results and Discussion

We compare the original third-order design with its approximation design by the even-order expansion. To obtain a diffraction-limited design, the third and even orders up to the 16th are needed. By decomposing the third coefficients by Eq. (17) and adding to the original even-order coefficients, we obtain the diffraction-limited design by expanded shape as well, in which even orders up to the 20th are needed. The aspherical coefficients are shown in Table 4.

Table 4

Odd order and approximation coefficients.

Aspherical coefficients	Original design	Approximation
$A 2$		$- 1.009321 E - 05$
$A 3$	$- 2.923781 E - 06$
$A 4$	$2.192910 E - 07$	$- 4.641273 E - 08$
$A 6$	$- 4.867503 E - 10$	$5.052101 E - 10$
$A 8$	$1.023192 E - 12$	$- 2.492068 E - 12$
$A 10$	$- 1.418425 E - 15$	$7.259246 E - 15$
$A 12$	$1.199493 E - 18$	$- 1.314277 E - 17$
$A 14$	$- 5.588738 E - 22$	$1.495951 E - 20$
$A 16$	$1.096804 E - 25$	$- 1.041559 E - 23$
$A 18$		$4.054177 E - 27$
$A 20$		$- 6.760936 E - 31$

The error in shape is shown in Fig. 4. If we use more Zernike terms, the dominant error at the origin will be gradually reduced.

Fig. 4

Fitting error in shape.

The peak-to-valley error in Fig. 4 is $3.88 \times 10^{- 5} mm$ , which corresponds to $0.034 λ$ . The comparison of wavefront map and the RMS wavefront $Δ W$ of the two designs is shown in Fig. 5. Note that the RMS wavefront is calculated without considering the central obscuration.

Fig. 5

Design comparison.

The expanded shape is a fair approximation of the original surface shape. From the viewpoint of practical optical design, the wavefront aberrations in both designs are almost identical. Thus, the third-order surface is almost completely approximated by a finite sum of even-order terms.

In addition, since there is the $A_{2}$ term in the expanded shape, the third-order surface acts as an aspherical surface containing the power, which is necessary in conventional designs. Although the power term induces slight change in the paraxial relation, the third-order surface does not affect the paraxial relation formally. Consequently, we can conclude that this property is one reason for the effectiveness of odd-order surfaces. Even though we have proven that the third-order aspherical surface is fairly approximated by a finite sum of even-order terms, we infer that there are some other reasons for the effectiveness. This is our future issue.

5. Conclusion

In this paper, we have shown that odd-order surfaces are fairly approximated by a finite number of Zernike polynomials. Rearranging the order of monomials, odd-order surfaces are fully approximated by a finite number of even-order power series.

To show this, we have proven the uniform convergence of Zernike expansion of odd-order surfaces both in surface shapes and slopes. By estimating the decreasing speed of the expansion coefficients, any odd-order surface is precisely approximated to be an optical surface by a finite number of Zernike polynomials. In other words, odd-order surfaces are fairly approximated by ordinary even-order aspherical surfaces. For the first-order surface, the Zernike expansion approximates the shape. However, the expansion of the slope does not make sense because the coefficients give a diverging series.

We have demonstrated the result by a design of a Schmidt surface and confirmed that the effect of the third-order surface can be fully expressed by the 20th even-order surface. The power term is necessary in classical design of Schmidt cameras. However, in the calculation of paraxial quantity, odd-order surfaces do not affect values such as focal lengths, paraxial magnifications, and so on. This is one reason why odd-order surfaces are effective in optical design. Since we infer that there are some other reasons for the odd order’s effectiveness, this will be our future problem.

Our analysis method is useful for not only analysis of aberration characteristics of aspherical surfaces but also optical measurement especially in the accuracy estimation of approximation of generic surfaces by finite number of terms. Analyzing another expression of odd-order surfaces such as Qcon polynomials is a future issue.

Appendices

Appendix A:

Conversion of the Radial Coordinate in the Derivative Formulae

In the definition of Zernike polynomials Eq. (1) and expansion Eq. (4), the converted coordinate $t = r^{2}$ is used instead of the original coordinate $r$ . Considering that the surface is given by $z = F (x, y)$ , due to the law of refraction, the direction of the exit ray is determined by the surface slopes $\partial F / \partial x$ and $\partial F / \partial y$ . If $F (x, y)$ is rotationally symmetric, $F (x, y)$ should be a function of $t = r^{2}$ , i.e., $F (x, y) = f (t)$ .

Consider the Zernike expansion of $f (t) = \sum a_{n} Q_{n} (t)$ , the slopes are rewritten as

Eq. (40)

\frac{\partial F}{\partial x} = \frac{\partial t}{\partial x} \frac{d f}{d t} = 2 x {(\sum a_{n} Q_{n} (t))}^{'}, \frac{\partial F}{\partial y} = \frac{\partial t}{\partial y} \frac{d f}{d t} = 2 y {(\sum a_{n} Q_{n} (t))}^{'} .

Thus, in the context of Zernike expansion and its derivatives, it is enough to prove the equality of the termwise derivative ${(\sum a_{n} Q_{n} (t))}^{'} = \sum a_{n} {Q_{n}}^{'} (t)$ and its convergence.

Appendix B:

Derivation of the Coefficients $b_{k}$

We present a straightforward derivation of the expansion coefficients $b_{k}$ . By differentiating Eq. (4) formally, we obtain

Eq. (41)

{(a_{0} Q_{0} (t) + a_{1} Q_{1} (t) + a_{2} Q_{2} (t) + \dots)}^{'} = 2 \sum_{n = 0}^{\infty} [{2 \times 2 n + 1} Q_{2 n} (t) \sum_{k = n + 1}^{\infty} a_{2 k - 1} + {2 (2 n + 1) + 1} Q_{2 n + 1} (t) \sum_{k = n + 1}^{\infty} a_{2 k}] .

At this stage, the convergence of Eq. (41) is not guaranteed.

The coefficients $b_{m}$ in Eq. (27) are calculated by the orthogonality of $Q_{m} (t)$ ,

Eq. (42)

b_{m} = (2 m + 1) \int_{0}^{1} f^{'} (t) Q_{m} (t) d t = (2 m + 1) [f (t) Q_{m} (t) |_{0}^{1} - \int_{0}^{1} f (t) {Q_{m}}^{'} (t) d t] = (2 m + 1) [1 - \int_{0}^{1} f (t) Q_{m}^{'} (t) d t] .

This calculation differs according to the parity of $m$ . If $m$ is even,

Eq. (43)

b_{2 n} = (4 n + 1) [1 - \int_{0}^{1} f (t) {Q_{2 n}}^{'} (t) d t] = (4 n + 1) [1 - 2 \int_{0}^{1} f (t) \sum_{k = 1}^{n} (4 k - 1) Q_{2 k - 1} (t) d t] .

Since

\int_{0}^{1} f (t) (4 k - 1) Q_{2 k - 1} (t) d t = {2 (2 k - 1) + 1} \int_{0}^{1} f (t) Q_{2 k - 1} (t) d t = a_{2 k - 1},

the expansion coefficient $b_{2 n}$ is rewritten as

Eq. (44)

b_{2 n} = (4 n + 1) [1 - 2 \sum_{k = 1}^{n} a_{2 k - 1}] .

By Eq. (25),

Eq. (45)

b_{2 n} = 2 (4 n + 1) \sum_{k = n + 1}^{\infty} a_{2 k - 1} .

This is identical to the coefficient of $Q_{2 n} (t)$ in Eq. (41). By the same procedure, the case of $m = odd$ is proven, namely,

Eq. (46)

b_{2 n + 1} = 2 (4 n + 3) \sum_{k = n + 1}^{\infty} a_{2 k} .

Acknowledgments

The authors would like to thank Mr. Makoto Fujino and Dr. Akira Takada for their valuable suggestions. Also, the authors would like to thank Mr. Shota Obigane, Mr. Yusuke Araki, Mr. Takashi Toriyabe, and Mr. Takayuki Nomura for their useful optical designs.

References

1.

I. A. Neil, “Optical design dependence on technology development,” Opt. Eng., 50 (12), 121706 (2011). http://dx.doi.org/10.1117/1.3631888 Google Scholar

2.

G. W. Forbes and C. P. Brophy, “Asphere, O Asphere, how shall we describe thee?,” Proc. SPIE, 7100 710002 (2008). http://dx.doi.org/10.1117/12.797770 PSISDG 0277-786X Google Scholar

3.

K. P. Thompson et al., “The Forbes polynomial surface: the road forward; a return to manufacturable aspheres,” SPIE Optifab Panel Discussion, (2011) http://spie.org/Documents/ConferencesExhibitions/2011_Optifab_panel_final.pdf May 2016). Google Scholar

4.

A. Yabe, “Sensitivity control to surface irregularity,” Proc. SPIE, 6342 634225 (2006) http://dx.doi.org/10.1117/12.692246 Google Scholar

5.

G. Shu et al., “Trapped ion imaging with a high numerical aperture spherical mirror,” J. Phys. B: At. Mol. Opt. Phys., 42 154005 (2009). http://dx.doi.org/10.1088/0953-4075/42/15/154005 JPAPEH 0953-4075 Google Scholar

6.

W. T. Plummer, “Viewfinder for a reflex camera,” U. S. Patent No. 3,810,221 (1974).

7.

G. Carsten, “Aberration-corrected microscope,” (2012). Google Scholar

8.

K. Sato, “Wide-angle zoom lens including at least one aspheric lens surface,” U.S. Patent No. 6,982,834 (2006).

9.

M. Amano, “Projection optical system and projection display apparatus,” U.S. Patent No. 8,967,812 (2015).

10.

J. P. McGuire, “EUV condenser with non-imaging optics,” U.S. Patent No. 7,046,338 (2006).

11.

M. Shibuya et al., “Theoretical investigation of the meaning of odd-order aspherical surface and numerical confirmation of effectiveness in rotational-symmetric but off-axis optics,” Opt. Eng., 49 (7), 073003 (2010). http://dx.doi.org/10.1117/1.3461991 Google Scholar

12.

M. Tanikawa et al., “Effectiveness of odd order aspherical surface,” Jpn. J. Opt., 36 (11), 646 –660 (2007). KOGAD5 0389-6625 Google Scholar

13.

M. Born and E. Wolf, Principles of Optics, 2nd (Revised) ed.Cambridge University Press, Cambridge, United Kingdom (1999). Google Scholar

14.

F. Zernike, “Beugungstheorie des schneidenver-fahrens und seiner verbesserten form, der phasenkontrastmethode,” Physica, 1 (7–12), 689 –704 (1934). http://dx.doi.org/10.1016/S0031-8914(34)80259-5 Google Scholar

15.

G. Arfken and H. Weber, Mathematical Methods for Physicists, 5th ed.Academic Press, New York (2001). Google Scholar

16.

A. B. Bhatia and E. Wolf, “On the circle polynomials of Zernike and related orthogonal sets,” Math. Proc. Cambridge Philos. Soc., 50 40 –48 (1954). http://dx.doi.org/10.1017/S0305004100029066 PCPSA4 0008-1981 Google Scholar

17.

G. Szegö, Orthogonal Polynomials, 4th ed.American Mathematical Society, Providence, Rhode Island (1975). Google Scholar

18.

J. J. M. Braat and A. J. E. M. Janssen, “Double Zernike expansion of the optical aberration function from its power series expansion,” J. Opt. Soc. Am. A, 30 (6), 1213 –1222 (2013). http://dx.doi.org/10.1364/JOSAA.30.001213 JOAOD6 0740-3232 Google Scholar

19.

B. R. A. Nijboer, “The diffraction theory of aberrations,” Rijksuniversiteit Groningen, (1942). Google Scholar

20.

A. J. E. M. Janssen, “Zernike expansion of derivatives and Laplacians of the Zernike circle polynomials,” J. Opt. Soc. Am. A, 31 (7), 1604 –1613 (2014). http://dx.doi.org/10.1364/JOSAA.31.001604 JOAOD6 0740-3232 Google Scholar

21.

A. Wünsche, “Generalized Zernike or disk polynomials,” J. Comput. Appl. Math., 174 135 –163 (2005). http://dx.doi.org/10.1016/j.cam.2004.04.004 Google Scholar

22.

M. Laikin, Lens Design, 2nd ed.Dekker, New York (1995). Google Scholar

23.

H. Gross, F. Blechinger and A. Achtner, Handbook of Optical Systems, Wiley-VCH, Weinheim (2008). Google Scholar

24.

“Optics and photonics—preparation of drawings for optical elements and systems—Part 12: aspheric surfaces,” Geneva (2007). Google Scholar

Biography

Takao Tanabe is a lens designer for Topcon Corporation. He received his BS degree in physics from the University of Tokyo in 2001 and his MS degree in mathematics from Tokyo Institute of Technology in 2005. He has designed lithographic optics, medical optics, and laser optics. He is studying optical design, metrology, and lens performance evaluation in Tokyo Polytechnic University. He is a professional engineer, Japan (mechanical). He is a member of SPIE.

Masato Shibuya is a professor in the Department of Media and Image Technology, Faculty of Engineering, Tokyo Polytechnic University. He graduated from Tokyo Institute of Technology in 1977 with a master’s degree in physics and joined Nikon Corporation. He received his PhD from the University of Tokyo in 1996. He joined Tokyo Polytechnic University in 2001, and his research interests include fundamental lens optics, optical lithography, super-resolution, and phase-shifting mask. He is a SPIE fellow.

Kazuhisa Maehara is a part-time lecturer in the center of Basic Education and Research, Tokyo Polytechnic University. He graduated from the University of Tokyo in 1974 with a master’s degree in mathematics. He joined Tokyo Polytechnic University in 1976, and he received his PhD from the University of Tokyo in 1984. He was an associate professor at Tokyo Polytechnic University until 2011. His research interest is classification of algebraic varieties from the viewpoint of algebraic varieties.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

Takao Tanabe, Masato Shibuya, and Kazuhisa Maehara "Convergence and differentiation of Zernike expansion: application for an analysis of odd-order surfaces," Optical Engineering 55(3), 035101 (1 March 2016). https://doi.org/10.1117/1.OE.55.3.035101

Published: 1 March 2016

Access the abstract

JOURNAL ARTICLE
9 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 4 scholarly publications.

Explore citations on Lens.org

KEYWORDS

Zernike polynomials

Optical design

Optical engineering

Cameras

Wavefront aberrations

Aspheric lenses

Wavefronts

1.

Introduction

2.

Definition of Zernike Polynomials and Their Basic Properties

2.1.

Definition of Rotationally Symmetrical Zernike Polynomials

Eq. (1)

Eq. (2)

Eq. (3)

Eq. (4)

Eq. (5)

2.2.

Theoretical Estimation of the Decreasing Speed of Zernike Coefficients

Eq. (6)

Eq. (7)

Eq. (8)

Eq. (9)

Eq. (10)

Eq. (11)

Eq. (12)

Eq. (13)

Eq. (14)

2.3.

Another Representation of Coefficients for Odd-Order Surfaces and a More Strict Theoretical Estimation of Their Decreasing Speed

Eq. (15)

Eq. (16)

Eq. (17)

Eq. (18)

2.4.

Derivatives of Zernike Polynomials

Eq. (19)

Eq. (20)

Eq. (21)

Eq. (22)

Eq. (23)

2.5.

Sums of Expansion Coefficients

Eq. (24)

Eq. (25)

3.

Uniform Convergence of Zernike Expansion in Shape and Slope for Odd-Order Aspherical Surface

3.1.

Convergence of Surface Shapes

Eq. (26)

3.2.

Convergence of Surface Slopes

Eq. (27)

Eq. (28)

Eq. (29)

Eq. (30)

3.3.

Numerical Confirmation for f(t)=t3/2=r3

Eq. (31)

Eq. (32)

Eq. (33)

Fig. 1

Table 1

Fig. 2

3.4.

Approximation for the Cone f(t)=t1/2=r

Eq. (34)

Eq. (35)

Eq. (36)

Eq. (37)

Eq. (38)

4.

Design Example

4.1.

Schmidt Surface

Fig. 3

Eq. (39)

Table 2

Table 3

4.2.

Design Results and Discussion

Table 4

Fig. 4

Fig. 5

5.

Conclusion

Numerical Confirmation for $f (t) = t^{3 / 2} = r^{3}$

Approximation for the Cone $f (t) = t^{1 / 2} = r$

Derivation of the Coefficients $b_{k}$