Vol. 17 -ray imaging enables the nondestructive investigation of interior structures in otherwise opaque samples. In particular the use of computed tomography Göttingen Series in X(CT) allows for arbitrary virtual slices through the object and 3D information about X-ray Physics intricate structures can be obtained. However, when it comes to image very small structures like single cells, the classical CT approach is limited by the weak ab- sorption of soft-tissue. The use of phase information, encoded in measureable intensity images by free-space propagation of coherent x-rays, allows a huge in- crease in contrast, which enables 3D reconstructions at higher resolutions. In this work the application of propagation-based phase-contrast tomography to lung tissue samples is demonstrated in close to in vivo conditions. Reconstructions of Martin Krenkel the lung structure of whole mice at down to 5 μm resolution are obtained at a self- built CT setup, which is based on a liquid-metal jet x-ray source. To reach even Cone-beam x-ray phase-contrast tomography higher resolutions, synchrotron radiation in combination with suitable holographic phase-retrieval algorithms is employed. Due to optimized cone-beam geometry, for the observation of single cells in whole organs fi eld of view and resolution can be varied over a wide range of parameters, so that information on different length scales can be achieved, covering several millime- ters fi eld of view down to a 3D resolution of 50 nm. Thus, the sub-cellular 3D imaging of single cells embedded in large pieces of tissue is enabled, which paves the way for future biomedical research. Martin Krenkel Cone-beam x-ray phase-contrast tomography for the observation of single cells in whole organs
ISBN 978-3-86395-251-8 ISSN 2191-9860N: Universitätsverlag Göttingen Universitätsverlag Göttingen
Martin Krenkel Cone-beam x-ray phase-contrast tomography for the observation of single cells in whole organs
This work is licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.
Published in 2015 by Universitätsverlag Göttingen as volume 17 in the series „Göttingen series in X-ray physics“
Martin Krenkel
Cone-beam x-ray phase-contrast tomography for the observation of single cells in whole organs
Göttingen series in X-ray physics Volume 17
Universitätsverlag Göttingen 2015 Bibliographic information published by the Deutsche Nationalbibliothek The Deutsche Nationalbibliothek lists this publication in the Deutsche Nationalbibliografie; detailed bibliographic data are available on the Internet at http://dnb.dnb.de.
Address of the Author Martin Krenkel e-mail: [email protected]
Dissertation zur Erlangung des mathematisch-naturwissenschaftlichen Doktorgrades „Doctor rerum naturalium“ der Georg-August-Universität Göttingen
Mitglieder des Betreuungsausschusses: Referent: Prof. Dr. Tim Salditt Institut für Röntgenphysik, Georg-August-Universität Göttingen 1. Korreferent: Prof. Dr. Dr. Detlev Schild Abtl. Neurophysiologie und Zelluläre Biophysik, Universitätsmedizin Göttingen 2. Korreferent: PD Dr. Timo Aspelmeier Institut für Mathematische Stochastik, Georg-August-Universität Göttingen
This work is protected by German Intellectual Property Right Law. It is also available as an Open Access version through the publisher’s homepage and the Göttingen University Catalogue (GUK) at the Göttingen State and University Library (http://www.sub.uni-goettingen.de). The conditions of the license terms of the online version apply.
Layout: Martin Krenkel Cover: Jutta Pabst Cover image: Martin Krenkel
© 2015 Universitätsverlag Göttingen http://univerlag.uni-goettingen.de ISBN: 978-3-86395-251-8 ISSN: 2191-9860 Preface of the series editors
The Göttingen series in x-ray physics is intended as a collection of research monographs in x-ray science, carried out at the Institute for X-ray Physics at the Georg-August-Universität in Göttingen, and in the framework of its related research networks and collaborations.
It covers topics ranging from x-ray microscopy, nano-focusing, wave propagation, image reconstruction, tomography, short x-ray pulses to applications of nanoscale xray imaging and biomolecular structure analysis.
In most but not all cases, the contributions are based on Ph.D. dissertations. The individual monographs should be enhanced by putting them in the context of related work, often based on a common long term research strategy, and funded by the same research networks. We hope that the series will also help to enhance the visibility of the research carried out here and help others in the field to advance similar projects.
Prof. Dr. Sarah Köster Prof. Dr. Tim Salditt Editors Göttingen June 2014
For this book
Depicting single labelled cells in the functional surrounding of an organ or tissue has remained a persisting challenge for modern microscopies. Using phase-contrast x-ray tomography, this work now brings cells to light which are deeply buried in tissue, without any cutting or slicing. With x-rays so famous for their penetration power, tremendous challenges had to be mastered: from issues of coherence to phase front stability, from the realisation of focusing optics to optical zoom, and from robust solutions of the phase problem to efficient numerical implementation. The efforts are rewarded: in full quantitative account of the respective three dimensional structure, one can now study how macrophages migrate and partition in lung tissue. The single cell as a part of the whole.
Prof. Dr. Tim Salditt
Contents
Introduction 1
1 Concepts of x-ray propagation imaging 5 1.1 Propagation of electromagnetic waves ...... 6 1.1.1 Wave equations in vacuum ...... 6 1.1.2 Paraxial wave equation ...... 8 1.1.3 Angular spectrum approach ...... 9 1.1.4 Huygens-Fresnel principle ...... 11 1.1.5 Numerical implementation of propagation ...... 13 1.2 X-ray interactions with matter ...... 15 1.2.1 Wave equations in the presence of matter ...... 15 1.2.2 The refractive index ...... 16 1.2.3 Projection approximation ...... 17 1.3 From wave fields to images ...... 20 1.3.1 Contrast transfer function ...... 21 1.3.2 Transport of intensity equation ...... 23 1.3.3 Imaging regimes ...... 25 1.3.4 Fresnel scaling theorem ...... 27 1.4 Coherence ...... 30 1.5 Dose and resolution ...... 31
2 Phase-retrieval approaches 35 2.1 Direct-contrast regime ...... 36 2.1.1 Pure phase objects ...... 36 2.1.2 Single material objects ...... 39 2.1.3 Bronnikov-aided cyorrection ...... 40 2.2 Holographic regime ...... 43 2.2.1 Holographic reconstruction ...... 43 2.2.2 CTF-based phase-retrieval ...... 45 2.2.3 Holo-TIE ...... 50 2.2.3.1 The ideal case ...... 50 viii Contents
2.2.3.2 Treatment of noisy data ...... 53 2.2.4 Iterative methods ...... 56 2.2.4.1 Projection methods ...... 56 2.2.4.2 Iterative reprojection phase-retrieval ...... 60 2.2.4.3 Iterative Newton methods ...... 62
3 Computed tomography 67 3.1 Mathematical background ...... 67 3.1.1 Radon transform ...... 67 3.1.2 Fourier slice theorem ...... 70 3.1.3 Filtered back projection ...... 71 3.1.4 Algebraic reconstruction technique ...... 73 3.1.5 Cone-beam reconstruction ...... 75 3.2 Artifacts in x-ray tomography ...... 76 3.2.1 Ring artifacts and ring-removal ...... 76 3.2.2 Beam hardening ...... 78 3.2.3 Interior tomography ...... 79 3.2.4 Sampling artifacts ...... 81 3.2.5 Geometry artifacts ...... 82 3.2.6 Motion artifacts ...... 84 3.3 Alignment procedures and algorithms ...... 84
4 Experimental realization 89 4.1 X-ray generation ...... 89 4.2 X-ray detectors ...... 91 4.2.1 Direct detectors ...... 91 4.2.2 Indirect detectors ...... 92 4.3 Laboratory setup JuLiA ...... 96 4.4 Using synchrotron-radiation ...... 99 4.4.1 KB based setups at ID22 & ID16a ...... 99 4.4.2 The empty-beam problem ...... 102 4.4.3 X-ray waveguides ...... 104 4.4.4 Waveguide based setup GINIX ...... 107
5 Tomographic imaging of the mouse lung 111 5.1 State of the art in small-animal imaging ...... 112 5.2Methods...... 113 5.2.1 Sample preparation ...... 113 Contents ix
5.2.2 Experimental setup and measurements ...... 114 5.3 Results ...... 117 5.3.1 Phase-contrast in mouse lungs with laboratory sources . . . 117 5.3.2 The influence of the energy spectrum ...... 123 5.4 Summary and outlook ...... 125
6 Single-cell imaging 127 6.1 High resolution imaging of macrophages ...... 128 6.2Methods...... 129 6.2.1 Sample preparation ...... 129 6.2.2 Staining procedures ...... 132 6.2.3 Experimental parameters ...... 134 6.2.4 Waveguide-stripe removal ...... 136 6.3 2D Imaging of single macrophages ...... 137 6.3.1 Dried macrophage cells ...... 137 6.3.1.1 Single-distance phase-retrieval ...... 138 6.3.1.2 Multi-distance phase-retrieval ...... 141 6.3.2 Stained and unstained cells embedded in resin ...... 144 6.3.3 Towards in-vivo x-ray imaging of macrophages ...... 146 6.4 3D imaging of single macrophages ...... 151 6.4.1 Algorithmic motion correction ...... 151 6.4.1.1 Vertical correction ...... 151 6.4.1.2 Horizontal correction ...... 153 6.4.2 Stained macrophages in aqueous gels ...... 154 6.4.3 Cells embedded in resin ...... 157 6.4.3.1 Osmium stained ...... 157 6.4.3.2 Unstained cells ...... 162 6.4.3.3 Osmium and barium stained ...... 163 6.5 Resolution of the 3D reconstructions ...... 166 6.6 Summary ...... 168
7 Tracking macrophages in the mouse lung 169 7.1 The medically relevant question ...... 169 7.2 Sample preparation ...... 171 7.3 Measurements and phase retrieval ...... 173 7.4 Experimental results ...... 178 7.4.1 Phase-contrast imaging of asthmatic mouse lungs ...... 178 7.4.2 Zoom tomography enabled by x-ray waveguides ...... 183 x Contents
7.4.3 Dose and resolution ...... 188 7.5 Summary ...... 189
8 Conclusions 191
Appendix 195 A.1 The Fourier transform and its properties ...... 195 A.2 Frechét derivative of the propagation operator ...... 196 A.3 List of Matlab functions ...... 197 A.3.1 Image preprocessing ...... 198 A.3.2 Phase retrieval in the direct-contrast regime ...... 199 A.3.3 Holographic phase-retrieval ...... 199 A.3.4 Functions useful for tomography ...... 200
Bibliography 201
Own publications 219
Danksagung 221
Lebenslauf 225 Introduction
For medical diagnostics, x-ray imaging is of great importance as the interactions with matter are so small that internal structures in oblique objects like the human body can be observed in a non-destructive way. Even shortly after the discovery of x-rays in 1895, the technique was used to image bones inside the hand [150], paving the way for a new kind of medical examination. As x-ray images are projections, overlapping structures can not unambiguously be discerned in a single x-ray image. One method to overcome this limitation is computed tomography (CT), in which several recordings under different angles are used to obtain a three-dimensional (3D) volume by means of numerical reconstruction using a computer [77]. In this way, the precise 3D location of internal structures is retrieved, which enables, e.g., the recognition of tumors and its treatment by radiation therapy [80]. In medical research, the precise location of single cells in relation to anatomical structures can help to understand fundamental biological processes. Data with sub-cellular resolution for, e.g., a whole animal would imply handling unreasonable amounts of data. Therefore, the ideal imaging technique should allow the observation of large samples for orientation with the possibility to zoom to specific regions of interest. In classical x-ray imaging the image formation is based on the partial absorption of the radiation, i.e. dense structures like the bones absorb more radiation than lighter elements comprising soft tissue. However, if structures with a weak absorp- tion are observed, the transparency for x-rays may be limiting, as the absorption contrast may be insufficient to obtain a clear image. As x-rays are electromagnetic waves, not only their absorption but also the relative phase difference can carry information about an object. Over the past two decades, x-ray phase-contrast imaging has been developed, for which the underlying physical constants are up to 1000 times larger than for absorption imaging [124, 133]. In particular if very small structures on the micro- and nanometer-scale are of interest, low absorption coefficients become even more restrictive, as reasonable absorption levels build up only over longer path lengths [89]. Thus, phase contrast is an essential tool, to obtain sufficient contrast for high-resolution hard x-ray imaging of biomedical specimen. There are several ways of implementing x-ray phase-contrast, like Zernike phase- 2 Introduction contrast in zone-plate microscopy [34, 156], grating interferometry [39, 111, 138], scanning diffraction microscopy [41, 167] and phase contrast based on coherent free-space propagation behind the sample [37, 125, 185]. Each method has its advantages and drawbacks and can be applied depending on the length scales and applications [99, 193]. For example, if Zernike phase-contrast in zone-plate based microscopes is utilized, a resolution down to some ten nanometers has been achieved in 2D [99], which however can not be achieved on specimen exceeding a thickness of about 300 micron [34]. Grating based phase-contrast imaging offers a very high sensitivity [193] and is compatible with low-brilliance laboratory sources [138], but the best resolution achieved so far is in the range of some micrometers [139], limited to the quality of the necessary gratings. Additionally, the need for optical elements behind the sample increases the radiation dose and imperfections of these elements may limit the quality of the reconstructions. Coherent lensless imaging techniques are potentially most dose efficient, as no optical element behind the sample is needed for the image formation. In scan- ning transmission x-ray microscopy (STXM) with ptychographic phase-retrieval and in coherent diffractive imaging (CDI) resolutions down to 10 nm have been achieved for strongly diffracting test-structures [122, 158, 186] and to about 50 nm on biological specimen [184]. Successful high-resolution reconstructions have been demonstrated in 3D on the organelle [122] and single cell level [82, 119]. However, the typically very small beam size and the scanning overhead are a major drawback of these methods. The time required to perform a typical ptychographic tomog- raphy measurement is in the order of 10 hours [61], which limits the usefulness of this method for specimen with a size of several 100 μm. Propagation-based phase-contrast enables a lensless full-field imaging approach compatible with a wide range of experimental parameters. In particular if a mag- nifying cone-beam illumination is employed, the field of view (FOV) and resolution can be tuned by changing the geometric magnification. The approach is concep- tually simple: The standard radiographic exposure is extended by a free-space propagation and enhanced coherence properties of the illumination . Based on self-interference of the radiation, an in-line hologram can be recorded, which de- pending on the actual propagation distance still shows resemblance to the original object. Especially for small propagation distances, the approach is also compatible with partially coherent laboratory sources [15, 132]. The main difficulty in prop- agation imaging is the phase-retrieval step, which serves to provide quantitative information about the specimen, needed e.g. for tomographic applications. Typ- ically the phase retrieval is based on additional information about the object like 3 measurements in several planes [2, 37] or a priori assumptions to the object (e.g. weak absorption, compact support or positivity of physical constants [17, 54, 151]) which also can be combined [88]. To enable 3D imaging of entire organs at large fields of view and of regions of in- terest at nanoscale resolution with the same technique, propagation-based phase- contrast imaging is the method of choice. Using this approach, reconstructions of whole organs with a low-resolution in the range of some micrometer [15, 126] as well as reconstructions of isolated cells with a high-resolution below 50 nm have been demonstrated [16, 17]. However, to provide a tool for biomedical research, the question arises if high-resolution reconstructions can also be obtained for very large samples. In particular, the observation of cells in relation to anatomical structures in their native environment can help to answer many biomedically rele- vant questions. The goal of the present work is to investigate if propagation-based phase-contrast can be used to obtain 3D structural information with sub-cellular resolution in whole organs. As the lung with its intricate structure is an ideal example how the 3D structure enables a biological function, different samples ranging from whole mice over thick lung-tissue slices down to single alveolar cells (macrophages) are used as speci- mens for proof-of-concept demonstrations in this work. Although macrophages are known to be involved in processes of allergic asthma [114], their exact role is still not well understood [10, 11, 43, 190]. Not only it is interesting whether single cells inside large surrounding tissues can be observed, but an important question is if macrophages are able to migrate through epithelial cells inside the lung. To provide a deeper insight and to answer the raised questions, the present thesis is organized as follows: After a brief introduction of the necessary concepts of Fourier optics, free-space propagation and image formation in chapter 1, several phase-retrieval approaches are presented, which provide a tool to obtain quanti- tative projection images from measured data. In chapter 3 the fundamentals of computed tomography are presented and typical artifacts that may occur in tomo- graphic imaging and their treatments are embraced. A large part of this work was the experimental realization and improvement of existing setups for propagation imaging, so that in chapter 4 a brief summary of the latest synchrotron-based and laboratory setups is given. Experimental applications of propagation based phase- contrast tomography to biomedical specimens are presented in chapter 5, where micrometer-scale 3D reconstructions of whole mouse lungs are obtained with a laboratory source. To achieve reconstructions with nanoscale resolution, following experiments are performed on synchrotron storage rings. 3D reconstructions of 4 Introduction isolated cells are presented at a resolution down to 50 nm and the influence of different contrast agents is investigated in chapter 6. Finally, the question about migration properties of macrophages in mouse lungs is addressed in chapter 7, us- ing high-resolution phase-contrast measurements on thick mouse lung-tissue slices. The thesis is closed with a short conclusion in chapter 8, where further experiments are proposed. 1 Concepts of x-ray propagation imaging
A generic propagation-imaging experiment can be described by the sketch in Fig. 1.1. First, a well-defined, divergent x-ray beam is produced that illuminates the sample. The sample modifies the beam and after a finite propagation distance z2, the phase-shift introduced by the object can be measured as it will transform into intensity variations. This chapter will deal with the theory needed to understand the concept of propa- gation-based phase-contrast tomography. Using the formalism of classical electro- dynamics and Fourier optics, wave propagation will be explained as an essential tool to understand the formation of a measurable intensity image, which plays a crucial role in the derivation and application of phase-retrieval algorithms. A large part of this work was the implementation and evaluation of different ap- proaches, so that chapter 2 will give a detailed overview of available phase-retrieval approaches. Following the introduction of the different propagation concepts, pos-
Figure 1.1: Sketch of the generic propagation-based phase-contrast setup: A well-defined, divergent x-ray beam, emerging from a source, propagates over a distance z1 to the object plane. The object will change the phase distribution of the illuminating beam, leading to a disturbed wave-front. At a propagation distance z2 the phase changes are transformed into intensity changes, which can be measured by a detector. 6 Concepts of x-ray propagation imaging sible interactions between x-ray radiation and biological samples are discussed. The projection approximation will be introduced that simplifies the interac- tion of object and illumination to a single plane, the so-called exit plane behind the object. Equipped with these tools, the formation of a propagation-based phase- contrast image can be described, for which the contrast transfer function plays an important role. The different imaging regimes used for experiments in this work are presented using numerical simulations, showing advantages and challenges for the application to relevant samples.
1.1 Propagation of electromagnetic waves
1.1.1 Wave equations in vacuum
At the beginning of the 20th century it was under discussion if x-rays behave like the previously discovered electromagnetic waves or if they can be fully described by rays without wave properties. Laue could experimentally show that x-rays are indeed electromagnetic waves [50] and therefore have to obey Maxwell’s equations. In vacuum1 these can be written as follows [25, 133]:
∇·E(r,t)=0, (1.1) ∇·B(r,t)=0, (1.2)
∇×E(r,t)+∂tB(r,t)=0 , (1.3)
∇×B(r,t) − 0μ0∂tE(r,t)=0 . (1.4)
Here E is the electric field and B the magnetic induction dependent on the three- T dimensional (3D) coordinate vector r =(x, y, z) ∈ R3 and the time t ∈ R. The T partial derivative in terms of t is denoted by ∂t, and ∇ =(∂x,∂y,∂z) denotes the gradient operator, so that ∇· is the divergence and ∇× the curl operator. −12 −1 −1 −7 −1 −1 0 =8.854 · 10 AsV m is the electric and μ0 =4π · 10 VsA m the magnetic field constant. The Maxwell equations describe the free-space evolution of a general electromagnetic field in space and time. We allow the solutions to these equations to be complex, being aware of the fact that only the real part of a quantity has a physical meaning. By using the Grassman vector identity two
1 Note that SI units are used here and throughout the thesis and vector quantities are indicated using bold letters. 1.1 Propagation of electromagnetic waves 7 independent wave equations can be derived using Eq. (1.1) to (1.4) 2 2 μ ∂t −∇ E(r,t)=0, (1.5) 0 0 2 2 0μ0∂t −∇ B(r,t)=0. (1.6)
Equations (1.5) and (1.6) state that each component of the electric field and mag- netic induction has to obey a wave equation of similar form, which motivates the introduction of a single scalar wave equation, also known as the d’Alembert equation 1 2 2 ∂t −∇ Ψ(r,t)=0. (1.7) c2 Here Ψ may be any component of the electric field E or the magnetic induction B, √ and the definition of the speed of an electromagnetic wave in vacuum c =1/ 0μ0 is introduced. One special solution to this differential equation is given by a monochromatic plane wave
Ψ(PW)(r,t)=e(ik·r−ωt) , (1.8)
T if the wave vector k =(kx,ky,kz) ∈ R3, with modulus |k| =: k =2π/λ and λ being the wavelength, fulfills k = ω/c. In these waves, points with the same phase travel with velocity c, which was experimentally determined to be the speed of light and therefore led to the discovery that light is an electromagnetic wave [133]. In order to find a time-independent equation let us consider an arbitrary wave field Ψ(r,t) that can be described by its spectral Fourier expansion2
ˆ∞ 1 −iωt Ψ(r,t)=√ ψω(r)e dω. (1.9) 2π 0
This decomposition can be understood such that the wave field Ψ(r,t) consists of monochromatic waves that are weighted with ψω(r). We insert Eq. (1.9) into the d’Alembert wave-equation (1.7), change the order of integration and differentiation and perform the time derivative, yielding
ˆ∞ 2 2 ω −iωt ∇ + ψω(r) e dω =0. (1.10) c2 0
To fulfill Eq. (1.10) for general ψω(r), the term in square brackets has to vanish.
2 Note that negative frequencies are ignored here, resulting in “a certain useful analyticity” [133]. 8 Concepts of x-ray propagation imaging
This condition for the time-independent part ψω(r) is known as the Helmholtz equation [133] 2 2 ∇ + k ψω(r)=0 , (1.11) with k = ω/c as in the case of plane waves. To obtain solutions that are not monochromatic, ω-dependent solutions ψω(r) can be put into the ansatz of spectral decomposition (1.9). For the sake of convenience we will ignore the subscript ω in the following equations. Analogously to the time-dependent plane waves Ψ(PW)(r,t), which are solutions to the d’Alembert equation (1.7), one can find monochromatic stationary plane waves ψ(PW)(r) to be solutions to the Helmholtz equation
ψ(PW)(r)=eik·r =ei(kxx+ky y+kz z) , (1.12) where the wave vector k determines the propagation direction. Other fundamental solutions are given by time-independent spherical waves
1 ψ(SW)(r)= eik|r−r0| , (1.13) |r − r0| for r = r0. These waves do not have a distinct propagation direction. Instead, concentric spheres of constant phase around an origin r0 will propagate equally in all directions.
1.1.2 Paraxial wave equation
A way to further simplify the wave equation is to use the so-called small-angle, Fresnel or paraxial approximation. A wave is called paraxial, if the normal vectors of the wavefront form small angles with respect to a given optical axis. This is a good description for “beam-like” waves that propagate in one direction and thus vary mainly perpendicular to that direction [133]. To derive the paraxial wave equation we start with the separation ansatz
ψ(r)=ψ(r) · eikz , (1.14) where, without loss of generality, we define the z-axis to be the optical axis. This wave consists of a fast oscillating plane-wave term eikz along the optical axis and an envelope ψ(r) that changes slowly along the optical axis. Inserting the ansatz (1.14) into the Helmholtz equation (1.11), performing the derivative and dividing 1.1 Propagation of electromagnetic waves 9 all terms by eikz yields 2 2 ∇⊥ + ∂z +2ik∂z ψ (r)=0. (1.15)
2 2 2 2 Here the Laplacian ∇ is split up into the lateral part ∇⊥ := ∂x +∂y and the axial 2 part ∂z . The paraxial approximation manifests itself in neglecting the second 2 derivative ∂z ψ (r). The resulting equation is called the homogenous paraxial equation 2 ∇⊥ +2ik∂z ψ (r)=0. (1.16)
One simple solution to this equation is given by a so-called parabolic beam
2 1 ikr /2z ψ (PB)(r)= e ⊥ , (1.17) z
T where r⊥ =(x, y) is introduced. This solution can be obtained by assuming a spherical wave originating at r = 0 and performing a small-angle approximation for the modulus of r [14]. Another important solution is the Gaussian beam, described in detail e.g. in [182]. It can be used to approximate beams exiting x-ray waveguides [90], which play an important role in this thesis.
1.1.3 Angular spectrum approach
Let us consider a solution to the Helmholtz equation ψ(r⊥,z) known in a plane perpendicular to the optical axis at position z. To understand how this field will change during propagation along the optical axis, we decompose it in lateral Fourier-components ˆ 1 ik⊥r⊥ ψ(r⊥,z)= ψ˜(k⊥,z)e dk⊥ , (1.18) 2π R2 with k⊥ =(kx,ky) and ψ˜(k⊥,z) denoting the lateral Fourier transform of ψ(r⊥,z). In the following, Fourier-transformed quantities will be marked with a tilde (~). Inserting Eq. (1.18) into the Helmholtz equation (1.11) yields a second-order, scalar differential equation in terms of z for ψ˜(k⊥,z). This initial value problem can be solved to yield the general solution [55] √ z k2−k2 iΔ ⊥ ψ˜(k⊥,z1)=ψ˜(k⊥,z0)e , (1.19) 10 Concepts of x-ray propagation imaging
2 2 2 2 with k⊥ = |k⊥| = kx + ky,Δz = z1 − z0 and z1 >z0 . Using this solution Eq. (1.18) reads ˆ √ 1 z k2−k2 ik r iΔ ⊥ ⊥ ψ(r⊥,z1)= ψ˜(k⊥,z0)e e dk⊥ . (1.20) 2π R2
This equation conveys how to propagate an arbitrary wave field, known at a dis- tance z0, to another distance z1 = z0 +Δz. Note that the propagation is only dependent on Δz and that the absolute value of z0 does not matter. Hence, we will set z0 = 0 in the following and regard Eq. (1.20) as a rule to propagate a given field by a relative distance Δz = z − 0=z. With this convention and with the operator F⊥ denoting the lateral Fourier transform, Eq. (1.20) can be written as √ − z k2−k2 1 i ⊥ ψ(r⊥,z)=F⊥ e F⊥ (ψ(r⊥, 0)) . (1.21)
Equation (1.21), also known as the free-space propagator [133], states that for the propagation of a wave field over a distance z, the initial wave field is laterally Fourier transformed, multiplied with a complex phase factor and back transformed to real space. In particular, the lateral dimensions defined by r⊥ do not change, so this propagator is well suited for near-field propagation.3 2 2 Note that for k⊥ >k the square-rooted term in Eq. (1.21) is negative. Hence, the exponential term will be real and the wave will be damped out within a range of a few wavelengths instead of being an oscillating, propagating wave. In this case one speaks of an evanescent wave. If there are no evanescent waves, the operator defined by Eq. (1.21) is unitary [112] and the equation can be used for z<0 to perform a back propagation, also known as inverse diffraction [160]. For experiments in this work, the x-ray beams are well described by paraxial beams. Therefore, we now derive a paraxial operator in analogy to Eq. (1.21). For the wave vector of such a beam it is [53]
2 2 k⊥ k . (1.22)
This motivates the Taylor expansion of the phase factor in Eq. (1.21)