Phase Only Pupil Filter Design Using Zernike Polynomials
 Author: Liu Jiang, Miao Erlong, Sui Yongxin, Yang Huaijiang
 Publish: Journal of the Optical Society of Korea Volume 20, Issue1, p101~106, 25 Feb 2016

ABSTRACT
A pupil filter is a useful technique for modifying the light intensity distribution near the focus of an optical system to realize depth of field (DOF) extension and superresolution. In this paper, we proposed a new design of the phase only pupil filter by using Zernike polynomials. The effect of design parameters of the new filters on DOF extension and superresolution are discussed, such as defocus Strehl ratio (S.R.), superresolution factor (G) and relative first side lobe intensity (M). In comparison with the other two types of pupil filters, the proposed filter presents its advantages on controlling both the axial and radial light intensity distribution. Finally, defocused imaging simulations are carried out to further demonstrate the effectiveness and superiority of the proposed pupil filter on DOF extension and superresolution in an optical imaging system.

KEYWORD
Wavefront coding , Noninteger phase mask , Depth of field , Imaging system

I. INTRODUCTION
In numerous applications it is necessary to extend the depth of field (DOF) and improve the resolution for an optical system, such as microlithography [1], optical data storage [2], and microscopy imaging [3, 4]. One of the most popular techniques to achieve DOF extension and superresolution includes introducing a pupil filter at the pupil plane of an optical system to modify the radial and axial spot size of the point spread function (PSF).
Since Toraldo made a seminal study on a radialsymmetric filter [5], a lot of work has been done on the design of a pupil filter. Although an amplitude only filter can provide effective superresolution and DOF extension, it still has two main difficulties. The first one is the reduction of image brightness and the other refers to the fabrication, especially for continuously varying amplitude functions [6, 7]. Therefore, in recent years phase only filter design has achieved a rapid development based on the parabolic approximation of the PSF for the general filters introduced by Sheppard [8, 9]. Typically, the design strategies of pupil filters include circular obstructions [1012], continuously varying functions, such as Bessel functions [13, 14] and Gaussian functions [15], and special forms, such as a spoke wheel filter [16].
In this paper, a phase only filter is proposed by using Zernike polynomials. Radial and axial spot sizes of an optical system with and without the designed filter are analyzed relatively based on scalar diffraction theory. Then, compared with the other two typical phase filters (threezones and Gaussian functions), we show its advantages on DOF extension and superresolution. Finally, we further demonstrate the availability and superiority of our proposed phase only filter through defocus imaging simulation.
II. BASIC PRINCIPLE
According to the Born and Wolf’s theory [17], at the focal region, the amplitude distribution of a converging monochromatic spherical wavefront passing through the pupil plane can be written as:
Where
ρ is the normalized radial coordinate of the pupil plane, represents the pupil function,P represents Bessel function of the first kind of order zero,J_{0} andu are the radial and axial dimensionless optical coordinates, given by:v Where
is the numerical aperture of the optical system,NA is the wavelength of the incident light,λ andr are the usual radial and axial distances.z Since most practical imaging systems use incoherent illumination, the light intensity distribution near the focal region can be calculated by its amplitude distribution, written as:
Where * denotes the complex conjugate.
After introducing the phase pupil filter into the optical system, the light intensity distribution will be modulated by the phase changing of the pupil function to realize the DOF extension and superresolution. These can be evaluated by the following three parameters:
1. Defocus Strehl ratio S.R.(0, zdefocus), defined as the ratio of the intensity at the focal point to that corresponding to a clear pupil, is a relevant parameter for analyzing the imaging quality. It is used to evaluate the ability of DOF extension for a pupil filter. Typically, the defocused imaging quality can be considered acceptable when S.R.(0, zdefocus) ≥ 0.8. Thus, the effective depth of focus can be written as DS.R. ≥ 0.8. 2. Superresolution factor G, defined as the ratio of the radius of the central lobe of the intensity distribution to that corresponding to a clear pupil, is a parameter for calculating the spot size reduction. G < 1 means the resolution enhancement for the optical system. 3. Relative first side lobe intensity M, defined as the ratio of between the central lobe intensity and the maximum intensity at the first side lobe, is also a parameter used to evaluate the ability of superresolution for a pupil filter. Typically, M must be limited to less than 0.2 in order not to reduce the resolution and cause energy loss.
These parameters can be written as:
Where
is the defocus distance,z_{defocus} (0,I_{max} ) is the intensity at the focal point with defocus distancez_{defocus} ,z_{defocus} is the intensity at the focal point of without defocus,I_{max}(0, 0) is the radius of central lobe of intensity,R_{central lobe} is the radial optical coordinate of maximum intensity at the first side lobe.v_{firstside lobe} III. PUPIL FILTER DESIGN
The pupil functions of an optical system with phase only pupil filter can be written as:
Where
is phase function of the pupil filter.Ψ(ρ) Compared with conventional power functions or Gaussian functions, there are some advantages for pupil filter design by using Zernike functions since they are closely related with optical aberrations. The phase functions with Zernike functions now can be written as:
Where
is the diffraction order, here set to 1,K is the number of Zernike coefficients,N is theZ_{i} Zernike standard polynomial andi^{th} is the corresponding coefficient, ρ is the normalized radial coordinate of the pupil plane andA_{i} is the angular coordinate.φ The first 16 polynomials of the Zernike functions are shown in Table 1 [18].
Here
polynomials of the Zernike functions are chosen to build the phase functions of the pupil filter by considering both availability for spot size modulation and the ease of fabrication. According to the relationships between Zernike functions and Seidel aberrations, roles that these polynomials play in the phase function can be analyzed reasonably.1^{st}, 4^{th}, 5^{th}, 9^{th}, 12^{th}, 16^{th} 1. 1st polynomial is a constant, representing the thickness of the pupil filter. Usually, the ratio between thickness and aperture of an optical element is not less than 0.2. Otherwise, the element will be fragile and hard to assemble. 2. 4th, 9th, 16th polynomials represent the defocus, primary spherical aberration and highorder spherical aberration respectively, which are useful to modulate the radius of the central lobe of the intensity distribution. In addition, these three have the ability to reduce the impact of defocus. 3. 5th, 6th, 12th, 13th polynomials represent the primary astigmatic aberration and highorder astigmatic aberration, which can limit the intensity distribution of the side lobe effectively. The 6th term is rotated 45 degrees from the 5th term and the 13th term is rotated 45 degrees from the 12th term. But they make no contribution to the DOF extension and superresolution. In order to simplify the optimization, only 5th, 12th terms are selected to suppress the energy of the side lobe.
The optimization process is described as follows:
A series of coefficients of chosen Zernike polynomials is used to calculate the PSF of the optical system with the phase pupil filter, and the
M can be obtained from the calculated PSF.If the value of
M is larger than the threshold (here the threshold is 0.2), select the next series of coefficients to recalculate the PSF since the optical system is no longer suitable for imaging with a biggerM .If the value of
M is not larger than the threshold, continue calculating defocus Strehl ratioS.R . and superresolution gainG of the optical system. Then, choose proper weight coefficients andκ_{1} to calculate the comprehensive evaluation result of this series of coefficients of Zernike polynomials. The coefficientsκ_{2} andκ_{1} are used to balance superresolution and DOF extension. In this paper, we are more concerned about DOF extension, so theκ_{2} equals 0.7 andκ_{1} equals 0.3.κ_{2} Repeat the above steps, selecting the best comprehensive evaluation result. The corresponding coefficients are the optimization results.
It can be presented as follows by Ctype codes:
Iteration:
for{(
A_{i} ) = Coef(i ),i ++}end
(
A_{i} )_{opt} =Coef (Maxima (Value ))In the optimization,
is a matrix containing all the combinations of coefficients within a certain range,Coef andValue_sr are used to keep the defocus Strehl ratioValue_g S.R . and superresolution gainG with different coefficients underM constraint, andκ_{1} are the weight coefficients.κ_{2} Thus, optimized coefficients of the proposed Zernike phase only pupil filter can be obtained by a global optimization algorithm, as shown below:
A _{1}= 0.250,A _{4}= −0.233,A _{5}= −0.126,A _{9}= 0.317,A _{12}= −0.071,A _{16}= 0.042.IV. DISCUSSION AND COMPARISON
By introducing the designed pupil filter into a conventional optical system with clear pupil, the improvement for imaging quality can be analyzed numerically through calculating the intensity distribution near the focal plane.
The radial intensity distributions at the focal point are shown in Fig. 1. The coordinates of the figure are in wavelengths. As Fig. 1 shows, the size of the central lobe is compressed by the designed pupil filter compared with the optical system of a clear pupil, which means there is resolution enhancement.
The ability of extending DOF of the Zernike phase only pupil filter can be reflected on the axial intensity distribution across the focal region, as shown in Fig. 2. Similarly, the coordinates are in wavelength. As can be seen from Fig. 2, the length of the bright spot is much longer than that of a conventional optical system, which characterizes the DOF extension.
Then, through comparison with the other two pupil filters, advantages of the Zernike phase only pupil filter are shown as follows. Here, forms of the two pupil filters are three zones [10] and Gaussian function [15], their phase functions can be respectively written as:
Three zones:
Where
andρ_{1} are the radiuses of the first and second rings, here,ρ_{2} = 0.64,ρ_{1} = 0.82 [10].ρ_{2} Gaussian functions:
Where
is the width of half height, andw is order of Gaussian functions, here,n = 0.12,w = 8 [15].n The comparison of radial intensity distributions at the focal point are shown in Fig. 3. It is clearly seen that the intensity distribution with the Zernike phase only pupil filter has the smallest radius of central lobe, which indicates its advantages on superresolution for a conventional optical system. Besides, the maximum intensity at the first side lobe of the proposed pupil filter is constrained quite low, which benefits from the
polynomials of the Zernike functions.5^{th}, 12^{th} The comparison of axial intensity distributions across the focal plane are shown in Fig. 4. The value of the transverse dotted line equals 0.8, which is used to indicate effective defocus Strehl ratio. That means the imaging quality is acceptable at a certain defocus distance if the corresponding normalized intensity is greater than this value. It can be seen in Fig. 4, the optical system with Zernike phase only pupil filter can keep its normalized intensity greater than 0.8 over the longest range of defocus distance.
The summary of the ability of DOF extension and superresolution for these three pupil filters is shown in the table below. In Table 1, effective depth of focal
≥ 0.8 of the clear pupil is normalized and superresolution factorD_{S.R.} is inverted to show the resolution enhancement more intuitively. Compared with conventional optical system, DOF is enlarged 4.15 times and resolution is enhanced 1.32 times when introducing the Zernike phase only pupil filter.G In order to demonstrate the ability of DOF extension and superresolution of our designed pupil filter further, imaging simulations are carried out. The virtual object used in the simulations is a resolution testing panel of line pairs with different spatial frequencies. Meanwhile, the defocus distance added in the simulations is from 0λ to 3λ. Results are shown in Table 2. From the first line of these figures in Table 2, the resolution of the optical system is clearly enhanced with all the three pupil filters compared to that with the clear pupil. In performance on DOF extension, the Zernike phase only pupil filter show its advantages. Its imaging results keep quite good quality over the entire defocus distance, which is very useful in numerous applications.
V. CONCLUSION
In this paper, we have proposed an optical phase only pupil filter using Zernike functions, used to realize DOF extension and superresolution. Based on analyzing the relationship between Seidel aberrations and Zernike polynomials, design parameters of the proposed pupil filter are optimized. Compared with other two typical pupil filters, the Zernike pupil filter shows its advantages on both DOF extension and superresolution, enlarging DOF 4.15 times and enhancing resolution 1.32 times. Finally, imaging simulations with different defocus distances further demonstrate its effectiveness and superiority on improving the performance of the conventional optical system.

[]

[]

[]

[]

[]

[]

[TABLE1.] Comparison of evaluation parameters for different phase pupil filters

[]

[]

[FIG. 1.] Radial intensity distributions of optical system with (b) and without (a) pupil filter.

[FIG. 2.] Axial intensity distributions of optical system with (b) and without (a) pupil filter.

[]

[]

[FIG. 3.] Comparison of normalized radial intensity distribution at focus plane with different phase pupil filters and clear pupil.

[FIG. 4.] Comparison of normalized axial intensity distribution across focus point with different phase pupil filters and clear pupil.

[TABLE 2.] Comparison of defocused imaging for different phase pupil filters