Abstract
We examine the interplay between spectral bandwidth and illumination curvature in ptychography. By tailoring the divergence of the illumination, broader spectral bandwidths can be tolerated without requiring algorithmic modifications to the forward model. In particular, a strong wavefront curvature transitions a far-field diffraction geometry to an effectively near-field one, which is less affected by temporal coherence effects. The relaxed temporal coherence requirements allow for leveraging wider spectral bandwidths and larger illumination spots. Our findings open up new avenues towards utilizing pink and broadband beams for increased flux and throughput at both synchrotron facilities and lab-scale beamlines.
Published by Optica Publishing Group under the terms of the Creative Commons Attribution 4.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.
1. Introduction
Ptychography [1,2] is a scanning version of coherent diffractive imaging (CDI), in which the object is illuminated in overlapping positions by a confined illuminating beam. It has found applications in a wide spectral range, including terahertz [3], near-infrared (NIR) [4], visible spectrum (VIS) [5], extreme ultraviolet (XUV) [6], soft X-ray [7], and hard X-rays [2,8]. While single-shot CDI requires a highly coherent illumination [9], ptychography can be used to robustly account for and mitigate decoherence effects [10–12]. For both single-shot CDI and ptychography, algorithmic and experimental techniques have been developed to account for reduced temporal coherence. These include physical modeling of dispersion effects in free-space propagation [13–18], tailoring or accounting for wavelength-dependent wavefronts to retrieve multiple spectral components in the illumination [12,19–21], or using techniques from Fourier transform spectroscopy [22–24]. While the aforementioned techniques allow for utilizing broadband radiation in ptychographic beamlines, they come at the cost of either algorithmic overhead or interferometric hardware modifications. Consequently, the question arises under which conditions a polychromatic forward model or hardware modifications are needed at all. Early work on single-shot CDI suggested a trade-off between spectral bandwidth and spatial resolution [25]. This criterion has been utilized for ptychography [26] as a rule-of-thumb. Furthermore, recent work [27] has empirically demonstrated that ptychography exhibits a higher tolerance for broadband conditions in the near-field as compared to far-field (Fraunhofer) diffraction geometries. However, switching to a near-field geometry in the short-wavelength regime such as XUV or x-rays is challenging due to physical constraints that limit the minimum object to detector distance. In our work, we explore how beam curvature can serve as an additional degree of freedom that relaxes coherence requirements, thereby allowing for larger spectral bandwidths and potentially faster scanning and higher throughput in ptychography experiments utilizing sources with low temporal coherence.
The purpose of this work is threefold: First, We examine the spectral bandwidth criterion for single-shot CDI and investigate its validity for ptychography under varied experimental conditions and the assumption of a monochromatic forward model. Second, We continuously transition between near- and far-field diffraction geometries by controlling the illumination curvature – a mechanism that can be harnessed in short-wavelength experimental setups without any hardware modifications. Finally, We examine implications for dispersive specimens. The latter point is critical: does tolerance to broad bandwidth conditions imply insensitivity to dispersion effects? We demonstrate that tolerance to broad bandwidth conditions and information retrieval about dispersion are not mutually conflicting goals. Our findings are supported by both simulation and experiment.
The paper is structured as follows: In Section 2, we revisit the geometrical arguments that lead to the single-shot CDI spectral bandwidth versus spatial resolution restriction, and provide an alternative theoretical viewpoint, where we discuss the stationary phase approximation, which provides a mechanism to understand how temporal coherence requirements are relaxed under increasing illumination curvature. In Section 3, we present experimental results that demonstrate improved tolerance to broadband conditions under curved illumination conditions. The results are generalized by a simulation study on dispersive specimens in Section 4. Our findings are discussed and summarized in Section 5.
2. Broadband diffraction considerations
2.1 Geometrical considerations
In this section we revisit an inequality between the bandwidth and the achievable lateral resolution in a single-shot CDI setup, originally described in [25]. With regard to Fig. 1, consider two laterally displaced point sources on a specimen. When sub-waves originating from these two points interfere on a distant detector, the contrast tends to vanish as the optical path length is increased. In fact, writing down an explicit expression for the optical path length difference and using the paraxial approximation the authors of [25] estimated that when the temporal coherence length is larger than the optical path length difference, the following inequality holds true:
Here $\delta$ is the smallest resolvable feature, $D$ is the lateral size of the illumination, and $\Delta \lambda$ and $\lambda _c$ are the spectrum’s bandwidth and its center wavelength respectively. Inequality (1) predicts a linear behavior between the spectral bandwidth and the smallest resolvable feature, namely that we should expect a twofold loss in spatial resolution if we double the bandwidth in a single-shot CDI experiment.
Since the smallest resolvable feature $\delta$ is bounded by the wavelength and numerical aperture of the setup geometry, Eq. (1) can be written in terms of the numerical aperture NA of the system,
where the right term of the inequality $\lambda _c^2/\Delta \lambda$ is an estimation of the temporal coherence length of the source [28].2.2 Stationary phase approximation
In this section, we analyse the oscillatory behavior of the Fresnel diffraction integral under curved wavefront illumination and apply the method of stationary phase to explain its implications for temporal coherence requirements. Consider a two-dimensional electric field $E$ diffracted from an aperture (see Fig. 1). Under the paraxial approximation, the diffracted field at a distance $z$ can be expressed in the form of the following Fresnel diffraction integral [29]
3. Experiments
3.1 Setup
The experimental setup is shown in Fig. 3(a)). A super-continuum laser source was spectrally filtered by means of bandpass filters, resulting in four different center wavelengths ($\lambda _c=$[632.8 nm, 650 nm, 643.95 nm, 636.57 nm]) and 4 different relative bandwidth ($\Lambda = \Delta \lambda / \lambda _{c}$) conditions (0.18%, 0.95%, 3.3%, 8.1%). The beam was spatially filtered (L1-P1-L2) to ensure effects due to reduced spatial coherence can be ignored. To this end P1 was set to $\approx$30% larger than the diffraction-limited resolution of L1. A 500 µm pinhole (P2) was imaged onto the specimen using a 4f-system. The detector downstream of the object had a pixel-size of 3.45 µm $\times$ 3.45 µm, an array size of 4438 $\times$ 4438, and a 12-bit dynamic range. The divergence of the illumination was controlled by a movable lens (L3) upstream of the pinhole (P2), which resulted in an illumination NAillu=[0.005, 0.02, 0.027] for its corresponding curvature $f_i$=[5 cm, 1.3 cm, 0.9 cm]. An example of the raw data for $f=$0.9 cm showing the spectral blurring due to the increased spectral bandwidth is shown in Figs. 3(b),(c). Four different bandwidth conditions and three different phase curvatures were used [Fig. 4(a)-(c)]. The sample-detector distance was found to be 55.28 mm using the zPIE axial position calibration algorithm [33], which resulted in an NA$\approx$0.14 and a lateral resolution limit of 2.3 µm. To generate comparable photon flux conditions between all 12 combinations for curvatures and bandwidths, the total number of intensity counts for each dataset was held constant. This resulted in a lower signal to noise ratio for the case of the more divergent beam. Each dataset consisted of 100 scan positions arranged in a Fermat spiral [34] with an average linear overlap [35] of 90% between consecutive positions. A USAF resolution test target was chosen as a specimen to allow for comparing the lateral resolution between each reconstruction.
3.2 Results
The open-source software PtyLab [36] was used to run 300 iterations of mPIE [37] using the same initial parameters for the reconstruction of all datasets. The reconstruction results are shown in Fig. 4. For all reconstructions an inset of the same region was inserted for comparison. Two general tendencies are seen from the reconstructed micrographs: First, the lateral resolution decreases with increasing bandwidth (left to right columns). Second, the loss of lateral resolution can be compensated by increasing the wavefront curvature (bottom to top row). Thus our main finding is that a larger spectral bandwidth can be tolerated when the curvature of the illumination increases. Comparing, for instance, conditions a-$\Lambda _3$ with c-$\Lambda _2$, we see that a comparable imaging resolution is achieved, albeit using a 3 times larger relative bandwidth in the former. To evaluate the achieved resolution, we applied the Rayleigh criterion to a lineout taken across the smallest features of the reconstructed USAF target, specifically where the visibility of the features exceeded 20%. Next, we identified the group and number of these features and computed the full pitch resolution using the USAF target formula. An overview of the lateral resolution $\delta$ for each panel is summarized in Table 1 and Fig. 5. The columns and rows in Table 1 are organized in the same manner as in Fig. 4. For easy comparison, a relative resolution $\delta /D$ represented in percentage was added. The blue shaded cells mark where Eq. (1) holds, while the red cells indicate that this inequality was violated. Thus we find both circumstances where Eq. (1) is outperformed and other situations where the lower resolution bound is not reached.
4. Dispersive sample considerations
The last section showed that a beam with a curved wavefront can increase the tolerance to a broader spectrum, provided that we assume a monochromatic forward model. However, ptychography is well-known to allow for extracting polychromatic wavefronts as well as individual color channels of dispersive specimens [12,17,19–21]. Thus we must ask the question how curvature affects our ability to extract individual spectral components if we assume a polychromatic forward model. Does increasing wavefront curvature entail a loss of sensitivity to polychromatic information? While in the last section the specimen was non-dispersive and the wavefront for individual spectral components is not expected to change significantly, the situation would certainly be different if the specimen was dispersive. Under such circumstances we would want to be able to spectrally resolve the specimen [12,38] and extract the individual spectral wavefronts [17,20,21]. To this end, we simulated a ptychographic experiment using a dispersive sample and a trichromatic, discrete spectrum.
The simulated beam consisted of three spectral lines at wavelengths 600 nm, 650 nm, and 700 nm. For the dispersive sample, we used a Siemens star resolution test target with the refractive index values of NBK7 glass and a binary thickness profile of 632 nm that results in a phase delay of OPD $\approx \pi$ radians at $\lambda$ = 650 nm. We also included an amplitude spectral response from this sample with a custom coating layer, whose transmissivity increases towards longer wavelengths such that we expect 50% transmission at the central wavelength, as shown in Figs. 6(a)-(c). With the aforementioned parameters we expected a distinct object spectral response (amplitude and phase shift) for each spectral line, as depicted in Table 2. The setup parameters were chosen to be similar to our experimental setup (i.e., beam diameter D=500 µm, distance to detector z=55.28 mm, detector size $L$=12.8 mm, 200 scans positions with an average linear overlap of $\sim$ 80% between consecutive positions).
The simulations were carried twice, once with a curved beam ($f$=20 mm), and then with a flatter illumination ($f$=100 mm). The reconstructed results for both of these conditions are shown in Fig. 7. On the columns (a)-(c), (g)-(i) we see the object spectral response, where the color represent the phase shift and amplitude is encoded in the brightness. Next to each object is the beam profile for each spectral line plotted in the same color scale as the object. At first sight there is no clear difference between the curved and flat beam object reconstructions. However, to assess the differences closely, we analyzed the histogram of the reconstructed object’s amplitude and phase from each spectral line (see Fig. 8(a)), and fitted a Gaussian mixture model to evaluate the shape and location of the reconstructed amplitude and phase values against the expected values in Table 2. We found that the histogram’s peaks of both reconstructions closely match the ground truth spectral response (black triangles). However, the histogram from the data set with a flatter beam (bottom row) shows wider lobes in both the amplitude- and phase-distributions, indicating a larger uncertainty in the reconstruction for the flat beam. On the contrary, the histograms retrieved with a curved beam (top row) show no significant artifacts, which is important for quantitative results and material specific imaging [39]. In order to exclude spurious oscillations from difficult to resolve high spatial frequencies, these histograms were extracted from a region of interest mainly containing mid-spatial frequencies, as shown in Fig. 8(b). Thus the main finding from this section is that curved beams do not entail a loss of sensitivity to polychromatic information. As in the monochromatic case, they help improve the reconstruction quality also for multi-spectral wavefronts.
It should be noticed that for the experimental results were obtained with a broad continuous spectrum while the numerical simulations presented here were performed with three discrete wavelengths – which is a simplification of a typical spectrum provided with HHG sources [19,21]. By refraining from simulating continuum spectra we circumvent another open question in ptychography, namely: how many discrete spectral points are required to represent a continuum spectrum? Previous work has already shown the capabilities to use ptychography with such spectra and reconstruct it a certain number of colors [17,18,40]. However, the criteria for choosing n-number of wavelengths to represent the input spectrum is still not clear. Thus, we restricted the simulations to discrete colors, such that the reconstructed object spectral response was not influenced by the aforementioned interrogation.
5. Discussion and conclusion
In this work, we expanded upon the experimental observations of [27], and observed an improved tolerance of ptychography to a broad spectrum under a monochromatic forward model for increasingly divergent beams. Furthermore, we found that the commonly used temporal coherence limit for single shot CDI [25] can be outperformed in ptychography, which is consistent with previous work [17,26]. This is clear evidence that the spectral bandwidths constraints in ptychography are more relaxed than in single shot CDI, and that curved beams affect this tolerance. However, we are unable to provide a theoretical bound for a rigorous spectral bandwidth constraint that holds for ptychography. Further investigations in this direction will be required in future work.
We presented an intuitive explanation for the increased tolerance to broadband illumination under increasing wavefront curvature via the stationary phase approximation. In particular, the effective region of integration in the Fresnel diffraction integral shrinks as curvature increases. This translates into an effectively smaller beam, which explains why the simplified single shot CDI rule of thumb (Eq. (2), which assumes a fixed beam diameter), is not valid. In other words, larger optical path differences do not contribute to the integral and the temporal coherence requirement is relaxed.
We extended our investigation to polychromatic specimens and wavefronts. A naive extrapolation of the previous finding could suggest that curved beams are harmful for polychromatic experiments, as they would inhibit sensitivity to dispersion effects under a monochromatic forward model. However, our simulations indicate that that beam curvature does not result in a loss of spectral sensitivity, in spite of the previously mentioned improved tolerance under a monochromatic forward model. The simulations results in section 4. showed that a polychromatic ptychography experiment likewise benefits from a curved illumination.
In essence, larger beam curvatures allow for more relaxed spectral bandwidth requirements. We believe these finding have important bearings on the achievable throughput in photon-limited ptychography experiments, including inherently broadband tabletop high-harmonic generation [39] and x-ray sources [17].
Funding
Fraunhofer-Gesellschaft (Cluster of Excellence Advanced Photon Sources); Helmholtz Association (Ptychography 4.0, ZT-I-PF-4-018 (AsoftXm)); Thüringer Ministerium für Bildung, Wissenschaft und Kultur (2018 FGR 0080); European Social Fund.
Disclosures
The authors declare no conflicts of interest.
Data availability
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
References
1. J. M. Rodenburg, “Ptychography and related diffractive imaging methods,” Adv. Imaging Electron Phys. 150, 87–184 (2008). [CrossRef]
2. P. Thibault, M. Dierolf, A. Menzel, O. Bunk, C. David, and F. Pfeiffer, “High-resolution scanning x-ray diffraction microscopy,” Science 321(5887), 379–382 (2008). [CrossRef]
3. L. Valzania, T. Feurer, P. Zolliker, and E. Hack, “Terahertz ptychography,” Opt. Lett. 43(3), 543–546 (2018). [CrossRef]
4. M. Du, L. Loetgering, K. S. E. Eikema, and S. Witte, “Measuring laser beam quality, wavefronts, and lens aberrations using ptychography,” Opt. Express 28(4), 5022–5034 (2020). [CrossRef]
5. T. Wang, S. Jiang, P. Song, R. Wang, L. Yang, T. Zhang, and G. Zheng, “Optical ptychography for biomedical imaging: recent progress and future directions,” Biomed. Opt. Express 14(2), 489–532 (2023). [CrossRef]
6. L. Loetgering, S. Witte, and J. Rothhardt, “Advances in laboratory-scale ptychography using high harmonic sources,” Opt. Express 30(3), 4133–4164 (2022). [CrossRef]
7. K. Giewekemeyer, P. Thibault, S. Kalbfleisch, A. Beerlink, C. Kewish, M. Dierolf, F. Pfeiffer, and T. Salditt, “Quantitative biological imaging by ptychographic x-ray diffraction microscopy,” Proc. Natl. Acad. Sci. 107(2), 529–534 (2010). [CrossRef]
8. J. Rodenburg, A. Hurst, A. Cullis, B. Dobson, F. Pfeiffer, O. Bunk, C. David, K. Jefimovs, and I. Johnson, “Hard-x-ray lensless imaging of extended objects,” Phys. Rev. Lett. 98(3), 034801 (2007). [CrossRef]
9. J. Spence, U. Weierstall, and M. Howells, “Coherence and sampling requirements for diffractive imaging,” Ultramicroscopy 101(2-4), 149–152 (2004). [CrossRef]
10. P. Thibault and A. Menzel, “Reconstructing state mixtures from diffraction measurements,” Nature 494(7435), 68–71 (2013). [CrossRef]
11. R. Röhrich, A. F. Koenderink, S. Witte, and L. Loetgering, “Spatial coherence control and analysis via micromirror-based mixed-state ptychography,” New J. Phys. 23(5), 053016 (2021). [CrossRef]
12. D. J. Batey, D. Claus, and J. M. Rodenburg, “Information multiplexing in ptychography,” Ultramicroscopy 138, 13–21 (2014). [CrossRef]
13. R. A. Dilanian, B. Chen, G. J. Williams, H. M. Quiney, K. A. Nugent, S. Teichmann, P. Hannaford, L. V. Dao, and A. G. Peele, “Diffractive imaging using a polychromatic high-harmonic generation soft-x-ray source,” J. Appl. Phys. 106(2), 023110 (2009). [CrossRef]
14. B. Abbey, L. W. Whitehead, H. M. Quiney, D. J. Vine, G. A. Cadenazzi, C. A. Henderson, K. A. Nugent, E. Balaur, C. T. Putkunz, A. G. Peele, G. J. Williams, and I. McNulty, “Lensless imaging using broadband x-ray sources,” Nat. Photonics 5(7), 420–424 (2011). [CrossRef]
15. J. Huijts, S. Fernandez, D. Gauthier, M. Kholodtsova, A. Maghraoui, K. Medjoubi, A. Somogyi, W. Boutu, and H. Merdji, “Broadband coherent diffractive imaging,” Nat. Photonics 14(10), 618–622 (2020). [CrossRef]
16. W. Eschen, S. Wang, C. Liu, R. Klas, M. Steinert, S. Yulin, H. Meißner, M. Bussmann, T. Pertsch, J. Limpert, and J. Rothhardt, “Towards attosecond imaging at the nanoscale using broadband holography-assisted coherent imaging in the extreme ultraviolet,” Commun. Phys. 4(1), 154 (2021). [CrossRef]
17. Y. Yao, Y. Jiang, J. Klug, Y. Nashed, C. Roehrig, C. Preissner, F. Marin, M. Wojcik, O. Cossairt, Z. Cai, S. Vogt, B. Lai, and J. Deng, “Broadband x-ray ptychography using multi-wavelength algorithm,” J. Synchrotron Radiat. 28(1), 309–317 (2021). [CrossRef]
18. D. Goldberger, D. Schmidt, J. Barolak, B. Ivanic, C. G. Durfee, and D. E. Adams, “Spatiospectral characterization of ultrafast pulse-beams by multiplexed broadband ptychography,” Opt. Express 29(20), 32474–32490 (2021). [CrossRef]
19. B. Zhang, D. F. Gardner, M. H. Seaberg, E. R. Shanblatt, C. L. Porter, R. Karl, C. A. Mancuso, H. C. Kapteyn, M. M. Murnane, and D. E. Adams, “Ptychographic hyperspectral spectromicroscopy with an extreme ultraviolet high harmonic comb,” Opt. Express 24(16), 18745–18754 (2016). [CrossRef]
20. L. Loetgering, X. Liu, A. C. De Beurs, M. Du, G. Kuijper, K. S. Eikema, and S. Witte, “Tailoring spatial entropy in extreme ultraviolet focused beams for multispectral ptychography,” Optica 8(2), 130–138 (2021). [CrossRef]
21. M. Du, X. Liu, A. Pelekanidis, F. Zhang, L. Loetgering, P. Konold, C. L. Porter, P. Smorenburg, K. S. Eikema, and S. Witte, “High-resolution wavefront sensing and aberration analysis of multi-spectral extreme ultraviolet beams,” Optica 10(2), 255–263 (2023). [CrossRef]
22. S. Witte, V. T. Tenner, D. W. Noom, and K. S. Eikema, “Lensless diffractive imaging with ultra-broadband table-top sources: from infrared to extreme-ultraviolet wavelengths,” Light: Sci. Appl. 3(3), e163 (2014). [CrossRef]
23. G. Jansen, D. Rudolf, L. Freisem, K. Eikema, and S. Witte, “Spatially resolved fourier transform spectroscopy in the extreme ultraviolet,” Optica 3(10), 1122–1125 (2016). [CrossRef]
24. M. Du, L. Loetgering, K. S. Eikema, and S. Witte, “Ptychographic optical coherence tomography,” Opt. Lett. 46(6), 1337–1340 (2021). [CrossRef]
25. F. Van der Veen and F. Pfeiffer, “Coherent x-ray scattering,” J. Phys.: Condens. Matter 16, 5003 (2004). [CrossRef]
26. B. Enders, M. Dierolf, P. Cloetens, M. Stockmar, F. Pfeiffer, and P. Thibault, “Ptychography with broad-bandwidth radiation,” Appl. Phys. Lett. 104(17), 171104 (2014). [CrossRef]
27. W. Xu, S. Ning, and F. Zhang, “Numerical and experimental study of partial coherence for near-field and far-field ptychography,” Opt. Express 29(24), 40652–40667 (2021). [CrossRef]
28. J. W. Goodman, Statistical Optics (John Wiley & Sons, 2015).
29. J. W. Goodman, Introduction to Fourier Optics (McGraw-Hill, 1996).
30. G. J. Gbur, Mathematical Methods for Optical Physics and Engineering (Cambridge University, 2011).
31. M. Guizar-Sicairos, M. Holler, A. Diaz, J. Vila-Comamala, O. Bunk, and A. Menzel, “Role of the illumination spatial-frequency spectrum for ptychography,” Phys. Rev. B 86(10), 100103 (2012). [CrossRef]
32. X. Ji, X. He, Z. Jiang, Y. Kong, S. Wang, and C. Liu, “Resolution enhancement with highly curved illumination in ptychography,” Appl. Opt. 61(34), 10150–10158 (2022). [CrossRef]
33. L. Loetgering, M. Du, K. S. Eikema, and S. Witte, “zpie: an autofocusing algorithm for ptychography,” Opt. Lett. 45(7), 2030–2033 (2020). [CrossRef]
34. X. Huang, H. Yan, R. Harder, Y. Hwu, I. K. Robinson, and Y. S. Chu, “Optimization of overlap uniformness for ptychography,” Opt. Express 22(10), 12634–12644 (2014). [CrossRef]
35. O. Bunk, M. Dierolf, S. Kynde, I. Johnson, O. Marti, and F. Pfeiffer, “Influence of the overlap parameter on the convergence of the ptychographical iterative engine,” Ultramicroscopy 108(5), 481–487 (2008). [CrossRef]
36. L. Loetgering, M. Du, D. B. Flaes, T. Aidukas, F. Wechsler, D. S. P. Molina, M. Rose, A. Pelekanidis, W. Eschen, J. Hess, T. Wilhein, R. Heintzmann, J. Rothhardt, and S. Witte, “PtyLab.m/py/jl: a cross-platform, open-source inverse modeling toolbox for conventional and Fourier ptychography,” arXivarXiv:2301.06595 (2023). [CrossRef]
37. A. Maiden, D. Johnson, and P. Li, “Further improvements to the ptychographical iterative engine,” Optica 4(7), 736–745 (2017). [CrossRef]
38. P. Song, R. Wang, J. Zhu, T. Wang, Z. Bian, Z. Zhang, K. Hoshino, M. Murphy, S. Jiang, C. Guo, and G. Zheng, “Super-resolved multispectral lensless microscopy via angle-tilted, wavelength-multiplexed ptychographic modulation,” Opt. Lett. 45(13), 3486–3489 (2020). [CrossRef]
39. W. Eschen, L. Loetgering, V. Schuster, R. Klas, A. Kirsche, L. Berthold, M. Steinert, T. Pertsch, H. Gross, M. Krause, J. Limpert, and J. Rothhardt, “Material-specific high-resolution table-top extreme ultraviolet microscopy,” Light: Sci. Appl. 11(1), 117 (2022). [CrossRef]
40. B. Enders, “Development and application of decoherence models in ptychographic diffraction imaging,” Ph.D. thesis, Technische Universität München (2016).