Abstract
1. INTRODUCTION
Isotropic resolution is progressively desirable in laser fabrication and microscopy imaging [1–15], while unfortunately, conventional focal region is inherently of anisotropic configuration subjected to optical lens design and the diffraction nature of light. For example, a plane wave is commonly transformed into an elliptical focal spot, with its axial size more than 2 times its lateral size, after a lens. Tremendous attention has been paid to increase the axial resolution in recent decades. The total internal reflection fluorescence microscopy uses evanescent wave illumination and provides high-contrast axial superresolution image of a thin layer of object within the evanescent wave penetration depth. In the case of acquiring a larger axial field of view, 3D scanning methods with a compressed point spread function are dominating currently, including stimulated emission depletion (STED) microscopy, two-photon/multi-photon microscopy, confocal microscopy, technique, and the simultaneous spatial temporal focusing approach [16–24]. Implementations of these methods rely on their unique optical configurations. STED incorporates two parallel laser beams in illumination, one of which is modulated by a vortex phase plate so as to form a hollow focus to deactivate the fluorescent signal generated by the other beam. Confocal configuration blocks the signal out of focus by placing a pinhole in the conjugate position of the illumination focus. technique utilizes two well aligned head-to-head objectives, which is particularly suited to transparent and thin samples. The simultaneous spatial temporal focusing approach uses a chromatically dispersed femtosecond laser in illumination.
In practice, an axial superresolution method with compact configuration, extraordinary compatibility, and flexibility that can be easily transferred to the other microscope systems to enable a better performance, is immensely preferred. Beam modulation by directly inserting a mask in the pupil plane to squeeze the focal spot axially possesses the features above. Different kinds of modulation masks have been reported [25–28], including the binary phase-only mask and the annular amplitude mask. These pupil masks, demonstrated previously, were designed either empirically or by search algorithms, as far as we know. The Toraldo concept [28] has established a Fourier transform relationship between the axial amplitude of the focus and the pupil transmittance distribution, which nonlinearly relates to the pupil radius and thus makes the pupil mask difficult to fabricate in practice. Lack of a universal theoretical model to analyze the pupil modulation-based axial superresolved focusing has largely burdened mask design and hindered achieving the optimal resolution.
In this paper, we analyze the axial superresolved focusing process from an interference perspective and find that the axial resolution closely relates to the gradient of the Gouy phase, which is used to describe the phase anomaly of a focused field with respect to an ideal spherical wave. Based on this finding, we propose an axial interference model and develop a Gouy phase steering method. This method can be universally applied to varied types of incident polarization and objective numerical apertures (NAs) and can be transplanted to optimize multiple kinds of pupil function for axial superresolution. It also shows the potential to attain a superresolved focal length that can even be extremely shortened, similar to the lateral resolution modulation based on superoscillation [29,30]. We use this method to optimize a shaded-ring pupil mask, i.e., a binary amplitude modulation mask. The approach has been implemented on the linearly polarized beam and radially polarized beam, respectively, to achieve the axial superresolution. The results have clearly evidenced its superiority to previous studies on pupil masks designed by experience or by a searching algorithm. We show that both in simulations and experiments, the axial resolution of a radially polarized beam can be 27% enhanced from the diffraction limit and can be 11% finer than the previously reported record without increasing the sidelobe intensity [25]. Additionally, by introducing more rings in the pupil mask based on this method, the energy residing in the sidelobes can be further suppressed and precisely redistributed according to their influence in the specific applications. A spherical focus, with a ratio of 1:1:1.06 between the full widths at half-maximum (FWHMs) in the , , and directions, is also demonstrated in simulation for a radially polarized beam, which evidences the isotropic resolution by using a compact and simple single objective configuration. Finally, by applying the method to design a phase ring mask, a superoscillation-like behavior can be observed in the axial dimension with the conceptual possibility toward extremely high axial resolution.
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
2. MECHANISM AND SIMULATIONS
In the axial superresolution method concerning the Gouy phase gradient (GPG), the pupil mask [] is decomposed to a clear pupil (constant distribution of 1) and an inverse mask []. Optimization of the inverse mask [] is equivalent to the optimization of the pupil mask []. Based on the Huygens–Fresnel principle, the final focus can be treated as an interference between the focus formed with the clear pupil (denoted as “clear focus” in the following) and the focus formed with the inverse pupil mask (denoted as “inverse focus” in the following). The completely destructive interference between the clear focus and the inverse focus along the axis can shorten the axial focal depth, which provides a theoretical constraint condition for the parameter optimization in the pupil plane. To demonstrate this interference concept and the mask design mechanism, we mainly take use of a shaded-ring mask model as an example, in which a ring zone is amplitude-attenuated [25,27,31] [Fig. 1(a)], owing to its usually enabling axially superresolved focus accompanied with sidelobes of low intensity. It should be emphasized that this method does not merely work for the shaded-ring mask design.
Figure 1.Schematic diagram of the pupil mask design for axially superresolved focusing based on GPG tuning. (a) Decomposition of a shaded-ring pupil mask [
For a shaded-ring pupil mask, the inverse pupil mask [] is an amplitude-attenuated, -phase ring-pass mask [Fig. 1(a)]. Based on the proposed interference model, the parameters, including the ring radius and the amplitude attenuation, can be designed and calculated, rather than using a search algorithm. To graphically show the effect of the interference on axial superresolved focusing, we plot the amplitude and phase distribution of the clear focus and the inverse focus along the axis [Figs. 1(b) and 1(c)]. Coordinates of the focus center are . To ensure the completely destructive interference between the clear focus and the inverse focus, phase difference of and equal amplitude near should be satisfied simultaneously between the two foci. At the focus center where , the phases of the two foci differ by [Fig. 1(c)] because the phases at the clear pupil and the inverse pupil differ by , and then the light propagates an equal distance from the pupil plane to the focus center. While the amplitude of the clear focus is larger than that of the inverse focus, at a certain distance from the focus center, equal amplitude can be satisfied [shown by gray dotted lines in Fig. 1(b)]. Only if the phases of the two foci still differ by at that position apart from the focus center, can the completely destructive interference occur. Consequently, the slope of the Gouy phase distribution [see Eq. (1)] along the axis should be the same between the clear focus and the inverse focus, to enable the unchanged phase difference of [shown in Fig. 1(c)]. We call it the equal Gouy phase gradient (EGPG) condition, which is a necessary condition for the completely destructive interference between the clear focus and the inverse focus to have an axially superresolved focus. Figure 6 in Appendix B shows the simulated microscopy images formed by the clear focus and superresolution focus, respectively, for a clearer demonstration. The EGPG condition is satisfied in Figs. 1(b)–1(d), and a completely destructive interference can be observed at the positions labeled by the two dark dotted lines, leading to a narrower width of the main lobe than the unmodulated focus [Fig. 1(d)]. Inevitably, the focusing efficiency in proportion to the energy residing in the main lobe would be decreased due to the destructive interferences in the main focal area.
Based on this understanding, we proposed a GPG tuning model to design the pupil mask for axial superresolved focusing. The GPG of the inverse focus at [denoted as ] can be calculated by Eq. (1). The calculation approach here can be understood as the tilted wave average method. Different from the vectorial-diffraction-theory-based approach to calculate the exact complex amplitude, the phase is analyzed by averaging the phase retardations of the decomposed tilted waves by their normalized amplitude [32] (see Appendix I for the detailed deduction). To conveniently tune the GPG of the inverse focus, we introduce a tool function of the convergence angle, and then the GPG of the inverse focus can be calculated by Eq. (2). The tool function is shown in Eq. (3), whose expression relates to the target value (denoted as ) to which the GPG of the inverse focus is to be tuned,
In order to tune the GPG of the inverse focus to the value of needs to be equal to , which ensures in Eq. (2). The EGPG condition can be satisfied by making , where denotes the GPG of the clear focus.
There are three steps in the GPG tuning-based method to calculate the parameters of a shaded-ring mask.
Figure 2.Calculation of the inner and outer radii of the inverse mask to satisfy EGPG condition. (a) Schematic of the inverse mask for a shaded-ring pupil mask; (b) radius determination tool function; (c) GPG of the inverse focus at focus center, which maintains at GPG of the clear focus, i.e.,
When the GPG of the inverse focus is tuned to be equal to the GPG of the clear focus, the completely destructive interference could occur within the Abbe diffraction length, and the minimum axial FWHM tends to appear. The GPG tuning method can be applied to various types of polarization. Here, we take a radially polarized beam as an example to illustrate the influence of the inverse focus GPG on the final axial resolution, as well as to show the mask design process in detail. The axial FWHMs of a radially polarized beam focus with the GPG of the inverse focus tuned at different values are plotted for comparison, as shown in Fig. 1(e). The power transmittance of the shaded-ring mask is set at 10%, 20%, 40%, and 60%, respectively (corresponding to different colors), and the sidelobe tolerance of all the calculated superresolved foci is set as 20%. The data set plotted in each color presents the axial FWHMs with the same power efficiency and sidelobe intensity, but for different GPGs of the inverse focus. When the power transmittance is large, the minimum axial FWHM is achieved with the GPG of the inverse focus vaguely less than the clear focus. In such a case, the axial FWHM changes slightly, which reveals that in a condition with high transmittance, the GPG of the inverse focus equal to the clear focus can give rise to an FWHM rather close to the minimum value. When the power transmittance decreases from 60% to 10%, a smaller axial FWHM can be achieved, and the GPG of the inverse focus that corresponds to the minimum axial FWHM converges to the GPG of the clear focus. It is noteworthy that, in the case of a smaller final focus, which is expected, the distance between the two completely destructive interfered points [labeled by gray dotted lines in Figs. 1(b)–1(d)] is closer, and the Gouy phase can be better approximated to a linear distribution; then the EGPG condition should be satisfied more critically as a consequence.
Based on the analysis above, we can tune the inverse focus GPG to the value of the clear focus GPG during the mask design. Figure 2(b) plots the constructed function value versus the convergence angle, by making in Eq. (3). The and should be determined corresponding to the equal function value, as illustrated by the gray dotted lines. Figure 2(c) shows the GPG of the inverse focus, which is maintained at the GPG of the clear focus with these calculated and . We randomly select two pairs of inner and outer convergence angles in Figs. 2(b) and 2(c), as shown by “Ring I” and “Ring II” in the figure, and plot the intensity and phase distribution of the focus, as shown in Figs. 2(d) and 2(e). The FWHM of the focus intensity distribution in the axial direction decreases with a wider ring, while the phase distribution on the optical axis changes slightly, and the phase gradient is unchanged at the focal point.
For each pair of and (or and ), the amplitude attenuation can be correspondingly obtained, given the power transmittance of the mask. Low transmittance of the total pupil can be overcome in practice by increasing the input power. In the meantime, the intensity of sidelobes, which is customarily pursued at a low level in microscope applications, should be more carefully controlled with a sophisticated pupil modulation design. The axial FWHMs and sidelobes of a radially polarized beam focus corresponding to 30% pupil power transmittance are shown in Fig. 3(a). Smaller FWHM comes with higher sidelobe intensity. A 0.75--FWHM can be achieved with 30% sidelobe intensity, as shown in Figs. 3(b) and 3(c) (pupil mask parameters: , ; and the amplitude attenuation of the ring is ). The focusing efficiency in the modulation is 6.87%. The lateral resolution could be slightly degraded when accompanied with the axial resolution improvement, as shown in Fig. 3(d). The axial FWHM is 27% improved from the 1.03--FWHM of the diffraction-limited focus and 11% finer than the FWHM record of 0.84-, with 30% sidelobe reported in the previous simulation work [25].
Figure 3.(a) Axial FWHM and sidelobe of the radially polarized beam focus with the shaded-ring pupil mask whose parameters are calculated based on EGPG; energy transmittance of the mask: 30%; (b)–(d) 2D intensity distribution in
To achieve a minimized quasi-spherical focus, we further used a lower power transmission (20%) and a higher sidelobe intensity (50%) for the radially polarized beam in our simulations. An axial superresolution of , 36% enhanced from the diffraction limit can be realized under the EGPG condition (pupil mask parameters: ; ; amplitude attenuation of the ring is ). Figure 3(e) shows the intensity distribution of the focus in the cross section, and Fig. 3(f) plots the intensity profile along and axis, respectively. The ratio between the FWHM () and the FWHM () is 1.06:1. The isotropy property of the focus is superior to the traditionally elliptical focus (with a typical aspect ratio between axial and lateral ). The maximum sidelobe intensity aside the optical axis is 64%, which is mainly contributed by the transversely polarized component. Figure 3(g) shows the pure longitudinal (i.e., -polarized) component intensity distribution with a much cleaner background. With a special photosensitive layer prepared to be sensitive only to this longitudinal field component, the strong transversely polarized sidelobes can be suppressed [33–35]. On the other hand, the influence of the strong sidelobe distant from the focus center could be conceivably eliminated to some extent by using confocal configuration, or by multiphoton absorption and threshold effect.
3. EXPERIMENTS
In the experiment, we utilize a phase spatial light modulator (SLM) (HED6010-NIR-011-C, HOLOEYE) to achieve the amplitude-only modulation in the pupil plane. A laser (MDL-C-405-50 mW, CNI CW Laser or Coherent Chameleon Ultra II Femtosecond Laser) is first filtered and expanded by an assembly of a lens (), a pinhole (10 μm in diameter, P10D, Thorlabs), and another lens () to obtain a collimated laser beam. The polarization direction of the beam is 45° with respect to the orientation of the liquid crystal in the SLM. The phase on the SLM is set to 0 or . The polarization direction of the beam is rotated by 90° in the region. By adjusting the angle of a polarizer after the SLM, the ratio between the amplitude of the regions corresponding to the 0 and phase on the SLM can be modulated to achieve the simulated shaded-ring mask, as illustrated in Figs. 4(a) and 4(b). A radial polarization converter (RADPOL4L, ARCoptix) is utilized to generate the radially polarized beam. The beam is tightly focused by a 1.4 NA oil-immersed objective (UPlanSApo, Oil, Olympus). An Au nanoparticle (80 nm or 150 nm in diameter) on a 3D nanometer-precision motorized positioning system (P-563.3CD, PI) is used as a point object to measure the intensity distribution of the modulated focus. The scattered light is collected by a multimode fiber in a confocal configuration connected with an avalanche photodiode (APD) (SPCM-AQRH-13-FC, EXCELITAS) and the wide-field images are captured by a CCD camera (MER-130-30UM, Daheng Optics). The overall view of the experimental setup can be found in Fig. 8 in Appendix D.
Figure 4.Experiment on the axially superresolved focusing based on the GPG tuning. (a) Schematic diagram of converting the phase modulation to polarization direction modulation. The amplitude modulation can then be realized by adjusting the direction of polarizer after the SLM. (b) Image of the intensity distribution in the pupil plane; (c), (d) measured
To show the applicability of this mask design method on various polarization beams, the experiments have been implemented for two types of CW laser beams of radial and linear polarization, and the experimental results turn out to be very close to the simulation results. The shaded-ring mask with parameters corresponding to the simulations in Figs. 3(b) and 3(c) have been realized by our experimental setup (, , and the amplitude attenuation of the ring is 0.855). The intensity distribution in the pupil plane is shown in Fig. 4(b). Figure 4(c) shows a measured focus intensity profile of a radially polarized beam (the raw data are in supplementary
4. APPLICATION OF EGPG METHOD ON DIFFERENT TYPES OF PUPIL FUNCTION
Figure 5.Application of the EGPG method on a
We categorize axial focusing behaviors via the EGPG method into three zones: I, diffraction limited; II, superresolved; and III, deep superresolved, depending on whether the axial FWHM of the main lobe is finer than the axial diffraction limit () or the limit of the lateral resolution with a uniform illumination, i.e., half of the media wavelength (). Although the infinitesimal main lobe is possible in Zone III, the cost is the exponentially increased sidelobe intensity, similar to the lateral resolution modulation by superoscillation in an extreme case. For practical applications, a better trade-off is gained in Zone II. Figure 5(d) presents the 3D intensity isosurface distribution of focus 3 in Zone II, in which the axial FWHM is shortened by half from the diffraction limit. Its 2D and cross section intensity distribution can be found in Fig. 11 in Appendix G. The lateral resolution is anisotropic in such a modulated linearly polarized tight-focus model. The axial sidelobe intensity is 7.65 times that of the main lobe, while the influence of the sidelobe may be conceivably minimized when the focus is utilized to deterministically manipulate the thickness of ultrathin 2D nanomaterials and fabricate functional integrated flat optoelectronic devices [36]. The vector property distribution of a deep superresolved focus with an ultrasmall FWHM of and accompanied by ultralarge sidelobes is presented in Fig. 12 in Appendix H to show the extremity of the modulation from a conceptual perspective. More details about the applications of the EGPG method to design a phase mask to modulate a radially polarized beam (with laterally isotropic focus shape) can be found in Fig. 13 in Appendix I.
5. CONCLUSION AND OUTLOOK
In conclusion, we have proposed a GPG tuning-based method with an interference model for axial superresolved focusing via beam modulation in the pupil, which possesses the merits of compact configuration and general applicability on the existing (superresolved or traditional) microscope system such as the confocal or the two-photon systems to upgrade their axial resolution. The pupil mask is decomposed to a clear pupil and an inverse mask, and the optimal axial resolution can be attained by utilizing the completely destructive interference between the two corresponding foci. The axial resolution of a shaded-ring mask modulated radially polarized beam can be 27% improved from the diffraction limit and 11% improved from the previous resolution record. A quasi-spherical focus with an aspect ratio of 1:1:1.06 for , , and FWHMs has been demonstrated in simulation, showing the potential to achieve the isotropic resolution using a simple single-objective configuration. Moreover, a superoscillation-like behavior has been observed in the axial dimension with the conceptual possibility toward extremely high axial resolution. The method has been proved to be efficient for multiple types of polarizations and pupil masks. By introducing more rings in the pupil mask, the sidelobe intensity distribution could be manipulated for adaptation to specific applications. In the future, the SLM could be substituted by a plate fabricated using photolithography and film coating techniques, to improve the axial resolution of a microscope system simply by inserting the thin plate in the pupil plane of the objective. The costs of the implementation include lowered focusing efficiency, slightly degraded lateral resolution, and increased sidelobe intensity, which could be overcome by increasing the source energy, introducing more modulation freedoms, and using an optical system that can suppress the defocused energy. The specific applications that can be explored in the future may include the 3D omnidirectionally superresolved microscopy with a confocal system, and the omnidirectional laser direct writing of a circular waveguide by combining the threshold effect of the photoresist, so as to suppress the strong sidelobes, and so forth. In a word, the research may benefit a wide spectrum of applications, such as 3D laser fabrication, superresolution imaging, and light trapping, which increasingly require axial superresolved and isotropic focusing properties.
Acknowledgment
Acknowledgment. We thank Dr. Dejiao Hu, Dr. Sicong Wang, Dr. Fei Qin, and Dr. Shichao Song for helpful discussions.
APPENDIX A: DEDUCTION OF THE GOUY PHASE GRADIENT OF THE INVERSE FOCUS
Considering the vector focusing theory, the complex field along the optical axis of the inverse focus can be written as Eq. (
We simplify Eqs. (
In the vicinity of the focus center, , and can be approximated to . So we have Eq. (
The Gouy phase of the inverse focus can be calculated by Eq. (
The GPG of the inverse focus can be calculated by Eq. (
APPENDIX B: UTILIZATION OF EGPG METHOD TO UPGRADE THE AXIAL RESOLUTION OF A CONFOCAL SYSTEM
Figure 6.Images in the
APPENDIX C: SIDELOBE INTENSITY SUPPRESSION AND MODULATION BY INTRODUCING MORE RINGS
Figure 7.Normalized intensity profile along
APPENDIX D: SCHEMATIC OF THE EXPERIMENTAL SETUP
Figure 8.Schematic of experimental setup. BS, beam splitter; SLM, spatial light modulator.
APPENDIX E: AXIAL SUPERRESOLVED FOCUS MODULATION FOR A FEMTOSECOND LASER
Figure 9.2D focus intensity distribution in the
APPENDIX F: SUPEROSCILLATION-LIKE BEHAVIOR OF THE AXIAL RESOLUTION MODULATION BASED ON EGPG
Figure 10.(a) Transient electric field and (b) the local wave vector distribution along the
APPENDIX G: INTENSITY DISTRIBUTION OF A LINEARLY POLARIZED FOCUS WITH TWO-FOLD IMPROVED AXIAL RESOLUTION
Figure 11.Intensity distribution in (a)
APPENDIX H: VECTOR PROPERTY OF VECTOR PROPERTY DISTRIBUTION OF A DEEP SUPERRESOLVED FOCUS WITH AN ULTRASMALL FWHM OF 0.1λ
Figure 12.Simulation of a deep superresolved focus of linearly polarized beam with a
APPENDIX I: APPLICATIONS OF THE EGPG METHOD TO DESIGN A π PHASE MASK TO MODULATE A RADIALLY POLARIZED BEAM
Figure 13.Applications of the EGPG method to design a
References
Set citation alerts for the article
Please enter your email address