Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Compact structured illumination microscopy with high spatial frequency diffractive lattice patterns

Open Access Open Access

Abstract

Structured illumination microscopy (SIM) enables live-cell super-resolution imaging with wide field of view (FOV) and high imaging speed, but the illumination system is usually bulky. With the advantages of small structure and high efficiency, lattice patterns assisted by diffractive optical elements (DOEs) have been used for structured illumination in SIM. But it is still challenging to raise the spatial frequency of diffractive lattice patterns when using traditional DOE design method, and thus the super-resolution imaging performance is restricted. In this paper, we propose a novel design method for DOE to generate lattice patterns with spatial frequency close to the cut-off frequency. It is the first time to obtain a lattice pattern with such high spatial frequency by diffractive optics. Finally, the proposed SIM achieves a lateral resolution of 131 nm at 519 nm fluorescent light while maintaining an original size as a standard inverted fluorescence microscope by only inserting a single well-designed DOE in the illumination optical path, which may promote this compact SIM applied in super-resolution imaging field.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

The spatial resolution of optical microscopy has been restricted to about half of the wavelength because of the diffraction limit [1]. However, several super-resolution methods have developed in recent years, bypassing this barrier and extending the lateral resolution to tens of nanometers [25]. Although not providing better resolution than other super-resolution approaches, structured illumination microscopy (SIM) is considered as one of the most promising candidates in live-cell imaging area for its advantages of wide field of view (FOV), high imaging speed, compatibility with standard labels and low phototoxicity [610].

In SIM, when the sample is illuminated by a periodic pattern, in mathematical way it represents the sample’s combining with the illumination pattern multiplicatively in spatial space, as well as their convolution operation in spatial frequency space. After convolved by a series of delta functions that are Fourier transforms of the periodic patterns, the high spatial frequency information unresolved before will be shifted into the low spatial frequency scope and then is made available within the band-limited imaging system. Because different spatial frequency components are mixed in the images acquired with detector, a series of low-resolution images with different angles and phases are collected to synthesize the super-resolution image. The frequency down conversion determined by the pattern period is also constrained by the diffraction limit, so linear SIM doubles the resolution of conventional microscopy [5].

Besides conventional line grid patterns, two-dimensional (2D) lattices have been employed for structured illumination because the 2D lattice patterns only require lateral phase shift and provide a better signal-to-noise ratio theoretically [1114]. Such lattice patterns are generated mainly through interference method and a bulky 4f projection system is inevitable. The diffraction orders of a grating are diffracted onto a median plane where zero and other orders are filtered out and only the first orders are projected onto the rear aperture of the objective lens. The spatial frequency of the pattern gets higher when the first diffraction orders approach nearer opposite edges of the rear aperture.

Diffractive optical elements (DOEs) are also suitable for generating lattice patterns due to its high efficiency and design flexibility [1517], however the spatial frequency is not easy to raise in the same way as interference method. Iterative methods such as modified Gerchberg-Saxton (G-S) algorithms as well as other non-iterative methods also have been developed for generating various lattice patterns [1820]. Nevertheless, these methods including our previous work [21] all restrict the period of patterns no smaller than the size of Airy spot, which means the spatial frequency of lattice patterns cannot exceed 0.58 fco, where fco is the cutoff frequency and can be seen as the frequency of the fringes interfered by two ideal point sources at the edge of the entrance pupil. Although some methods have been proposed to generate sub-diffraction patterns [2224], they only generated a single very small spot or a few periods rather than a lattice pattern. Therefore, it is necessary to propose a new method for designing DOE to raise the spatial frequency of the lattice patterns and then increase the imaging resolution.

Here, we propose a new method for DOE design to achieve high spatial frequency lattice patterns. Firstly, through a mathematical analysis we explain how to choose a propriate initial phase in detail to keep the subsequent iterative process away from divergence. Moreover, we modified the G-S algorithm to optimize the uniformity and efficiency of the lattice patterns. Using this method, we can obtain the illumination pattern with a high spatial frequency of 0.9fco. Generating such high spatial frequency lattice pattern with DOE, first reported to our knowledge, is the key factor to realize high resolution for SIM. We also further establish the system with the epi-illumination where the phase shifts are achieved through the light source movement and the DOE is mounted on the filter cube to provide a better support for oil immersion objective lens. As a result, the size of this compact illumination system is 75 mm × 75 mm × 60 mm, which is about 2 orders of magnitude smaller than that of existing commercial SIM [25]. Lastly, imaging comparisons, using both 100 nm fluorescence beads and biological specimens, indicate obvious resolution enhancement with the compact SIM system.

2. Design of DOE

Here we present an iterative method to generate high spatial frequency lattice patterns. Considering high spatial frequency of the target pattern in this case, we should reduce the sampling interval in the output plane. But as the sampling interval decreases, the number of points to be optimized increases, which makes the optimization difficult to be convergent with the GS algorithm and its existing modified versions [26]. The iterative process is sensitive to the initial phase choice, so it is possible to improve the optimization process by setting different starting points [27,28]. Instead of random phase in most iterative algorithms for designing DOE, we choose an explicit quadratic phase as the initial value to keep the subsequent iterative process away from encountering a local stagnation and assure a higher convergence speed. In what follows, we will present how to choose the propriate initial phase.

Assuming N represents the number of frequencies in the lattice pattern, for example, N is two for a square lattice and three for a hexagonal lattice. The ideal amplitude of lattice pattern can be expressed as

$${a_{pt}}({\mathbf r}) = \sum\limits_{n = 1}^N {\cos (} 2\mathrm{\pi }{{\mathbf f}_{ptn}} \cdot {\mathbf r} - {\alpha _0}^{(n)}),$$
where fptn (n = 1, 2, …, N) represents the n-th spatial frequency vector, r is the spatial position vector, and α0(n) is the constant relative phase. And the corresponding intensity pattern ipt(r) = [apt(r)]2. In coordinate form spatial frequency vector f can be written as (fx, fy) and spatial position vector r can be written as (x, y). Here we denote the light field distribution in the sample plane as ep(r) = ap(r)d(r)exp [iφ(r)], where ap(r) and φ(r) are respectively the amplitude and phase distribution, d(r) is a binary mask function associated with FOV. The DOE could be considered as recording the Fourier transform of the complex field ep(r). Thus, the distribution in the DOE plane can be derived via the Fourier transform as
$$\begin{array}{l} E({\mathbf f}) = \mathrm{{\cal F}}\{ {e_p}({\mathbf r})\} = \mathrm{{\cal F}}\{ {a_p}({\mathbf r}) \times d({\mathbf r})\exp [i\varphi ({\mathbf r})]\} \\ \;\;\;\;\;\;\; = \mathrm{{\cal F}}\{ {a_p}({\mathbf r})\} \otimes \mathrm{{\cal F}}\{ d({\mathbf r})\exp [i\varphi ({\mathbf r})]\} = P({\mathbf f}) \otimes D({\mathbf f}) \end{array}$$
where⊗denotes the convolution operation, P(f) and D(f) are the Fourier transform of ap(r) and $d({\mathbf r})\exp [i\varphi ({\mathbf r})]$ respectively. Considering the amplitude in the sample plane as the ideal pattern, we further derive the distribution in the DOE plane as
$$E({\mathbf f}) = \sum\limits_{n = 1}^N {[D({\mathbf f} - {{\mathbf f}_{ptn}}) + D({\mathbf f} + {{\mathbf f}_{ptn}})} ].$$

The fact that E(f) is phase-only and constrained by the system cut-off frequency, leads to D(f) also phase-only and frequency constrained. From this point of view, our purpose is to design a proper phase φ(r) superimposing on the pattern and keep the Fourier transform of the complex field ep(r) phase-only and band-limited. It is consistent with the flat-top beam shaping studied by some researchers before [29,30].

As described in Ref. [31], the Fourier transform of the quadratic phase is relative smooth and meanwhile band-limited when suitable parameters are set. Derived from the perspective of geometric optics, the field $\exp (i{\alpha _1}{x^2})$ with size Δx and the field $\exp (i{\alpha _2}{f_x}^2)$ with size Δfx are approximate Fourier transform pair if their parameters satisfy

$${\alpha _1} = \pi \frac{{\Delta {f_x}}}{{\Delta x}}\;\;\;\;{\alpha _2} = \pi \frac{{\Delta x}}{{\Delta {f_x}}}\;. $$

Supposing fp is the target frequency for the lattice pattern, and FOV is circular with radius of r0, then the size Δfx is set as fcofp to make those sperate Fourier spectrums fill the whole DOE region while not overlap with each other as far as possible. According to Eq. (4), the initial phase can be obtained as

$${\varphi _0}(x,y) = \exp [i\mathrm{\pi }\frac{{{f_{co}} - {f_p}}}{{2{r_0}}}({x^2} + {y^2})]. $$

The iterative process based on GS algorithm is depicted in Fig. 1. Note that some modifications are applied and will be explained in detail. The first one is the initial phase and it is clarified below. Next one is that the sampling rate of amplitude in the sample plane is reduced for describing target pattern here more precisely. The last modification is the constraint in the sample plane, compared with the direct replacement with the target amplitude, the amplitude constraint is defined as

$${a_p}^{(k + 1)} = d({\mathbf r}){a_p}^{(k)}{(\gamma {a_{pt}}({\mathbf r})/{a_p}^{\prime}({\mathbf r}))^{{\beta _k}}} + (1 - d({\mathbf r})){a_p}^{\prime}({\mathbf r})$$

 figure: Fig. 1.

Fig. 1. Flowchart of our iterative algorithm for designing phase-only DOE.

Download Full Size | PDF

It means the sample plane is separated into signal and noise windows by binary mask function d(r). γ is a constant between 0 to 1 for controlling the signal energy ratio. And ${\beta _{k + 1}} = \sqrt {{\beta _k}} \;({\beta _1} = 1{\textrm{0}^{\textrm{ - 9}}})$ is used for obtaining next iteration adaptive weight to accelerate the convergence speed.

With this method, we have generated both square and hexagonal lattices. The simulation result of square lattice and the corresponding spectrum are shown in Fig. 2. In the simulation, we use the same wavelength and effective numerical aperture as in the experiment, and the phase in the DOE plane is quantized with eight-step. Square lattice has nine orders, and the highest frequency in the diagonal direction we can design is very close to the cutoff frequency while in previous work this frequency is not higher than 0.58fco plotted by solid line [21].

 figure: Fig. 2.

Fig. 2. Simulation result and the spectrum of the square lattice. (a) Simulation of the square lattice and scale bar 1 µm, (b) Frequency orders in Fourier spectrum, the cut-off frequency is marked with dash line, the obtained highest frequency and the frequency realized by previous work are marked with dotted line and solid line, respectively.

Download Full Size | PDF

Theoretically, for SIM, the hexagonal lattice can provide an isotropic resolution and a better signal-to-noise ratio than the square one [11,32]. Therefore, we choose hexagonal lattice as the illumination pattern in our SIM system. To evaluate this method, we design many DOEs to generate hexagonal lattice patterns with different spatial frequencies.

The diffraction efficiency is defined as the ratio of the signal energy within the target area over the total energy. The lost energy is shifted to the region away from the target area with our method. And the modulation depth of illumination pattern to evaluate the uniformity, is defined as the energy ratio of the highest frequency order over zero frequency order. High modulation depth assures high-quality construction for SIM, otherwise it will cause severe artifacts due to low signal-to-noise ratio raw images [33]. Therefore, we design our DOE to generate patterns while keeping the modulation depth over 80%. As shown in Fig. 3, the diffraction efficiency will remarkably decline with the increasing frequency when the ratio fp over fco (fp/fco) is larger than 0.5.

 figure: Fig. 3.

Fig. 3. Simulation results of hexagonal lattices with different spatial frequencies.

Download Full Size | PDF

Although we have fabricated the DOE whose fp/fco is set as 0.9, the efficiency is low and the super-resolution construction is of poor quality because of high zero order intensity and inferior modulation in the experimental pattern resulting from the fabrication error. We have to use the DOE with the spatial frequency of 0.75fco to carry out the experiments, and Fig. 4 demonstrates the simulation and experimental patterns realized by such a DOE. Despite the intensity deviation, the experimental result is consistent with the simulation, especially from the detail comparison in the magnification view. Another convincing point lies that they have the same spatial frequency orders. Considering the NA, the cut-off frequency fco = 5.08 µm-1, and the highest frequency of pattern fp = 3.81 µm-1, which is consistent with the design target. The experimental efficiency is 41.4% by calculating the ratio of the energy in lattice area over the total incidence energy.

 figure: Fig. 4.

Fig. 4. Verification for DOE design. (a) Simulation result of hexagonal pattern and enlarged ROI in box. (b) experiment results of hexagonal pattern and enlarged region of interest (ROI) in box; Scale bar 1µm. (c, d) Fourier spectrums of (a) and (b) show all the frequency orders, and the cut-off frequency is marked with white dash line.

Download Full Size | PDF

But the experimental pattern is deviated from the ideal one if we use a phase-determined coefficient matrix to separate different spectrum components similar in linear SIM, there will be artifact in the reconstruction image. Therefore, we use an iterative algorithm, termed as pattern-illuminated Fourier ptychography (pi-FP) detailed in Ref. [34], to avoid the artifact and robustly construct the super-resolution image.

3. SIM system with DOE

For the purpose of implementing the lattice as the illumination pattern, a beam modulation experimental set is established as shown in Fig. 5, and the DOE consists of 1250 × 1250 square pixels of 4µm size. In the illumination system the 488 nm laser beam coupled into a single mode fiber is collimated with a lens L1. After being modulated by the DOE, which is mounted on the filter cube (B-2A fluorescence filter cube, Nikon), the laser beam reflects from the dichroic mirror, passes through the objective (100×/1.45, Nikon) and forms a hexagonal lattice to illuminate the sample plane. The illumination system is highly compact because it generates the illumination pattern just through a single ∼1 mm thickness DOE mounted on the filter cube.

 figure: Fig. 5.

Fig. 5. The experimental setup. Illumination and fluorescence are shown in blue and green, respectively. (a) the phase in the DOE plane. (b) Simulation of hexagonal pattern in the sample plane. (c) DOE mounted on the filter cube.

Download Full Size | PDF

Epi-fluorescence illumination is used here for further simplifying the whole system. The image system shares the same objective (100×/1.45, Nikon) with the illumination system. An emission filter is used for filtering out the excitation light (Sapphire 488, Coherent), and the fluorescence raw images are captured by an sCMOS camera (edge 4.2, pco.). Furthermore, the phase shift unit is established specially to provide a better support for oil immersion objective. The fiber port is stick to a motorized stage (C-483, PI), scanning through 25 (5 × 5) different positions to provide phase shifts. Because the objective and collimation lens L1 form a telescope system, the displacement of fiber port can be transferred with a certain scale to the sample plane. The step size of phase shift in the sample plane is 1/5fp in x direction and 2/(${5}\sqrt {3} $fp) in y direction. Considering that the highest spatial frequency of hexagonal lattice used in the experiment is 3.81µm-1 and the magnification of the telescope system is 20×, the scanning step size is 2.10µm in x direction and 2.42µm in y direction.

4. Experimental results

4.1 Quantitative resolution measurement with fluorescence beads

The resolution performance of this system is testified quantitatively by imaging a cluster of 100 nm fluorescence microspheres (F36924, ThermoFisher). Comparing to the diffraction-limited images, the resolution enhancement is clearly visible with the SIM system shown in Fig. 6. In this paper, the widefield image is the average of all scanned raw images. In addition, those raw images could be used to estimate the modulation depth. We compute the modulation depth per pixel as shown in Fig. S3. Then the modulation depth of 20 pixels, mostly the center of the beads, are averaged to estimate the modulation depth of the lattice pattern as 0.712.

 figure: Fig. 6.

Fig. 6. Imaging results of fluorescence beads. (a) Widefield image of fluorescence beads and enlarged ROI (lower right) in the boxed region, and Fourier spectrum (upper right). (b) SIM image of fluorescence beads and enlarged ROI (lower right) in the boxed region, and Fourier spectrum (upper right). (c) Intensity profiles of fluorescence beads on the selected line in (a) and (b). (d) Profiles of PSFs in widefield microscope and SIM system. Scale bar 1 µm.

Download Full Size | PDF

As shown in the enlarged images that two adjacent fluorescence beads with a central distance of 177 nm can be resolved in the SIM system but the situation is not in widefield microscope. Intensity distributions of 20 beads are fitted to Gaussian model and averaged to approximately represent the point spread function (PSF) of the image system, and the full-width-at-half-maximum (FWHM) value of the function also backs up the resolution enhancement. Here the average size of the FWHM in the SIM system amounts 131 nm, while 240 nm in widefield microscope, meeting the diffraction limit considering the fluorescence label (Ex488/Em519) and objective (100×/1.45 NA, Nikon) used in this experiment. In addition, the Fourier spectrum of the construction image is expanded in all direction compared to that of the widefield image, which verifies the isotropic resolution improvement with hexagonal lattice pattern.

4.2 Super-resolution imaging with fixed cells

To demonstrate the SIM on biological sample we use a FluoCell Prepared Slide #1 (F36924, Thermofisher), which contains bovine pulmonary artery endothelial (BPAE) cells stained with green-fluorescent Alexa Fluor 488 phalloidin. Figure 7 presents the obtained diffraction-limited image and the SIM reconstruction. The intensity profiles on the yellow selected lines present that the SIM visualizes the separation of two neighbor filaments, which are not resolved in the diffraction-limited image. The Fourier spectrum of the construction image is expanded compared to that of the widefield image. Furthermore, we evaluate the results through analysis of circular averaged power spectral density (PSDca) [35]. This method estimates the resolvable frequency at which the spectral power density decreases to noise level, and it also supports the resolution enhancement of the SIM system, from 230 nm to 137 nm.

 figure: Fig. 7.

Fig. 7. Imaging results for tubulin of BPAE cells. (a) diffraction-limited image of tubulins, enlarged ROI (lower left), and Fourier spectrum (upper right). (b) SIM image of tubulins, enlarged ROI (lower left), and Fourier spectrum (upper right). (c) Intensity profiles on the yellow lines in (a) and (b). (d) Circular average power spectral density of image in (a) and (b). Scale bar 1 µm.

Download Full Size | PDF

5. Discussion and conclusion

In conclusion, we use a well-designed DOE to achieve a lattice pattern with high spatial frequency approaching to the cut-off spatial frequency for illumination. The SIM system, mounting the DOE on the filter cube, avoids extra room compared to an inverted fluorescence microscope. The experiment lateral resolution is quantitatively testified by measuring FWHM of fluorescence beads and circular average power spectral density of BPAE cells. This advantage of miniaturization may help this SIM system more suitable for those imaging applications where both resolution and portability are needed.

We have experimentally demonstrated the pattern with spatial frequency of 0.75fco, and the simulation have shown spatial frequency of 0.9fco, but this high spatial frequency DOE is not applied here because of the low efficiency and poor pattern quality resulting from the fabrication error of the DOE. It may be improved by optimizing the super-resolution construction algorithm to alleviate the influence of illumination pattern. Furthermore, an ultra-high spatial frequency over the cut-off frequency still deserves to be explored for more than two-fold resolution enhancement. It suggested that multi-lobe superoscillating pattern with spatial frequency over the cut-off frequency is potential to be used in SIM [36]. Theoretically, this multi-lobe superoscillating pattern may be generated with DOE, but the key is to separate the superoscillating area with those strong lobes and increase the number of the oscillations, which will be further explored.

Funding

National Natural Science Foundation of China (Grant No. 62075112).

Acknowledgments

We acknowledge Assistant Professor Lingfeng Pan of Institute of Semiconductors, Chinese Academy of Sciences, for fabricating the designed DOEs.

Disclosures

The authors declare that there are no conflicts of interest related to this paper.

Data Availability

Raw data presented in this paper and code for analyzing the data are available in Code 1, Ref. [37].

Supplemental document

See Supplement 1 for supporting content.

References

1. J. W. Goodman, Introduction to Fourier Optics, 4th edition (Macmillan Learning, 2017).

2. E. Betzig, G. H. Patterson, R. Sougrat, O. W. Lindwasser, S. Olenych, J. S. Bonifacino, M. W. Davidson, J. Lippincott-Schwartz, and H. F. Hess, “Imaging intracellular fluorescent proteins at nanometer resolution,” Science 313(5793), 1642–1645 (2006). [CrossRef]  

3. M. J. Rust, M. Bates, and X. Zhuang, “Sub-diffraction-limit imaging by stochastic optical reconstruction microscopy (STORM),” Nat. Methods 3(10), 793–796 (2006). [CrossRef]  

4. S. W. Hell and J. Wichmann, “Breaking the diffraction resolution limit by stimulated emission: stimulated-emission-depletion fluorescence microscopy,” Opt. Lett. 19(11), 780–782 (1994). [CrossRef]  

5. M. G. L. Gustafsson, “Surpassing the lateral resolution limit by a factor of two using structured illumination microscopy,” J. Microsc. 198(2), 82–87 (2000). [CrossRef]  

6. A. Markwirth, M. Lachetta, V. Mönkemöller, R. Heintzmann, W. Hübner, T. Huser, and M. Müller, “Video-rate multi-color structured illumination microscopy with simultaneous real-time reconstruction,” Nat. Commun. 10(1), 4315 (2019). [CrossRef]  

7. R. Lin, E. T. Kipreos, J. Zhu, C. H. Khang, and P. Kner, “Subcellular three-dimensional imaging deep through multicellular thick samples by structured illumination microscopy and adaptive optics,” Nat. Commun. 12(1), 3148 (2021). [CrossRef]  

8. D. Li, L. Shao, B.-C. Chen, X. Zhang, M. Zhang, B. Moses, D. E. Milkie, J. R. Beach, J. A. Hammer, M. Pasham, T. Kirchhausen, M. A. Baird, M. W. Davidson, P. Xu, and E. Betzig, “Extended-resolution structured illumination imaging of endocytic and cytoskeletal dynamics,” Science 349(6251), aab3500 (2015). [CrossRef]  

9. Y. Guo, D. Li, S. Zhang, Y. Yang, J.-J. Liu, X. Wang, C. Liu, D. E. Milkie, R. P. Moore, U. S. Tulu, D. P. Kiehart, J. Hu, J. Lippincott-Schwartz, E. Betzig, and D. Li, “Visualizing intracellular organelle and cytoskeletal interactions at nanoscale resolution on millisecond timescales,” Cell 175(5), 1430–1442.e17 (2018). [CrossRef]  

10. A. Sandmeyer, M. Lachetta, H. Sandmeyer, W. Hübner, T. Huser, and M. Müller, “Cost-Effective Live Cell Structured Illumination Microscopy with Video-Rate Imaging,” ACS Photonics 8(6), 1639–1648 (2021). [CrossRef]  

11. M. Schropp and R. Uhl, “Two-dimensional structured illumination microscopy,” J. Microsc. 256(1), 23–36 (2014). [CrossRef]  

12. R. Heintzmann, “Saturated patterned excitation microscopy with two-dimensional excitation patterns,” Micron 34(6-7), 283–291 (2003). [CrossRef]  

13. C.-W. Chen, P.-H. Wang, L.-J. Chou, Y.-Y. Lee, B.-J. Chang, and S.-Y. Chiang, “High-resolution light-scattering imaging with two-dimensional hexagonal illumination patterns: system implementation and image reconstruction formulations,” Opt. Express 25(18), 21652 (2017). [CrossRef]  

14. “ZEISS Elyra 7 with Lattice SIM Introduced,” https://www.labbulletin.com/articles/ZEISS-Elyra-7-Lattice-SIM-introduced.

15. E. Betzig, “Excitation strategies for optical lattice microscopy,” Opt. Express 13(8), 3021 (2005). [CrossRef]  

16. J.-Z. Wang, H. Pang, M. Zhang, L.-F. Shi, A.-X. Cao, Q.-L. Deng, and S. Hu, “Diffractive element design for generating multi-channel structured light field,” in 8th International Symposium on Advanced Optical Manufacturing and Testing Technologies: Design, Manufacturing, and Testing of Micro- and Nano-Optical Devices and Systems; and Smart Structures and Materials (SPIE, 2016), Vol. 9685, pp. 156–161.

17. I. Moreno, J. A. Davis, D. M. Cottrell, N. Zhang, and X.-C. Yuan, “Encoding generalized phase functions on Dammann gratings,” Opt. Lett. 35(10), 1536 (2010). [CrossRef]  

18. C. Wu, H. Gu, Z. Zhou, and Q. Tan, “Design of diffractive optical elements for subdiffraction spot arrays with high light efficiency,” Appl. Opt. 56(31), 8816–8821 (2017). [CrossRef]  

19. B.-C. Chen, W. R. Legant, K. Wang, L. Shao, D. E. Milkie, M. W. Davidson, C. Janetopoulos, X. S. Wu, J. A. Hammer, Z. Liu, B. P. English, Y. Mimori-Kiyosue, D. P. Romero, A. T. Ritter, J. Lippincott-Schwartz, L. Fritz-Laylin, R. D. Mullins, D. M. Mitchell, J. N. Bembenek, A.-C. Reymann, R. Böhme, S. W. Grill, J. T. Wang, G. Seydoux, U. S. Tulu, D. P. Kiehart, and E. Betzig, “Lattice light-sheet microscopy: Imaging molecules to embryos at high spatiotemporal resolution,” Science 346(6208), 1257998 (2014). [CrossRef]  

20. L. Zhu, J. Yu, D. Zhang, M. Sun, and J. Chen, “Multifocal spot array generated by fractional Talbot effect phase-only modulation,” Opt. Express 22(8), 9798–9808 (2014). [CrossRef]  

21. G. Liu, N. Xu, H. Yang, Q. Tan, and G. Jin, “Miniaturized structured illumination microscopy with diffractive optics,” Photonics Res. 10(5), 1317–1324 (2022). [CrossRef]  

22. N. Xu, G. Liu, and Q. Tan, “Adjustable super-resolution microscopy with diffractive spot array illumination,” Appl. Phys. Lett. 116(25), 254103 (2020). [CrossRef]  

23. Y. Ogura, M. Aino, and J. Tanida, “Design and demonstration of fan-out elements generating an array of subdiffraction spots,” Opt. Express 22(21), 25196 (2014). [CrossRef]  

24. L. C. Thomson, Y. Boissel, G. Whyte, E. Yao, and J. Courtial, “Simulation of superresolution holography for optical tweezers,” New J. Phys. 10(2), 023015 (2008). [CrossRef]  

25. “N-SIM E,” https://www.microscope.healthcare.nikon.com/products/super-resolution-microscopes/n-sim-e.

26. W. Qu, H. Gu, Q. Tan, and G. Jin, “Precise design of two-dimensional diffractive optical elements for beam shaping,” Appl. Opt. 54(21), 6521–6525 (2015). [CrossRef]  

27. L. Chen, S. Tian, H. Zhang, L. Cao, and G. Jin, “Phase hologram optimization with bandwidth constraint strategy for speckle-free optical reconstruction,” Opt. Express 29(8), 11645 (2021). [CrossRef]  

28. F. Wyrowski and O. Bryngdahl, “Iterative Fourier-transform algorithm applied to computer holography,” J. Opt. Soc. Am. A 5(7), 1058 (1988). [CrossRef]  

29. S. Zhang, G. Neil, and M. Shinn, “Single-element laser beam shaper for uniform flat-top profiles,” Opt. Express 11(16), 1942–1948 (2003). [CrossRef]  

30. D. Wang, “Adaptive-weight iterative algorithm for flat-top laser beam shaping,” Opt. Eng. 51(7), 074301 (2012). [CrossRef]  

31. H. Pang, J. Wang, M. Zhang, A. Cao, L. Shi, and Q. Deng, “Non-iterative phase-only Fourier hologram generation with high image quality,” Opt. Express 25(13), 14323 (2017). [CrossRef]  

32. M. Schropp, C. Seebacher, and R. Uhl, “XL-SIM: Extending superresolution into deeper layers,” Photonics 4(4), 33 (2017). [CrossRef]  

33. M. Li, Y. Li, W. Liu, A. Lal, S. Jiang, D. Jin, H. Yang, S. Wang, K. Zhanghao, and P. Xi, “Structured illumination microscopy using digital micro-mirror device and coherent light source,” Appl. Phys. Lett. 116(23), 233702 (2020). [CrossRef]  

34. S. Dong, P. Nanda, R. Shiradkar, K. Guo, and G. Zheng, “High-resolution fluorescence imaging via pattern-illuminated Fourier ptychography,” Opt. Express 22(17), 20856–20870 (2014). [CrossRef]  

35. K. Fliegel, M. Klíma, and J. Pospíšil, “Assessing resolution in live cell structured illumination microscopy,” in Photonics, Devices, and Systems VII (SPIE, 2017), p. 39.

36. N. Shapira, Z. Deng, R. Remez, D. Singh, E. Katzav, and A. Arie, “Multi-lobe superoscillation and its application to structured illumination microscopy,” Opt. Express 27(24), 34530–34541 (2019). [CrossRef]  

37. C. Zhang, “Reconstruction code,” figshare (2022) https://doi.org/10.6084/m9.figshare.21201001.

Supplementary Material (2)

NameDescription
Code 1       construction code used in the paper
Supplement 1       supplemental document

Data Availability

Raw data presented in this paper and code for analyzing the data are available in Code 1, Ref. [37].

37. C. Zhang, “Reconstruction code,” figshare (2022) https://doi.org/10.6084/m9.figshare.21201001.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. Flowchart of our iterative algorithm for designing phase-only DOE.
Fig. 2.
Fig. 2. Simulation result and the spectrum of the square lattice. (a) Simulation of the square lattice and scale bar 1 µm, (b) Frequency orders in Fourier spectrum, the cut-off frequency is marked with dash line, the obtained highest frequency and the frequency realized by previous work are marked with dotted line and solid line, respectively.
Fig. 3.
Fig. 3. Simulation results of hexagonal lattices with different spatial frequencies.
Fig. 4.
Fig. 4. Verification for DOE design. (a) Simulation result of hexagonal pattern and enlarged ROI in box. (b) experiment results of hexagonal pattern and enlarged region of interest (ROI) in box; Scale bar 1µm. (c, d) Fourier spectrums of (a) and (b) show all the frequency orders, and the cut-off frequency is marked with white dash line.
Fig. 5.
Fig. 5. The experimental setup. Illumination and fluorescence are shown in blue and green, respectively. (a) the phase in the DOE plane. (b) Simulation of hexagonal pattern in the sample plane. (c) DOE mounted on the filter cube.
Fig. 6.
Fig. 6. Imaging results of fluorescence beads. (a) Widefield image of fluorescence beads and enlarged ROI (lower right) in the boxed region, and Fourier spectrum (upper right). (b) SIM image of fluorescence beads and enlarged ROI (lower right) in the boxed region, and Fourier spectrum (upper right). (c) Intensity profiles of fluorescence beads on the selected line in (a) and (b). (d) Profiles of PSFs in widefield microscope and SIM system. Scale bar 1 µm.
Fig. 7.
Fig. 7. Imaging results for tubulin of BPAE cells. (a) diffraction-limited image of tubulins, enlarged ROI (lower left), and Fourier spectrum (upper right). (b) SIM image of tubulins, enlarged ROI (lower left), and Fourier spectrum (upper right). (c) Intensity profiles on the yellow lines in (a) and (b). (d) Circular average power spectral density of image in (a) and (b). Scale bar 1 µm.

Equations (6)

Equations on this page are rendered with MathJax. Learn more.

a p t ( r ) = n = 1 N cos ( 2 π f p t n r α 0 ( n ) ) ,
E ( f ) = F { e p ( r ) } = F { a p ( r ) × d ( r ) exp [ i φ ( r ) ] } = F { a p ( r ) } F { d ( r ) exp [ i φ ( r ) ] } = P ( f ) D ( f )
E ( f ) = n = 1 N [ D ( f f p t n ) + D ( f + f p t n ) ] .
α 1 = π Δ f x Δ x α 2 = π Δ x Δ f x .
φ 0 ( x , y ) = exp [ i π f c o f p 2 r 0 ( x 2 + y 2 ) ] .
a p ( k + 1 ) = d ( r ) a p ( k ) ( γ a p t ( r ) / a p ( r ) ) β k + ( 1 d ( r ) ) a p ( r )
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.