TY - JOUR AU1 - Yashiro,, Wataru AB - Abstract The self-imaging phenomenon referred to as the Talbot effect in the field of optics was discovered by H.F. Talbot in the 1830s, and is now widely used for imaging using not only visible light but also X-rays, electrons, neutrons, and matter waves. In this review, the author introduces the current progress being made in hard-X-ray imaging microscopy based on the self-imaging phenomenon. Hard-X-ray imaging microscopy is a promising technique for non-destructively visualizing internal structures in specimens with a spatial resolution up to a few tens of nanometers. The use of the self-imaging phenomenon makes it possible to realize highly sensitive phase-contrast X-ray imaging microscopes. These approaches have several advantages over conventional X-ray imaging microscopes, including the widely used Zernike X-ray phase-contrast microscopes, and can provide a powerful way of quantitative visualization with a high spatial resolution and a high sensitivity even for thick specimens. X-ray imaging, X-ray microscopy, grating, Talbot effect, interferometry Introduction Non-destructive and quantitative high-spatial-resolution visualization of internal structures of materials consisting of light elements will bring about significant progress in biological and materials sciences. Hard-X-ray imaging microscopy is a technique for non-destructively visualizing internal structures of specimens with a spatial resolution up to a few tens of nanometers [1]. It has an advantage over other X-ray microscopies such as X-ray scanning and diffraction microscopies in that it does not require time-consuming measurements. However, similar to the other X-ray microscopies, it is less sensitive to materials consisting of lighter elements because hard X-rays weakly interact with them. One way to improve the sensitivity of hard-X-ray imaging microscopy is by using the phase-shift of X-rays passing through a specimen [2–4]. In principle, this so-called X-ray phase-contrast imaging microscopy provides three orders of magnitude higher sensitivity than conventional ones based on absorption contrast. Several X-ray phase-contrast imaging microscopes have been reported since the mid-1990s, but they require a high-brilliance synchrotron X-ray source or specially designed optics. Zernike phase-contrast microscopy for visible light [5,6] is attainable even for hard X-rays [7–10] and it works with a compact low-brilliant laboratory X-ray source, but it is quantitative only for a week-phase object. In the last decade, X-ray phase-contrast imaging techniques using the self-images (referred to as ‘Fourier images’ [11] in the field of electron microscpy) of a grating generated by the Talbot effect [12–16], which is called X-ray grating interferometry [17–52], have attracted much attention. Because X-ray grating interferometry enables highly sensitive X-ray imaging with a continuous-spectrum low-brilliant laboratory X-ray source, it is suitable for medical diagnostics. Besides, it works with spherical-wave X-rays, which cannot be used other X-ray phase-contrast imaging techniques such as crystal interferometry and diffraction enhanced imaging (DEI). This fact enables the use of an X-ray lens, and X-ray phase-contrast imaging microscopes can be constructed on the basis of X-ray grating interferometry [27,29,30,32,34,40,50]. In this review, X-ray phase-contrast imaging microscopy is reviewed. The next section is devoted to the theoretical description of the self-imaging phenomenon. In Section 3, X-ray phase-contrast imaging using the X-ray grating interferometry is reviewed. In Section 4, several phase-contrast X-ray imaging microscopes using the self-imaging phenomenon are introduced. Finally, in Section 5, the future prospects of the X-ray imaging microscopes are discussed and conclusions are made. Self-imaging phenomenon The self-imaging phenomenon called the Talbot effect was discovered by H.F. Talbot [12]. A typical setup where the phenomenon can be observed is shown in Fig. 1, where spatially coherent quasi-monochromatic light emanated from a point source illuminates a one-dimensional periodic pattern (diffraction grating) located on the (x1,y1)-plane at a distance R1 from the source. When the pitch d1 of the grating is much larger than the wavelength λ of the light, the images of the grating, called self-images or Talbot images, are formed at specific distances downstream of the grating. Fig. 1. View largeDownload slide (Color online) Illustration of setup where Talbot effect occurs. Fig. 1. View largeDownload slide (Color online) Illustration of setup where Talbot effect occurs. This self-imaging phenomenon can be attributed to the Fresnel diffraction by the grating as follows. For simplicity, we consider the electric field of light close to the optical axis (z-axis). In addition, we assume that the electric field is polarized in the y-direction, but the following theoretical description is approximately correct even that polarized in the x-direction because only the Fresnel diffraction close to the optical axis (in the paraxial approximation) contributes to the self-imaging phenomenon. The electric field E2 on the (x2, y2)-plane at a distance z12 downstream from the grating can be written in the paraxial approximation by E2(x2,y2)≈−iE10λexp[2πiz12λ]z12×∬T(x1+x′,y1+y′)×exp[iπ(x′2+y′2)λzeff,12]dx′dy′. (1) Here, E10≡E0exp[2πiR1λ]exp [πi(x12+y12)λR1]R1 corresponds to the electric field just in front of the grating, T(x1,y1) is the complex transmission function of the grating, zeff,12 is the effective distance defined by zeff,12≡(1R1+1z12)−1 ⁠, and (x2,y2)≡R2R1(x1,y1) ⁠, where R2≡R1+z12 ⁠. If we assume that the size of the grating is regarded to be infinite, T(x1,y1) can be expanded by a Fourier series ∑nanexp[2πinx1d1] ⁠, where an is the n th Fourier component of the transmission function. Using an ⁠, Eq. (1) can be rewritten by E2(x2,y2)≈E20∑nbnexp[2πinx2d2], (2) where E20≡E0exp[2πiR2λ]exp [πi(x22+y22)λR2]R2 corresponding to the electric field on the (x2,y2)-plane when no grating would be placed, d2≡R2R1d1 ⁠, and bn is defined by bn≡anexp[−iπn2λzeff,12d12]. (3) Equation (3) means that if zeff,12 is an even-integer multiple of d12λ ⁠, the electric field just behind the grating is completely reproduced with a magnification of R2R1 ⁠. This reproduction of the electric field forms the self-image. Historically, the effective distance where the first self-image appears, i.e. 2d12λ ⁠, which was first found for nearly plane-wave illumination by Rayleigh [53], is called the Talbot distance. Note that even at an odd-integer multiple of d12λ ⁠, the electric field just behind the grating is also reproduced but with a half-pitch shift. It is convenient to define the Talbot order p by p≡λzeff,12d12. (4) By using p ⁠, Eq. (3) can be rewritten as bn(p)≡anexp[−iπn2p], (5) In the cases where an absorption grating is used, the intensity of light just behind the grating is reproduced with a magnification of R2R1 when P is a positive integer (the intensity is shifted by a half pitch when P is an odd positive integer). In the cases where a phase grating is used, an intensity distribution that is the same as the shape of the grating, which is also called a self-image, often appears at a fractional Talbot order. This phenomenon is referred to as the fractional Talbot effect. For example, when a rectangular π/2-phase grating with a duty cycle of 0.5 is used, its self-images appear at odd-integer multiples of 12 [54], whereas when a rectangular π-phase grating with a duty cycle of 0.5 is used, its self-images with a compression ratio α of 2 (meaning that the pitch of a self-image is given by d22 ⁠) appear at odd-integer multiples of 18 [55]. At a Talbot order other than such specific Talbot orders, various intensity distribution, which is generally referred to as the Fresnel image [56], appears. Note that in the case of two dimensional grating, the Talbot distance is given by the minimum distance where λzeff,12|Na1⁎+Ma2⁎| is equal to an even integer for all N and M [56]. Here, N and M are integers and a1⁎ and a2⁎ are reciprocal lattice vectors defined by the primitive lattice vectors a1 and a2 of the two dimensional grating. We can interpret Eq. (2) in geometrical wave optics [52]. In fact, we can rewrite Eq. (2) into E2(x2,y2)≈E20∑nTn(x1−Δxn)exp[2πiΔlnλ], (6) where Tn(x1)≡anexp[2πinx1d1], (7) Δxn≡nλd2z12, (8) =npd1, (9) Δln≡Δxn22zeff,12. (10) Equation (6) can be interpreted as follows: the nth order diffracted wave propagates along the path passing through the point Pn shown in Fig. 2, and Δln corresponds to the optical path difference between the nth and 0th orders, passing through the points Pn and P0, respectively. Fig. 2. View largeDownload slide (Color online) Schematic illustration of geometrical-optics interpretation of Talbot effect for monochromatic X-rays. Electric field at distance z12 downstream of grating (G1) is formed from interference of waves passing through points Pn on G1 (n=0,±1,±2,⋯) ⁠, where Pn is separated at distance of Δxn from P0. Here, Δln is optical path difference between nth and 0th order. Fig. 2. View largeDownload slide (Color online) Schematic illustration of geometrical-optics interpretation of Talbot effect for monochromatic X-rays. Electric field at distance z12 downstream of grating (G1) is formed from interference of waves passing through points Pn on G1 (n=0,±1,±2,⋯) ⁠, where Pn is separated at distance of Δxn from P0. Here, Δln is optical path difference between nth and 0th order. The intensity I2(x2,y2) of the image at the distance of z12 from the grating is given by I2(x2,y2)≈I20∑mTm+n(x1−Δxm+n)Tm⁎(x1−Δxm)×exp[2πiΔlm+n−Δlmλ], (11) where I20≡|E20|2 ⁠. Equation (11) is convenient when we consider the effect of spatial and temporal coherences of light: the spatial- and temporal-coherence lengths on the grating that are required for mth and (m+n)th orders to interfere are given by Δxn and Δln ⁠, respectively. For a finite size of quasi-monochromatic light source, the effect of the spatial coherence can be described by the van Cittert–Zernike theorem [6,15] and can be represented by a complex coherence factor μn [6]; using μn ⁠, I2(x2,y2) can be expressed in the following form [47]: I2(x2,y2)≈I20∑mμnTm+n(x1−Δxm+n)Tm⁎(x1−Δxm)×exp[2πiΔlm+n−Δlmλ]. (12) For the self-image to be observed, it is necessary for the mth and (m+1)th orders to interfere. Because the mth and (m+1)th order waves propagate along two different paths separated by pd1 on the grating (see Fig. 2), a spatial coherence length that is comparable or larger than pd1 is necessary on the grating. On the other hand, the self-image can be observed even for a broad spectrum because, from the Wiener–Khintchine’s theorem [6], the temporal coherence length of the light used is roughly given by λ2Δλ ⁠, where Δλ is the bandwidth of the light, and the temporal coherence that is necessary for the mth and (m+1)th orders to interfere is Δl1 ⁠, which is comparable to λ when d1≫λ ⁠. X-ray grating interferometry In the X-ray grating interferometry based on the Talbot effect, the phase shift caused by a sample is detected by the deformation of the self-image. When a sample is located at a distance of Rs from the X-ray source in Fig. 1, the electric field on the (x2,y2)-plane is expressed from Eqs. (2) and (6) by [47] E2(x2,y2)≈E20∑nbn(p)T's(xs−npsd1,ys)×exp[2πinx2d2]. (13) Here, (xs,ys)≡RsR1(x1,y1) ⁠, ps is the effective Talbot order given by ps≡{pRsR1(0≤Rs≤R1)pR2−RsR2−R1(R1≤Rs≤R2), (14) and E20T's(xs,ys) corresponds to the electric field on the (x2,y2)-plane that formed because of the free-space propagation of X-rays from the sample to the plane when no grating is located: T's(xs,ys) is expressed using the complex transmission function Ts(xs,ys) of the sample by T's(xs,ys)≈∬Ts(xs+x′s,ys+y′s)exp[iπ(x′s2+y′s2)λzeff,s2]dx′dy′, (15) where zeff,s2≡(1Rs+1R2−Rs)−1 ⁠. The complex transmission function Ts(xs,ys) of the sample represents the effects of the absorption and the phase-shift of the sample: Ts(xs,ys)=exp[−μs(xs,ys)2]exp[−iΦ(xs,ys)], (16) where μs(xs,ys) corresponds to the linear absorption coefficient and Φ(xs,ys) is the phase shift caused by the sample. The intensity of the self-image is expressed from Eqs (12) and (13) by I2(x2,y2)≈I20∑mμnbm+n(p)bm⁎(p)T's(xs−(m+n)psd1,y1)×T's⁎(xs−mpsd1,y1)exp[2πinx2d2]. (17) Equation (17) means that this optical configuration works as a multi-wave interferometer: the deformation of the self-image due to the insertion of a sample includes the effect of interference between the (m+1)th and mth orders. In other words, we can realize a shearing interferometer with a shearing distance of psd1 by retrieving the 1st order Fourier component of the intensity of the self-image. From geometrical consideration (see Fig. 2), it can be shown that the spatial coherence length on the grating that is necessary for the shearing interferometer is given by pd1. This is the reason why the interferometer does not require high-spatial-coherence X-rays. Note that, if the shearing distance is larger than the size of the sample, the interferometer works as a Mach–Zehnder interferometer (an array of Mach–Zehnder interferometers), which provides a phase image (see Fig. 6 (a)) as X-ray and electron holographies do [57,58]. A fringe scanning [59–61] or Fourier transform method [62] is commonly used to retrieve the first order Fourier component q1(x2,y2) of the intensity of the self-image. If we neglect the effect of the free-space propagation, i.e. T's(xs,ys)≈Ts(xs,ys) ⁠, the phase difference between Φ(xs−(m+1)psd1,ys) and Φ(xs−mpsd1,ys) is obtained from the argument of q1 ⁠. Note that when a rectangular grating with a duty cycle of 0.5 is used and the compression ratio of the self-image is 1, only the three terms a−1 ⁠, a0 ⁠, a1 contribute to q1(x2,y2) ⁠, so that the optical configuration works as a three-wave interferometer. In fact, from Eq. (17), the first order Fourier component of the intensity of the self-image is proportional to P(x2,y2)=μ1(b0b−1⁎Ts(xs,ys)Ts⁎(xs+psd1,ys)+b1b0⁎Ts(xs−psd1,ys)Ts⁎(xs,ys))×exp[2πix2d2], (18) and for a weakly absorbing sample, the argument of Eq. (18) is proportional to Φ(xs−psd1,ys)−Φ(xs+psd1,ys)2, (19) which gives a twin image with an opposite sign separated by 2psd1[29,30,32], as later shown in Fig. 6(a) and (c). Here, we used (b0b−1⁎)⁎=b1b0⁎ ⁠. Since the configuration gives the twin image, it can be regarded as a quasi-two-wave interferometer, but at points where ΔP(x2,y2) defined by ΔP(x2,y2)≡arg[b0b−1⁎Ts(xs,ys)Ts⁎(xs−psd1,ys)]−arg[b1b0⁎Ts(xs+psd1,ys)Ts⁎(xs,ys)], (20) is close to π ⁠, the photon number that contributes to form the twin image becomes very small to make the image have a low signal-to-noise ratio. In addition, when ΔP(x2,y2) exceeds an odd-integer multiple of π ⁠, π-wrapping, i.e. π-phase jump, occurs [32] (see Fig. 7 (a)). This kind of wrapping is different from that caused by the fact that the range of the value of arg[b0b−1⁎Ts(xs,ys)Ts⁎(xs+psd1,ys)+b1b0⁎Ts(xs−psd1,ys)Ts⁎(xs,ys)] is limited from –π to π ⁠; unwrapping algorithms for the latter do not work well for the former because we can experimentally obtain π-wrapped images after they are blurred by the finite size of the X-ray source and/or point spread function of the X-ray image detector. It is often assumed that Φ(xs,ys) is a slowly varying function of xs ⁠. In this case, Φ(xs,ys) can be expanded by a Taylor series and approximated up to the first order: Φ(xs,ys−mpsd1)≈Φ(xs,ys)−mpsd1∂Φ(xs,ys)∂xs. (21) Thus, by taking the argument of the 1st order Fourier component of the intensity of the self-image, arg[q1] ⁠, we can obtain a differential-phase image ∂Φ(xs,ys)∂xs ⁠. Note that there is a tradeoff between the magnification of the sample and the sensitivity to ∂Φ(xs,ys)∂xs ⁠: increasing magnification of the sample makes the slope of the wavefront gentle and, as a result, reduces the sensitivity to ∂Φ(xs,ys)∂xs ⁠. In other words, the effective Talbot order ps becomes smaller when the sample is located closer to the X-ray source to increase the magnification of the sample. There are several setups that have so far been proposed for obtaining the differential-phase image. Figure 3 shows four typical setups of the X-ray grating interferometry for this purpose. In the setup shown in Fig. 3(a) [24,38], the self-image generated downstream of a grating (G1) is resolved by a high-spatial-resolution X-ray image detector, and arg[q1] was obtained by the Fourier transform method [62], where the spatial resolutions of the obtained images are limited by the pitch of the self-image. Because the setup of Fig. 3(a) requires a high-spatial resolution X-ray image detector, which generally has a small field-of-view, it is suitable only in the case where a small field-of-view is sufficient. The interferometer shown in Fig. 3(b) is called a Talbot interferometer [19,63,64], where an absorption grating (G2) is overlaid on a self-image to form an image with moiré fringes (a moiré image). Generating moiré fringes makes it possible to detect the deformation of the self-image without using a high-resolution X-ray image detector and to realize a large field of view. The setups of Fig. 3(c) [35] and (d) [23] enable the use of using a low-brilliance X-ray source. As shown in Section 2, a spatial coherent length that is comparable or larger than the pitch of the grating used is required for the Talbot effect to occur. Because of the weak interaction of hard X-rays with materials, the grating has to have a high aspect-ratio, which makes it difficult to reduce the pitch of the grating to less than 1 μm. Thus, from the van Cittert–Zernike theorem, a micro-focus X-ray source or a very large distance from a large-size X-ray source is necessary to meet the requirement of the spatial coherence length. The setups shown in Fig. 3(c) and (d) meet the requirement with even a low-brilliance normal-size X-ray source: the source grating (G0) located close to the X-ray source works as an array of individually coherent but mutually incoherent X-ray sources. In this case, another condition is necessary for the self-images that is generated by the mutually incoherent X-ray sources to constructively be added: 1d0=αd1−1d2, (22) where α is the compression ratio of the self-image. The use of a normal-size X-ray source drastically reduces the time that is required for obtaining a moiré image because the X-ray flux that is generated by an X-ray source is generally proportional to the size of the source. Fig. 3. View largeDownload slide (Color online) Four typical setups for X-ray grating interferometry where one (a), two ((b) and (c)), and three (d) gratings are used. Setups illustrated in (b), (c), and (d) are referred to as Talbot, Lau, and Talbot–Lau types. Fig. 3. View largeDownload slide (Color online) Four typical setups for X-ray grating interferometry where one (a), two ((b) and (c)), and three (d) gratings are used. Setups illustrated in (b), (c), and (d) are referred to as Talbot, Lau, and Talbot–Lau types. It should be noted that, in addition to the merit that continuous-spectrum spherical-wave X-rays are available, the X-ray grating interferometry has multi-modality: it can provide three independent images called transmittance, differential-phase, and visibility-contrast images [25,31,36,37,41–47], which can be given by |q0| ⁠, arg[q1] ⁠, and 2|q1||q0| ⁠, respectively. The first image corresponds to the conventional absorption-contrast image, while the contrast of the third image has several origins including mainly ultra-small-angle X-ray scattering from the sample, which enables structure analysis at each pixel. X-ray imaging microscopy using self-imaging phenomenon In Section 2, we saw that the Talbot effect can occur even when a grating is illuminated by continuous-spectrum spherical-wave X-rays and that it enables to realize X-ray phase-contrast imaging. Since a spherical wave can magnify an image, we can construct an X-ray phase-contrast microscope by using the Talbot effect. The simplest type of X-ray phase-contrast microscopy using the Talbot effect is the projection types that have been shown in Fig. 3. In these cases, their spatial resolutions are determined not only by magnification but also by the total size of the X-ray source and the spatial resolution of the X-ray image detector. To attain a spatial resolution of, e.g. 1 μm using the projection type X-ray phase-contrast microscopy, it is necessary to use an X-ray source with a size <1 μm and/or an X-ray image detector with a spatial resolution <1 μm, both of which reduce the photon number available for X-ray imaging. In addition, free-space propagation of X-rays from the sample to the X-ray image detector blurs the image of the sample projected on the detector because of refraction caused by the sample, although this blurring causes an edge-enhancement effect, which is positively used for the so-called propagation-based imaging technique. The setup of the X-ray imaging microscope shown in Fig. 4(a), consisting of an X-ray lens and a grating, can solve the problem of the defocusing. For simplicity, we assume here that the X-ray lens (objective lens) is ideal, i.e. the amplitude spread function (ASF) is given by a delta function. In fact, from detailed analytical calculations, it was shown that the electric field created on the image plane (the (x'2,y'2)-plane) can be expressed by E2(x′2,y′2)≈E′20∑nbn(p′)Ts(x′s−np′d′1,−y′s)×exp[2πinx′2d′2], (23) which is the same as that in the case of the X-ray projection microscope shown in Fig. 3(a) except that there is no defocusing effect. Here, E′20 corresponds to the electric field on the image plane without the sample and the grating, expressed by E′20=−exp[2πi(a+b)λ]exp[πi(x2′2+y2′2)λR2′]M, (24) where a ⁠, b and R′2 are geometrically defined distances from the object plane to the X-ray lens, from the X-ray lens to the image plane, and from the focal spot of the X-ray source to the image plane, respectively, M is the magnification of the sample, which is given by ba ⁠, and the (x′s,y′s)-plane is defined on the object plane, defined by (x′s,y′s)≡−1M(x′2,y′2) ⁠. Equation (23) includes both the cases where the grating is located upstream and downstream of the X-ray lens if we define the Talbot order p' and the pitch d'2 of the Fresnel image on the image plane as follows: p′≡{−λzeff,s1d′12(Rs≤R1