International
Tables for
Crystallography
Volume C
Mathematical, physical and chemical tables
Edited by E. Prince

International Tables for Crystallography (2006). Vol. C, ch. 5.3, pp. 516-534

Section 5.3.3. Methods with counter recording

E. Gałdeckaa

aInstitute of Low Temperature and Structure Research, Polish Academy of Sciences, PO Box 937, 50-950 Wrocław 2, Poland

5.3.3. Methods with counter recording

| top | pdf |

5.3.3.1. Introduction

| top | pdf |

Although, theoretically, the limit of accuracy in all methods based on the Bragg law [equation (5.3.1.1)[link]] is given by the accuracy of the wavelength measurement [(\delta\lambda/\lambda\sim10^{-6})], with photographic recording this limit is not attained. Surprisingly high accuracy may be offered by accurately applied Kossel or divergent-beam techniques. In practice, however, even in this case the accuracy achieved is poorer by an order of magnitude.

The use of Geiger–Müller, proportional, or scintillation counters together with a step-scanning motor makes it possible to record the diffraction profile in a quantitative numerical form convenient for data processing, to locate it with better accuracy and precision and, as a consequence, to obtain better accuracy and precision for the Bragg angle and thus for the lattice parameter. To make the most of this possibility, theoretical papers concerning methods of peak location, estimation of systematic and statistical errors, and optimization of the measurement were developed in parallel with constructional and experimental methods.

Methods of lattice-parameter determination with counter recording form a large and heterogeneous group. As well as measurements on two- or four-circle standard diffractometers, a separate method developed by Bond (1960[link]) and a variety of non-dispersive (X-ray and optical interferometry) and pseudo-non-dispersive methods (two- and three-circle spectrometers, multiple-beam techniques, and combined methods) are included in this group.

5.3.3.2. Standard diffractometers

| top | pdf |

The determination of lattice parameters by the use of a standard diffractometer is based, as in the case of photographic methods, on (5.3.1.1)[link] and (5.3.1.2)[link], and the main task is to measure a sufficient number of reflections (the [\theta] values for various hkl indices) for determining and solving the equations and for calculating the unknown parameters. The reflections can be chosen arbitrarily or in a special way (high [\theta] angle, axial or non-axial reflections).

The characteristic feature of measurements performed on a diffractometer is, however, that to satisfy the Ewald condition for a given reflection the crystal and the detector are rotated or, depending on the geometry (equatorial or inclination), shifted round their axes as well. Basic and more detailed information about the geometry of diffractometers is given elsewhere (Arndt & Willis, 1966[link], Chap. 3; Stout & Jensen, 1968[link], Section 6.3; Kheiker, 1973[link], Chap. 4; Luger, 1980[link], Chap. 4; Section 2.2.6[link] of this volume). For calculating the setting angles for given hkl reflections, the lattice parameters (at least preliminary values) have to be known, and conversely, if the setting angles are known, it is possible to calculate or to refine lattice parameters. Therefore, not only the [\theta] values (given by the angle [2\theta] of rotation of the detector about the goniometer axis) but also the values of the remaining setting angles (i.e. [\omega], [\varphi], and χ of the crystal rotation in equatorial geometry, or μ and [\varphi] for the crystal and v for the detector in inclination geometry) can be used for lattice-parameter determination. This problem can be treated by a matrix analysis.

5.3.3.2.1. Four-circle diffractometer

| top | pdf |

In the case of an automated four-circle (equatorial geometry) diffractometer, the setting angles are calculated by means of the orientation matrix U, i.e. a matrix such that [{\bi A}^*={\bi UA}_G, \eqno (5.3.3.1)]where [{\bi A}^*= \left[\matrix{ a^* \cr b^* \cr c^*}\right] \eqno (5.3.3.1a)]is the reciprocal-axis system with metric [{\bi G}^{-1}=\left[\matrix{a^{*2}&a^*b^*\cos\gamma^*&a^*c^*\cos\beta^* \cr a^*b^*\cos\gamma^*&b^{*2}&b^*c^*\cos\alpha^* \cr a^*c^*\cos\beta^*&b^*c^*\cos\alpha^*&c^{*2}}\right] \eqno (5.3.3.1b)]and [{\bi A}_G=\left[\matrix{a_G \cr b_G \cr c_G}\right] \eqno (5.3.3.1c)]is the crystal-fixed orthonormal system. As can be proved (Busing & Levy, 1967[link]; Hamilton, 1974[link]; Luger, 1980[link], Section 4.1.1; Gabe, 1980[link]), the reciprocal-cell parameters are related to the orientation matrix by the following equation: [{\bi A^*A^*{'}=U\cdot U'}, \eqno (5.3.3.2)]where [{\bi A^*A^*{'}=G}^{-1}] is given by (5.3.3.1b[link]). It is thus possible to calculate the lattice parameters from the terms of the orientation matrix.

The determination of the orientation matrix is usually the first step in measurements performed on the four-circle diffractometer. This task can be accomplished when the preliminary lattice-parameter values are known, and even when they are unknown. In the first case, the setting angles of two reflections, and, in the second, of three reflections, have to be determined. The procedure (Busing & Levy, 1967[link]; Hamilton, 1974[link]) is usually accomplished by the software of the four-circle diffractometer. Least-squares refinement of the lattice and orientation parameters may be performed when the setting angles of several reflections have been observed (Clegg, 1984[link]). Appropriate constraints, resulting from the presence of symmetry elements in the given crystal structure, to be introduced during the refinement, are discussed by Bolotina (1989[link]).

In a particular case, the four-circle diffractometer can be used for lattice-parameter measurements performed in the plane perpendicular to the main goniometer axis (say, the horizontal plane), for which χ = 0°, so that, in practice, only [2\theta] and [\omega] values are used for lattice-parameter determination (see also §5.3.3.4.1[link]). The equations to be solved can be simplified if only axial reflections are taken into account. In an example described by Luger (1980[link], Section 4.2.2), the [{\bf b}^*] axis of a monoclinic crystal is oriented in the direction of the main axis. Then each of the two axial lengths, [a^*] and [c^*] (see Fig. 5.3.3.1[link] ), can be obtained from only one measurement: [\eqalignno{ a^*&={2\sin\theta \over|h|\lambda}, & (5.3.3.3a) \cr c^*&={2\sin\theta \over|l|\lambda}, &(5.3.3.3b)}%fd5.3.3.3b]whereas [\varphi] values of two reflections are used to determine the [\beta^*] angle between [{\bf a}^*] and [{\bf c}^*] axes, since [\beta^*=\varphi_{h00}-\varphi_{00l}. \eqno (5.3.3.3c)]This method is more suitable for orthogonal systems than for non-orthogonal ones, because of the difficulties in obtaining the proper orientation in the case of the monoclinic and, particularly, the triclinic system. In the latter case, the crystal has to be set three times.

[Figure 5.3.3.1]

Figure 5.3.3.1 | top | pdf |

Determination of reciprocal-lattice angles on the θ circle (after Luger, 1980[link]).

5.3.3.2.2. Two-circle diffractometer

| top | pdf |

Lattice-parameter determination by the use of the two-circle (inclination) diffractometer, the so-called `Weissenberg diffractometer', is more troublesome than by means of the four-circle one, because only two rotations [[\omega] (or [\varphi]) of the crystal, and [2\theta] (or γ) of the detector] are motor-driven under computer control, while two inclination angles (μ for the crystal and ν for the detector) must be set by hand.

The problem of application of the popular two-circle (Eulerian-cradle) diffractometer for measurements similar to those presented in §5.3.3.2.1[link] was discussed by Clegg & Sheldrick (1984[link]). The main idea of their paper was to introduce equations combining setting angles, obtained for selected reflections, with reciprocal-cell parameters, for calculating the latter. The authors started with zero-layer reflections for which, for a crystal mounted about the c axis, [\eqalignno{ \sin\theta& = (x^2+y^2){}^{1/2}, &(5.3.3.4a) \cr \omega &= \omega_0+\theta-\tan^{-1} (\,y,x), &(5.3.3.4b)}%fd5.3.3.4b]where [\eqalignno{ x &= \lambda(ha^*+kb^*\cos\gamma^*)/2, &(5.3.3.4c) \cr y&= (\lambda kb^*\sin\gamma^*)/2, &(5.3.3.4d)}%fd5.3.3.4d]and [\omega_0] is a zero-point correction.

The remaining parameter c had to be determined from the inclination angle μ, measured by hand. The use of zero-layer reflections was advantageous, apart from the simplicity of the formulae (5.3.3.4a[link],b[link],c[link],d[link]), because they were less affected by crystal misalignment than were upper-layer reflections. However, a zero-point correction [\omega_0] for [\omega] had to be performed. For this purpose, the [\omega_0] value was treated as an additional parameter in off-line least-squares refinement.

As the next step, the authors introduced equations for a general crystal orientation instead of an aligned crystal (cf. §5.3.3.2.1[link]) and derived equations defining the setting angles for an arbitrary reflection useful for data collection from a randomly oriented crystal if preliminary lattice-parameter values had been assumed. This made possible measurements of reflections on a range of layers; only one crystal mounting was required. The matrix formulae suitable for Eulerian-geometry diffractometers are also given by Kheiker (1973[link], Chap. 3, Section 9) and Gabe (1980[link]).

In order to perform precise refinement of all six cell parameters, Clegg & Sheldrick (1984[link]) used least squares with empirical weights: [W_{hkl}=1/\sqrt {\omega_{hkl}}, \eqno (5.3.3.5)]where [\omega_{hkl}] is the width of the hkl reflection. An additional (third) motor to control the μ circle was proposed.

The authors point out that the two-circle diffractometer, owing to its simpler construction in comparison with the four-circle one, is well suited to operations that require additional attachments; for example, for low-temperature operation.

5.3.3.3. Data processing and optimization of the experiment

| top | pdf |

5.3.3.3.1. Models of the diffraction profile

| top | pdf |

Every measurement is based on a certain model of its object. By `model' we understand here2 all the systematized a priori knowledge concerning the given measurement, necessary for planning and performing the experiment and for estimating parameters being determined. The use of an incorrect model results in a bias, i.e. an additional systematic error that may appear aside from physical and geometric aberrations. Therefore, the choice of a well founded model is essential in accurate measurements.

In the case of lattice-parameter determination, the object of direct measurements is a diffraction profile, already mentioned in Subsection 5.3.1.1[link], and the quantity that is directly determined from the experiment is the Bragg angle [\theta].

The a priori information about the diffraction profile should define: (i) the way in which the Bragg angle [\theta] is related to the measured profile [h(\omega)], i.e. a measure of location; (ii) the mean values of the measured intensities within the profile; and (iii) their variances.

  • (i) In traditional photographic methods, the Bragg angle is determined from the measurement of distance on the film, where points or lines of the most intense blackening are taken into account. The blackening, which corresponds to the recorded intensity, may be estimated qualitatively (`by eye') or quantitatively, by means of a special device. In the second case, the intensity is determined as a function of the coordinates on the photograph, which, in turn, are related to the angular positions of diffracted beams. The distribution so obtained, i.e. the line profile or the diffraction profile, allows more precise measurements of the distances and the determination of [\theta] angles, if a definition of the point [(\theta_0, h_0)] of the profile [h(\theta)], corresponding to the Bragg angle, i.e. a measure of location, is accepted. The analogous situation appears when the diffraction profile is recorded by means of the counter diffractometer. Then the intensities are measured by a counter, while the angular positions of the detector ([2\theta] scan) or the sample ([\theta] scan), or both ([\omega][2\theta] scan), are controlled by stepping motor. The device is normally combined with a computer, which facilitates the data processing.

    There are various measures of location of the diffraction profile (Wilson, 1965[link]; Thomsen & Yap, 1968[link]). The most popular are:

    • (1) the centroid or the centre of gravity, defined as [\theta_c=\textstyle\int\limits^{\Omega_2}_{\Omega_1} \theta h(\theta)\,{\rm d}\theta\Big/\int\limits^{\Omega_2}_{\Omega_1}\,h(\theta)\,{\rm d}\theta, \eqno (5.3.3.6)]where [\Omega_1] and [\Omega_2] are the selected truncation limits;

    • (2) the median, the value [\theta_m] that equally divides some specified portion of the line profile, i.e. [\textstyle\int\limits^{\theta_m}_{\Omega_1} h(\theta)\,{\rm d}\theta=\int\limits^{\Omega_2}_{\theta_m} h(\theta)\,{\rm d}\theta; \eqno (5.3.3.7)]

    • (3) the geometrical peak – the abscissa value [\theta_p] for which the maximum occurs, i.e. [[{\rm d} h(\theta)/{\rm d}\theta]_{\theta=\theta_p}=0; \eqno (5.3.3.8)]

    • (4) the extrapolated peak or the midchord peak, introduced by Bearden (1933[link]) – the point θep of intersection of two curves, one of them approximating the midpoints of chords drawn through the profile parallel to the abscissa axis (or to the background) and the other approximating the data points (Fig. 5.3.3.2[link] );

      [Figure 5.3.3.2]

      Figure 5.3.3.2 | top | pdf |

      The extrapolated-peak procedure (after Bearden, 1933[link]).

    • (5) the single midpoint of a chord θmc drawn horizontally at the defined height, [\alpha H], where H is the peak height and α is the truncation level, [0\lt\alpha\lt1].

    The advantages and disadvantages of these measures of location have been widely discussed (Wilson, 1965[link], 1967[link], 1968[link], 1969[link]; Thomsen & Yap, 1968[link]; Segmüller, 1970[link]; Kirk & Caulfield, 1977[link]; Grosswig, Jäckel & Kittner, 1986[link]; Gałdecka, 1994[link]), the errors, both systematic (biases) and statistical (variances), resulting from each of these definitions being taken into account. The dependence of these errors on the scanning range (truncation limits) is of great importance. Such features of the definitions as their simplicity or current usage were also considered.

    The geometrical peak of the least-squares parabola, approximating the data points near the top of the profile, distinguishes itself with the best precision but rather large bias (because of the asymmetry of the profiles met in practice); the extrapolated peak – commonly used in the case of the Bond (1960[link]) method (definition 4) – permits location of the peak with better accuracy and omitting the dispersion error (cf. §5.3.3.4.3.2[link]). The centre of gravity, very useful in theoretical considerations (Wilson, 1963[link]), is strongly dependent on the truncation limits and requires a rather large scanning range. The choice of the definition of the measure of location is the first step of lattice-parameter calculations and also of systematic and statistical error estimation.

    In the papers that appeared in the mid-1950's, and which were mainly concerned with powder samples, the centre of gravity as a measure of location was more often used than the peak, probably owing to its property of additivity (the total systematic error in the Bragg angle is a sum of the partial errors related to various physical and apparatus factors) and the estimated errors were consequently referred to this point. The papers were reviewed by Wilson (1963[link], 1980[link]), one of the authors, in the form of a homogeneous mathematical theory of X-ray powder diffractometry. Some of the formulae describing corrections for displacements of the centroid caused by physical and geometrical factors (collected in convenient tables) proved to be useful for single-crystal methods as well (Smakula & Kalnajs, 1955[link]; Kheiker & Zevin, 1963[link]). Wilson (1963[link]) derived the general formula for calculations of the peak displacements due to various factors. As results from this, the displacements are not additive and, in the case when at least one of the partial distributions is asymmetric, the convolution of the curves [see equation (5.3.1.6)[link]] may lead to an appreciable peak shift, if the distributions are not known. The problem has been treated by Berger (1984[link], 1986a[link]), who used computer modelling.

    In later single-crystal methods, in particular in the Bond (1960[link]) method, the peak position of the profile was determined rather than the centroid and the respective corrections referred to the peak (§5.3.3.4.3.2[link]). As a rule, the corrections that related to the peak position were treated as being independent. In practice, this simplifying assumption can be sufficient in measurements with moderate and even high accuracy. However, if the highest accuracy, say of 1 part in 107, is required, the joint effect of all the aberrations should be considered (the so-called `cross terms' are used besides the main terms). Such considerations [Härtwig & Grosswig, 1989[link]; cf. §5.3.3.4.3.2[link], point (7)] must be based on a well-founded physical model of the diffraction profile.

  • (ii) As already mentioned in Subsection 5.3.1.1[link], the diffraction profile can be described as a convolution of several factors (distributions), namely the wavelength distribution, crystal profile and certain aberration profiles. To the so-obtained net profile [equation (5.3.1.6)[link]], a background should be added – constant in the case of an [\omega] scan (as in one-crystal spectrometers, for example), and more complex (but usually approximated with a straight line within a narrow angular range) in other cases. Thus, to describe accurately the distribution of the mean values of measured intensities, all individual distributions must be given.

    Such complete syntheses of the diffraction profile are rarely performed, and only for the highest-accuracy absolute measurements (Härtwig, Hölzer, Förster, Goetz, Wokulska & Wolf, 1994[link]). Since one of the basic factors of the convolution model is the wavelength distribution that characterizes a given source of radiation, its accurate determination and proper scaling in metric units is of primary importance in high-accuracy lattice-parameter measurements. At present, only a few such measurements are reported, which relate to the Cu Kα emission line (Berger, 1986b[link]; Härtwig, Hölzer, Wolf & Förster, 1993[link]; Härtwig, Bąk-Misiuk, Berger, Brühl, Okada, Grosswig, Wokulska & Wolf, 1994[link]) and to the Cu Kβ line (the latter paper). Owing to a relatively simple analytical model proposed by Berger (1986b[link]) to describe the [K\alpha_{1,2}] doublet, the measurement results are easy to handle.

    Profiles connected with individual apparatus factors (collimation, for example) can also be, in principle, described analytically, under some simplifying assumptions. Examples of such profiles are distributions related to the vertical divergence of the beam (Eastabrook, 1952[link]) and to the horizontal (in-plane) divergence (Urbanowicz, 1981a[link]). These are general enough, so can be calculated for given apparatus parameters. While performing high-accuracy measurements, however, the validity of all respective accompanying assumptions must be carefully considered (Urbanowicz, 1981b[link]; Härtwig & Grosswig, 1989[link]; Härtwig et al., 1993[link]).

    In wider practice, there is a tendency towards using simpler descriptions of the diffraction profile. Often, one of the factors, apart from the spectral distribution, is dominant, and the influence of the other ones can be neglected. Berger (1986b[link]), for example, neglecting small effects of both the vertical divergence and the crystal profile, obtained an analytical model of the measured Cu Kα emission spectrum, with several adjusted parameters, and so managed to determine the pure Cu Kα emission-spectrum profile without the necessity of calculating the deconvolution of the measured spectrum in relation to the horizontal-divergence profile.

    The choice of model of the shape of the diffraction profile depends, of course, on the purpose for which it is applied. The simplest possible descriptions are used in low- or medium-accuracy measurements, in which first the measured values of Bragg angles are determined by approximation of the measured profiles with simple analytical functions (polynomials or so-called shape functions), the parameters of which have no physical meaning, and then all necessary corrections are calculated and subtracted from the measured Bragg angles – under the assumption of their additivity, mentioned in (i) – to obtain their true values. Another application of the simple models is just the estimation of systematic and statistical errors of the Bragg-angle determination. The choice and use of such simple models will be shown in §5.3.3.3.2[link].

  • (iii) The knowledge of variances (and covariances) of recorded counts is needed to evaluate the goodness of fit while approximating the measured profile with a given model function (appropriate criteria have been formulated by Gałdecka, 1993a[link],b[link]) and to estimate the precision of the Bragg-angle determination.

    Most often, one assumes that the variances of measured intensities are defined by the Poisson statistic, i.e. [\sigma^2(h)=h, \eqno (5.3.3.9)]where h is the intensity in number of counts.

    Other factors affecting the statistics of recorded counts and the validity of the assumption [equation (5.3.3.9)[link]] have been taken into consideration by Bačkovský (1965[link]) [see also equations (5.3.3.17)[link] and (5.3.3.18)[link] and the comments on these], Wilson (1965[link]), and Gałdecka (1985[link]). The factors are mostly errors in the angle setting and reading and also fluctuations of the primary-beam intensity, of the counting time, and of the temperature of the sample. The use of automatic scanning can cause correlations between intensities measured at different points in the profile (Gałdecka, 1985[link]).

5.3.3.3.2. Precision and accuracy of the Bragg-angle determination; optimization of the experiment

| top | pdf |

The analysis of the variance [\sigma^2(\theta_0)] of a chosen measure of location permits a combination of the precision of the Bragg-angle determination, and so of the lattice-parameter determination [equation (5.3.1.4)[link]], with the scanning range [2\Omega] = [\Omega_2-\Omega_1] [see definition (1), §5.3.3.3.1[link]] or truncation level α [see definition (5)], the number of measuring points n (usually [n=2p+1]), the parameters of the profile (number of counts H in the peak position, the half-width [\omega_h]), and its shape. It is convenient to present the profile [h(\theta)] in a standardized form (Thomsen & Yap, 1968[link]) as: [h(\theta)=Hv[x (\theta)], \eqno (5.3.3.10)]where [x(\theta)=2{\theta-\theta_0\over\omega_h} \eqno (5.3.3.10a)]are standardized angle values and [v(x)=h/H \eqno (5.3.3.10b)]is the shape function, not dependent on the parameters H and [\omega_h]. For each measure of location [definitions (1)–(5) of §5.3.3.3.1(i)[link]], there is the dependence: [\sigma^2(\theta_0)=F{\omega^2_h \over I_pT}, \eqno (5.3.3.11)]where [I_p] is the peak intensity, T is the total counting time, and F is a dimensionless factor that depends on the measure of location and the shape of the profile.

Since, in the case of fixed-time counting, the total counting time T is proportional to the number n of measuring points: [T=n\Delta t, \eqno (5.3.3.12)]where [\Delta t] is the counting time, and since the number of counts h is proportional to the intensity I: [h=I\Delta t, \eqno (5.3.3.13)]and, in particular, the number of counts H in the peak position is proportional to the peak intensity [I_p]: [H=I_p\Delta t, \eqno (5.3.3.13a)]the dependence (5.3.3.11)[link] can be presented as [\sigma^2(\theta_0)= {F\over n} {\omega\,^2_h \over H}. \eqno (5.3.3.14)]Thus, for a given measure of location and given shape of the profile [equations (5.3.3.10)[link], (5.3.3.10b)[link]], the variance [\sigma^2(\theta_0)] depends on the ratio [\omega^2_h/H] of the profile parameters [(\omega_h, H)] and decreases with an increase of the number of points, n.

In particular, the variance [\sigma^2(\theta_p)] of the peak [definition (3), §5.3.3.3.1[link]] of the least-squares parabola has been estimated (Wilson, 1965[link]) as [\sigma^2(\theta_p) = {3H\over2p\Omega{^2}[h''(\theta_p)]{}{^2}}, \eqno (5.3.3.15)]where [h''(\theta_p)] is the second derivative of [h(\theta)] in the peak position and p is a number such that n = 2p + 1 (n [\approx ] 2p, if p is sufficiently large).

Taking into account the standardization performed [equations (5.3.3.10)[link], (5.3.3.10a[link],b[link])], equation (5.3.3.15)[link] can be rewritten in the form: [\sigma^2(\theta_p) = {1\over n}\, {3\over4X{^2}[v''(0)]{}{^2}} \,{\omega^2_h \over H}, \eqno (5.3.3.16)]where X is the standardized scanning range [X=2\Omega/\omega_h \eqno (5.3.3.16a)]and [v''(0)] is the second derivative of the shape function in the peak positions. By comparing (5.3.3.16)[link] and (5.3.3.14)[link], we find the factor F in this case to be [F = {3\over 4X{^2}[v''(0)]{}{^2}}. \eqno (5.3.3.16b)]From (5.3.3.14)[link] and (5.3.3.16b[link]), the variance of the peak of the least-squares parabola decreases with an increase of the scanning range. On the other hand, the bias of the peak position, resulting from the asymmetry of the profile, is proportional to [\Omega^2] (Wilson, 1965[link]): [\Delta\theta_p = 2(\Omega{^2}/\omega_h) [v'''(u_p)/v''(u_p)], \eqno (5.3.3.16c)]where [v''(u_p)] and [v'''(u_p)] are the second and third derivatives of a function describing the profile at its peak position [u_p]. These two aspects should be taken into account in choosing the scanning range. Yet, as shown by Gałdecka (1993b[link]; Section 5), (5.3.3.16c[link]) may be applied to reduce the bias by extrapolating to [\Omega=0] the results obtained within various scanning ranges.

In the case of polynomials of higher (and even) degrees (m = 4, 6, 8) and [0.5\ll X\ll1], the factor F can be expressed by a semi-empirical dependence (Thomsen, 1974[link]; Gałdecka, 1993b[link]): [F=0.0017\,m^2\,(\tan^{-1}X)/X^3, \eqno (5.3.3.16d)]but it is difficult to evaluate the bias. Therefore, as shown by Gałdecka (1993b[link]), polynomials of higher degrees have no advantage over a least-squares parabola.

To minimize the bias, a reasonable shape function may be used rather than a polynomial (Gałdecka, 1993a[link], b[link]). The function should be continuous (including its derivatives), not negative and closely related to known physical models of the diffraction profiles. Since the measured diffraction profiles are, as a rule, asymmetric, the proper selection of a description of asymmetry is of primary importance. The use of the so-called `split functions', consisting of two `half' functions of the same (or different) shape and different half-widths, leads to a noticeable bias, so such functions must not be used for accurate lattice-parameter determination.

The variance [\sigma^2(\theta_{\rm mc})] of a single midpoint of a chord [definition (5), §5.3.3.3.1[link]] has been estimated by Bačkovský (1965[link]) as [\sigma^2(\theta_{\rm mc})=\sigma^2(\theta_i)+\sigma^2(h_i)/[h'(\theta_i)]{}^2, \eqno (5.3.3.17)]where [\sigma^2(\theta_i)] and [\sigma^2(h_i)] are the variances of the coordinates [\theta] and h, respectively, and [h'(\theta)] is the first derivative at the ith point. If it is assumed that [\sigma^2(\theta_i)] is small in relation to the second component of (5.3.3.17)[link] and if (5.3.3.9)[link] and the standardizations (5.3.3.10)[link], (5.3.3.10a[link],b[link]) are taken into consideration, (5.3.3.17)[link] can be rewritten in the form: [\sigma^2(\theta_{\rm mc})={v_i \over 4[v'(x_i)]{^2}}\, {\omega^2_h\over H}, \eqno (5.3.3.18)]where [v'(x_i)] is the first derivative of the shape function in the ith position.

Comparison of (5.3.3.18)[link] and (5.3.3.14)[link], with n = 2, leads to [F= {v_i \over 2[v'(x_i)]{}{^2}}. \eqno (5.3.3.18a)]For an arbitrary shape function [v(x)] describing the diffraction profile, it is thus possible to find such a truncation level [\alpha=\alpha_{\rm opt}] [§5.3.3.3.1[link], definition (5)[link]], for which F is a minimum. If the shape function is the Cauchy function, [v={ 1\over 1+x^2}, \eqno (5.3.3.19)]the optimum truncation level is [\alpha_{\rm opt}=2/3], and the resulting F factor, [F=F_{\rm min}=0.84].

In spite of a large bias introduced by the midpoint of a single chord (the difference between its position and the peak position), this measure of location is preferred by Barns (1972[link]), because the calculations are less time-consuming than those for other points of the profile. Barns takes α = 0.5 [F = 1 for the Cauchy function; equations (5.3.3.18a[link]), (5.3.3.19)[link]] and compensates the bias at this level by determining an effective value of the wavelength based on a silicon standard.

The estimators of the variance for the centroid and the median given by Wilson (1967[link], 1968[link], 1969[link]), or estimators of both the variance and the bias of the extrapolated-peak position given by Gałdecka (1994[link]) can also be the basis of the choice of the scanning range if these measures of location are applied.

The other possibility of affecting the precision of the measurements is to change the shape and the parameters of the profile [see equations (5.3.3.14)[link], and (5.3.3.16b[link]) or (5.3.3.18a[link])] by changing the apparatus parameters [the influence on [h_A(\theta)], equation (5.3.1.6)[link]], or the X-ray source profile [h_\lambda(\theta)], or the crystal profile [h_C(\theta)].

An example of the first possibility is the optimization of the parameters of in-plane collimation in the case when the peak of the least-squares parabola is used as the measure of the location (Urbanowicz, 1981a[link]). Since both the shape and the parameters of the profile depend on the collimation parameters, the task is to choose collimator-slit dimensions to minimize the value [(\omega^2_h/H)\{1/[v''(0)]{^2}\}] [cf. equation (5.3.3.16)[link]]. As a result of detailed considerations, under the assumption given by (5.3.3.9)[link], the optimum exists and is defined by the following formula: [d_1=d_2=0.565L\omega_\lambda, \eqno (5.3.3.20)]where [d_1] and [d_2] are the widths of the slits, L is the collimator length, and [\omega_\lambda] is the half-width of the original profile [h_\lambda(\theta)] [cf. equations (5.3.1.6)[link], (5.3.1.7)[link], and (5.3.1.8)[link]]. Systematic errors connected with collimation have been discussed separately (Urbanowicz, 1981b[link]).

The width of the original profile [h_\lambda(\theta)] can be reduced by means of spectrally narrow sources or by the use of additional crystal(s) in multiple-crystal methods (Subsection 5.3.3.7[link]). The latter also affects the crystal profile [h_C(\theta)].

5.3.3.4. One-crystal spectrometers

| top | pdf |

5.3.3.4.1. General characteristics

| top | pdf |

A diffractometer in which both [2\theta] and [\omega] scans are available, intended for precise and accurate lattice-parameter determination, is sometimes called a one-crystal spectrometer, by analogy with a similar device used for wavelength determination. This name has been used by Lisoivan (1982[link]), who in his review paper described various properties and applications of such a device.

Bragg-angle determination with the one-crystal spectrometer can be performed in an asymmetric as well as in a symmetric arrangement (Arndt & Willis, 1966[link], pp. 262–264). In the asymmetric arrangement (Fig. 5.3.3.3a[link] ), the angle [2\theta] is the difference between two detector positions, related to the maximum intensity of the diffracted and the primary beam, respectively. Bragg-angle determination in such an arrangement is subject to several systematic errors; among these zero error, eccentricity, and absorption are of great importance. As shown by Berger (1984[link]), the latter two errors can be eliminated when Soller slits are used.

[Figure 5.3.3.3]

Figure 5.3.3.3 | top | pdf |

Determination of the Bragg angle by means of the one-crystal spectrometer using (a) an asymmetric or (b) a symmetric arrangement. The zero position of the detector arms must be known in (a), but not in (b). After Arndt & Willis (1966[link]).

To eliminate the zero error, a symmetric diffractometer may be used, in which each measurement of the Bragg angle is performed twice, for two equivalent diffracting positions of the sample, symmetrical in relation to the primary-beam direction (Fig. 5.3.3.3b[link]). The respective positions of the counter (or counters, since sometimes two counters are used) are also symmetrical. Such an arrangement may be considered to be (Beu, 1967[link]), in some ways, the diffractometer counterpart of the Straumanis film method (Straumanis & Ieviņš, 1940[link]). From geometric considerations, the absolute value of the angle between the two counter positions is [4\theta] and the absolute value of the angle between the two sample positions, [\omega_1] and [\omega_2], is 180° − [2\theta], so that both [2\theta] and [\omega] scans can be used for the Bragg-angle determination.

As was mentioned in §5.3.2.3.4(vi)[link], the idea of calculating the [\theta] angle from the two sample positions has been used with photographic methods (Bragg & Bragg, 1915[link]; Weisz, Cochran & Cole, 1948[link]). Bond (1960[link]), in contrast, was the first to apply this to measurements on the counter diffractometer, and proved that, owing to the geometry, not only the zero error but also the eccentricity, absorption, and several other errors can be reduced.

5.3.3.4.2. Development of methods based on an asymmetric arrangement and their applications

| top | pdf |

Although the Bond (1960[link]) method, based on a symmetric arrangement presented in §5.3.3.4.3[link], makes possible higher accuracy than that obtained by means of a standard diffractometer, an asymmetric arrangement proves to be more suitable for certain tasks connected with lattice-parameter measurement, because of its greater simplicity. The more detailed arguments for the use of such a device result from some disadvantages of the Bond method, discussed in §5.3.3.4.3.4[link].

One of the earliest and most often cited methods of lattice-parameter determination by means of the counter single-crystal diffractometer (in an asymmetric arrangement) is that of Smakula & Kalnajs (1955[link]). The authors reported unit-cell determinations of eight cubic crystals. The systematic errors due to seven factors were analysed according to the formulae derived by Wilson (1950[link]) and Eastabrook (1952[link]) for powder samples, and valid also for single crystals. The lattice parameters computed for various diffraction angles were plotted versus [\cos^2\theta]; extrapolation to [2\theta=180]° gave the lattice parameters corrected for systematic errors. Accuracy of 4 parts in 105, limited by the uncertainty of the X-ray wavelength, and precision of 1 part in 106 were achieved.

A more complete list of factors causing broadening and asymmetry of the diffraction profile, and so affecting statistical and systematic errors of lattice-parameter determination, has been given by Kheiker & Zevin (1963[link], Tables IV, IVa, and IVb). Since the systematic errors due to the factors causing asymmetry (specimen transparency, axial divergence, flat specimen) are, as a rule, dependent on the Bragg angle and proportional to [\cos\theta], [\cos^2\theta], [\cot\theta] or [\cot^2\theta], they can be removed or reduced – as in the method of Smakula & Kalnajs (1955[link]) – by means of extrapolation to [\theta=90]°. The problem has also been discussed by Wilson (1963[link], 1980[link]) in the case of powder diffractometry [cf. §5.3.3.3.1(i)[link]]. When comparing the considerations of Kheiker & Zevin and Wilson [the list of references concerning the subject given by Kheiker & Zevin (1963[link]) is, with few exceptions, contained in that given by Wilson (1963[link])], it will be noticed that some differences in the formulae result from differences in the geometry of the measurement rather than from the different nature of the samples (single crystal, powder).

As in the photographic methods, the accurate recording of the angular separation between [K\alpha] and [K\beta] diffraction lines can be the basis for lattice-parameter measurements with a diffractometer (Popović, 1971[link]). The method allows one to reduce the error in the zero setting of the [2\theta] scale and the error due to incorrect positioning of the sample on the diffractometer, since the angular separations are independent of the zero positions of the [2\theta] and [\omega] scales.

An example of a contemporary method of lattice-parameter determination is given by Berger (1984[link]). As has been mentioned in §5.3.3.4.1[link], the characteristic feature of the device is the Soller slits, which limit the divergence of both primary and diffraction beams and, at the same time, eliminate errors due to eccentricity and absorption. On the other hand, systematic errors due to refraction, vertical inclination, vertical divergence, and Soller-slit inaccuracy, as well as asymmetry of profiles and crystal imperfection, have to be analysed.

Since, in this case, the angle between the incident and the reflected beam is measured, the inclinations of both beams must be considered. As a result of the analysis [analogous to that of Burke & Tomkeieff, 1969[link]; referred to in §5.3.3.4.3.2(4)[link]], the following expression for the angular correction [\Delta\theta_t] (to be added to the measured value of [\theta]) is obtained: [\Delta\theta_t={\alpha\gamma \over2\sin2\theta} + {\alpha^2+\gamma^2 \over 4\tan 2\theta}, \eqno (5.3.3.21)]where α and γ are the vertical inclinations of the incident and reflected beams, respectively. The correction for vertical divergence is presented in §5.3.3.4.3.2(3)[link].

The Soller-slit method, the accuracy and precision of which are comparable to those obtained with the Bond method, is suitable both for imperfect crystals, since only a single diffracting position of the sample is required, and for perfect samples, when an exactly defined irradiated area is required. It is applicable to absolute and to relative measurements. Examples are given by Berger, Rosner & Schikora (1989[link]), who worked out a method of absolute lattice-parameter determination of superlattices; by Berger, Lehmann & Schenk (1985[link]), who determined lattice-parameter variations in PbTe single crystals; and by Berger (1993[link]), who examined point defects in II–VI compounds.

An original method, based on determining the Bragg angle from a two-dimensional map of the intensity distribution (around the reciprocal-lattice point) of high-angle reflections as a function of angular positions of both the specimen and the counter, was described by Kobayashi, Yamada & Nakamura (1963[link]) and Kobayashi, Mizutani & Schmidt (1970[link]). A finely collimated X-ray beam, with a half-width less than 3′, was used for this purpose. The accuracy of the counter setting was ±0.1°, the scanning step [\Delta\theta=0.01]°. Systematic errors depending on the depth of penetration and eccentricity of the specimen were reported, and were corrected both experimentally (manifold measurements of the same planes for different diffraction ranges, and rotation of the crystal around its axis by 180°) and by means of extrapolation. The correction for refraction was introduced separately. The method was used in studies of the antiparallel 180° domains in the ferroelectric barium titanate, which were combined with optical studies.

The determination of variations in the cell parameter of GaAs as a function of homogeneity, effects of heat treatments, and surface defects has been presented by Pierron & McNeely (1969[link]). Using a conventional diffractometer, they obtained a precision of 3 parts in 106 and an accuracy better than 2 parts in 106. The systematic errors were removed both by means of suitable corrections (Lorentz–polarization factor and refraction) and by extrapolation.

A study of the thermal expansion of α-LiIO3 over a wide range of temperatures (between 20 and 520 K) in the vicinity of the phase transition has been reported by Abrahams et al. (1983[link]). Lattice-parameter changes were examined by means of a standard diffractometer (CAD-4); absolute values at separate points were measured by the use of a Bond-system diffractometer.

An apparatus for the measurement of uniaxial stress based on a four-circle diffractometer has been presented by d'Amour et al. (1982[link]). The stress, produced by turning a differential screw, can be measured in situ, i.e. without removing the apparatus from the diffractometer. An example of lattice-parameter measurement of Si stressed along [111] is given, in which the stress parameter ζ is calculated from intensity changes of the chosen 600 reflection.

5.3.3.4.3. The Bond method

| top | pdf |

5.3.3.4.3.1. Description of the method

| top | pdf |

By the use of the symmetric arrangement presented in §5.3.3.4.1[link] (Fig. 5.3.3.3b[link]), it is possible to achieve very high accuracy, of about 1 part in 106 (Bond, 1960[link]), and high precision (Baker, George, Bellamy & Causer, 1968[link]) but, to make the most of this, some requirements concerning the device, the sample, the environmental conditions, the measurement itself, and the data processing have to be fulfilled; this problem will be continued below.

Bond (1960[link]) in his notable work used a large, highly pure and perfect single crystal (zone-refined silicon) in the shape of a flat slab. The scheme of the method is given in Fig. 5.3.3.4[link] . The crystal was mounted with the reflecting planes accurately parallel to the axis of the shaft on a graduated circle (clinometer), the angular position of which could be read accurately (to 1′′). The X-ray beam travelling from the tube through a collimator (two 50 µm slits, 215 mm apart, so that the half-width of the primary beam was 0.8′) fell directly upon the crystal, set in one of the two diffracting positions. The diffracted beam was intercepted by one of two detectors [Geiger–Müller (G–M) counters], which were fixed in appropriate positions. The detectors were wide open, so that their apertures were considerably wider than the diffracted beam, which eliminated some systematic errors depending on the counter position. The crystal was rotated step by step through the reflecting position to record the diffraction profile. Next, the peak positions of both profiles were determined by the extrapolated-peak procedure [§5.3.3.3.1[link], definition (4)] to find the accurate positions of the sample, [\omega_1] and [\omega_2], from which the Bragg angle was calculated by use of a formula that can be written in a simple form as [\theta=|180^\circ- |\omega_1-\omega_2|/2|. \eqno (5.3.3.22)]Before calculating the interplanar distance [equation (5.3.1.1)[link]] or, in the simplest case, the lattice parameter directly, the systematic errors have to be discussed and evaluated. Sometimes, corrections are made to the parameters themselves rather than to the [\theta] values. The reader is referred to §5.3.3.4.3.2[link], in which present knowledge is taken into account, rather than to Bond's original paper.

[Figure 5.3.3.4]

Figure 5.3.3.4 | top | pdf |

Schematic representation of the Bond (1960[link]) method.

Bond performed measurements at room temperature (298 K) for reflections 444, 333, and 111 and, after detailed discussion of errors, reported the [a_0] values (in kXU), which related to these measurements (standard deviations are given in parentheses), as 5.419770 (0.000019), 5.419768 (0.000031), 5.419790 (0.000149). These values are referred to λ = 1.537395 kXU. These results were then tested by Beu, Musil & Whitney (1962[link]) by means of the likelihood-ratio method to test the hypothesis of `no remaining systematic errors'. They proved that the estimate for this sample of silicon is accurate within the stated precision (1 part in 390 000).

The results reported in Bond (1960[link]) – very high accuracy and remarkable reproducibility (low standard deviation), obtained by use of a relatively simple device, which can be realized on the basis of a standard diffractometer – encourage experimenters to perform similar measurements. However, many problems arise with the adaptation of the Bond method to other kinds of samples and/or to other purposes than those described by Bond (1960[link]) in his original paper. Both theoretical and experimental work have increased the accuracy and the precision of the method during the last 35 years.

5.3.3.4.3.2. Systematic errors

| top | pdf |

As mentioned above (§5.3.3.4.1[link]), some systematic errors that affect the asymmetric diffractometer are experimentally eliminated in the Bond (1960[link]) arrangement. According to Beu (1967[link]), who has supplemented the list of errors given by Bond, the following systematic errors are eliminated at the [0.001^\circ\theta] level:

  • (a) absorption, source profile, radial divergence and surface flatness; removed since the detectors are used only to measure intensities and not angular positions;

  • (b) zero, eccentricity, misalignment and diffractometer radius; eliminated since [\theta] depends only on the difference in the crystal-angle positions and not on these geometrical factors;

  • (c) ratemeter recording does not affect the measurements since the detectors are used only for point-by-point counting;

  • (d) 2:1 tracking error is eliminated because the 2:1 tracking used in most commercial asymmetric diffractometers is not used;

  • (e) dispersion, if the peak position of the profile is determined rather than the centroid or the median, and the wavelength has been determined for the peak position also.

As well as these errors there are other systematic errors, due to both physical and apparatus factors, which should be eliminated by suitable corrections.

  • (1) Lorentz–polarization error. The Lorentz–polarization factor Lp influences the shape of the profile as follows: [v_1(x)={\rm Lp} v(x), \eqno (5.3.3.23)]where v(x) and [v_1(x)] are the shape functions [see equations (5.3.3.10)[link], (5.3.3.10a[link],b[link])] of the undistorted and distorted profiles, respectively. It therefore produces a shift [(\Delta\theta_{\rm Lp})] in the peak position.

    The correction for the Lp factor was estimated, assuming that [v(x)] is the Cauchy function [equation (5.3.3.19)[link]], by Bond (1960[link], 1975[link]), Segmüller (1970[link]), and Okazaki & Ohama (1979[link]) for two cases. For perfect crystals, when the Lp factor has the form (James, 1967[link], p. 59; Segmüller, 1970[link]; Okazaki & Ohama, 1979[link]) [{\rm Lp}=(1+|\cos2\theta|)/\sin2\theta, \eqno (5.3.3.24)]the correction is given by [\theta-\theta_p=(\omega_h/2){}^2[\cot2\theta_p+\sin2\theta_p/(1+|\cos2\theta_p|)], \eqno (5.3.3.25)]where [\omega_h] is the half-width of the profile, [\theta_p] is the Bragg angle related to the distorted profile, and [\theta] is the corrected Bragg angle. In contrast, the following formulae are valid for mosaic crystals: [{\rm Lp}=(1+\cos^22\theta)/(2\sin2\theta), \eqno (5.3.3.26)]and [\theta-\theta_p=(\omega_h/2){^2}\cot2\theta_p(2+\sin{^2}2\theta_p)/(2-\sin{^2}\theta_p). \eqno (5.3.3.27)]Because of a notable difference between the values calculated from (5.3.3.25)[link] and (5.3.3.27)[link], the problem is to choose the formulae to be used in practice. However, the Lorentz–polarization error is usually smaller than the rest.

  • (2) Refraction. In the general case, when the crystal surface is not parallel to the reflecting planes but is rotated from the atomic planes around the measuring axis by the angle [\varepsilon], the correction, which relates directly to the determined interplanar distance, has the form (Bond, 1960[link]; Cooper, 1962[link]; Lisoivan, 1974[link], 1982[link]) [d=d_p\left[1+{\delta\cos^2\varepsilon \over \sin(\theta+\varepsilon)\sin(\theta-\varepsilon)}\right], \eqno (5.3.3.28)]where δ is unity minus the refractive index of the crystal for the X-ray wavelength used, and [d_p] and d are the uncorrected and corrected interplanar distances, respectively.

  • (3) Errors due to axial and horizontal (in-plane) divergence. The axial divergence of the primary beam, given by an angle [2\Delta_p] depending on the source and collimator dimensions, causes the angle [\theta'], formed by a separate ray of the beam with a given set of crystallographic planes, to differ from the proper Bragg angle. In general, if the plane of diffraction is not sufficiently perpendicular to the axis of rotation but lacks perpendicularity by an angle [\Delta], the measured Bragg angle [\theta'] can be described, according to Bond (1960[link]), as [\sin\theta'=\sec\Delta\sin\theta. \eqno (5.3.3.29)]Let us assume that both the crystal and the collimator have been accurately adjusted so that the lack of perpendicularity results from axial divergence only. By averaging the expression (5.3.3.29)[link] over the limits [\pm\Delta_p], the mean value of [\sin\theta'] can be found and, as a consequence, the following formula describing the correct d spacing can be obtained: [d=d'(1+\Delta^2_p/6), \eqno (5.3.3.30)]where [d'] is the apparent d spacing.

    According to Berger (1984[link]), this correction is valid only for the case of infinitely small focus, when all rays have the same intensity. Taking into consideration the shift of the centroid caused by vertical divergence when the focus emits uniformly within the axial limits (−F, F), he proposes an alternative correction for [\theta]: [\Delta\theta_d=\textstyle{1\over6}\tan\theta(P^2+F^2), \eqno (5.3.3.31)]where 2P is the sample height.

    As tested using computer modelling (Urbanowicz, 1981b[link]) and estimated analytically (Härtwig & Grosswig, 1989[link]), the effect of the horizontal divergence on the peak position of the recorded profiles cannot be neglected, contrary to suggestions of Bond (1960[link]). The respective systematic error is dependent on asymmetries of both the focus-tube emissivity and the spectral line, and so it is difficult to express it with a simple formula [cf. point (7) below]. In practice (Härtwig, Grosswig, Becker & Windisch, 1991[link]), it proves to be the second largest error. (The first is the one caused by refraction.)

  • (4) Specimen-tilt and beam-tilt error. Since the three main sources of systematic error in diffractometer measurements, i.e. zero, eccentricity, and absorption, have been eliminated in the Bond method, two errors due to misalignment of the crystal and the collimator can strongly influence results of lattice-parameter determination. They are difficult to control because of the random character; numerous authors analysing the Bond method have tried to cope with them. A review is given by Nemiroff (1982[link]).

    Bond (1960[link]) considered the crystal-tilt error separately from the collimator tilt. However, in subsequent papers on this subject it was shown that the errors connected with the crystal tilt and the collimator tilt, i.e. with the angles that the normals to the crystal and collimator make, respectively, with the plane of angular measurement, are dependent and should be treated jointly.

    Foreman (in Baker, George, Bellamy & Causer, 1968[link]) derived a formula for the real value of the angle between two reflecting positions [i.e. [\omega_1] and [\omega_2] in equation (5.3.3.22)[link]] when affected by both tilts. Burke & Tomkeieff (1968[link], 1969[link]), in contrast, have found a dependence between the crystal tilt α and the beam tilt β and the relative error Δa/a in lattice parameter a in the form [\Delta a/a=\alpha\beta/\sin\theta-(\alpha^2+\beta^2)/2. \eqno (5.3.3.32)]A separate analysis is given by Gruber & Black (1970[link]) and by Filscher & Unangst (1980[link]).

    Two approaches are used to eliminate the systematic errors considered, based on the above formula:

    • (i) The error resulting from the crystal tilt and the collimator tilt can be reduced experimentally. Baker, George, Bellamy & Causer (1968[link]) have given a simple procedure that allows a collimator tilt of small but unknown magnitude to be tolerated and, at the same time, the tilt of the crystal to be adjusted to its optimum value. Burke & Tomkeieff (1968[link], 1969[link]) propose a method for setting the crystal so that α = β, since, as is obvious from (5.3.3.32)[link], the error has then its minimum value; α and β have to be of the same sign. Then the influence of crystal tilt and beam tilt on the accuracy of lattice-parameter determination is negligible at the level of 1 part in 106.

    • (ii) Equation (5.3.3.32)[link] permits calculation of the exact correction due to both crystal and collimator tilts, if the respective values of α and β are known. Halliwell (1970[link]) proposed a method for determining the beam and the crystal tilt that requires measuring reflections from both the front and back surfaces of the crystal. In a method described by Nemiroff (1982[link]), the two tilts are measured and adjusted independently within ±0.5 mrad.

  • (5) Errors connected with angle reading and setting. Errors in angle reading and angle setting depend both on the class of the device and on the experimenter's technique. Some practical details are discussed by Baker, George, Bellamy & Causer (1968[link]). Since the angles are measured by counting pulses to a stepping motor connected to a gear and worm, the errors due to angle setting and reading depend on the fidelity with which the gear follows the worm. To diminish errors affected by the gearwheel (notably eccentricity), the authors propose a closed error-loop method, which involves using each part of the gear in turn to measure the angle and averaging the results. In the diffractometer reported in the above paper, there was, originally, an angular error of about +15′′ around the gearwheel, and this can be corrected by means of a cam so that the residual error is reduced to about ±5′′.

    Another example of a high-precision drive mechanism is given by Pick, Bickmann, Pofahl, Zwoll & Wenzl (1977[link]). In the diffractometer described in their paper (see also §5.3.3.7.2[link]), the gear was shown to follow the worm with fidelity even down to 0.01′′ steps, and a drift of ±10% per step was traced to insufficient stability of temperature (±0.15 K).

  • (6) Temperature correction. An error [\Delta d_T] in the lattice parameter d owing to the uncertainty ΔT of the temperature T can be estimated from the formula (Łukaszewicz, Pietraszko, Kucharczyk, Malinowski, Stępień-Damm & Urbanowicz, 1976[link]): [\Delta d_T=d\,\alpha_d\,\Delta T, \eqno (5.3.3.33)]if the thermal-expansion coefficient αd in the required direction is known.

    In the case of the 111 reflection of silicon, for which [\alpha_d\approx2.33\times10^{-6}], to obtain a relative accuracy (precision) of 1 part in 106, the temperature has to be controlled with accuracy (precision) not worse than ±0.05 K if the temperature correction is to be neglected (Segmüller, 1970[link]; Hubbard & Mauer, 1976[link]; Łukaszewicz et al., 1976[link]).

  • (7) Remarks. The above list of corrections, sufficient when the Bond (1960[link]) method is applied under the conditions similar to those described by him (large, perfect, specially cut single crystal; well collimated primary beam; large open detector window) has to be sometimes complemented in the case of different specimens and/or different measurement conditions (§5.3.3.4.3.3[link]). When an asymmetric diffractometer is used, all the systematic errors listed in this section (see also §5.3.3.4.1[link]) must be taken into account.

    Using a complete convolution model of the diffraction profile, Härtwig & Grosswig (1989[link]) were able to derive all known aberrations (and so respective corrections) in a rigorous, analytical way. The analytical expressions given by the authors, though based on some simplifying assumptions, are usually much more complex than the ones shown in points (1)–(6) above. Some coefficients in their equations depend on physical parameters characterizing the particular device and experiment. So, to follow the idea of Härtwig & Grosswig, one must individually consider all preliminary assumptions. As shown by the authors, to achieve the accuracy of 1 part in 107, all aberrations mentioned by them must be taken into account. The most important aberrations prove to be those related to refraction and to horizontal divergence.

5.3.3.4.3.3. Development of the Bond method and its applications

| top | pdf |

The Bond (1960[link]) method, in its first stage, was meant for large, specially cut and set samples. In principle, only one lattice parameter can be determined in one measuring cycle. As has been shown, the method can also be adapted to other samples, with non-cubic symmetry, and to geometries of the illuminated area, different from those used by Bond. This task needs, however, some additional operations and often some additional corrections for systematic errors.

The basic application of the Bond (1960[link]) method, because its geometry reduced several systematic errors, was to absolute lattice-parameter measurements. The method also proved useful in precise investigations of lattice-parameter changes.

Bond-system diffractometers were most often realized in practice on the basis of standard diffractometers under computer control (Baker, George, Bellamy & Causer, 1968[link]; Segmüller, 1970[link]; Pihl, Bieber & Schwuttke, 1973[link]; Kucharczyk, Pietraszko & Łukaszewicz, 1993[link]). Some were designed for special investigations, such as high-precision measurements, [\sigma(d)/d=10^{-7}] (Baker, George, Bellamy & Causer, 1966[link]; Grosswig, Härtwig, Alter & Christoph, 1983[link]; Grosswig et al., 1985[link]; Grosswig, Härtwig, Jäckel, Kittner & Melle, 1986[link]); local measurements at chosen points of a specimen (Lisoivan & Dikovskaya, 1969[link]; Lisoivan, 1974[link], 1982[link]); examination of lattice-parameter changes over a wide temperature range (Łukaszewicz et al., 1976[link], 1978[link]; Okada, 1982[link]); or the effect of high pressure on lattice parameters (Mauer, Hubbard, Piermarini & Block, 1975[link]; Leszczyński, Podlasin & Suski, 1993[link]).

By introduction of synchrotron radiation to a Bond-system diffractometer (Ando et al., 1989[link]), a highly collimated and very narrow beam has been obtained, so lattice-parameter measurements can be accomplished reliably and quickly with a routinely achieved precision of 2 parts in 106; these can be combined with X-ray topographs made in selected areas of the sample.

(1) Crystals with different symmetry. Cooper (1962[link]) used the Bond (1960[link]) diffractometer and method for absolute measurements of lattice parameters of several crystals belonging to various orthogonal systems. Special attention was paid to preparing the samples, i.e. cutting and polishing, to obtain crystal surfaces parallel to the planes of interest. One sample of a given substance was sufficient to find the lattice in the case of cubic crystals but two samples were required for tetragonal and hexagonal systems, and three were necessary for the orthorhombic system. This difficulty increases when non-orthogonal lattices have to be examined. This problem was resolved by Lisoivan (1974[link], 1982[link]), who used very thin single-crystal slabs, which made possible measurements both in reflection and in transmission. Lisoivan (1981[link], 1982[link]), developing his first idea, derived the requirements for a precision determination of all the interaxial angles for an arbitrary system. The coplanar lattice parameters can also be determined in one crystal setting when only reflection geometry is used (Grosswig et al., 1985[link]).

Superlattices can be determined using the system proposed by Bond; a simple method for this purpose was derived by Kudo (1982[link]).

(2) Different sample areas. A separate problem is to adapt the Bond method for measurement of small spherical crystals, commonly used in structure investigations. A detailed analysis of this problem is given by Hubbard & Mauer (1976[link]), who indicate that the effect of absorption and horizontal divergence has to be taken into account if the sample dimensions are less than the cross section of the primary beam. As has been mentioned above (§§5.3.3.4.1[link], 5.3.3.4.3.2[link]), these factors, as well as eccentricity and uncertainty of the zero point, could be neglected in Bond's (1960[link]) experiment. Kheiker (1973[link]) considered systematic errors resulting from the latter two factors when small crystals are used. He proposed a fourfold measurement of the sample position (rather than a twofold one used by Bond), in which `both sides' of a given set of planes are taken into account, so that measurement by the Bond method is performed for two pairs of specimen positions: [\omega_1] and [\omega_2=\omega_1 - 2\theta], and [\omega_3=180^\circ+\omega_1] and [\omega_4=180^\circ+\omega_2]. The corresponding positions of the counter are also determined and used in calculations of the Bragg angle (cf. §5.3.3.4.1[link]). The mean value of the [\theta] angle is not subject to the errors mentioned. A similar idea has been presented by Mauer et al. (1975[link]).

In many practical cases, it is necessary to determine lattice parameters of thin superficial layers. One of the possibilities is to use the Bond method for this purpose. Wołcyrz, Pietraszko & Łukaszewicz (1980[link]) used asymmetric Bragg reflections with small angles of incidence, to reduce the penetration depth of X-rays. This rather simple method permits high accuracy if proper corrections (the formulae are given by the authors) resulting from the dynamical theory of diffraction of X-rays are carefully determined. This method was used to estimate the gradient of the lattice parameter inside diffusion layers. The penetration depth was changed by rotation of the sample. Golovin, Imamov & Kondrashkina (1985[link]) achieved a penetration depth as small as about 1 to 10 nm, using X-ray total-reflection diffraction (TRD) from the planes normal to the surface of the specimen. The sample was oriented in such a way that the conditions for total external reflection were satisfied when the X-ray beam fell on the sample at a small angle of incidence, about 0.5°.

The homogeneity of the crystal in a direction parallel to its surface may be examined by means of local measurements, described by Lisoivan & Dikovskaya (1969[link]) and Lisoivan (1974[link]), in which the goniometer head was specially designed so that the sample could be precisely set and displaced.

(3) Determination of lattice-parameter changes. Baker, George, Bellamy & Causer (1968[link]) have shown that a carefully manufactured and adjusted Bond-system diffractometer (mentioned above) with good stability of environmental conditions (temperature, pressure, power voltage) may be a suitable tool for the investigation of lattice-parameter changes. A static method of thermal-expansion measurement is proposed, in which changes in angle of an in situ specimen due to changes in the lattice parameter with temperature are quickly determined. If it is assumed that the intensity and the shape of the peak have not altered with the change of conditions (cf. the method based on double-crystal diffractometers in §5.3.3.7.1[link]), the change in angle can be determined by intensity measurement alone. The reported precision of the relative measurement is 1 part in 107. Since the shape of the profile may change with the change of conditions, the whole profile must be determined accurately and precisely, so that the whole experiment, consisting of a series of measurements, is time-consuming. The optimization problems resulting from this inconvenience have been discussed above (§5.3.3.3.2[link]; Barns, 1972[link]; Urbanowicz, 1981a[link],b[link]).

In particular, thermal-expansion studies can detect phase transitions and the resulting changes in crystal symmetry (Kucharczyk, Pietraszko & Łukaszewicz, 1976[link]; Kucharczyk & Niklewski, 1979[link]; Pietraszko, Waśkowska, Olejnik & Łukaszewicz, 1979[link]; Horváth & Kucharczyk, 1981[link]; Pietraszko, Tomaszewski & Łukaszewicz, 1981[link]; Keller, Kucharczyk & Küppers, 1982[link]; Åsbrink, Wołcyrz & Hong, 1985a[link],b[link]).

Another group of applications of the Bond method is connected with single-crystal characterization problems (homogeneity, doping, stoichiometry) resulting from technological operations (epitaxy, diffusion, ion implantation) producing changes in lattice spacings, δd/d = 10−2 to 10−5. The examples cited below show a variety of applications.

Stępień, Auleytner & Łukaszewicz (1972[link]) and Stępień-Damm, Kucharczyk, Urbanowicz & Łukaszewicz (1975[link]) examined γ-irradiated NaClO3. The effect of X-ray irradiation on the lattice parameter of TGS crystals in the vicinity of the phase transition was studied by Stępień-Damm, Suski, Meysner, Hilczer & Łukaszewicz (1974[link]). Pihl, Bieber & Schwuttke (1973[link]) dealt with ion-implanted silicon, using a Bond-system diffractometer for local measurements. The effect of silicon doping on the lattice parameters of gallium arsenide was studied by Fewster & Willoughby (1980[link]). Crystal-perfection studies by the Bond method were reported by Grosswig, Melle, Schellenberger & Zahorowski (1983[link]), and Wołcyrz & Łukaszewicz (1982[link]). In the latter paper, the measurements were performed on a superficial single-crystal layer by the use of the geometry described above [paragraph (2)] (Wołcyrz, Pietraszko & Łukaszewicz, 1980[link]). Lattice distortion in LiF single crystals was examined by Dressler, Griebner & Kittner (1987[link]), who used the method of Grosswig et al. (1985[link]) [cf. paragraph (1)]. The use of anomalous dispersion in studies of microdefects was considered by Holý & Härtwig (1988[link]).

5.3.3.4.3.4. Advantages and disadvantages of the Bond method

| top | pdf |

The significant advantages of the Bond (1960[link]) method, such as:

  • (a) very high accuracy;

  • (b) rather high precision;

  • (c) well elaborated analysis of errors;

  • (d) a simple arrangement, which may be realized on the basis of a standard diffractometer with computer control and, if necessary, supplemented with suitable attachment; and

  • (e) variety of applications;

make this method one of the most popular at present.

The method, however, has the following limitations:

  • (i) Special requirements concerning the sample are difficult to satisfy in some cases.

  • (ii) Problems arise with determination of all the lattice parameters of non-cubic crystals. Multiple-sample preparation or a special approach is needed in such cases.

  • (iii) Lattice-spacing determination from small spherical crystals requires additional corrections or fourfold measurements.

  • (iv) Displacement of the irradiated area on the sample surface (Wołcyrz, Pietraszko & Łukaszewicz, 1980[link]; Berger, 1984[link]) complicates examination of the real structure (for example, by local measurements).

  • (v) The method is rather time-consuming, since twofold scanning of the profile is required for determination of a single [\theta] value.

  • (vi) Because two detectors, or a wide range of rotations of only one detector, are required, measurement with additional attachments is more difficult than on an asymmetric diffractometer.

Nevertheless, the geometry proposed by Bond (1960[link]), owing to its advantages, is commonly used in precise and accurate multiple-crystal spectrometer methods (§§5.3.3.7.1[link], 5.3.3.7.2[link]).

Other limitations concerning the precision and accuracy of the method are common to it and to all the `traditional' methods (Subsection 5.3.3.5[link]).

5.3.3.5. Limitations of traditional methods

| top | pdf |

As `traditional' are considered the methods that depend on a comparison of the lattice spacings to be determined with the wavelength values of characteristic X-radiation that comes directly from laboratory (Bremsstrahlung) sources. The emission lines are wide and asymmetric, which limits both the accuracy and precision of lattice-parameter measurements (as discussed in Subsection 5.3.1.1[link]). One of the limiting factors is the uncertainty of the wavelength value. For many years, the wavelength values determined by Bearden (1965[link], 1967[link]) with an accuracy of 5 parts in 106 were widely used. At present, owing to remarkable progress in the measurement technique, it is possible to achieve an accuracy in wavelength of an order better, and nowadays remeasurements of some characteristic emission X-ray wavelengths are reported [cf. §5.3.3.3.1(iii)[link] and Subsection 5.3.3.8[link]]. Yet, even after reducing the uncertainty in wavelength, and after introducing all necessary corrections for systematic errors, the highest accuracy of traditional methods does not exceed 1 part in 106 (cf. Subsection 5.3.3.8[link]).

The accuracy of an order better is possible with X-ray and optical interferometry. This non-dispersive method (cf. Sub­section 5.3.3.8[link]) is used for accurate lattice-spacing determination of highly perfect standard crystals; the standards are next used for both lattice-parameter determination with a double-beam comparison technique (Baker & Hart, 1975[link]; see also §5.3.3.7.3[link]) and for the accurate wavelength determination mentioned above.

Another problem is the limited precision attainable by traditional methods. As was discussed in Subsection 5.3.1.1[link], the width of the diffraction profile depends on the spectral distribution of the wavelength, (5.3.1.6)[link], (5.3.1.7)[link], (5.3.1.8)[link], and cannot be less than this owing to the wavelength dispersion. However, much has been done to approach this limit and to attain the precision and accuracy of the diffraction profile location (cf. Subsection 5.3.3.3[link]). The highest precision of lattice-parameter determination that it is possible to achieve with traditional methods is about 1 part in 107. For some problems connected with single-crystal characterization, such as the effect of irradiation, stress, defect concentration, including local measurement (topography), better precision is required.

From (5.3.1.9)[link], the other possibility of increasing precision, besides choosing optimum parameters for the measurement and improvement of profile-location methods, is to influence the original profile [h_\lambda(\omega)]. This aim can be attained either by applying spectrally narrower X-ray sources or by reducing the width of the original profile by means of arrangements with additional crystals playing the role of monochromator and reference crystal. This second possibility is applied in double- or triple-crystal spectrometry, in multiple-beam methods, or in combined methods. These methods are called `pseudo-non-dispersive' methods, since the width of the diffraction profile is considerably limited in them owing to considerable limitation of the width of the original profile. A similar situation to that in n-crystal spectrometers, in which the beam reflected from one set of crystal planes is the source of radiation for the second (or the next) diffraction phenomena, arises in multiple-diffraction methods; this is described in Subsection 5.3.3.6[link].

A systematic and well illustrated review of pseudo-non-dispersive and other differential methods is given by Hart (1981[link]), who is the author of numerous papers on this subject.

5.3.3.6. Multiple-diffraction methods

| top | pdf |

Multiple diffraction occurs when two or more sets of planes simultaneously satisfy the Bragg law for a single wavelength λ. The beam diffracted from one set of planes becomes the incident beam within the crystal for the next diffraction. In the reciprocal-space representation, this means that three or more reciprocal-lattice points lie simultaneously on the Ewald sphere (Fig. 5.3.3.5[link] ). These points can be detected by successive rotations of the crystal, as described below. This phenomenon, known also as simultaneous reflection or (after Renninger, 1937[link]) `Umweganregung', may be observed in both X-ray and neutron experiments. In the first case, it occurs both in methods with counter recording, initiated by Renninger (1937[link]), and in methods with photographic recording in which a highly divergent beam is used (§5.3.2.4.2[link]). The intersections of conic sections encountered in the methods developed by Kossel (1936[link]) and Lonsdale (1947[link]) are the cases of multiple-diffraction phenomena in photographic methods.

[Figure 5.3.3.5]

Figure 5.3.3.5 | top | pdf |

Schematic representation of multiple diffraction in reciprocal space (after Post, 1975[link]).

Simultaneous reflection, undesirable in some cases (`forbidden' reflections in measurements of intensities) can be very useful in others. Its various applications have been reviewed by Terminasov & Tuzov (1964[link]) and Chang (1984[link]). Only the utility of multiple diffraction in lattice-parameter determination will be discussed here.

The principle of Renninger's (1937[link]) experiment, which is also the basis of the method described by Post (1975[link]), is shown in Fig. 5.3.3.6[link] . The crystal in the shape of a slab is at first set in a position to diffract the primary X-ray beam. A primary reflection whose intensity is very low or which is forbidden by the space group of the crystal is usually selected. Its intensity determines the background intensity of the pattern, which should be low. The detector, with a wide-open window, is situated in the appropriate position and remains fixed throughout the experiment while the crystal is rotated around the axis perpendicular to the crystal planes (and its surface) to record successive reflections.

[Figure 5.3.3.6]

Figure 5.3.3.6 | top | pdf |

Schematic representation of the multiple-diffraction method. (a) Experimental set-up (after Cole, Chambers & Dunn, 1962[link]; Post, 1975[link]). (b) Geometric representation in reciprocal space.

The multiple-diffraction pattern (an example is shown in Fig. 5.3.3.7[link] ) has next to be indexed. The principle of the method of indexing, known as the reference-vector method (Cole, Chambers & Dunn, 1962[link]; Post, 1975[link]; Chang, 1984[link]), is shown in Fig. 5.3.3.6(b)[link]. Directions of the primary and diffracted beams are marked by vectors [{\bf K}_0] and K. The ends of the vectors lie on the Ewald sphere, the radius of which is equal to 1/λ. The reciprocal vector [{\bf P}=h_0{\bf a}^*+k_0{\bf b}^*+l_0{\bf c}^*], being the difference between the vectors K and [{\bf K}_0], represents the first diffraction phenomenon, which is observed for setting angles equal to [\varphi_0] and μ (usually [\mu=\theta]). Let us assume that the reciprocal vector [{\bf H}=h{\bf a}^*+k{\bf b}^*+l{\bf c}^*], observed for setting angles [\varphi_0+\beta] and ν, represents the next diffraction. The vector components of H, parallel and normal to P, are denoted by [{\bf H}_p] and [{\bf H}_n], respectively. For a given wavelength, the lengths [P,H, H_p,H_n] of respective vectors [{\bf P},{\bf H},{\bf H}_p, {\bf H}_n] are functions of lattice parameters and diffraction indices.

[Figure 5.3.3.7]

Figure 5.3.3.7 | top | pdf |

The multiple-diffraction pattern at the 222 position in germanium (Cole, Chambers & Dunn, 1962[link]).

The task is to find the relationship between the difference β of angles of rotation (or between two values of the setting angles, β and ν) and the lengths of the reciprocal vectors. The following relations result from Fig. 5.3.3.6(b)[link]: [\eqalign{ \cos\beta &={(C'A'){^2}+H^2_n-R'^2 \over 2H_n(C'A')}, \cr (C'A'){^2} &=R^2-P^2/4, \cr H_p &=P/2-R\sin \nu, \cr R' &=R\cos \nu.}]Taking these into consideration, and remembering that [H^2=H^2_p+H^2_n], we finally obtain [\cos\beta ={H^2-H_pP \over 2H_n(R^2-P^2/4){}^{1/2}}. \eqno (5.3.3.34)]

Since cos(−β) = cosβ, the appearance of successive reflections does not depend on the direction of rotation. A detailed discussion of (5.3.3.34)[link] is given by Cole, Chambers & Dunn (1962[link]) and Chang (1984[link]). When preliminary values of the lattice parameters are known, (5.3.3.34)[link] can be applied for indexing multiple-diffraction patterns. A computer program (Rossmanith, 1985[link]) can be very useful in rather complicated calculations and in the graphical representation of the multiple-diffraction pattern. Formula (5.3.3.34)[link], since [C'A'=R\cos\mu], can be presented in another form: [R={H^2-H_pP \over 2H_n\cos \mu\cos\beta}, \eqno (5.3.3.35)]where two setting angles, μ and β, are taken into account. When the indices are known, both (5.3.3.34)[link] and (5.3.3.35)[link] can be used for the determination or refinement of lattice parameters.

Another analytical method for indexing multiple-diffraction patterns, based on the determination of the Lorentz point, has been described by Kshevetsky, Mikhailyuk, Ostapovich, Polyak, Remenyuk & Fomin (1979[link]).

Formulae (5.3.3.34)[link] and (5.3.3.35)[link] are valid for all crystal systems. In practice, however, the rather complicated method is used mainly for cubic crystals, and a special approach proved to be needed in order to adapt the method to other (rectangular) systems (Kshevetsky, Mikhalchenko, Stetsko & Shelud'ko, 1985[link]). In the case of a cubic lattice, it is convenient to substitute [R=a/\lambda \eqno (5.3.3.35a)]into (5.3.3.34)[link] and (5.3.3.35)[link] rather than R = 1/λ used in the general case, so that the lengths of the reciprocal vectors, being now functions of the indices only, are: [\eqalignno{ P &=(h^2_0+k^2_0+l^2_0){}^{1/2}, &(5.3.3.35b) \cr H&=(h^2+k^2+l^2){}^{1/2}, & (5.3.3.35c) \cr H_p P&={\bf HP} = h_0h+k_0k+l_0l. &(5.3.3.35d)}%fd5.3.3.35c,fd5.3.3.35d]The lengths of the components of H can be determined from (5.3.3.35b,c,d[link]) taking [H_p=H_pP/P] and [H_n=(H^2-H^2_p){}^{1/2}]. After introducing the alterations [equation (5.3.3.35a[link]) and the resulting equations (5.3.3.35b,c,d[link])], (5.3.3.35)[link] now describes a simple dependence between the ratio a/λ, the indices, and the setting angles.

The accuracy of the lattice-parameter determination resulting from (5.3.3.35)[link] in the cubic case can be assumed to be: [{\Delta a\over a}=\tan\mu\Delta\mu+\tan\beta\Delta\beta; \eqno (5.3.3.36)]this thus depends on the values of the setting angles β, μ and their accuracies [\Delta\beta], [\Delta\mu]. The latter depend on various systematic errors.

Since the differences between the two angular settings at which a given set of planes diffracts are measured rather than their absolute values, the systematic errors due to absorption, specimen displacement, and zero-setting are eliminated. In contrast, errors due to vertical divergence, refraction and the change of wavelength of the incident radiation (when it enters the crystal), alignment, and dynamical effects should be taken into account. In the case described by Post (1975[link]), when a fine focus (effective size 0.4 × 0.5 mm) and collimation limiting the beam divergence to 2′ were used, the vertical divergence causing the relative error in d of about [5\times10^{-8}] could be ignored.

The errors due to the real structure (inhomogeneity, mosaicity and internal stress) were discussed by Kshevetsky et al. (1979[link]).

The accuracy possible by this method (from 1 to 4 parts in 106) is comparable with that obtained with the Bond (1960[link]) method. The advantages of this method from the point of view of lattice-parameter determination are as follows:

  • (a) a large number of reflections can be measured without realigning or removing the crystal;

  • (b) all the lattice parameters can be determined and not only one, as in the Bond (1960[link]) method;

  • (c) the narrow diffraction profiles can be located with very high accuracy and precision;

  • (d) the arrangement makes it possible to remove some systematic errors;

  • (e) the high accuracy resulting from (a)–(d), which is comparable with that obtained by means of the Bond (1960[link]) method;

  • (f) the high precision that results from (a) and (c).

A disadvantage, on the other hand, is the complicated interpretation (indexing) of multiple-diffraction patterns, so that this method is less popular than the Bond (1960[link]) method.

The Post (1975[link]) method has been applied to the accurate lattice-parameter determination of germanium, silicon, and diamond single crystals (Hom, Kiszenick & Post, 1975[link]).

5.3.3.7. Multiple-crystal – pseudo-non-dispersive techniques

| top | pdf |

5.3.3.7.1. Double-crystal spectrometers

| top | pdf |

Detailed information concerning the double-crystal spectrom­eter, which consists of two crystals successively diffracting the X-rays, can be found in James (1967[link], pp. 306–318), Compton & Allison (1935[link]), and Azároff (1974[link]). This device, usually used for wavelength determination, may also be applied to lattice-parameter determination, if the wavelength is accurately known. The principle of the device is shown in Fig. 5.3.3.8[link] . The first crystal, the monochromator, diffracts the primary beam in the direction defined by the Bragg law for a given set of planes, so that the resulting beam is narrow and parallel. It can thus be considered to be both a collimator (or an additional collimator, if the primary beam has already been collimated) and a wavelength filter. The final profile [h(\theta)], obtained as a result of the second diffraction by the specimen when the first crystal remains stationary and the second is rotated, is narrower than that which would be obtained with only one crystal. The final crystal profile [h_C(\theta)] [cf. equation (5.3.1.6)[link]] is due to both crystals, which, if it is assumed that they are cut from the same block, can be described by the autocorrelation function (Hart, 1981[link]): [h_C(\theta)=K\textstyle\int\limits^\infty_{-\infty}\,R(\theta')R(\theta'-\theta)\,{\rm d}\theta', \eqno (5.3.3.37)]where [R(\theta)] is an individual reflectivity function of one crystal and K is a coefficient of proportionality. Its half-width is 1.4 times larger than that related to only one crystal. In spite of this, the recorded profile can be as narrow as, for example, 2.6′′ (Godwod, Kowalczyk & Szmid, 1974[link]), since the profile due to the wavelength [h_\lambda(\theta)], modified by the first crystal, is extremely narrow. Additional advantages of the diffraction profile are: its symmetry, because [h_C(\theta)] is symmetric as an effect of autocorrelation, and smoothness, as an effect of additional integration. The profile can thus be located with very high accuracy and precision.

[Figure 5.3.3.8]

Figure 5.3.3.8 | top | pdf |

Schematic representation of the double-crystal spectrometer.

When there is a small difference in the two lattice spacings, so that one has a value d and the other [d+\delta d], if [\delta\lambda/\lambda] is small enough, it can be assumed that the profile does not alter in shape but in its peak position [cf. §5.3.3.4.3.3[link], paragraph (3)]. If for two identical crystals this were located at [\theta_0], the peak position shifts to [\theta_0-\tan\theta\delta d/d]. The measurement of this shift rather than the absolute position of the rocking curve is the basis of all the double-crystal methods. An example of the application of a double-crystal spectrometer with photographic recording has been given in §5.3.2.3.5[link] (Bearden & Henins, 1965[link]).

The basic requirements that should be fulfilled to make the most of the double-crystal spectrometer are: limitation of the primary beam by means of a collimator, parallelism of the two axes [precision as high as 1′′ obtained by Godwod, Kowalczyk & Szmid (1974[link])], and high thermal stability (0.1 K; Godwod, Kowalczyk & Szmid, 1974[link]). Alignment procedure, errors, and corrections valid for the double-crystal spectrometer have been considered by Bearden & Thomsen (1971[link]).

The double-crystal diffractometer, because of the small width of the diffraction profile, is a very suitable tool for local measurements of lattice-parameter differences, for example between an epitaxic layer and its substrate. Hart & Lloyd (1975[link]) carried out such a measurement on a standard single-axis diffractometer (APEX) to which a simple second axis, goniometer head, and detector were added (Fig. 5.3.3.9[link] ). The diffracted beam was recorded simultaneously by three detectors. A symmetric arrangement with two detectors, [D_1] and [D_2], with no layer present, makes possible the determination of the absolute value of the lattice parameter of the substrate, as in the Bond (1960[link]) method. The third detector makes it possible to record the double-crystal rocking curve, which usually fully resolves the layer and substrate profiles. The changes in the lattice parameter between the two components can be used for determination of strain (at 1 part in 104).

[Figure 5.3.3.9]

Figure 5.3.3.9 | top | pdf |

Schematic representation of the double-crystal arrangement of Hart & Lloyd (1975[link]) for the examination of epitaxic layers. (a) Experimental set-up. (b) Diffraction profiles recorded by detectors D1, D2, and D3.

The very important advantage of this method, from the point of view of local measurements, is that single- or double-crystal diffraction can be selected, simultaneously if needed, on exactly the same specimen area. Other examples of strain measurements by means of a double-crystal spectrometer are given by Takano & Maki (1972[link]), who measured lattice strain due to oxygen diffusing into a silicon single crystal; by Fukahara & Takano (1977[link]), who compared experimental rocking curves and theoretical ones computed within the frame of the dynamical theory; and Barla, Herino, Bomchil & Pfister (1984[link]), who examined the elastic properties of silicon.

The standard double-crystal technique does not allow determination of relatively small strains, i.e. ones that affect the lattice parameter by, for example, less than 2–3 parts in 105, as in the case of (004) Si reflection and Cu Kα radiation. To overcome this difficulty, Zolotoyabko, Sander, Komem & Kantor (1993[link]) propose a new method that combines double-crystal X-ray diffraction with high-frequency ultrasonic excitation. Since ultrasound has a wavelength a little less than the X-ray excitation length, it affects the diffraction profile close to the Bragg position and so permits the detection of very small profile broadenings caused by lattice distortions. With this method, lattice distortion as small as 5 parts in 106 can be measured.

As has been shown in the case of the device used by Hart & Lloyd (1975[link]), the symmetric arrangement due to Bond (1960[link]) proves to be very useful when the double-crystal spectrometer is to be used for absolute lattice-parameter determination, since such an arrangement combines the high precision and sensitivity of a double-crystal spectrometer with the high absolute accuracy of the Bond method. Other examples of a similar idea are presented by Kurbatov, Zubenko & Umansky (1972[link]), who report measurements of the thermal expansion of silicon; Godwod, Kowalczyk & Szmid (1974[link]), who also discuss the theoretical basis of their arrangement; Ridou, Rousseau & Freund (1977[link]), who examine a phase transition; Sasvári & Zsoldos (1980[link]), and Fewster (1982[link]). The latter two papers are concerned with epitaxic layers. A rapid method is proposed by Sasvári & Zsoldos (1980[link]) for deconvoluting the overlapping peaks due to the layer and the substrate. A particular feature of the arrangement proposed in the first of these papers (Kurbatov, Zubenko & Umansky, 1972[link]) is the use of a germanium-crystal monochromator with anomalous transmission, to obtain a nearly parallel primary beam (the horizontal divergence is 28′′ and the vertical 14′′).

The error analyses given by Godwod, Kowalczyk & Szmid (1974[link]) and Sasvári & Zsoldos (1980[link]) show that systematic errors due to eccentricity, absorption, and zero position are eliminated experimentally, owing to the symmetric arrangement, as in the Bond (1960[link]) method. In contrast, the errors due to crystal tilt, refraction and the Lorentz–polarization factor [their uncertainties in lattice parameters, as evaluated by Sasvári & Zsoldos (1980[link]), are 10−6 Å each], axial divergence (2 × 10−6 Å), angle reading (10−4 Å), and instrument correction and calculations (each to 5 × 10−5 Å) should be taken into account. The effect of absorption, discussed by Kurbatov, Zubenko & Umansky (1972[link]), proved to be negligible. The final accuracy achieved for silicon single crystals by Godwod, Kowalczyk & Szmid (1974[link]) is comparable with that obtained by Bond (1960[link]).

A specific group of double-crystal arrangements is formed by those in which white X-radiation is used instead of characteristic. Such an arrangement makes possible very large values of the Bragg angle (larger than about 80°), which increases the accuracy, precision, and sensitivity of measurement of the lattice parameters and their change with change of temperature. This task is rather difficult to realize by means of traditional methods, in which both the wavelengths and the lattice parameters are fixed, and it is difficult to find a suitable combination of their values.

The principle of the method presented by Okazaki & Kawaminami (1973a[link]) is shown in Fig. 5.3.3.10[link] . The first crystal (the specimen to be measured) remains fixed during a single measurement, the second (the analyser) is mounted on the goniometer of an X-ray diffractometer and can be operated with either an [\omega] or a [\theta][2\theta] scan. As diffraction phenomena appear for both the specimen and the analyser (in general of different materials) whose interplanar spacings are equal to [d_s] and [d_A], respectively, the following relation results from Bragg's law: [d_s\sin\theta_s=d_A\sin \theta_A, \eqno (5.3.3.38)]where [\theta_s] and [\theta_A] are the respective Bragg angles. Since [d_A] and [\theta_s] are kept constant, a change in [d_s] as a function of temperature is determined from a change in [\theta_A]. The relative error [\delta d/d] resulting from (5.3.3.38)[link] with [\theta_A\approx90^\circ] is [\eqalignno{ {\delta d_s \over d_s} = \cot\theta_A\delta\theta_A &= \tan(\pi/2-\theta_A)\delta\theta_A \cr &\approx (\pi/2-\theta_A)\delta\theta_A. &(5.3.3.39)}]The method initiated by Okazaki & Kawaminami (1973a[link]) has been developed by Okazaki & Ohama (1979[link]), who constructed the special diffractometer HADOX (the positions of the specimen and the analyser were interchanged) and discussed systematic errors. Precision as high as 1 part in 107 was reported. Examples of the application of such an arrangement for measuring the temperature dependence of lattice param­eters were given by Okazaki & Kawaminami (1973b[link]) and Ohama, Sakashita & Okazaki (1979[link]). Various versions of the HADOX diffractometer are still reported. By introducing two slits (Soejima, Tomonaga, Onitsuka & Okazaki, 1991[link]) – one to limit the area of the specimen surface to be examined and the other to define the resolution of [2\theta] – it is possible to combine [\omega] and [2\theta] scans and obtain a two-dimensional intensity distribution in the plane parallel to the plane of the diffractometer, and to determine the temperature dependence of lattice parameters on a selected area of the specimen (avoiding the effects of the surroundings). The HADOX diffractometer may work with both a rotating-anode high-power X-ray source (examples reported above) and a sealed-tube X-ray source. In the latter case (Irie, Koshiji & Okazaki, 1989[link]), to increase the efficiency of the X-ray tube, the distance between the X-ray source and the first crystal has been shortened by a factor of five. As is implied by (5.3.3.39)[link], one can increase the relative precision of the method by using the analyser angle close to π/2. This idea has been realized by Okazaki & Soejima (2001[link]), who achieved the relative accuracy of determination of lattice-parameter changes as high as 1 part in [10^{9}][10^{10}] by extending the Bragg angle from 78° (previous versions) to 89.99° and by elimination of systematic errors due to crystal tilt, crystal displacement, temperature effects and radiation damage.

[Figure 5.3.3.10]

Figure 5.3.3.10 | top | pdf |

Schematic representation of the double-crystal arrangement of Okazaki & Kawaminami (1973a[link]); white incident X-rays are used.

An original method for the measurement of lateral lattice-parameter variation by means of a double-crystal arrangement with an oscillating slit was proposed by Korytár (1984[link]). This method permitted simultaneous recording of two rocking curves from two locations on a crystal. Precision of 3 parts in 107 was reported. The method has been applied for the measurement of growth striations in silicon.

The main disadvantage of double-crystal spectrometers, in their basic form (Fig. 5.3.3.8[link]), is that they cannot be used for measurements on an absolute scale. Combination of the double-crystal arrangement with the system proposed by Bond (1960[link]) makes it possible to recover the origin of the angular scale and thus such an absolute measurement, but the reported precision is rather moderate.

There are two other ways to overcome this difficulty in pseudo-non-dispersive methods: addition either of a third crystal (more accurately, a third reflection) (§5.3.3.7.2[link]) or of a second source (a second beam) (§5.3.3.7.3[link]). Such arrangements require additional detectors. Combinations of both techniques are also available (§5.3.3.7.4[link]).

5.3.3.7.2. Triple-crystal spectrometers

| top | pdf |

Higher precision than that obtained with the double-crystal arrangements (§5.3.3.7.1[link]) can be achieved by means of triple-crystal diffractometers. Arrangements specially designed for the determination of lattice-parameter changes are described by Buschert (1965[link]) and Skupov & Uspeckaya (1975[link]), and reviewed by Hart (1981[link]).

The principle of the triple-axis spectrometer is shown in Fig. 5.3.3.11[link] . The arrangement consists of one standard crystal S, ultimately replaced by the sample under investigation, and two reference crystals [R_1] and [R_2]. The principle of the measurement is as follows. First, the crystals S, [R_1], and [R_2] are set to their diffraction (peak) positions using two detectors [D_1] and [D_2]. Then the standard crystal S is replaced by the sample and the new peak position is found by means of [D_1] when the sample is turned from its original position to its reflecting position. The angle of rotation of the sample [\Delta\omega_s] depends on the lattice-parameter difference [\Delta d] between the sample and the standard. The relation is given by (Hart, 1981[link]) [\Delta\omega_s=-\tan\theta\Delta d/d. \eqno (5.3.3.40)]Next, the second reference crystal [R_2] is turned through the angle [\Delta\omega_R] to its diffracting position, the intensity being controlled with the second detector [D_2]. From the geometry of the arrangement, [\Delta\omega_R=2\Delta\omega_s. \eqno (5.3.3.41)]Because the origin of the [\omega_s] scale is lost during the crystal exchange, this second angle of rotation [(\Delta\omega_R)] is used to determine [\Delta d] rather than the first one [(\Delta\omega_s)], by using (5.3.3.41)[link] and (5.3.3.40)[link].

[Figure 5.3.3.11]

Figure 5.3.3.11 | top | pdf |

Schematic representation of the triple-crystal spectrometer developed by Buschert (1965[link]) (after Hart, 1981[link]).

The diffraction profiles observed in the second detector, described by Hart (1981[link]), [h(\theta)_R=\int\limits^\infty _{-\infty}\,R^2(\theta') R(\theta'-\theta)\,{\rm d}\theta, \eqno (5.3.3.42)]are not symmetric but can be as narrow as 0.1–1′′, so that a precision of 2 parts in 108 is possible.

The main experimental problem here is to adjust the tilts of the crystals. The errors resulting both from the crystal tilts and from the vertical divergence were discussed by Skupov & Uspeckaya (1975[link]).

Triple-crystal spectrometers are often applied as lattice-spacing comparators, when very small changes of lattice parameters [(10^{-8}\le|\Delta d|/][d\le10^{-6})] are to be detected, in particular for the examination of a correlation between lattice parameter and the dopant or impurity concentration (Baker, Tucker, Moyer & Buschert, 1968[link]). Such an arrangement can also be a very suitable tool in deformation studies, since it allows the separation of the effect of deformation on the Bragg angle from that due to lattice-parameter change (Skupov & Uspeckaya, 1975[link]).

The basis of the accurate lattice-parameter comparison proposed by Bowen & Tanner (1995[link]) is the use of a high-purity silicon standard (cf. §5.3.3.9[link] below) with a well known lattice parameter. To compensate an error that may result from a slight misalignment of crystal planes in relation to the axes of the instrument, the authors recommend a twofold measurement of the diffraction-peak position of the reference crystal (for a given diffraction position and after rotating the specimen holder through 180° about the axis normal to its surface) and a similar twofold measurement of the diffraction-peak position of the sample – after replacing the reference crystal by the sample. The mean positions of the reference crystal and of the sample are used in calculations of the Bragg-angle difference and then of the unknown interplanar spacing. The method uses a standard double-crystal diffractometer fitted with a monochromator (therefore, a third crystal), which provides a well defined wavelength, and with a specimen rotation stage. The measurement is accompanied by a detailed error analysis. The accuracy of absolute lattice-parameter determination as high as a few tens of parts in [10^{6}], and a much greater relative sensitivity are reported.

By combining a triple-axis spectrometer with the Bond (1960[link]) method, the device can be used for absolute measurements (Pick, Bickmann, Pofahl, Zwoll & Wenzl, 1977[link]). The device described in the latter paper is an automatic triple-crystal diffractometer that permits intensity measurement to be made in any direction in reciprocal space in the diffraction plane with step sizes down to 0.01′′ and therefore can be used for very precise measurements [see also §5.3.3.4.3.2[link], paragraph (5)].

5.3.3.7.3. Multiple-beam methods

| top | pdf |

The other possibility of recovering the crystal-angle scale in differential measurements with a double-crystal spectrometer (cf. §§5.3.3.7.1[link], 5.3.3.7.2[link]) is to obtain reflections from two crystal planes [for example, from (hkl) and [(\bar h\bar k\bar l)] planes] by means of a double-beam arrangement and to measure them simultaneously.

The second X-ray beam may come from an additional X-ray source (Hart, 1969[link]) or may be formed from a single X-ray source by using a beam-splitting crystal (Hart, 1969[link], second method; Larson, 1974[link]; Cembali, Fabri, Servidori, Zani, Basile, Cavagnero, Bergamin & Zosi, 1992[link]). In particular, two beams with different wavelengths [(K\alpha_1,K\beta_1)] separated with a slit system can be used for this purpose (Kishino, 1973[link], second technique). The principle of the double-beam method is shown in Fig. 5.3.3.12[link] . The beams are directed at the first crystal (the reference crystal) so that the Bragg condition is simultaneously fulfilled for both beams, and they then diffract from the second crystal (the specimen). As the second crystal is rotated, a double-crystal diffraction profile is recorded first in one detector and then in the other. The angle [\Delta\theta] of crystal rotation between the two rocking curves is given by (Baker & Hart, 1975[link]): [\Delta\theta=(\theta_1-\theta_2)=\tan\theta\Delta d/d. \eqno (5.3.3.43)]This formula leads to the lattice-parameter changes Δd.

[Figure 5.3.3.12]

Figure 5.3.3.12 | top | pdf |

Schematic representation of the double-beam comparator of Hart (1969[link]).

A double-beam diffractometer can be used for the examination of variations in lattice parameters of about 10 parts in 106 within a sample in a given direction. An example was reported by Baker, Hart, Halliwell & Heckingbottom (1976[link]), who used Larson's (1974[link]) arrangement for this task.

The highest reported sensitivity (1 part in 109) can be achieved in the double-source double-crystal X-ray spectrometer proposed by Buschert, Meyer, Stuckey Kauffman & Gotwals (1983[link]). The device can be used for the investigation of small concentrations of dopants and defects.

The method can also be applied for the absolute determination of a lattice parameter, if that of the reference crystal is accurately known and the difference between the two parameters is sufficiently small. Baker & Hart (1975[link]), using multiple-beam X-ray diffractometry (Hart, 1969[link], first technique), determined the d spacing of the 800 reflection in germanium by comparing it with the d spacing of the 355 reflection in silicon. The latter had been previously determined by optical and X-ray interferometry (Deslattes & Henins, 1973[link]; the method is presented in Subsection 5.3.3.8[link]).

In the case of two different wavelengths and diffraction from two different diffraction planes [(h_1k_1l_1)] and [(h_2k_2l_2)], the lattice parameter [a_0] of a cubic crystal can be determined using the formula (Kishino, 1973[link]) [a_0=\textstyle{1\over2}\{(L\lambda_1)^2+[(M\lambda_2-L\lambda_1\cos\theta_{1- 2})/\sin\theta_{1 - 2}]^2\}^{1/2}, \eqno (5.3.3.44)]where [L=(h^2_1+k^2_1+l^2_1){}^{1/2}], [M=(h^2_2+k^2_2+l^2_2){}^{1/2}], and [\theta_{1 - 2}] is the difference between the two Bragg angles for the specimen crystal, estimated from the measurement of [\Delta\theta=|\theta_{1 - 2}-\theta'_{1 - 2}|] if the difference [\theta_{1 - 2}'] for the first (reference crystal) is known beforehand. The idea of Kishino was modified by Fukumori, Futagami & Matsunaga (1982[link]) and Fukumori & Futagami (1988[link]), who used the Cu Kα doublet instead of [K\alpha_1] and [K\beta_1] radiation. Owing to the change, they could use only one detector (Kishino's original method needs two detectors), but a special approach is sometimes needed to resolve two peaks that relate to the components of the doublet. A similar problem of separation of two peaks (recorded by two detectors) is reported by Cembali et al. (1992[link]). By introducing a computer simulation of the reflecting curves (using a convolution model), the authors managed to determine the separation with an error of 0.01′′ and to achieve a precision of some parts in 107. The same precision is reported by Fukumori, Imai, Hasegawa & Akashi (1997[link]), who introduced a precise positioning device and a position-sensitive proportional counter to their instrument.

As in the other multiple-crystal methods, the most important experimental problem is accurate crystal setting. Larson (1974[link]), as a result of detailed analysis, gave the dependence between the angular separation of two peaks and angles characterizing misalignment of the first and second crystals.

5.3.3.7.4. Combined methods

| top | pdf |

The idea of multiple-beam measurement (§5.3.3.7.3[link]) can be applied to other arrangements that combine the features of the double-beam comparator with those of the triple-crystal spectrometer; there are additional advantages in such a system. The application of the double-beam technique makes it possible to realize a triple-reflection scheme for comparing lattice param­eters on the basis of a double-axis spectrometer. The arrangement proposed by Ando, Bailey & Hart (1978[link]), shown in Fig. 5.3.3.13[link] , consists of a sample and a reference crystal, which are made from the same material but differ in purity (or strain, stoichiometry, vacancy concentration, etc). The angular difference [\Delta\theta] in the Bragg angles of the sample and the reference crystal, [\theta_S] and [\theta_R], respectively, [\Delta\theta=\theta_S-\theta_R, \eqno (5.3.3.45a)]is measured as the sample angle [\Delta\omega] between the double-reflected peak D and the triply diffracted peak T: [\eqalignno{ \Delta\theta &= \Delta\omega, &(5.3.3.45b) \cr \Delta\omega&=\omega_D-\omega_T. &(5.3.3.45c)}%fd5o3o3o45c]Assuming that [\Delta\theta] is entirely due to changes [\Delta d] in atomic spacings, the authors use the following relation for determination of the latter: [\Delta d/d=-\cot\theta\Delta\theta. \eqno (5.3.3.46)]The experimental requirements are simple and inexpensive, owing to simple shapes of both the reference crystal and the sample crystal, so that the measurement can be made quickly. By combining the two reference crystals into a single monolithic reference crystal, excellent stability, difficult to achieve with triple-axis arrangements (cf. §5.3.3.7.2[link]), is obtained at the same time. The disadvantage of the method is that it covers a smaller range of lattice parameters than the other double-beam methods (Hart, 1969[link]; Larson, 1974[link]) described in §5.3.3.7.3[link]. A new version of the double-crystal triple-reflection scheme (Häusermann & Hart, 1990[link]) allows one to achieve a precision of 1 part in 108 in 2 min of measurement time, which includes the data analysis; 30 min are needed to change the sample. Errors due to the crystal tilt and thermal drifts are considered.

[Figure 5.3.3.13]

Figure 5.3.3.13 | top | pdf |

The double-axis lattice-spacing comparator of Ando, Bailey & Hart (1978[link]); a triple-diffracted beam is used.

Another example of the triple-reflection scheme realized by means of the double-beam technique has been presented by Kovalchuk, Kovev & Pinsker (1975[link]), who realized the triple-crystal arrangement on the basis of a double-crystal spectrometer by parallel mounting of the two crystals to be compared (the sample and the reference crystal) on one common axis. The advantage of this system is that Bragg angles as high as 80° are available. The device can be applied in studies of the real structure of a single crystal.

High-sensitivity (Δd/d up to [\pm3\times10^{-8}]) lattice-parameter-comparison measurement over a wide range of temperatures can be performed by means of the triple-crystal (more accurately, triple-axis) X-ray spectrometer realized by Buschert, Pace, Inzaghi & Merlini (1980[link]). The arrangement (Fig. 5.3.3.14[link] ) consists of four crystals. The first is used for obtaining a very wide but extremely parallel exit beam, which is incident on both the standard crystal S and an unknown crystal X, placed side by side on a common axis in the cryostat. The reflected beams from S and X are recorded by partially transmitting detectors [DA_2] and [DB_2], so that the beams reflect from the third crystal and are detected by the counters [DA_3] and [DB_3]. There is a small, sensitive, angle adjustment to rotate the crystal X with respect to the standard S and it is used to bring the peaks of S and X into approximate coincidence. The angular difference in the peak positions on the third axis is used for determination of lattice-parameter changes from (5.3.3.46)[link], so that [\Delta\theta=\Delta\theta_3/2-\Delta\theta_2, \eqno (5.3.3.47)]where [\Delta\theta_2] and [\Delta\theta_3] are the differences in peak positions at axes (2) and (3), respectively. The device was used, for example, to study the effect of isotope concentration on the lattice parameter of germanium perfect crystals (Buschert, Merlini, Pace, Rodriguez & Grimsditch, 1988[link]). The measured differences in the lattice parameter, of the order of 1 part in 105, were compared with those evaluated theoretically, and a very good agreement was obtained.

[Figure 5.3.3.14]

Figure 5.3.3.14 | top | pdf |

Schematic representation of the double-beam triple-crystal spectrometer of Buschert et al. (1980[link]).

Another variant of a multiple-beam arrangement, based on a triple-crystal spectrometer, was proposed by Kubena & Holý (1988[link]). The authors compared the distances of lattice planes in a direction perpendicular to the surface of the sample while studying the growth striations. One well collimated and monochromated beam coming from the first crystal was directed into the sample, and then two beams – one transmitted and one diffracted in the sample – diffracted in the reference crystal. Intensities of the diffracted beams were measured by two detectors. The difference of lattice spacings of the sample and the reference crystal was determined from the difference in positions of respective peaks. The accuracy of the lattice-spacing comparison of 2 parts in 107 and the precision of 1 part in 107 were obtained.

A four-crystal six-reflection diffractometer (Fewster, 1989[link]) was built to study crystals distorted by epitaxy and defects in nearly perfect crystals. Fig. 5.3.3.15[link] is a schematic diagram of this device. The two-crystal four-reflection Bartels monochromator (Bartels, 1983[link]) defines a narrow reflectivity profile. The analyser selects the angular range diffracted from the sample. The device may be used for recording both near-perfect rocking curves from distorted crystals (when rotations of the sample and the analyser are coupled) and a diffraction-space map for studying the diffuse scattering (when the two rotations are uncoupled). Various applications of such high-sensitivity multiple-crystal X-ray spectrometers for reciprocal-space mapping and imaging (topography), which are outside the scope of the present paper, are reviewed by Fewster (1993[link], and references therein).

[Figure 5.3.3.15]

Figure 5.3.3.15 | top | pdf |

The geometry of the diffractometer used by Fewster & Andrew (1995[link]). The scattering angle, 2ω′, is the fundamental angle for determination of the interplanar spacing and P is the analyser-groove entrance.

As was shown a few years later by Fewster & Andrew (1995[link]), the device can also be used for absolute lattice-parameter measurements of single-crystal and polycrystalline materials with a relative accuracy of a few parts in [10^6]. The authors checked the angular resolution and the sample centring of their instrument, and discussed systematic errors due to refraction, the Lorenz and polarization factor, the diffracting-plane tilt and the peak-position determination.

5.3.3.8. Optical and X-ray interferometry – a non-dispersive technique

| top | pdf |

The accuracy of an absolute measurement can be improved, in relation to that obtained in traditional methods (cf. Subsection 5.3.3.5[link]), either if the wavelength of the radiation used in an experiment is known with better accuracy [cf. equation (5.3.1.3)[link]] or if a high-quality standard single crystal is given, whose lattice spacing has been very accurately determined (Baker & Hart, 1975[link]; mentioned in §5.3.3.7.3[link]). The two tasks, i.e. very accurate determination of both lattice spacings and wavelengths in metric units, can be realized by use of combined optical and X-ray interferometry. This original concept of absolute-lattice-spacing determination directly in units of a standard light wavelength has been proposed and realized by Deslattes (1969[link]) and Deslattes & Henins (1973[link]).

The principle of the method is presented in Fig. 5.3.3.16[link] . The silicon-crystal X-ray interferometer is a symmetric Laue-case type (Bonse & te Kaat, 1968[link]). The parallel translation device consists of the stationary assembly (a) formed by two specially prepared crystals, and a moveable one (b), to which belongs the third crystal. One of the two mirrors of a high-resolution Fabry–Perot interferometer is attached to the stationary assembly and the second to the moving assembly. A stabilized He–Ne laser is used as a source of radiation, the wavelength of which has been established relative to visible standards. The first two crystals produce a standing wavefield, which is intercepted by the third crystal, so that displacement of the third crystal parallel to the diffraction vector (as suggested by the large arrow) produces alternate maxima and minima in the diffracted beams, detected by X-ray detector (c). Resonant transmission maxima of the optical interferometer are detected simultaneously by the photomultiplier indicated at (d). Analysis of the fringes (shown in Fig. 5.3.3.17[link] ) is the basis for the calculation of the lattice-spacing-to-optical-wavelength ratio (d/λ), which is given by [{2d\over\lambda} = {n\over m}\,{\cos\alpha \over \cos\beta}, \eqno (5.3.3.48)]where n and m are the numbers of optical and X-ray diffraction fringes, respectively, and α and β are the measured angular deviations of the optical and X-ray diffraction vectors from the direction of motion. The measurements are carried out in two steps. First, the lattice parameter of silicon along the [110] crystallographic direction was measured in the metric system, independently of the X-ray wavelength used in the experiment. As the next step, a specimen of known lattice spacing, treated as a reference crystal, was used for the accurate wavelength determination of Cu Kα1 and Mo Kα1. Accuracy better than 1 part in 106 was reported (see Section 4.2.2[link] ).

[Figure 5.3.3.16]

Figure 5.3.3.16 | top | pdf |

Optical and X-ray interferometry. Schematic representation of the experimental set-up (after Deslattes & Henins, 1973[link]; Becker et al., 1981[link]).

[Figure 5.3.3.17]

Figure 5.3.3.17 | top | pdf |

Portion of a dual-channel recording of X-ray and optical fringes (Deslattes, 1969[link]).

The above experiment was a turning point in accurate measurements of both wavelengths and lattice parameters. Owing to the idea of Deslattes & Henins, it became possible to determine the wavelength in nanometres rather than in troublesome XU or Å* units (cf. §4.2.1.1.1[link] ). However, the results obtained and the method itself needed verification and some adjustments. These were performed by another group of experimenters with a similar but different measuring device (Becker, Seyfried & Siegert, 1982[link], and references therein; Siegert, Becker & Seyfried, 1984[link]).

5.3.3.9. Lattice-parameter and wavelength standards

| top | pdf |

An extended series of measurements performed by means of the optical and X-ray interferometry (cf. §5.3.3.8[link]) led, among other things, to evaluation of the lattice spacing of a highly perfect silicon sample WASO 4.2.A (Becker et al., 1981[link]). Such silicon samples may be used as reference crystals in successive lattice-spacing comparison measurements – with a double-source double-crystal spectrometer (Windisch & Becker, 1990[link]), for example. The latter measurements provided new excellent lattice-spacing standards (WASO 9, for example) of the well known lattice-parameter values. As shown by the authors, the differences in lattice parameters of different samples of float-zone silicon (due to oxygen or carbon content) were not greater than a few parts in [10^{8}]. Finally, the lattice parameter of silicon, a = 5.43102088 (16) Å, has been accepted as the atomic scale length standard (Mohr & Taylor, 2000[link]).

Another reference material reported is crystals of pure rhombohedral corrundum (α-Al2O3), i.e. of ruby or sapphire (Herbstein, 2000[link], and references therein; Shvyd'ko et al., 2002[link]).

With silicon standards, measurements or remeasurements of [K\alpha _{1,2}] and/or [K\beta _{1,3}] X-ray emission lines and absolute wavelength determinations of most of the 3d transition metals (Cr, Mn, Fe, Co, Ni and Cu) have been performed [Härtwig, Grosswig, Becker & Windisch, 1991[link]; Hölzer, Fritsch, Deutsch, Härtwig & Förster, 1997[link] (see §4.2.2[link] )].

The standard crystals may also be used for determination of such physical quantities as the Avogadro constant (Deslattes et al., 1974[link]; Deslattes, Henins, Schoonover, Caroll & Bowman, 1976[link]). The single accurate wavelength values, on the other hand, may be used both in simple measurements of lattice parameters [based directly on the Bragg law, equation (5.3.1.1)[link]] and for accurate scaling of the wavelength spectra, in order to use them, for example, in high-accuracy lattice-parameter measurements based on complete convolution models [cf. §5.3.3.3.1, point(ii)[link]].

Unlike the X-rays emitted from an X-ray tube, for which the spectral line and the characteristic wavelength are known, there are no such characteristic features in synchrotron radiation. Therefore, special energy-selective monochromators should be applied in relative lattice-spacing measurements using synchrotron radiation. Obaidur (2002[link]) proposes two measurement schemes, using two types of high-resolution channel-cut monolithic monochromators. The first scheme (see Fig. 5.3.3.18[link] ) is a modification of the Bond method. The second one (see Fig. 5.3.3.19[link] ) uses the simultaneous Bragg condition for the indices (5, 1, 3), [(\bar5,\bar1,\bar3)], (1, 5, 3) and [(\bar1,\bar5,\bar3)]. The lattice-spacing differences in Si wafers were determined in the sub-parts in [10^6] range of 0.6 parts in [10^6] (in the first scheme) and of 0.2 parts in [10^6] (second scheme).

[Figure 5.3.3.18]

Figure 5.3.3.18 | top | pdf |

Synchrotron radiation, SR, from the bending magnet incident on the Si(111) double-crystal monochromator and, after four reflections from the monolithic monochromator (0.1410 nm), impinges on sample Si(444). Two diffractions are recorded at the photodiode detectors, PIN1 and PIN2. The ω1 and ω2 values of the crystal positions are recorded using a Heiden height encoder.

[Figure 5.3.3.19]

Figure 5.3.3.19 | top | pdf |

Synchrotron radiation, SR, from the bending magnet incident on the Si(111) double-crystal monochromator and, after four reflections from the monolithic monochromator (0.1612 nm), impinges on sample Si(153). Two diffractions are recorded at the photodiode detectors, PIN1 and PIN2. The differences between the two peaks' D values are recorded using a Heiden height encoder.

Recently, a new atomic scale wavelength standard was proposed by Shvyd'ko et al. (2000[link]), instead of the wavelength of the Cu [K\alpha _{1}] emission line or of the lattice parameter of a silicon standard. It is the wavelength, λM, of the 57Fe Mössbauer radiation, i.e. of γ radiation of natural linewidth from nuclear transitions. It has been measured to the sub-parts in [10^6] accuracy: λM = 0.86025474 (16) Å (relative accuracy 0.19 parts in [10^6]). Its advantage, in relation to the previous standards, is the high spectral sharpness of the Mössbauer radiation of [3.5\times 10^{-13}] in relative units, which makes its wavelength λM extremely well defined. This standard wavelength value, which lies a little outside of scope of the present review (X-ray methods), was next used for the lattice-parameter determination of sapphire single crystals with a relative accuracy of about 0.5 parts in [10^6] (Shvyd'ko et al., 2002[link]). Fig. 5.3.3.20[link] is a diagram of the measurement arrangement.

[Figure 5.3.3.20]

Figure 5.3.3.20 | top | pdf |

Experimental set-up for measuring lattice parameters. X-rays after a high-heat-load monochromator (not shown) pass through the vertical slits S1 and S2 at a distance of 26.6 m. λ is a `λ-meter'; F is a 57Fe foil used as a source of the Mössbauer radiation of high brightness; D is a semi-transparent avalanche photodiode with 0.7 ns time resolution; Al2O3 is a sapphire single crystal in a furnace on a four-circle goniometer.

References

Abrahams, S. C., Liminga, R., Marsh, P., Schrey, F., Albertsson, J., Svensson, C. & Kvick, Å. (1983). Thermal expansivity of α-LiIO3 between 20 and 520 K. J. Appl. Cryst. 16, 453–457.
d'Amour, H., Denner, W., Schulz, H. & Cardona, M. (1982). A uniaxial stress apparatus for single-crystal X-ray diffraction on a four-circle diffractometer: Application to silicon and diamond. J. Appl. Cryst. 15, 148–153.
Ando, M., Bailey, D. & Hart, M. (1978). A simple Bragg-spacing comparator. Acta Cryst. A34, 484–489.
Ando, M., Hagashi, Y., Usuda, K., Yasuami, S. & Kawata, H. (1989). A precision Bond method with SR. Rev. Sci. Instrum. 60, 2410–2413.
Arndt, U. W. & Willis, B. T. M. (1966). Single crystal diffractometry. Cambridge University Press.
Åsbrink, S., Wołcyrz, M. & Hong, S.-H. (1985a). X-ray Bond-type diffractometric investigations on V3O5 in the temperature interval 298 to 480 K including the phase transition temperature Tt = 428 K. Phys. Status Solidi A, 87, 135–140.
Åsbrink, S., Wołcyrz, M. & Hong, S.-H. (1985b). X-ray Bond-type diffractometric investigations on V3O5 in the temperature interval 298 to 480 K including the phase transition temperature Tt = 428 K. Erratum. Phys. Status Solidi A, 89, 415.
Azároff, L. V. (1974). X-ray spectroscopy, Chap. 2. New York: McGraw-Hill.
Bačkovský, J. (1965). On the most accurate measurements of the wavelengths of X-ray spectral lines. Czech. J. Phys. B15, 752–759.
Baker, J. A., Tucker, T. N., Moyer, N. E. & Buschert, R. C. (1968). Effects of carbon on the lattice parameter of silicon. J. Appl. Phys. 39, 4365–4368.
Baker, J. F. C. & Hart, M. (1975). An absolute measurement of the lattice parameter of germanium using multiple-beam X-ray diffractometry. Acta Cryst. A31, 364–367.
Baker, J. F. C., Hart, M., Halliwell, M. A. G. & Heckingbottom, R. (1976). Precise lattice parameter determination of dislocation-free gallium arsenide. I. X-ray measurements. Solid-State Electron. 19, 331–334.
Baker, T. W., George, J. D., Bellamy, B. A. & Causer, R. (1966). Very high precision X-ray diffraction. Nature (London), 210, 720–721.
Baker, T. W., George, J. D., Bellamy, B. A. & Causer, R. (1968). Fully automated high-precision X-ray diffraction. Adv. X-ray Anal. 11, 359–375.
Barla, K., Herino, R., Bomchil, G. & Pfister, J. C. (1984). Determination of lattice parameter and elastic properties of porous silicon by X-ray diffraction. J. Cryst. Growth, 68, 727–732.
Barns, R. L. (1972). A strategy for rapid and accurate (p.p.m.) measurement of lattice parameters of single crystals by Bond's method. Adv. X-ray Anal. 15, 330–338.
Bartels, W. J. (1983). Characterization of thin layers on perfect crystals with a multipurpose high-resolution X-ray diffractometer. J. Vac. Sci. Technol. B1, 338–345.
Bearden, J. A. (1933). The wavelengths of the silver, molybdenum, copper, iron and chromium Kα1 lines. Phys. Rev. 43, 92–97.
Bearden, J. A. (1965). Selection of W Kα1 as the X-ray wavelength standard. Phys. Rev. 137, BY55–BY61.
Bearden, J. A. (1967). X-ray wavelengths. Rev. Mod. Phys. 39, 78–124.
Bearden, J. A. & Henins, A. (1965). Precision measurement of lattice imperfections with a photographic two-crystal method. Rev. Sci. Instrum. 36, 334–338.
Bearden, J. A. & Thomsen, J. S. (1971). The double-crystal X-ray spectrometer: corrections, errors, and alignment procedure. J. Appl. Cryst. 4, 130–138.
Becker, P., Dorenwendt, K., Ebeling, G., Lauer, R., Lucas, W., Probst, R., Rademacher, H.-J., Reim, G., Seyfried, P. & Siegert, H. (1981). Absolute measurement of the (220) lattice plane spacing in silicon crystal. Phys. Rev. Lett. 46, 1540–1543.
Becker, P., Seyfried, P. & Siegert, H. (1982). The lattice parameter of highly pure silicon single crystals. Z. Phys. B, 48, 17–21.
Berger, H. (1984). A method for precision lattice-parameter measurement of single crystals. J. Appl. Cryst. 17, 451–455.
Berger, H. (1986a). Systematic errors in precision lattice parameter determination of single crystals caused by asymmetric line profiles. J. Appl. Cryst. 19, 34–38.
Berger, H. (1986b). Study of the Kα emission spectrum of copper. X-ray Spectrom. 15, 241–243.
Berger, H. (1993). X-ray diffraction studies on point defects in II–VI compounds. Cryst. Res. Technol. 28, 795–801.
Berger, H., Lehmann, A. & Schenk, M. (1985). Lattice parameter variations in PbTe single crystals. Cryst. Res. Technol. 20, 579–581.
Berger, H., Rosner, B. & Schikora, D. (1989). Lattice parameter determination of superlattices. Cryst. Res. Technol. 24, 437–441.
Beu, K. E. (1967). The precise and accurate determination of lattice parameters. Handbook of X-rays, edited by E. F. Kaelble, Chap. 10. New York: McGraw-Hill.
Beu, K. E., Musil, F. J. & Whitney, D. R. (1962). Precise and accurate lattice parameters by film powder methods. I. The likelihood ratio method. Acta Cryst. 15, 1292–1301.
Bolotina, N. B. (1989). Refinement of unit-cell parameters and orientation of specimen in diffractometer, taking account of symmetry of single crystal. Kristallografiya, 34, 598–601. (English transl: Sov. Phys. Crystallogr. 34, 355–357.)
Bond, W. L. (1960). Precision lattice constant determination. Acta Cryst. 13, 814–818.
Bond, W. L. (1975). Precision lattice constant determination. Erratum. Acta Cryst. A31, 698.
Bonse, U. & te Kaat, E. (1968). A two-crystal X-ray interferometer. Z. Phys. 214, 16–21.
Bowen, D. K. & Tanner, B. K. (1995). A method for the accurate comparison of lattice parameters. J. Appl. Cryst. 28, 753–760.
Bragg, W. H. & Bragg, W. L. (1915). X-rays and crystal structure, Chap. 2. London: G. Bell and Sons.
Burke, J. & Tomkeieff, M. V. (1968). Specimen and beam tilt errors in Bond's method of lattice parameter determination. Acta Cryst. A24, 683–685.
Burke, J. & Tomkeieff, M. V. (1969). Errors in the Bond method of lattice parameter determinations. Further considerations. J. Appl. Cryst. 2, 247–248.
Buschert, R. C. (1965). X-ray lattice parameter and linewidth studies in silicon. Bull. Am. Phys. Soc. 10, 125.
Buschert, R. C., Merlini, A. E., Pace, S., Rodriguez, S. & Grimsditch, M. H. (1988). Effect of isotope concentration on the lattice parameter of germanium perfect crystals. Phys. Rev. B, 38, 5219–5221.
Buschert, R. C., Meyer, A. J., Stuckey Kauffman, D. & Gotwals, J. K. (1983). A double-source double-crystal X-ray spectrometer for high-sensitivity lattice-parameter difference measurements. J. Appl. Cryst. 16, 599–605.
Buschert, R. C., Pace, S., Inzaghi, D. & Merlini, A. E. (1980). A high-sensitivity double-beam triple-crystal spectrometer for lattice parameter and topographic measurements. J. Appl. Cryst. 13, 207–210.
Busing, W. R. & Levy, H. A. (1967). Angle calculations for 3- and 4-circle X-ray and neutron diffractometers. Acta Cryst. 22, 457–464.
Cembali, F., Fabri, R., Servidori, M., Zani, A., Basile, G., Cavagnero, G., Bergamin, A. & Zosi, G. (1992). Precise X-ray relative measurement of lattice parameters of silicon wafers by multiple-crystal Bragg-case diffractometry. Computer simulation of the experiment. J. Appl. Cryst. 25, 424–431.
Chang, S.-L. (1984). Multiple diffraction of X-rays in crystals, Chap. 7.2 in particular. Berlin: Springer-Verlag.
Clegg, W. (1984). Orientation matrix refinement during four-circle diffractometer data collection. Acta Cryst. A40, 703–704.
Clegg, W. & Sheldrick, G. M. (1984). The refinement of unit cell parameters from two-circle diffractometer measurements. Z. Kristallogr. 167, 23–27.
Cole, H., Chambers, F. W. & Dunn, H. M. (1962). Simultaneous diffraction: indexing Umweganregung peaks in simple cases. Acta Cryst. 15, 138–144.
Compton, A. H. & Allison, S. K. (1935). X-rays in theory and experiment. New York: Van Nostrand.
Cooper, A. S. (1962). Precise lattice constants of germanium, aluminium, gallium arsenide, uranium, sulphur, quartz and sapphire. Acta Cryst. 15, 578–582.
Deslattes, R. D. (1969). Optical and X-ray interferometry of a silicon lattice spacing. Appl. Phys. Lett. 15, 386–388.
Deslattes, R. D. & Henins, A. (1973). X-ray to visible wavelength ratios. Phys. Rev. Lett. 31, 972–975.
Deslattes, R. D., Henins, A., Bowman, H. A., Schoonover, R. M., Caroll, C. L., Barnes, I. L., Machlan, L. A., Moore, L. J. & Shields, W. R. (1974). Determination of the Avogadro constant. Phys. Rev. Lett. 33, 463–466.
Deslattes, R. D., Henins, A., Schoonover, R. M., Caroll, C. L. & Bowman, H. A. (1976). Avogadro constant – correction to an earlier report. Phys. Rev. Lett. 36, 898–900.
Dressler, L., Griebner, U. & Kittner, R. (1987). Precision measurement of lattice parameters in LiF monocrystals. Cryst. Res. Technol. 22, 1431–1435.
Eastabrook, J. N. (1952). Effect of vertical divergence on the displacement and breadth of X-ray powder diffraction lines. Br. J. Appl. Phys. 3, 349–352.
Fewster, P. F. (1982). Absolute lattice-parameter measurements of epitaxial layers. J. Appl. Cryst. 15, 275–278.
Fewster, P. F. (1989). A high-resolution multiple-crystal multiple-reflection diffractometer. J. Appl. Cryst. 22, 64–69.
Fewster, P. F. (1993). Structural characterisation of materials by combining X-ray diffraction space mapping and topography. Philips J. Res. 47, 235–245.
Fewster, P. F. & Andrew, N. L. (1995). Absolute lattice-parameter measurement. J. Appl. Cryst. 28, 451–458.
Fewster, P. F. & Willoughby, A. F. W. (1980). The effect of silicon doping on the lattice parameter of gallium arsenide grown by liquid-phase epitaxy, vapour-phase epitaxy and gradient-freeze techniques. J. Cryst. Growth, 50, 648–653.
Filscher, G. & Unangst, D. (1980). Bond-method for precision lattice constant determination. Dependence of lattice constant error on sample adjustment and collimator tilt. Krist. Tech. 15, 955–960.
Fukahara, A. & Takano, Y. (1977). Determination of strain distributions from X-ray Bragg reflexion by silicon single crystals. Acta Cryst. A33, 137–142.
Fukumori, T. & Futagami, K. (1988). Measurements of lattice parameters and half-widths of the rocking curve on GaAs crystal by the X-ray double-crystal method using a Cu Kα doublet. Jpn. J. Appl. Phys. 27, 442–443.
Fukumori, T., Futagami, K. & Matsunaga, K. (1982). X-ray double-crystal method for crystal lattice parameter measurements using Cu Kα doublet. Jpn. J. Appl. Phys. 21, 1525.
Fukumori, T., Imai, K., Hasegawa, T. & Akashi, Y. (1997). Precision lattice spacing measurement using X-ray Cu Kα doublet. J. Phys. Soc. Jpn, 66, 1976–1978.
Gabe, E. J. (1980). Diffractometer control with minicomputers. Computing in crystallography, edited by R. Diamond, S. Ramaseshan & K. Venkatesan, pp. 1.01–1.18. Bangalore: Indian Academy of Sciences.
Gałdecka, E. (1985). The variances and covariances of measured intensities in precise lattice-constant determination by the Bond method. Structure & statistics in crystallography, edited by A. J. C. Wilson, pp. 137–149. New York: Adenine Press.
Gałdecka, E. (1993a). Description and peak-position determination of a single X-ray diffraction profile for high-accuracy lattice-parameter measurements by the Bond method. I. An analysis of descriptions available. Acta Cryst. A49, 106–115.
Gałdecka, E. (1993b). Description and peak-position determination of a single X-ray diffraction profile for high-accuracy lattice-parameter measurements by the Bond method. II. Testing and choice of description. Acta Cryst. A49, 116–126.
Gałdecka, E. (1994). The extrapolated-peak method for the peak-position determination of an X-ray diffraction profile, and the accuracy of the Bragg-angle measurements. Sci. Bull. Tech. Univ. Łódź Branch Bielsko-Biała, 22(4), 1–24.
Godwod, K., Kowalczyk, R. & Szmid, Z. (1974). Application of a precise double X-ray spectrometer for accurate lattice parameter determination. Phys. Status Solidi A, 21, 227–234.
Golovin, A. L., Imamov, R. M. & Kondrashkina, E. A. (1985). Absolute measurements of lattice spacings in surface layers of crystals. Phys. Status Solidi A, 89, K5–K7.
Grosswig, S., Härtwig, J., Alter, U. & Christoph, A. (1983). Precision lattice parameter determination of coloured quartz monocrystals. Cryst. Res. Technol. 18, 501–511.
Grosswig, S., Härtwig, J., Jäckel, K.-H., Kittner, R. & Melle, W. (1986). A novel arrangement for the absolute measurement of geometric lattice parameters of monocrystals with high precision. Nauch. Apparat. Sci. Instrum. 1, 29–44.
Grosswig, S., Jäckel, K.-H. & Kittner, R. (1986). Peak position determination of X-ray diffraction profiles in precision lattice parameter measurements according to the Bond-method with help of the polynomial approximation. Cryst. Res. Technol. 21, 133–139.
Grosswig, S., Jäckel, K.-H., Kittner, R., Dietrich, B. & Schellenberger, U. (1985). Determination of the coplanar geometric lattice parameters of monocrystals with high precision. Cryst. Res. Technol. 20, 1093–1100.
Grosswig, S., Melle, W., Schellenberger, U. & Zahorowski, W. (1983). High precision lattice parameter determination of KDP with different crystal perfection. Cryst. Res. Technol. 18, K28–K30.
Gruber, E. E. & Black, R. E. (1970). Analysis of the axial misalignment error in precision lattice parameter measurement by the Bond technique. J. Appl. Cryst. 3, 354–357.
Halliwell, M. A. G. (1970). Measurement of specimen tilt and beam tilt in the Bond method. J. Appl. Cryst. 3, 418–419.
Hamilton, W. C. (1974). Angle settings for four-circle diffractometers. International tables for X-ray crystallography, Vol. IV, pp. 274–284. Birmingham: Kynoch Press.
Hart, M. (1969). High precision lattice parameter measurements by multiple Bragg reflexion diffractometry. Proc. R. Soc. London Ser. A, 309, 281–296.
Hart, M. (1981). Bragg angle measurement and mapping. J. Cryst. Growth, 55, 409–427.
Hart, M. & Lloyd, K. H. (1975). Measurement of strain and lattice parameter in epitaxic layers. J. Appl. Cryst. 8, 42–44.
Härtwig, J., Bąk-Misiuk, J., Berger, H., Brühl, H.-G., Okada, Y., Grosswig, S., Wokulska, K. & Wolf, J. (1994). Comparison of lattice parameters obtained from an internal monocrystal standard. Phys. Status Solidi A, 142, 19–26.
Härtwig, J. & Grosswig, S. (1989). Measurement of X-ray diffraction angles of perfect monocrystals with high accuracy using a single crystal diffractometer. Phys. Status Solidi A, 115, 369–382.
Härtwig, J., Grosswig, S., Becker, P. & Windisch, D. (1991). Remeasurement of the Cu Kα1 emission X-ray wavelength in the metrical system (present stage). Phys. Status Solidi A, 125, 79–89.
Härtwig, J., Hölzer, G., Förster, E., Goetz, K., Wokulska, K. & Wolf, J. (1994). Remeasurement of characteristic X-ray emission lines and their application to line profile analysis and lattice parameter determination. Phys. Status Solidi A, 143, 23–34.
Härtwig, J., Hölzer, G., Wolf, J. & Förster, E. (1993). Remeasurement of the profile of the characteristic Cu Kα emission line with high precision and accuracy. J. Appl. Cryst. 26, 539–548.
Häusermann, D. & Hart, M. (1990). A fast high-accuracy lattice-parameter comparator. J. Appl. Cryst. 23, 63–69.
Herbstein, F. H. (2000). How precise are measurements of unit-cell dimensions from single crystals? Acta Cryst. B56, 547–557.
Holý, V. & Härtwig, J. (1988). The role of diffuse scattering on microdefects in the precise lattice parameter measurement. Phys. Status Solidi B, 145, 363–372.
Hölzer, G., Fritsch, M., Deutsch, M., Härtwig, J. & Förster, E. (1997). 1,2 and Kβ1,3 X-ray emission lines of the 3d transition metals. Phys. Rev. A, 56, 4554–4568.
Hom, T., Kiszenick, W. & Post, B. (1975). Accurate lattice constants from multiple reflection measurements. II. Lattice constants of germanium, silicon and diamond. J. Appl. Cryst. 8, 457–458.
Horváth, J. & Kucharczyk, D. (1981). Temperature dependence of lattice parameters of PbHPO4 and PbDPO4 single crystals. Phys. Status Solidi A, 63, 687–692.
Hubbard, C. R. & Mauer, F. A. (1976). Precision and accuracy of the Bond method as applied to small spherical crystals. J. Appl. Cryst. 9, 1–8.
Irie, K., Koshiji, N. & Okazaki, A. (1989). High-angle double-crystal X-ray diffractometry (HADOX): combination with a sealed-tube X-ray source. Jpn. J. Appl. Phys. 28, 1504–1506.
James, R. W. (1967). The optical principles of the diffraction of X-rays. London: Bell.
Keller, H. L., Kucharczyk, D. & Küppers, H. (1982). The ferroelastic monoclinic low temperature modification of ammonium hydrogen oxalate hemihydrate. Z. Kristallogr. 158, 221–232.
Kheiker, D. M. (1973). Rentgenowskaya diffraktometriya monokristallow, Chaps. 3, 4, 5. Leningrad: Mashinostroyenie.
Kheiker, D. M. & Zevin, L. S. (1963). Rentgenowskaya diffraktometriya (X-ray diffractometry), Chap. 4. Moscow: Fizmatgiz.
Kirk, D. & Caulfield, P. B. (1977). Location of diffractometer profiles in X-ray stress analysis. Adv. X-ray Anal. 20, 283–289.
Kishino, S. (1973). Improved techniques of lattice parameter measurements using two X-ray beams. Adv. X-ray Anal. 16, 367–378.
Kobayashi, J., Mizutani, I. & Schmidt, H. (1970). X-ray study on the lattice strains of ferroelectric iron iodine boracite Fe3B7O13. Phys. Rev. B, 1, 3801–3808.
Kobayashi, J., Yamada, N. & Nakamura, T. (1963). Origin of the visibility of the antiparallel 180° domains in barium titanate. Phys. Rev. Lett. 11, 410–414.
Korytár, D. (1984). Lateral lattice parameter variation measurement by means of a double crystal X-ray method with oscillating slit. Czech. J. Phys. B34, 1277–1281.
Kossel, W. (1936). Messungen am vollständigen Reflexsystem eines Kristallgitters. Ann. Phys. (Leipzig), 26, 533–553.
Kovalchuk, M. V., Kovev, E. K. & Pinsker, Z. G. (1975). The X-ray triple-crystal spectrometer and precision determination of Δdhkl. (In Russian.) Kristallografiya, 20, 142–148.
Kshevetsky, S. A., Mikhailyuk, I. P., Ostapovich, M. V., Polyak, M. I., Remenyuk, P. I. & Fomin, V. G. (1979). Application of multiple diffraction to determination of lattice parameters. (In Russian.) Ukr. Fiz. Zh. 24, 1480–1485.
Kshevetsky, S. A., Mikhalchenko, V. P., Stetsko, Yu. P. & Shelud'ko, S. A. (1985). The lattice parameter refinement of single crystals by means of multiple-wave diffractometry. (In Russian.) Ukr. Fiz. Zh. 30, 1843–1848.
Kubena, J. & Holý, V. (1988). Precise relative X-ray measurement of the lattice parameter of silicon crystals with growth striations. J. Appl. Cryst. 21, 245–251.
Kucharczyk, D. & Niklewski, T. (1979). Accurate X-ray determination of the lattice parameters and the thermal expansion coefficients of VO2 near the transition temperature. J. Appl. Cryst. 12, 370–373.
Kucharczyk, D., Pietraszko, A. & Łukaszewicz, K. (1976). Temperature dependence of lattice parameters of NaNO2 single crystals. Phys. Status Solidi A, 37, 287–294.
Kucharczyk, D., Pietraszko, A. & Łukaszewicz, K. (1993). An automatic four-circle diffractometer designed for precise lattice-parameter determination. J. Appl. Cryst. 26, 467.
Kudo, S. (1982). X-ray determination of incommensurate superlattices in K2SeO4 and (NH4)2BeF4. Jpn. J. Appl. Phys. 21, 255–258.
Kurbatov, B. A., Zubenko, V. V. & Umansky, M. M. (1972). The use of the monochromator crystal with anomalous transmission of X-rays in precise lattice parameter measurements. (In Russian.) Kristallografiya, 17, 1058–1060.
Larson, B. C. (1974). High-precision measurements of lattice parameter changes in neutron-irradiated copper. J. Appl. Phys. 45, 514–518.
Leszczyński, M., Podlasin, S. & Suski, T. (1993). A 109 Pa high-pressure cell for X-ray and optical measurements. J. Appl. Cryst. 26, 1–4.
Lisoivan, V. I. (1974). Local determination of all the lattice parameters of single crystals. (In Russian.) Appar. Methody Rentgenovskogo Anal. 14, 151–157.
Lisoivan, V. I. (1981). Experimental refinement of the angles between unit-cell axes. (In Russian.) Kristallografiya, 26, 458–463.
Lisoivan, V. I. (1982). Measurements of unit-cell parameters on one-crystal spectrometer. (In Russian.) Novosibirsk: Nauka.
Lisoivan, V. I. & Dikovskaya, R. R. (1969). Local precision determination of lattice constants of a single crystal. Prib. Tech. Eksp. No. 4, pp. 164–166; English transl: Instrum. Exp. Tech. (USSR), 4, 992–994.
Lonsdale, K. (1947). Divergent-beam X-ray photography of crystals. Proc. R. Soc. London Ser. A, 240, 219–250.
Luger, P. (1980). Modern X-ray analysis of single crystals. In particular, Chap. 4 and Section 4.2.2. Berlin: de Gruyter.
Łukaszewicz, K., Kucharczyk, D., Malinowski, M. & Pietraszko, A. (1978). New model of the Bond diffractometer for precise determination of lattice parameters and thermal expansion of single crystals. Krist. Tech. 13, 561–567.
Łukaszewicz, K., Pietraszko, A., Kucharczyk, D., Malinowski, M., Stępień-Damm, J. & Urbanowicz, E. (1976). Precyzyjne pomiary stałych sieciowych kryształów metoda Bonda (Precision measurements of lattice constants of crystals by the Bond method). Wrocław: Instytut Niskich Temperatur i Badań Strukturalnych PAN.
Mauer, F. A., Hubbard, C. R., Piermarini, G. J. & Block, S. (1975). Measurement of anisotropic compressibilities by a single crystal diffractometer method. Adv. X-ray Anal. 18, 437–453.
Mohr, P. J. & Taylor, B. N. (2000). CODATA recommended values of the fundamental physical constants. Rev. Mod. Phys. 72, 351–495.
Nemiroff, M. (1982). Precise lattice-constant determinations using measured beam and crystal tilts. J. Appl. Cryst. 15, 375–377.
Obaidur, R. M. (2002). Energy-selective (+,+) monolithic monochromator and relative lattice-spacing measurement of Si wafers with synchrotron radiation. J. Synchrotron Rad. 9, 28–35.
Ohama, N., Sakashita, H. & Okazaki, A. (1979). Improvement of high-angle double-crystal X-ray diffractometry (HADOX) for measuring temperature dependence of lattice constants. II. Practice. J. Appl. Cryst. 12, 455–459.
Okada, Y. (1982). A high-temperature attachment for precise measurement of lattice parameters by Bond's method between room temperature and 1500 K. J. Phys. E, 15, 1060–1063.
Okazaki, A. & Kawaminami, M. (1973a). Accurate measurement of lattice constants in a wide range of temperature: use of white X-ray and double-crystal diffractometry. Jpn. J. Appl. Phys. 12, 783–789.
Okazaki, A. & Kawaminami, M. (1973b). Lattice constant of strontium titanate at low temperatures. Mater. Res. Bull. 8, 545–550.
Okazaki, A. & Ohama, N. (1979). Improvement of high-angle double-crystal X-ray diffractometry (HADOX) for measuring temperature dependence of lattice constants. I. Theory. J. Appl. Cryst. 12, 450–454.
Okazaki, A. & Soejima, Y. (2001). Ultra-high-angle double-crystal X-ray diffractometry (U-HADOX) for determining a change in the lattice spacing: theory. Acta Cryst. A57, 708–712.
Pick, M. A., Bickmann, K., Pofahl, E., Zwoll, K. & Wenzl, H. (1977). A new automatic triple-crystal X-ray diffractometer for the precision measurement of intensity distribution of Bragg diffraction and Huang scattering. J. Appl. Cryst. 10, 450–457.
Pierron, E. D. & McNeely, J. B. (1969). Precise cell parameters of semiconductor crystals and their applications. Adv. X-ray Anal. 12, 343–353.
Pietraszko, A., Tomaszewski, P. E. & Łukaszewicz, K. (1981). X-ray and optical study of the phase transition in LiCsSO4. Phase Transit. 2, 131–150.
Pietraszko, A., Waśkowska, A., Olejnik, S. & Łukaszewicz, K. (1979). X-ray study of the phase transition in RbHSeO4. Phase Transit. 1, 99–106.
Pihl, C. F., Bieber, R. L. & Schwuttke, G. H. (1973). Precision lattice parameter studies of ion-implanted silicon. Phys. Status Solidi A, 17, 359–369.
Popović, S. (1971). An X-ray diffraction method for lattice parameter measurements from corresponding Kα and Kβ reflexions. J. Appl. Cryst. 4, 240–241.
Post, B. (1975). Accurate lattice constants from multiple diffraction measurements. I. Geometry, techniques and systematic errors. J. Appl. Cryst. 8, 452–456.
Renninger, M. (1937). `Umweganregung', eine bisher unbeachtete Wechselwirkungserscheinung bei Raumgitterinterferezen. Z. Phys. 106, 141–176.
Ridou, C., Rousseau, M. & Freund, A. (1977). Détermination précise des paramètres cristallins au voisinage de changement de phase cubique quadratique dans RbCaF3. J. Phys. (Paris), 38, L-359–L-363.
Rossmanith, E. (1985). UMWEG – a computer program for calculation and graphical representation of Umweganregung patterns. Z. Kristallogr. 171, 253–254.
Sasvári, J. & Zsoldos, É. (1980). Accurate lattice parameter measurements of epitaxial layers. Invited paper at International Symposium on Industrial Applications of X-ray Spectrometry and Diffractometry, Turawa, Poland, 15–18 April 1980.
Segmüller, A. (1970). Automated lattice parameter determination on single crystals. Adv. X-ray Anal. 13, 455–467.
Shvyd'ko, Yu. V., Lerche, M., Jäschke, J., Lucht, M., Gerdau, E., Gerken, M., Rüter, H. D., Wille, H.-C., Becker, P., Alp, E. E., Sturhahn, W., Sutter, J. & Toellner, T. S. (2000). γ-ray wavelength standard for atomic scales. Phys. Rev. Lett. 85, 495–498.
Shvyd'ko, Yu. V., Lucht, M., Gerdau, E., Lerche, M., Alp, E. E., Sturhahn, W., Sutter, J. & Toellner, T. S. (2002). Measuring wavelengths and lattice constants with the Mössbauer wavelength standard. J. Synchrotron Rad. 9, 17–23.
Siegert, H., Becker, P. & Seyfried, P. (1984). Determination of silicon unit-cell parameters by precision measurements of Bragg plane spacings. Z. Phys. B, 56, 273–278.
Skupov, V. D. & Uspeckaya, G. I. (1975). The combined X-ray spectrometer for deformation measurements in single crystals. (In Russian.) Prib. Tekh. Eksp. No. 2, pp. 210–213.
Smakula, A. & Kalnajs, J. (1955). Precision determination of lattice constants with a Geiger-counter X-ray diffractometer. Phys. Rev. 99, 1737–1743.
Soejima, Y., Tomonaga, N., Onitsuka, H. & Okazaki, A. (1991). Two-dimensional intensity distribution in high-angle double-crystal X-ray diffractometry (HADOX). Z. Kristallogr. 195, 161–168.
Stępień, J. A., Auleytner, J. & Łukaszewicz, K. (1972). X-ray examination of the real structure of γ-irradiated NaClO3 single crystals. Phys. Status Solidi A, 10, 631–638.
Stępień-Damm, J. A., Kucharczyk, D., Urbanowicz, E. & Łukaszewicz, K. (1975). Effect of γ-irradiation on the thermal expansion of sodium chlorate NaClO3. Bull. Acad. Pol. Sci. Ser. Sci. Chim. Geol. Geogr. 23, 985–988.
Stępień-Damm, J. A., Suski, T., Meysner, L., Hilczer, B. & Łukaszewicz, K. (1974). Effect of X-ray irradiation on the lattice constant of TGS crystal in the vicinity of phase transition. Bull. Acad. Pol. Sci. Ser. Sci. Chim. Geol. Geogr. 22, 685–688.
Stout, G. H. & Jensen, L. H. (1968). X-ray structure determination. London: Macmillan.
Straumanis, M. & Ieviņš, A. (1940). Die Präzizionsbestimmung von Gitterkonstanten nach der asymmetrischen Methode. Berlin: Springer. [Reprinted by Edwards Brothers Inc., Ann Arbor, Michigan (1948).]
Takano, Y. & Maki, M. (1972). X-ray measurement of lattice strain of oxygen diffused silicon. Acta Cryst. A28, S171.
Terminasov, Yu. S. & Tuzov, L. V. (1964). The double-diffraction of X-rays in crystals. (In Russian.) Usp. Phys. Nauk, 83, 223–258. (English transl: Sov. Phys. Usp. 7, 734.)
Thomsen, J. S. (1974). High-precision X-ray spectroscopy. X-ray spectroscopy, edited by L. V. Azároff, pp. 26–132. New York: McGraw-Hill.
Thomsen, J. S. & Yap, Y. (1968). Effect of statistical counting errors on wavelengths criteria for X-ray spectra. J. Res. Natl Bur. Stand. Sect. A, 72, 187–205.
Urbanowicz, E. (1981a). The influence of in-plane collimation on the precision and accuracy of lattice-constant determination by the Bond method. I. A mathematical model. Statistical errors. Acta Cryst. A37, 364–368.
Urbanowicz, E. (1981b). The influence of in-plane collimation on the precision and accuracy of lattice-constant determination by the Bond method. II. Verification of the mathematical model. Systematic errors. Acta Cryst. A37, 369–373.
Weisz, O., Cochran, W. & Cole, W. F. (1948). The accurate determination of cell dimensions from single-crystal X-ray photographs. Acta Cryst. 1, 83–88.
Wilson, A. J. C. (1950). Geiger-counter X-ray spectrometer – influence of size and absorption coefficient of specimen on position and shape of powder diffraction maxima. J. Sci. Instrum. 27, 321–325.
Wilson, A. J. C. (1963). Mathematical theory of X-ray powder diffractometry. Philips Technical Library. Eindhoven: Centrex Publishing Company.
Wilson, A. J. C. (1965). The location of peaks. Br. J. Appl. Phys. 16, 665–674.
Wilson, A. J. C. (1967). Statistical variance of line-profile parameters. Measures of intensity, location and dispersion. Acta Cryst. 23, 888–898.
Wilson, A. J. C. (1968). Statistical variance of line-profile parameters. Measures of intensity, location and dispersion: Corrigenda. Acta Cryst. A24, 478.
Wilson, A. J. C. (1969). Statistical variance of line-profile parameters. Measures of intensity, location and dispersion: Addendum. Acta Cryst. A25, 584–585.
Wilson, A. J. C. (1980). Accuracy in methods of lattice-parameter measurement. Natl Bur. Stand. (US) Spec. Publ. No. 567. Proceedings of Symposium on Accuracy in Powder Diffraction, NBS, Gaithersburg, MD, USA, 11–15 June 1979.
Windisch, D. & Becker, P. (1990). Silicon lattice parameters as an absolute scale of length for high precision measurements of fundamental constants. Phys. Status Solidi A, 118, 379–388.
Wołcyrz, M. & Łukaszewicz, K. (1982). The evaluation of crystal perfection by means of the asymmetric Bragg reflections. J. Appl. Cryst. 15, 406–411.
Wołcyrz, M., Pietraszko, A. & Łukaszewicz, K. (1980). The application of asymmetric Bragg reflections in the Bond method of measuring lattice parameters. J. Appl. Cryst. 13, 12–16.
Zolotoyabko, E., Sander, B., Komem, Y. & Kantor, B. (1993). Improved strain analysis in semiconductor crystals by X-ray diffractometry enhanced with ultrasound. Appl. Phys. Lett. 63, 1540–1542.








































to end of page
to top of page