Optical Design and Engineering

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

[+] Author Affiliations
Takao Tanabe

Tokyo Polytechnic University, Faculty of Engineering, 1583, Iiyama, Atsugi-shi, Kanagawa 243-0297, Japan

Topcon Corporation, 75-1, Hasunuma-cho, Itabashi-ku, Tokyo 174-8580, Japan

Masato Shibuya, Kazuhisa Maehara

Tokyo Polytechnic University, Faculty of Engineering, 1583, Iiyama, Atsugi-shi, Kanagawa 243-0297, Japan

Opt. Eng. 55(3), 035101 (Mar 01, 2016). doi:10.1117/1.OE.55.3.035101
History: Received December 1, 2015; Accepted January 26, 2016
Text Size: A A A

Open Access Open Access

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.

Figures in this Article

New types of aspheric surfaces and their mathematical expressions have been proposed in optical design.15 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,6 microscopes,7 camera lenses,8 projection display optics,9 EUV optical systems,1012 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.11 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.11,12 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 Zernike13,14 as the eigenfunctions of a rotationally symmetrical and self-adjoint differential equation on the unit disk. The orthogonality and the L2-completeness of Zernike polynomials are direct consequences of this definition.15 Another definition was given by Bhatia and Wolf13,16 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ö,17 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 Janssen18 gave an example of explicit Zernike expansion coefficients by factorial functions or the gamma function. For derivatives of Zernike polynomials, many researchers such as Nijboer19 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. Janssen20 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ünsche21 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.

Definition of Rotationally Symmetrical Zernike Polynomials

Rotationally symmetrical Zernike polynomials Z2n0(r)   are given as13Display Formula

Z2n0(r)=Qn(t)=(1)nn!dndtntn(1t)n,(1)
where t=r2, 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.15 Explicitly, one obtains the first six polynomials as Display Formula

Q0(t)=1,Q1(t)=2t1,Q2(t)=6t26t+1,Q3(t)=20t330t2+12t1,Q4(t)=70t4140t3+90t220t+1,Q5(t)=252t5630t4+560t3210t2+30t1.(2)
Their orthogonality is expressed as Display Formula
01Qn(t)Qm(t)dt=12n+1δm,n,(3)
where δm,n is Kronecker’s delta.

Since {Qn(t)} is a complete orthogonal set on the unit interval of [0, 1], any L2 function f(t) is expanded by {Qn(t)}. Namely, Display Formula

f(t)=a0Q0(t)+a1Q1(t)+a2Q2(t)+,(4)
where the coefficients are Display Formula
an=(2n+1)01f(t)Qn(t)dt.(5)

The expansion by orthogonal polynomials such as Eq. (4) is also expressed as a polynomial f(t)=c0+c1t+c2t2+cNtN 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)=c0+c1t+c2t2+ is a Taylor expansion, the coefficients ck s do not depend on the number of terms. However, in this case, the constant term c0=k=0N(1)kak 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)=t3/2 as an example.

Theoretical Estimation of the Decreasing Speed of Zernike Coefficients

For a generic monomial f(t)=r2α=tα  (α>0), the coefficients of Eq. (5) are expressed by the gamma function.18 That is, Display Formula

an=(2n+1)Γ(α+1)2Γ(α+n+2)Γ(αn+1).(6)
Hence, the expansion becomes Display Formula
f(t)=tα=n=0(2n+1)Γ(α+1)2Γ(α+n+2)Γ(αn+1)Qn(t).(7)

Braat and Janssen18 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, Display Formula

an=(2n+1)(1)nn!01tαdndtntn(1t)ndt=2n+1n!α(α1)(αn+1)01tαntn(1t)ndt=2n+1n!Γ(α+1)Γ(αn+1)B(α+1,n+1)=(2n+1)Γ(α+1)2Γ(αn+1)Γ(α+n+2).(8)
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 an. For representing an odd-order aspherical surface, let α be a noninteger and α1/2. The minimum α=1/2 corresponds to the first-order surface or the cone. By Euler’s reflection formula Γ(z)Γ(1z)=sinπzπ, Display Formula

Γ(αn+1)=  sinπ(nα)π1Γ(nα).(9)
By the definition of the gamma function of Γ(z+1)=zΓ(z), Display Formula
Γ(α+n+2)=(α+n+1)(α+n)Γ(n+α).(10)
The denominator of Eq. (6) is Display Formula
Γ(α+n+2)Γ(αn+1)=sinπ(nα)π(α+n+1)(α+n)Γ(n+α)Γ(nα).(11)
Since α1/2, the inequality n+α>n+α1nα holds. Since Γ(x)>Γ(y) for positive numbers, x>y>0. Then, Γ(n+α) is evaluated as Display Formula
Γ(n+α)=(n+α1)Γ(n+α1)>(n+α1)Γ(nα)where  n>α.(12)
Therefore, Display Formula
|Γ(α+n+2)Γ(αn+1)|>  |sinπ(nα)π|(α+n+1)(α+n)(n+α1).(13)
Thus, Display Formula
|an|=|(2n+1)Γ(α+1)2Γ(α+n+2)Γ(αn+1)|<M|(2n+1)(α+n+1)(α+n)(α+n1)|,(14)
where M=  |sinπ(nα)π| is a positive number. Hence, |an| decreases faster than or equally to n2.

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 (r2k1=t2k12). As far as we know, this concrete estimation has not been studied. If nk, substituting the equalities of the gamma function of half-integer15Display Formula

Γ(12+m)=(2m)!4mm!π,Γ(12m)=(4)mm!(2m)!π(15)
into Eq. (6), we obtain Display Formula
an=4n+k+1(n+k+1)!{2(n+k+1)}!{2(nk)}!(4)nk(nk)!(2n+1)((2k)!4kk!)2=(1)nk42  k+1(n+k+1)(nk+1)(2n+2k+2)(2n2k+1)(2n+1)((2k)!4kk!)2.(16)
Since Eq. (16) includes 2k+1 terms of n in the numerator and 4k+1 terms of n in the denominator, an decreases as fast as n2  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 Display Formula

an=(1)n(2n+1)916×210(n+3)(n+2)(n+1)n(n1)(2n+6)(2n+5)(2n+4)(2n+3)(2n+2)(2n+1)2n(2n1)(2n2)(2n3).(17)
Thus, the coefficient an decreases as fast as n4.

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

an=(1)n1(2n+1)×(n+2)(n+1)n(n1)(2n+4)(2n+3)(2n+2)(2n+1)2n(2n1).(18)

This decreases as fast as n2, which corresponds to the minimum speed of decrease described in Eq. (14).

Derivatives of Zernike Polynomials

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

According to Appendix 1, it is enough to discuss the derivatives by t=r2 instead of the radial coordinate r. Since Qn(t)   is a polynomial to the n’th power, its derivative Qn(t) must be described as a linear combination of {Qk(t)}k=0n1. Thus, Display Formula

Qn(t)=c0Q0(t)+c1Q1(t)+cn1Qn1(t).(19)

By the orthogonality shown in Eq. (3), the coefficient cm is Display Formula

cm=(2m+1)01Qn(t)Qm(t)dt=(2m+1)[Qn(t)Qm(t)]01(2m+1)01Qn(t)Qm(t)dt.(20)

Since Qn(0)=(1)n and Qn(1)=1, the first term of the right side is (2m+1){(1)n+m1}. The second term is zero because the order of Qm(t) is less than Qn(t). Thus Display Formula

cm=(2m+1){(1)n+m1}.(21)
Therefore, Eq. (19) is represented as Display Formula
Q2k+1(t)=2j=0k(4j+1)Q2j(t),Q2k(t)=2j=1k(4j1)Q2j1(t).(22)

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

Q0=0,Q1=2Q0,Q2=6Q1,Q3=2Q0+10Q2,Q4=6Q1+14Q3,Q5=2Q0+10Q2+18Q4.(23)

Sums of Expansion Coefficients

When we consider f(t)=tα, then f(0)=0 and f(1)=1. Since Qn(0)=(1)n and Qn(1)=1, when substituting t=0,1 into Eq. (4), the following relations hold: Display Formula

f(0)=a0a1+a2=n=1(1)nan=0,f(1)=a0+a1+a2+=n=1an=1.(24)
By adding and subtracting the two equations of Eq. (24) Display Formula
k=1a2k=k=1a2k1=12.(25)

In general, let us consider a function g(t) and its expansion g(t)=hnWn(t), where {Wn(t)} is an arbitrary sequence of functions. Even though g(t)=hnWn(t) converges uniformly, its termwise differential hnWn(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)=anQn(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 anQn(t) converges to the surface slope f(t). In this section, we prove the convergence of both f(t)=anQn(t) and f(t)=anQn(t).

Convergence of Surface Shapes

We present a proof of the uniform convergence of Eq. (7). For an odd-order surface, |an| decreases faster than or equally to n2. Since |Qn(t)|1, there exists a positive number C and an integer N such that Display Formula

|a0Q0(t)+a1Q1(t)+a2Q2(t)|n=0|an|n=0N|an|+Cn=N+11n2<n=0N|an|+Cπ26.(26)

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.15

Convergence of Surface Slopes

To discuss the convergence of the surface slope of an odd-order surface, let us consider the monomial f(t)=t2k12=r2k1 on the interval [0, 1]. To guarantee the continuity of the derivative, we suppose that the parameter k2 (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)=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 Display Formula

f(t)=b0Qo(t)+b1Q1(t)+b2Q2(t)+=k=0bnQn(t),(27)
By the analogy of Fourier analysis of derivatives, the convergence in slopes is proven if the following two conditions are satisfied:

  • Convergence of slope: k=0bnQn(t) converges to f(t).
  • Consistency of termwise differential: the equality Display Formula
    k=0anQn(t)=k=0bnQn(t)(28)

holds.

Since f(t) is also the odd-order surface of order (2k3)/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 Display Formula

(a0Q0(t)+a1Q1(t)+a2Q2(t)+)=2(a1+a3+a5+)Q0+6(a2+a4+a6+)Q1+10(a3+a5+a7+)Q2+14(a4+a6+a8+)Q3+18(a5+a7+a9+)Q4(29)
By the procedure described in Appendix 2 and Eq. (25), the coefficients bk are Display Formula
b0=1=2(a1+a3+a5+),b1=6(12a0)=6(a2+a4+a6+),b2=10(12a1)=10(a3+a5+a7+),b3=14(12a0a2)=14(a4+a6+a8+),b4=18(12a1a3)=18(a5+a7+a9+),b5=22(12a0a2a4)=22(a6+a8+a10+).(30)

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

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

By Eq. (17), we obtain Display Formula

t3225Q0(t)+1835Q1(t)+221Q2(t)2165Q3(t)+185005Q4(t)21365Q5(t)=42145+28143t+224143t2224143t3+168143t42465t5.(31)

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

b0=1,b1=6×[12a0]=35,b2=10×[12a1]=17,b3=14×[12(a0+a2)]=115,b4=18×[12(a1+a3)]=377.(32)

Thus, the expansion of the derivatives is Display Formula

ddt(t3/2)Q0(t)+35Q1(t)17Q2(t)+115Q3(t)377Q4(t).(33)

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

Graphic Jump Location
Fig. 1
F1 :

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   k=0n(1)nan for the shape and k=0n(1)nbn for the slope. Both of them converge to zero when n.

Table Grahic Jump Location
Table 1Errors in shape and slope for the third order surface (N15).

Since f(t)=t3/2 and f(t)=(3/2)t1/2, the decreasing speed of the coefficients is as fast as n4   in shapes by Eq. (17) and n2 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.

Graphic Jump Location
Fig. 2
F2 :

Convergence of the maximum errors for f(t)=t3/2.

Approximation for the Cone f(t)=t1/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)=t1/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 Display Formula

bn=1/2(2n+1)Γ(1/2)2Γ(3/2+n)Γ(1/2n).(34)
By Eq. (15), the denominator is Display Formula
Γ(32+n)  Γ(12n)=(2n+2)!4n+1(n+1)!π×(4)nn!(2n)!π=(1)n(2n+2)(2n+1)π4(n+1)=(1)n2n+12π.(35)
Since Γ(1/2)=π, Display Formula
bn=(1)n.(36)
Thus, Display Formula
k=0bnQn(0)=k=0(1)n(1)n=k=01=.(37)
Hence, the sum |k=0bnQn(t)| does not converge at the origin.

Another explanation is possible as follows. Since the cone is described as   f(t)=t1/2=|r| on the interval [1,1], the differential is Display Formula

dfdr={1  (r<0)1  (r>0).(38)
This conclusion intuitively meets the fact that the cone has a singular point at the vertex.

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,22,23 the design with odd-order terms has been known as the Schmidt surface24.

Graphic Jump Location
Fig. 3
F3 :

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: Display Formula

z=A2r2+A3r3+A4r4+A6r6+(39)
Note that all optical design and analysis is achieved with Code-VTM.

Table Grahic Jump Location
Table 2Design specification.
Table Grahic Jump Location
Table 3Basic lens data.
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 Grahic Jump Location
Table 4Odd order and approximation coefficients.

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.

Graphic Jump Location
Fig. 4
F4 :

Fitting error in shape.

The peak-to-valley error in Fig. 4 is 3.88×105  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.

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 A2 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.

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.

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=r2 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 F/x and F/y. If F(x,y) is rotationally symmetric, F(x,y) should be a function of t=r2, i.e., F(x,y)=f(t).

Consider the Zernike expansion of f(t)=anQn(t), the slopes are rewritten as Display Formula

Fx=txdfdt=2x(anQn(t)),Fy=tydfdt=2y(anQn(t)).(40)

Thus, in the context of Zernike expansion and its derivatives, it is enough to prove the equality of the termwise derivative (anQn(t))=anQn(t) and its convergence.

Derivation of the Coefficients bk

We present a straightforward derivation of the expansion coefficients bk. By differentiating Eq. (4) formally, we obtain Display Formula

(a0Q0(t)+a1Q1(t)+a2Q2(t)+)=2n=0[{2×2n+1}Q2n(t)k=n+1a2k1+{2(2n+1)+1}Q2n+1(t)k=n+1a2k].(41)
At this stage, the convergence of Eq. (41) is not guaranteed.

The coefficients bm in Eq. (27) are calculated by the orthogonality of Qm(t), Display Formula

bm=(2m+1)01f(t)Qm(t)dt=(2m+1)[f(t)Qm(t)|0101f(t)Qm(t)dt]=(2m+1)[101f(t)Qm(t)dt].(42)

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

b2n=(4n+1)[101f(t)Q2n(t)dt]=(4n+1)[1201f(t)k=1n(4k1)Q2k1(t)dt].(43)

Since Display Formula

01f(t)(4k1)Q2k1(t)dt={2(2k1)+1}01f(t)Q2k1(t)dt=a2k1,

the expansion coefficient b2n is rewritten as Display Formula

b2n=(4n+1)[12k=1na2k1].(44)

By Eq. (25), Display Formula

b2n=2(4n+1)k=n+1a2k1.(45)

This is identical to the coefficient of Q2n(t) in Eq. (41). By the same procedure, the case of m=odd is proven, namely, Display Formula

b2n+1=2(4n+3)k=n+1a2k.(46)

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.

Neil  I. A., “Optical design dependence on technology development,” Opt. Eng.. 50, (12 ) 121706  (2011).CrossRef
Forbes  G. W., and Brophy  C. P., “Asphere, O Asphere, how shall we describe thee?,” Proc. SPIE. 7100, , 710002  (2008). 0277-786X CrossRef
Thompson  K. P.  et al., “The Forbes polynomial surface: the road forward; a return to manufacturable aspheres,” SPIE Optifab Panel Discussion. , 11  May  2011, http://spie.org/Documents/ConferencesExhibitions/2011_Optifab_panel_final.pdf (16  January  2016).
Yabe  A., “Sensitivity control to surface irregularity,” Proc. SPIE. 6342, , 634225  (2006).CrossRef
Shu  G.  et al., “Trapped ion imaging with a high numerical aperture spherical mirror,” J. Phys. B: At. Mol. Opt. Phys.. 42, , 154005  (2009). 0953-4075 CrossRef
Plummer  W. T., “Viewfinder for a reflex camera,” U. S. Patent No. 3,810,221 (1974).
Carsten  G., “Aberration-corrected microscope,” WO2012025349 A1 (2012).
Sato  K., “Wide-angle zoom lens including at least one aspheric lens surface,” U.S. Patent No. 6,982,834 (2006).
Amano  M., “Projection optical system and projection display apparatus,” U.S. Patent No. 8,967,812 (2015).
McGuire  J. P., “EUV condenser with non-imaging optics,” U.S. Patent No. 7,046,338 (2006).
Shibuya  M.  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).CrossRef
Tanikawa  M.  et al., “Effectiveness of odd order aspherical surface,” Jpn. J. Opt.. 36, (11 ) 646 –660 (2007) (in Japanese). 0389-6625 
Born  M., and Wolf  E., Principles of Optics. , 2nd (Revised) ed.,  Cambridge University Press ,  Cambridge, United Kingdom  (1999).
Zernike  F., “Beugungstheorie des schneidenver-fahrens und seiner verbesserten form, der phasenkontrastmethode,” Physica. 1, (7–12 ), 689 –704 (1934).CrossRef
Arfken  G., and Weber  H., Mathematical Methods for Physicists. , 5th ed.,  Academic Press ,  New York  (2001).
Bhatia  A. B., and Wolf  E., “On the circle polynomials of Zernike and related orthogonal sets,” Math. Proc. Cambridge Philos. Soc.. 50, , 40 –48 (1954). 0008-1981 CrossRef
Szegö  G., Orthogonal Polynomials. , 4th ed.,  American Mathematical Society ,  Providence, Rhode Island  (1975).
Braat  J. J. M., and Janssen  A. J. E. M., “Double Zernike expansion of the optical aberration function from its power series expansion,” J. Opt. Soc. Am. A. 30, (6 ), 1213 –1222 (2013). 0740-3232 CrossRef
Nijboer  B. R. A., “The diffraction theory of aberrations,” PhD Thesis, Rijksuniversiteit Groningen (1942).
Janssen  A. J. E. M., “Zernike expansion of derivatives and Laplacians of the Zernike circle polynomials,” J. Opt. Soc. Am. A. 31, (7 ), 1604 –1613 (2014). 0740-3232 CrossRef
Wünsche  A., “Generalized Zernike or disk polynomials,” J. Comput. Appl. Math.. 174, , 135 –163 (2005).CrossRef
Laikin  M., Lens Design. , 2nd ed., Revised and Expanded ,  Dekker ,  New York  (1995).
Gross  H., , Blechinger  F., and Achtner  A., Handbook of Optical Systems. , Vol. 4,  Wiley-VCH ,  Weinheim  (2008).
ISO-10110-12:2007, “Optics and photonics—preparation of drawings for optical elements and systems—Part 12: aspheric surfaces,”  International Organization for Standardization ,  Geneva  (2007)

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.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.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

Takao Tanabe ; Masato Shibuya and Kazuhisa Maehara
"Convergence and differentiation of Zernike expansion: application for an analysis of odd-order surfaces", Opt. Eng. 55(3), 035101 (Mar 01, 2016). ; http://dx.doi.org/10.1117/1.OE.55.3.035101


Figures

Graphic Jump Location
Fig. 1
F1 :

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

Graphic Jump Location
Fig. 2
F2 :

Convergence of the maximum errors for f(t)=t3/2.

Graphic Jump Location
Fig. 3
F3 :

Layout of Schmidt camera.

Graphic Jump Location
Fig. 4
F4 :

Fitting error in shape.

Tables

Table Grahic Jump Location
Table 1Errors in shape and slope for the third order surface (N15).
Table Grahic Jump Location
Table 2Design specification.
Table Grahic Jump Location
Table 3Basic lens data.
Table Grahic Jump Location
Table 4Odd order and approximation coefficients.

References

Neil  I. A., “Optical design dependence on technology development,” Opt. Eng.. 50, (12 ) 121706  (2011).CrossRef
Forbes  G. W., and Brophy  C. P., “Asphere, O Asphere, how shall we describe thee?,” Proc. SPIE. 7100, , 710002  (2008). 0277-786X CrossRef
Thompson  K. P.  et al., “The Forbes polynomial surface: the road forward; a return to manufacturable aspheres,” SPIE Optifab Panel Discussion. , 11  May  2011, http://spie.org/Documents/ConferencesExhibitions/2011_Optifab_panel_final.pdf (16  January  2016).
Yabe  A., “Sensitivity control to surface irregularity,” Proc. SPIE. 6342, , 634225  (2006).CrossRef
Shu  G.  et al., “Trapped ion imaging with a high numerical aperture spherical mirror,” J. Phys. B: At. Mol. Opt. Phys.. 42, , 154005  (2009). 0953-4075 CrossRef
Plummer  W. T., “Viewfinder for a reflex camera,” U. S. Patent No. 3,810,221 (1974).
Carsten  G., “Aberration-corrected microscope,” WO2012025349 A1 (2012).
Sato  K., “Wide-angle zoom lens including at least one aspheric lens surface,” U.S. Patent No. 6,982,834 (2006).
Amano  M., “Projection optical system and projection display apparatus,” U.S. Patent No. 8,967,812 (2015).
McGuire  J. P., “EUV condenser with non-imaging optics,” U.S. Patent No. 7,046,338 (2006).
Shibuya  M.  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).CrossRef
Tanikawa  M.  et al., “Effectiveness of odd order aspherical surface,” Jpn. J. Opt.. 36, (11 ) 646 –660 (2007) (in Japanese). 0389-6625 
Born  M., and Wolf  E., Principles of Optics. , 2nd (Revised) ed.,  Cambridge University Press ,  Cambridge, United Kingdom  (1999).
Zernike  F., “Beugungstheorie des schneidenver-fahrens und seiner verbesserten form, der phasenkontrastmethode,” Physica. 1, (7–12 ), 689 –704 (1934).CrossRef
Arfken  G., and Weber  H., Mathematical Methods for Physicists. , 5th ed.,  Academic Press ,  New York  (2001).
Bhatia  A. B., and Wolf  E., “On the circle polynomials of Zernike and related orthogonal sets,” Math. Proc. Cambridge Philos. Soc.. 50, , 40 –48 (1954). 0008-1981 CrossRef
Szegö  G., Orthogonal Polynomials. , 4th ed.,  American Mathematical Society ,  Providence, Rhode Island  (1975).
Braat  J. J. M., and Janssen  A. J. E. M., “Double Zernike expansion of the optical aberration function from its power series expansion,” J. Opt. Soc. Am. A. 30, (6 ), 1213 –1222 (2013). 0740-3232 CrossRef
Nijboer  B. R. A., “The diffraction theory of aberrations,” PhD Thesis, Rijksuniversiteit Groningen (1942).
Janssen  A. J. E. M., “Zernike expansion of derivatives and Laplacians of the Zernike circle polynomials,” J. Opt. Soc. Am. A. 31, (7 ), 1604 –1613 (2014). 0740-3232 CrossRef
Wünsche  A., “Generalized Zernike or disk polynomials,” J. Comput. Appl. Math.. 174, , 135 –163 (2005).CrossRef
Laikin  M., Lens Design. , 2nd ed., Revised and Expanded ,  Dekker ,  New York  (1995).
Gross  H., , Blechinger  F., and Achtner  A., Handbook of Optical Systems. , Vol. 4,  Wiley-VCH ,  Weinheim  (2008).
ISO-10110-12:2007, “Optics and photonics—preparation of drawings for optical elements and systems—Part 12: aspheric surfaces,”  International Organization for Standardization ,  Geneva  (2007)

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

PubMed Articles
Advertisement


 

  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.