The reduction in visibility in x-ray grating interferometry based on the Talbot effect is formulated by the autocorrelation function of spatial fluctuations of a wavefront due to unresolved micron-size structures in samples. The experimental results for microspheres and melamine sponge were successfully explained by this formula with three parameters characterizing the wavefront fluctuations: variance, correlation length, and the Hurst exponent. The ultra-small-angle x-ray scattering of these samples was measured, and the scattering profiles were consistent with the formulation. Furthermore, we discuss the relation between the three parameters and the features of the micron-sized structures. The visibility-reduction contrast observed by x-ray grating interferometry can thus be understood in relation to the structural parameters of the microstructures.
© 2010 Optical Society of America
Since their discovery in 1895 by Röntgen, x-rays have been utilized for visualizing internal structures non-destructively in a wide range of fields including those of physics, chemistry, biology, and medicine. The absorption of x-rays has been conventionally used for imaging applications, but in the early 1990s several imaging techniques using x-ray phase shift were proposed. These so-called x-ray phase imaging techniques [1, 2] have been highlighted because they can have a sensitivity to light elements three orders of magnitude higher than that of absorption-contrast imaging. However, most of these techniques have essentially required a highly brilliant synchrotron x-ray source. A great deal of attention has recently been paid to x-ray phase imaging techniques using transmission gratings [3, 4, 5, 6, 7, 8, 9, 10, 11, 12], represented by Talbot and Talbot-Lau interferometry, because they work even with polychromatic and cone beam from a compact laboratory source.
We can retrieve two kinds of quantitative images with x-ray Talbot (-Lau) interferometry, i.e. absorption and differential-phase, from a series of experimentally obtained moiré images. Pfeiffer et al. has recently proposed another approach to forming image contrast, where a relative reduction in the visibility of the moiré image is quantified by defining normalized visibility . They reported that this visibility contrast was formed through the mechanism of small-angle x-ray scattering from microstructures with a scale much smaller than the spatial resolution of the imaging system. Their approach is fascinating because it can provide structural information that is inaccessible from the absorption and differential-phase images, and it shows promise of offering a broad range of applications [14, 15, 16, 17, 18, 19, 20]. However no general formulation of the phenomenon, which is essential for quantitative structure analysis, has thus far been provided.
This paper explains how the reduction in the visibility can generally be formulated by an autocorrelation function describing the spatial fluctuations of a wavefront due to the microstructures. Visibility reduction obtained experimentally with microspheres and a melamine sponge was successfully explained by the formula. The experimental results were also consistent with results from ultra-small angle x-ray scattering measurements in the far field. Three parameters are required to characterize the wavefront fluctuations: variance, correlation length, and the Hurst exponent. We demonstrate from analytical and numerical calculations that the correlation length represents the average size of the microstructures, while the Hurst exponent is relevant to their shape. The results of least-squares fitting to the experimental data were well explained by this interpretation. Thus, the visibility contrast is quantitatively related to the structural parameters of the microstructures. Our approach can directly provide two- and three-dimensional structural information in real space, a scale that is only accessible by ultra-small-angle x-ray scattering (the near-field speckle region ), and it is expected to be broadly applied to medical, biological, and material sciences.
2. Theoretical description of the reduced visibility
First, we derive a general formula for the reduced visibility on the basis of the spatial fluctuations of a wavefront. Let us consider the x-ray Talbot interferometer system schematically illustrated in Fig. 1. Here, we assume that a plane wave with wavelength λ is illuminating a sample located just in front of the first grating; the following can easily be extended to the case of polychromatic- and spherical-wave illumination .
The complex refractive index of the sample is expressed by 1 − δ + iβ, and the electric field just behind the sample can be written in the projection approximation by E 0 exp[−α(x,y)/2]exp[−iΦ(x,y)], where E 0 is the amplitude of the incident x-rays, and α and Φ are given by (4π/λ) ∫ β (x,y,z)dz and (−2π/λ) ∫ δ (x,y,z)dz. To discuss the effect of phase fluctuations, we present α and Φ as superpositions of smooth (resolvable) and fine (unresolvable) features , α = α s + α f and Φ = Φs + Φf. In the following, we assume that α f is negligible. We also assume, for the sake of simplicity, that the size of the first grating is infinite.
The electric field just in front of the second grating, located at distance z T from the first grating, can be given in the paraxial approximation by 
Here, a′n is the nth Fourier coefficient of the electric field generated in front of the second grating (n = 0,±1, …), given by an exp(−iπpn 2), where an is the nth Fourier coefficient of the amplitude transmission function of the first grating and p is the Talbot order , d is the pitch of the gratings, and Φn(x,y) is defined by Φ(x − npd,y).
The intensity just behind the second grating detected by a pixel located at point (x,y) is approximated by
Here, I 0 ≡ ∣E 0∣2, cN is the Nth Fourier coefficient of the intensity transmission function of the second grating (N = 0,±1, …), µm is the complex coherence factor of x-rays  at two points separated by distance mpd on the first grating (m = 0,±1, …), χ is the relative displacement of the second grating to the first grating in the x direction, and PSF(x,y) is the normalized point spread function (PSF) of the detector. Note that, in Eq. (2), we have assumed that the two gratings are parallel to the y axis without loss of generality and that the scale resolved by the image system is sufficiently larger than the pitch of the gratings.
The visibility is approximately proportional to the ratio of the modulus of the 1st order Fourier coefficient term, q 1, in Eq. (2) (corresponding to N = 1) to that of the 0th order, q 0, (corresponding to N = 0), both of which are experimentally obtained by using the fringe-scanning technique . The 0th order Fourier coefficient, q 0, is given by
The factor, Ã n,m,0, is rewritten by
where Φf,n(x,y) ≡ Φf(x − npd,y). Here, we have considered that the resolvable feature, Φs, is a slowly varying function of x and that the derivative of Φs higher than the first can be neglected. If we assume that the unresolvable features, Φf, are distributed randomly and that the width of PSF is sufficiently larger than the characteristic scale of the unresolvable features, Ã n,m,0(x,y) can be further approximated by 
Here, the bar means the averaging around (x,y):
where D is the width of the PSF. Since the phase Φf,n(x,y) − Φf,n−m(x,y) in Eq. (7) is random, the terms of Ã n,m,0(x,y) for m ≠ 0 are much smaller than Ã n,0,0(x,y). Hence,
Thus, ∣q 0∣ is finally given by
where C 0 ≡ ∑nμ 0∣a′n∣2 c 0.
We can similarly obtain an approximate form of ∣q 1∣. From Eq. (2),
Since the phase, Φf,n(x,y) − Φf,n−m(x,y), in Eq. (15) is random, the terms, Ã n,m,1(x,y), for m ≠ −1 are much smaller than Ã n,−1,1(x,y). Hence,
If a grating with a Ronchi ruling is used as the first grating, the even order Fourier coefficients of the amplitude transmission function of the grating except for the 0th order vanish. Then
Expressing a′0 a′* 1 by |a′0 a′* 1|exp[iΩ], and using the fact a′−1 a′* 0 = (a′0 a′* 1)*, we can finally obtain an approximate form of ∣q 1∣:
where C 1 ≡ 2∣μ −1∣∣c 1∣∣a′0 a′* 1∣ cosΩ and ∆Φf(x,y; pd) ≡ Φf,0(x,y) − Φf,1(x,y). Here, we have assumed that the average phase, , is zero, and that
because D was assumed to be much larger than pd.
Hence, the normalized visibility, which is the ratio of visibility with the sample (V) to that without the sample (V 0), is given by
In the derivation of Eq. (22), we assumed that Φf could be modeled as a random Gaussian process [23, 27], the width of which is given by σ. The term, γ, is the normalized autocorrelation function given by
As wavefront fluctuations can be treated as surface roughness, we can use the simplest general model for γ that was proposed by Sinha et al. to deal with height fluctuations on surfaces :
where H is the Hurst exponent (0<H <1) and ξ is the correlation length of phase fluctuations. Thus, we can generally formulate the reduced visibility in terms of the autocorrelation function of spatial fluctuations of a wavefront. Note that, from Eq. (22), V/V 0 approaches zero when σ increases. In addition, for a given sample (with a given ‘scattering power’), V/V 0 is a function of pd.
3. Verification of formulation through experiments
3.1. Dependence of V/V0 on pd
We measured the dependence of V/V 0 on pd to corroborate Eq. (22) by using monochromatic x-rays. The experiment was carried out with synchrotron x-rays at the beamline 14C, Photon Factory (PF), Japan. Cross-linked PMMA microspheres (EPOSTAR MA, Nippon Shokubai Co., Ltd.) with various radii and a melamine sponge were used as the samples. Microspheres with a fixed radius were dispersed in glycerin liquid, which was dispensed into a 10 mm-thick plastic cell. The volume fraction of the microspheres was fixed at 0.057. The melamine sponge was 1.6 mm thick. Gold gratings with a Ronchi ruling (a π/2-phase grating for the first grating and a 30 µm-thick absorption grating) were used. They were aligned parallel to each other. The images were recorded using a charge-coupled device (CCD)-based x-ray image detector (Spectral Instruments), where the CCD (4096 × 4096 pixels) was connected to a 40 µm GOS screen with a 2:1 fiber coupling. The effective pixel size was 18 µm, and the width of the line spread function was 70 µm.
The plots in Fig. 2 (a) represent the experimental results for the microspheres. The wavelength of the x-rays was fixed at 0.71 Å. Here, 5.3 µm-pitch gratings were used and the Talbot order, p, was changed by changing the distance between the first and second gratings. The curves in the figure plot the results of least-squares fitting to the experimental data by using Eqs. (22) and (24). The experimental data are in good agreement with the fitting curves. This means that the description of the reduced visibility obtained by using Eqs. (22) and (24) is valid.
Figure 2(b) plots the experimental results for the melamine sponge, where the results for three wavelengths and two pitches of the gratings are shown. The experiments were performed at the beamline 14C, PF, and at the beamline 20XU, SPring-8, Japan. At SPring-8, an x-ray camera consisting of a phosphor screen (10 µm-P43, Gd2O2S:Tb+ fine powders), a relay lens, and a CCD camera (Hamamatsu Photonics C4742-98-24A, 1344 × 1024 pixels) was used. The effective pixel size of the detector was 3.14 µm. The open and filled circles in Fig. 2 (b) plot the results for λ = 1.0 and 0.5 Å (SPring-8), while the triangles and squares plot those for another sample with the same thickness for d = 5.3 and 8.0 µm at 0.7 Å (PF). At SPring-8, 20 × 20 pixel binning was implemented to make D almost the same as that of the detector used at PF. The good agreement between the experimental results for the 5.3 µm and 8.0 µm pitches corroborates that Eq. (22) is correct. In addition, the experimental data were well fitted by using Eqs. (22) and (24) for all the wavelengths. Note that the results did not depend on the direction of the sample. The fitting results in Figs. 2(a) and 2(b) are summarized in Table 1.
3.2. Ultra-small-angle scattering measurement
The angular distribution of small-angle scattering in the far field should be related to the Fourier transform of Eq. (21). Hence, we can also check the validity of Eq. (22) through small-angle x-ray scattering measurements. The experiment was performed using Mo Kα characteristic x-rays from a laboratory x-ray source. The non-dispersive double-crystal arrangement seen in Fig. 3(a) was employed. The intensity, Ix, integrated in the y-direction was measured with this setup, which should agree with the Fourier transform of Eq. (22) in the x-direction. The plots in Fig. 3(b) represent the measured intensity of the (ultra) small-angle scattering for the microspheres samples. The curves in Fig. 3(b) are the Fourier transforms of the V/V 0 of the fitting curves in Fig. 2(a) convoluted by the measured rocking curve, the width of which was 7 µrad. The good agreement between the plots and curves substantiates the validity of Eq. (21).
4. Relation between fitting parameters and structural parameters of sample
Finally, we will explain how the parameters that characterize the spatial fluctuations of a wavefront are related to the structural parameters of a sample. First, let us discuss σ and γ. For simplicity, we approximate the sample by a matrix containing a large number of unresolved microscopic particles. If the positions and directions of the particles are uncorrelated in Eq. (23), there is also no correlation between the projections of the particles. Then, γ(x,y;∆x) in Eq. (22) can simply be given by the autocorrelation function of the projection of a particle. In addition, since γ is normalized as it satisfies γ(x,y;0) = 1, γ is determined only by the shape of a particle. Other parameters such as the number density, N, of particles and the thickness, T, of the sample are contained in σ 2. For example, in microspheres with a radius of a, γ and σ 2 can be analytically given as:
where ∆x′ ≡ ∆x/(2a), ∆ρ is the number density of electrons, and r e is the classical electron radius. Here we assumed that TNπa 2 ≫ 1. These results were also confirmed by numerical calculations for a system where microspheres disperse. It should be noted that the first factor on the right hand side of Eq. (26) corresponds to the number of trials in the random walk problem, while the second factor originates from the phase shift by a microsphere. Since the number of trials is proportional to N and T and the phase shift by a particle is proportional to λ far from an absorption edge, σ 2 and ln(V/V 0) are proportional to TNλ 2 for any homogeneous sample. It can be seen that σ 2 is proportional to λ 2 for the melamine sponge (see Table 1). We also confirmed that ln(V/V 0) is proportional to the thickness of the melamine sponge.
The correlation length, ξ, and the Hurst exponent, H, are related to the average-structure parameters of the microscopic particles. In fact, Eq. (25) can also be fitted by using Eq. (24), resulting in ξ/a = 0.95 and H = 0.97. The former suggests that correlation length ξ represents the average size of the particles. On the other hand H should be relevant to the average shape of the particles. The results of numerical calculations of H for cylinder particles with random directions are indicated in Fig. 4 by the crosses, where ξ and H are plotted against the ratio of the diameter (2a) to the height (h). The plot indicates that H represents the asymmetry of the shape of the particle, i.e., if we define asymmetry factor β by abs (log (2a/h)), H should be approximately given by exp[−0.5β]. Note that, if we define the average radius, ā, of a cylinder by the radius of a sphere with the same volume as the cylinder, the deviation of ξ/ā from unity increases when β increases (open circles in the figure).
Thus, we can relate normalized visibility to the structural parameters of microstructures. Our experimental results showed that H ~ 1 for microspheres, which was consistent with our calculations. On the other hand H ~ 0.7 for the melamine sponge, which corresponds to β ~ 0.7. This was also consistent with its real shape; it has a thin fiber structure (see Fig. 2 (b)). Our results of ξ for microspheres indicate that the possibility of microspheres aggregating increases when their radius is smaller. This is also supported by the results from ultra-smallangle scattering in Fig. 3 because otherwise the experimentally obtained angular broadening, representing the inverse size of microstructures, cannot be explained. The results of ξ for the melamine sponge were also consistent with observations through an optical microscope.
Finally, we discuss the case where the size and shape of the microscopic particles are not uniform in a sample. Let us consider a sample consisting of uniform domains, each of which contains particles with a different size and shape from the others. For simplicity, we assume that each domain size is larger than the size that can be resolved by the detector, and each domain has no structural correlation with the others. From Eqs. (21) and (22), V/V 0 at a point (x,y) on the detector is given by
where j represents the contribution from the jth domain on the path to the point (x,y) along the z-axis. Because σ 2 j should be proportional to the thickness of the jth domain, Eq. (29) has a similar form to the Beer-Lambert law. Hence, −ln [V/V 0] can be given by
where σ and γ are expressed as a function of z. Note that ∂ (σ 2)/∂z is proportional to λ 2 and (∆ρ)2. Thus, we can also carry out tomography and determine the three-dimensional distribution of (1 − γ)∂ (σ 2)/∂z. Obtaining tomograms for ∂(σ 2)/∂z, ξ, and H is also possible by performing scans for at least three pds.
In conclusion, we formulated the reduction in visibility in x-ray Talbot interferometry [Eqs. (22) and (24)]. The formulae quantitatively related the visibility contrast to the spatial fluctuations of a wavefront due to unresolved microstructures. The experimental results for the microspheres and the melamine sponge were consistent with our formulation. Three parameters were required to characterize the wavefront fluctuations: variance, correlation length, and the Hurst exponent. We showed from analytical and numerical calculations that the correlation length represents the average size of microstructures, while the Hurst exponent is relevant to their shape. The results of least-squares fitting to the experimental data (Table 1) were well explained by this interpretation. Thus, the visibility contrast is quantitatively related to the structural parameters of the microstructures.
Defining (1 − γ)∂ (σ 2)/∂z, we can also carry out tomography if γ is not regarded as being dependent on z in each domain. Obtaining tomograms for ∂ (σ 2)/∂z, ξ, and H is also possible by performing scans for at least three pds. Furthermore, we can use two-dimensional gratings to obtain structural parameters in both x and y directions. Our approach can directly provide twoand three-dimensional structural information on unresolved microstructures, and we expect it to be broadly applied to medical, biological, and material sciences.
We appreciate the assistance given by Professor T. Hattori and Dr. D. Noda in the fabrication of the gratings. The experiments using synchrotron radiation were performed at Photon Factory and SPring-8. This study was financially supported by the Japan Science and Technology Agency (JST).
References and links
1. R. Fitzgerald, “Phase-sensitive x-ray imaging,” Phys. Today 53, 23–26 (2000). [CrossRef]
2. A. Momose, “Recent advances in x-ray phase imaging,” Jpn. J. Appl. Phys. 44, 6355–6367 (2005). [CrossRef]
3. A. Momose, S. Kawamoto, I. Koyama, Y. Hamaishi, K. Takai, and Y. Suzuki, “Demonstration of x-ray Talbot interferometry,” Jpn. J. Appl. Phys. 42, L866–L868 (2003). [CrossRef]
4. F. Pfeiffer, T. Weitkamp, O. Bunk, and C. David, “Phase retrieval and differential phase-contrast imaging with low-brilliance x-ray sources,” Nat. Phys. 2, 258–261 (2006). [CrossRef]
5. A. Olivo and R. Speller, “A coded-aperture technique allowing x-ray phase contrast imaging with conventional sources,” Appl. Phys. Lett. 91, 074106 (2007). [CrossRef]
6. F. Pfeiffer, M. Bech, O. Bunk, P. Kraft, E.F. Eikenberry, C.H. Brönnimann, C. Grünzweig, and C. David, “Hard-x-ray dark-field imaging using a grating interferometer,” Nat. Mat. 7, 134–137 (2008). [CrossRef]
8. Y. Takeda, W. Yashiro, T. Hattori, A. Takeuchi, Y. Suzuki, and A. Momose, “Differential phase x-ray imaging microscopy with x-ray Talbot interferometer,” Appl. Phys. Express 1, 117002 (2008). [CrossRef]
9. W. Yashiro, Y. Takeda, and A. Momose, “Efficiency of capturing a phase image using conebeam x-ray Talbot interferometry,” J. Opt. Soc. Am. A 25, 2025–2039 (2008). [CrossRef]
10. A. Momose, W. Yashiro, and Y. Takeda, “X-ray phase imaging with Talbot interferometry,” in Biomedical Mathematics: Promising Directions in Imaging, Therapy Planning and Inverse Problems, edited by Y. Censor, M. Jiang, and G. Wang, (Medical Physics Publishing, Madison, Wisconsin, USA, 2009).
11. Z.-F. Huang, K.-J. Kang, L. Zhang, Z.-Q. Chen, F. Ding, Z.-T. Wang, and Q.-G. Fang, “Alternative method for differential phase-contrast imaging with weakly coherent hard x rays,” Phys. Rev. A 79, 013815 (2009). [CrossRef]
12. A. Olivo, S. E. Bohndiek, J. A. Griffiths, A. Konstantinidis, and R. D. Speller, “A non-free-space propagation x-ray phase contrast imaging method sensitive to phase effects in two directions simultaneously,” Appl. Phys. Lett. 94, 044108 (2009). [CrossRef]
13. W. Yashiro, Y. Takeda, A. Takeuchi, Y. Suzuki, and A. Momose, “Hard x-ray phase-difference microscopy using a Fresnel zone plate and a transmission grating,” Phys. Rev. Lett. 103, 180801 (2009). [CrossRef] [PubMed]
14. C. Grünzweig, C. David, O. Bunk, M. Dierolf, G. Frei, G. Kühne, J. Kohlbrecher, R. Schäfer, P. Lejcek, H. M. R. Ronnow, and F. Pfeiffer, “Neutron decoherence imaging for visualizing bulk magnetic domain structures,” Phys. Rev. Lett. 101, 025504 (2008). [CrossRef] [PubMed]
16. C. Grünzweig, C. David, O. Bunk, M. Dierolf, G. Frei, G. Kuhne, R. Schafer, S. Pofahl, H. M. R. Ronnow, and F. Pfeiffer, “Bulk magnetic domain structures visualized by neutron dark-field imaging,” Appl. Phys. Lett. 93, 112504 (2009). [CrossRef]
17. F. Pfeiffer, M. Bech, O. Bunk, T. Donath, B. Henrich, P. Kraft, and C. David, “X-ray dark-field and phase-contrast imaging using a grating interferometer,” J. Appl. Phys. 105, 102006 (2009). [CrossRef]
18. S.A. MacDonald, F. Marone, C. Hintermüller, G. Mikuljan, C. David, F. Pfeiffer, and M. Stampanoni, “Advanced phase-contrast imaging using a grating interferometer,” J. Synchrotron Radiation 16, 562–572 (2009). [CrossRef]
19. H. Wen, E. E. Bennett, M. M. Hegedus, and S. Rapacchi, “Fourier x-ray scattering radiography yields bone structural information,” Radiography 251, 910–918 (2009).
20. Z.-T. Wang, K.-J. Kang, Z.-F. Huang, and Z.-Q. Chen, “Quantitative grating-based x-ray dark-field computed tomography,” Appl. Phys. Lett. 95, 094105 (2009). [CrossRef]
21. R. Cerbino, L. Peverini, M. A. C. Potenza, A. Robert, P. Bosecke, and M. Giglio, “X-ray-scattering information obtained from near-field speckle,” Nat. Phys. 4, 238–243 (2008). [CrossRef]
22. K. Patorski, “Self-imaging and its applications,” in Progress in Optics XXVII, edited by E. Wolf (ELSEVIER SCIENCE PUBLISHERS B.V., Amsterdam, 1989).
23. Y.I. Nesterets, “On the origins of decoherence and extinction contrast in phase-contrast imaging,” Opt. Commun. 281, 533–542 (2008). [CrossRef]
24. Here, we defined the pth Talbot order as it expresses the pth position where the electric field just behind the first grating is perfectly reproduced for any grating. For plane-wave illumination, the pth Talbot order corresponds to the position of zT = pd2/λ. This definition is convenient because it is independent of what kind of grating is used as the first grating. Note that, once zT is given by pd2/λ, we can use p for specifying any position behind the first grating instead of zT. The analytical calculations presented in this paper can be applied to any periodic image generated behind the first grating at any position of p (> 0).
25. J. W. Goodman, Statistical Optics, (A Wiley-Interscience Publication, New York, 2000).
26. J. H. Bruning, D. R. Herriott, J. E. Gallagher, D. P. Rosenfeld, A. D. White, and D. J. Brangaccio, “Digital Wavefront Measuring Interferometer for Testing Optical Surfaces and Lenses,” Appl. Opt. 13, 2693–2703 (1974). [CrossRef] [PubMed]
27. I. A. Vartanyants and I. K. Robinson, “Origins of decoherence in coherent x-ray diffraction experiments,” Opt. Commun. 222, 29–50 (2003). [CrossRef]
28. S.K. Sinha, E.B. Sirota, and S. Garoff, “X-ray and neutron scattering from rough surfaces,” Phys. Rev. B 38, 2297–2312 (1988). [CrossRef]