# Spherical Aberration Correction Using Refractive-Diffractive Lenses with an Analytic-Numerical Method

Authors: Sergio Vázquez-Montiel; Omar García-Liévanos; Juan Alberto Hernández-Cruz
Journal: Advances in Optical Technologies
 (2010)

---

## Abstract

We propose an alternative method to design diffractive lenses free of spherical aberration for monochromatic light. Our method allows us to design diffractive lenses with the diffraction structure recorded on the last surface; this surface can be flat or curved with rotation symmetry. The equations that we propose calculate the diffraction profiles for any substratum, for any f-number, and for any position of the object. We use the lens phase coefficients to compensate the spherical aberration. To calculate these coefficients, we use an analytic-numerical method. The calculations are exact, and the optimization process is not required.


---

## Body

## 1. Introduction

Spherical aberration is, in many cases, the most important of all primary aberrations, because it affects the whole field of the lens, including the vicinity of the optical axis. It is due to different focus positions for a marginal ray, meridional ray, and paraxial rays. An alternative to minimize the spherical aberration is to use diffractive optical elements (DOE). Diffractive lenses are essentially gratings with a variable spacing groove which introduces a chromatic aberration that is worse than conventional refractive/reflective optical elements. In some applications, an optical component may require a diffractive surface combined with a classic lens element. By using the diffractive properties, it is possible to design hybrid elements to obtain an achromatically corrected element [1]. In other cases, the requirements can be satisfied by just using a diffractive element.  In general, iterative methods are used to design these lenses [2]. Also, some people have used analytical third-order and numerical integrator methods to design diffractive lenses [1, 3]. The diffractive lenses we describe in this paper are limited to monochromatic applications; however, our proposed method is valid for a wide range of wavelengths.We use lens phase coefficients to compensate spherical aberration. To calculate these coefficients, we use an analytic-numerical method. The calculations are exact, simple and quick. A process of optimization is not required.The manufacturing problem of diffractive lenses is not considered here; to solve this problem you can read Castro-Ramos et al. [4].First, we describe the diffractive lenses theory. Also, we give a brief derivation of the general grating equation to trace a couple of light rays through a rotationally symmetrical surface. Then, we establish the analytic-numerical method to minimize spherical aberration.  We propose some heights to correct the spherical aberration. Finally, we conclude by providing a design example.


## 2. Theory of the Diffractive Lenses

Diffractive lenses can be described by a polynomial phase function [5](1)ϕ(x,y)=2πλ∑m∑namnxmyn,
where λ is the design wavelength; amn are the lens phase coefficients; x, y are the coordinates in the diffractive lens.We will consider that the diffractive lens is rotationally symmetrical, so (1) is rerewritten as(2)ϕ(y)=2πλ(a0+a2y2+a4y4+a6y6+a8y8+⋯).
Here, the longitudinal displacement of the reference sphere is a0=0 because we have assumed it is in the ideal focus. The coefficient a2 are implicit lens paraxial properties; it is equal to -1/2f, where f is the focal length. The remainder coefficients in (2) give the amount of spherical aberrations of the first- second- and higher-order [6, 7].Designers usually use some commercial optical design programs to obtain the lens phase coefficients by using an optimization process. We will describe an analytical method to obtain these coefficients.To trace a pencil of rays through the diffractive optical surface, we use the grating equation. For a flat surface, the grating equation is given by(3)n′sinI′-nsinI=mλf,
where n′ and n are refractive indexes for two different mediums, I′ and I are the diffractive and incident angles; f is the grating frequency; m is the diffracted order.To analyze the light propagation through a diffractive curved surface, we have to change the form of the last equation. After some algebra, we obtain the general grating equation(4)(n′M2-nM1)cosθN+(nN1-n′N2)sinθN=mλcosθNf,
where the direction of refracted and diffracted rays is given by the direction cosines M1,M2,N1,  and N2 as are shown in the Figure 1. This analysis considers that the diffractive lens is rotationally symmetrical, and then the direction cosines L0=0,L1=0,L2=0. θN is the angle between the normal at surface and optical axis, it is given by (5)sinθN,cosθN=(∂F/∂y,∂F/∂z){(∂F/∂x)2+(∂F/∂y)2+(∂F/∂z)2}1/2,
here, F is the surface function in which the diffractive lens will be recorded, and x and y are the surface coordinates.Figure 1
Lens parameter. The diffractive surface is on the second surface.The grating frequency of (3) and (4) can be calculated in one dimension by(6)fy=12π∂ϕ∂y,fy=1λ[∑k=1∞2ka2ky2k-1],
where ϕ is the phase function given by (2) k is an integer equal to 1, 2, 3, 4, …, and the diffracted order m=1. With this, (4) can be rewritten as(7)[(n′M2-nM1)cosθN+(nN1-n′N2)sinθN]=cosθN[∑k=1∞2ky2k-1a2k].
Using (5) and (7), we can trace n rays through the surface at different heights on the pupil. Then they can arrange a k×k equations system. The number of equations depends of the number of coefficients that we want to find.Then, if we want to findk coefficients, we need to solve an equation system similar to(8)(8)[w11w12w13⋯w1kw21w22w23⋯w2kw31w32w33⋯w3k⋮⋮⋮⋯⋮wk1wk2wk3⋯wkk][-12fa4a6⋮a2k]=[A0A1A2⋮Ak],
where w represents the different constants of the right side of (7), and A are the constants of the left side of the same equation, all for different height rays on the pupil.


## 3. How Many and What Heights Should Be Corrected

The spherical aberration of the ray in any optical system can be expressed as(9)∂W(0,y)∂y=4b1(y3)+6b2(y5)+8b3(y7)+⋯.
Considering only big f/numbers, the spherical aberration can be represented only for the first and second terms, and combining these terms, the spherical aberration of the edge can be corrected. Then, the peak of the spherical aberration residual occurs when y is equal to the marginal ym multiplied by 3/5=0.7746. This analysis is similar to Kingslake [8]. The difference is that the defocus term is not considered here.It is possible to correct the residual spherical aberration by using the third term of the expansion (9), but now the ray aberration curve has two opposite peaks above and below the 0.7746 zone. The zones with maximum and minimum residuals fall at values of y given by y/ym=0.5550 or 0.9137 (see Figure 2).Figure 2
Transversal spherical aberration of the example one.If we consider f/numbers to be small, we should correct the spherical aberration residual, and its peaks fall at valuesy=ym(0.5550) or ym(0.9137); then we need fourth and fifth term to correct these other y’s, now the ray aberration curve has two opposite peaks above and two below of 0.5550 and 0.9137 zones. The zones with maximum and minimum residuals fall at values of y given by y/ym=0.9681, 0.8505, 0.6661, or 0.3740, Figure 7. This analysis can continue because the expansion (9) is infinity.The points for Kingslake analysis arey/ym=1, 0.8880, 0.7071, and 0.4597, and for our analysis y/ym=1, 0.9137, 0.7746, and 0.5550.The number ofy that must be corrected for each optical system depends on the optical system tolerances, for example, with one value of y ((1)(ym)), we correct a lens with f/number bigger than f/5; with two different values of y((1)(ym) and (0.7746)(ym)), we correct a lens with f/number bigger than f/2; with four different values of y((1)(ym), (0.7746)(ym), (0.5550)(ym), and (0.9137)(ym)), we correct a lens with f/number bigger than f/1, but only the designer should decide the correction that he needs.


## 4. Results

We have proposed a general expression to compute the phase coefficients. Now, we will show how theses coefficients minimize the spherical aberration with some numerical examples. All examples considered in this section have the diffracted orderm=1.


### 4.1. Example  1

In this example, we consider that the diffractive surface is on a spherical surface (the last surface of the system) with 50 mm of diameter aperture, numerical aperture 0.375, object distance 200 mm, andλ=0.587μm.In Table1, other characteristics of the refractive-diffractive lens are shown.Table 1
Example1 data.
SurfacesRadius (mm)Thickness (mm)Radius aperture (mm)Glass1∞200Air2101.9548.13825BK73DOE−101.95466.05925Air4∞0AirWe must trace light rays until the last surface, and then we can calculate all constants of (8).The number of rays traced depends on the number of coefficients. In this example, we use two coefficients, and we get the next equations system(10)[6.340530E+46.130961E+72.924391E+41.6731E+7][a4a6]=[0.0670530.032394].We have solved (8) to compute the phase coefficients for two different pupil positions on the surface; they are shown in Table 2.Table 2
Coefficients value for the example1.
CoefficientsAperture height (mm)Valuea2Paraxial−0.005 mm-1a4(25)(1)=251.180409 × 10-6 mm-3a6(25)(0.7746)=19.363−1.270732 × 10-10 mm-5Figure2 shows the spherical aberration of the refractive-diffractive lens; the graphics were obtained using the commercial optical design program “OSLO” [9].We can see in the graphic a maximum transversal spherical aberration of about 0.004 mm, having zeros on two pupil positions. This is because we had computed two coefficients for the system. The corresponding Strehl Ratio is of about 0.151.In the Figure2, FBY and FBX are the fractional object coordinates, and WV1 is the wavelength (λ=0.587μm) for the evaluation.


### 4.2. Example  2

We consider the same optical system but now using four phase coefficients. Solving the next equations system(11)[6.34053E+46.130961E+75.269617E+104.246204E+134.821205E+43.869253E+72.760235E+101.84602E+132.924391E+41.6731E+78.50856E+94.056593E+121.071387E+43.119059E+68.071391E+81.958143E+11][a4a6a8a10]=[0.0670530.0519530.0323940.012299], we obtain phase coefficients which are shown in Table 3.Table 3
Coefficients value for the example2.
CoefficientsAperture height (mm)Valuea2Paraxial−0.005 mm-1a4(25)(1)=251.194430 × 10-6 mm-3a6(25)(0.9137)=22.842−1.684803 × 10-10 mm-5a8(25)(0.7746)=19.3633.602457 × 10-14 mm-7a10(25)(0.555)=13.875−5.856954 × 10-14 mm-9In Figure3, we can see a maximum traversal spherical aberration of the refractive-diffractive lens of about 0.00005 mm, having zeros on four pupil positions. The reason is that we had computed four coefficients for this system. The corresponding Strehl Ratio is of about 1. Figure 3 also shows the difference between the points proposed by Kingslake [6] and our selected points. It can be seen that the points suggested in this paper to correct the spherical aberration are slightly better than the Kingslake points.Transversal spherical aberration of the example2.
(a)(b)


### 4.3. Example  3

Now we consider the same optical system but the diffractive surface on a hyperbolic surface (last surface) with conic constant, diameter apertureK=-4.654, 50 mm, numerical aperture of 0.375, object distance of 200 mm, and λ=0.587μm. We must trace rays to the hyperbolic surface because in this way we can calculate all constants of (8) for this example. We use four phase coefficients to solve the following equations system:(12)[6.422962E+46.209924E+75.336845E+104.299859E+134.86619E+43.905076E+72.785591E+101.862844E+132.939127E+41.681484E+78.550953E+94.076689E+121.072902E+43.123454E+68.082727E+81.960884E+11][a4a6a8a10]=[0.0487920.0380010.0239069.203867E-3]In the Table 4 are the new coefficients for this optical system.Table 4
Coefficients value for the example 3.
CoefficientsAperture height (mm)Valuea2Paraxial−0.005 mm-1a4(25)(1)=259.106091 × 10-7 mm-3a6(25)(0.9137)=22.842−1.941142 × 10-10 mm-5a8(25)(0.7746)=19.3635.218183 × 10-14 mm-7a10(25)(0.555)=13.875−9.916744 × 10-18 mm-9Figure4 shows the aberration of this refractive-diffractive lens.Figure 4
Transversal spherical aberration of the example 3.We can see again a very small spherical aberration, and its maximum value is of around 2 × 10-5 mm. It has 4 zeros because we have used 4 phase coefficients. The irradiance distribution corresponding to this system is shown in Figure 5.Figure 5
The point spread function of the example 3.Our proposed method is also for flat surfaces. We only use zero for the angle between the normal to the surface and optical axis in (4), and then we obtain the grating (3) for a flat surface. Then, we can use the procedure that we used in the previous examples.If the designer wants to use the first surface, the conjugates must be changed, and then the method proposed can be applied.


## 4.1. Example  1

In this example, we consider that the diffractive surface is on a spherical surface (the last surface of the system) with 50 mm of diameter aperture, numerical aperture 0.375, object distance 200 mm, andλ=0.587μm.In Table1, other characteristics of the refractive-diffractive lens are shown.Table 1
Example1 data.
SurfacesRadius (mm)Thickness (mm)Radius aperture (mm)Glass1∞200Air2101.9548.13825BK73DOE−101.95466.05925Air4∞0AirWe must trace light rays until the last surface, and then we can calculate all constants of (8).The number of rays traced depends on the number of coefficients. In this example, we use two coefficients, and we get the next equations system(10)[6.340530E+46.130961E+72.924391E+41.6731E+7][a4a6]=[0.0670530.032394].We have solved (8) to compute the phase coefficients for two different pupil positions on the surface; they are shown in Table 2.Table 2
Coefficients value for the example1.
CoefficientsAperture height (mm)Valuea2Paraxial−0.005 mm-1a4(25)(1)=251.180409 × 10-6 mm-3a6(25)(0.7746)=19.363−1.270732 × 10-10 mm-5Figure2 shows the spherical aberration of the refractive-diffractive lens; the graphics were obtained using the commercial optical design program “OSLO” [9].We can see in the graphic a maximum transversal spherical aberration of about 0.004 mm, having zeros on two pupil positions. This is because we had computed two coefficients for the system. The corresponding Strehl Ratio is of about 0.151.In the Figure2, FBY and FBX are the fractional object coordinates, and WV1 is the wavelength (λ=0.587μm) for the evaluation.


## 4.2. Example  2

We consider the same optical system but now using four phase coefficients. Solving the next equations system(11)[6.34053E+46.130961E+75.269617E+104.246204E+134.821205E+43.869253E+72.760235E+101.84602E+132.924391E+41.6731E+78.50856E+94.056593E+121.071387E+43.119059E+68.071391E+81.958143E+11][a4a6a8a10]=[0.0670530.0519530.0323940.012299], we obtain phase coefficients which are shown in Table 3.Table 3
Coefficients value for the example2.
CoefficientsAperture height (mm)Valuea2Paraxial−0.005 mm-1a4(25)(1)=251.194430 × 10-6 mm-3a6(25)(0.9137)=22.842−1.684803 × 10-10 mm-5a8(25)(0.7746)=19.3633.602457 × 10-14 mm-7a10(25)(0.555)=13.875−5.856954 × 10-14 mm-9In Figure3, we can see a maximum traversal spherical aberration of the refractive-diffractive lens of about 0.00005 mm, having zeros on four pupil positions. The reason is that we had computed four coefficients for this system. The corresponding Strehl Ratio is of about 1. Figure 3 also shows the difference between the points proposed by Kingslake [6] and our selected points. It can be seen that the points suggested in this paper to correct the spherical aberration are slightly better than the Kingslake points.Transversal spherical aberration of the example2.
(a)(b)


## 4.3. Example  3

Now we consider the same optical system but the diffractive surface on a hyperbolic surface (last surface) with conic constant, diameter apertureK=-4.654, 50 mm, numerical aperture of 0.375, object distance of 200 mm, and λ=0.587μm. We must trace rays to the hyperbolic surface because in this way we can calculate all constants of (8) for this example. We use four phase coefficients to solve the following equations system:(12)[6.422962E+46.209924E+75.336845E+104.299859E+134.86619E+43.905076E+72.785591E+101.862844E+132.939127E+41.681484E+78.550953E+94.076689E+121.072902E+43.123454E+68.082727E+81.960884E+11][a4a6a8a10]=[0.0487920.0380010.0239069.203867E-3]In the Table 4 are the new coefficients for this optical system.Table 4
Coefficients value for the example 3.
CoefficientsAperture height (mm)Valuea2Paraxial−0.005 mm-1a4(25)(1)=259.106091 × 10-7 mm-3a6(25)(0.9137)=22.842−1.941142 × 10-10 mm-5a8(25)(0.7746)=19.3635.218183 × 10-14 mm-7a10(25)(0.555)=13.875−9.916744 × 10-18 mm-9Figure4 shows the aberration of this refractive-diffractive lens.Figure 4
Transversal spherical aberration of the example 3.We can see again a very small spherical aberration, and its maximum value is of around 2 × 10-5 mm. It has 4 zeros because we have used 4 phase coefficients. The irradiance distribution corresponding to this system is shown in Figure 5.Figure 5
The point spread function of the example 3.Our proposed method is also for flat surfaces. We only use zero for the angle between the normal to the surface and optical axis in (4), and then we obtain the grating (3) for a flat surface. Then, we can use the procedure that we used in the previous examples.If the designer wants to use the first surface, the conjugates must be changed, and then the method proposed can be applied.


## 5. Conclusions

We have established a new exact method to correct the spherical aberration for any optical system using diffractive lenses; this method makes use of the general grating equation and exact ray trace. With our method, we can decide how many zeros the spherical aberration should have and fix its position in the exit pupil. The method can only be applied to the first and last surface of the optical system.We also have proposed some heights to correct the spherical aberration and how many rays must be traced depending on the f/number.In the first and second examples, we have shown that we can have a high control of spherical aberration, minimized at points on the surface where we have wanted. Also, we have shown that our method is valid for any rotationally symmetrical surface.In general, spherical aberration will have as many zeros as the coefficients we calculate. It is very important to see that in order to minimize spherical aberration, we use only as many coefficients as necessary.Finally, to calculate the coefficients, we only use the analytic-numerical method. The calculations are exact, simple, and quick. A process of optimization is not required.



---

*Source: 783206-2010-08-17.xml*