research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Multiple-beam X-ray interferometry for phase-contrast microtomography

CROSSMARK_Color_square_no_text.svg

aInstitute of Physics, University of Dortmund, D-44221 Dortmund, Germany, and bHamburger Synchrotronstrahlungslabor (HASYLAB) at Deutsches Elektronen-Synchrotron (DESY), Notkestrasse 85, D-22603 Hamburg, Germany
*Correspondence e-mail: bonse@physik.uni-dortmund.de

(Received 7 September 2000; accepted 2 November 2000)

The first successful operation of an X-ray interferometer under conditions of non-planar three-beam diffraction for phase-contrast X-ray microtomography is reported. Intrinsic phase differences of the reflections used cancel from the three-dimensional phase image of the specimen. With simultaneous hkl and [hk\bar{l}] reflections of a synchrotron radiation beam in a side-by-side geometry, the size of the usable field of view is doubled and the investigated specimen volume is increased by a factor of four. As an example, the reconstructed slice of a mouse kidney is shown in phase contrast at 71 keV. Optimized choices of three-beam reflections and matching interferometer geometries useful for applications are presented.

1. Introduction

In recent years, with the increasing use of synchrotron radiation, microtomography (µCT) has become a powerful tool for performing three-dimensional studies of structures on the mesoscopic scale with spatial resolutions ranging from 0.5 µm to 20 µm. Results have been summarized, for instance, in reviews by Graeff & Engelke (1991[Graeff, W. & Engelke, K. (1991). Handbook of Synchrotron Radiation, edited by S. Ebashi, M. Koch & E. Rubinstein, pp. 361-405. Amsterdam: Elsevier Science.]), Kinney & Nichols (1992[Kinney, J. H. & Nichols, M. C. (1992). Annu. Rev. Mater. Sci. 22, 121-152.]), Graeff (1996[Graeff, W. (1996). Proceedings of the International School of Physics Enrico Fermi, Course CXXVII, edited by. E. Burattini & A. Balerna, pp. 155-170. Amsterdam: IOS Press.]) and Bonse & Busch (1996[Bonse, U. & Busch, F. (1996). Prog. Biophys. Molec. Biol. 65, 133-169.]). In most cases it has been sufficient to simply employ X-ray attenuation for achieving the necessary image contrast. However, attenuation contrast is very poor with specimens consisting predominantly of light elements like biological tissues or resin-based materials. To reveal structures in light specimens, X-ray phase-contrast microtomography (PµCT), generated by placing the specimen in one of the interfering beams of an X-ray interferometer (Bonse & Hart, 1965a[Bonse, U. & Hart, M. (1965a). Appl. Phys. Lett. 6, 155-156.]), is an excellent alternative (Momose et al., 1995[Momose, A., Takeda, T. & Y. Itai, Y. (1995). Rev. Sci. Instrum. 66(2), 1434-1436.]; Beckmann et al., 1995[Beckmann, F., Bonse, U., Busch, F., Günnewig, O. & Biermann, T. (1995). HASYLAB Annual Rep. II, pp. 691-692. HASYLAB, Hamburg, Germany.], 1997[Beckmann, F., Bonse, U., Busch, F., Günnewig, O. & Biermann, T. (1997). J. Comput. Assist. Tomogr. 21, 539-553.]) to attenuation contrast.

State-of-the-art PµCT using synchrotron radiation is based on generating parallel projections of the cylindrical specimen by employing a low-divergence synchrotron radiation beam of sufficient width X0 to illuminate the entire diameter of the specimen. Since there are increasing requests to study specimens with diameters exceeding the beam width of most synchrotron radiation sources, some means of beam widening are required. It is the purpose of this investigation to demonstrate that, by establishing three-beam diffraction in the phase-sensing interferometer, widths X0 given by the synchrotron radiation source can be doubled. The principle of this method and its actual use in an experiment performed at the synchrotron radiation source DORIS at DESY, Hamburg, will be described.

2. Principle of multiple-beam interferometry

The general set-up of the PµCT camera is shown in Fig. 1[link]. Synchrotron radiation with energies E in the range 60–200 keV generated by wiggler BW5 is made monochromatic to a relative bandwidth of ΔE/E ≃ 10−3 by Bragg reflection in transmission (Laue) geometry from a single-crystal monochromator LM. Fig. 2[link] shows, in a side and a top view, the skew-symmetric X-ray interferometer (Bonse & Hart, 1965b[Bonse, U. & Hart, M. (1965b). Appl. Phys. Lett. 7, 99-100.]) featuring a long section of parallel interfering beams between mirrors M1 and M2. Transmitted beams at mirrors M1 and M2 have been omitted in the drawing. The upper beam (Fig. 2[link], side view) propagates through the sample, Sa, which normally is a cylinder with axis vertical and located inside a fluid-containing measuring cell, LC (Fig. 1[link]). Wafers M2 and A serve to recombine both beams to form the interference pattern which carries phase-contrast encoded information about the structure of Sa. The pattern is recorded by an area detector consisting of the fluorescent screen FS (Fig. 1[link]), the optical lens, Le, and the CCD detector, CC. The phase shifter, PS, is used to determine the projected phase shift modulo 2π. Details of the phase extraction and the elimination of the 2π ambiguity from the interference pattern have been described by Beckmann et al. (1997[Beckmann, F., Bonse, U., Busch, F., Günnewig, O. & Biermann, T. (1997). J. Comput. Assist. Tomogr. 21, 539-553.]).

[Figure 1]
Figure 1
Phase-contrast microtomography (PµCT) using ≥70 keV synchrotron radiation. e+: synchrotron radiation source; LM: Laue monochromator; Co: collimators; In: interferometer manufactured of a perfect silicon crystal; Sa: sample; SM: sample manipulator; LC: liquid cell; PS: phase shifter; FS fluorescent screen; CC: CCD detector; Le: lens imaging FS to CC; LG: lead glass.
[Figure 2]
Figure 2
X-ray interferometer using simultaneous three-beam case diffraction hkl, [hk\bar{l}]. S, M1, M2, A: Bragg diffracting crystal wafers for coherent splitting and superposition of interfering beams. Sa: sample; Δz(W): separation of beams in vertical (horizontal) direction. Note doubling of incident beam width X0 when W = X0.

Equation (18)[link] below gives the separation Δz of the interfering beams in the vertical plane (Fig. 2[link], side view). As may be seen, in order to accommodate at high energies a sufficient portion of a longer sample within only the upper beam, higher-order reflections have to be used and the distance y1 (Fig. 2[link]) must be increased. To give an example, with 70 keV X-rays a good choice of high-order reflection would be 660 with Bragg angle ΘB = 7.953°. Making y1 = 27.8 mm, sufficient vertical beam separation of 7.69 mm is achieved.

Three-beam diffraction by a wavevector triplet of type {Khkl, K[_{hk\bar{l}}], K000} can be employed for widening the field of view in the horizontal direction. For the upper (lower) beam in Fig. 2[link] the widening occurs through reflection by wafer S (M2). In order to avoid at M1 and A the direct generation of (hkl) beams by ([hk\bar{l}]) beams (and vice versa), we exclude even values of h, k, l which guarantees the difference reflection (0 0 2l) to be forbidden in the diamond lattice. The geometrical situation in reciprocal space is illustrated in Fig. 3[link]. The horizontal plane in real space is given by the incident wave K000 together with the diffractometer axis [001], and a vertical plane is indicated as a rectangular net. K000 is orthogonal to [001] and has been drawn from the center of the Ewald sphere, M, to the origin of the reciprocal lattice. Some reflection pairs of h k ±l type, i.e. {551, [55\bar1]}, {553, [55\bar3]} and {733, [73\bar3]} are shown. Starting from the 660 point, only small changes of direction and/or energy of the incident K000 will be sufficient to excite a three-beam diffraction with one or the other of the pairs. As an example, the wave vectors of the three-beam diffraction {K553, K[_{55\bar{3}}], K000} are explicitly drawn in Fig. 3[link]. Beams reflected as pairs separate from each other in the (horizontal) [001] direction (Fig. 2[link], top view). Equation (19)[link] below gives the separation W acquired when the beam hits the next crystal wafer.

[Figure 3]
Figure 3
Three-beam diffraction hkl, [hk\bar{l}] explained in reciprocal space. M: center of Ewald sphere. For 553, [55\bar3] the waves are shown. Further possible cases, e.g. 551, [55\bar1] and 733, [73\bar3], are indicated.

3. Pilot experiment

To demonstrate the feasibility of three-beam diffraction in PµCT, we investigated a mouse kidney specimen in a pilot experiment using the X-ray interferometer shown in Fig. 2[link] working with the {K551, K[_{55\bar1}], K000} triplet and an X-ray energy of 71 keV. We made y1 = 27.8 mm, the distance between M1 and M2 be 34.5 mm, and the thickness of crystals S, M1, M2 and A be 0.7 mm. All wafers were 15 mm high and about 10 mm wide.

Fig. 4[link] shows projections with the kidney specimen inserted in the beam (left-hand side) and three projections without the specimen (right-hand side). Water was used as the cell fluid. Overall phase shifts are φ = 0, π/3 and 2π/3. The phase shift induced by the sample is obtained by evaluating the change between projections on the left-hand side and projections on the right-hand side. Input data for tomographic reconstruction were collected by repeating the same set of exposures at 1° incremental rotations of the sample for a total of 180 angular positions. For reconstruction we used the method of filtered back-projections (Kak & Slaney, 1987[Kak, A. C. & Slaney, M. (1987). Principles of Computerized Tomographic Imaging. New York: IEEE Press.]). Fig. 5[link] shows one 16.2 µm-thick slice of the reconstructed specimen cylinder which was 2.0 mm in height. The kidney tissue is the bright structure in the center. Phase-shift modulations inside the specimen are recognizable. The tissue is surrounded by the embedding epoxy (grey). The black area is the cell fluid water. The bar at the bottom left indicates in a quantitative way the phase advance which the X-ray wave gains per µm path length travelled, relative to that when travelling in water. Since phase shift per path length is proportional to electron density and the latter is approximately proportional to mass density, we can re-calibrate the phase-shift scale directly to densities if we input the density value 1.18 g cm−3 of epoxy. From this we find the kidney density to range from 1.182 g cm−3 to 1.224 g cm−3, which is reasonable.

[Figure 4]
Figure 4
Set of projections taken at fixed angular sample position and at three different values of overall phase shift φ, with and without the specimen (mouse kidney) in the beam. E = 71 keV. X0: width of incident beam. The incomplete horizontal separation W of beams (W < X0) did not hamper the tomographic reconstruction (in this pilot experiment) (Fig. 5[link]).
[Figure 5]
Figure 5
One reconstructed 16.2 µm-thick slice of mouse kidney specimen. Center bright structure: tissue. Note phase shift modulations. Grey ring: embedding epoxy; black area: water (cell fluid). Bar: phase advance of wave per µm path relative to water. The phase advance is approximately proportional to mass density. From this the spatial mass density distribution of this slice of kidney is found to range between 1.182 g cm−3 and 1.224 g cm−3. E = 71 keV.

On one of the projections in Fig. 4[link], the widths X0(551), X0([55\bar1]) and locations of the pictures due to the outgoing waves K551, K[_{55\bar{1}}] are indicated. Obviously both beams overlap in the center since their relative sidewise shift W < X0. To make optimum use of the widening effect one should make W = X0 and thereby double the beam width. This is straightforward by choosing y1 larger and/or use a reflection with a larger l value. From the general calculations given in the final section[link] we find that, by employing 553 reflections instead of the used 551, a field width of 10 mm would be obtained. For a particular synchrotron radiation beamline, the value of X0 usually has to stay within narrow bounds. At BW5, X0 is 5–6 mm so that a usable beam of width 10–12 mm can be obtained for PµCT by employing three-beam diffraction. It should be noted, however, that the reconstruction was not hampered by the partial overlap of beams in this pilot experiment.

4. Optimizing the choice of reflections and geometry

In order to survey the favorable choices of geometry and reflections for three-beam PµCT, we include a brief derivation of the underlying formalism. To be specific, we consider the three-beam reflection {Khkl, K[_{hk\bar{l}}], K000} on the skew-symmetric interferometer shown in Fig. 2[link], but with a slightly more general choice of orientation. Our aim is to calculate the vertical distance Δz between interfering beams and the horizontal widening W of beams for a variety of reflections and crystallographic orientations of the interferometer and for different energies. In Fig. 2[link] we see that the unit normal nc on the surfaces of the interferometer wafer crystals is defined by postulating that it has to be normal to the `main' reflection gmain = (h, k, 0)/a0 (generating the lower and upper interfering beams) and to the pair of `sidewise splitting' reflections ±gsplit = (0, 0, ±l)/a0. Referring to crystallographic coordinates, this means that

[{\bf n}_c=(k,\bar{h},0)/C^{1/2},\eqno(1)]

with

[C=h^2+k^2.\eqno(2)]

Again from Fig. 2[link], we see that the unit normal nb parallel to the interferometer's wafers and base is

[{\bf n}_b=(h,k,0)/C^{1/2}.\eqno(3)]

The diffraction vectors are

[{\bf g}^+=(h,k,l)/a_0,\quad\quad\quad{\bf g}^-=(h,k,\bar{l})/a_0,\eqno(4)]

where a0 = 0.543102 nm is the lattice parameter of silicon. The incident wave K000 = (K1, K2, K3) has to be symmetric to the splitting reflections ±gsplit given above which means that K3 = 0. Diffracted waves are determined by Laue equations

[{\bf K}_{hkl}=\left(K_1^+,K_2^+,K_3^+\right)={\bf K}_{000}+{\bf g}^+,\eqno(5)]

[{\bf K}_{hk\bar{l}}=\left(K_1^-,K_2^-,K_3^-\right)={\bf K}_{000}+{\bf g}^-.\eqno(6)]

We combine equations (1)[link]–(6)[link] by elementary vector algebra under conditions of elastic scattering, i.e. |Khkl| = |K[_{hk\bar{l}}]| = |K000| = λ−1, where λ [nm] = 1.23985/E [keV] is the wavelength to be calculated from the used photon energy E. We arrive at equations quadratic in K1, K2, from which we obtain

[K_1=\lambda^{-1}\left[C^{-1}\left(kF\cos\Theta_{\rm B}-hG\sin\Theta_{\rm B}\right)\right],\eqno(7)]

[K_2=\lambda^{-1}\left[C^{-1}\left(-hF\cos\Theta_{\rm B}-kG\sin\Theta_{\rm B}\right)\right],\eqno(8)]

[F=\left(C-l^{\,2}\tan^2\Theta_{\rm B}\right)^{1/2},\eqno(9)]

[G=\left(C+l^{\,2}\right)^{1/2},\eqno(10)]

and ΘB is the Bragg angle. The condition for the quadratic equation to yield real-valued roots K1 and K2 is found to be

[\lambda\le\lambda_{\rm max}=2a_0C^{1/2}/G^2.\eqno(11)]

For a wavelength λ > λmax, no three-beam diffraction is possible. When λ = λmax, the three-beam case is planar as employed for beam splitting in the interferometer by Graeff & Bonse (1977[Graeff, W. & Bonse, U. (1977). Z. Phys. B, 27, 19-32.]). For λ < λmax, we have non-planar three-beam cases which solely are of interest here. It is obvious from Fig. 3[link] that, for the beam widening to function, a wavelength λ considerably smaller than λmax has to be used.

Of importance is also the unit normal n[_\perp] = (n[_{\perp1}], n[_{\perp2}], 0) in the vertical direction for which we find

[n_{\perp1}=C^{-1}\left(hF\cos\Theta_{\rm B}+kG\sin\Theta_{\rm B}\right),\eqno(12)]

[n_{\perp2}=C^{-1}\left(kF\cos\Theta_{\rm B}-hG\sin\Theta_{\rm B}\right).\eqno(13)]

From equations (4)[link]–(10)[link] we find the unit vector s+ = (s1, s2, s3) in the direction of a diffracted wave

[s_1=C^{-1}\left(kF\cos\Theta_{\rm B}-hG\sin\Theta_{\rm B}\right)+\lambda{h}/a_0,\eqno(14)]

[s_2=C^{-1}\left(-hF\cos\Theta_{\rm B}-kG\sin\Theta_{\rm B}\right)+\lambda{k}/a_0,\eqno(15)]

[s_3=\lambda{l}/a_0.\eqno(16)]

Next we consider the ray which travels in the direction of the diffracted wave bridging the gap between wafers S and M1 of width y1,

[{\bf r}_{\rm SM1}={\bf s}^+y_1/\left({\bf s}^+\!\cdot{\bf n}_c\right).\eqno(17)]

Its projection in the vertical direction, rSM1·n[_\perp], is the distance Δz between interfering beams. With equations (1)[link]–(2)[link], (9)[link]–(10)[link] and (12)[link][link][link][link][link]–(17)[link] we obtain

[\Delta{z}=2y_1\sin\Theta_{\rm B}\left(1+l^{\,2}/C\right)^{-1/2}.\eqno(18)]

The amount W of horizontal separation of beams hkl, [hk\bar{l}] is twice the projection of the ray of equation (17)[link] along the horizontal [001]. We find

[W=2y_1\sin\varepsilon\,\left(\cos^2\Theta_{\rm B}-\sin^2\Theta_{\rm B}l^{\,2}/C\right)^{-1/2}.\eqno(19)]

where sin = {lλ/a0}. 2 is the angle between beams Khkl and K[_{hk\bar{l}}].

Δz and W are important criteria for choosing three-beam reflections suitable for use with given sample dimensions, sample absorption properties and X-ray energy E. The latter is often determined mainly by sample absorption. Using the described algorithm, we surveyed a large number of reflections for energies ranging from 10 keV to 105 keV. In Table 1[link] we list some typical and possibly useful choices of reflections at 80 keV (20 keV) for an interferometer with y1 = 31.7 mm (y1 = 27.8 mm).

Table 1
Exemplary choices of three-beam reflections h k ±l yielding optimal vertical beam separation [Δz of equation (18)[link]] and horizontal beam widening [W of equation (19)[link]]

E: photon energy; λ: wavelength, θB: Bragg angle.

E (keV) λ (nm) hkl θB (°) Δz (mm) W (mm)
80 0.0155 533 5.368 5.274 5.464
80 0.0155 553 6.292 6.396 5.743
80 0.0155 733 6.706 6.889 5.477
20 0.0620 131 10.91 10.03 6.489
20 0.0620 331 14.40 13.64 6.578
20 0.0620 511 17.25 16.18 6.672

5. Conclusion and outlook

We have for the first time employed non-planar multiple-beam diffraction in an X-ray interferometer and were able to successfully apply it to X-ray interferometric PµCT. The pilot experiment presented proves that the phase-sensing abilities and the reconstruction process are not affected by multiple diffraction inside the optical components S, M1, M2 and A of the interferometer. We show that PµCT is well feasible up to photon energies of at least 71 keV. Furthermore, we have illustrated that multiple diffraction can very profitably be employed to widen narrow beams delivered by present synchrotron radiation sources and to adapt the beam geometry to growing sample dimensions. Multiple-beam diffraction can also be used in interferometry with thermal neutrons (Bauspiess et al., 1974[Bauspiess, W., Bonse, U., Rauch, H. & Treimer, W. (1974). Z. Phys. 271, 177-182.]). Finally, we have presented a simple algorithm for surveying the properties of three-beam cases of the h k ±l type over a wide range of energies and include exemplary results useful for PµCT using synchrotron radiation.

Acknowledgements

We thank J. R. Schneider and T. Lippmann of DESY for fruitful discussions regarding the use of beamline BW5 of HASYLAB, and T. Biermann for manufacturing the interferometer and for help during the PµCT measurements. Financial support by the Minister for Wissenschaft und Forschung, NRW (Düsseldorf), is gratefully acknowledged.

References

First citationBauspiess, W., Bonse, U., Rauch, H. & Treimer, W. (1974). Z. Phys. 271, 177–182.  CrossRef CAS Google Scholar
First citationBeckmann, F., Bonse, U., Busch, F., Günnewig, O. & Biermann, T. (1995). HASYLAB Annual Rep. II, pp. 691–692. HASYLAB, Hamburg, Germany.  Google Scholar
First citationBeckmann, F., Bonse, U., Busch, F., Günnewig, O. & Biermann, T. (1997). J. Comput. Assist. Tomogr. 21, 539–553.  CrossRef CAS PubMed Web of Science Google Scholar
First citationBonse, U. & Hart, M. (1965a). Appl. Phys. Lett. 6, 155–156.  CrossRef Web of Science Google Scholar
First citationBonse, U. & Hart, M. (1965b). Appl. Phys. Lett. 7, 99–100.  CrossRef Web of Science Google Scholar
First citationBonse, U. & Busch, F. (1996). Prog. Biophys. Molec. Biol. 65, 133–169.  CrossRef CAS Web of Science Google Scholar
First citationGraeff, W. & Engelke, K. (1991). Handbook of Synchrotron Radiation, edited by S. Ebashi, M. Koch & E. Rubinstein, pp. 361–405. Amsterdam: Elsevier Science.  Google Scholar
First citationGraeff, W. (1996). Proceedings of the International School of Physics Enrico Fermi, Course CXXVII, edited by. E. Burattini & A. Balerna, pp. 155–170. Amsterdam: IOS Press.  Google Scholar
First citationGraeff, W. & Bonse, U. (1977). Z. Phys. B, 27, 19–32.  CrossRef CAS Google Scholar
First citationKak, A. C. & Slaney, M. (1987). Principles of Computerized Tomographic Imaging. New York: IEEE Press.  Google Scholar
First citationKinney, J. H. & Nichols, M. C. (1992). Annu. Rev. Mater. Sci. 22, 121–152.  CrossRef CAS Google Scholar
First citationMomose, A., Takeda, T. & Y. Itai, Y. (1995). Rev. Sci. Instrum. 66(2), 1434–1436.  Google Scholar

© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775
Follow J. Synchrotron Rad.
Sign up for e-alerts
Follow J. Synchrotron Rad. on Twitter
Follow us on facebook
Sign up for RSS feeds