X-ray phase imaging : Demonstration of extended conditions with homogeneous objects

We discuss contrast formation in a propagating x-ray beam. We consider the validity conditions for linear relations based on the transportof-intensity equation (TIE) and on contrast transfer functions (CTFs). From a single diffracted image, we recover the thickness of a homogeneous object which has substantial absorption and a phase-shift of −0.37radian. © 2004 Optical Society of America OCIS codes: (340.7480) X-rays; (340.7440) X-ray imaging; (340.7460) X-ray microscopy References and links 1. G. Schmahl, D. Rudolph and P. Guttmann, “Phase contrast x-ray microscopy – experiments at the BESSY storage ring,” in X-ray Microscopy II, D. Sayre, M. Howells, J. Kirz and H. Rarback, eds., Vol 56 in Springer Series in Optical Sciences (Springer-Verlag, Berlin, 1988), pp. 228–232. 2. D. Sayre and H. N. Chapman, “X-ray microscopy,” Acta. Cryst. A51, 237–252 (1995). 3. K. A. Nugent, T. E. Gureyev, D. F. Cookson, D. Paganin and Z. Barnea, “Quantitative phase imaging using hard x rays,” Phys. Rev. Lett. 77, 2961–2964 (1996). 4. A. Pogany, D. Gao and S. W. Wilkins, “Contrast and resolution in imaging with a microfocus x-ray source,” Rev. Sci. Instrum. 68, 2774–2782 (1997). 5. U. Bonse and M. Hart, “An x-ray interferometer,” Appl. Phys. Lett. 6, 155–157 (1965). 6. A. Momose, T. Takeda and Y. Itai, “Phase-contrast x-ray computed tomography for observing biological specimens and organic materials,” Rev. Sci. Instrum. 66, 1434–1436 (1995). 7. Y. Kohmura, H. Takano, Y. Suzuki and T. Ishikawa, “Shearing x-ray interferometer with an x-ray prism and its improvement,” in Proc. 7th Intern. Conf. on X-ray Microscopy, D. Joyeux, F. Polack, eds., J. de Physique IV, Vol 104, J. Susini (EDP Sciences, Les Ulis, 2003), pp. 571–574 . 8. T. Wilhein, B. Kaulich, E. Di Fabrizio, F. Romanato, M. Altissimo, J. Susini, B. Fayard, U. Neuhäusler, S. Cabrini and F. Polack, “Differential interference contrast x-ray microscopy with twin zone plates at ESRF beamline ID21,” in Proc. 7th Intern. Conf. on X-ray microscopy, ibid. 9. E. M. Di Fabrizio, D. Cojoc, S. Cabrini, B. Kaulich, J. Susini, P. Facci and T. Wilhein,“Diffractive optical elements for differential interference contrast x-ray microscopy,” Opt. Express 11, 2278–2288 (2003), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-11-19-2278. 10. E. Forster, K. Goetz and P. Zaumseil, “Double crystal diffractometry for the characterization of targets for laser fusion experiments,” Krist. Tech. 15, 937–945 (1980). 11. J. R. Palmer and G. R. Morrison, “Differential phase contrast imaging in the scanning transmission x-ray microscope,” in Short Wavelength Coherent Radiation, P. H. Bucksbaum and N. M. Ceglio, eds., Vol. 11 of OSA Proceedings Series (Optical Society of America, Washington, D. C., 1991), pp. 141–145. 12. C. Jacobsen, M. Howells, J. Kirz and S. Rothman, “X-ray holographic microscopy using photoresists,” J. Opt. Soc. Am. A 7, 1847–1861 (1990). 13. L. J. Allen, W. McBride and M. P. Oxley, “Exit wave reconstruction using soft x-rays,” Opt. Commun. 233, 77–82 (2004). 14. M. R. Teague, “Deterministic phase retrieval: Greens function solution,” J. Opt. Soc. Am. A 73, 1434–41 (1983). 15. D. Paganin and K. A. Nugent, “Noninterferometric phase imaging with partially coherent light,” Phys. Rev. Lett. 80, 2586–2589 (1998). (C) 2004 OSA 28 June 2004 / Vol. 12, No. 13 / OPTICS EXPRESS 2960 #4387 $15.00 US Received 20 May 2004; revised 19 June 2004; accepted 20 June 2004 16. T. E. Gureyev, A. Pogany, D. M. Paganin and S. W. Wilkins, “Linear algorithms for phase retrieval in the Fresnel region,” Opt. Commun. 231, 53–70 (2004). 17. X. Wu and H. Li, “A general theoretical formalism for x-ray phase contrast imaging,” J. X-ray Sci. Tech. 11, 33–42 (2003). 18. J.-P. Guigay, R. H. Wade and C. Delpha, “Optical diffraction of Lorentz microscope images,” in Proceedings of the 25th meeting of the Electron Microscopy and Analysis Group, W. C. Nixon, ed. (The Institute of Physics, London, 1971), pp. 238–239. 19. J.-P. Guigay, “Fourier transform analysis of Fresnel diffraction patterns,” Optik 49, 121–125 (1977). 20. D. Paganin, S. C. Mayo, T. E. Gureyev, P. R. Miller and S. W. Wilkins, “Simultaneous phase and amplitude extraction from a single defocused image of a homogeneous object,” J. Microsc 206, 33–40 (2001). 21. P. Cloetens, W. Ludwig, J. Baruchel, D. Van Dyck, J. Van Landuyt, J.-P. Guigay and M. Schlenker, “Holotomography: Quantitative phase tomography with micrometer resolution using hard synchrotron radiation x-rays,” Appl. Phys. Lett. 75, 2912–2914 (1999). 22. M. H. Maleki and A. J. Devaney, “Noniterative reconstruction of complex-valued objects from two intensity measurements,” Opt. Eng. 33, 3243–3253 (1994). 23. A. N. Tikhonov and V. Y. Arsenin, “Solutions of Ill-posed Problems” (V. H. Winston, Washington D.C., 1977). 24. E. C. Harvey and P. T. Rumsby, “Fabrication techniques and their application to produce novel micromachined structures and devices using excimer laser projection,” in Micromachining and Microfabrication Process Technology III, S. Chang and S. W. Pang, eds., Proc. SPIE 3223, 26–33 (1997). 25. L. D. Turner, K. P. Weber, D. Paganin and R. E. Scholten, “Off-resonant defocus-contrast imaging of cold atoms,” Opt. Lett. 29, 232–234 (2004). 26. L. D. Turner, K. F. E. M. Domen, W. Rooijakkers and R. E. Scholten, School of Physics, University of Melbourne 3010, Australia are preparing a manuscript to be called “Holographic imaging of cold atoms”. 27. M. Centurion, Y. Pu, Z. Liu, D. Psaltis and T. W. Hänsch, “Holographic recording of laser-induced plasma,” Opt. Lett. 29, 772–774 (2004).


Introduction
Absorption contrast has been the principal imaging mode for x-rays for over 100 years, nonetheless there has been considerable development recently in implementing phase contrast techniques.Phase contrast can be strong when absorption contrast is minimal; for instance for low-Z materials or for high energy x-rays [1,2,3].Phase contrast can also be used without delivering a high dose to the sample [2].And certain phase methods require no additional optics leading to source-limited, rather than optics-limited, resolution [4].
Demonstrated methods for obtaining the phase from an x-ray wavefield are now legion.Interferometric methods include the use of Bonse and Hart [5] type interferometers [6], shearing interferometers [7], differential interference contrast arrangements using two zone plates [8] or diffractive optical elements [9].Zernike phase contrast [1] has also been demonstrated.Other methods include refraction measurement using crystal diffraction [10], and segmented detectors [11].Propagation-based methods have also been developed and involve recovery of the object phase and/or amplitude from one or more measurements of the object diffraction pattern.Methods include in-line holography [12], iterative schemes [13] and approaches based on solution of the equations governing the propagation of the wavefield [14,15,3].
In this paper we consider aspects of the latter method.Here the Fresnel integral describing the diffracted intensity is linearized allowing a straightforward retrieval of object phase and transmission.A first order Taylor expansion obtains the transport of intensity equation (TIE) solution [16,3].Alternatively a Born-type approximation [17,16] gives a solution identical in form to that derived by Guigay [18,19] in the context of electron microscopy.We will refer to this as the contrast transfer function (CTF) solution.Both the TIE and the CTF solution can be further simplified under the assumption of a homogeneous object [20,16].This permits the thickness distribution of an object to be retrieved from a single diffracted image.
In the TIE it is the first-order Taylor expansion that restricts the validity of the solution.Other than a requirement for paraxiality, there is no limitation on the magnitude of the phase or the absorption.On the other hand, the Born-type approximation previously used in deriving the CTF method can be quite restrictive.Guigay showed [19] that a less restrictive requirement applies for a pure phase object.Like the TIE condition, this condition depends on the feature sizes present as well as the wavelength and propagation distance.
Here we show that the Fresnel-diffracted intensity can be linearized at a later stage of the derivation.As a consequence we find that the less restrictive phase condition extends to weaklyabsorbing objects.In Section 2 we outline the derivations of the TIE and the CTF solutions and their validity conditions and show how the assumption of a homogeneous object allows retrieval of the object thickness from a single diffracted image.In Section 3 we present experimental results demonstrating quantitative thickness retrieval under our derived validity conditions.

TIE and CTF
We begin with the Fourier transform, F[•], of the intensity I(r, z) obtained under Fresnel diffraction of an object-plane wavefield f (r) = f (r, z = 0) with wavelength λ [18]: The transverse spatial coordinates and their corresponding Fourier conjugates are given by r and u respectively.To obtain the TIE solution we Taylor expand the wavefield to first order, where ∇ is the gradient operator in the transverse plane.Substituting into Eq.( 1) gives where φ is the phase of the wavefunction f .The validity condition for the TIE solution is therefore that the higher order Taylor expansion terms can be disregarded: This condition can always be satisfied by choosing a sufficiently small propagation distance, while requiring no approximation regarding the magnitude of the amplitude or phase.
The previous approach to obtain the CTF solution was to write the object wavefunction in terms of its absorption µ and phase φ components such that: where the intensity of the plane wavefield incident on the object.Then the Borntype approximation of µ 1 and |φ | 1 was made so that Substituting Eq. ( 6) into Eq.( 1) and retaining µ and φ to first order, obtained Guigay's result in which δ (u) denotes the Dirac delta distribution.When the object is pure phase (µ = 0), Guigay [19] also showed, by substituting Eq. ( 5) into Eq.( 1), that the corresponding form of Eq. ( 7) can be obtained if, for all r, This is sometimes referred to as the slowly-varying phase condition.It should be noted that the displacement vector λ zu, over which points in the phase should be similar, is a function of propagation distance z and of spatial frequency u.

TIE for a homogeneous object
Consider the case of an optically thin and homogeneous object for which µ(r) = kβ T (r) and φ (r) = −kδ T (r), where k = 2π/λ , the refractive index is n = 1 − δ + iβ and T is the thickness of the object.We substitute into Eq.( 1), re-factor and make a first-order Taylor expansion of T yielding If we assume that |λ zu • ∇φ (r)| 1 then we can expand the second exponential to first order, and then applying the Fourier derivative theorem ] we obtain: This may be solved for the thickness T : The validity condition on the thickness T is similar to that for the TIE in Eq. (4).

Derivation of extended validity CTF solution
Here we begin by substituting Eq. ( 5) into Eq.( 1): Assuming both real and imaginary parts of the exponential are small, we expand, noting the Fourier transforms: Applying the Fourier shift theorem to each term and rearranging recovers Eq. (7).If the object is homogeneous then, substituting Eq. ( 9), we can retrieve the thickness T from a single diffracted image: The linearizing assumption made in obtaining Eq. ( 15) is that 2µ(r) 1 and |φ (r These conditions on the validity of Eq. ( 7) are much less stringent than previously realised.While the object must be weakly absorbing it need not be non-absorbing, which is less restrictive than the pure phase assumption made by Guigay.Furthermore, the weak phase condition required in the Born-type approach is here relaxed to the slowly-varying condition Eq. ( 8).
Assuming object homogeneity allows us to incorporate the effect of object absorption.Consequently the contrast transfer function, which is the denominator term of Eq. ( 16), is non-zero at u = 0. Nulls in the CTF at frequencies above u = 1/ √ 2λ z may be avoided by incorporating additional images at different propagation distances [21,22], or simply by Tikhonov regularisation [23].We find that Tikhonov regularisation leaves only minor artefacts at high frequencies whereas assuming a pure phase object introduces intractable low-frequency instability.Fig. 1.Inverse of the contrast transfer functions for the TIE (blue dotted line) and CTF (red line) forms calculated for an infinite grid with spatial feature size of 1.9 µm and the experimental conditions.The green dashed line is at the experimental distance.The embedded movie shows the inverse CTFs as shown here with a slider indicating the z position corresponding to diffraction distance.Also shown in the movie is a plot of an input amplitude (green) and the amplitude retrieved using either the CTF (red) or TIE (blue) methods for the indicated propagation distance.The CTF method correctly accounts for the contrast reversals that arise on propagation.The TIE method should only be applied for z closer than the first contrast reversal; it may retrieve inverted amplitudes if applied at greater z.

Experimental results
Laser ablation [24] was used to etch a grid of lines in a polyimide film (composition C 22 H 10 N 2 O 4 and density 1.45 g cm −3 ).The lines were measured by atomic force microscope (AFM) to be approximately 1.9 µm apart and 90 nm high and were superimposed on an 80 µm square of 650 nm in height.
The experiments were performed at a wavelength of 0.436 nm and with z = 0.0475 m, where z is as defined above.Through the thickest part of the object the transmission exp(−2µ) is 98.7% and the phase-shift φ is −0.37 radian.Accordingly, neither the Born-type approximation nor the phase-only requirement are met.However, for most positions and spatial frequencies the slowly-varying phase condition Eq. ( 8) is obeyed.
Figure 1 shows a plot of the inverse of the contrast transfer function from the CTF solution Eq. ( 16) for the 1.9 µm grid features.A contrast transfer function can also be defined, in the weak absorption limit, for the TIE solution Eq. ( 13) and its inverse is also shown in Fig. 1.Inspection shows that the TIE solution should yield grid features reversed in contrast compared to the CTF solution.No contrast reversal is expected for the 80 µm square.
The experiments were performed at beamline 2-ID-B at the Advanced Photon Source.A beam of 2844 eV (∆E/E 10 −3 ) x-rays with FWHM size 1.5 × 0.5 mm illuminated a 160 µm diameter gold zone plate with an outer zone width of 50 nm.The focal length was 18.4 mm and a 10 µm diameter order sorting aperture was placed at a distance of 17.2 mm from the zone plate.A 30 µm wide gold beamstop blocked the zero order beam.This zone plate configuration provided a point source of illumination at a distance of R 1 = 54.7 mm from the object.After passing through the object the beam was allowed to propagate through a He-filled flight tube a distance of R 2 = 363 mm onto a crystal scintillator which was imaged through a 20× objective by a CCD camera with 13.5 µm pixels.This expanding beam mode is analyzed using the parallel beam mode derivations described above using the conversion [4] where AFM measurements also confirm the presence of grid lines outside the square.These are not a retrieval artefact, unlike the circular fringes around the contaminant at centre right.
The contaminating material presumably violates the assumption of an homogeneous object.
pattern due to a point source a distance R 1 behind the object and I ∞ is parallel beam diffractionpattern.The object thickness map was retrieved from the measured diffraction pattern using both the CTF solution Eq. ( 16) with Tikhonov regularization, and the TIE solution Eq. ( 13), and is shown in Fig. 2(a).We line average the essentially one-dimensional image in the region outlined in white.Both CTF and TIE methods retrieve the mean height of the 80 µm square as 610 ± 50 nm, concurring with the AFM result of 650 nm.However, the TIE retrieval grossly underestimates the height of the grid pattern and, as predicted in Fig 1, reverses the contrast (Fig. 2(b)).As expected from the validity conditions for Eq. ( 16), the CTF retrieval of the grid pattern (Fig. 2(b)) is in excellent agreement with the AFM measurement of 90 nm.

Conclusions
We have explored the validity conditions for the linear CTF expression relating object phaseshift and absorption to the contrast of the Fresnel diffraction-pattern.The linear expression is found not to be restricted to weakly phase-shifting objects: it applies to a substantially wider class of objects which show weak absorption and slowly-varying phase.If an object is made of one material with known complex refractive index, the CTF expression may be inverted to retrieve the object thickness from a single diffracted image.We demonstrated an example where the CTF solution could correctly retrieve thickness features of a weaklyabsorbing object with large, but slowly-varying, phase-shift.Thickness features at two wellseparated spatial frequencies were retrieved by the CTF solution while the TIE retrieval was valid only at the lower spatial frequency.These results augur well for wider applications of the CTF technique such as imaging cold atom clouds [25,26] and plasmas [27].

R 2 )Fig. 2 .
Fig. 2. (a) CTF-retrieved thickness map for the square with grid lines.(b)Column-average of retrieved thickness for the grid pattern in the region shown in (a) for the TIE solution (blue) and the CTF solution (red).The AFM result (green) shows excellent agreement.AFM measurements also confirm the presence of grid lines outside the square.These are not a retrieval artefact, unlike the circular fringes around the contaminant at centre right.The contaminating material presumably violates the assumption of an homogeneous object.