Received 6 June 2012
Bent perfect crystals as X-ray focusing polychromators in symmetric Laue geometry
The focusing properties of cylindrically bent crystals in symmetric Laue geometry are discussed using the formalism of Fresnel diffraction and the analytical solution of the Takagi-Taupin equations for a point source on the entrance surface. The existence of a focal shift in the dynamical focusing effect is pointed out and discussed. The present theoretical framework is applied to experiments performed at the energy-dispersive X-ray absorption spectroscopy beamline of the European Synchrotron Radiation Facility concerning the position and the size of the focal spot obtained from a polychromatic source at a large distance from the bent crystal.
Bent crystal plates in Laue geometry are frequently used at synchrotron beamlines as focusing elements, forming a demagnified image of a distant source. The present theoretical framework was developed for the symmetric Laue case with the purpose of comparing theoretical predictions with experimental results obtained by some of the authors (D. Bhattacharyya, O. Mathon and S. Pascarelli).
Rays of different energies coming out of a point source and in exact Bragg incidence on a flat crystal in transmission geometry converge downstream of it in the point symmetric to the source with respect to the crystal mid-plane. In the case of a cylindrically bent crystal, this `polychromatic focusing' effect, which does not depend to a first approximation on the crystal thickness, occurs under the condition (Chukhovskii & Krisch, 1992)
where p is the source-to-crystal distance, q is the crystal-to-focus distance, is the Bragg angle and R is the radius of curvature of the crystal, which is taken as positive in the present investigation, that is to say the source faces the convex side of the crystal. In the design of an energy-dispersive X-ray absorption spectroscopy beamline, one uses a position-sensitive detector downstream of a sample placed in the focal position. In this way, the absorption spectrum of the sample can be recorded over the whole energy range of the reflected beam (Hagelstein et al., 1995).
The optical performance of such a crystal (polychromator), in particular the size of the focal spot, depends on the crystal thickness via the `dynamical focusing' effect, which was described decades ago by Afanasev & Kohn (1977) and Aristov et al. (1978, 1980) in the flat crystal case, and by Kushnir & Suvorov (1982) in the bent crystal case; see also the book by Pinsker (1978). According to Mocella et al. (2004, 2008), the focal size should be minimized when the dynamical and polychromatic foci are made to coincide by the choice of the crystal thickness (see also Sanchez del Rio et al., 1994).
A novel aspect of the theory of dynamical focusing is reported in the present study and regards a theoretical framework explaining a non-negligible displacement of the focus position (focal shift) as well as the influence of anomalous absorption; this is briefly presented in the flat crystal case for the sake of clarity in §2. Our main contribution is, however, to show that the wavefunction downwards from a bent crystal is not obtained by successive convolutions, owing to the analytical form of the phase factor by which the propagator of a bent crystal differs from that of the flat crystal case. This finding has lead us not to use the Fourier transform formalism which is convenient for flat crystals, as shown for instance by Kohn et al. (2000).
The calculations are instead carried out following the geometric scheme depicted in Fig. 1. Let us consider a monochromatic point source S and let O be the point of the crystal entrance surface such that SO is a ray in exact Bragg incidence, and let be the point at the crystal exit surface of this ray. With being the wavelength and SO = p, the incident amplitude along the axis O perpendicular to SO is
which becomes a delta function as (point source on the crystal surface). The crystal propagator is expressed in such a form as to obtain the Bragg-reflected amplitude D(u) along the axis u perpendicular to the Bragg direction q; the reflected amplitude at any point downstream from the crystal can be then calculated by a Fresnel diffraction integral.
| || Figure 1 |
Illustration of the optical setup under consideration, showing the different systems of coordinates used in the text. The O axis is perpendicular to Oso; the u axis is perpendicular to q and to Osh. The shape of the caustic curve is schematically shown on the right-hand side.
Let us initially consider the case of a point source on the entrance surface of a flat crystal. The reflected amplitude D(u) along the axis u in vacuum is obtained by setting in the function ; ; so,h are oblique coordinates along the incident and Bragg directions, respectively, with a common origin in O; are the Fourier coefficients of the crystal susceptibility; t is the crystal thickness; and is the half-width of the reflected beam (see Fig. 1). The Bessel function argument becomes then Z(a2 - u2)1/2, with , the real part of Z being chosen as positive: . For and , by using the asymptotic expression
and also , one obtains
This phase modulation is responsible for the u-dependent angular deviation:
the modulus of which is larger than . The rays around u = 0 cut q at distances shorter than q0 and are tangential to the caustic curve shown schematically in Fig. 1, its cusp being the paraxial focus at the distance q0; this is indeed a cylindrical aberration effect. The point of maximum intensity along the optical axis q, considered as the `effective focus', is expected to be at an `effective focal distance' , with . Numerical calculations of the intensity distribution along the optical axis , according to the integral expression
allow the determination of the parameter , which is a function of the complex parameter . In the case of a non-absorbing crystal, aZ is real and equal to , where is the Pendellösung distance. We have found numerically that increases from for aZ = 10 to for aZ = 40. The focusing condition can be written as in order to include the virtual focus corresponding to the divergent term of equation (4).
Formula (4) also gives the opportunity to discuss the influence of anomalous absorption (Borrmann effect) in a new and straightforward manner. Clearly, due to the phase term , if the real and imaginary parts of Z are positive, the modulus of the convergent term in equation (4) is decreased, whereas that of the divergent term is increased. This is an unfavourable condition for achieving a good focusing of the reflected wave. A favourable condition would be instead and for the first term in equation (4) or and for the second term in equation (4), namely that the real and the imaginary parts of should have different signs. This condition can be most simply formulated as .
One can write , where and are, respectively, the Fourier coefficients of the real and imaginary parts of the complex susceptibility and are written in the form and . In the case of a centrosymmetric crystal, the cosine term is equal to and is always equal to 1 for a mono-atomic crystal. A counterexample would be the 222 reflection of calcite, CaCO3 (Authier, 2005). It would be interesting to investigate experimentally the dynamical focusing effect in such a case, i.e. under conditions of strong anomalous absorption.
Let us use two mutually orthogonal axes (O, Oz) tangential and normal, respectively, to the crystal surface. The diffraction vector of modulus is antiparallel to O, as shown in Fig. 1. Upon bending, the length of the neutral surface of radius R is not changed, whereas the lengths of the entrance and exit surfaces are multiplied respectively by . The component of the lattice displacement along O is . In terms of the oblique coordinates we obtain
Using the notation r for the position vector of an arbitrary point in the crystal, the total wavefield can be expressed as , being the vacuum wavevector parallel to Oso. The amplitudes Do,h(so,sh) are solutions of the Takagi-Taupin equations (TTEs) (Takagi, 1962, 1969; Taupin, 1964, 1967)
the incident amplitude transforming into Einc(so,sh) = .
In the case of a point source in O on the crystal entrance surface, Dinc(so,sh) and Einc(so,sh) are both equal to the delta function and Eh(so,sh) reduces to the well known expression for the flat crystal case (Kato, 1961; Pinsker, 1978; Authier, 2005). This means that
In the more general case of a distant source, corresponding to a distribution of coherent elementary point sources on the entrance surface, we need the analytical form of Dh(so,sh) for a point source on the entrance surface such that its oblique coordinates () are not equal to 0. We then set and consequently Einc(so,sh) = . We thus obtain
where the notation is used. Relations (8) and (9) are used in the next section to obtain the crystal propagator linking the reflected amplitude D(u) along the u axis with the incident amplitude along the O axis; we shall omit the factor , which is a constant for a crystal plate of uniform thickness.
Consider the points N and Q on the same line parallel to Osh, N being on the axis Ou, Q on the curved exit surface (see Fig. 1). These points have the same coordinate so = , where u denotes the position of N on the axis u. It is convenient to write the product [sosh]Q as a polynomial function of degree 2 in u. The point Q being the intersection of the crystal boundary curve z = by the straight line = a + u, we can calculate the rectangular coordinates
and consequently the oblique coordinates
The last term in the numerator can be neglected since . This approximation indicates that the effects of the curvature of the exit surface are negligibly small in the calculation of [sosh]Q.
The linear phase term accounts for the uniform angular deviation of the reflected beam. The focusing condition
is derived by the same approach used in §2, by taking into account the additional phase term . It should be noted that does not depend on the crystal curvature. A real focus is obtained if the value of q as given by equation (11) is positive. Since q0 [see equation (5)] and are both positive, we obtain:
(a) If , there is one real focus at a distance .
(b) If , there are two real foci at distances .
For a point M on the O axis and the point P on the curved entrance surface, with MP parallel to SO, one can easily find , where denotes the position of the point M on the axis O (see Fig. 1). Using equation (6), we obtain - u(u + a)]. Bearing in mind equation (2) for the incident amplitude, the reflected amplitude along the u axis is the integral expression
where the product of the two last factors in the integrand is the crystal propagator. We point out that this is not a function of the difference only; in other words, this propagator is not space-invariant; the propagation through the bent crystal is not a convolution, in contrast to the propagation through a flat crystal. This can be understood physically by regarding the fact that if the incident monochromatic wave is a plane one, the Bragg-reflected wave is not a plane wave, because the departure from the Bragg condition varies continuously along the entrance surface of the bent crystal.
Using the integration variable and introducing the definition , we obtain
We also define qe to be such that and note that pe - 1 + qe - 1 = p - 1 + q - 1. The amplitude distribution at the distance q is obtained by means of the Fresnel diffraction integral
having denoted L = p + q; we also define Le = pe + qe and we obtain the general formula
and the simple focusing condition
This novel (to the best of our knowledge) formula is a clear generalization of the focusing condition in the case of a flat crystal, which under the same conditions would be . The focal distances are thereby
Since a real focus corresponds to a positive value of q1 or of q2, we obtain:
(a) no real focus if (note that );
(b) one real focus if has a value between pe and ; and
(c) two real foci if .
The intensity profile is symmetric with respect to x = - a(p - q)/2p. The focusing condition is . There is one real focus at the distance
where the last expression is valid for , a condition generally satisfied by synchrotron experimental setups. Note that the crystal curvature is such that .
In the approach of Mocella et al. (2004, 2008), the illumination on the crystal entrance surface is considered to be incoherent, as a consequence of polychromaticity. There would therefore be no influence of the source-to-crystal distance on the dynamical focusing behaviour. Propagation in the bent crystal is regarded to be the same as in a flat crystal, as if the displacement vector were perpendicular to the diffraction vector, which is true only along the neutral surface. The effect of bending is a phase factor related to the curvature of the exit surface. Using the dynamical focusing condition , which is the same as in §4.1 with the sign `+' and with , we obtain the focusing condition 2pq(p + q) - 1 = q0 for coincidence of the dynamical and polychromatic focusing. The focus position is then obtained as
which is in strong disagreement with formula (17) of the present formulation. Actually, the assumption of effective incoherence, based on considering the crystal propagator to be nearly independent of , is justified for the calculation of the reflected intensity distribution close to the crystal exit surface. This is, however, not the case for the in-vacuum propagation after the crystal.
We have applied the present theoretical framework and in particular equation (16) to analyse in a computational manner a series of experiments performed at the energy-dispersive EXAFS beamline ID24 of the ESRF, using the Si 111 reflection of symmetrically cut polychromators in Laue geometry with crystal thicknesses of 200 or 300 µm and mean photon energies 8.3 and 20.0 keV. The horizontally diffracting crystal was placed in an adjustable bender at a fixed distance p = 29.7 m from a secondary source produced by a horizontally focusing mirror acting on the primary beam. A detector was placed at a variable distance q downstream from the crystal, the crystal curvature being adjusted experimentally to minimize the beam size in order to get the focal spot matching the polychromatic focusing conditions for each value of q. No auxiliary measurement of the crystal curvature was made. The experiment's purpose was to detect a minimum of the adjusted beam size for a value of q corresponding to the polychromatic and dynamical focusing coincidence: following the approach by Mocella et al. (see §4.3), this was expected to occur for , as predicted by equation (18). The FWHMs shown in Figs. 2 and 3 were obtained via a Gaussian-Lorentzian curve fitting of the signal recorded while scanning the focal spot with a slit of approximately 10 µm aperture. In order to obtain the actual FWHM of the beam one should deconvolve the recorded signal taking into account the intrinsic width of the slit. The initial goal of the experiment could not be met since, for each combination of the photon energy and the crystal thickness, a continuous decrease of the measured beam size was detected for q decreasing from about 1500 to 600 or 700 mm (see Figs. 2 and 3). Further reduction of q was not possible due to the mechanical constraints of the beamline setup. In the following, we interpret the experimental results within the present theoretical framework.
| || Figure 2 |
Beam-size measurements for bent Si 111 crystals of two different thicknesses and a polychromatic beam of mean energy 20.0 keV, p = 29 700 mm.
| || Figure 3 |
Beam-size measurements for a bent Si 111 crystal and a polychromatic beam of mean energy 8.3 keV, p = 29 700 mm.
Let us first consider the experiment at 20 keV. For the 200 µm thick crystal (Fig. 2), the experimental findings show that the minimum spot size is found at a distance of less than 600 mm, i.e. a trend which appears to be in qualitative agreement with the calculations predicting a distance of 500 mm (cf. Fig. 4 and Table 1). For the 300 µm crystal, the calculations (Fig. 5) predict a distance of 820 mm, which is inside the experimental range, so that the minimum should have been observed. The failure to detect the minimum can be explained by an additional random deformation of the perfect bent crystal due to internal stresses induced by the bending action and by the heat load; this may result in random fluctuations in the wavefront converging from the crystal to the focus and consequently cause a focus broadening which would increase with the distance q. This is likely to produce a shift of the minimum spot size towards smaller q values; it can also explain the fact that the measured minimum spot sizes are about twice as large as those predicted by the calculations shown in Figs. 4 and 5. It should be mentioned that the energy-dependent numerical values of the crystal susceptibility used in the calculations were taken from the database of optical constants included in the software package XOP (Sanchez del Rio, 2011).
| || Figure 4 |
Calculated intensity curves for a bent Si 111 crystal, energy = 20 keV, crystal thickness = 200 µm, p = 29 700 mm, , 1.684×10 - 12 + i 1.666 ×10- 14, q0 = 2997 mm. (a) Intensity at the symmetry centre versus q showing the intensity maximum at q = 500 mm, (b) intensity profiles for q = 500 mm (solid curve) and for q = q0/4 = 750 mm (dashed curve). Note that = 0.67.
| || Figure 5 |
Calculated intensity curves for a bent Si 111 crystal, energy = 20 keV, crystal thickness = 300 µm, q0 = 4500 mm, the remaining parameters are the same as for Fig. 4. (a) Intensity at the symmetry centre versus q showing the intensity maximum at q = 820 mm, (b) intensity profiles for q = 820 mm (solid curve) and for q = q0/4 = 1125 mm (dashed curve). Note that = 0.73.
There is apparently no significant spot broadening effect associated with the polychromatic focusing process: the width of the beam hitting the crystal is several millimetres, so the crystal acts as a polychromator diffracting an energy band as large as about 1000 eV; when a narrow slit placed close to the crystal in order to select a limited energy band is moved across the incident beam, the recorded focal spots do not show any significant lateral displacement; this shows that the bending device allows one to approach the following ideal conditions for polychromatic focusing in symmetric Laue geometry: a focusing bent crystal in the form of a hyperbola with the source at one of its foci (Hrdy, 1990).
For an energy of 8.3 keV, the calculated intensities at the symmetry centre show a peculiar behaviour as a function of q (Fig. 6) when compared to the calculations for the 20.0 keV measurements (Figs. 4a and 5a). The related lateral beam profiles are rapidly oscillating and also differ significantly from the corresponding profiles at 20 keV. This is related to anomalous absorption, as discussed in §2. The interference of the convergent and divergent wave components is a minor effect at 20 keV, because the spreading out of the divergent component causes its amplitude to be much smaller than that of the convergent component in the focal region; at 8.3 keV, the amplitude of the divergent component is enhanced with respect to that of the convergent one, as a result of anomalous absorption, and their interference effect becomes significant, in spite of the spreading out of the divergent component.
| || Figure 6 |
Calculated intensity at the symmetry centre versus q for a bent Si 111 crystal, energy 8.3 keV, p = 29 700 mm, , 5.781×10 - 11 + i 3.260×10 - 12. (a) Crystal thickness = 200 µm, (b) crystal thickness = 300 µm.
The present approach to focusing via a bent crystal in symmetric Laue geometry is mainly analytical, in contrast e.g. to the approach of Nesterets & Wilkins (2008), which is based on a numerical solution of the TTEs.
Equation (4), which results directly from the spherical-wave theory (Kato, 1961), allows one to explain in a very simple way the dynamical focusing effect and the influence of anomalous absorption on it. It is shown that the bent crystal propagator is not space-invariant and a new formulation of the dynamical focusing is given in terms of only two quantities, which are the effective focal distance of the flat crystal case and , which is the product of the radius of curvature and the cosine of the Bragg angle. This novel formulation is applied to the case in which dynamical and polychromatic focusing coincide and provides a comprehensive interpretation of the experimental results, ensuring a far better agreement compared to former computational approaches.
In a future study we intend to propose a generalization of the present theoretical approach to the case of asymmetric crystal reflections, involving the use of confluent hypergeometric functions as the solutions of the TTEs (Petrashen', 1974; Chukhovskii & Petrashen', 1977).
We are greatly indebted to the referees for their comments. Following their recommendations, we have stressed which results are novel compared to the existing literature on this subject.
Afanasev, A. M. & Kohn, V. G. (1977). Fiz. Tverd. Tela, 19, 1775.
Aristov, V. V., Polovinkina, V. I., Afanas'ev, A. M. & Kohn, V. G. (1980). Acta Cryst. A36, 1002-1013.
Aristov, V. V., Polovinkina, V. I., Shmytko, I. M. & Shulakov, E. V. (1978). JETP Lett. 28, 4-7.
Authier, A. (2005). Dynamical Theory of X-ray Diffraction, 3rd ed. IUCr Monographs on Crystallography. Oxford University Press.
Chukhovskii, F. N. & Krisch, M. (1992). J. Appl. Cryst. 25, 211-213.
Chukhovskii, F. N. & Petrashen', P. V. (1977). Acta Cryst. A33, 311-319.
Goodman, J. W. (1968). Introduction to Fourier Optics. San Francisco: McGraw-Hill.
Hagelstein, M., Ferrero, C., Hatje, U., Ressler, T. & Metz, W. (1995). J. Synchrotron Rad. 2, 174-180.
Hrdy, J. (1990). Czech. J. Phys. 40, 1086-1089.
Kato, N. (1961). Acta Cryst. 14, 526-532.
Kohn, V. G., Snigireva, I. & Snigirev, A. (2000). Phys. Status Solidi B, 222, 407-423.
Kushnir, V. I. & Suvorov, E. V. (1982). Phys. Status Solidi A, 69, 483-490.
Mocella, V., Ferrero, C., Hrdý, J., Wright, J., Pascarelli, S. & Hoszowska, J. (2008). J. Appl. Cryst. 41, 695-700.
Mocella, V., Guigay, J. P., Hrdý, J., Ferrero, C. & Hoszowska, J. (2004). J. Appl. Cryst. 37, 941-946.
Nesterets, Y. I. & Wilkins, S. W. (2008). J. Appl. Cryst. 41, 237-248.
Petrashen', P. V. (1974). Sov. Phys. Solid State, 15, 2096-2097.
Pinsker, Z. G. (1978). Dynamical Scattering of X-rays in Crystals. Springer Series in Solid-State Sciences. Berlin: Springer Verlag.
Sanchez del Rio, M. (2011). Proc. SPIE, 8141, 814115.
Sanchez del Rio, M., Ferrero, C., Chen, G.-J. & Cerrina, F. (1994). Nucl. Instrum. Methods Phys. Res. A, 347, 338-343.
Takagi, S. (1962). Acta Cryst. 15, 1311-1312.
Takagi, S. (1969). J. Phys. Soc. Jpn, 26, 1239-1253.
Taupin, D. (1964). Bull. Soc. Fr. Mineral. Cristallogr. 87, 69.
Taupin, D. (1967). Acta Cryst. 23, 25-35.