In a number of previous studies on light focusing, the asymmetric axial intensity distribution with intensity peak shifted away from the paraxial focal plane was demonstrated for lenses working with a low Fresnel number. Here, the axial asymmetry of the three-dimensional point spread function (PSF) and the aberration effects are examined in a magnified phase-shifting holographic imaging achieved by the mismatch of reference and reconstruction waves. In the analysis, an optimal combination of experimental parameters and the range of applicable lateral magnifications are found for which the axial asymmetry of the PSF is not apparent and the aberration effects are acceptable. The focal shift and the axial asymmetry of the PSF and the effects of holographic aberrations are evaluated by approximate quantitative criteria whose validity is verified in exact numerical models and experiments. The optimal design of in-line holographic geometry is demonstrated by reconstructing the three-dimensional PSFs and the image of the resolution target recorded in the experimental setup using a spatial light modulator.
© 2017 Optical Society of America
In optical imaging theory, the basic imaging functions and criteria for evaluation of the imaging performance were preferably defined for two-dimensional objects. The three-dimensional (3D) imaging of volume samples became possible after the invention of confocal microscopy , and it has been further developed in modern imaging techniques including digital holographic microscopy  or super-resolution localization microscopy [3,4]. Properties of the 3D imaging have been of particular importance also in laser manipulation  and particle tracking . The point spread function (PSF) representing the blurred image of a point object was traditionally used as an indicator of the optical performance of imaging systems. To determine the 3D PSF in lens imaging, diffraction of a convergent spherical wave bounded by lens aperture was solved. In exact treatment, Fresnel-Kirchhoff and Rayleigh-Sommerfeld diffraction theories were formulated and the paraxial and Debye approximations proposed for the investigation of optical imaging. To apply these approximations correctly, the specific imaging conditions must be carefully considered. In optical systems with a low numerical aperture, where apodization, polarization and aberration effects are not essential, the paraxial approximation is well applicable. In high-aperture optical systems, the Debye approximation is used successfully. It is based on the representation of the focused light field by a superposition of plane waves whose propagation vectors fall inside the cone given by the aperture angle. A paraxial elaboration of the Debye approximation has also been proposed. It provides the 3D PSF with a symmetrical distribution of the axial intensity. This result is correct for higher apertures but fails in systems with a low Fresnel number, where the axial intensity exhibits asymmetry and its maximum is shifted away from the geometric focus [7, 8]. These effects are not apparent in lens imaging implemented by optical systems with commonly used parameters but can significantly affect the holographic imaging. In optical systems with higher numerical apertures, the axial asymmetry caused by diffraction is suppressed but the attention must be devoted to optical aberrations deteriorating the image quality and reducing the image resolution .
Here we examine the 3D PSF reconstructed from phase-shifted point holograms recorded in the experimental setup for digital in-line holographic microscopy . The main attention is paid to the distortion of the axial profile of the PSF whose intensity peak is shifted out of the paraxial image plane. This undesirable effect is studied as a result of the diffraction or holographic aberrations and its connection with the geometric parameters of experiments and the lateral image magnification is clarified. To implement quantitative criteria for the assessment of the diffraction and aberration deterioration of the 3D PSF, suitable approximations are used whose validity is verified by precise numerical simulations and experiments. The analysis provides an optimal combination of experimental parameters and an applicable range of the lateral magnifications for which the diffraction and aberration deterioration of the PSF is maintained at an acceptable level. Theoretical findings are verified in experiments using a spatial light modulator (SLM) that changes the tested parameters flexibly and permits implementation of the phase-shifting.
2.1. Three-dimensional PSF in digital in-line holography
The performance of optical imaging realized by digital in-line holography is studied using Fig. 1(a), illustrating both recording and numerical reconstruction of a point hologram. Laser light of the wavelength λ is focused on a pinhole of few microns in diameter, which acts as a source from which a reference wave ψr emanates. This wave illuminates a point-like scatterer that is located just behind the pinhole and generates a divergent signal wave ψs. The interference pattern captured by the CCD represents a point hologram given as H = |ψs + ψr|2. To determine the 3D PSF, the hologram is illuminated by a virtual reconstruction wave and the image is obtained by calculating free propagation of light diffracted at the hologram. In the exact treatment, Fresnel-Kirchhoff and Rayleigh-Sommerfeld diffraction integrals are solved numerically to reconstruct the image spot . For small numerical apertures, the Fresnel approximation can be used to solve the diffraction integral with a sufficient accuracy .
In digital holography, the reconstruction is carried out using a virtual monochromatic wave ψc having the same wavelength as that used in the recording. The light diffracted on the hologram placed in the plane x, y can be described using the first Rayleigh-Sommerfeld diffraction integral . The complex amplitude of the image ψ evaluated at the lateral position xi, yi located in the distance zi from the hologram can be written as(2), the remaining terms cause undesirable blurring. In optical holography, this problem has been solved by off-axis arrangement ensuring sufficient spatial separation of the dc term and the holographic twin image from the reconstructed image. In digital in-line holography various techniques have been proposed, in which the dc term was subtracted , or both the dc term and the twin image were completely eliminated by processing three or more holograms with different phase shifts between signal and reference waves . In experiments, the required phase shifts were precisely implemented using a piezoelectric transducer , a spiral phase plate , geometric phase optical elements , a liquid crystal variable retarder  or a spatial light modulator [19–22]. Applying the phase-shifting technique , the complex hologram is created allowing to reconstruct either virtual or real image. Here the 3D PSF is examined for a virtual image reconstructed from the complex function t given as (1) is solved numerically for the complex point hologram (3) to obtain the 3D PSF characterizing the reconstructed image. The axial asymmetry and the focal shift of the diffraction-limited PSF and its deterioration caused by aberrations are studied in appropriate approximations introduced in recording and reconstruction of the holograms.
2.2. Focal shift and axial asymmetry of the diffraction-limited PSF
Description of the experiments implemented with a low numerical aperture can be successfully performed using the paraxial approximation. The signal, reference and reconstruction waves then become paraboloidal waves whose phase (4) is simplified to the form(1) the Fresnel approximation is used, while in the amplitude term a less accurate replacement |ri| ≈ |zi| is applied. In the calculations, the optical sign convention is respected, in which zs < 0, zr < 0, zc < 0 and zi < 0 are assigned to divergent signal, reference and reconstruction waves and a virtual reconstructed image, respectively. To simplify the discussion, the pinhole located on the optical axis is considered, |r⊥0r| = 0. Neglecting terms that do not affect the shape of the 3D PSF, the complex amplitude Eq. (1) using the complex function Eq. (3) can be rewritten as the Fourier transform of the hologram aperture function T multiplied by a quadratic phase term,
In order to examine the axial asymmetry of the PSF, the reconstruction of the image is performed at positions shifted away from the focal plane, zi = f + Δzi. When the hologram is bounded by a circular aperture of radius ρh, its aperture function is given by a circle function, T = circ(|r⊥|/ρh). In this case, the intensity distribution providing the 3D PSF must be calculated numerically. If the axial profile of the PSF is investigated for a point scatterer positioned on the optical axis (|R⊥| = 0), the diffraction integral has an analytical solution providing the axial intensity defined as I = |ψ (|R⊥| = 0, Δzi)|2. The normalized axial intensity that is unitary in the paraxial focal plane, I(Δzi = 0) = 1, then can be written asFig. 1(b). These effects are influenced by the geometry of experiments and for their quantification the Fresnel number can be used that is defined as Ni = |ni|. The axial intensity varies periodically and takes zero values when q = 2mπ, where m = ±1, ±2, ···. The asymmetry of the axial intensity profile is assessed by the positions of the nearest zero points and related to m = ±1 [Fig. 1(b)], 8], equation (11) must be determined numerically. The position of the intensity maximum is then given by Fig. 1(b). When examining the axial asymmetry of the PSF in the image space, the coefficient Qi can be used, which is defined as
The optical performance achieved in digital in-line holography is influenced by three basic geometric parameters, ρp, Δ and zr, which denote the pinhole radius and distances from the pinhole to the sample and CCD, respectively. The pixel size of the CCD, pCCD, determines the conditions of the discretization of the record and is also an important parameter affecting the image quality. The lateral image magnification is dependent on Δ and zr but is also influenced by the reconstruction geometry. If the reconstruction wave is exactly matched to the reference wave, zc = zr, an aberration-free image is reconstructed at the object position with the unitary lateral magnification, f = zs and β = 1 [Eqs. (7) and (8)]. To get the magnified image, the reconstruction wave originating from a virtual source placed at the distances zc ∈ (−∞, zr) from the hologram is used. For the specified parameters Δ and zr, the largest lateral magnification is achieved for a plane wave reconstruction (zc → ∞). The mismatch between the reference and reconstruction waves results in optical aberrations and diffraction effects causing the axial asymmetry of the PSF. Manifestations of these effects vary with the lateral image magnification and will be examined both theoretically and experimentally.
When the point hologram is recorded using a low numerical aperture, the diffraction-limited image is reconstructed. The radius of the central Airy spot is given as Δri = αλ/NAi, where α = 0.61 for the circular hologram aperture and NAi ≈ ρh/|f|. The lateral resolution in the object space then can be written as Δr = Δri/|β|, hence Δr being inversely proportional to NAi|β|. As this product is independent of zc, NAi|β| = ρh/(zr − Δ), the theoretical lateral resolution in the object space remains unchanged when the reconstruction wave is changed. If the hologram aperture is defined by the light spot created by the reference wave on the CCD, the radius ρh is determined by the pinhole radius ρp as ρh = αλ|zr|/ρp. The resolution in the object space can then be rewritten to the form
In this paper, the image deterioration is investigated for the plane wave reconstruction providing the highest lateral magnification for given Δ and zr. With the increasing distance between the pinhole and the sample Δ, both the object space resolution and the axial symmetry of the PSF are improved [Eqs. (15) and (16)], but the detection requirements increase strongly. Considering the Nyquist-Shannon sampling theorem for the reconstruction using zc → ∞, the largest specimen to pinhole distance Δ permissible for the pixel size pCCD is determined by the condition(15), the holographic aberrations and diffraction effects causing the axial asymmetry of the PSF must be included in the analysis.
2.3. Influence of holographic aberrations
The virtual image is reconstructed from complex function (3) determined by the product of complex amplitudes of signal, reference and reconstruction waves, whose amplitudes and phases are given by (4). If the hologram is reconstructed by a monochromatic wave of the same wavelength as that used in the recording, and the paraxial approximation operating with paraboloidal waves is applied, the phase of the reconstructed hologram matches the phase of the wave converging at the paraxial image point with the coordinates xi = βxs, yi = βys and zi = f. When the experiment is performed with parameters beyond the paraxial approximation, the hologram phase determined by the product of the complex amplitudes of the spherical waves no longer corresponds to the phase of the spherical wave converging at the paraxial image point. The differences that appear are holographic aberrations. The third-order holographic aberrations are described by Taylor expansion applied to the phase Φj given by Eq. (4), in which only first three terms are included . To simplify the discussion, the pinhole and the source of the reconstruction wave located on the optical axis are considered, xj = yj = 0, j = r, c. Following the procedure specified in Appendix, the third-order holographic aberrations are obtained in the formEq. (3) and the signal wave is fully recovered. In this case, the aberration-free image with the unitary magnification is obtained (zi = zs, β = 1, hence S = C = A = F = D = 0).
When optimizing the magnified lensless holographic imaging achieved by the mismatch of the reconstruction and reference waves, the focus shift and the axial asymmetry of the PSF caused by diffraction are taken into account together with the PSF deterioration due to the spherical aberration. The off-axis aberrations are then assessed to determine the acceptable field of view. The analysis is performed for the plane wave reconstruction (zc → ∞) providing the highest lateral magnification. If the spherical aberration is evaluated in the paraxial image plane (zi = f), the coefficient S can be written as(15) and result in the increased Fresnel number (16) that ensures elimination of the axial asymmetry and the focal shift of the diffraction-limited PSF. Unfortunately, since Δ is positive and zr negative, the spherical aberration increases significantly with larger Δ as shown in (20). An accurate assessment of the effect of the holographic spherical aberration can be done numerically but a good estimate is even possible with the Strehl ratio K accessible by simplified calculations . Assuming a circular aperture with the radius ρh and a holographic wave aberration W, the Strehl ratio is defined as (22) results in (23) and (24), the decrease in Strehl ratio caused by spherical aberration S is after optimal refocusing approximately 16× smaller than in the paraxial image plane. As an indicator of the high imaging performance, the Strehl ratio K ≥ 0.8 is applied resulting in the requirement of a weak holographic spherical aberration limited by S ≤ 6/k. To achieve the desired values of K in the optimal image plane, the condition for the combination of basic experimental parameters can be obtained using (20), (25) and (26) ensure that the reconstructed image is not deteriorated by the spherical aberration at an optimal image plane.
3. Results and discussion
The magnified lensless holographic imaging realized by the mismatch of the reconstruction and reference waves is significantly affected by the pinhole radius ρp and the distances zr and Δ determining the pinhole position relative to the CCD and the sample, respectively. By combining these parameters inappropriately, the axial asymmetry and the focal shift or the holographic aberrations deteriorate the 3D PSF. Theoretical predictions of these effects resulting in an optimal choice of geometrical parameters were verified in experiments exploiting versatility of a SLM.
3.1. Experimental setup
The optical setup used is shown in Fig. 2. The light beam from He-Ne laser (10 mW, 632.8 nm) is spatially filtered by the single-mode fiber (Thorlabs P1-630A-FC-2, mode field diameter 3.6 – 5.3 μm) and collimated by the lens L1 (achromatic doublet Thorlabs AC254-150-A, diameter 25.4 mm, focal length 150 mm). The created beam passes through a beam splitter BS towards the reflective SLM (Hamamatsu X10468, 800 × 600, pixel size 20 μm). By a random pixel selection [19, 20], two different phase functions are simultaneously displayed on the SLM generating converging and diverging lenses with the focal lengths ± fm. Three different constant phase shifts ϑj are added to the quadratic phase of the converging lens so the phase functions tj sequentially sent to the SLM are given as15,25,], a complex hologram corresponding to the virtual image was obtained and the 3D PSF reconstructed numerically. When reconstructing image of the resolution target, only one convergent lens was displayed on the SLM that created the reference wave in combination with lens L2. The signal wave was created by diffraction of the reference wave on the resolution target.
3.2. Experimental results
In preparing the measurement, a theoretical dependence of the axial asymmetry and the aberration deterioration of the PSF on geometric parameters of the setup was examined using established quantitative indicators. In all demonstrated results, the monochromatic plane wave was used to reconstruct the recorded holograms. In Fig. 3, the deterioration of the reconstructed image of a point object caused by spherical aberration is demonstrated by the color-coded Strehl ratio evaluated in the image plane obtained by an optimal refocusing. The changes in the Strehl ratio are evaluated for different pinhole settings zr relative to the CCD and different positions of the point object from the pinhole Δ. In Fig. 3(a), the pinhole radius ρp = 2.75 μm is used, and for monitored settings zr, the radius of the detected point hologram is determined by the size of the diffraction spot created by the reference wave, ρh = αλ|zr|/ρp. In this case, the image deterioration does not change significantly with zr and the unacceptable decrease in the Strehl ratio occurs for Δ > 4 mm. In Fig. 3(b), the Strehl ratio is evaluated for the constant radius of the hologram ρh = 5.6 mm, corresponding to the diffraction spot created by the pinhole at the closest position in front of the CCD.
The dependence of the Strehl ratio K on Δ is shown in Fig. 4 (black line) for the constant pinhole position zr = −43.5 mm and the pinhole radius ρp = 2.75 μm, together with the diffraction-limited resolution in the object space Δr (green line), the resolution limit given by sampling conditions (thin black line), the coefficient of asymmetry Q (red line) and the coefficient of spherical aberration S (blue line). The coefficient Q is evaluated in the object space using a varying longitudinal magnification to transform the longitudinal distances from the image to object space. The dependence of K on Δ corresponds to the cross section of the color-coded graph in Fig. 3(a) indicated by a white dashed line. The range of the evaluated distances between the pinhole and the sample Δ ∈ 〈0, 10〉 mm is limited by the sampling conditions determined for CCD QImaging Retiga 4000R. For the greatest value of Δ, the best lateral resolution Δr ≈ 3.3λ is obtained. The lateral resolution is demonstrated by the radius of the diffraction spot Δr determined in the object space. While Δr decreases linearly with increasing Δ, the asymmetry and aberration coefficients Q and S show just opposite changes. For small values of Δ, the coefficient Q rapidly decreases indicating a strong asymmetry and shift of the axial PSF profile. In this case, the spherical aberration given by the coefficient S is very low so the Strehl ratio K is almost unitary. The situation is just opposite for the large distances Δ. The axial symmetry of the PSF with the coefficient Q ≈ 1 is achieved, while the spherical aberration increases and the Strehl ratio drops sharply. Axial symmetry of the PSF and nearly aberration-free imaging are guaranteed if both the asymmetry coefficient and the Strehl ratio are close to unit value. With the conditions Q ≥ 0.8 and K ≥ 0.8, the range of applicable positions Δ can be determined for the fixed parameters zr and ρp defining three different areas in Fig. 4. The areas identified as I, II and III determine the parameters resulting in diffraction-limited imaging with strong axial asymmetry of the PSF, low-aberration imaging with the nearly symmetrical 3D PSF and imaging strongly deteriorated by the spherical aberration, respectively. For ρp = 2.75 μm and zr = −43.5 mm, the applicable area II is defined by the limit values of the sample to pinhole distance Δ = 0.17 mm and 3.5 mm. When the position of the sample is changed within this range, the lateral magnification varies from β = 257 to 12.5, respectively. The maximum usable value of Δ is determined by the decrease of the Strehl ratio to K = 0.8 and follows from Eq. (26). The minimum value of Δ is given by the decrease of the asymmetry coefficient to Q = 0.8 and must be determined numerically. The experiments carried out with the parameters belonging to the individual areas of Fig. 4 (marked by violet points) provide results which closely match the theoretical predictions.
Verification measurements were performed in the setup shown in Fig. 2, in which the light spots created at the focal points Fr and Fs replaced the role of the pinhole and the point object in the lensless holographic configuration in Fig. 1. The radius of the diffraction spot created at the focal point Fr was approximately equal to the pinhole radius used in the analysis presented in Fig. 4. The distance Δ was operatively altered by displaying various converging and diverging lenses on the SLM. Each point hologram was recorded repeatedly with three different phase shifts between the signal and reference waves (0, 2π/3, 4π/3) that were set by the SLM. By processing the records using the phase-shifting technique , a complex hologram was created allowing to reconstruct the PSF undisturbed by the nondiffracted light and the holographic twin image.
The results obtained are demonstrated in Figs. 5, 6 and 7. In all three cases, the same position of the reference focal spot zr = −43.5 mm was used and the plane wave reconstruction of the PSF was carried out using the point holograms recorded with the same radius ρh = 6.1 mm. In the individual cases, the distance between the reference and signal focal spots Fr and Fs was set to Δ = 0.1 mm, 1.5 mm and 4.9 mm, respectively. The settings were carried out by the SLM and resulted in different focal lengths of the reconstructed complex hologram. With the parameters used in the individual cases, nearly the same diffraction-limited lateral resolution was achieved in the object space, specified by the theoretical values Δr = 2.7 μm, 2.65 μm and 2.4 μm, respectively. If the axial PSF profile is calculated using the paraxial Debye approximation, a symmetrical shape described by the sinc function is obtained. The half-width of the longitudinal spot determined by the first zero point of the sinc function remains approximately the same for all three cases provided the axial PSFs are compared in the object space. The theoretical values obtained for the examined cases are Δz = 64 μm, 60 μm and 56 μm. Although the demonstrated image reconstructions provide nearly the same theoretical values of the lateral and longitudinal resolution in the paraxial Debye approximation, significant differences are found when the 3D PSFs are experimentally implemented and correctly evaluated. The 3D PSFs shown in Figs. 5, 6 and 7 can be assigned, according to the effects influencing their shape, to the areas I, II and III in Fig. 4. Significant differences in the reconstructed PSFs are due to the fact that nearly the same lateral and longitudinal resolutions were achieved in a different manner in each case.
The PSF illustrated in Fig. 5 was reconstructed from the holograms recorded with ρh = 6.1 mm, zr = −43.5 mm and Δ = 0.1 mm. The used parameters correspond to the left violet dot in Fig. 4. With these parameters, the paraxial lateral and longitudinal resolutions Δr = 2.7 μm and Δz = 64 μm were obtained by combining a very low numerical aperture of the reconstructed hologram (NAi = 0.0003) with the extremely large lateral and longitudinal magnifications (β = 435, γ = β2). As the numerical aperture was low, the image was not impaired by aberrations, but the diffraction effects caused a strong asymmetry of the axial PSF due to the low Fresnel number (Ni = 3.1). The asymmetric axial PSF (blue line) was reconstructed from experimentally acquired holograms by algorithms based on the Kirchhoff transform . To demonstrate a good consistency between experiment and theory, a numerical model was created in which both hologram recording and reconstruction were simulated with parameters used in the real experiment. Results obtained in the precise nonparaxial approach are illustrated by a red line while the use of the paraxial approximation corresponds to a black dashed line.
The PSF shown in Fig. 6 belongs to the area II in Fig. 4 and was obtained by reconstructing the records acquired with optimally selected parameters ρh = 6.1 mm, zr = −43.5 mm and Δ = 1.5 mm (middle violet dot in Fig. 4). In this case, the paraxial lateral and longitudinal resolutions Δr = 2.65 μm and Δz = 60 μm were achieved by a balanced combination of the hologram aperture (NAi = 0.005) and the lateral and longitudinal magnifications (β = 29, γ = β2). The Fresnel number Ni = 48.3 was high enough for the axial PSF to be nearly symmetric and the aberrations were so low that the imaging was still almost diffraction-limited. The combination of parameters used complied with the condition (25), hence the Strehl ratio K ≥ 0.8 was ensured in the optimal image plane.
The PSF in Fig. 7 belongs to the area III in Fig. 4 (ρh = 6.1 mm, zr = −43.5 mm, Δ = 4.9 mm, right violet dot in Fig. 4), and corresponds to the paraxial lateral and longitudinal resolutions Δr = 2.4 μm and Δz = 56 μm obtained with the high image space aperture (NAi = 0.017) and the low lateral and longitudinal magnifications (β = 8.9, γ = β2). The Fresnel number is extremely high (Ni = 172) so the diffraction image asymmetry does not appear, but the lateral and longitudinal PSF profiles are strongly impaired by the spherical aberration due to the high image space aperture. The experimental PSF (blue line) is again in a good agreement with the PSF obtained in the numerical model (red line) in which both the recording and reconstruction of the hologram were simulated in the nonparaxial approach.
The impact of the geometric parameters on the imaging performance was also investigated using a resolution target. Experiments were carried out in the setup shown in Fig. 2, in which the distance of the reference focus Fr from the CCD was maintained constant, zr = −46.4 mm. In the examined cases, the resolution target was placed in three different positions relative to the reference focus given by Δ = 0.7 mm, 3 mm and 10 mm, respectively. The reconstruction of acquired holograms was performed using the monochromatic plane wave. The various positions of the target resulted in significantly different combinations of the image space aperture and the magnification. The image space apertures change as a result of different distances in which the holograms are reconstructed. In the settings under consideration, the reconstruction distances and the lateral magnifications get values zi = −3029 mm, −671 mm and −169 mm, and β = 66.3, 15.5 and 4.6, respectively. The imaging performance of individual settings can be assessed using the object space resolution estimated as the radius of the Airy diffraction spot divided by the lateral magnification, Δr = 0.61λ|zi|/(ρhβ). Hence, the resolution in the object space is improved with reduced |zi|/β ratio for constant ρh. This ratio decreases in each of the examined cases, giving values |zi|/β = 45.7, 43.4 and 36.7, respectively.
In the reconstructions of the resolution target in Fig. 8, the trend predicting improvement of the object space resolution with the increasing Δ is not fully confirmed because the aberration effects occur in these settings, as previously shown by the PSF reconstructions. The experimental results shown in Fig. 8 were obtained with a positive resolution target in which the line width in the finest group was 2.5 μm. In each reconstruction, an enlarged cut of the finest target group is shown along with the cross profiles and the average visibility.
In the image reconstruction in Fig. 8(a), performed with the low image space aperture and the high lateral magnification, aberration effects are negligible, but the resolution is limited by the |zi|/β ratio, which is highest of all three cases evaluated. In addition, axial asymmetry of the PSF occurs with these parameters.
The target reconstruction shown in Fig. 8(b) represents an optimal combination of geometric parameters resulting in imaging only weakly affected by holographic aberrations and providing resolution approaching the resolution limit determined by the sampling conditions applied to the CCD used. The PSF assigned to this case provides almost perfect axial symmetry.
Although the parameters used for the target reconstruction shown in Fig. 8(c) give the smallest |zi|/β ratio, the best resolution is not achieved in this case. The reason for the lower image contrast is spherical aberration which deteriorates both the transverse and axial profiles of the PSF as previously demonstrated.
We present theoretical and experimental analysis of the lensless holographic microscopy focused on a diffraction asymmetry of the axial PSF, conditions of correct sampling of experimental data and effects of holographic aberrations that occur when the magnified image is created by the mismatch of reconstruction and reference waves. We show in the imaging using the plane wave reconstruction that the same theoretical resolution can be achieved with various combinations of experimental parameters having significantly different impact on the imaging performance. In extreme cases of the design, the same theoretical lateral and longitudinal resolutions are achieved by combining a very low numerical aperture of the reconstructed hologram with an extremely high lateral and longitudinal magnifications, or conversely, a high hologram aperture is used with low magnifications. While in the first case an extremely asymmetric axial PSF arises with its maximum shifted outside the paraxial focus, in the latter case the spherical aberration unacceptably deteriorates both lateral and longitudinal profiles of the PSF. In suitably designed experiments, a symmetrical, almost diffraction-limited PSF is obtained by a balanced combination of the numerical aperture of the reconstructed hologram and the magnification. The choice of the geometric parameters that ensure optimal imaging performance is governed by the conditions that are found and verified in this paper.
Effects examined theoretically are verified experimentally in the setup utilizing the versatility of a SLM. The experimental results presented by reconstructions of the 3D PSFs and the resolution target are in a good agreement with theoretical predictions.
To examine the holographic aberrations, the square roots in the phase of signal and reference waves (4) are expanded to Taylor series. In the same way, the phase of the wave corresponding to the paraxial image is expanded for the point of convergence (xi, yi, zi). Considering only the first three terms of the series and respecting the introduced sign convention, the phase of waves can be written as(7) and (8) can be obtained. Holographic wave aberrations defined as the optical path difference between the reference wavefront and the equiphase surface given by the complex hologram are determined by the terms of the third order in zj, (30) into (29) and using the normalized radial coordinate ρN = ρ/ρh, where ρh is the radius of the circular hologram aperture, the total holographic aberration given by Eq. (18) is obtained.
Grant Agency of the Czech Republic (No. 15-14612S); Palacký University in Olomouc (IGA-PrF-2017-002).
P. Bouchal from Brno University of Technology, the Czech Republic, is acknowledged for his contribution to the analysis of axially asymmetric holographic imaging made at the early stage of the research.
References and links
1. M. Gu, Principles of Three-Dimensional Imaging in Confocal Microscopes (World Scientific, 1996). [CrossRef]
2. P. Marquet, B. Rappaz, P.J. Magistretti, E. Cuche, Y. Emery, T. Colomb, and C. Depeursinge, “Digital holographic microscopy: a non invasive contrast imaging technique allowing quantitative visualization of living cells with sub wavelength accuracy,” Opt. Lett. 30(5), 468–470 (2005). [CrossRef] [PubMed]
3. G. Grover, K. DeLuca, S. Quirin, J. DeLuca, and R. Piestun, “Super-resolution photon-efficient imaging by nanometric double-helix point spread function localization of emitters (SPINDLE),” Opt. Express 20(24), 26681–26695 (2012). [CrossRef] [PubMed]
4. H. I. D. Lee, S. J. Sahl, M. D. Lew, and W. E. Moerner, “The double-helix microscope super-resolves extended biological structures by localizing single blinking molecules in three dimensions with nanoscale precision,” Appl. Phys. Lett. 100, 153701 (2012). [CrossRef] [PubMed]
5. K. Dholakia and T. Čižmár, “Shaping the future of manipulation,” Nature Photon. 5, 335–342 (2011). [CrossRef]
6. X. Yu, J. Hong, Ch. Liu, and M. K. Kim, “Review of digital holographic microscopy for three-dimensional profiling and tracking,” Opt. Eng. 53(11), 112306 (2011). [CrossRef]
7. J. H. Erkkila and M. E. Rogers, “Diffracted fields in the focal volume of a converging wave,” J. Opt. Soc. Am. A 71(7), 904–905 (1981). [CrossRef]
8. Y. Li and E. Wolf, “Focal shifts in diffracted converging spherical waves,” Opt. Commun. 39(4), 211–215 (1981). [CrossRef]
11. U. Schnars and W. Juptner, “Digital recording and numerical reconstruction of holograms,” Meas. Sci. Technol. 13(9) R85–R101 (2002). [CrossRef]
12. T. Kreis, Handbook of Holographic Interferometry: Optical and Digital Methods (Wiley-VCH Verlag, 2005), Chap. 4.
13. J. W. Goodman, Introduction to Fourier optics (McGraw-Hill, 1996).
18. P. Bouchal, R. Čelechovský, and Z. Bouchal, “Polarization sensitive phase-shifting Mirau interferometry using a liquid crystal variable retarder,” Opt. Lett. 40(19) 4567–4570 (2015). [CrossRef] [PubMed]
20. B. Katz, J. Rosen, R. Kelner, and G. Brooker, “Enhanced resolution and throughput of Fresnel incoherent correlation holography (FINCH) using dual diffractive lenses on a spatial light modulator (SLM),” Opt. Express 20(8) 9109–9121 (2012). [CrossRef] [PubMed]
23. W. Reinhard, “Magnification and Third-Order Aberrations in Holography,” J. Opt. Soc. Am. A 55(8), 987–992 (1965). [CrossRef]
24. V. N. Mahajan, “Strehl ratio for primary aberrations: some analytical results for circular and annular pupils,” J. Opt. Soc. Am. 72(9), 1258–1266 (1982). [CrossRef]
25. V. Micó, Z. Zalevsky, and J. Garcia, Superresolved common-path phase-shifting digital inline holographic microscopy using a spatial light modulator,” Opt. Lett. 37(23), 4988–4990 (2012). [CrossRef] [PubMed]
26. P. Ferraro, A. Wax, and Z. Zalevsky, Coherent Light Microscopy: Imaging and Quantitative Phase Analysis (Springer, 2011), Chap. 1.