Tight focusing of a higher-order radially polarized beam transmitting through multi-zone binary phase pupil filters

When the pupil filters are used to improve the performance of the imaging system, the conversion efficiency is a critical characteristic for real applications. Here, in order to take full advantage of the subwavelength focusing property of the radially polarized higher-order Laguerre-Gaussian (LG) beam, we introduce the multi-zone binary phase pupil filters into the imaging system to deal with the problem that the focal spot is split along the z axis for the small size parameter of the incident LG beam. We provide an easy-to-perform procedure for the design of multi-zone binary phase pupil filters, where the zone numbers of π phase are uncertain when the optimizing procedure starts. Based on this optimizing procedure, we successfully find the set of optimum structures of a seventeen-belt binary phase pupil filters and generate the excellent focal spot, where the depth of focus, the focal spot transverse size, the Strehl ratio, and the sidelobe intensity are 9.53λ , 0.41λ , 41.75% and 16.35% in vacuum, respectively. Most importantly, even allowing the power loss of the incident LG beam truncated by the pupil of the imaging system, the conversion efficiency is still as high as 37.3%. Theoretical calculations show that we succeed to have sufficient conversion efficiency while utilizing the pupil filters to decrease the focal spot and extend the depth of focus. ©2013 Optical Society of America OCIS codes: (110.2990) Image formation theory; (050.1960) Diffraction theory; (260.5430) Polarization. References and links 1. B. Richards and E. Wolf, “Electromagnetic diffraction in optical systems II. Structure of the image field in an aplanatic system,” Proc. R. Soc. Lond. A Math. Phys. Sci. 253(1274), 358–379 (1959). 2. H. Guo, J. Chen, and S. Zhuang, “Resolution of aplanatic systems with various semiapertures, viewed from the two sides of the diffracting aperture,” J. Opt. Soc. Am. A 23(11), 2756–2763 (2006). 3. R. Dorn, S. Quabis, and G. Leuchs, “Sharper focus for a radially polarized light beam,” Phys. Rev. Lett. 91(23), 233901 (2003). 4. Y. J. Yoon, W. C. Kim, N. C. Park, K. S. Park, and Y. P. Park, “Feasibility study of the application of radially polarized illumination to solid immersion lens-based near-field optics,” Opt. Lett. 34(13), 1961–1963 (2009). 5. H. Guo, J. Chen, and S. Zhuang, “Vector plane wave spectrum of an arbitrary polarized electromagnetic wave,” Opt. Express 14(6), 2095–2100 (2006). 6. K. J. Moh, X. C. Yuan, J. Bu, S. W. Zhu, and B. Z. Gao, “Radial polarization induced surface plasmon virtual probe for two-photon fluorescence microscopy,” Opt. Lett. 34(7), 971–973 (2009). 7. Y. Ku, C. Kuang, X. Hao, Y. Xue, H. Li, and X. Liu, “Superenhanced three-dimensional confinement of light by compound metal-dielectric microspheres,” Opt. Express 20(15), 16981–16991 (2012). 8. Q. Zhan, “Trapping metallic Rayleigh particles with radial polarization,” Opt. Express 12(15), 3377–3382 (2004). #183308 $15.00 USD Received 10 Jan 2013; revised 18 Feb 2013; accepted 18 Feb 2013; published 25 Feb 2013 (C) 2013 OSA 11 March 2013 / Vol. 21, No. 5 / OPTICS EXPRESS 5363 9. N. Hayazawa, Y. Saito, and S. Kawata, “Detection and characterization of longitudinal field for tip-enhanced Raman spectroscopy,” Appl. Phys. Lett. 85(25), 6239–6241 (2004). 10. K. Yoshiki, K. Ryosuke, M. Hashimoto, T. Araki, and N. Hashimoto, “Second-harmonic-generation microscope using eight-segment polarization-mode converter to observe three-dimensional molecular orientation,” Opt. Lett. 32(12), 1680–1682 (2007). 11. M. Meier, V. Romano, and T. Feurer, “Material processing with pulsed radially and azimuthally polarized laser radiation,” Appl. Phys., A Mater. Sci. Process. 86(3), 329–334 (2007). 12. G. Terakado, K. Watanabe, and H. Kano, “Scanning confocal total internal reflection fluorescence microscopy by using radial polarization in the illumination system,” Appl. Opt. 48(6), 1114–1118 (2009). 13. L. Novotny, M. R. Beversluis, K. S. Youngworth, and T. G. Brown, “Longitudinal field modes probed by single molecules,” Phys. Rev. Lett. 86(23), 5251–5254 (2001). 14. K. S. Youngworth and T. G. Brown, “Focusing of high numerical aperture cylindrical-vector beams,” Opt. Express 7(2), 77–87 (2000). 15. H. Wang, L. Shi, B. Lukyanchuk, C. Sheppard, and C. T. Chong, “Creation of a needle of longitudinally polarized light in vacuum using binary optics,” Nat. Photonics 2(8), 501–505 (2008). 16. K. Kitamura, K. Sakai, and S. Noda, “Sub-wavelength focal spot with long depth of focus generated by radially polarized, narrow-width annular beam,” Opt. Express 18(5), 4518–4525 (2010). 17. C. Kuang, X. Hao, X. Liu, T. Wang, and Y. Ku, “Formation of sub-half-wavelength focal spot with ultra long depth of focus,” Opt. Commun. 284(7), 1766–1769 (2011). 18. H. Dehez, A. April, and M. Piché, “Needles of longitudinally polarized light: guidelines for minimum spot size and tunable axial extent,” Opt. Express 20(14), 14891–14905 (2012). 19. Y. Kozawa and S. Sato, “Sharper focal spot formed by higher-order radially polarized laser beams,” J. Opt. Soc. Am. A 24(6), 1793–1798 (2007). 20. Y. Kozawa, T. Hibi, A. Sato, H. Horanai, M. Kurihara, N. Hashimoto, H. Yokoyama, T. Nemoto, and S. Sato, “Lateral resolution enhancement of laser scanning microscopy by a higher-order radially polarized mode beam,” Opt. Express 19(17), 15947–15954 (2011). 21. Y. Kozawa and S. Sato, “Focusing of higher-order radially polarized Laguerre-Gaussian beam,” J. Opt. Soc. Am. A 29(11), 2439–2443 (2012). 22. B. Tian and J. Pu, “Tight focusing of a double-ring-shaped, azimuthally polarized beam,” Opt. Lett. 36(11), 2014–2016 (2011). 23. C. J. R. Sheppard and A. Choudhury, “Annular pupils, radial polarization, and superresolution,” Appl. Opt. 43(22), 4322–4327 (2004). 24. M. Yun, L. Liu, J. Sun, and D. Liu, “Three-dimensional superresolution by three-zone complex pupil filters,” J. Opt. Soc. Am. A 22(2), 272–277 (2005). 25. T. G. Jabbour and S. M. Kuebler, “Vector diffraction analysis of high numerical aperture focused beams modified by twoand three-zone annular multi-phase plates,” Opt. Express 14(3), 1033–1043 (2006).


Introduction
In an imaging system with high numerical aperture (NA), because of the asymmetry of polarization with respect to the imaging system, the focal spot is not a circle but ellipsoidal elongating in the polarization direction for a linearly polarized beam [1,2]. If a circularly polarized beam is used, the focal spot is circularly symmetric, but the focal spot size is larger than that for a linearly polarized beam [2]. Compared with a linearly polarized beam, a radially polarized beam can produce a smaller focal spot with circularly symmetric intensity distribution for a high NA imaging system [3], which means that it cannot only improve the resolutions, but also make the resolutions in all transverse directions same of the imaging system [4]. From the view of vector plane wave spectrum, as there are only TM-wave (the polarized direction of magnetic field perpendicular to the plane formed by wave vector and longitudinal direction) [5], the surface plasmon mode is excited more efficiently by the radially polarized beam [6,7]. Another attractive characteristic of the radially polarized beam is that it can create the focal spot with a strong longitudinal electric field. This strong longitudinal electric field can be effectively used for particle manipulation [8], Raman spectroscopy [9], second-harmonic generation [10], material processing [11], and fluorescence imaging [12]. However, most of the studies on the radially polarized beam and its applications have concerned the so-called doughnut-shaped beam, which is represented as a fundamental radially polarized mode ( 01 TM mode) [3,4,[6][7][8][9][10][11][12][13][14][15][16][17][18]. Recently, Kozawa et. al report that the radially polarized higher-order Laguerre-Gaussian beam (LG) has higher performance than the low-order one and can generate the focal spot beyond the Abbe's diffraction limit 0.5λ  [19][20][21]. Tian and Pu also utilize the double-ring-shaped azimuthally polarized LG beam to form a subwavelength focal hole near the focus [22]. This kind of higher-order LG beam has concentric multi ring-shaped intensity pattern with π phase shift between the rings. Therefore, the higher-order LG beam is equivalent to the zero-order LG beam transmitting through the combination of a multi-zone binary phase and a successive amplitude pupil filters. However, for the higher-order LG generated directly by the laser, as the concentric multi ring-shaped intensity pattern is the inherent properties of the higher-order LG beam, the equivalent successive amplitude pupil filters have no loss of light energy when we cannot consider how to generate the higher-order LG. In other words, for the higher-order LG generated directly by the laser, the equivalent successive amplitude pupil filters not only can reduce the size of the focal spot, but also have no loss of light energy unlike the amplitude pupil filters placed at the pupil of the imaging system illuminated by the zero-order LG beam. Unfortunately, the equivalent combination of a multi-zone binary phase and a successive amplitude pupil filters and the relationship between the two cannot be adjusted arbitrarily for the higher-order LG beam. The natural question arises whether we can further reduce the focal spot size or extend the depth of focus by utilizing the multi-zone binary phase pupil filters to redistribute the phase shift between the rings of the radially polarized higher-order LG beam.
The binary phase pupil filters are widely used to reduce the focal spot size or extend the depth of focus, but, because of the complexity of optimizing multi-zone binary phase pupil filters, most of authors only focus on the design for the two-or three-zone binary phase pupil filters [23][24][25]. Wang et al. have an optimized design on the five zone binary phase filters, but the optimization method is not a strictly defined procedure [15]. To our knowledge, there is not an easy-to-perform procedure for the design of multi-zone binary pupil filters until now. However, as the focusing properties of the radially polarized higher-order LG beam are already excellent, i.e., beyond the Abbe's diffraction limit 0.5λ  , and the intensity of each ring is successive and inhomogenous, we think that the multi-zone binary phase pupil filters are more likely to adjust finely the phase shift between the multi ring-shaped intensity patterns of the higher-order LG beam so that one might further improve its focusing properties.
In addition, when the pupil filters are used to improve the performance of the imaging system, the conversion efficiency is a critical characteristic for real applications. Wang et. al, characterize the conversion efficiency as the ratio of the total energy within the focal volume after using the pupil filters to its original energy within the focal volume before using the pupil filters [15]. How to have sufficient conversion efficiency while utilizing the pupil filters to decrease the focal spot and extend the depth of focus is still a big challenge until now. For example, the narrow-width annular beams are used to decrease the focal spot and extend the depth of focus, but all have insufficient conversion efficiency (at the level of a few per cent) [16][17][18].
In this paper, we introduce the multi-zone binary phase pupil filter into the imaging system while taking full advantage of the subwavelength focusing property of the radially polarized higher-order LG beam. We provide an easy-to-perform procedure for the design of multi-zone binary phase pupil filters, where the zone numbers of π phase are uncertain when the optimizing procedure starts. Based on this optimizing procedure, we successfully find the set of optimum structures of a seventeen-belt binary phase pupil filters and generate the excellent focal spot, where the depth of focus, the focal spot transverse size, the Strehl ratio, and the sidelobe intensity are 9.53λ , 0.41λ , 41.75% and 16.35% in vacuum, respectively. Most importantly, even allowing the power loss of the incident LG beam truncated by the pupil of the imaging system, the conversion efficiency is still as high as 37.3%. To the best of our knowledge, this is the first report of the focal spot with subwavelength transverse size ( 0.41λ ), ultra long depth of focus ( 9.53λ ), and high conversion efficiency (37.3%) until now.  Fig. 1, we assume that the time dependence is exp(

Shown in
and the incident light propagating along z + axis is focused by an aplanatic imaging system L obeying the sine condition. The ideal focus O is the origin of rectangular coordinates ( , , ) x y z and its corresponding cylindrical coordinates ( , , ) r z ϕ . On the basis of the classical vector diffraction theory of polarized beams [1], the electric fields in the focal region for the high-NA imaging system illuminated by a radially beam can be expressed as [14,19] 1 2 where the maximum aperture angle a sin(NA / ) n α = , n is the refractive index in image space, k is the wave number in image space, A is a constant, m J is the first kind of Bessel functions of orders m , and ( ) T θ is transmission function of the pupil filters. 0 ( ) l θ is the relative amplitude of the electric field of the incident beam at the pupil plane. For the higherorder LG incident beam, its relative amplitude can be expressed as [19] where the size parameter 0 β of the incident beam is defined as the ratio of the pupil radius to the beam waist, 1 p L is the generalized Laguerre polynomial, and p is the radial mode number.
In this paper, the multi-zone binary phase pupil filters is used (see Fig. 1). As the transmission function ( ) T θ of the multi-zone binary phase pupil filters is irrelevant to the azimuthal angle ϕ and only dependent on the angle θ with 0 θ α ≤ < , ( ) T θ is a onedimensional function of the variable θ . Therefore, in order to clearly introduce the optimizing procedure for the multi-zone binary phase pupil filters, we consider the multi circular belts shown in Fig. 1 as the multi rectangular belts with variation of angle θ (see Fig.   2), namely each rectangular belt denoting one circular belt. Shown in Fig. 2(a), the integral interval 0 θ α ≤ < of angle θ in Eqs. (1a) and (1b) is divided into N rectangular belts and each rectangular belt corresponds to the angle interval The widths of each rectangular belt, i.e., the angle interval 1 i i θ θ θ − ≤ < , may be same or different. We number all rectangular belts from left to right and utilize these numbers to denote the positions of the rectangular belts. Our optimizing procedure for the multi-zone binary phase pupil filters are divided into four steps. In the first step, we assume the binary phases 0 and π are used in this paper. We let no E and yes E denote the electric fields in the focal region for the imaging system without pupil filters and with circular belts of phase π , respectively. Namely no E is the results of the integrand of Eqs. (1a) and (1b)   In the second step, shown in Fig. 2(a), we first assume all N rectangular belts have not phase π and the numbers of rectangular belts that will be added phase π are M ( 1 M N ≤ ≤ ). Then we start to add phase π into M rectangular belts of the total N rectangular belts. This adding phase π procedure are always carried out in sequence, from left to right, and can be described by the following recursive function where M still denotes the numbers of rectangular belts that will be added phase π indicated before. m denotes that we are adding the th m phase π into the pupil filters. In the third step, when the recursive function (3) stops, we start to calculate the electric fields yes E in terms of the actual positions of the M rectangular belts with phase π recorded by ℜ . Then we examine whether whole E meets the requirements determined by the merit parameters. The merit parameters we have selected are the Strehl ratio S , the sidelobe intensity r I , and the superresolution gain factor G , respectively. S is a relevant parameter for estimating conversion efficiency [15] and is defined as the ratio of the intensity at the focal point of the imaging system with to that corresponding to the one without pupil filters. r I is defined as the second maximum intensity relative to the central peak intensity at the focal plane. G is defined as the ratio of the focal spot size of the imaging system with to that corresponding to the one without pupil filters, where the focal spot size is defined as the full width at half-maximum (FWHM) of the focal spot at the focal plane. As indicated by Wang et. al, the conversion efficiency is a critical characteristic for real applications [15]. So, during optimizing procedure, 30% S > , 20% r I < , and 1 G < must be met because 30% S > can basically assure a conversion efficiency above 20%. Through the above three steps, we can obtain the optimum structures of the multi-zone binary pupil filters when M rectangular belts are added phase π .
In the final step, we further search the final optimum structures of the multi-zone binary phase pupil filters by changing M . As the phases of the belts are either zero or π , we only consider the case of M varying from one to the minimum integer bigger than 2 N because the optimizing results are same for the other half. After the optimizing procedure finished, the final optimum structures (i.e., the zone numbers and their radii) of the multi-zone binary phase pupil filters can be calculated by ℜ . It is noted that M is not the zone numbers of the binary phase pupil filters because many rectangular belts with phase zero or π might be adjacent (see Fig. 2(b)).
The above optimizing procedure is not only a strictly defined procedure, but also an easyto-perform. Moreover, the zone numbers of binary phase pupil filters are uncertain when the optimizing procedure starts. The final optimization structures of the multi-zone binary phase pupil filters are majorly determined by the merit parameters.

Simulations and discussions
As the higher the order of LG beam is, the smaller the focal spot size is in the case of the same NA [19], we will only focus on the focusing properties of the radially polarized fiveorder LG beam, namely 5 p = in the generalized Laguerre polynomial 1 p L (hereafter called the 5,1 LG beam). In addition, the NA 0.95 = , the refractive index 1 n = in image space, and the constant 1 A = in Eqs. (1a) and (1b) are used in our calculations.

Effects of the size parameter 0 β of the incident beam on the distribution of the focal spot
In this Subsection, the multi-zone binary phase pupil filters are not used in all calculations, i.e., ( ) 1 T θ = in Eqs. (1a) and (1b). Figure 3(a) and 3(b) show the intensity ( 2 E ) profiles along the transverse axis ( r ) at the focal plane ( 0 z = ) and the z axis for various size parameters 0 β , respectively. The intensity profiles are normalized by the intensity of the focal point ( 0 x y z = = = ). In this paper, the focal spot size and the depth of focus are defined as the FWHMs of the focal spot at the focal plane and at the z axis, respectively. Shown in Figs. 3(a) and 3(b), the smaller the 0 β is, the smaller the focal spot size is, whereas the longer the depth of focus is, which means that the smaller 0 β is beneficial to generate the needles of longitudinally polarized light [15]. However, when the 0 β is smaller than 3.58 for the radially polarized 5,1 LG beam, the focal spot has central peak intensity at the focal plane (see Fig. 3(a)), but start to be split along the z axis (see Fig. 3(b)). Namely, the intensity of the focal point is no longer the peak intensity along the z axis. Meanwhile, in terms of the definition of the focal plane, i.e., the plane with intensity being the peak intensity along the z axis, the actual focal plane will shift the original focal plane ( 0 z = ) when the 0 β is smaller than 3.58. In the actual focal plane, the focal spot size become bigger. Therefore, there is a discrepancy between the focal spot size reduced and the focal spot split along the z axis when decreasing the 0 β . When the imaging system illuminated by the radially polarized 5,1 LG beam, the 0 β cannot be smaller than 3.58. As the incident beam is truncated by the pupil of the imaging system, we define the effective power of the incident beam as the ratio of the power of the incident beam after truncated by the pupil to that of the incident beam before truncated by the pupil. Shown in , the focal spot size is about 0.45λ and slightly smaller than that for 0 3.68 β = . We note that the focal spot size is about 0.414λ and the effective power of the incident beam also reaches 90% for 0 3.36 β = (see Fig. 4). Therefore, when the radially polarized 5,1 LG beam is used to generate the small focal spot, one need not to focus on enhancing the effective power of the incident beam because the effective power of above 90% is enough for the real applications. Compared with the case of 0 3.68 β = , the focal spot size corresponding to 0 3.36 β = is decreased by 10%, which is attractive for the high resolution imaging and lithography. Unfortunately, the 0 3.36 β = cannot be used because the focal spot will be split along the z axis.
Based on the above discussions, we conclude that for the imaging system with NA 0.95 = illuminated by the radially polarized 5,1 LG beam, the smallest focal spot that can be achieved is about 0.45λ and the optimum size parameter 0 β of the incident beam should be 3.58. The smaller the 0 β is, the smaller the focal spot. However, when 0 3.58 β < , the focal spot will be split along the z axis. The natural question arises whether we can deal with the problem that the focal spot is split for the smaller 0 β by utilizing the multi-zone binary phase pupil filters.
In the following, we will deal with it. LG beam.

Generate the excellent focal spot with multi-zone binary phase pupil filters
In our optimizing procedure of the multi-zone binary phase pupil filters given before, the zone numbers of binary pupil filters are uncertain when the optimizing procedure starts. Before the optimizing procedure runs, we need to set the numbers N of rectangular belts and calculate the angles i θ that are determined by the used method of dividing the integral interval 0 θ α ≤ < of angle θ . The simplest way is to divide the integral interval 0 θ α ≤ < into N equal parts. We do not use this way in this paper. In order to get some rectangular belts with narrower width, we first divide the integral interval 0 θ α ≤ < into fifteen equal parts, and then divide the odd rectangular belts into four equal parts. So, the numbers of rectangular belts are 39 N = . In addition, as indicated before, during optimizing procedure, the merit parameters 30% S > , 20% r I < , and 1 G < must be met. Basing on the above optimizing procedure of the multi-zone binary phase pupil filters, for the 0 3.36 β = , we find the set of optimum structures of a seventeen-belt binary phase pupil filters: The transmission function ( ) T θ of this binary phase pupil filters is shown in Fig. 5(e).
Figures 5(a)-5(c) describe the intensity distribution of the focal spot in the rz plane for the total intensity, the transverse component, and the longitudinal component, respectively. Figure 5(d) describes the total intensity distribution of the focal spot along the r (red dashed curve) and z (blue solid curve) axes, respectively. The focal spot size (i.e., the FWHM of the total intensity of the focal spot at the focal plane) is 0.41λ and the depth of focus is 9.53λ . The Strehl ratio S and the sidelobe intensity r I are 41.75% and 16.35%, respectively. Compared with the case of no binary phase pupil filters (the focal spot size 0.414λ ), the focal spot size is almost unchanged, but the binary phase pupil filters designed by us successfully deal with the problem that the focal spot is split for 0 3.36 β = and make the depth of focus reach 9.53λ . Most importantly, even allowing the effective power of 90% for the radially polarized 5,1 LG beam with 0 3.36 β = , the conversion efficiency is still as high as 37.3%. In addition, the longitudinal component is dominant in this focal spot. To the best of our knowledge, this is the first report of the focal spot with subwavelength transverse size ( 0.41λ ), ultra long depth of focus ( 9.53λ ), and high conversion efficiency (37.3%) until now. We also do some optimizing designs for various 0 β and find that for the radially polarized 5,1 LG beam, the multi-zone binary phase pupil filters can extend effectively the depth of focus, but have no obvious effects on the decrease of the focal spot size. The focal spot size is majorly determined by the radially polarized 5,1 LG beam itself (i.e., the 0 β when