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. 4.2, pp. 245-248
Section 4.2.6.2.3. Relativistic theories
D. C. Creaghb
|
It is necessary to consider relativistic effects for atoms having all but the smallest atomic numbers. Cromer & Liberman (1970) produced a set of tables based on a relativistic approach to the scattering of photons by isolated atoms that was later reproduced in IT IV (1974
). Subsequent experimental determinations drew attention to inaccuracies in these tables in the neighbourhood of absorption edges owing to the poor convergence of the Gaussian integration technique, which was used to evaluate the real part of the dispersion correction. In a later paper, Cromer & Liberman (1981
) recalculated 34 instances for which the incident radiation lay close to the absorption edges of atoms using a modified integration procedure. Care should be exercised when using the Cromer & Liberman computer program, especially for calculations of
for high atomic weight elements at low photon energies. As Creagh (1990
) and Chantler (1994
) have shown, incorrect values of
can be calculated because an insufficient number of values of
are calculated prior to performing the Kramers–Kronig transform. In a new tabulation, Chantler (1995
) presents the Cromer & Liberman data using a finer integrating grid. It should be noted that the relativistic correction is the same as that used in this tabulation.
These relativistic calculations are based on the scattering formula developed by Akhiezer & Berestetsky (1957) for the scattering amplitude for photons by a bound electron, viz:
Here the angular frequencies of the incident and scattered photons are
and
, respectively, and the initial and final energy states of the atom are
and
, respectively. The scattering factor f is a complicated expression that includes the initial and final polarization states of the photon, the Dirac velocity operator, and the phase factors
and
for the incident and scattered waves, respectively. Summation is over all positive and negative intermediate states except those positive energy states occupied by other atomic electrons. The form of this expression is not easily related to the form-factor formalism that is most widely used by crystallographers, and a number of manipulations of the formula for the scattering factor are necessary to relate it more directly to the crystallographic formalism. In doing so, a number of assumptions and simplifications were made. Cromer & Liberman restricted their study to coherent, forward scatter in which changes in photon polarization did not occur. With these approximations, and using the electrical dipole approximation
, they were able to show that
In equation (4.2.6.25)
, f(0) is the atomic form factor for the case of forward scatter
, and the term
arises from the application of the dipole approximation to determine the contribution of bound electrons to the scattering process. The term
is related to the photoelectric scattering cross section expressed as a function of photon energy
by
and
These equations may be compared with equations (4.2.6.23)
and (4.2.6.22)
, respectively. But equation (4.2.6.25)
differs from equation (4.2.6.21)
by the term
, which is constant for each atomic species, and is related to the total Coulomb energy of the atom. Evidently, to keep the formalism the same, one must write
In Table 4.2.6.1
, values of
are set out as a function of atomic number for elements ranging in atomic number from 3 to 98.
|
To develop their tables, Cromer & Liberman (1970) used the Brysk & Zerby (1968
) computer code for the calculation of photoelectric cross sections, which was based on Dirac–Slater relativistic wavefunctions (Liberman, Waber & Cromer, 1965
). They employed a value for the exchange potential of 0.667
and experimental rather than computed values of the energy eigenvalues for the atoms.
The wide use of their tables by crystallographers inevitably meant that criticism of the accuracy of the tables was forthcoming on both theoretical and experimental grounds. Stibius-Jensen (1979) drew attention to the fact that the use of the dipole approximation too early in the argument caused an error of
in the tabulated values. More recently, Cromer & Liberman (1981
) include this term in their calculations. Some experimental deficiencies of the tabulated values of
have been discussed by Cusatis & Hart (1977
), Hart & Siddons (1981
), Creagh (1980
, 1984
, 1985
, 1986
), Deutsch & Hart (1982
), Dreier, Rabe, Malzfeldt & Niemann (1984
), Bonse & Hartmann-Lotsch (1984
), and Bonse & Henning (1986
).
In the latter two cases, the Kramers–Kronig transformation of photoelectric scattering results has been performed without taking into account the term that arises in the relativistic case for the total Coulomb energy of the atom. Although good agreement with the Cromer & Liberman tables is claimed, their failure to include this term is an implied criticism of the Cromer & Liberman tables. That this is unjustified can be seen by references to Fig. 4.2.6.2 taken from Bonse & Henning (1986
), which shows that their interferometer results [which measure
directly] and the Kramers–Kronig results differ from one another by
in the neighbourhood of the K-absorption edge of niobium in the compound lithium niobate.
![]() | The relativistic correction in electrons per atom for: (a) the modified form-factor approach; (b) the relativistic multipole approach; (c) the relativistic dipole approach. |
Further theoretical objections have been made by Creagh (1984) and Smith (1987
), who has shown that the Stibius-Jensen correction is not valid, and that, when higher-order multipolar expansions and retardation are considered, the total self-energy correction becomes
rather than
. Fig. 4.2.6.1
shows the variation of the self-energy correction with atomic number for the modified form factor (Creagh, 1984
; Smith, 1987
; Cromer & Liberman, 1970
).
For the imaginary part of the dispersion correction , which depends on the calculation of the photoelectric scattering cross section, better agreement is found between theoretical results and experimental data. Details of this comparison have been given elsewhere (Section 4.2.4
). Suffice it to say that Creagh & Hubbell (1990
), in reporting the results of the IUCr X-ray Attenuation Project, could find no rational basis for preferring the Scofield (1973
) Hartree–Fock calculations to the Cromer & Liberman (1970
, 1981
) and Storm & Israel (1970
) Dirac–Hartree–Fock–Slater calculations.
Computer programs based on the Cromer & Liberman program (Cromer & Liberman, 1983) are in use at all the major synchrotron-radiation laboratories. Many other laboratories have also acquired copies of their program. This program must be modified to remove the incorrect Stibius–Jensen correction term, and, as will be seen later, the energy term should be modified to be
.
Kissel, Pratt & Roy (1980) have developed a computer program based on the second-order S-matrix formalism suggested by Brown, Peierls & Woodward (1955
). Their aim was to provide a prescription for the accurate (∼1%) prediction of the total-atom Rayleigh scattering amplitudes.
Their model treats the elastic scattering as the sum of bound electron, nuclear, and Delbrück scattering cross sections, and treats the Rayleigh scattering by considering second-order, single-electron transitions from electrons bound in a relativistic, self-consistent, central potential. This potential was a Dirac–Hartree–Fock–Slater potential, and exchange was included by use of the Kohn & Sham (1965) exchange model. They omitted radiative corrections.
In principle, the observables in an elastic scattering process are momentum and polarization
. The complex polarization vectors
satisfy the conditions
In quantum mechanics, elastic scattering is described in terms of a differential scattering amplitude, M, which is related to the elastic cross section by equation (4.2.6.16).
If polarization is not an observable, then the expression for the differential scattering cross section takes the form of equation (4.2.6.17). If polarization is taken into account, as may be the case when a polarizer is used on a beam scattered from a sample irradiated by the linearly polarized beam from a synchrotron-radiation source, the full equation, and not equation (4.2.6.17)
, must be used to compute the differential scattering cross section.
The principle of causality implies that the forward-scattering amplitude should be analytic in the upper half of the
plane, and that the dispersion relation
should hold, with the consequence that
This may be rewritten as
with the value of
defined by equation (4.2.6.15)
. Using the conservation of probability,
which is to be compared with equation (4.2.6.23)
.
Starting with Furry's extension of the formalism of quantum mechanics proposed by Feynman and Dyson, the total Rayleigh amplitude may be written as where
and
The
are the complete set of bound and continuum states in the external field of the atomic potential. Singularities occur at all photon energies that correspond to transitions between bound
and bound state
. These singularities are removed if the finite widths of these states are considered, and the energies E are replaced by
, where
is the total (radiative plus non-radiative) width of the state (Gavrila, 1981
). By use of the formalism suggested by Brown et al. (1955
), it is possible to reduce the numerical problems to one-dimensional radial integrals and differential equations. The required multipole expansions of
and the specification of the radial perturbed orbitals that are characterized by angular-momentum quantum numbers have been discussed by Kissel (1977
). Ultimately, all the angular dependence on the photon scattering angle is written in terms of the associated Legendre functions, and all the energy dependence is in terms of multipole amplitudes.
Solutions are not found for the inhomogeneous radial wave equations, and Kissel (1977) expressed the solution as the linear sum of two solutions of the homogeneous equation, one of which was regular at the origin and the other regular at infinity.
Because excessive amounts of computer time are required to use these direct techniques for calculating the amplitudes from all the subshells, simpler methods are usually used for calculating outer-shell amplitudes. Kissel & Pratt (1985) used estimates for outer-shell amplitudes based on the predictions of the modified form-factor approach. A tabulation of the modified relativistic form factors has been given by Schaupp, Schumacher, Smend, Rullhusen & Hubbell (1983
).
Because of the generality of their approach, the computer time required for the calculation of the scattering amplitudes for a particular energy is quite long, so that relatively few calculations have been made. Their approach, however, does not confine itself solely to the problem of forward scattering of photons as does the Cromer & Liberman (1970) approach. Using their model, Kissel et al. (1980
) have been able to show that it is incorrect to assign a dependence of the dispersion corrections on the scattering vector
. This is at variance with some established crystallographic practices, in which the dispersion corrections are accorded the same dependence on
as
, and also at variance with the predictions of Wagenfeld's (1975
) non-relativistic model.
References
































