**π**/

**2**, arbitrary-length super-Gaussian optical needles are created with consistent beam size of

**0.36λ**(full width at half maximum) and the electric field being pure longitudinally polarized (polarization conversion efficiency greater than 99%). Numerical calculations show that the on-axis intensity distributions are super-Gaussian, and the peak-valley intensity fluctuations are all within 1% for

**4λ , 6λ , 8λ ,**and 10λ long light needles. The method remarkably improves the nondiffraction beam quality, compared with the subwavelength Gaussian light needle, which is generated by a narrow-width annular paraboloid mirror. Such a light beam may suit potential applications in particle acceleration, optical trapping, and microscopy.

## 1.

## Introduction

Extensive researches have been conducted in the past decade on the focusing^{1}2.3.4.5.6.7.^{–}^{8} and generation^{1}^{,}^{3}^{,}^{9}^{,}^{10} of cylindrical vector beams both theoretically^{1}^{,}^{2}^{,}^{4}5.6.7.^{–}^{8} and experimentally.^{1}^{,}^{3}^{,}^{9}^{,}^{10} Many applications have been reported, e.g., focusing of radially polarized vector beams in probing a tight focal spot^{1}2.3.^{–}^{4} and creating longitudinally polarized nondiffraction beams.^{5}6.7.^{–}^{8} Particularly, generation of subwavelength nondiffraction beams (called “needle of light” by Wang et al.^{5}) has attracted wide interest,^{6}7.^{–}^{8}^{,}^{11}12.13.14.15.16.17.^{–}^{18} because such a light beam suits a variety of applications in optical microlithography,^{19} high density optical data storage,^{20} or microscopic imaging.^{21} More importantly, the subwavelength nondiffraction beam with the electric field being longitudinally polarized has potential applications especially in particle acceleration,^{5}^{,}^{22} optical trapping,^{5}^{,}^{23} and near-field scanning optical microscopy.^{1}^{,}^{5}

Following the work by Wang et al.^{5} and again using radially polarized vector beams in free space, Kitamura et al. reported creating a subwavelength beam by focusing a narrow-width annular beam under a high numerical aperture (NA) aplanatic lens system.^{16} Full width at half maximum (FWHM) of the transverse beam width in the focal plane is about $0.4\lambda $ and the extended depth of focus is roughly $4\lambda $.^{16} Under a high aperture ($\mathrm{NA}\approx 1$) paraboloid mirror system, Dehez et al. further reported suppressing the nondiffraction beam to be $0.36\lambda $ in the focal plane and with an ultra-long focal depth of about $40\lambda $.^{17} The tight focusing performance of a high aperture paraboloid mirror was demonstrated by Meixner et al.^{24} These two methods^{16}^{,}^{17} generate longitudinally polarized light needles under lens or mirror systems, which are sharper than the result ($0.43\lambda $) in Ref. 5; however, the compressed light needles are nonuniform (approximately in a Gaussian shape) within the extended focal depth compared with Ref. 5, and the minimum beam size has only been localized near the focal plane, implies that the created light needles stringently diffracting propagate along the axial direction. Other researches have created either a flattop light needle of beam size broader than $0.41\lambda $^{6} or ultra-long light needles with much broader beam size and nonuniform intensity distribution.^{7}^{,}^{8} All these light needles are of dominantly longitudinal polarization, not purely especially at an out of focal plane. Besides, it is still possible to use some unconventional methods to generate light needles, e.g., using either a plasmonic lens^{14} or a super-oscillatory lens^{18}; however, these light needles are not longitudinally polarized and nonuniform, either within the near-field region^{14} or tens of wavelengths away from the element surface.^{18} So far, a light needle with radial beam size of $0.36\lambda $ (FWHM) might have approached the minimum size in free space under conventional far-field lens- or mirror-based systems ^{5}6.7.^{–}^{8}^{,}^{13}14.15.16.^{–}^{17}^{,}^{25}; however, at an out of focal plane, the transverse beam size varies largely and the nondiffraction beam oscillates or decays. Above all, generation of pure longitudinally polarized, super-Gaussian light needles of consistent beam size of $0.36\lambda $ and arbitrary length has not yet been achieved. As indicated by Wang et al.,^{5} a balance is difficult to draw among sharpening the beam size, extending the focal depth, homogenizing the axial intensity distribution, and purifying the polarization state. This problem still exists in Refs. 6–8 and 13–18. In this paper, we aim to solve this problem and describe a method to retain the beam size of $0.36\lambda $ over an arbitrarily long focal depth. We use an annular paraboloid mirror^{16}^{,}^{17} and modulate the incident radially polarized beam by the cosine synthesized filter (CSF). CSF is an amplitude filter, which can be designed quite flexibly to reshape the three-dimensional light field. We provide, e.g., light needles with consistent beam size of $0.36\lambda $ and super-Gaussian intensity distribution within an axial range of $4\lambda $, $6\lambda $, $8\lambda $, or over $10\lambda $, respectively.

## 2.

## Method

Suppose a radially polarized beam is incident upon a high NA reflection paraboloid mirror system (Fig. 1), and the radially and longitudinally polarized components of the electric field in the focal region are described, according to the vectorial Debye–Wolf diffraction integral,^{17}^{,}^{26}^{,}^{27} as

## (1)

$${E}_{r}^{(r)}(r,z)={A}_{r}{\int}_{0}^{\alpha}{l}_{0}(\theta )\frac{\mathrm{sin}(2\theta )}{1+\mathrm{cos}\theta}{J}_{1}(kr\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta )\phantom{\rule{0ex}{0ex}}\mathrm{exp}(-jkz\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\theta )\mathrm{d}\theta $$## (2)

$${E}_{z}^{(r)}(r,z)=-j{A}_{r}{\int}_{0}^{\alpha}{l}_{0}(\theta )\frac{2{\text{\hspace{0.17em}}\mathrm{sin}}^{2}\text{\hspace{0.17em}}\theta}{1+\mathrm{cos}\theta}{J}_{0}(kr\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta )\phantom{\rule{0ex}{0ex}}\mathrm{exp}(-jkz\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\theta )\mathrm{d}\theta ,$$^{3}4.

^{–}

^{5}

## (3)

$${E}_{r}^{(t)}(r,z)=-{A}_{t}{\int}_{0}^{\alpha}{l}_{0}(\theta ){\mathrm{cos}}^{1/2}\text{\hspace{0.17em}}\theta \text{\hspace{0.17em}}\mathrm{sin}(2\theta ){J}_{1}(kr\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta )\phantom{\rule{0ex}{0ex}}\mathrm{exp}(jkz\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\theta )\mathrm{d}\theta $$## (4)

$${E}_{z}^{(t)}(r,z)=-2j{A}_{t}{\int}_{0}^{\alpha}{l}_{0}(\theta ){\mathrm{cos}}^{1/2}\text{\hspace{0.17em}}\theta \text{\hspace{0.17em}}{\mathrm{sin}}^{2}\text{\hspace{0.17em}}\theta {J}_{0}(kr\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta )\phantom{\rule{0ex}{0ex}}\mathrm{exp}(jkz\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\theta )\mathrm{d}\theta ,$$^{2}

^{,}

^{26}

^{,}

^{27}

Specifically, consider a radially polarized Bessel–Gaussian beam with the waist plane at the pupil plane of an aplanatic lens system, and the corresponding amplitude of the Bessel–Gaussian beam is expressed as^{2}^{,}^{5}^{,}^{7}

## (5)

$${l}_{0}(\theta )=\mathrm{exp}[-{\beta}_{0}^{2}{\left(\frac{\mathrm{sin}\text{\hspace{0.17em}}\theta}{\mathrm{sin}\text{\hspace{0.17em}}\alpha}\right)}^{2}]{J}_{1}\left(2{\beta}_{0}\frac{\mathrm{sin}\text{\hspace{0.17em}}\theta}{\mathrm{sin}\text{\hspace{0.17em}}\alpha}\right),$$## (6)

$${l}_{0}(\theta )=\mathrm{exp}\{-{\beta}_{0}^{2}{\left[\frac{\mathrm{tan}(\theta /2)}{\mathrm{tan}(\alpha /2)}\right]}^{2}\}{J}_{1}\left[2{\beta}_{0}\frac{\mathrm{tan}(\theta /2)}{\mathrm{tan}(\alpha /2)}\right]$$^{5}

^{,}

^{7}Other radially polarized vector beams, e.g., Laguerre–Gaussian, could be similarly modeled and analyzed.

We now turn to reshaping the light field distribution in the focal region of a high aperture paraboloid mirror or aplanatic lens system. In order to create a subwavelength, arbitrarily long light needle with uniform intensity distribution, the CSF is formulized as

where $C(\theta )$ and $F(\theta )$ are constructed to separately suppress the radial beam size and elongate the axial focal depth. A variety of super-resolving pupil functions can, in principle, be used to replace $C(\theta )$^{28}; for a high aperture paraboloid mirror system, one might simply adopt the most widely used annular pupil filter, equivalently, using an annular beam, or an annular paraboloid mirror. The annular pupil function (or an annular beam) for a paraboloid mirror system, $C(\theta )$, is mathematically expressed as

## (8)

$$C(\theta )=\mathrm{circ}\left[\frac{\mathrm{tan}(\theta /2)}{\mathrm{tan}(\alpha /2)}\right]-\mathrm{circ}\left[\frac{\mathrm{tan}(\theta /2)}{\mathrm{tan}(\gamma /2)}\right],$$^{7}

^{,}

^{28}

## (9)

$$F(\theta )=\sum _{n=0}^{N}{a}_{n}\text{\hspace{0.17em}}\mathrm{cos}(kn\mu \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\theta ),$$^{7}

^{,}

^{28}though with much broader beam size and nonuniform axial intensity distribution.

^{7}The shift factor in $F(\theta )$, $\mu $, is to longitudinally shift the decomposed light segments in order to arbitrarily extend the focal depth by coherent superimposition. The balance factor, ${a}_{n}$, is to homogenize the superimposed beam such that a light needle with uniform intensity distribution might be obtained. $\tau $ is a normalization constant to maximize the light throughput.

The light field in the focal volume can be calculated by replacing ${l}_{0}(\theta )$ with ${l}_{0}(\theta ){T}_{\text{CSF}}(\theta )$ in Eqs. (1)(2)(3)–(4). The electric energy densities are calculated using ${|{E}_{r}|}^{2}$, ${|{E}_{z}|}^{2}$, and ${|\mathbf{E}|}^{2}={|{E}_{r}|}^{2}+{|{E}_{z}|}^{2}$. We further calculate the polarization conversion efficiency,^{5} $\eta {=\mathrm{\Phi}}_{z}{/(\mathrm{\Phi}}_{r}{+\mathrm{\Phi}}_{z})$, defined as the ratio of the longitudinally polarized electric energy to the total energy in the focal plane ($z=0$), with ${\mathrm{\Phi}}_{r,z}=\phantom{\rule{0ex}{0ex}}2\pi {\int}_{0}^{{r}_{0}}{|{E}_{r,z}(r,0)|}^{2}r\mathrm{d}r$, and ${r}_{0}$ denoting the radius of the central main lobe of ${|\mathbf{E}|}^{2}$. $\eta >50\%$ implies that the longitudinally polarized component predominates in the total energy. We compare the radial beam size between a high aperture reflection paraboloid mirror system and a transmission aplanatic lens system with respect to the obstruction factor, $\epsilon $, as shown in Fig. 2(a). The radial focal spot has been dramatically compressed for both mirror and lens systems with a large central obstruction, and is much sharper in the paraboloid mirror system for the same focusing angle. FWHM is used to characterize both the radial beam size and axial focal depth, denoted as ${r}_{\mathrm{FWHM}}$ and ${z}_{\mathrm{FWHM}}$, respectively. For NA 0.95 ($\alpha =71.8\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$) aplanatic lens system with clear pupil, ${r}_{\mathrm{FWHM}}=0.68\lambda $, ${z}_{\mathrm{FWHM}}=1.46\lambda $, and $\eta =45.0\%$. (Hence, the radially polarized parasitic energy dominates, as $1-\eta =55.0\%$.) While for a reflection paraboloid mirror system of $\alpha =71.8\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$, ${r}_{\mathrm{FWHM}}=0.56\lambda $, ${z}_{\mathrm{FWHM}}=1.53\lambda $, and $\eta =58.2\%$; for $\alpha =90\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$, parameters are calculated as ${r}_{\mathrm{FWHM}}=0.43\lambda $, ${z}_{\mathrm{FWHM}}=1.07\lambda $, and $\eta =79.9\%$, respectively. Specifically, the radial focal spot (${r}_{\mathrm{FWHM}}=0.37\lambda $) is much sharper with $\epsilon =0.7$, and compared with $\alpha =71.8\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{deg}$ lens in Fig. 2(b) and mirror system in Fig. 2(c), respectively. The polarization conversion efficiency is much higher for the paraboloid mirror system, i.e., the radially polarized component has been highly suppressed, as shown in Fig. 2(e); as a result, the electric field in the focal region is much more longitudinally polarized.

## 3.

## Numerical Calculations and Comparisons

For a high aperture paraboloid mirror system with maximum focusing semi-angle $\alpha $ being $\pi /2$, light needle with radial beam size of $0.36\lambda $ (FWHM) could be uniformly elongated with CSF, e.g., from $4\lambda $ to over $10\lambda $. Within the flattop range of the modulated light needles, the peak-valley intensity fluctuations are all less than 1%. Equations (1), (2), and (6)(7)(8)–(9) are used in the following examples. In order to highly compress the radially polarized component, the obstruction factor is set to be 0.7 for all cases, viz., almost half the area of the incident beam is blocked. It is found through calculations that the radial beam size of the elongated light needle remains $\sim 0.36\lambda $, even further increasing $\epsilon $. In order to show the proposed method is powerful and quite flexible, variable length light needles over $4\lambda $, $6\lambda $, $8\lambda $, and $10\lambda $, are designed, respectively, and the total energy densities, ${|\mathbf{E}|}^{2}$, are plotted in Fig. 3. Corresponding parameters in Fig. 3 ($\epsilon =0.7$) are tabulated in Table 1. The on-axis intensity distributions of the flattop light needles in Fig. 3 can be characterized by a series of super-Gaussian functions with different orders and widths,^{15} in contrast to the Gaussian shape, which is generated by an annular paraboloid mirror.^{16}^{,}^{17}

## Table 1

Parameters of cosine synthesized filter for various length flattop light needles (ε=0.7).

N | a0 | a1 | a2 | a3 | a4 | a5 | a6 | μ/λ |
---|---|---|---|---|---|---|---|---|

2 | 0.731 | 0.988 | 1 | —a | — | — | — | 1.240 |

4 | 0.530 | 0.900 | 1 | 0.673 | 0.635 | — | — | 0.805 |

5 | 0.531 | 0.690 | 1 | 0.673 | 0.610 | 0.550 | — | 0.950 |

6 | 0.568 | 0.993 | 0.971 | 1 | 0.826 | 0.606 | 0.605 | 0.916 |

## a

— denotes that an is invalid for n larger than N.

For the light needle in Fig. 3(d), evaluation parameters are calculated as ${r}_{\mathrm{FWHM}}=0.36\lambda $, ${z}_{\mathrm{FWHM}}=14\lambda $, $\eta =99.9\%$, and the flattop peak-valley intensity fluctuation is less than 0.7% within an ultra-long distance of over $10\lambda $ (Fig. 4), which is propagating without any divergence, as shown in Fig. 4(c) and 4(d), and localized into a small beam area of $\sim 0.10{\lambda}^{2}$ (area at half intensity). If $\lambda =405\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$, e.g., the radial beam size and the beam area are merely 145.8 nm and 0.0164 $0.0164\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mu \mathrm{m}}^{2}$, respectively. The corresponding amplitude distribution of ${T}_{\text{CSF}}$ is plotted in Fig. 4(a), which is valid within [69.984 deg, 90 deg]. Electric energy densities, ${|{E}_{r}|}^{2}$, ${|{E}_{z}|}^{2}$, and ${|\mathbf{E}|}^{2}$, are plotted in Fig. 4(b) in the focal plane ($z=0$), where the radially polarized component almost disappears (see the inset figure). When compared with the previous reports, the light needle (with radial beam size of $0.36\lambda $ and as long as $\sim 14\lambda $), as shown in Figs. 3(d) and 4(c), is much sharper and remarkably longer than the result (with radial beam size of $0.43\lambda $ and as long as $\sim 4\lambda $)^{5}; again, the light needle is 2.5 times sharper than the report (with radial beam size of $0.9\lambda $),^{8} and much more uniform over the whole axial range of $\sim 14\lambda $ (FWHM)^{8}; moreover, the light needle is pure longitudinally polarized ($\eta >99\%$) within the flattop range of the light needle in contrast to Refs. 5–8 and 13–18. Lastly, when compared with Ref. 17, in the focal plane, both light needles are localized to be a small beam size of $0.36\lambda $; however, at an out of focal plane, the light needle, as shown in Fig. 4, remains the beam size of $0.36\lambda $. Such a good balance has not been achieved in Ref. 17. The detail of the polarization state is further plotted in Fig. 4(d). Above all, the modulated light needle keeps the minimum possible beam size among the known reports ^{5}6.7.^{–}^{8}^{,}^{13}14.15.16.17.^{–}^{18} and draws a good balance among sharpening the beam size, extending the focal depth, homogenizing the axial intensity distribution, and purifying the polarization state.

It should be indicated that it is impossible to modulate an ultra-long optical needle from the proposed method. As the theoretical basis, i.e., the vectorial Debye–Wolf interal,^{26} it is only valid to accurately describe the light field in the vicinity of the focus (far-field approximation condition). According to the previous researches,^{2}^{,}^{5}6.7.^{–}^{8}^{,}^{16}^{,}^{17} the length of the light needle, which can be modulated through the proposed method, could vary from several wavelengths to tens of wavelengths. When the amplitude distribution of the incident radially polarized vector beam changes, the structural parameters of CSF must be re-optimized in the same manner.

Practically, there might be two challenges for creating such light needles:^{24}^{,}^{29} first, a high NA paraboloid mirror with sufficiently high surface accuracy, e.g., peak-valley deviation better than $\lambda /10$, is required;^{24} however, this requirement might be alleviated if adaptive optics is introduced to correct the wavefront aberration induced by the imperfection of the fabricated paraboloid mirror; second, accurate realization of the amplitude transmission (${T}_{\mathrm{CSF}}$) depends on the state-of-the-art coating and photolithography techniques;^{29} alternatively, a method based on binary computer holography might be promising, where a ferroelectric liquid spatial light modulator with accurate binary amplitude modulation is used, and as a result, the accurate vector pupil apodization might be generated.^{30}

## 4.

## Conclusion

In conclusion, a method is presented to generate an arbitrary-length super-Gaussian (flattop) optical needle with consistent beam size of $0.36\lambda $ (FWHM) and the electric field being pure longitudinally polarized (polarization conversion efficiency greater than 99%). It is realized through modulating the Bessel–Gaussian radially polarized vector beam by the CSF under a reflection paraboloid mirror system with maximum focusing semi-angle of $\pi /2$. Numerical calculations show that the on-axis intensity distribution of the modulated light needle is super-Gaussian, in contrast to the subwavelength Gaussian light needle, which is generated by a narrow-width annular paraboloid mirror. Such a light beam may suit several potential applications in particle acceleration, optical trapping, and near-field scanning microscopy.^{1}^{,}^{5}^{,}^{10}

## Acknowledgments

The authors are grateful to the funding from National Natural Science Foundation of China (NSFC) under Grants Nos. 51275121 and 61008039.

## References

## Biography

**Tao Liu** obtained his BS degree in 2007 from Harbin Institute of Technology. Currently, he is a PhD candidate at Harbin Institute of Technology. His research interests are superresolution focusing theory and superresolution confocal microscopy.

**Jiubin Tan** was awarded the BS, MS, and PhD degrees in 1982, 1988, and 1991, respectively, from the Department of Automation Measurement and Control, Harbin Institute of Technology. In 1995, he became a professor of measurement and instrumentation. He was the former dean of the Department of Automation Measurement and Control. He is currently the director of the Center of Ultra-precision Optoelectronic Instrument Engineering, Harbin Institute of Technology. His research is broad in the fields of ultra-precision measurement and instrumentation.

**Jie Lin** received his doctorate in optics from Harbin Institute of Technology in 2007. He currently works as an associate professor in the discipline of instrument science and technology at Harbin Institute of Technology. His research interests include optical beam shaping, super-resolution focusing, two-dimensional grating for lithography, binary optical elements, and the optical performance of metal micro/nanostructure.

**Jian Liu** obtained his BS degree from Southwest Jiaotong University in 1997. He was then awarded the MS and PhD degrees from Harbin Institute of Technology in 2002 and 2009, respectively. In 2012, he became a professor in the Department of Automation Measurement and Control, Harbin Institute of Technology. His research interests are superresolution light microscopy, precision measurement, and metrology.