Simulation and optimization of volume holographic imaging systems in ...

Simulation and optimization of volume holographic imaging systems in Zemax  Patrick Wissmann1,2† , Se Baek Oh1† , and George Barbastathis 1 1 Department

of Mechanical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, , Cambridge, MA 02139 USA wissmann, sboh, [email protected]

2 Department

of Mechanical Engineering, RWTH Aachen University of Technology, Aachen, Germany [email protected] † Authors contributed equally to this work.

Abstract: We present a new methodology for ray-tracing analysis of volume holographic imaging (VHI) systems. Using the k-sphere formulation, we apply geometrical relationships to describe the volumetric diffraction effects imposed on rays passing through a volume hologram. We explain the k–sphere formulation in conjunction with ray tracing process and describe its implementation in a Zemax  UDS (User Defined Surface). We conclude with examples of simulation and optimization results and show proof of consistency and usefulness of the proposed model. © 2008 Optical Society of America OCIS codes: (090.7330) Volume gratings;(220.4830) Optical systems design.

References and links 1. A. Sinha, W. Sun, T. Shih, and G. Barbastathis, “Volume Holographic Imaging in Transmission Geometry,” Appl. Opt. 43, 1533–1551 (2004). 2. www.zemax.com. 3. J. W. Goodman, Introduction to Fourier Optics, 2nd ed. (McGraw Hill, 1996). 4. W. Sun, K. Tian and G. Barbastathis, “Hyper-spectral imaging with volume holographic lenses,” in Conference on Laser and Electro-Optics (CLEO), Paper CFP2 (2005). 5. W. Liu, D. Psaltis and G. Barbastathis, “Real-time spectral imaging in three spatial dimensions,” Opt. Lett. 27, 854–856 (2002). 6. G. Barbastathis, ”The transfer function of volume holographic optical systems,” Photorefractive Materials and Their Applications,(Springer-Verlag, 2007) Vol. 3. 7. A. Sinha and G. Barbastathis, “Volume holographic telescope,” Opt. Express 27, 1690–1692 (2002). 8. G. Barbastathis, M. Balberg and D. J. Brady, “Confocal microscopywith a volume holographic filter,” Opt. Lett. 24, 811–813 (1999). 9. G. Barbastathis and D. Psaltis, ”Volume Holographic Multiplexing Methods,” in Holographic Data Storage (Springer, 2000). 10. C. Kitel, Introduction to Solid-State Physics, 6th ed. (John Wiley, 1986). 11. F. A. Jenkins and H. E. White, Fundamentals of Optics, 4th ed. (McGraw-Hill, 1976). 12. J. R. Meyer-Arendt, Introduction to Classical and Modern Optics, 4th ed. (Prentice-Hall, 1989). 13. R. K. Kostuk, Multiple grating reflection volume holograms with application to optical interconnects, Ph. D. Thesis, (Stanford University, 1986) 14. R. R. A. Syms and L. Solymar, ”Analysis of Volume Holographic Cylindrical Lenses,” J. Opt. Soc. Am. 72, 179–186 (1982)

#92941 - $15.00 USD

(C) 2008 OSA

Received 22 Feb 2008; revised 25 Apr 2008; accepted 28 Apr 2008; published 9 May 2008

12 May 2008 / Vol. 16, No. 10 / OPTICS EXPRESS 7516

1. Introduction Volume holographic imaging (VHI) is a multi-dimensional imaging technique utilizing wavefront information stored in volume holograms [1]. Volume holograms are recorded by two or more mutually coherent wavefronts, whose interference pattern modulates the refractive index of a photo sensitive material. The modulation of the refractive index generates volumetric diffraction when the volume hologram is exposed to a probe wavefront, reaching maximum diffraction efficiency if the entire probe wavefront satisfies the Bragg condition [3]. VHI systems make use of the volume hologram as a 3D optical element, exploiting the strong shift variance due to the Bragg condition to achieve depth selective imaging and hyperspectral imaging [4, 5]. We have evaluated the system response of various VHI systems [6] and demonstrated multi–dimensional imaging systems such as the volume holographic telescope and microscope [1, 7, 8]. To design a VHI system, Bragg diffraction from volume holograms has to be analyzed. Analytical models of volume holograms have already been reported [1, 6, 9]. However, they are based on wave optics and not convenient to use in design process. In this paper, we propose a ray–tracing implementation of volume hologram modeling. In modern optical system design, ray tracing software packages are frequently used because they provide more functionality and flexibility with various analysis and optimization features, which makes the design process much more efficient and reliable. In section 2, we revisit the geometrical representation of volumetric diffraction, the k-sphere. In section 3, we describe its implementation in a Zemax  [2] UDS. In section 4, we investigate the consistency and accuracy of the approach and compare the simulation results with our conventional method. We present illustrative examples to demonstrate the design capabilities, focusing on the optimization of aberrated systems. Section 5 discusses advantages and limitation of our ray-tracing implementation and addresses future developments. 2. k–sphere formulation A volume hologram is considered as a phase grating, but being significantly thicker than conventional diffraction gratings, it is subject to volumetric diffraction and shows strong wavefront selectivity when exposed to a probe beam. We have previously derived a mathematical model of volumetric diffraction of scalar fields assuming Born’s first order approximation [9]. The k–sphere formulation is the geometrical representation of the above model; it originates from a method used in crystallography and solid-state physics [10] and has the potential to simplify the calculation of volumetric diffraction by geometric construction.

(a) Construction of the grating vector from the recording wave vectors

(b) Bragg matched reconstruction

(c) Bragg mismatch

(d) Wavelength degeneracy

Fig. 1. Geometric construction of volumetric diffraction using the k–sphere

#92941 - $15.00 USD

(C) 2008 OSA

Received 22 Feb 2008; revised 25 Apr 2008; accepted 28 Apr 2008; published 9 May 2008

12 May 2008 / Vol. 16, No. 10 / OPTICS EXPRESS 7517

In Fig. 1, we explain the concept of the k–sphere formulation. In this particular case, we assume that the volume hologram is recorded by plane waves for simplicity. Later we will describe how to use the k–sphere formulation for volume holograms recorded by non–plane waves. We denote by λ f and λp the wavelengths of the recording and probe beams, respectively. Let kf , ks denote the reference and signal   wavevectors, respectively (|k f | = |ks | = 2π /λf ); while kp denotes the probe wavevector ( kp  = 2π /λp). As shown in Fig. 1(a), the grating vector K g is obtained as Kg = ks − kf . (1) To determine the diffracted field in the Bragg regime, we need to examine the vector k p + Kg [9]. The direction of the wavevector k d of the diffracted field is specified by demanding that the lateral component of k d equals the lateral component of k p + Kg, and that the tip of k d lies on the k–sphere (i.e., the condition for propagating field); that is, kd · xˆ =

(kp + Kg) · xˆ ,

kd · yˆ = (kp + Kg) · yˆ , |kd | = 2π /λp.

(2) and

(3) (4)

These constraints are obtained assuming that the lateral size of the hologram is much larger than the hologram thickness. If k p + Kg satisfies the above conditions, then k d ≡ kp + Kg , and the probe is said to be Bragg–matched as shown in Fig. 1(b). If the Bragg matching condition is not satisfied, kd is still obtained from Eqs. (3)–(4) but the diffraction efficiency is decreased according to     L(kp + Kg − kd ) · zˆ Lδ kd · zˆ η ∝ sinc2 = sinc2 , (5) 2π 2π where L is the hologram thickness. The vector

δ kd ≡ kp + Kg − kd

(6)

is referred to as the Bragg mismatch vector; clearly, Bragg matching requires δ k d = 0, whereas if δ kd = 0 the diffraction efficiency is determined directly by the magnitude of δ k d (Fig. 1(c)). It is evident from the geometrical construction that the Bragg condition can be fulfilled at an infinite number of combinations of probe beam wavevectors k p and probe wavelengths λ p . The locus of the kp vector’s tip as function of λ p constitutes the wavelength degeneracy curve of the volume hologram (Fig. 1(d)). It should be noted that the k–sphere formulation can be applied either to slanted (θ s = θf ) or unslanted holograms (θ s = θf ). Moreover, reflection holograms (|θ s − θf | > 90◦) as well as transmission holograms (|θ s − θf | < 90◦ ) can be described by the k–sphere formulation. Typically, the k–sphere formulation is not frequently used to analyze Bragg diffraction of holograms recorded by non–plane waves, because the direction of the grating vector (K g in Fig. 1) changes locally over the hologram surfaces. However, ray tracing software can handle this situation. Once the local grating vector where the incident ray hits the volume hologram is known, the diffracted ray can be computed by the k–sphere formulation. Thus we implement two options to compute grating vectors; 1) holograms recorded by either planes waves or spherical waves and 2) holograms that have arbitrary grating vectors. In the first case, once the angle and/or the position of a point source are known, then the grating vector can be computed analytically. In the second case, we discretize the hologram and assign local grating vectors. A similar approach has been used to model non–plane wave reference holograms using coupled mode theory [13, 14]. Our implementation interpolates the assigned grating vectors to compute local grating vectors. #92941 - $15.00 USD

(C) 2008 OSA

Received 22 Feb 2008; revised 25 Apr 2008; accepted 28 Apr 2008; published 9 May 2008

12 May 2008 / Vol. 16, No. 10 / OPTICS EXPRESS 7518

(a) Non-slanted plane–wave erence/signal transmission (θf = θs = 20◦ )

refVH

(b) Slanted plane wave reference/signal transmission VH (θf = 0◦ , θs = 20◦ )

(c) Plane wave reference/signal reflection VH (θf = 0◦ , θs = 240◦ )

(d) Spherical wave reference–plane wave signal transmission VH (dr = −60 mm, θs = 20◦ )

(e) Spherical wave reference–plane wave signal reflection VH (dr = −60 mm, θs = 220◦ )

(f) Spherical wave reference–spherical wave signal reflection VH (dr = −60 mm, ds = −60 mm, θs = 40◦ )

Fig. 2. Layouts of several VHI systems with different volume hologram. Note that images are produced by Zemax ) with a realistic bi–convex lens (BK7, f = 30 mm). dr and ds are distances from a volume hologram to point sources. They are used to specify the origin of spherical wavefronts.

Figure 2 shows various geometries of VHI systems with different volume holograms. These figures are generated by Zemax  with our volume hologram model. Figures 2(a), (b), and (d) are transmission holograms and Figs. 2(c), (e) and (f) are reflection holograms. Figure 2(a) uses an unslanted hologram and the rest of them use slanted ones. Not only plane waves but also spherical waves can be used as shown in Figs. 2(d), (e), and (f). 3. Implementation of the k–sphere formulation in a Zemax  User Defined Surface We implemented the k–sphere formulation in a User Defined Surface (UDS) of Zemax  . The UDS feature allows making a user’s own special optical elements, where users can define special refraction, reflection, diffraction laws. A UDS is represented by a Microsoft Windows  Dynamic Link Library (DLL), the file which contains functions invoked by Zemax  at design and simulation time. Besides Zemax  , other commercial ray tracing software suites also provide customizable features. Hence, our approach can be easily integrated in other software packages with proper modifications. Figure 3 describe the overall process. First, we consider a volume hologram as a surface. The thickness of the volume hologram is taken into account when the diffraction efficiency is computed. Zemax  chooses a single ray arriving on the hologram UDS. Then Zemax  computes

#92941 - $15.00 USD

(C) 2008 OSA

Received 22 Feb 2008; revised 25 Apr 2008; accepted 28 Apr 2008; published 9 May 2008

12 May 2008 / Vol. 16, No. 10 / OPTICS EXPRESS 7519

Fig. 3. Schematic of volume hologram UDS implementation in Zemax , where Kg is the local grating vector and η0 is the local diffraction efficiency at (x, y, z), kp and kd are k– vectors of the incident and diffracted ray, respectively. kd is computed by the k–sphere formulation. The intensity I of the diffracted ray is η0 × η , where η is computed by Eq. (10).

the local grating vector K g and the k–vector of the incident ray by using direction cosines (l, m, n) and wavelength as ⎛ ⎞ l 2π (7) kp = ⎝ m ⎠ . λp n Using the geometrical relationships of the k–sphere, the diffracted ray vector k d and the Bragg mismatch vector δ k d are identified to be ⎞ ⎛ kp,x + Kg,x ⎟ ⎜ kp,y + Kg,y kd = ⎝

(8) ⎠ 2π 2 2 2 ( λp ) − kp,x + Kg,x − kp,y + Kg,y and

δ kd = kp + Kg − kd .

(9)

Continuing with the implementation of the UDS, we normalize the diffracted ray vector k d and directly write the result into the ray data structure of diffracted rays. The attenuation of diffraction efficiency as a result of Bragg mismatch is taken into account by δ kd in (5). Since the normalized diffraction efficiency is equivalent to the relative transmittance of the surface for the given ray, the corresponding parameter of the data structure is rel surf tran, defined as rel surf tran = η0 (x, y, z)sinc2

Lδ k · zˆ d , 2π

(10)

where η0 (x, y, z) is the local diffraction efficiency. This feature allows changing absolute efficiency locally. Typically manufacturers provide hologram efficiency maps, hence the data can be substituted into η 0 (x, y, z). Note that kp , Kg , kd and rel surf tran are repeatedly computed for all rays being traced.

#92941 - $15.00 USD

(C) 2008 OSA

Received 22 Feb 2008; revised 25 Apr 2008; accepted 28 Apr 2008; published 9 May 2008

12 May 2008 / Vol. 16, No. 10 / OPTICS EXPRESS 7520

4. Examples In this section, we present simulated results generated by Zemax  and our volume hologram model. First we reproduce the experimental data [1] by simulation, and then we address new directions of using Zemax  in VHI system design. 4.1. Depth selectivity In the case of the VHI system with a plane wave reference hologram as shown in Figs. 2(a) and (b), maximum Bragg diffraction is achieved if the light from a point source is collimated and the resulting plane wave enters the volume hologram at the Bragg matched angle. If the point source leaves the front focal plane and moves along the optical axis, due to the defocus, then the wavefront entering the volume hologram is no longer a plane wave, causing Bragg mismatch except for a Bragg matched location on the optical axis. This example demonstrates how the Bragg condition can be exploited to achieve depth selectivity with volume holograms. We used a 4–f VHI system as shown in Fig. 2(b) (focal lengths f 1 = f2 = 50.2 mm) with a slanted volume hologram (radius a = 3.5 mm, θ f = 0◦ , θs = 12◦ , L=2 mm) at the Fourier plane, and the wavelength of recording and probe beam is 532 nm. Figure 4(a) shows the ray-tracing simulated plot of the total diffracted intensity at the detector plane versus the defocus of the point source. The longitudinal PSF is defined[1] as the normalized diffraction efficiency with respect to defocus. Figure 4(b) shows an analytical simulation, overlayed with an experimental result [1] for the comparison.

(a) Ray-tracing in Zemax

(b) Analytical model and experiment (from Ref. [1])

Fig. 4. Longitudinal PSF of a 4–f VHI system (see text for parameters)

4.2. Effect of geometrical aberrations in volumetric diffraction An important advantage of the ray-tracing approach over analytical methods is the straightforward fashion in which aberrations can be taken into account. In the analytical model [1], the paraxial approximation was used and rotational symmetry was assumed for the wavefront inside the volume hologram. This prohibited aberrations other than spherical to be taken into account. Generally, we can assume that aberration reduces the Bragg diffraction efficiency, and consequently the depth and wavelength selectivity of VHI systems may be degraded. Therefore, the ability to model aberration effects is essential for making full use of analysis and optimization tools in the design process. Figure 5 presents a simple example showing the effect of aberration to the longitudinal PSF. In this case, we used the ideal volume hologram represented by one K g vector and two different objective lenses: first, an ideal thin lens (which, observing the paraxial approximation does not introduce any aberration in Zemax  ); and second, a realistic bi-convex lens made of BK7. The effective focal lengths of the two lenses are identical at the wavelength of 532 nm. #92941 - $15.00 USD

(C) 2008 OSA

Received 22 Feb 2008; revised 25 Apr 2008; accepted 28 Apr 2008; published 9 May 2008

12 May 2008 / Vol. 16, No. 10 / OPTICS EXPRESS 7521

Fig. 5. Longitudinal PSFs with paraxial and aberrated lens. ( f1 =50.2 mm, L=2 mm, θs = 23.1◦ , λf = λp = 532 nm )

The aberrated lens yields slightly lower Bragg diffraction than the paraxial lens; however, unexpectedly the FWHM(Full Width Half Maximum) of the longitudinal PSF is improved by 10% after the normalization. This implies that the depth resolution (defined as FWHM of the longitudinal PSF) can be engineered by managing aberrations. We now present another simple example to demonstrate the capabilities of our ray tracing model in conjunction with the built-in analysis features of Zemax  . In the paraxial 4–f system, we replace the first lens with an aberrated lens of constant focal length and variable q factor defined by R1 + R2 , (11) q= R1 − R2 where R1 , R2 are the radii of curvature of the front and back surfaces of a lens, respectively [11]. We conducted a simulation where we changed the both radii of curvature of the surfaces while maintaining its focal length. As is well known, the severity of the spherical aberration is a function of q in this case. The following movie shows the transformation of the first lens through different stages of shape factor, the resulting diffraction pattern and a ray-fan diagram to quantify the severity of aberration. 4.3. Optimization example As in any optical system, VHI requires careful optimization of the overall optical elements, including the objective lens used to transform the object wavefront before it is diffracted by the hologram. In this section, we demonstrate a simple optimization of a VHI system with respect to aberration. We first run the optimization for both a regular 4–f system and VHI system with the same merit function. Intuitively, the optimization would yield different result because the characteristics of both systems are somewhat different. Then we try to obtain similar result with different merit function, which additional constraints are required in VHI optimization. We first used the default merit function, which minimizes the overall optical path length for the evaluated rays from 18 predefined points. As the parameters to be optimized, we choose both radii of curvature of the aberrated lens. Column “DFM” (Default Merit Function) of Table 1 shows the result after 1, 000 iterations. In this particular example, the point source is on-axis so that spherical aberration is dominant. With a singlet, the minimum spherical aberration is achieved when the shape factor is qmin.SA = #92941 - $15.00 USD

(C) 2008 OSA

−2(n2 − 1) , n+2

(12)

Received 22 Feb 2008; revised 25 Apr 2008; accepted 28 Apr 2008; published 9 May 2008

12 May 2008 / Vol. 16, No. 10 / OPTICS EXPRESS 7522

Diffraction Image

Lens Shape

Ray Fan

Lens Shape q = 1.597400

Fig. 6. Diffraction image, lens shape and ray fan. (1.4MB)

Parameters R1 [mm] R2 [mm] Shape factor q [N.A.] Merit function value [N.A.] Diffraction efficiency [%]

4–f system Initial DFM 30.5 123.892 -30.5 -19.774 0 0.724 20.132 2.231 – –

Initial 30.5 -30.5 0 21.167 3.566

VHI system DFM DFM+IMAE 34.416 68.815 -31.362 -21.910 0.0464 0.5170 10.869 11.891 23.746 33.264

Table 1. Optimization of 4–f and VHI system using the default merit function. The last column is computed by the IMAE operand in the merit function.

where n is the refractive index [12]. In the example, we use BK7 glass (n=1.5195 @ 532.8 nm), qmin.SA is 0.7437. Note that in the regular 4–f system, the optimized shape factor (0.7240) is almost same as qmin.SA , while the shape factor of VHI is only 0.0464. This is because of the merit function we chose, all rays diffracted by the hologram have the same weighting values with respect to the optimization goal. It is obvious that the optimization results in different shape factor for the VHI system. Instead of using the default merit function, a more practical choice is to increase the intensity of the diffracted beams. Hence, we modified the merit function to give less weight to the rays attenuated due to Bragg mismatch. In Zemax  , this can be done by defining the merit function to maximize the IMAE operand, which measures the integrated intensity at the detector plane. The rightmost column labeled “DFM+IMAE” of Table 1 shows how the modified merit function leads to an improved overall diffraction efficiency of the VHI system. Note the optimized shape factor (0.5170) with IMAE operand is much closer to q min.SA than the value (0.0464) computed with the default merit function only. 5. Discussion and conclusions We described a ray-tracing approach for the simulation and optimization of VHI systems. Our model for volumetric diffraction is based on the k–sphere formulation, whose applicability is verified by comparison with results of our previous model [9]. We implemented a Ze-

#92941 - $15.00 USD

(C) 2008 OSA

Received 22 Feb 2008; revised 25 Apr 2008; accepted 28 Apr 2008; published 9 May 2008

12 May 2008 / Vol. 16, No. 10 / OPTICS EXPRESS 7523

max User Defined Surface to demonstrate how the model can be easily integrated into current commercial ray-tracing software. Our method can simulate various volume holograms including reflection/transmission and slanted/unslated ones. Moreover, our model computes the local grating vector analytically or numerically from given data, which allows simulating not only plane/spherical wave reference holograms but also holograms with arbitrary grating vectors. Using the local efficiency parameter η 0 , we also account for absolute diffraction efficiency changes resulting from material properties or index variation,etc. Through several examples, we demonstrated how ray-tracing based analysis and optimization tools are used in conjunction with our model to assist in the design process of VHI systems. We showed that our model is ready to simulate aberrations in VHI systems and provided examples of how optimization tools can be applied to mitigate the effects of aberration or even benefit from them. We also want to mention that our approach can be integrated within other ray tracing software packages with proper modification. However, there is one limitation: the k–sphere formulation assumes that the lateral extent of the volume hologram is infinite. Currently, we have verified the model for spatially uniform volumetric gratings, e.g. an ideal plane wave reference volume hologram. It is left for future research to extend the model for non-uniform gratings, such as spherical reference wave holograms, holograms recorded with aberrated beams, etc. Acknowledgment The authors would like to thank Wojciech Matusik and Wenyang Sun for useful discussions. This research was partially funded by Mitsubishi Electric Research Laboratories (Cambridge, MA) and the DARPA Montage research program.

#92941 - $15.00 USD

(C) 2008 OSA

Received 22 Feb 2008; revised 25 Apr 2008; accepted 28 Apr 2008; published 9 May 2008

12 May 2008 / Vol. 16, No. 10 / OPTICS EXPRESS 7524