X-ray wavefront characterization using a rotating shearing interferometer technique Hongchang Wang,1 Kawal Sawhney,1 Sébastien Berujon,1,2 Eric Ziegler,2 Simon Rutishauser,3 and Christian David3 1
Diamond Light Source, Harwell Science and Innovation Campus, Didcot, Oxfordshire, OX11 0DE, UK 2 European Synchrotron Radiation Facility, BP-220, F-38043 Grenoble, France 3 Laboratory for Micro- and Nanotechnology, Paul Scherrer Institut, 5232 Villigen PSI, Switzerland *
[email protected]
Abstract: A fast and accurate method to characterize the X-ray wavefront by rotating one of the two gratings of an X-ray shearing interferometer is described and investigated step by step. Such a shearing interferometer consists of a phase grating mounted on a rotation stage, and an absorption grating used as a transmission mask. The mathematical relations for X-ray Moiré fringe analysis when using this device are derived and discussed in the context of the previous literature assumptions. X-ray beam wavefronts without and after X-ray reflective optical elements have been characterized at beamline B16 at Diamond Light Source (DLS) using the presented X-ray rotating shearing interferometer (RSI) technique. It has been demonstrated that this improved method allows accurate calculation of the wavefront radius of curvature and the wavefront distortion, even when one has no previous information on the grating projection pattern period, magnification ratio and the initial grating orientation. As the RSI technique does not require any a priori knowledge of the beam features, it is suitable for routine characterization of wavefronts of a wide range of radii of curvature. ©2011 Optical Society of America OCIS codes: (050.2770) Gratings; (340.7450) X-ray interferometry; (120.3940) Metrology; (120.4640) Optical instruments.
References and links 1.
E. Ziegler, L. Peverini, I. V. Kozhevnikov, T. Weitkamp, and C. David, “On-line mirror surfacing monitored by X-ray shearing interferometry and X-ray scattering,” AIP Conf. Proc. 879, 778–781 (2007). 2. C. David, B. Nöhammer, H. H. Solak, and E. Ziegler, “Differential x-ray phase contrast imaging using a shearing interferometer,” Appl. Phys. Lett. 81(17), 3287–3289 (2002). 3. T. Weitkamp, B. Nöhammer, A. Diaz, C. David, and E. Ziegler, “X-ray wavefront analysis and optics characterization with a grating interferometer,” Appl. Phys. Lett. 86(5), 054101–054103 (2005). 4. A. Momose, S. Kawamoto, I. Koyama, Y. Hamaishi, K. Takai, and Y. Suzuki, “Demonstration of X-Ray Talbot Interferometry,” Jpn. J. Appl. Phys. 42(Part 2, No. 7B), L866–L868 (2003). 5. 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(4), 258–261 (2006). 6. T. Weitkamp, A. Diaz, C. David, F. Pfeiffer, M. Stampanoni, P. Cloetens, and E. Ziegler, “X-ray phase imaging with a grating interferometer,” Opt. Express 13(16), 6296–6304 (2005). 7. F. Pfeiffer, O. Bunk, C. Schulze-Briese, A. Diaz, T. Weitkamp, C. David, J. F. van der Veen, I. Vartanyants, and I. K. Robinson, “Shearing interferometer for quantifying the coherence of hard x-ray beams,” Phys. Rev. Lett. 94(16), 164801 (2005). 8. T. Weitkamp, A. Diaz, B. Nöhammer, F. Pfeiffer, M. Stampanoni, E. Ziegler, and C. David, “Moire interferometry formulas for hard x-ray wavefront sensing,” Proc. SPIE 5533, 140–144 (2004). 9. A. Diaz, C. Mocuta, J. Stangl, M. Keplinger, T. Weitkamp, F. Pfeiffer, C. David, T. H. Metzger, and G. Bauer, “Coherence and wavefront characterization of Si-111 monochromators using double-grating interferometry,” J. Synchrotron Radiat. 17(3), 299–307 (2010). 10. J.-P. Guigay, S. Zabler, P. Cloetens, C. David, R. Mokso, and M. Schlenker, “The partial Talbot effect and its use in measuring the coherence of synchrotron X-rays,” J. Synchrotron Radiat. 11(6), 476–482 (2004). 11. M. Born and E. Wolf, Principles of Optics (Pergamon Press, 1993).
#149463 - $15.00 USD
(C) 2011 OSA
Received 17 Jun 2011; revised 25 Jul 2011; accepted 10 Aug 2011; published 12 Aug 2011
15 August 2011 / Vol. 19, No. 17 / OPTICS EXPRESS 16550
12. K. J. S. Sawhney, I. P. Dolbnya, M. K. Tiwari, L. Alianelli, S. M. Scott, G. M. Preece, U. K. Pedersen, R. D. Walton, R. Garrett, I. Gentle, K. Nugent, and S. Wilkins, “A test beamline on diamond light source,” AIP Conf. Proc. 1234, 387–390 (2010). 13. C. David, J. Bruder, T. Rohbeck, C. Grünzweig, C. Kottler, A. Diaz, O. Bunk, and F. Pfeiffer, “Fabrication of diffraction gratings for hard X-ray phase contrast imaging,” Microelectron. Eng. 84(5-8), 1172–1177 (2007). 14. I. Zanette, T. Weitkamp, T. Donath, S. Rutishauser, and C. David, “Two-dimensional x-ray grating interferometer,” Phys. Rev. Lett. 105(24), 248102 (2010). 15. H. Itoh, K. Nagai, G. Sato, K. Yamaguchi, T. Nakamura, T. Kondoh, C. Ouchi, T. Teshima, Y. Setomoto, and T. Den, “Two-dimensional grating-based X-ray phase-contrast imaging using Fourier transform phase retrieval,” Opt. Express 19(4), 3339–3346 (2011).
1. Introduction The combination of at-wavelength (i.e. using x-rays) and in situ metrology methods is considered the best pathway for surpassing the present optics performance and to realize the true nano-focusing optical elements needed for future synchrotron experiments [1]. Nevertheless, the need to push the resolution and sensitivity to nanometer scale still requires improvements that are the object of significant attention from the X-ray community. Of the many methods that exist in visible optics for optical components metrology, only a few are applied to X-Rays. One such method is the grating based shearing interferometer (SI), nowadays widely used by the synchrotron community for X-Ray wavefront characterization, tomography and radiography [1–6]. Indeed, the X-ray shearing interferometer presents the advantages of compactness, robustness, and a low requirement on longitudinal coherence [7]. The instrument consists of a first phase grating as a beam splitter and a second absorption grating acting as a transmission mask. The combined effects of the two gratings on an X-ray beam produce Moiré fringes that can be resolved by a 2D detector, and the calculation of their distortion from the straight lines allows one to characterize wavefront aberrations in the order of the wavelength [3]. A common method used with the SI is the phase stepping method in which at least one period of a grating is scanned while the two gratings are set at the same angle with respect to the beam. To reduce the artifact effects when using the phase stepping method for radiography and tomography, it is important to consider the divergence of the beam and accordingly use gratings with perfectly matching periods [6]. However, for routine wavefront characterization, changing the gratings several times to select the adequate one is time consuming, and then, this becomes a limiting factor when working with different optics or different configurations on a beamline. Moreover, the careful calibration needed in the phase stepping method to set the gratings parallel to each other is a stringent requirement that can lead to non- negligible systematic errors. The Moiré fringes analysis method for wavefront characterization is based on the assumption that the two gratings are inclined symmetrically in the opposite direction [8], leading usually to unexpected variations in the wavefront curvature measurements [9]. We show that these limitations are overcome when applying the rotation method: from several images collected at different rotation angles, the relevant grating angles are accurately extracted, and then the X-ray wavefront radius of curvature and its distortions can be obtained from any of these images. 2. Theoretical considerations The pitch of the phase grating and the absorption grating are respectively noted d0 and d2, and the period of the interference fringes created by the phase grating, at the position of the absorption grating, is noted d1. One can define the ratio η of the two grating pitches by:
η=
d0 d2
(1)
Working with a collimated beam, the Talbot distance [10] beyond the π-shifting phase grating is defined by
#149463 - $15.00 USD
(C) 2011 OSA
Received 17 Jun 2011; revised 25 Jul 2011; accepted 10 Aug 2011; published 12 Aug 2011
15 August 2011 / Vol. 19, No. 17 / OPTICS EXPRESS 16551
Lm = m
d 02 8λ
(2)
where λ is the wavelength, and m an odd integer which corresponds to the order of the fractional Talbot distance. For a diverging beam with a radius of curvature R, these distances rescale to [11]:
L*m =
R Lm R − Lm
(3)
The interference pattern produced by the phase grating at the distance Lm* where the second grating is located, is magnified from the expected self-referenced Talbot periodic grating pattern, by a factor:
M =
L*m + R R
(4)
And so the relation linking d1 to d2 can be written as
d1 =
η d2 2
M =
d2
κ
(5)
where κ is the revised demagnification factor defined as:
κ=
2R ( L*m + R)η
(6)
So far, the mathematical relations used in the literature for hard X-ray wavefront analysis of Moiré fringes were derived assuming that the small inclination angle β of the phase grating from the detector axis is equal to the inclination angle α of the absorption grating [8]. As illustrated in Fig. 1, the observed Moiré fringes on the detector make an angle θ with the vertical detector axis. In the illustration, the ratio η of the two grating pitch is set to 1.3, and the two gratings are inclined to α = 17.2° and β = −5.7°, respectively. The periods of the recorded Moiré fringe along the x and y axis directions of the detector are noted dx and dy.
#149463 - $15.00 USD
(C) 2011 OSA
Received 17 Jun 2011; revised 25 Jul 2011; accepted 10 Aug 2011; published 12 Aug 2011
15 August 2011 / Vol. 19, No. 17 / OPTICS EXPRESS 16552
Fig. 1. Moiré fringes formation. The Moiré fringes will have an inclination angle θ with respect to the vertical axis when the two gratings are tilted from the horizontal axis by angles β and α, respectively
The fundamental harmonics of the interference pattern intensity created by the phase grating I1(x, y), and the intensity transmitted through the absorption grating I2(x, y) can be written as:
2 π 2 I1 ( x, y ) = cos ( y cos β − x sin β ) ≡ cos B d1 2 π 2 I 2 ( x, y ) = cos d ( y cos α − x sin α ) ≡ cos A 2
(7)
where B and A respectively represent the quantities in square brackets in the above equation. Following the methodology of reference [8], we calculate the intensity distribution on the detector which is equal to the product of I1(x, y) and I2(x, y).
I ( x, y ) = I1 ( x, y ) × I 2 ( x, y ) = cos 2 B × cos 2 A
1 2 2 = cos(2 B + 2 A) + cos(2 B − 2 A) + 4 cos + 4 cos − 2 B A 8 I1 ( x , y ) I2 ( x, y ) Im ( x, y )
(8)
The Moiré fringes information recorded on the detector is generated by the term Im(x, y):
2π 2π I m ( x, y ) = cos ( y cos β − x sin β ) − ( y cos α − x sin α ) d2 d1
(9)
Substituting the expression of d1 from Eq. (5), Eq. (9) can be written as:
#149463 - $15.00 USD
(C) 2011 OSA
Received 17 Jun 2011; revised 25 Jul 2011; accepted 10 Aug 2011; published 12 Aug 2011
15 August 2011 / Vol. 19, No. 17 / OPTICS EXPRESS 16553
2π 2π I m ( x, y ) = cos ( sin α − κ sin β ) x − ( cos α − κ cos β ) y d2 d2
(10)
As described in reference [8], the inclination angle θ of the Moiré fringes (see Fig. 1) can be written as:
tan θ =
cos α − κ cos β sin α − κ sin β
(11)
and the components dx and dy of the Moiré fringe period along x and y are:
d2 d2 d x = sin α − κ sin β ≡ γ x d d 2 d = ≡ 2 y cos α − κ cos β γ y
(12)
γ x = sin α − κ sin β γ y = cos α − κ cos β
(13)
where
For each horizontal line j in one image, the intensity of the fundamental component of the Moiré fringes along the horizontal x direction can be written as:
2π I j ( x) = a0 + b0 cos x +ϕ dx
(14a)
Using the Fourier method [3], the phase φ as well as the period dx can be extracted for each row j. The phase along the same fringe is constant c, that is
2π x +ϕ = c dx
(14b)
Therefore the inclination θ of the fringes can be derived as:
tan θ =
d ∂ϕ ∂x =− x ∂y 2π ∂y
(15)
Once θ and dx are derived, the parameters γx and γy can be calculated using Eqs. (11) and (12). Although the parameters α and κ can be determined from Eq. (13) using only two different images with two different angles β, a much better way to proceed is to extract both at the same time by fitting γy as a function of γx when varying the phase grating angle β. One can deduce the relation between γx and γy from Eq. (13):
γ y = cos α − κ 2 − ( γ x − sin α )
2
(16)
Here, the parameters γx, γy and κ are average values over one image. Once the angle α has been calculated, Eq. (16) can be rewritten as:
κ ( y) =
#149463 - $15.00 USD
(C) 2011 OSA
(γ
( y ) − cos α ) + ( γ x ( y ) − sin α ) 2
y
2
(17)
Received 17 Jun 2011; revised 25 Jul 2011; accepted 10 Aug 2011; published 12 Aug 2011
15 August 2011 / Vol. 19, No. 17 / OPTICS EXPRESS 16554
The revised demagnification ratio κ(y) is determined for each line j. Once η and κ are known, the radius of curvature can be determined from the rearranged form of Eq. (6):
L*mηκ ( y ) 2 − ηκ ( y )
R( y ) =
(18)
Finally, the wavefront slope can be recovered by integration of the inverse of the radius of curvature R−1 along the vertical direction y [8]. y
S y = ∫ R ( y )−1 dy 0
(19)
One can see that both γx and γy depend on the phase grating angle β. The numerical extraction of β is however not required for the wavefront slope and radius of curvature reconstruction process. 3. Experiments and results The experiments were carried out at the Diamond Light Source beamline B16 [12] where Xrays are produced by a bending magnet on a 3 GeV storage ring. A schematic of the rotating shearing interferometer at B16 beamline is given in Fig. 2. The phase grating G1 was made from a silicon wafer with a designed pitch of d0 = 4.0 µm. The absorption grating G2 with a design pitch d2 = 2.0 µm was made by covering the lines of a silicon grating with electroplated gold [13]. The rotating shearing interferometer was mounted on a versatile optics test bench, with each grating placed on a motorized tower capable of independent three axis translations. The distance between the two gratings was changed by moving the absorption grating stages. The distance from the source to the first grating G1 was 46.5 m and the angle of this grating was changed using a rotation motor. The phase grating was designed to produce a phase shift of π rad at the X-ray energy of 14.8 keV, thus allowing maximum fringe visibility. The exact energy value was selected by a silicon double-crystal monochromator (DCM). An X-ray 2D CCD detector with indirect illumination and with an effective pixel size of 6.4 µm was used to record the Moiré fringe pattern. The wavefronts before and after insertion into the beam of a focusing curved test mirror were characterized using the presented RSI method.
Fig. 2. A schematic of the rotating shearing interferometer at B16 beamline.
In the first instance, the direct beam wavefront was characterized by retracting the test mirror from the beam path. An inter-grating distance L*m = 0.075m, corresponding to the 3rd order Talbot distance, was employed, which gave an average visibility [9] of 0.45 for the Moiré fringes. Figure 3 shows four interferograms recorded at different phase grating tilt angles β. As expected, the inclination angle θ and the horizontal period dx of the Moiré fringes vary with this angle. Interferograms were taken at various odd Talbot orders (distance between gratings), and up to 25th Talbot orders were used. The absorption grating tilt angle α was kept the same for all data sets. For each inter-grating distance L*m, the average γx and γy values were extracted and fitted using Eq. (16). #149463 - $15.00 USD
(C) 2011 OSA
Received 17 Jun 2011; revised 25 Jul 2011; accepted 10 Aug 2011; published 12 Aug 2011
15 August 2011 / Vol. 19, No. 17 / OPTICS EXPRESS 16555
Fig. 3. Interferograms recorded at four different angles: (a) β = 0.019 rad, (b) β = 0.025 rad, (c) β = 0.041 rad, (d) β = 0.048 rad.
Figure 4 shows some data collected at the 3rd, 9th, 15th and 21st Talbot order distances. For each Talbot order, Eq. (16) was used to calculate the inclination angle α of the absorption grating and the revised demagnification ratio κ. The good fitting of the experimentally extracted points γy as a function of γx confirm the validity of the model described by Eq. (16). Only some values close to γx = 0 do not fit well. The reason for this is that for γx ≈0 the Moiré fringes are almost horizontal, rendering the phase extraction complex and inaccurate. γx ≈0 corresponds to the two gratings being parallel to each other i.e. β = α. RSI therefore provides a quantitative method to more accurately find the required zero angle of β and α, for use in the phase stepping method.
Fig. 4. Experimental data and fitted model from Eq. (16) for γy as function of γx at the 3rd, 9th, 15th and 21st order Talbot distances. The points corresponding to the four images shown in Fig. 2 are circled.
#149463 - $15.00 USD
(C) 2011 OSA
Received 17 Jun 2011; revised 25 Jul 2011; accepted 10 Aug 2011; published 12 Aug 2011
15 August 2011 / Vol. 19, No. 17 / OPTICS EXPRESS 16556
The calculated α and κ values for the 3rd, 9th, 15th and 21st Talbot orders are tabulated in Table 1. The approximations d0 = 4.0 µm and d2 = 2.0 µm (η = d0/d2 = 2) only give approximate radii of curvatures R’, as derived from Eq. (18). The R’ values are also tabulated in Table 1 and are seen to depend on the Talbot order used. For instance, a difference of 4.6 m is observed between the values derived from measurements at the 3rd and 21st orders. Thus, there is a need to know the period ratio η very accurately. This has been done by fitting the model for κ described by Eq. (6), with the κ values calculated using Eq. (16), as shown in Fig. 5. A value of η = 1.99975, with an accuracy of 10−5, has been extracted from the fit. It is therefore possible to compare the periods of the two gratings with sub-nanometer accuracy with the RSI technique. The average radii of curvatures R calculated for each image using this new value of η are given in Table 1. The difference in the radii of curvatures between 3rd and 21st orders is only 0.6m now, which comes about from the measurement error (~1mm) of inter-grating distance L*m. For instance, if the inter-grating distance is changed by 1 mm, then both 3rd and 21st orders return a value of 49.7m for the radius of curvature of the wavefront. It can be seen from Table 1 that the α values calculated from different Talbot orders are non-zero but identical, as they should be, since the absorption grating is held fixed, and only β is scanned. If, however, one had assumed the absorption grating angle α to be equal to zero, the radius of curvature R’ would have varied depending on the β angle. For instance, for the 3rd Talbot order, R’ would have changed from 39.2m to 44.1m, as β changed from 0.025 rad to 0.019 rad (images a and b in Fig. 3), thus demonstrating the sensitivity of the radius of curvature calculation to the correct knowledge of the α angle value. The value extracted from the regression of κ (Fig. 5) for the radius of curvature is 49.8 m, which is 5% higher compared to the expected value of 47.5 m. This discrepancy may be due to the heat bump caused by the non-uniform X-ray illumination of the water-cooled DCM. Moreover, we will like to emphasize that these detailed measurements at different Talbot distances, which might seem time consuming, have only to be performed once, in order to determine the exact period ratio η. This value can be then directly used for all subsequent wavefront measurements. Table 1. Radius of Curvature of DCM direct beam at B16 with RSI Order
L*m (m)
α (rad)
κ
R’ (m)
R (m)
21st
0.507
0.0351
0.9900
50.4
49.8
15th
0.363
0.0350
0.9929
50.7
49.8
9th
0.219
0.0350
0.9958
51.3
49.8
3rd
0.075
0.0350
0.9986
55.0
50.4
3rd Fig. 3(a)
0.075
0.0000
0.9981
39.2
36.7
3rd Fig. 3(b)
0.075
0.0000
0.9983
44.1
41.1
3rd
0.074
0.0350
0.9986
54.2
49.7
21st
0.506
0.0351
0.9900
50.3
49.7
#149463 - $15.00 USD
(C) 2011 OSA
Received 17 Jun 2011; revised 25 Jul 2011; accepted 10 Aug 2011; published 12 Aug 2011
15 August 2011 / Vol. 19, No. 17 / OPTICS EXPRESS 16557
Fig. 5. Linear regression of the κ factors calculated at the different Talbot orders. The 3rd, 9th, 15th and 21st Talbot orders are circled.
Once η and α are accurately known, the local radii of curvature R(y) can be recovered using Eqs. (17) and 18. As shown in Fig. 6(1), the calculated radii of curvature of the wavefronts, extracted from the images (a-d) of Fig. 3, are consistent. The mean value of R is 51.4m for image (a) and 51.1m for image (b). The plots show that the wavefronts are not spherical but their curvatures vary with a standard deviation of 9.2 m and 7.4 m respectively. Next, the wavefront slope S(y) was calculated by integration of R−1 along y, as per Eq. (19). The slope error ∆S(y), defined as the difference between S(y) and its best linear fit, is plotted as a function of the vertical position in Fig. 6(2). The standard deviations of the slope error are 0.18 µrad for image (a) and 0.20 µrad for image (b), and hence are in good agreement with each other.
Fig. 6. Data extracted from the images presented in Fig. 3: (1) Radius of curvature as function of y, (2) Wavefront slope error obtained after subtraction of the linear fit. Different line types correspond to calculated data from different images (a-d) of Fig. 3.
#149463 - $15.00 USD
(C) 2011 OSA
Received 17 Jun 2011; revised 25 Jul 2011; accepted 10 Aug 2011; published 12 Aug 2011
15 August 2011 / Vol. 19, No. 17 / OPTICS EXPRESS 16558
Figure 7 illustrates a second application of the RSI method, using the same pair of gratings as for the direct beam study, but now for measuring the beam wavefront reflected from a curved focusing mirror. The average values of γx and γy were extracted from different images taken during a rotation scan of the phase grating and the same procedure as described previously was applied. This time the demagnification ratio κ was determined to be 0.9723 and the wavefront radius of curvature to be 3.3 m. The wavefront from the focusing mirror has a strong curvature, and its precise retrieval by the RSI therefore shows that the wavefronts can be characterized by RSI even when the pitches of the two gratings are far from satisfying the criteria of divergence matching.
Fig. 7. Experimental data and its fitting to Eq. (16): γy as function of γx at the 3rd Talbot order distance (L = 93 mm) for a wavefront reflected from a focusing mirror.
7. Conclusions The rotating shearing interferometer method has been established at DLS B16 beamline enabling different wavefronts to be accurately characterized. The experimental results are in good agreement with the physical expectations. Full wavefront characterization can be determined from a single image, once the transmission grating angle offset (α) is known (derived from a simple rotation scan). The RSI technique allows accurate derivation of the projected fringes period, without any previous knowledge on the beam characteristics. Unlike the phase stepping method, the RSI technique with the 1D interferometer provides averaged information in one direction, it however relaxes the requirement of divergence matching of the two gratings. The RSI method has been demonstrated to be routinely usable for X-ray wavefront characterization and at-wavelength metrology for beams with various radii of curvatures. Recently, genuine 2D gratings have been successfully fabricated and used for twodimensional Talbot interferometry [14,15]; allowing the Fourier analysis to be performed in two orthogonal planes simultaneously. The presented RSI method will also work with such 2D-interferometers. Acknowledgments This work was carried out with the support of the Diamond Light Source Ltd UK. We thank Timm Weitkamp from Synchrotron Soleil for useful discussions. We also gratefully acknowledge Andrew Malandain for his support during the setting up of the experiment. The research leading to these results has received partial funding from the European Community’s Seventh Framework Programme (FP7/2007-2013) under the grant agreement nº 226716.
#149463 - $15.00 USD
(C) 2011 OSA
Received 17 Jun 2011; revised 25 Jul 2011; accepted 10 Aug 2011; published 12 Aug 2011
15 August 2011 / Vol. 19, No. 17 / OPTICS EXPRESS 16559