## Abstract

We present a theoretical investigation of an optical microscope design that achieves wide-field, multiphoton fluorescence microscopy with finer axial resolution than confocal microscopy. Our technique creates a thin plane of excitation light at the sample using height-staggered microlens arrays (HSMAs), wherein the height staggering of microlenses generate temporal focusing to suppress out-of-focus excitation, and the dense spacing of microlenses enables the implementation of structured illumination technique to eliminate residual out-of-focus signal. We use physical optics-based numerical simulations to demonstrate that our proposed technique can achieve diffraction-limited three-dimensional imaging through a simple optical design.

© 2013 OSA

## 1. Introduction

Fluorescence is one of the most important optical contrasts in biological imaging, for it offers high molecular sensitivity and specificity. The general desire for obtaining three-dimensional (3D) information of dynamic biological systems has spurred the development of several fluorescence microscopies [1–4] with the ability to selectively image thin sections of an intact sample. This ability of optical sectioning can be categorized into two general regimes - intrinsic and extrinsic. Intrinsic optical sectioning utilizes an axially peaked excitation profile along the imaging path to obtain optical sectioning; examples include two-photon excitation microscopy (2PE) [2] and selective plane illumination microscopy (SPIM) [3]. Extrinsic optical sectioning is performed post-excitation and primarily relies on a far-field imaging principle: the higher spatial-frequency components of the images decay more rapidly with defocusing [5]. As a result, one can create a high-spatial-frequency excitation pattern and use appropriate spatial filters to selectively extract the in-focus signal. Confocal microscopy (CFM) [5] and structured illumination microscopy (SIM) are two examples that use physical and numerical spatial filters, respectively, to achieve equivalent optical sectioning [6].

Present imaging techniques, however, have important drawbacks. For CFM and 2PE, phototoxicity and photothermal damages are noted [1, 7], and are exacerbated if video-rate acquisition is required. To reduce these effects, spatially/temporally parallelized CFM/2PE systems, such as spinning disk confocal [8], multifocal multiphoton microscopy (MMM) [9], and time-multiplexed multifocal multiphoton microscopy (TM-MMM) [10] were developed. In these techniques, multiple foci are created as independent, parallel channels for excitation in and detection from the sample. To reduce crosstalk among parallel channels, however, the foci spacing has to be sufficiently sparse (typically <1 foci per 10^{2}-focus area, see Appendix A, B), which limits the degree of parallelization. SIM and SPIM do not suffer from parallelization issues; however both techniques have their own limitations. SIM can lead to significant photobleaching and low signal-to-noise ratio (SNR) in the processed images, because it excites fluorophores and receives fluorescence over a wide axial range at each acquisition. SPIM requires a compromise between axial resolution and the size of the field of view. Moreover, SPIM raises design issues and challenges in sample handling and manipulation for its separate excitation and imaging objectives in close proximity.

In principle, intrinsic and extrinsic optical sectioning can be combined to take the advantages of both. For example, it has been shown that a SPIM system can incorporate SIM to improve its optical sectioning when a reasonably large field of view is desired [11]. Here, we propose a hybrid optical microscope design (Fig. 1) incorporating SIM and temporal focusing [4], a wide-field, intrinsic sectioning technique. Our design utilizes height-staggered microlens arrays (HS-MAs) having periodic patterns (Fig. 1(b)) to precisely control the time staggering for temporal focusing, and to provide an excitation pattern for SIM to further suppress residual out-of-focus signal. We herein refer to this approach as Microlens Array and Structured Illumination-based Plane-Projection Multiphoton Microscopy (MASI-PPMP).

## 2. MASI-PPMP: integrating (intrinsic) temporal focusing and (extrinsic) SIM

We have previously demonstrated a new technique, plane-projection multiphoton microscopy (PPMP) [12], in which temporal focusing with randomized excitation patterns is achieved through the randomly staggered light pulses (in the time domain) created by the surface roughness of a ground-glass optical diffuser. The randomized excitation patterns, as suggested by a previous report [13], could be used for SIM. However, mitigating the excitation patterns in the processed images requires considerable effort. For simplicity, periodic excitation patterns are commonly used [6], and chosen for the present study. Further, the design control afforded in the HSMA allows for optimization of the intrinsic sectioning.

The combination of the microlens arrays and temporal masks, as used in TM-MMM [14], can provide equivalent functions of the proposed HSMAs, although additional precision alignment and positioning are required. Notably, the step sizes of time staggering Δ*t* in TM-MMM and MASI-PPMP are different. In TM-MMM, Δ*t* is set to be much larger than the pulse duration *τ*_{0} to reduce temporal inter-foci interferences [14], by which individual foci can be considered as independent, parallel channels. In contrast, MASI-PPMP allows Δ*t* < *τ*_{0} for practical reasons detailed below and in Supplementary Materials.

We begin the design of MASI-PPMP with consideration of the physical limitations of HS-MAs. The dimensions of HSMAs are limited by fabrication technologies as well as diffractive losses. Existing techniques cannot straightforwardly engineer micro-optics of depth variation >1 mm [15]. In addition, when we consider the HSMA as an array of time-delay channels, a light pulse propagating in one channel can leak into its neighboring channels due to the nature of diffraction. A simple estimation of the appropriate dimensions of HSMAs is given in Appendix A. In short, to have negligible inter-channel light leakage, the spatial parameters of the HSMAs should satisfy

where*d*is the aperture diameter of a single time-delay channel,

*λ*

_{0}is the central wavelength of the excitation light (here ∼ 800 nm), and Δ

*h*is the height staggering between the longest and shortest time-delay channels. Accordingly, we restrict Δ

_{max}*h*to ∼ 300

_{max}*μ*m and 36

*μ*m ≥

*d*≥ 18

*μ*m, via Eq. (1). These dimensions make the HSMA fabrication feasible through existing techniques, and can create up to ∼ 10

^{6}foci within a 2-inch aperture, thereby providing the periodic, high-spatial-frequency pattern for SIM. For a simple analysis, we assume that the total amount of time delay Δ

*t*is separated equally into

_{tot}*N*distinct time-delay steps (with step size Δ

_{t}*t*, Eq. (2)). These time-delay steps are then arranged in a prescribed periodic pattern in the HSMA (Fig. 1(b)). Considering the propagation speed of light in a material, we have

*c*and

*n*are the speed of light in vacuum and the refractive index of the material of the HSMAs (set to be 1.5), respectively.

Because of the limitation of Δ*h _{max}*, the estimated Δ

*t*in Eq. (2) can be around or shorter than the pulse duration of conventional ultrafast oscillators or amplifiers above certain values of

*N*. In such cases, we should take into account the temporal interferences among light pulses of different time delays. Notably, temporal masks with much larger Δ

_{t}*t*have been proposed to avoid temporal interferences and to achieve scanningless TM-MMM [14]. However, Eqs. (1) and (2) indicate that the required Δ

_{tot}*h*falls far beyond the limits of existing fabrication techniques, and the aperture sizes of these temporal masks could be too large for standard biomedical microscopes (see Appendix B).

_{max}To have fabrication-feasible HSMAs and to improve the SNR of the SIM post-processed images, we turn to optimizing the parameters of HSMAs. In the following sections, we first construct a model that considers the temporal interferences (§3); through this model we investigate how intrinsic optical sectioning depends on the choice of *N _{t}* (and Δ

*t*) for a given Δ

*t*(§4), and the arrangement of time delays in the HSMA (Appendix C). Then, we couple HSMAs with SIM to enhance optical sectioning (§5).

_{tot}## 3. Construct a physical optics-based model taking into account temporal interferences

Given the temporal focusing effect in the proposed technique, and with the non-negligible temporal interferences among the light pulses, a time-independent model such as that previously used for TM-MMM [14] is no longer sufficient for analyzing the performances of our setup. Thus, we develop a new model taking into account the time-dependent optical phase of multiple spectral components in an ultrafast pulse. For simplicity, we use a Gaussian-pulse approach, i.e., in the excited area, the electric field *E*(**r**, *z*, *t*) at position (**r**, *z*) (*z* = 0 at the specimen plane) and time *t* is approximated as the Gaussian-weighted sum of a series of constant-interval (in *k*–space), in-phase light waves,

*k*

_{0}is the central wavenumber of the pulse spectrum, and

*E*

_{kj}is the scalar field of the light wave of wavenumber

*k*. To approximate the ultrafast pulse train generated by the amplified laser system we used previously [12], we set

_{j}*k*

_{0}≈ 7.85 × 10

^{4}cm

^{−1}and a pulse duration

*τ*

_{0}of ≈ 30 fs (by using an appropriate

*σ*). We then employ the amplitude point spread function (PSF) derived previously for high numerical-aperture (NA) lenses [14, 16],

_{k}*E*

_{kj}, as

*α*is the maximal focusing angle

*θ*of the objective lens, and

*J*

_{0}is the zero-order Bessel function of the first kind. The objective lens used in all the simulations presented here is a 60X oil-immersion lens of NA 1.42.

Through Eqs. (3) and (4) we can numerically evaluate the time-dependent amplitude PSF of an ultrafast pulse focused by a well-corrected objective lens, *E _{PSF}* (

**r**,

*z*,

*t*). Having solved

*E*numerically, we estimate the electric field near the specimen plane,

_{PSF}*E*(

_{SP}**r**,

*z*,

*t*), as the linear superposition of the

*E*from the individual microlenses,

_{PSF}**r**and Δ

_{m}*t*

**are the central position and time delay of the ultrafast pulse going though the**

_{m}**m**-th HSMA channel, respectively. An example of numerically simulated

*E*(

_{SP}**r**,

*z*,

*t*) is shown in Media 1, which also reveals the dynamic process of temporal focusing (see Appendix D). Notably, to fulfill the wide-field illumination condition and to simplify the simulations, our model assumes that a ’unit’ microlens array is infinitely replicated in the transverse coordinates, as shown in Fig. 1(b). Under such a periodic condition, the physical optics properties in the projected region of one unit microlens array is sufficient to represent the entire system (the validity and details of this model are given in Appendix E).

Through Eq. (5), the excitation intensity profiles *I _{SP}*(

**r**,

*z*) can be derived by integrating the excitation intensity over time, as

*n*is the number of photons required in single excitation event (here

_{p}*n*= 2).

_{p}## 4. Optimize intrinsic optical sectioning through tuning *N*_{t} and Δ*t*

_{t}

To quantify the intrinsic optical sectioning created by a particular design of HSMA, we calculate the total fluorescence signal at a depth *z*, *S*(*z*), through integration of *I _{SP}*(

**r**,

*z*) in Eq. (6) over the transverse coordinates, as

*S*(

*z*) corresponds to the detected fluorescence signal from a thin fluorescent film placed at a depth

*z*.

To determine the choice of (*N _{t}*, Δ

*t*) that produces the most efficient optical sectioning, we evaluate

*S*(

*z*) and the ratio of out-of-focus and in-focus signal,

*S*/

_{out}*S*, for various sets of

_{in}*N*and Δ

_{t}*t*(constrained by Eq. (2)) with two inter-foci spacings,

*d*(≡

_{foci}*d*/

*M*, where

*M*is the magnification of the microscopy system) ≈ 0.4 and 0.8

*λ*

_{0}, corresponding to ∼ 56 and ∼ 14 foci per 10

^{2}-focus area, respectively (Fig. 2). Because of the square geometry of our HSMAs, we examine

*N*= 2

_{t}^{2}, 3

^{2}, 4

^{2}, ... and 9

^{2}. Analyzing

*S*/

_{out}*S*reveals that the decay of the out-of-focus excitation significantly slows down between

_{in}*N*= 25 and 64, corresponding to Δ

_{t}*t*≈ 20-8 fs (i.e., 2/3

*τ*

_{0}-1/4

*τ*

_{0}) (Fig. 2(c)). We observe similar trends of

*S*/

_{out}*S*for various geometrical arrangements of the distinct time-delay steps (see Appendix C). This result suggests that the intrinsic optical sectioning for a fixed Δ

_{in}*t*is optimized when Δ

_{tot}*t*is slightly smaller than

*τ*

_{0}. Further reducing Δ

*t*(equivalent to increasing

*N*(Eq. (2))) can complicate the fabrication of the HSMAs without major improvement of the intrinsic optical sectioning. In addition, increasing

_{t}*d*leads to weaker out-of-focus excitation and a less complex axial excitation profile

_{foci}*S*(

*z*) (Fig. 2(b)).

## 5. Obtain extrinsic optical sectioning via implementing SIM

Next, we implement SIM to extrinsically remove the residual out-of-focus signal of the optimized HSMA. SIM relies on the post processing of multiple images, each noted as *I _{IM}*. Because the emitted fluorescence (wavelength assumed to be ∼ 0.56

*λ*

_{0}) from the specimen is generally incoherent, we can estimate

*I*from a convolution of the excitation intensity profile

_{IM}*I*and the intensity PSF of the microscopy system,

_{SP}*I*[5], as

_{SYS}*f*is the concentration distribution of the fluorophore in the specimen. To quantify the optical-sectioning effect, we assume that the specimen is an ideal thin fluorescent film placed at

*z*=

*z*, i.e.,

_{f}*f*(r,

*z*) =

*δ*(

*z*−

*z*). Eq. (8) then becomes

_{f}*I*is 2-dimensional. We thus use 3-by-3 shifts (

_{SP}*I*

_{IM}_{1},

*I*

_{IM}_{2}, ...,

*I*

_{IM}_{9}denote the obtained images) and apply SIM post-processing to extract the optically sectioned images, We then substitute

*I*in Eq. (7) with

_{SP}*I*to evaluate the overall strength of optical sectioning (Fig. 3).

_{SIM}The results show that MASI-PPMP can achieve a better optical sectioning than CFM. Such results are reasonable because the fundamental mechanisms of intrinsic and extrinsic optical sectioning are independent from each other and can be combined. We further simulate the images of a virtual 3D object obtained by the different imaging techniques (Fig. 4) to illustrate the optical sectioning of MASI-PPMP. Compared with conventional epifluorescence microscopy (Fig. 4(b)), the images obtained before SIM post-processing (denoted as MA-PPMP) shows the ability of intrinsic optical sectioning. In the log-scale intensity plot (Fig. 4(c)), MASI-PPMP has the highest signal contrast between fluorescent and non-fluorescent areas, consistent with the results in Fig. 3. Moreover, the reconstructed 3D views show that MASI-PPMP successfully reproduces the details of the object (Fig. 4(d)).

We should note that the SIM process introduced here, i.e., the 9-frame imaging procedure and Eq. (10), is not the only way to remove out-of-focus signal from the obtained images *I _{IM}*. Other methods using a high-spatial-frequency illumination pattern to distinguish in-focus and out-of-focus signals, such as HiLo microscopy [17, 18], may also be applied to the proposed optical setup. These methods use different imaging procedures and post-processing algorithms, which may lead to different imaging properties such as speed and spatial resolutions, such that different methods may be most suitable for different imaging applications. For example, the HiLo method, whose algorithms are more complicated than Eq. 10, requires only 2 frames to retrieve in-focus signal, one frame each with uniform and non-uniform illumination; the reduced number of frames may therefore shorten the imaging time. However, such an advantage exists only when the frame rate of the camera is much lower than the repetition rate of the ultrafast pulse train. If these two rates fall in the same order of magnitude, one will need to increase the exposure time to shift the structural illumination pattern around the sample so as to mimic the effect of uniform illumination, and the overall imaging time will eventually be similar to SIM.

## 6. Conclusion

In conclusion, our numerical modeling shows that MASI-PPMP can efficiently provide intrinsic optical sectioning via temporal focusing effect generated by the HSMAs, and achieve an axial resolution finer than CFM via the implementation of extrinsic optical-sectioning technique. Our analysis provides the design guide of HSMAs with small time-delay increments; the dimensions of the proposed HSMAs are compatible with standard biomedical microscopes, and feasible for existing fabrication methods. In contrast to conventional SIM, the intrinsic sectioning of our method can reduce photobleaching and increase the SNR of the processed images.

Compared with parallelized CFM and 2PE, MASI-PPMP has a denser spacing of foci (> 10 foci per 10^{2}-focus area), and thus has greater potential in high-speed imaging. In addition to the finer axial resolution than CFM, the lateral resolution can be enhanced via utilization of the periodic excitation structures and Fourier analysis of the obtained images [19]. Moreover, MASI-PPMP can use the pulse train generated by ultrafast amplifiers, thereby yielding a sufficient signal level and moderate impact on the specimen at high acquisition rates, as demonstrated in our previous PPMP setup [12]. Taken together, we have shown that MASI-PPMP has the potential to achieve (better-than-conventional) diffraction-limited, high-frame-rate 3D imaging with a simple, wide-field optical design.

## Appendix A: Determining the dimensions of the time-delay channels

In this section we discuss the appropriate dimensions of the height-staggered microlens arrays (HSMAs) utilized in our technique. As mentioned in the main article, the light waves propagating along one particular time-delay channel can leak into the neighboring channels due to diffraction. In the microscope, this leakage appears as stray light, which induces multiphoton excitation at all depths and thus increases out-of-focus excitation.

To quantify such an effect, we use a simple model consisting of two adjacent microlenses with height staggering Δ*h*, as shown in Fig. 5. When the light of wavelength *λ*_{0} exits from channel 1 (Ch1) and propagates to distance *z*, diffraction theory predicts the beam size, *d _{B}*, as a function of

*λ*

_{0}and

*z*. With this regard, we consider the inter-channel leakage to be non-negligible if

*d*before propagating to

_{B}*z*= Δ

*h*is larger than the diameter of the microlens,

*d*. In other words, the inter-channel leakage is negligible if

*z*= Δ

*h*, so that the

*d*decreases upon exiting Ch1 and reaches its minimum,

_{B}*d*, at

_{f}*z*= Δ

*h*. Based on the diffraction-limited spot size of an ideal thin lens, we can estimate

*d*as where

_{f}*f*

_{#}is defined as

*f*/

*d*[20]. Combining Eqs. (11), (12) and the geometry, we have and the criterion of negligible inter-channel light leakage is thus To ensure the inter-channel light leakage is negligible regardless of the geometrical arrangements of distinct time-delay steps, we can substitute the largest height difference in the HSMA, Δ

*h*, for Δ

_{max}*h*in Eq. (14), and thus derive

The argument here also applies to the design of the combination of temporal masks and microlens arrays used in TM-MMM [14]. In such a case, the ends of the time-delay channels are flat surfaces, and one can use the single-slit diffraction to approximate the beam size with propagation. Such an approach yields a similar criterion.

## Appendix B: Size of the microlens arrays to achieve scanningless imaging using conventional time-multiplexing approach

In this section we estimate the size of the HSMAs for scanningless TM-MMM, using the criterion derived in the previous section. To achieve a complete, gap-less plane of illumination for a truly scanningless wide-field imaging requires the placement of adjacent foci no more than one focal spot diameter apart. However, to preserve the intrinsic optical sectioning of conventional multiphoton microscopy, the foci (of the same timing) have to be spatially separated by a certain distance, *d*_{0}, which is commonly much larger than the size of the focal spot of the objective lens, *d _{fs}*. As a result, there will be substantial un-illuminated area between the foci. For example, take an oil-immersion, NA 1.42 objective lens: the resultant

*d*

_{0}and

*d*is approximately 5

_{fs}*λ*[14] and 0.3

*λ*, which makes the fraction of un-illuminated area approximately $1-{\left(\frac{{d}_{fs}}{{d}_{0}}\right)}^{2}\approx 99.6\%$, equivalent to ∼0.4 foci per 10

^{2}-focus area. To cover the un-illuminated area, Egner et al. [14] proposed to use foci that are largely separated in time, so that the interferences among these foci is negligible even though they partially overlap in space. Consequently, the number of distinct time-delay steps required to cover the un-illuminated area,

*N*, can be estimated as

_{t}Given *N _{t}*, and using the propagation speed of light in the HSMA material, we can estimate the difference of height between the longest and shortest time channels, Δ

*h*, as:

_{max}*t*,

*c*and

*n*are the step size of the distinct time-delay steps, speed of light in vacuum and the refractive index of the material of the HSMAs, respectively. The appropriate values of Δ

*t*for negligible temporal interferences, as suggested in the previous study of TM-MMM [14], are equal to or larger than twice the pulse duration of the light pulse,

*τ*

_{0}, i.e., Δ

*t*≥ 2

*τ*

_{0}(where

*τ*

_{0}was set as ≈ 100 fs for conventional ultrafast oscillators) [14]. Accordingly, we find that Δ

*h*≈ 33 mm and

*d*≈ 0.15 mm using Eqs. (15)–(17), wherein

*λ*and

*n*are assumed to be ∼ 800 nm and 1.5, respectively. As a result, if one needs 1000-by-1000 foci in the field of view (FOV), the typical aperture of the entire HSMA will be as large as 1000×

*d*≈ 150 mm, which is considerably larger than the optical elements of a standard biomedical microscope. Furthermore, when an ultrafast pulse propagates in a long, narrow glass cylinder, both group velocity dispersion and modal dispersion broaden the pulse width, and thus reduce the axial resolution as well as excitation efficiency. From Eqs. (15), (17) and the requirement of Δ

*t*≥ 2

*τ*

_{0}, we can derive that Δ

*h*∝

*τ*

_{0}and $d\propto \sqrt{{\tau}_{0}}$, indicating that Δ

*h*and

*d*can be scaled down via using small

*τ*

_{0}. However, commonly available ultrafast systems typically achieve

*τ*

_{0}≥ 10 fs, such that the resultant Δ

*h*and

*d*remain considerably large. Moreover, the ultrafast pulses of smaller

*τ*

_{0}(and thus greater spectral bandwidth) will suffer more from group velocity dispersion and modal dispersion while propagating in the HSMA materials.

## Appendix C: Various geometrical arrangements of distinct time-delay steps show a similar trend of *S*_{out} / *S*_{in}

_{out}

_{in}

Figure 6 shows the simulation results of the ratio between out-of-focus and in-focus signal, *S _{out}* /

*S*, using various geometrical arrangements of distinct time-delay steps. As concluded in the main article, these

_{in}*S*/

_{out}*S*curves show that the decays of the out-of-focus excitation slow down significantly when

_{in}*N*= 25–64, (Δ

_{t}*t*≈ 2/3

*τ*

_{0}-1/4

*τ*

_{0}), although the curves from different patterns of time delays show slight differences in the values of

*S*/

_{out}*S*.

_{in}## Appendix D: Simulation of the dynamic *E*_{SP} in space reveals temporal focusing

_{SP}

In
Media 1 (a static frame shown in Fig. 7), we show an example of the numerically simulated *E _{SP}*. In particular, we plot the distributions of the real part of the electric field (Re(

*E*), upper panel) and intensity (|

_{SP}*E*|

_{SP}^{2}, lower panel). The geometrical arrangement of the distinct time-delay steps used for this simulation is a spiral pattern (Fig. 1(b) of main article) with

*N*= 49 and

_{t}*d*≈ 0.4

_{foci}*λ*

_{0}. Media 1 shows that at most of the area, except for the in-focus region (

*z*≈ 0), the intensity is moderate during nearly the entire time course (100–700 fs). On the other hand, at in-focus region the intensity is high in a relatively short period of time (450–540 fs). Such a phenomenon is commonly recognized as temporal focusing [4].

## Appendix E: Validity of the assumption of infinitely replicated unit HSMA in a wide-field microscope

In this section, we validate the model of infinitely replicated unit HSMAs for conventional wide-field biomedical microscopes. As described in the main article, the analysis of the properties of the proposed system can become simpler if we assume that the whole HSMA consists of infinitely, laterally repeated unit HSMAs. This is because the system is laterally periodic under such an assumption (referred to as the inf-HSMA model herein) and thus the physical properties (ex. the distribution of electric field) of one unit region can be used to describe the entire system. Although the infinite replication of unit HSMAs is not realistic for either experiments or numerical simulations, the inf-HSMA model is made tractable by assuming that, for the multiphoton excitation intensity *I _{SP}* at arbitrary positions, the contributions from very distant foci are negligible. In other words, one can approximate the results of using the infinitely replicated unit HSMAs via consideration of only the contributions from the foci within a certain range. To justify such an approach, we need to examine if the computed

*I*becomes stationary (i.e., converging to a certain value) as the range of foci considered to contribute, defined by

_{SP}*r*(Fig. 8(a)), is expanded.

_{tot}To proceed, we compute time-integrated intensity *I* using Eq. (6) of the main article, with only the contributions to *E _{SP}* from the foci within a distance

*r*considered (Fig. 8(a)). Here we use

_{tot}*I*, the intensity at a given position as derived with a relatively large

_{end}*r*(≈ 375

_{tot}*λ*

_{0}, where

*λ*

_{0}is the central wavelength of the light pulse), as a reference value, and plot

*I*/

*I*as

_{end}*r*increases from 0 to 125

_{tot}*λ*

_{0}. Fig. 8(b) shows that values of

*I*at all of the 100 randomly chosen positions converge to their corresponding

*I*as

_{end}*r*increases. In particular, we note that when

_{tot}*r*≥ 62.5

_{tot}*λ*

_{0}(equivalent to 50

*μ*m in physical dimensions), the error of

*I*, defined as |

*I*−

*I*|/

_{end}*I*, is less than 0.1%. Such results suggest that, for a given error tolerance in numerical simulations, we can use the inf-HSMA model by considering only the contributions of foci within a certain distance

_{end}*r*. Here, we have

_{inf}*r*≈ 50

_{inf}*μ*m for a 0.1% error tolerance. In the physical microscopy system, we can also use

*r*to determine the region wherein the inf-HSMA assumption is valid, as shown in Fig. 8(c). For conventional biomedical microscopes using

_{inf}*M*= 60X objective lenses, the diameter of the full FOV is typically lager than 300

*μ*m. Thus, the inf-HSMA model is valid in the central region of diameter larger than 200

*μ*m. At the image plane, this region corresponds to a disk of diameter ∼ 12 mm (200

*μ*m×

*M*) or larger, which is able to cover most conventional imaging devices.

## Acknowledgment

We gratefully acknowledge Dr. Alexander Egner of Laser-Laboratorium Gttingen e.V. for his kind advices on algorithm development. We are indebted to Mr. Yu-Hang Chen of National Taiwan University (NTU) for his assistance with the simulations, and we thank Dr. Chien-Cheng Chen of NTU for helpful comments. The support from Ellison Medical Foundation is acknowledged.

## References and links

**1. **J. B. Pawley, *Handbook of Biological Confocal Microscopy*, 3rd ed. (Springer, 2006). [CrossRef]

**2. **W. Denk, J. Strickler, and W. Webb, “2-photon laser scanning fluorescence microscopy,” Science **248**, 73–76 (1990). [CrossRef]

**3. **J. Huisken, J. Swoger, F. Del Bene, J. Wittbrodt, and E. Stelzer, “Optical sectioning deep inside live embryos by selective plane illumination microscopy,” Science **305**, 1007–1009 (2004). [CrossRef]

**4. **D. Oron, E. Tal, and Y. Silberberg, “Scanningless depth-resolved microscopy,” Opt. Express **13**, 1468–1476 (2005). [CrossRef]

**5. **T. Wilson, *Confocal Microscopy* (Academic Press, 1990).

**6. **M. Neil, R. Juskaitis, and T. Wilson, “Method of obtaining optical sectioning by using structured light in a conventional microscope,” Opt. Lett. **22**, 1905–1907 (1997). [CrossRef]

**7. **B. Masters, P. So, C. Buehler, N. Barry, J. Sutin, W. Mantulin, and E. Gratton, “Mitigating thermal mechanical damage potential during two-photon dermal imaging,” J. Biomed. Opt. **9**, 1265–1270 (2004). [CrossRef]

**8. **I. Akira, T. Takeo, I. Katsumi, S. Yumiko, K. Yasuhito, M. Kenta, A. Michio, and U. Isao, “High-speed confocal fluorescence microscopy using a nipkow scanner with microlenses for 3-d imaging of single fluorescent molecule in real time,” Bioimaging **4**, 57–62 (1996-06).

**9. **J. Bewersdorf, R. Pick, and S. Hell, “Multifocal multiphoton microscopy,” Opt. Lett. **23**, 655–657 (1998). [CrossRef]

**10. **V. Andresen, A. Egner, and S. Hell, “Time-multiplexed multifocal multiphoton microscope,” Opt. Lett. **26**, 75–77 (2001). [CrossRef]

**11. **P. J. Keller, A. D. Schmidt, A. Santella, K. Khairy, Z. Bao, J. Wittbrodt, and E. H. K. Stelzer, “Fast, high-contrast imaging of animal development with scanned light sheet-based structured-illumination microscopy,” Nat. Methods **7**, 637–U55 (2010). [CrossRef]

**12. **J.-Y. Yu, C.-H. Kuo, D. B. Holland, Y. Chen, M. Ouyang, G. A. Blake, R. Zadoyan, and C.-L. Guo, “Wide-field optical sectioning for live-tissue imaging by plane-projection multiphoton microscopy,” J. Biomed. Opt. **16**, 116009 (2011). [CrossRef]

**13. **C. Ventalon and J. Mertz, “Quasi-confocal fluorescence sectioning with dynamic speckle illumination,” Opt. Lett. **30**, 3350–3352 (2005). [CrossRef]

**14. **A. Egner and S. Hell, “Time multiplexing and parallelization in multifocal multiphoton microscopy,” J. Opt. Soc. Am. A **17**, 1192–1201 (2000). [CrossRef]

**15. **J. Jahns and K.-H. Brenner, *Microoptics: from Technology to Applications*, vol. v. 97 (Springer, 2004).

**16. **M. Born and E. Wolf, *Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light*, 6th ed. (Cambridge University Press, 1997).

**17. **D. Lim, K. K. Chu, and J. Mertz, “Wide-field fluorescence sectioning with hybrid speckle and uniform-illumination microscopy,” Opt. Lett. **33**, 1819–21 (2008). [CrossRef]

**18. **D. Lim, T. N. Ford, K. K. Chu, and J. Mertz, “Optically sectioned in vivo imaging with speckle illumination hilo microscopy,” J. Biomed. Opt. **16**, 016014 (2011). [CrossRef]

**19. **R. Heintzmann and P. A. Benedetti, “High-resolution image reconstruction in fluorescence microscopy with patterned excitation,” Appl. Opt. **45**, 5037–5045 (2006). [CrossRef]

**20. **E. Hecht, *Optics*, 4th ed. (Addison-Wesley, 2002).