International
Tables for Crystallography Volume C Mathematical, physical and chemical tables Edited by E. Prince © International Union of Crystallography 2006 |
International Tables for Crystallography (2006). Vol. C. ch. 7.4, pp. 654-656
Section 7.4.2.2. TDS correction factor for X-rays (single crystals)
B. T. M. Willisd
|
The differential cross section, representing the intensity per unit solid angle for Bragg scattering, is where N is the number of unit cells, each of volume V, and F(h) is the structure factor. H is the scattering vector, defined by
with k and k0 the wavevectors of the scattered and incident beams, respectively. (The scattering is elastic, so k = k0 = 2π /λ, where λ is the wavelength.) 2πh is the reciprocal-lattice vector and the delta function shows that the scattered intensity is restricted to the reciprocal-lattice points.
The integrated Bragg intensity is given by where the integration is over the solid angle Ω subtended by the detector at the crystal and over the time t spent in scanning the reflection. Using
with dH = H dθ, equation (7.4.2.2)
reduces to the familiar result (James, 1962
)
where
is the angular velocity of the crystal and 2θ the scattering angle.
The differential cross section for one-phonon scattering by acoustic modes of small wavevector q is [see Section 4.1.1
of IT B (2001
)]. Here, ej(q) is the polarization vector of the mode (jq) , where j is an index for labelling the acoustic branches of the dispersion relations, m is the mass of the unit cell and
is the mode energy. The delta function in (7.4.2.4)
shows that the scattering from the mode (jq) is confined to the points in reciprocal space displaced by ±q from the reciprocal-lattice point at q = 0. The acoustic modes involved are of small wavenumber, for which the dispersion relation can be written
where vj is the velocity of the elastic wave with polarization vector ej(q). Substituting (7.4.2.5)
into (7.4.2.4)
shows that the intensity from the acoustic modes varies as 1/q2, and so peaks strongly at the reciprocal-lattice points to give rise to the TDS error.
Integrating the delta function in (7.4.2.4) gives the integrated one-phonon intensity
with ρ the crystal density. The sum over the wavevectors q is determined by the range of q encompassed in the intensity scan. The density of wavevectors is uniform in reciprocal space [see Section 4.1.1
of IT B (2001
)], and so the sum can be replaced by an integral
Thus, the correction factor (E1/E0) is given by
where
The integral in (7.4.2.6)
is over the range of measurement, and the summation in (7.4.2.7)
is over the three acoustic branches. Only long-wavelength elastic waves, with a linear dispersion relation, equation (7.4.2.5)
, need be considered.
The frequencies ωj(q) and polarization vectors ej(q) of the elastic waves in equation (7.4.2.7) can be calculated from the classical theory of Voigt (1910
) [see Wooster (1962
)]. If
,
,
are the direction cosines of the polarization vector with respect to orthogonal axes x, y, z, then the velocity vj is determined from the elastic stiffness constants
by solving the following equations of motion.
Here, Akm is the km element of a 3 × 3 symmetric matrix A; if
,
,
are the direction cosines of the wavevector q with reference to x, y, z, the km element is given in terms of the elastic stiffness constants by
The four indices klmn can be reduced to two, replacing 11 by 1, 22 by 2, 33 by 3, 23 and 32 by 4, 31 and 13 by 5, and 12 and 21 by 6. The elements of A are then given explicitly by
The setting up of the matrix A is a fundamental first step in calculating the TDS correction factor. This implies a knowledge of the elastic constants, whose number ranges from three for cubic crystals to twenty one for triclinic crystals. The measurement of elastic stiffness constants is described in Section 4.1.6
of IT B (2001
).
For each direction of propagation , there are three values of
(j = 1, 2, 3), given by the eigenvalues of A. The corresponding eigenvectors of A are the polarization vectors ej(q). These polarization vectors are mutually perpendicular, but are not necessarily parallel or perpendicular to the propagation direction.
The function J(q) in equation (7.4.2.7) is related to the inverse matrix A−1 by
where
,
,
are the x, y, z components of the scattering vector H, and classical equipartition of energy is assumed [Ej(q) = kBT] . Thus A−1 determines the anisotropy of the TDS in reciprocal space, arising from the anisotropic elastic properties of the crystal.
Isodiffusion surfaces, giving the locus in reciprocal space for which the intensity J(q) is constant for elastic waves of a given wavelength, were first plotted by Jahn (1942a,b
). These surfaces are not spherical even for cubic crystals (unless
), and their shapes vary from one reciprocal-lattice point to another.
Inserting (7.4.2.8) into (7.4.2.6)
gives the TDS correction factor as
where Tmn, an element of a 3 × 3 symmetric matrix T, is defined by
Equation (7.4.2.9)
can also be written in the matrix form
with
representing the transpose of H.
The components of H relate to orthonormal axes, whereas it is more convenient to express them in terms of Miller indices hkl and the axes of the reciprocal lattice. If S is the 3 × 3 matrix that transforms the scattering vector H from orthonormal axes to reciprocal-lattice axes, then where hT = (h, k, l) . The final expression for α, from (7.4.2.11)
and (7.4.2.12)
, is
This is the basic formula for the TDS correction factor.
We have assumed that the entire one-phonon TDS under the Bragg peak contributes to the measured integrated intensity, whereas some of it is removed in the background subtraction. This portion can be calculated by taking the range of integration in (7.4.2.10) as that corresponding to the region of reciprocal space covered in the background measurement.
To evaluate T requires the integration of the function A−1 over the scanned region in reciprocal space (see Fig. 7.4.2.2).
Both the function itself and the scanned region are anisotropic about the reciprocal-lattice point, and so the TDS correction is anisotropic too, i.e. it depends on the direction of the diffraction vector as well as on
.
Computer programs for calculating the anisotropic TDS correction for crystals of any symmetry have been written by Rouse & Cooper (1969), Stevens (1974
), Merisalo & Kurittu (1978
), Helmholdt, Braam & Vos (1983
), and Sakata, Stevenson & Harada (1983
). To simplify the calculation, further approximations can be made, either by removing the anisotropy associated with A−1 or that associated with the scanned region. In the first case, the element Tmn is expressed as
where the angle brackets indicate the average value over all directions. In the second case,
where
is the radius of the sphere that replaces the anisotropic region (Fig. 7.4.2.2
) actually scanned in the experiment, and dS is a surface element of this sphere. qm can be estimated by equating the volume of the sphere to the volume swept out in the scan.
If both approximations are employed, the correction factor is isotropic and reduces to with vL representing the mean velocity of the elastic waves, averaged over all directions of propagation and of polarization.
Experimental values of α have been measured for several crystals by γ-ray diffraction of Mössbauer radiation (Krec & Steiner, 1984). In general, there is good agreement between these values and those calculated by the numerical methods, which take into account anisotropy of the TDS. The correction factors calculated analytically from (7.4.2.14)
are less satisfactory.
The principal effect of not correcting for TDS is to underestimate the values of the atomic displacement parameters. Writing , we see from (7.4.2.14)
that the overall displacement factor is increased from B to B + ΔB when the correction is made. ΔB is given by
Typically, ΔB/B is 10–20%. Smaller errors occur in other parameters, but, for accurate studies of charge densities or bonding effects, a TDS correction of all integrated intensities is advisable (Helmholdt & Vos, 1977
; Stevenson & Harada, 1983
).
References













