Wigner analysis of three dimensional pupil with finite lateral aperture

Hsi-Hsun Chen; Se Baek Oh; Xiaomin Zhai; Jui-Chang Tsai; Liang-Cai Cao; George Barbastathis; Yuan Luo

doi:10.1364/OE.23.004046

1. Introduction

The use of volume holographic (VH) pupils in imaging systems includes applications such as hyper-spectral image acquisition [1], profilometric imagers [2, 3], unhindered imaging capability under broadband illumination, and multi-focal microscopic imaging [4, 5]. Compared to conventional 2D imaging pupils such as thin gratings and lenses, the VH pupil at a 4-f imagers’ Fourier plane acts as a 3D optical element, providing unique opportunities to process the optical field in three spatial dimensions as well as wavelength [6, 7].

Utilizing holographic recording techniques and advances in material [5], the 3D VH imaging pupil can be further functionally engineered to achieve wavelength-coded [8] and phase-coded [9] VH imaging gratings. Several mathematical models, including coupled wave theory [10], k-sphere formulation [11], and weak diffraction approximation [12], have been used for 3D pupil analysis. Here we are interested in phase space (i.e. space-spatial frequency) information transport [13] between the object and image planes through the 3D pupil. Recent, prior work [14] introduced the Wigner Distribution Function (WDF) analysis method for that purpose; however, infinite lateral aperture was assumed and, therefore, several important features of the resulting diffraction patterns and point spread function (PSF) were missed.

In this paper, we analyze a weakly diffracting VH pupil of arbitrary thickness and width. We derive basic volume holographic properties of angular selectivity in phase space as well as phase space information of diffracted beams at various locations along the VH pupil, we also analyze of the relationship between the response of 3D pupils in phase space and point spread function as it impacts imaging performance we compare; the simulation results with experiments under different aperture conditions validating the theory.

2. Theory

Fig. 1 illustrates a 4-f imager consisting of a 3D pupil, recorded as a volume hologram by two mutually coherent plane waves with an inter-beam angle $θ_{s}$ . Figure 1(a) shows the recording geometry, whereas Fig. 1(b) is the corresponding k-sphere diagram [11] illustrating the Bragg-matching relationship between the diffracted and probe wavevectors ${\vec{k}}_{d}$ and ${\vec{k}}_{p}$ , respectively, and the grating vector $\vec{K}$ . The Bragg-matching condition allows the 3D VH pupil to perform its unique spatial-spectral filtering and imaging functions [8, 11, 12, 15].

Fig. 1 . Illustration of (a) volume holographic recording geometry, (b) k-sphere diagram probed and recorded using wavevectors kp and kd, respectively, for the recorded grating vector $\vec{K}$ ; (c) 4-f imaging system geometry with 3D pupil, where f is the focal length, a and L are the lateral and axial apertures, respectively, and red dash-line denotes the Fourier plane. To obtain a clear pupil, the VH would be removed from the system and be replaced by a simple aperture of the same size, located at the Fourier plane.

Download Full Size | PDF

To simplify the analysis of the 3D VH pupil onto Winger space, we rather consider the impulse response of a 4-f VH system based on Fig. 1(c). Although we are neglecting the un-diffracted (i.e. 0-order) beam for convenience, the geometry is convenient since its 3D spatial transfer function of the 4-f VH system has been reported previously without missing the general Bragg properties of the 3D pupil [12]. In Fig. 1(c), the 3D pupil is located at the Fourier plane and has finite aperture width $a$ , and thickness L. In the vicinity of the Fourier plane, the input and output WDF $W_{3}$ , $W_{4}$ respectively to the VH pupil can be related to input and output WDF $W_{1}$ , $W_{2}$ , respectively as [14, 16-17]

W_{3} (x_{3}, u_{3}) = W_{1} (- λ f u_{3}, \frac{L}{2 f} u_{3} + \frac{x_{3}}{λ f}), \begin{matrix}  \end{matrix} W_{4} (x_{4}, u_{4}) = W_{2} (λ f u_{4}, \frac{L}{2 f} u_{4} - \frac{x_{4}}{λ f}),

where

λ

is the operation wavelength. As known in Wigner space through coordinate transforms, the WDF of

W_{3}

and

W_{4}

are related as [18]:

W_{4} (x_{4}, u_{4}) = \int K_{3 D} (x_{4}, u_{4}; x_{3}, u_{3}) W_{3} (x_{3}, u_{3}) d x_{3} d u_{3} .

where

K_{3 D}

is the kernel describing the action of the 3D pupil.

For a volume hologram,

K_{3 D} (x_{4}, u_{4}; x_{3}, u_{3}) = K_{V H I S} (λ f u_{4}, \frac{L}{2 f} u_{4} - \frac{x_{4}}{λ f}; - λ f u_{3}, \frac{L}{2 f} u_{3} + \frac{x_{3}}{λ f}),

\begin{array}{l} K_{V H I S} (x_{2}, u_{2}; x_{1}, u_{1}) = \iint d {x^{'}}_{2} d {x^{'}}_{1} \exp (- i 2 π (u_{2} {x^{'}}_{2} - u_{1} {x^{'}}_{1})) \times \\ \begin{matrix} \begin{matrix}  \end{matrix} \end{matrix} h_{V H I S} (x_{2} + \frac{{x^{'}}_{2}}{2}; x_{1} + \frac{{x^{'}}_{1}}{2}) h_{V H I S}^{*} (x_{2} - \frac{{x^{'}}_{2}}{2}; x_{1} - \frac{{x^{'}}_{1}}{2}), \end{array}

and

h_{V H I S}

is the coherent PSF of the volume holographic imaging system. Assuming the 1st-order Born approximation [12],

h_{V H I S}

is obtained as

h_{V H I S} (x_{2}; x_{1}) = \sin c (\frac{a}{λ f} (x_{1} + x_{2} - f θ_{s})) \times \sin c (\frac{L}{2 λ f^{2}} (x_{1}^{2} - x_{2}^{2} + f^{2} θ_{s}^{2})),

Note that for L = 0, Eq. (5) becomes identical to the diffraction-limited PSF of a conventional 4-f imaging system with clear pupil. Substituting Eq. (5) into Eq. (4), we obtain

\begin{array}{l} K_{3 D} (x_{4}, u_{4}; x_{3}, u_{3}) = \iint d {u^{'}}_{3} d {u^{'}}_{4} \exp {- i 2 π [{u^{'}}_{4} (\frac{L λ}{2} u_{4} - x_{4}) + {u^{'}}_{3} (\frac{L λ}{2} u_{3} + x_{3})]} \times \\ \sin c {a [- u_{3} + u_{4} + \frac{1}{2} (- {u^{'}}_{3} + {u^{'}}_{4}) - \frac{θ_{s}}{λ}]} \sin c {a [- u_{3} + u_{4} - \frac{1}{2} (- {u^{'}}_{3} + {u^{'}}_{4}) - \frac{θ_{s}}{λ}]} \times \\ \sin c {\frac{L λ}{2} [{(u_{3} + \frac{{u^{'}}_{3}}{2})}^{2} - {(u_{4} + \frac{{u^{'}}_{4}}{2})}^{2} + {(\frac{θ_{s}}{λ})}^{2}]} \times \sin c {\frac{L λ}{2} [{(- u_{3} + \frac{{u^{'}}_{3}}{2})}^{2} - {(u_{4} - \frac{{u^{'}}_{4}}{2})}^{2} + {(\frac{θ_{s}}{λ})}^{2}]} . \end{array}

In the case of 2D clear pupil i.e. L = 0, the Wigner function

K_{2 D}

can be simplified as

K_{2 D} (x_{4}, u_{4}; x_{3}, u_{3}) = \frac{a}{2} δ (x_{3} - x_{4}) Λ (\frac{x_{4}}{a / 2}) \sin c {(2 a - 4 | x_{4} |) (u_{3} - u_{4})},

where

Λ

denotes the triangle function.

3. Simulation and analysis

Fig. 2(i-iv) shows the $x_{3} - u_{3}$ planar cross sections of the 4-dimensional $K_{3 D}$ at various coordinates $x_{4}$ at the back plane of 3D pupil. The width of the visible slit in the $x_{3} - u_{3}$ space is proportional to the pupil thickness L indicating the Fourier-conjugate relationship with the width of the visible slit in the input $x_{1}$ plane due to Bragg selectivity. On the contrary, this thickness-induced visible slit does not appear in the case of a clear pupil, as shown in Fig. 3, owing to $δ (x_{3} - x_{4})$ in $K_{2 D}$ . Moreover, there exists an offset between the center of probing location $x_{3}$ and the maximum diffracted location $x_{4}$ in VHIS; this is because of the Bragg-matching condition. For example, in case (i) of Fig. 2 its diffracted beam is from $x_{4} = 0$ with its corresponding $x_{3} = - 0.25$ , while in case (iv) its diffracted beam originates from $x_{4} = 0.8$ (outside the recorded area) with its corresponding $x_{3} = 0.375$ . This is true independent of the probe beam position, i.e. near the center or edge of the VH pupil. Due to the lateral extent along the $x_{3}$ plane as well as finite lateral pupil width, while the probe beam moves towards the pupil boundary, peak values of the $K_{3 D}$ kernel gradually reduce (right-hand side color-bar in Fig. 2). This indicates that Bragg diffraction effects become weaker at the vicinity of the pupil edge, which may be thought of as the grating having reduced effective thickness near the edges. Furthermore, Fig. 4 shows the cross section between $x_{3}$ and $x_{4}$ under Bragg match condition with $u_{3}$ and $u_{4} = θ_{s} / λ$ at various coordinates. The projection along $x_{4}$ in Fig. 4(c) provides Gaussian-like shape, whose width is wider than that along $x_{3}$ , and there is an offset between the peak values. These show agreement with previous finding in Fig. 2.

Fig. 2 Wigner function of input corresponding to different output location of diffracted beam. The output location $(x_{4})$ is (i) 0, (ii) 0.25,(iii), 0.5,(iv) 0.8mm with same diffracted angle $u_{4} = θ_{s} / λ$ , where $a = 1 m m$ , $L = 1 m m$ , $λ = 500 n m$ and $θ_{s} = 30^{\circ}$ Red dash frame denotes the edge of the recorded region inside the holographic material.

Download Full Size | PDF

Fig. 3 Wigner function of clear pupil, where $a = 1 m m$ , $L = 0 m m$ , $λ = 500 n m$ Note that the spatial coordinates is defined as $x_{3} = x_{4}$ owing to $δ (x_{3} - x_{4})$ in Eq. (7).

Download Full Size | PDF

Fig. 4 The illustration of Bragg match condition (a), (b) is the mapping against $x_{3}$ and $x_{4}$ , where $u_{3} = 0$ and $u_{4} = θ / λ$ as well as Bragg match. (c) is the projection on $x_{3}$ (orange curve) and $x_{4}$ (blue curve). The grid white lines in (b) denote the locations of peak value mapping against $x_{3}$ and $x_{4}$ .

Download Full Size | PDF

Rather than analyzing the Bragg selectivity of the 3D pupil at $x_{1} - x_{2}$ spatial coordinates [12], $u_{3} - u_{4}$ phase coordinates reveal a more detailed transformation relationship at various locations on the 3D pupil. Figs. 5(a-c) show the comparison among $K_{2 D}$ at the center, $K_{3 D}$ kernel at the center, and $K_{3 D}$ kernel near the edge of the 3D pupil in $u_{3} - u_{4}$ phase coordinates, respectively. The values along the diagonal lines $(u_{3} = u_{4})$ in Fig. 5 represent angular response under geometrical conjugate condition of $x_{2} = - x_{1}$ . Figure 5(a) shows a shift-invariant system with a 2D clear pupil so that the values along diagonal line are constant. Owing to the shift-variance property induced by strong Bragg angular selectivity in the 3D pupil, the values along diagonal line in Fig. 5(b) degrade significantly away from its centered Bragg-matched region. Note that the peak values along the output coordinate $u_{4}$ in both Figs. 5(b, c) shift to $u_{4} = θ_{s} / λ$ , depending on the Bragg-matched angle. In addition, the bandwidth of $K_{3 D}$ near the pupil center (in Fig. 5(b)) is narrower than that near the pupil edge (in Fig. 5 (c)), while maximum value on Fig. 5(b) (i.e. Bragg diffraction efficiency) of the kernel in the vicinity of the pupil center is much higher than that in Fig. 5(c). In other words, Bragg diffraction is weaker and angular selectivity is less sensitive near the edge, which show agreement with the previous findings in Fig. 2.

Fig. 5 The phase to phase representation as the beam illuminates on (a) 2D clear pupil, (b) vicinity of center of 3D pupil, where $x_{3} = - 0.25 m m$ and $x_{4} = 0 m m$ , and (c) edge of 3D pupil, where $x_{3} = 0.375 m m$ and $x_{4} = 0.8 m m$ . Here, $a = 1 m m$ , $λ = 500 n m$ and $θ_{s} = 30^{\circ}$ , and the thickness is (a) $L = 0 m m$ for the clear pupil, (b, c) $L = 1 m m$ for the volume hologram.

Download Full Size | PDF

Next, we investigate how the 3D pupil with thickness and finite width affects VHIS performance. Since WDF is a bilinear transform [19], its intensity projected along spatial frequency domain and the point spread function of a 4-f imager can be expressed as

I_{4} (x_{4}) = \int W_{4} (x_{4}, u_{4}) d u_{4}, \begin{matrix} {| P S F (x_{2}) |}^{2} = \int W_{4} (\frac{L}{2 f} x_{2} - λ f u_{2}, \frac{x_{2}}{λ f}) \end{matrix} d u_{2},

where

W_{4}

is calculated at the condition of

W_{3} (x_{3}, u_{3}) = δ (u_{3})

with coherent on-axis plane wave.

W_{4}

under different recording angles

θ_{s}

is plotted in Fig. 6 and Fig. 7, respectively. When

θ_{s}

increases,

I_{4}

decays faster and its full width at half maximum (FWHM) becomes narrower [Figs. 6(b), 7(b)]. This can be thought of as fringes near the 3D pupil edge with a larger slant, which may result in larger portion of smaller effective thickness as well as weaker Bragg diffraction effects inside the recorded area.

Fig. 6 (a) Wigner function of $x_{4} - u_{4}$ . (b) $I_{4}$ for a 4-f imager with 3D pupil (marked in red line), and with 2D clear pupil (marked in blue line). (c) PSF of the 4-f imager with 3D (marked in red line), and 2D pupil (marked in blue line), with $a = 1 m m$ , $L = 1 m m$ , $λ = 500 n m$ and $θ_{s} = 30^{\circ}$ .

Download Full Size | PDF

Fig. 7 (a) Wigner function of $x_{4} - u_{4}$ . (b) $I_{4}$ for a 4-f imager with 3D pupil (marked in red line), and with 2D clear pupil (marked in blue line). (c) PSF of the 4-f imager with 3D (marked in red line), and 2D pupil (marked in blue line), with $a = 1 m m$ , $L = 1 m m$ , $λ = 500 n m$ and $θ_{s} = 68^{\circ}$ .

Download Full Size | PDF

Compared to a conventional 4-f imager with 2D pupil in Figs. 6(b) and 7(b), $I_{4}$ of the 3D pupil in VHIS has wider cut-off range and decays gradually towards the pupil boundary, while $I_{4}$ with 2D pupil, as expected, has rectangular profile with narrower cut-off range. Therefore, Figs. 6(c) and 7(c) show that PSF in VHIS exhibits much fewer side-lobes than that in a conventional 4-f imager with 2D pupil. These PSF simulation results also indicate that in VHIS there exist fewer side-lobes along the degeneracy axis (i.e. red line in Figs. 6(c) and 7(c)), and more side-lobes along the non-degeneracy axis (i.e. blue line in Figs. 6(c) and 7(c)), which acts as a conventional 4-f imager with 2D pupil.

4. PSF measurements and conclusion

In addition, an experimental setup was built to further verify our theoretical analysis with PSF performance measurements. The 3D pupil implemented by VH was located in the Fourier plane based on Fig. 1(c). Thickness of the 3D pupil in the PSF measurements was ~1.2mm with $θ_{s} = 68^{\circ}$ , and Bragg-matched operation wavelength was 532nm (Millennia, Spectra-Physics). The 3D pupil was illuminated by a collimated beam, and an iris was placed in front of the VH to adjust the pupil width at will. The measured PSF performance at two different values of 3D pupil width is shown in Figs. 8(a) and 8(c). In the measurement, we deliberately saturated the intensity of main lobe in order to identify the clearer side lobe, and the simulated PSF was also saturated for fair comparison. Simulation results are shown respectively in Figs. 8(b) and 8(d), and show good agreement of the lobes’ distribution with the measurements. Indeed, both cases result in fewer side-lobes along the degeneracy axis ( $x_{2}$ axis). The correlation coefficients between simulated and measured 1-D PSFs passing through peak of main lobe along x and y axis were also calculated according to

C_{M,S} = \frac{\sum_{i = 1}^{n} (M_{i} - 〈 M 〉) (S_{i} - 〈 S 〉)}{\sqrt{\sum_{i = 1}^{n} {(M_{i} - 〈 M 〉)}^{2} \sum_{i = 1}^{n} {(S_{i} - 〈 S 〉)}^{2}}},

where M and S denotes the measurements and simulations of 1-D PSFs, respectively. The i and n is the index and total pixel number of M and S, respectively, and

〈 〉

is the average operator. For 1-D PSF along

x_{2}

axis, the correlation coefficients are 0.95 and 0.92, for

a = 0.57 m m

and

a = 1 m m

respectively. In case of 1-D PSF along

y_{2}

axis, the correlation coefficients are 0.93 and 0.97, for

a = 0.57 m m

and

a = 1 m m

respectively. The corresponding intensity profiles passing through the peak of the main lobe along

x_{2}

and

y_{2}

axis in both simulations and measurements are shown in Figs. 8(i-iv). The possible error could be occurred from neglecting aberration and absorption in our model. This proposed phase space derivation and simulation can be also potentially applied to analyze 3D pupils recorded by spherical reference wave. Future work will include a study of aberration in the imager, and absorption caused by the VH in VHIS.

Fig. 8 The measurement (a, c) and simulation (b, d) of PSF of VHIS, (a) and (b) is PSF for $a = 0.57 m m$ ; (c) and (d) is PSF for $a = 1 m m$ . (i, ii) intensity profile passing through the peak of the main lobe along $x_{2}$ and $y_{2}$ for $a = 0.57 m m$ , and (iii, iv) intensity profile passing through the peak of the main lobe along $x_{2}$ and $y_{2}$ for $a = 1 m m$ .

Download Full Size | PDF

Acknowledgments

We acknowledge financial support from Ministry of Science and Technology (100-2218-E-002-026-MY3, 102-2218-E-002-013-MY3, 103-2221-E-002-156-MY3), funded by National Health Research Institutes (NHRI-EX103-10220EC), National Taiwan University (102R7832), and National Taiwan University Hospital (MP03, UN102-17). The authors are also grateful to the following sponsors: the U.S. National Institutes of Health (NIH RO1CA134424), National Research Foundation, Singapore through the Singapore MIT Alliance for Research and Technology's BioSystems and Micromechanics Inter-Disciplinary Research Group (015824-039).

References and links

1. W. Liu, D. Psaltis, and G. Barbastathis, “Real-time spectral imaging in three spatial dimensions,” Opt. Lett. 27(10), 854–856 (2002). [CrossRef] [PubMed]

2. A. Sinha and G. Barbastathis, “Volume holographic imaging for surface metrology at long working distances,” Opt. Express 11(24), 3202–3209 (2003). [CrossRef] [PubMed]

3. W. Sun, “Profilometry with volume holographic imaging,” Ph.D. dissertation (Massachusetts Institute of Technology, 2006).

4. Y. Luo, J. Castro, J. K. Barton, R. K. Kostuk, and G. Barbastathis, “Simulations and experiments of aperiodic and multiplexed gratings in volume holographic imaging systems,” Opt. Express 18(18), 19273–19285 (2010). [CrossRef] [PubMed]

5. Y. Luo, P. J. Gelsinger, J. K. Barton, G. Barbastathis, and R. K. Kostuk, “Optimization of multiplexed holographic gratings in PQ-PMMA for spectral-spatial imaging filters,” Opt. Lett. 33(6), 566–568 (2008). [CrossRef] [PubMed]

6. C.-C. Sun, “Simplified model for diffraction analysis of volume holograms,” Opt. Eng. 42(5), 1184–1185 (2003). [CrossRef]

7. H.-T. Hsieh, W. Liu, F. Havermeyer, C. Moser, and D. Psaltis, “Beam-width-dependent filtering properties of strong volume holographic gratings,” Appl. Opt. 45(16), 3774–3780 (2006). [CrossRef] [PubMed]

8. Y. Luo, S. B. Oh, and G. Barbastathis, “Wavelength-coded multifocal microscopy,” Opt. Lett. 35(5), 781–783 (2010). [PubMed]

9. S. B. Oh, Z.-Q. J. Lu, J.-C. Tsai, H.-H. Chen, G. Barbastathis, and Y. Luo, “Phase-coded volume holographic gratings for spatial-spectral imaging filters,” Opt. Lett. 38(4), 477–479 (2013). [CrossRef] [PubMed]

10. H. W. Kogelnik, “Coupled Wave Theory for Thick Hologram Gratings,” Bell Syst. Tech. J. 48(9), 2909–2947 (1969). [CrossRef]

11. P. Wissmann, S. B. Oh, and G. Barbastathis, “Simulation and optimization of volume holographic imaging systems in Zemax®,” Opt. Express 16(10), 7516–7524 (2008). [CrossRef] [PubMed]

12. G. Barbastathis, “The Transfer Function of Volume Holographic Optical Systems,” in Photorefractive Materials and Their Applications 3, P. Günter, and J.-P. Huignard, eds., pp. 51–76 (Springer, 2007).

13. G. Situ and J. T. Sheridan, “Holography: an interpretation from the phase-space point of view,” Opt. Lett. 32(24), 3492–3494 (2007). [CrossRef] [PubMed]

14. S. B. Oh and G. Barbastathis, “Wigner distribution function of volume holograms,” Opt. Lett. 34(17), 2584–2586 (2009). [CrossRef] [PubMed]

15. Y. Luo, V. R. Singh, D. Bhattacharya, E. Y. S. Yew, J.-C. Tsai, S.-L. Yu, H.-H. Chen, J.-M. Wong, P. Matsudaira, P. T. C. So, and G. Barbastathis, “Talbot holographic illumination nonscanning (THIN) fluorescence microscopy,” Laser and Photonics Reviews 8(5), L71–L75 (2014). [CrossRef]

16. S. B. Oh, “Volume holographic pupils in ray, wave, statistical optics, and Wigner space,” Ph.D. dissertation (Massachusetts Institute of Technology, 2009).

17. K. H. Brenner and J. Ojeda-Castañeda, “Ambiguity Function and Wigner Distribution Function Applied to Partially Coherent Imagery,” Optica Acta: International Journal of Optics 31(2), 213–223 (1984). [CrossRef]

18. W. Pan, “Double Wigner distribution function of a first-order optical system with a hard-edge aperture,” Appl. Opt. 47(1), 45–51 (2008). [CrossRef] [PubMed]

19. M. J. Bastiaans, “Wigner distribution function and its application to first-order optics,” J. Opt. Soc. Am. 69(12), 1710–1716 (1979). [CrossRef]

Wigner analysis of three dimensional pupil with finite lateral aperture

Abstract

1. Introduction

2. Theory

3. Simulation and analysis

4. PSF measurements and conclusion

Acknowledgments

References and links

Cited By

Figures (8)

Equations (9)

Optics Express