## Abstract

We propose a three-dimensional phase contrast digital holographic microscopy. The object to be observed is a low-contrast transparent refractive index distribution sample, such as biological tissue. Low contrast phase objects are converted to high contrast images through the microscopy we propose. In order to gain high three-dimensional resolution, the direction of pump plane wave is scanned, and separate holographic images produced at each angle are acquired and decoded into complex amplitude in Fourier space. The three-dimensional image is reconstructed in a computer from all information acquired through the system. The resolution in the direction of the optical axis is increased by utilizing a 4π configuration of objective lenzes.

© 2007 Optical Society of America

## 1. Introduction

Recently, studies on biomedical science have evolved remarkably, and some novel microscopies for biological material have been developed energetically. For example, the confocal microscopy is widely used in biology, medical science, and pharmaceutical science. The development of confocal microscopy has resulted in the ability to observe living, unstained biological tissue [1–3]. While the image can be formed by autofluorescence from the specimen in confocal microscopy, the radiance of autofluorescence typically suffers from low signal values. The diffracted wave reflected or transmitted through the sample also enables one to obtain the image. With confocal reflection microscopy, however, the signal from refractive index boundaries is emphasized and the structure containing several phase layers cannot be visualized precisely because the 0-order light does not interfere with the diffracted wave at the pinhole right before the detector [1]. While confocal transmission microscopy can resolve an absorptive object that has an absorption coefficient distribution, it cannot form a high contrast image of a phase object, so that a detector with an excessively wide dynamic range is required [1]. Therefore, in most cases, the biological tissue is stained for fluorescence detection, which is, however, a major disadvantage in observing living tissue. Furthermore, it is difficult to separate the essential signal generated in the fluorescent dyes from background autofluorescence.

Various tools have been developed that do not require staining the specimen before observation. Some multi-photon microscopies detect the second harmonic generation (SHG) signal [4–6], the third harmonic generation (THG) signal [7,8], or the coherent anti-Stokes Raman scattering (CARS) signal [9,10], which are generated through nonlinear optical interaction between the molecules in the specimen and the light. Multi-photon microscopy, however, requires a relatively expensive high power pulse laser as a light source. While phase contrast microscopy, which can convert a phase object into a high contrast irradiance image, is one of the most powerful tools to observe biological tissue, the three-dimensional (3-D) character at a phase object cannot be resolved [11,12]. Digital holographic microscopy, which can provide the 3-D images by numerically calculating the complex amplitude of wavefront diffracted by the object using a hologram recorded through a digital camera, does not have a sufficient resolution [13–16]. Finally, although optical coherence tomography is an exquisite device for biological tissue, the development of a broadband light source is indispensable in order to gain high resolution in the depth direction [17,18].

In this paper, a powerful tool is proposed that can resolve an unstained specimen three-dimensionally and is operated easily. This technique permits 3-D visualization of objects at any viewing angle. 3-D image reconstruction is accomplished by using digital holographic technique combined with scanning the direction of the wave vector of the pump plane wave in order to increase the resolution. The digital holographic technique is used to acquire the complex amplitude of the object waves produced by the pump beams at every incident angle. The resolution of the microscopy system using the 3-D optical transfer function (*OTF*) and an estimate of the error in object reconstruction are also discussed.

## 2. Setup of 3-D phase contrast holographic microscopy

The object is a low contrast refractive index distribution. Most biological tissues satisfy this condition. Unlike ordinary microscopy, which comprises an image-formation optical system, the interference pattern between an object wave and a reference wave is recorded with a two-dimensional digital camera placed in the Fourier plane of the object. The complex amplitude of the object wave is calculated in a computer using a digital holographic technique [13–16]. In order to gain high 3-D resolution, the direction of the pump beam is scanned two-dimensionally by a scanning mirror system, and the 3-D image of the local scattering cross section in the object is reconstructed from the complex amplitude information of the object waves produced over the angle range of the pump beam.

Figure 1 shows a schematic of the 3-D holographic microscopy instrument. Any two successive lenses in Fig. 1 compose infinite image-formation systems. A stable near-infrared laser with a sufficiently long coherence length is used as the light source. Since the direction of the polarization of the laser is perpendicular to the plane of the drawing (S-polarization), the beam is reflected by polarized beam splitter (PBS1). The beam is collimated through a tube lens and a primary objective, and the object is illuminated by a plane wave whose angle is scanned two-dimensionally by the scanning mirror. While the maximum angle of the pump beam after propagating through the primary objective is approximately 64 degrees, which corresponds to a numerical aperture (NA) of 1.2 in water, the maximum beam angle after reflection by the scanning mirror is about 1.15 degree, since the focal length of the tube lens is 60 times as long as that of the primary objective. Two primary objectives, which have the identical specifications, face each other, and a sample holder is set between them. After the 0-order light is transmitted through the sample, it propagates through another primary objective and is transmitted through the polarized beam splitter (PBS1) after being changed into P-polarization by a half wave plate (*λ*/2 -A). The 0-order light is collimated by the tube lens and reflected by the scanning mirror again, which changes the angle of the 0-order light back into the original angle of the pump beam emerging from the light source. After reflection by the scanning mirror, the 0-order light converges onto a micromirror of about 200*μ*m in diameter. While a 50/50 beam splitter inevitably leads to light loss in the detection path, it can be reduced by using an 80/20 beam splitter. The majority of the 0-order light is reflected by the micromirror and is used as a reference wave for the hologram. The remainder of the 0-order light is transmitted through the micromirror and advances toward the digital camera. A certain quantity of 0-order light is required to form the image. The transmitted 0-order light is used for compensation of phase error that will be discussed later.

The transmitted component of the scattered wave emerging from the sample travels along the same optical path as the 0-order light and passes through the outside of the micromirror, whereas the low frequency component of the scattered wave is intercepted by the micromirror. The micromirror is set on the optical axis at the Fourier plane of the two primary objectives (the middle plane in the sample). The scattered wave passing through the outside of the microrirror is transmitted through PBS2 and propagated toward the digital camera set at the Fourier plane of the two primary objectives. The location of the micromirror is also conjugate to the digital camera except for the reference beam path.

The reflected component of the scattered wave emerging from the sample propagates back along the optical path and is reflected by PBS1. It then follows the identical optical path with the scattered transmitted wave. After reflected by the scanning mirror, the reflected scattered wave passes around the outside of the microrirror. The S-polarization is reflected by PBS2, and the polarization is converted the P-polarization by a half wave plate (*λ*/ 2 -B), and it propagates toward the digital camera. In the same way as the transmitted component, a few frequency components of the scattered reflected wave in the vicinity of the 0-order light are intercepted by the micromirror. While the well-known high-NA depolarization effect occurs in the primary objectives, two polarizers placed on the transmission (P-polarization) and reflection (S-polarization) sides ensure the linear polarization.

The 0-order light transmitted in the micromirror separates into two paths due to PBS2 after its polarization is converted to circular polarization by a micro quarter wave plate (*λ*/4 ) of the same size as the micromirror. A part of the 0-order light transmitted through the PBS2 propagates to the digital camera along with the transmitted component of the scattered wave and the other part of the 0-order light reflected by the PBS2 propagates with the reflected component of the scattered wave. The two divided 0-order light parts is used to compensate the phase error for the transmitted and reflected components, which will be described in detail in Section 4. Although the micro quarter wave plate can be removed for a simple setup, in this case the phase error of the reflected component cannot be minimized. The digital camera can be divided into two parts, which are used for the transmission and reflection waves. Since the transmitted and reflected scattered waves and the reference wave are P-polarized on the observation plane of the digital camera, the transmission and reflection waves interfere with the reference wave, and holographic images are recorded through the digital camera. Since the central point of the scanning mirror is placed at the conjugate position of the middle point of the sample on the optical axis, the incident angle of the transmission and reflection waves is perpendicular to the digital camera while the pump beam scanned. The circular area of the scattered wave at the digital camera is shifted during scanning. The 0-order light is fixed at the central positions of the observation planes in each of part, as shown in Fig. 2. The incident angle of the reference wave onto the digital camera never varies during the scanning, since the 0-order light is always reflected at the center of the micromirror after being reflected twice by the scanning mirror. Holographic images of the transmission and reflection waves are recorded for each incident direction of the pump wave.

The sample is set in a sample holder with two circular apertures fixed between the two primary objectives, as shown in Fig. 3. The primary objective in this system is a Fourier lens [19], in which the aberration of pupil is corrected properly and the condition *h* = *F* sin *θ* is satisfied, where *h* represents the height of the principal ray on the pupil plane, *F* denotes the focal length, and *θ* stands for the angle of the principal ray emitted from the central point on the object plane to the optical axis (see Fig. 3). In particular, a water immersion objective lens with numerical aperture (*NA*) 1.2 can be designed and produced easily for the Fourier lens, whose specification is that the angular field of view of the incident plane wave on the object side is 1.2 in *NA* in water and the *NA* on the detector side is 0.0075 in air. In this case, if the focal length of the tube lens is 200mm and the magnification of the primary objective is 60X, the focal length of the primary objective is about 4.4mm in water, and the field of view in the sample is approximately 50 *μ*m. The digital camera is set to be conjugate to the pupil plane of the primary objective. If we assume that the optical system between the pupil plane of the primary objective and the detector plane comprises the image-formation system of magnification *β*, we obtain the relation

where *μ* represents the size of a pixel of the digital camera and *N* denotes the number of the pixels along a side of the digital camera. The length of the side of detector plane is assumed to be equivalent to *μN*, which is twice as long as the diameter of a circular area of the scattered wave on the digital camera.

## 3. Diffraction due to low contrast phase object

Before the explanation of the algorithm for the object reconstruction, it is significant to consider diffraction due to a phase object that has a low contrast refractive index distribution. The first Born approximation is applied, in which only the 0-order light and the first-order scattered wave are taken into account, as the 0-order light is overwhelmingly intense and the higher-order scattered waves are negligible. The refractive index distribution in the specimen is given by *n*(*x*) = *n*
_{0} + *εg*(*x*), where *n*
_{0} is a scalar constant representing the average refractive index, *x* stands for the three-dimensional position in the specimen, ε is a constant that is sufficiently small *ε* ≫ 1, and *g*(*x*) is a scalar function related to the net refractive index distribution and its range is -1 < *g*(*x*) < 1. If the object is illuminated by a plane wave, the scattered wave appears. Since the digital camera is set conjugate to the pupil plane of the primary objective, one can presume that the wave vector of the scattered light is recorded.

#### 3.1 Electric field of the scattered wave

After the incident plane wave *E*
^{(i)}(*x*′,*t*) with wave vector *n*
_{0}
*k*
_{0} is scattered by the object, the total electric field including the 0-order light and the scattered wave on an infinite-radius reference sphere satisfies the following equation [20]

where *x*′ is the three-dimensional position on the reference sphere, *c* is the speed of the light in vacuum, *α*(*x*) is the polarizability relative to average refractive index *n*
_{0}, and *N*(*x*) denotes the number of molecules in a unit volume. The product *α*(*x*)*N*(*x*) is given by [20]

The effective electric field *E*′ (*x*,*t* - ∣*x*′-*x*∣*n*
_{0} / *c*) in Eq. (2) is

where *ω* is the angular frequency of the light used and *k*
_{0} = *ω*/*c* is the magnitude of the wave vector *k*
_{0}. Since the reference sphere is located far from the object, that is ∣*x*′-*x*∣ ≫ *λ*/2*π*, the integrand of the second term in Eq. (2) approximates [20]

where *λ* is the wavelength of light in vacuum, and *θ* is the angle between the direction of the effective electric field *E*
_{0}′(*x*) and the propagation direction of the scattered light exp[*i n*
_{0}
*k*
_{0} ∣ *x*′-*x*∣]/∣*x*′-*x*∣, and *φ* is azimuthal angle and *r* is radius in polar coordinates. While Eq. (5) implies the well-known high-NA depolarization effect, two polarizers placed on the transmission and reflection sides ensure the linear polarization of the scattered waves, as mentioned above. Hereafter, the polarization effects of scattering are not taken into account, namely, the vector electric field shown by Eq. (5) is considered to be the scalar field whose amplitude is equivalent to *E _{θ}* with an approximation of sin

*θ*= 1. Note that this is not the paraxial approximation. If only the scalar electric field is considered, the integrand of the second term in Eq. (2) is

Assuming the harmonic dependence of the scalar electric field results in *E*′(*x*′,*t*)→*E*
_{0}′(*x*′)*e*
^{-iωt} and *E*
^{(i)} (*x*′,*t*) →*E*
_{0}
^{(i)}(*x*′)*e*
^{-iωt}, and substitution of Eqs. (3) and (6) into Eq. (2) yields

Ignoring the second and higher order of *ε* (the first order Born approximation) and substitution of *E*
_{0}
^{(i)} (*x*′) = exp[*i n*
_{0}
*k*
_{0} · *x*] yields

$$\phantom{\rule{.5em}{0ex}}=\mathrm{exp}\left[i{n}_{0}{\mathit{k}}_{0}\xb7{\mathit{x}}^{\prime}\right]+\epsilon {k}_{0}^{2}\frac{{n}_{0}}{2\pi}\int g\left(\mathit{x}\right)\mathrm{exp}\left[i{n}_{0}{\mathit{k}}_{0}\xb7\mathit{x}\right]\frac{\mathrm{exp}\left[i{n}_{0}{k}_{0}\mid {\mathit{x}}^{\prime}-\mathit{x}\mid \right]}{\mid {\mathit{x}}^{\prime}-\mathit{x}\mid}{d}^{3}\mathit{x}.$$

Since the distance from any position ** x** in the sample to arbitrary position

**on the reference sphere is infinite, the second term of the right side in Eq. (8) is proportional to the Fourier transform of**

*x*′*g*(

*x*), which is measured at the digital camera.

#### 3.2 Diffraction efficiency

We will consider the diffraction efficiency in the case of the simple object that contains only a single spatial frequency, namely *g*(*x*) = cos(*Κ*· *x*) = exp[*i κ*·*x*]/2 + exp[-*i κ*·*x*]/2, where ** κ** is the three-dimensional grating vector. For simplicity, we will consider the diffraction due to only the positive frequency component exp[i

*κ*·

*x*]/2. Substitution of

*g*(

*x*) = exp[i

*κ*·

*x*]/2 into Eq. (8) yields

with

$$\phantom{\rule{.5em}{0ex}}=\frac{{n}_{0}}{4\pi}\epsilon {k}_{0}^{2}\int \mathrm{exp}\left[i\left\{\left({n}_{0}{k}_{x}+{\kappa}_{x}\right)\mathit{x}+\left({n}_{0}{k}_{y}+{\kappa}_{y}\right)y+\left({n}_{0}{k}_{z}+{\kappa}_{z}\right)z\right\}\right]\frac{\mathrm{exp}\left[i{n}_{0}{k}_{0}R\right]}{R}{d}^{3}\mathit{x},$$

where $R=\sqrt{{\left({x}^{\prime}-x\right)}^{2}+{\left({y}^{\prime}-y\right)}^{2}+{\left({z}^{\prime}-z\right)}^{2}}$ and *k _{x}*,

*k*and

_{y}*k*are the

_{z}*x*,

*y*, and

*z*components of

*k*

_{0}, respectively. A central position of the sample is defined as origin of the coordinate and

*z*direction is optical axis of the system. While the integral over

*x*and

*y*is conducted from negative infinity to positive infinity, the integral domain

**, which corresponds to the size of the sample, is assumed to be sufficiently small compared with the distance between the origin and the position**

*x***on the reference sphere,**

*x*′$$x\prime x+y\prime y+z\prime z\ll {x\prime}^{2}+{y\prime}^{2}+{z\prime}^{2},$$

which results in

Inserting Eq. (12) into the exponential function and $R\approx \sqrt{{x\prime}^{2}+{y\prime}^{2}+{z\prime}^{2}}$ into the dominator in Eq. (10) yields

$$\phantom{\rule{.5em}{0ex}}\times {\int}_{-\infty}^{\infty}\mathrm{exp}\left[i\left({n}_{0}{k}_{\mathrm{y}}+{\kappa}_{y}-{n}_{0}{k}_{y}\prime \right)y\right]\mathrm{exp}\left[i\frac{{n}_{0}{k}_{0}}{2r\prime}{y}^{2}\right]dy$$

$$\phantom{\rule{.5em}{0ex}}\times {\int}_{-\frac{L}{2}}^{\frac{L}{2}}\mathrm{exp}\left[i\left({n}_{0}{k}_{z}+{\kappa}_{z}-{n}_{0}{k}_{z}\prime \right)z\right]\mathrm{exp}\left[i\frac{{n}_{0}{k}_{0}}{2r\prime}{z}^{2}\right]\text{d}z,$$

where *L* represents a thickness of the sample in the *z* direction, and *k _{x}*′,

*k*′, and

_{y}*k*′ are

_{z}*x*,

*y*, and

*z*components of wave vector of the scattered wave

*k*

_{0}′, and

*k*

_{0}

*x*′/

*r*′=

*kx*′

*etc*

*. are used, and $r\prime =\sqrt{{x\prime}^{2}+{y\prime}^{2}+{z\prime}^{2}}\left(\gg \sqrt{{x}^{2}+{y}^{2}+{z}^{2}}\right).$*

*Now we will consider the integral over x in Eq. (13),*

*$${J}_{x}=\underset{r\prime \to \infty}{lim}{\int}_{-r\prime}^{r\prime}\mathrm{exp}\left[i\phantom{\rule{.2em}{0ex}}\mathrm{ax}\right]\mathrm{exp}\left[i\phantom{\rule{.2em}{0ex}}b{x}^{2}\right]\text{d}x,$$*

*with a = n
_{0}
k_{x} + κ_{x} - n
_{0}
k_{x}′ and b = n
_{0}
k
_{0} /(2r′). As r′ approaches infinity, b approaches zero. In order to execute the integral J_{x}, the integrand is multiplied by a factor exp[-βx
^{2}], where β is real and positive and approaches zero β → +0 (β ≪ b). Then Eq. (14) becomes*

*$${J}_{x}=\underset{\underset{\beta \to +0}{r\prime \to \infty}}{lim}{\int}_{-\mathrm{r\prime}}^{\mathrm{r\prime}}\mathrm{exp}\left[i\phantom{\rule{.2em}{0ex}}\mathrm{ax}\right]\mathrm{exp}\left[\mathrm{ib}{x}^{2}\right]\mathrm{exp}\left[-\beta {x}^{2}\right]dx$$*

$$\phantom{\rule{.5em}{0ex}}=\underset{\underset{\beta \to +0}{r\prime \to \infty}}{lim}\mathrm{exp}\left[-\frac{{a}^{2}}{4\left(\beta -\mathrm{ib}\right)}\right]{\int}_{-\mathrm{r\prime}}^{\mathrm{r\prime}}\mathrm{exp}\left[-{\left(\sqrt{\beta -\mathrm{ib}}x-\frac{\mathrm{ia}}{2\sqrt{\beta -\mathrm{ib}}}\right)}^{2}\right]dx.$$

$$\phantom{\rule{.5em}{0ex}}=\underset{\underset{\beta \to +0}{r\prime \to \infty}}{lim}\mathrm{exp}\left[-\frac{{a}^{2}}{4\left(\beta -\mathrm{ib}\right)}\right]{\int}_{-\mathrm{r\prime}}^{\mathrm{r\prime}}\mathrm{exp}\left[-{\left(\sqrt{\beta -\mathrm{ib}}x-\frac{\mathrm{ia}}{2\sqrt{\beta -\mathrm{ib}}}\right)}^{2}\right]dx.$$

*Note that $\sqrt{\beta -\mathrm{ib}}$ is a complex number with phase ${\phi}_{\beta ,b}=\frac{1}{2}{\mathrm{tan}}^{-1}\left(\frac{-b}{\beta}\right),$ which approaches -π/4 as β → +0 (β ≪ b). The change of variables $z=\sqrt{\beta -\mathrm{ib}}x-\frac{\mathrm{ia}}{\left(2\sqrt{\beta -\mathrm{ib}}\right)}$ allows one to write*

*$${J}_{x}=\underset{\underset{\beta \to +0}{r\prime \to \infty}}{lim}\mathrm{exp}\left[-\frac{{a}^{2}}{4\left(\beta -\mathrm{ib}\right)}\right]\frac{1}{\sqrt{\beta -\mathrm{ib}}}{\int}_{c}\mathrm{exp}\left[-{z}^{2}\right]dz,$$*

*where the contour C makes an angle φ_{βb} to the real axis. If β = (1/r′)^{α} (1 < α < 2), the
C and the real axis. Since there are no singularities of the integrand anywhere between C and the real axis, one can again deform the contour to the real axis,*

*$${J}_{x}=\underset{r\prime \to \infty}{lim}\mathrm{exp}\left[-\frac{{a}^{2}}{4\left(\beta -\mathrm{ib}\right)}\right]\frac{1}{\sqrt{\beta -\mathrm{ib}}}{\int}_{-\mathrm{r\prime}}^{\mathrm{r\prime}}\mathrm{exp}\left[-{x}^{2}\right]dx\text{}$$*

$$\phantom{\rule{.5em}{0ex}}=\underset{r\prime \to \infty}{lim}-\frac{\sqrt{\pi}}{\sqrt{\beta -\mathrm{ib}}}\mathrm{exp}\left[-\frac{i}{4b}{a}^{2}\right]\mathrm{exp}\left[-\frac{\beta}{4{b}^{2}}{a}^{2}\right]$$

$$\phantom{\rule{.5em}{0ex}}=\underset{\gamma \to \infty}{lim}\sqrt{i}\frac{\sqrt{\pi}}{\sqrt{b}}\mathrm{exp}\left[-\frac{i}{4b}{a}^{2}\right]\mathrm{exp}\left[-\gamma {a}^{2}\right].\text{}$$

$$\phantom{\rule{.5em}{0ex}}=\underset{r\prime \to \infty}{lim}-\frac{\sqrt{\pi}}{\sqrt{\beta -\mathrm{ib}}}\mathrm{exp}\left[-\frac{i}{4b}{a}^{2}\right]\mathrm{exp}\left[-\frac{\beta}{4{b}^{2}}{a}^{2}\right]$$

$$\phantom{\rule{.5em}{0ex}}=\underset{\gamma \to \infty}{lim}\sqrt{i}\frac{\sqrt{\pi}}{\sqrt{b}}\mathrm{exp}\left[-\frac{i}{4b}{a}^{2}\right]\mathrm{exp}\left[-\gamma {a}^{2}\right].\text{}$$

*Replacing a = n
_{0}
k_{x} + κ_{x} - n
_{0}
k_{x} and b = n
_{0}
k
_{0}/(2r′) yields*

*$${J}_{x}=\underset{\gamma \to \infty}{lim}\frac{\sqrt{i}\sqrt{\pi}\sqrt{2r\prime}}{\sqrt{{n}_{0}{k}_{0}}}\mathrm{exp}\left[-i\frac{r\prime {\left({n}_{0}{k}_{x}+{\kappa}_{x}-{n}_{0}{k}_{x}\prime \right)}^{2}}{2{n}_{0}{k}_{0}}\right]\mathrm{exp}\left[-\gamma {\left({n}_{0}{k}_{x}+{\kappa}_{x}-{n}_{0}{k}_{x}\prime \right)}^{2}\right],$$*

*and likewise the integral over y is*

*$${J}_{y}=\underset{\gamma \to \infty}{lim}\frac{\sqrt{i}\sqrt{\pi}\sqrt{2r\prime}}{\sqrt{{n}_{0}{k}_{0}}}\mathrm{exp}\left[-i\frac{\mathrm{r\prime}{\left({n}_{0}{k}_{y}+{\kappa}_{y}-{n}_{0}{k}_{y}\prime \right)}^{2}}{2{n}_{0}{k}_{0}}\right]\mathrm{exp}\left[-\gamma {\left({n}_{0}{k}_{y}+{\kappa}_{y}-{n}_{0}{k}_{y}\prime \right)}^{2}\right].$$*

*When r′ approaches infinity r′ → ∞, the integral over z in Eq. (13) is*

*$${J}_{z}=L\mathrm{sin}c\left[\frac{L}{2}({n}_{0}{k}_{z}+{\kappa}_{z}-{n}_{0}{k}_{z}\prime \right],$$*

*where sinc[ x] =sin x/x. Substitution of these integrals j_{x}, j_{y}, J_{z} into Eq. (13) yields*

*$${F}_{0}\prime (r\prime )=\underset{\gamma \to \infty}{lim}\frac{\mathrm{i\epsilon}{\phantom{\rule{.2em}{0ex}}k}_{0}}{2}\mathrm{exp}\left[i{n}_{0}{k}_{0}\mathrm{r\prime}\right]\mathrm{exp}\left[-i\frac{\mathrm{r\prime}}{2{n}_{0}{k}_{0}}\left\{{\left({n}_{0}{k}_{x}+{\kappa}_{x}-{n}_{0}{{k}_{x}}^{\prime}\right)}^{2}+{({n}_{0}{k}_{y}+{\kappa}_{y}-{n}_{0}{k}_{y}\prime )}^{2}\right\}\right]$$*

$$\times \mathrm{exp}\left[-\gamma \left\{{\left({n}_{0}{k}_{x}+{\kappa}_{x}-{n}_{0}{k}_{x}\prime \right)}^{2}+{({n}_{0}{k}_{y}+{\kappa}_{y}-{n}_{0}{k}_{y}\prime )}^{2}\right\}\right]\text{}$$

$$\times L\mathrm{sin}c\left[\frac{L}{2}({n}_{0}{k}_{\mathrm{z}}+{\kappa}_{z}-{n}_{0}{k}_{z}\prime )\right]$$

$$=\left\{\begin{array}{c}\frac{\mathrm{i\pi \epsilon L}}{\lambda}\mathrm{exp}\left[i{n}_{0}{k}_{0}\mathrm{r\prime}\right]\mathrm{sin}c\left[\frac{L}{2}\left({n}_{0}{k}_{z}+{\kappa}_{z}-{n}_{0}{k}_{z}\prime \right)\right],\\ \left(\mathrm{as}\phantom{\rule{.2em}{0ex}}{n}_{0}{k}_{x}+{\kappa}_{x}-{n}_{0}{k}_{x}\prime =0,{n}_{0}{k}_{y}+{\kappa}_{y}-{n}_{0}{k}_{y}\prime =0\right)\\ 0,\phantom{\rule{.3em}{0ex}}\left(\mathrm{others}\right).\end{array}\right.$$

$$\times \mathrm{exp}\left[-\gamma \left\{{\left({n}_{0}{k}_{x}+{\kappa}_{x}-{n}_{0}{k}_{x}\prime \right)}^{2}+{({n}_{0}{k}_{y}+{\kappa}_{y}-{n}_{0}{k}_{y}\prime )}^{2}\right\}\right]\text{}$$

$$\times L\mathrm{sin}c\left[\frac{L}{2}({n}_{0}{k}_{\mathrm{z}}+{\kappa}_{z}-{n}_{0}{k}_{z}\prime )\right]$$

$$=\left\{\begin{array}{c}\frac{\mathrm{i\pi \epsilon L}}{\lambda}\mathrm{exp}\left[i{n}_{0}{k}_{0}\mathrm{r\prime}\right]\mathrm{sin}c\left[\frac{L}{2}\left({n}_{0}{k}_{z}+{\kappa}_{z}-{n}_{0}{k}_{z}\prime \right)\right],\\ \left(\mathrm{as}\phantom{\rule{.2em}{0ex}}{n}_{0}{k}_{x}+{\kappa}_{x}-{n}_{0}{k}_{x}\prime =0,{n}_{0}{k}_{y}+{\kappa}_{y}-{n}_{0}{k}_{y}\prime =0\right)\\ 0,\phantom{\rule{.3em}{0ex}}\left(\mathrm{others}\right).\end{array}\right.$$

*The diffracted wave has a phase shift of π/2 relative to 0-order light. When the Bragg’s condition n
_{0}
k
_{0} + κ - n
_{0}
k
_{0}′=0 is satisfied, the diffraction efficiency η is*

*which corresponds to the well known formula [21]. The diffraction efficiency of the phase object composed of a single spatial frequency is deduced assuming that only single scattering is taken into account and multi-scattering can be ignored if ε is sufficiently small. This result is approximately true as long as πεL/λ < 1. For example, ε < 0.006 when L = 45μm and λ = 850 nm. If πεL/λ ≫1, η becomes the square of sine function and oscillates, since there is a tight coupling between the 0-order light and the diffraction wave and the energy transfer between the two modes occurs periodically due to multi-scattering. Even if ε is somewhat large ε > 0.006, an image can be reconstructed through this measurement system. However, the image might be deformed compared with the original object.*

*4. Algorithm of object reconstruction*

*As mentioned above, the digital camera is conjugate to the pupil plane of the primary objective. The circular area of the scattered wave shifts two-dimensionally as the pump beam is scanned, with the position of the 0-order light fixed at the center of the camera area (see Fig. 2). The number of pixels of the digital camera for each transmitted and reflected components is N×N pixels. It is useful to consider the detection position of the scattered wave on the camera to correspond to wave number space, since the entrance pupil of the primary objective is located at infinity and the radius of the reference sphere is also infinity. The complex amplitude of the scattered wave on the detection plane is calculated numerically by utilizing a digital off axis holography technique [14–16].*

*The observable volume of the specimen is limited by two circular apertures on both sides of the sample. The diameters of the apertures are three times as large as that of the cross section of the pump beam, which is three times as large as the thickness of the sample. The pump beam is converted into a plane wave and is incident on the sample after converging onto a point in the pupil plane of the primary objective (PO1). The cross section of the pump beam on the aperture plane is a circle that is independent of the incident angle, which is one of the features of the Fourier lens (see Fig. 3). The NA of the scattered wave in the pupil plane of the primary objective (PO2) is confined by the apertures to be three times as large as that of the 0-order light.*

*In order to acquire the digital holographic image, the incident angle ϕ of the reference wave to the normal of the detection plane is adjusted as λ/sin ϕ = 4μ, which means that a unit cycle of the phase of the reference wave corresponds to four pixels of the digital camera. The difference in spatial frequency between a top and a bottom or a right side and a left side of the digital camera corresponds to 4NA/λ. In this case, the size of the object D in real space, which is calculated by Fourier transforming the wave number space, is*

*from Eq. (1). The diameter of the aperture R and the size of the object D have the relation R = 3D/4.*

*Now the algorithm of the object reconstruction is considered. After the holographic images of the transmitted and reflected scattered waves produced by the pump beam at every incident angle are acquired through the digital camera, the holographic images of N×N pixels are Fourier transformed in the computer. As a result of that computation, the complex amplitude of the object wave is obtained on one side at the matrix and a conjugated wave appears on another side, as shown in Fig. 4. While the object wave and the conjugated wave overlap, due to the relation of R =3D/4, the central areas of the object and conjugated waves of the size of N/4 do not overlap. Because the irradiance of the reference wave is adjusted by an attenuator to be about a hundred times (ten times in amplitude) as intense as that of the diffraction wave at the most intense position, only a central bright spot produced by the dc term appears. The autocorrelation of the object wave, which is supposed to appear around the central bright spot, is negligible. A section of the object wave of N/4 × N/4 elements is cut out and the inverse Fourier transform is applied to this 2-D matrix, which implies that the size of the object to be reconstructed is restricted to D/4 . As a result, the complex amplitude of the scattered wave on the pupil plane is obtained.*

*Since information about the amplitude of the scattered wave is contained in the circle of diameter N/8 in Fig. 4, the circle is cut out and projected onto the “partial sphere” that corresponds to the reference sphere and is mapped into wave number space. The partial sphere lies in a three-dimensional matrix of N/4 × N/4 × N/4 elements. The complex amplitudes in all elements other than those on the partial sphere are set to zero (see Fig. 4). In this step, the complex amplitude is projected in the z direction onto the digital partial sphere composed of a set of voxels across which the analog partial sphere passes. If a (x, y) element of the digital partial sphere consists of two voxels in the z direction, the amplitude is distributed equally to the two voxels. While this step can lead to some artifacts in real space, the error is reduced by using a partial sphere convoluted with a Gaussian function containing a few elements in the z direction. After the resultant partial sphere is digitalized, the amplitude is projected onto the digital partial sphere with a weight of the Gaussian distribution. In this case, the peripheral intensity of the image in real space in the z direction becomes weaker.*

*The orientation of each partial sphere in frequency space is such that the 0-order light is placed at the frequency origin, as shown in Fig. 5. This construction allows one to figure the physical properties of the system optical transfer function ( OTF). If the system has no error in phase for the measurement of the hologram corresponding to each partial sphere, the argument of the origin in the frequency space is zero. While the amplitude corresponding to the 0-order light lies on the partial sphere in the transmitted component, the amplitude of the 0-order light for the reflected component is located at the origin apart from the partial sphere, which has a hole at the corresponding position where the 0-order light arrives on the digital camera. The 0-order light arriving on the digital camera for the reflected component is subtracted from the partial sphere in the reflected component and is used for the phase error correction of the reflected component. The partial sphere for the reflected component and the 0-order light lie on the same sphere. The two 3-D matrices for the transmitted and reflected components are added, which results in a matrix composed of two partial spheres facing each other. This matrix of two partial spheres, which share the same radius and center of curvature, is hereafter referred to as a twin partial sphere.*

*We will consider two calculation methods to reconstruct the object. These two methods are equivalent for object reconstruction, as long as the sample is low contrast, which will be discussed in detail in the next section. In the first method, the twin partial sphere is Fourier transformed after the phase shift of π/2 is added only to the 0-order light, and the square of the modulus is calculated. The same calculation is performed for every direction of the pump plane wave. The final reconstructed object with the size of D/4 × D/4 × D/4 in real space is obtained by adding all of these calculated matrices. Since the square of the modulus is independent of shifts in the Fourier plane, a positioning of the twin partial sphere is not required. This method of the 3-D object reconstruction corresponds to image formation of conventional microscopy with an incoherent Kohler illumination system.*

*The relative optical path length between the transmitted scattered wave and the reference wave may change due to the temperature variation or vibration during the measurement. A difference in phase among the measurements for each hologram could be generated, but the error of the reconstructed object vanishes if only the transmitted component is taken into account, since we calculate the square of the modulus of Fourier transform of each complex amplitude in the frequency domain. However, the phase error remains due to the reflected component. The relative phase difference between the transmitted and reflected components of the twin partial sphere exists for the changes in relative optical length among the three paths of the transmitted and reflected diffraction waves and the reference wave. In order to minimize the phase error, the partial sphere for the transmitted (reflected) component is divided by the phase term of the 0-order light arriving on the digital camera for the transmitted (reflected) component before the two partial spheres are added to form the twin partial sphere. The optical path of the scattered wave is different from that of the 0-order light in the looped part, as outlined in Fig. 1, which includes the sample holder and PBS1. This looped part requires an extremely steady and stable mechanical design. Relative phase error caused by the looped part is discussed in Section 6.*

*In the second method, the object is reconstructed in the following way. After the complex amplitudes of the twin partial spheres for every pump direction are calculated, each twin partial sphere is translated so that the position corresponding to the 0-order light is placed at the origin of the 3-D matrix. The partial sphere of the transmitted (reflected) component in the twin partial sphere is divided by the phase term of the 0-order light arriving on the transmission (reflection) side of the digital camera in the same way as the first method. Thus, also in the second method, phase correction for each hologram by means of the 0-order light is achieved, and it ensures that the phase at the origin of the 3-D matrices is zero. Even with this phase correction, some relative phase error can remain in the reflected component, as described for the first calculation method. The total amplitude in the frequency domain is obtained by coherently adding the 3-D matrices of the twin partial spheres for every pump direction. Finally, the object intensity is reconstructed through a Fourier transform of the total amplitude after the phase shift of π/2 is added to the origin of the 3-D matrix, and taking the square of the modulus. These two techniques are closely related to filtered backprojection. The filtering is later done by the application of the inverse total OTF.*

*5. Three-dimensional image formation feature*

*Consider the 3-D pupil functions P
_{T}(f) for the transmission side and P
_{R}(f) for the reflection side which have a relation P
_{T}(f) = P*_{R}(-f), as shown in Fig. 6. For simplicity, it is assumed that P
_{T} (f) (P
_{R} (f)) is unity on the partial sphere for the transmission (reflection) side and zero outside the partial sphere. That is,*

*$${P}_{T}\left(\mathit{f}\right)=\{\begin{array}{c}1,\left(\mathrm{on}\phantom{\rule{.2em}{0ex}}\mathrm{shell\; of\; transmission\; side}\right)\\ 0,\left(\mathrm{others}\right).\end{array}\text{}$$*

$${P}_{R}\left(\mathit{f}\right)=\{\begin{array}{c}1,\left(\mathrm{on}\mathrm{shell}\mathrm{of}\mathrm{reflection}\mathrm{side}\right)\phantom{\rule{.2em}{0ex}}\\ 0,\left(\mathrm{others}\right).\end{array}$$

$${P}_{R}\left(\mathit{f}\right)=\{\begin{array}{c}1,\left(\mathrm{on}\mathrm{shell}\mathrm{of}\mathrm{reflection}\mathrm{side}\right)\phantom{\rule{.2em}{0ex}}\\ 0,\left(\mathrm{others}\right).\end{array}$$

*In practical systems, the pupil function is not uniform and is a function of NA. As NA approaches higher edges, the value of the pupil function shrinks. While we assume a uniform apodization hereafter, one can replace it with an arbitrary apodization, and also in this case the equations can be developed similarly to the derivation in this section.*

*The wave number of the pump wave 2 πf
_{0}(=n
_{0}
k
_{0}) is scanned over the pupil function P
_{T} (f) and the scattered wave emerging from the object is transmitted though the both sides of the pupil P(f) = P
_{T}(f) + P
_{R}(f). The amplitude of the 0-order light, which can be assumed to be a real number after the phase is shifted by π/2, is attenuated by the micromirror before arriving at the digital camera by the factor of 0 < a < 1. If the amplitude of the object is given by O(x), the amplitude on the twin partial sphere for a certain wave number of the pump wave 2πf
_{0} is*

*$$F\left(\mathit{f},{\mathit{f}}_{0}\right)=\left\{\int O\left(\mathit{x}\right){e}^{-i2\pi \mathit{f}\xb7\mathit{x}}d\mathit{x}-a\delta \left(\mathit{f}\right)\right\}\left\{{P}_{T}\left(\mathit{f}+{\mathit{f}}_{0}\right)+{P}_{R}\left(\mathit{f}+{\mathit{f}}_{0}\right)\right\}{P}_{T}^{*}\left({\mathit{f}}_{0}\right),$$*

*where the asterisk stands for the complex conjugate. Based on the first method, the irradiance of the image of the reconstructed object I
_{I}(x′) is given by*

*$${I}_{I}\left(\mathit{x\prime}\right)=\int \mid \int F\left(\mathit{f},{\mathit{f}}_{0}\right){e}^{i2\pi \phantom{\rule{.2em}{0ex}}\mathit{f}\xb7\mathit{x}\prime}{d\mathit{f}\mid}^{2}d{\mathit{f}}_{0}$$*

$$=\int [\int \{\int O\left({\mathit{x}}_{1}\right){e}^{-i2\pi {\mathit{f}}_{1}\xb7{\mathit{x}}_{1}}d{\mathit{x}}_{1}-\mathrm{a\delta}\left({\mathit{f}}_{1}\right)\}\left\{{P}_{T}\left({\mathit{f}}_{1}+{\mathit{f}}_{0}\right)+{P}_{\text{R}}\left({\mathit{f}}_{1}+{\mathit{f}}_{0}\right)\right\}{P}_{T}^{*}\left({\mathit{f}}_{0}\right){e}^{i2\pi {\mathit{f}}_{1}\xb7\mathit{x\prime}}d{\mathit{f}}_{1}]$$

$$\times [\int \left\{\int {O}^{*}\left({\mathit{x}}_{2}\right){e}^{i2\pi {\mathit{f}}_{2}\xb7{\mathit{x}}_{2}}d{\mathit{x}}_{2}-\mathrm{a\delta}\left({\mathit{f}}_{2}\right)\right\}\left\{{P}_{T}^{*}\left({\mathit{f}}_{2}+{\mathit{f}}_{0}\right)+{P}_{R}^{*}\left({\mathit{f}}_{2}+{\mathit{f}}_{0}\right)\right\}{P}_{T}\left({\mathit{f}}_{0}\right){e}^{-i2\pi {\mathit{f}}_{2}\xb7\mathit{x\prime}}d{\mathit{f}}_{2}]d{\mathit{f}}_{0}$$

$$=\int [\int O\left({\mathit{x}}_{1}\right){P}_{T}^{*}\left({\mathit{f}}_{0}\right){e}^{-i2\pi {\mathit{f}}_{0}\xb7\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)}\left\{{U}_{T}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)+{U}_{R}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)\right\}d{\mathit{x}}_{1}-a\left\{{P}_{T}\left({\mathit{f}}_{0}\right)+{P}_{R}\left({\mathit{f}}_{0}\right)\right\}{P}_{T}^{*}\left({\mathit{f}}_{0}\right)]$$

$$\times [\int {O}^{*}\left({\mathit{x}}_{2}\right){P}_{T}\left({\mathit{f}}_{0}\right){e}^{-i2\pi {\mathit{f}}_{0}\xb7\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)}\left\{{U}_{T}^{*}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)+{U}_{\text{R}}^{*}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)\right\}d{\mathit{x}}_{2}-a\left\{{P}_{T}^{*}\left({\mathit{f}}_{0}\right)+{P}_{\text{R}}^{*}\left({\mathit{f}}_{0}\right)\right\}{P}_{\text{T}}\left({\mathit{f}}_{0}\right)]d{\mathit{f}}_{0,}$$

$$=\int [\int \{\int O\left({\mathit{x}}_{1}\right){e}^{-i2\pi {\mathit{f}}_{1}\xb7{\mathit{x}}_{1}}d{\mathit{x}}_{1}-\mathrm{a\delta}\left({\mathit{f}}_{1}\right)\}\left\{{P}_{T}\left({\mathit{f}}_{1}+{\mathit{f}}_{0}\right)+{P}_{\text{R}}\left({\mathit{f}}_{1}+{\mathit{f}}_{0}\right)\right\}{P}_{T}^{*}\left({\mathit{f}}_{0}\right){e}^{i2\pi {\mathit{f}}_{1}\xb7\mathit{x\prime}}d{\mathit{f}}_{1}]$$

$$\times [\int \left\{\int {O}^{*}\left({\mathit{x}}_{2}\right){e}^{i2\pi {\mathit{f}}_{2}\xb7{\mathit{x}}_{2}}d{\mathit{x}}_{2}-\mathrm{a\delta}\left({\mathit{f}}_{2}\right)\right\}\left\{{P}_{T}^{*}\left({\mathit{f}}_{2}+{\mathit{f}}_{0}\right)+{P}_{R}^{*}\left({\mathit{f}}_{2}+{\mathit{f}}_{0}\right)\right\}{P}_{T}\left({\mathit{f}}_{0}\right){e}^{-i2\pi {\mathit{f}}_{2}\xb7\mathit{x\prime}}d{\mathit{f}}_{2}]d{\mathit{f}}_{0}$$

$$=\int [\int O\left({\mathit{x}}_{1}\right){P}_{T}^{*}\left({\mathit{f}}_{0}\right){e}^{-i2\pi {\mathit{f}}_{0}\xb7\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)}\left\{{U}_{T}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)+{U}_{R}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)\right\}d{\mathit{x}}_{1}-a\left\{{P}_{T}\left({\mathit{f}}_{0}\right)+{P}_{R}\left({\mathit{f}}_{0}\right)\right\}{P}_{T}^{*}\left({\mathit{f}}_{0}\right)]$$

$$\times [\int {O}^{*}\left({\mathit{x}}_{2}\right){P}_{T}\left({\mathit{f}}_{0}\right){e}^{-i2\pi {\mathit{f}}_{0}\xb7\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)}\left\{{U}_{T}^{*}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)+{U}_{\text{R}}^{*}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)\right\}d{\mathit{x}}_{2}-a\left\{{P}_{T}^{*}\left({\mathit{f}}_{0}\right)+{P}_{\text{R}}^{*}\left({\mathit{f}}_{0}\right)\right\}{P}_{\text{T}}\left({\mathit{f}}_{0}\right)]d{\mathit{f}}_{0,}$$

*where*

*$${U}_{T}\left(\mathit{x}\right)=\int {P}_{T}\left(\mathit{f}\right){e}^{i2\pi \mathit{f}\xb7\mathit{x}}d\mathit{f}$$*

$${U}_{R}\left(\mathit{x}\right)=\int {P}_{R}\left(\mathit{f}\right){e}^{i2\pi \mathit{f}\xb7\mathit{x}}d\mathit{f}.$$

$${U}_{R}\left(\mathit{x}\right)=\int {P}_{R}\left(\mathit{f}\right){e}^{i2\pi \mathit{f}\xb7\mathit{x}}d\mathit{f}.$$

*The Fourier transform of the 3-D entrance pupil function P(f) = P
_{T}(f) + P
_{R}(f) is equivalent to the 3-D coherent point spread function U(x) = U
_{T}(x) + U
_{R}(x), which has the relations*

*$${U}_{T}\left(\mathit{x}\right)={U}_{R}^{*}\left(\mathit{x}\right)$$*

$${U}_{T}\left(\mathit{x}\right)={U}_{T}^{*}\left(-\mathit{x}\right)\text{}$$

$${U}_{R}\left(\mathit{x}\right)={U}_{R}^{*}\left(-\mathit{x}\right).$$

$${U}_{T}\left(\mathit{x}\right)={U}_{T}^{*}\left(-\mathit{x}\right)\text{}$$

$${U}_{R}\left(\mathit{x}\right)={U}_{R}^{*}\left(-\mathit{x}\right).$$

*Since the pupil function has rotational symmetry to the optical axis,*

*$${I}_{I}\left(\mathit{x}\prime \right)=\int \int \gamma \left({\mathit{x}}_{1}-{\mathit{x}}_{2}\right)O\left({\mathit{x}}_{1}\right){O}^{*}\left({\mathit{x}}_{2}\right)\left\{{U}_{T}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)+{U}_{R}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)\right\}\left\{{U}_{T}^{*}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)+{U}_{R}^{*}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)\right\}d{\mathit{x}}_{1}d{\mathit{x}}_{2}$$*

$$-a\int O\left({\mathit{x}}_{1}\right){U}_{R}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)\left\{{U}_{T}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)+{U}_{R}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)\right\}d{\mathit{x}}_{1}$$

$$-a\int {O}^{*}\left({\mathit{x}}_{2}\right){U}_{T}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)\left\{{U}_{T}^{*}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)+{U}_{R}^{*}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)\right\}d{\mathit{x}}_{2}\text{}$$

$$+{a}^{2}\int {\mid {P}_{T}\left({f}_{0}\right)\mid}^{4}d{f}_{0},$$

$$-a\int O\left({\mathit{x}}_{1}\right){U}_{R}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)\left\{{U}_{T}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)+{U}_{R}\left(\mathit{x}\prime -{\mathit{x}}_{1}\right)\right\}d{\mathit{x}}_{1}$$

$$-a\int {O}^{*}\left({\mathit{x}}_{2}\right){U}_{T}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)\left\{{U}_{T}^{*}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)+{U}_{R}^{*}\left(\mathit{x}\prime -{\mathit{x}}_{2}\right)\right\}d{\mathit{x}}_{2}\text{}$$

$$+{a}^{2}\int {\mid {P}_{T}\left({f}_{0}\right)\mid}^{4}d{f}_{0},$$

*with*

*$$\gamma \left({\mathit{x}}_{1}-{\mathit{x}}_{2}\right)=\int {\mid {P}_{T}\left({\mathit{f}}_{0}\right)\mid}^{2}{e}^{i2\pi {\mathit{f}}_{0}\xb7\left({\mathit{x}}_{1}-{\mathit{x}}_{2}\right)}d{\mathit{f}}_{0,}$$*

*which is referred to as the 3-D mutual intensity.*

*Consider a low contrast object O(x) = 1 + ε
_{0}
o(x) where o(x) is proportional to g(x) and ε
_{0} ≪ 1. Inserting O(x) = 1 + ε
_{0}
o(x) into Eq. (29) yields*

*$${I}_{I}(\mathit{x}\prime )=\int \left\{1+{\epsilon}_{0}o\left(\mathit{x}\right)+{\epsilon}_{0}^{*}{o}^{*}\left(\mathit{x}\right)\right\}{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}d\mathit{x}$$*

$$+\int {\epsilon}_{0}o\left(\mathit{x}\right){U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)d\mathit{x}+\int {\epsilon}_{0}^{*}{o}^{*}\left(\mathit{x}\right){U}_{T}^{2}\left(\mathit{x}\prime -\mathit{x}\right)d\mathit{x}$$

$$-a\int \left\{1+{\epsilon}_{0}o\left(\mathit{x}\right)\right\}\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}$$

$$-a\int \left\{1+{\epsilon}_{0}^{*}{o}^{*}\left(\mathit{x}\right)\right\}\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{T}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}$$

$$+{a}^{2}\int {\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}d\mathit{x},$$

$$+\int {\epsilon}_{0}o\left(\mathit{x}\right){U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)d\mathit{x}+\int {\epsilon}_{0}^{*}{o}^{*}\left(\mathit{x}\right){U}_{T}^{2}\left(\mathit{x}\prime -\mathit{x}\right)d\mathit{x}$$

$$-a\int \left\{1+{\epsilon}_{0}o\left(\mathit{x}\right)\right\}\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}$$

$$-a\int \left\{1+{\epsilon}_{0}^{*}{o}^{*}\left(\mathit{x}\right)\right\}\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{T}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}$$

$$+{a}^{2}\int {\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}d\mathit{x},$$

*where Parseval’s theorem ∫∣ P
_{T}(f
_{0})∣^{2}df
_{0} = ∫∣U
_{T}(x′-x)∣^{2}dx and ∣P
_{T}(f
_{0})∣^{4} = ∣P
_{T}(f
_{0})∣^{2} are used in the last term and the second order terms of ε
_{0} are ignored. A further simple calculation leads to*

*$${I}_{I}\left(\mathit{x}\prime \right)=\left(1-a\right)\lfloor \left(1-a\right)\int {\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}d\mathit{x}$$*

$$+\int {\epsilon}_{0}o\left(\mathit{x}\right)\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}\text{}$$

$$+\int {\epsilon}_{0}^{*}{o}^{*}\left(\mathit{x}\right)\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{T}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}]$$

$$=\left(1-a\right)\int {\mid {U}_{T}\left(\mathit{x}\right)\mid}^{2}d\mathit{x}[\left(1-a\right)$$

$$+{\epsilon}_{0}\int \tilde{o}\left(\mathit{f}\right)\mathrm{OTF}\left(\mathit{f}\right){e}^{i2\pi \mathit{f}\xb7\mathit{x}}d\mathit{f}$$

$$+{\epsilon}_{0}^{*}\int {\tilde{o}}^{*}\left(\mathit{f}\right){\mathrm{OTF}}^{*}\left(\mathit{f}\right){e}^{-i2\pi \mathit{f}\prime \mathit{x}}d\mathit{f}],$$

$$+\int {\epsilon}_{0}o\left(\mathit{x}\right)\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}\text{}$$

$$+\int {\epsilon}_{0}^{*}{o}^{*}\left(\mathit{x}\right)\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{T}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}]$$

$$=\left(1-a\right)\int {\mid {U}_{T}\left(\mathit{x}\right)\mid}^{2}d\mathit{x}[\left(1-a\right)$$

$$+{\epsilon}_{0}\int \tilde{o}\left(\mathit{f}\right)\mathrm{OTF}\left(\mathit{f}\right){e}^{i2\pi \mathit{f}\xb7\mathit{x}}d\mathit{f}$$

$$+{\epsilon}_{0}^{*}\int {\tilde{o}}^{*}\left(\mathit{f}\right){\mathrm{OTF}}^{*}\left(\mathit{f}\right){e}^{-i2\pi \mathit{f}\prime \mathit{x}}d\mathit{f}],$$

*where ∫ U
_{R}
^{2} (x′-x)dx = ∫U
_{T}
^{2} (x′-x)dx = 0 and õ(f) is the Fourier transform of o(x). The 3-D optical transfer function is defined as*

*$$\mathrm{OTF}\left(\mathit{f}\right)=\frac{\int \{{\mid {U}_{T}\left(\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\right)\}{e}^{-i2\pi \mathit{f}\xb7\mathit{x}}d\mathit{x}}{\int {\mid {U}_{T}\left(\mathit{x}\right)\mid}^{2}d\mathit{x}},$$*

*which is equivalent to the correlation function between the entrance pupil P(f) = P
_{T}(f) + P
_{R}(f) and the pupil on transmission side P
_{T}(f) for the pump wave, that is*

*$$\mathrm{OTF}\left(\mathit{f}\right)=\frac{\int P\left(\mathit{f}\prime \right){P}_{T}^{*}\left(\mathit{f}\prime -\mathit{f}\right)d\mathit{f}\prime}{\int {\mid {P}_{T}\left(\mathit{f}\prime \right)\mid}^{2}d\mathit{f}\prime}.\text{}$$*

*Figure 7 shows the OTF in the case for NA = 1.2 with the primary objective in water. The OTF has rotational symmetry in the f_{z} direction, which is the spatial frequency in the direction of the optical axis. A cross section involving the f_{z} axis is described in the figure. Note that spatial frequencies along the optical axis cannot be resolved for both the transmitted and reflected components, since the scattered wave propagated in the vicinity of the 0-order light is intercepted by the micromirror, which is similar to the conventional phase contrast microscopy. The depth resolution is gained from the reflected component, and a part of the OTF corresponding to the reflected component lies in the region known as the missing cone. Although the gap between the two portions of the OTF corresponding to the transmitted and reflected components exists and the spatial frequency in the gap cannot be resolved, it can be reduced by using higher NA objectives.*

*If the amplitude of the object is a real-valued function O(x) = O*(x), which holds for a 3-D phase object with a low contrast refractive index distribution, Eq. (32) becomes*

*$${I}_{I}\left(\mathit{x}\prime \right)={\left(1-a\right)}^{2}\int {\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}d\mathit{x}+\left(1-a\right)\int {\epsilon}_{0}o\left(\mathit{x}\right)\left\{2{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)+{U}_{T}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}$$*

$$=\left(1-a\right)\lfloor \int \left\{\left(1-a\right)+2{\epsilon}_{0}o\left(\mathit{x}\right)\right\}\mathrm{PSF}\left(\mathit{x}\prime -\mathit{x}\right)/2d\mathit{x}\rfloor ,$$

$$=\left(1-a\right)\lfloor \int \left\{\left(1-a\right)+2{\epsilon}_{0}o\left(\mathit{x}\right)\right\}\mathrm{PSF}\left(\mathit{x}\prime -\mathit{x}\right)/2d\mathit{x}\rfloor ,$$

*where*

*$$\mathrm{PSF}\left(\mathit{x}\prime -\mathit{x}\right)=2{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)+{U}_{T}^{2}\left(\mathit{x}\prime -\mathit{x}\right)}^{\text{}}$$*

$$={\mid U\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}$$

$$={\mid U\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}$$

*is the 3-D point spread function. Interestingly, intensity of the low contrast object ∣ O(x)∣^{2} = 1 + 2ε_{0}
o(x) is effectively converted into ∣O
_{a}(x)∣^{2} = (1 - a) + 2ε
_{0}
o(x), which means that the initial contrast of the object _{2ε0} is enhanced to _{2ε0/(1-a)} by attenuation of the 0-order light. Finally, a simple equation for object reconstruction is obtained,*

*$${I}_{I}\left(\mathit{x}\prime \right)=\left(1-a\right)\{\int {\mid {O}_{a}\left(\mathit{x}\right)\mid}^{2}\mathrm{PSF}\left(\mathit{x}\prime -\mathit{x}\right)/2d\mathit{x}\}.$$*

*The irradiance of the image reconstructed through this algorithm is represented by the convolution of the object intensity, which has the enhanced contrast, with the 3-D point spread function. This method requires a relatively long computing time, because the 3-D Fourier transform is calculated in the computer for every incident direction of the pump wave.*

*In the second method, the irradiance of the image is given by*

*$${I}_{\mathrm{II}}\left(\mathit{x}\prime \right)={\mid \int \{\int F\left(\mathit{f},{\mathit{f}}_{0}\right)d{\mathit{f}}_{0}\}{e}^{i2\pi \mathit{f}\xb7\mathit{x}\prime}d\mathit{f}\mid}^{2}$$*

$$\phantom{\rule{.5em}{0ex}}={\mid \int [\int \left\{\int O\left(\mathit{x}\right){e}^{-i2\pi \mathit{f}\xb7\mathit{x}\prime}d\mathit{x}-\mathrm{a\delta}\left(\mathit{f}\right)\right\}\left\{{P}_{T}\left(\mathit{f}+{\mathit{f}}_{0}\right)+{P}_{R}\left(\mathit{f}+{\mathit{f}}_{0}\right)\right\}{P}_{T}^{*}\left({\mathit{f}}_{0}\right)d{\mathit{f}}_{0}]{e}^{i2\pi \mathit{f}\xb7\mathit{x}\prime}d\mathit{f}\mid}^{2}$$

$$\phantom{\rule{.5em}{0ex}}={\mid \int O\left(\mathit{x}\right){U}_{R}\left(\mathit{x}\prime -\mathit{x}\right)\left\{{U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)+{U}_{R}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}-a\int {\mid {P}_{T}\left({\mathit{f}}_{0}\right)\mid}^{2}d{\mathit{f}}_{0}\mid}^{2}$$

$$\phantom{\rule{.5em}{0ex}}={\mid \int O\left(\mathit{x}\right)\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}-a\int {\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}d\mathit{x}\mid}^{2},$$

$$\phantom{\rule{.5em}{0ex}}={\mid \int [\int \left\{\int O\left(\mathit{x}\right){e}^{-i2\pi \mathit{f}\xb7\mathit{x}\prime}d\mathit{x}-\mathrm{a\delta}\left(\mathit{f}\right)\right\}\left\{{P}_{T}\left(\mathit{f}+{\mathit{f}}_{0}\right)+{P}_{R}\left(\mathit{f}+{\mathit{f}}_{0}\right)\right\}{P}_{T}^{*}\left({\mathit{f}}_{0}\right)d{\mathit{f}}_{0}]{e}^{i2\pi \mathit{f}\xb7\mathit{x}\prime}d\mathit{f}\mid}^{2}$$

$$\phantom{\rule{.5em}{0ex}}={\mid \int O\left(\mathit{x}\right){U}_{R}\left(\mathit{x}\prime -\mathit{x}\right)\left\{{U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)+{U}_{R}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}-a\int {\mid {P}_{T}\left({\mathit{f}}_{0}\right)\mid}^{2}d{\mathit{f}}_{0}\mid}^{2}$$

$$\phantom{\rule{.5em}{0ex}}={\mid \int O\left(\mathit{x}\right)\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}-a\int {\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}d\mathit{x}\mid}^{2},$$

*where Parseval’s theorem is used. If the sample is a low contrast object, substituting O(x) = 1 + ε
_{0}
o(x) into Eq. (38) yields*

*$${I}_{\mathrm{II}}\left(\mathit{x}\prime \right)={\mid \int \left\{1+{\epsilon}_{0}o\left(\mathit{x}\right)\right\}\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}-a\int {\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}d\mathit{x}\mid}^{2}$$*

$$\phantom{\rule{.5em}{0ex}}=A\left(1-a\right)[\left(1-a\right)\int {\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}d\mathit{x}$$

$$\phantom{\rule{1em}{0ex}}+\int {\epsilon}_{0}o\left(\mathit{x}\right)\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}$$

$$\phantom{\rule{1em}{0ex}}+\int {\epsilon}_{0}^{*}{o}^{*}\left(\mathit{x}\right)\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{T}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}]$$

$$\phantom{\rule{.5em}{0ex}}=A{I}_{I}\left(\mathit{x}\prime \right)$$

$$\phantom{\rule{.5em}{0ex}}=A\left(1-a\right)[\left(1-a\right)\int {\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}d\mathit{x}$$

$$\phantom{\rule{1em}{0ex}}+\int {\epsilon}_{0}o\left(\mathit{x}\right)\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{R}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}$$

$$\phantom{\rule{1em}{0ex}}+\int {\epsilon}_{0}^{*}{o}^{*}\left(\mathit{x}\right)\left\{{\mid {U}_{T}\left(\mathit{x}\prime -\mathit{x}\right)\mid}^{2}+{U}_{T}^{2}\left(\mathit{x}\prime -\mathit{x}\right)\right\}d\mathit{x}]$$

$$\phantom{\rule{.5em}{0ex}}=A{I}_{I}\left(\mathit{x}\prime \right)$$

*where A= ∫∣U
_{T}(x′-x)\2 dx and the second order terms of ε
_{0} are ignored. The image irradiance of the reconstructed object through the second method is proportional to that of the first algorithm, as long as the sample is a low contrast object. The second method can reduce the computing time, because it requires only one 3-D Fourier transform in the final stage. The images of the low contrast object reconstructed by the two methods show identical optical features.*

*The image can be deconvoluted with the OTF in the same way as a conventional optical system. In the first method, each shifted partial sphere is divided by the OTF before it is Fourier transformed. In the second method, the deconvolution is achieved by dividing the total amplitude in the frequency domain by the OTF before taking the square of the modulus. Both methods are phase error free in the transmitted component, so low-frequency objects can be resolved with almost no error in reconstruction. While the relative phase shifts between different reflected holograms can be minimized by analyzing the phase in the overlap areas, a slight dc phase difference between the transmitted and reflected components remains because of the looped part.*

*6. Error estimation of object reconstruction*

*In the looped part, the optical path lengths between the centers of the sample holder and PBS1 on the transmission and reflection sides are designed to be equivalent. The optical path lengths for the transmitted and reflected diffraction waves between the micromirror and the digital camera are also designed to be identical. As mentioned above, any errors in phase due to changes in relative optical length among the three paths for the transmitted and reflected diffraction waves and the reference wave between the micromirror and the digital camera are corrected in the computer. However, there remains another phase error for the reflected component due to the change in relative optical length between the transmission and reflection sides in the looped part. Even though extremely steady and stable mechanical design is achieved for this part, the change in optical path lengths could occur during scanning and every hologram could have a slightly different dc phase component.*

*Error in object reconstruction is now evaluated, assuming that the phase error of each hologram is a Gaussian-distributed random pattern whose statistical profile is (mean value, standard deviation) = (0, σ). The root mean square (RMS) error is given by*

*$$\mathrm{RMS}\left(\sigma \right)=\sqrt{\sum _{i=1}^{N}\frac{{\left({I}_{i}^{\sigma}-{I}_{i}^{0}\right)}^{2}}{N},}$$*

*where N represents the number of the elements in the 3-D matrix of the image, I_{i}^{σ} denotes the value of the i-th element in the image with phase error σ, and I_{i}^{0} is the value of the i-th element in the image with no error. Figure 8 shows an example of the calculation result of the normalized RMS for a cubic object of 3.0μ m in size. RMS is evaluated for σ between zero and 2π×0.65. Each data point in Fig. 8 is an average of five hundred images, where a different Gaussian white noise distribution is used for each image. The calculation is performed with deconvoluted images. The object is well reconstructed as long as the standard deviation of the error is less than 0.2λ, which can be achieved by an ordinary mechanical design. Since the phase error is only in the reflected component and the part of the OTF corresponding to the reflected component is located apart from the origin, objects composed of only low spatial frequencies can be resolved almost perfectly. Different object distributions were also calculated, and it turns out that the configuration of the object does not strongly affect the RMS. All results have the tendency to produce a similar RMS versus σ relationship.*

*7. Discussion*

*The setup of the system composed only of the transmitted component (T-type) is simple and phase error free, resulting in being more practical. While this T-type microscopy has a sufficient capability and offers a high performance, the entire system composed of transmitted and reflected components (TR-type) yields a higher optical resolution, particularly in the depth direction. For 3-D specimens, as is the case in biological studies, high 3-D resolution is often required to visualize the 3-D image at any viewing angle. T-type microscopy, however, is still functional for the 3-D specimens which are composed of low spatial frequencies.*

*In the practical experiment, a CCD (Charge Coupled Device Image Sensor) or a CMOS (Complementary Metal Oxide Semiconductor Image Sensor) can be used as the digital camera. One of the most reasonable sizes N
^{2} of the digital camera is 1024×1024 pixels for each of the transmitted and reflected components. In this case, the size of the 3-D matrix (N/4)^{3} for the image, which is Fourier transformed three-dimensionally in the computer, is 256^{3} elements. The size of the reconstructed object in real space D/4 is approximately 45μm (the resolution X/2NA is 354nm) if a pump beam of 850nm in wavelength is used as the light source and NA = 1.2.*

*While it is necessary to process each axial sectioning separately and sequentially with the confocal microscopy in order to acquire the 3-D image, holographic microscopy proposed here requires only two-dimensional scanning to reconstruct the 3-D specimen. The total time to acquire a 3-D image through this system is dependent on the number of scanning directions. Image acquisition time will be at least a couple of minutes using current technology for the detector and computer. If the size of the 3-D matrix for the image is 256 ^{3}, 128^{2} × π/4 scan directions are required for maximum resolution, where the coefficient π/4 implies the circular pupil. However, the number of scanning directions can be reduced by balancing the quality of the image with the scanning time.*

*Near-infrared rays, for example the wavelength of 850nm, are used as the light source in this system because most biological specimens are transparent in that wavelength region. Then, the specimen can be considered to be a phase object that can be reconstructed through this system. Even though some absorbing materials exist in the specimen, an image can be obtained that contains the information on both an absorbing part and the phase part. The image only of the phase part is provided by adding the phase shift of π/2 to the 0-order light, and an image only of the absorbing part is obtained separately from the phase part by not adding the phase shift. While the absorbing object can be resolved by other microscopy techniques, such as transmission confocal microscopy, one of the advantages of this system is to be able to visualize 3-D phase objects. As mentioned above, the first Born approximation is assumed in the generation of the scattered wave through the interaction between the object and the pump wave, which implies that the specimen must be considered to be a low contrast object. If the specimen has a high contrast refractive index distribution, the object cannot be well-reconstructed, and thus the system has this applicative limitation.*

*8. Conclusion*

*A novel 3-D holographic microscopy is proposed that can reconstruct low contrast phase objects and convert them into high contrast images. Two algorithms for object reconstruction are given. It is deduced that these two methods show identical optical features of the image if the object has a low contrast structure. The second algorithm is more useful because of its shorter computing time. Mechanical tolerance is also evaluated, which indicates the upper limit of the allowable error in phase during the scanning. The whole system, including the mechanical design, can be built using current technology. Further advances in two-dimensional detector and computer power will be able to improve the quality of the image and reduce the computing time.*

*References and links*

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

**2. **W. B. Amos, J. G. White, and M. Fordham, “Use of confocal imaging in the study of biological structures,” Appl. Opt. **26**, 3239 (1987). [CrossRef] [PubMed]

**3. **G. J. Brakenhoff, H. T. M. van der Voort, E. A. van Spronsen, and N. Nanninga, “3-Dimensional imaging of biological structures by high resolution confocal scanning laser microscopy,” Scanning Microsc. **2**, 33 (1988). [PubMed]

**4. **I. Freund and M. Deutsch, “2nd-harmonic microscopy of biological tissue,” Opt. Lett. **11**, 94 (1986). [CrossRef] [PubMed]

**5. **P. J. Campagnola, H. A. Clark, W. A. Mohler, A. Lewis, and L. M. Loew, “Second-harmonic imaging microscopy of living cells,” J. Biomed. Opt. **6**, 277 (2001). [CrossRef] [PubMed]

**6. **J. Mertz and L. Moreaux, “Second-harmonic generation by focused excitation of inhomogeneously distributed scatterers,” Opt. Commun. **196**, 25 (2001). [CrossRef]

**7. **Y. Barad, H. Eisenberg, M. Horowitz, and Y. Silberberg, “Nonlinear scanning laser microscopy by third-harmonic generation,” Appl. Phys. Lett. **70**, 922 (1997). [CrossRef]

**8. **M. Muller, J. Squier, K. R. Wilson, and G. J. Brakenhoff, “3D microscopy of transparent objects using third-harmonic generation,” J. Microsc. **191**, 266 (1998). [CrossRef] [PubMed]

**9. **M. D. Duncan, J. Reintjes, and T. J. Manuccia, “Scanning coherent anti-Stokes Raman microscope,” Opt. Lett. **7**, 350 (1982). [CrossRef] [PubMed]

**10. **A. Zumbusch, G. R. Holtom, and X. S. Xie, “Vibrational microscopy using coherent anti-Stokes Raman scattering,” Phys. Rev. Lett. **82**, 4014 (1999).

**11. **F. Zernike, “Das Phasenkontrastverfahren bei der mikroskopischen Beobachtung,“ Z. Tech. Phys. **16**, 454 (1935).

**12. **F. Zernike, “How I discovered phase contrast,” Science **121**, 345 (1955). [CrossRef] [PubMed]

**13. **W. S. Haddad, D. Cullen, J. C. Solem, J. W. Longworth, A. McPherson, K. Boyer, and C. K. Rhodes, “Fourier-transform holographic microscope,” Appl. Opt. **31**, 4973 (1992). [CrossRef] [PubMed]

**14. **U. Schnars and W. Jüptner, “Direct recording of holograms by a CCD target and numerical reconstruction,” Appl. Opt. **33**, 179 (1994). [CrossRef] [PubMed]

**15. **J. H. Massig, “Digital off-axis holography with a synthetic aperture,” Opt. Lett. **27**, 2179 (2002). [CrossRef]

**16. **S. Kostianovski, S. G. Lipson, and E. N. Ribak, “Interference microscopy and Fourier fringe analysis applied to measuring the spatial refractive-index distribution,” Appl. Opt. **32**4744 (1993). [CrossRef] [PubMed]

**17. **D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang,, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography,” Science , **254**1178 (1991). [CrossRef] [PubMed]

**18. **T. Dresel, G. Hausler, and H. Venzke, “Three-dimensional sensing of rough surfaces by coherence radar,” Appl. Opt. **31**919 (1992). [CrossRef] [PubMed]

**19. **M. Mansuripur, *Classical Optics and its Applications* (Cambridge University Press2002)

**20. **M. Born and E. Wolf, *Principles of Optics* 5th. ed., (Pergamon Press, 1974).

**21. **H. Kogelnik, “Coupled wave theory for thick hologram gratings,“ Bell Syst. Tech. J , **48**2909 (1969).