RTAB: the Rayleigh scattering database
PREVIOUS SECTION | CONTENTS | NEXT SECTION

RTAB data files are available at:
ftp://www-phys.llnl.gov/pub/rayleigh/RTAB

7. Anomalous scattering factors (tables_ASF, data_ASF)

Anomalous scattering factors (ASF) prepared by previous authors has been restricted to keV energies (1-70 keV for Cromer and Liberman, 1-30 keV for Henke et al.), but our studies have required the values over a wider energy range. We have computed ASF values for all atoms over an expanded energy range of 0-10 MeV. These values are tabulated on a variable grid that allows accurate interpolation to intermediate energies. More details of our ASF calculation are provided in Kissel et al. (1995).

Two separate tabulations (stored in the folders 'tables_ASF' and 'data_ASF') of our anomalous scattering factors are provided, that differ in how the energy-dependent part of the bound-bound resonances are stored. The values in 'data_ASF' separately tabulate the real and imaginary anomalous scattering factors on independent grids. Further, only the constant contribution of bound-bound transitions is included in the real anomalous scattering factors; an analytic expression and separately tabulated bound-bound oscillator strengths are needed to compute the full result. As a consequence, these values can be accurately interpolated to all intermediate energies using appropriate algorithms, and can be safely used as input for further computations to the ASFTAB and FFTAB codes, as examples.

The tables in 'tables_ASF' explicitly include the full contribution of the bound-bound transitions and have been prepared by ASFTAB from data in 'data_ASF' folder. While these tables are more readily accessible for immediate use without further computations, they cannot be accurately interpolated to intermediate energies in all cases, as one cannot tabulate the full energy dependence of the bound-bound resonances on a dense enough grid.

In summary, the values in 'data_ASF' have been prepared for subsequent use in further calculations, while the values in 'tables_ASF' have been prepared for direct use without interpolation to intermediate energies.

A variety of notations, phase coventions and normalizations are in common use for anomalous scattering factors. We designate

    (f',f") as the Kissel and Pratt corrections for the relativistic form factor, f(q),

    (g',g") as the Kissel and Pratt corrections for the modified relativistic form factor, g(q),

    (f'CL,f"CL) as the corrections defined by Cromer and Liberman (1970a,b),

    (f1,f2) as the corrections defined by Henke et al (1981,1993),

to indicate the phase and magnitude of the anomalous scattering factors published by these authors. With this notation, we note the following relationships [footnote - These equations are meant to indicate the phase and normalization relationships between the anomalous scattering factors of various authors. It would not be true, for example, that adding N to our value of f' would yield exactly the f1 value as that published by Henke et al. (1993). Instead, adding N to our value of f' would yield a quantity that could be compared directly to the f1 value published by Henke et al.]

g'(omega) = f'(omega) - f'(infinity) ,
g"(omega) = f"(omega) ,
f'CL(omega) = f'(omega) ,
f"CL(omega) = -f"(omega) ,
f1(omega) = N + f'(omega) ,
f2(omega) = f"CL(omega) = -f"(omega) .

Sample anomalous scattering factors extracted from the RTAB database are shown in Table 5.

An interesting feature of our ASF values that differs from other authors is the explicit inclusion of bound-bound resonant transitions. In our underlying model of single-electron transitions in a potential, a bound-bound resonant transition occurs at a single energy (our levels have no widths), the difference of the energies of the two orbitals involved in the transition. This infinitely narrow transition is manifested as a delta-function spike in the imaginary scattering amplitude and a resonance approaching infinity in the real scattering amplitude. Although these explicit spikes and infinities are unphysical, the underlying strength of the transition is important and contributes significantly to the scattering at low energies. The inclusion of bound-bound transitions in our anomalous scattering factors is important for satisfying the Thomas-Reiche-Kuhn sum rule (see, Kissel et al., 1995), wherein an appropriate integral over all energies of the imaginary scattering amplitude should equal the number of electrons in the atom. In many cases, bound-bound transitions contribute 30% or more of the contribution to the TRK sum rule.

A challenge for evaluating anomalous scattering factors via the relativistic dispersion relation and the optical theorem involves energies above 2mc^2. As discussed by Pratt et al. (1994), partitioning of the many particle scattering amplitude to yield R+NT+… introduces an additional contribution to the imaginary amplitude for Rayleigh scattering as computed by the optical theorem. Above 2mc^2, to lowest order, one needs to subtract the contribution of bound-electron pair production from that of photoeffect to correctly compute the total cross section. [footnote - By bound-electron pair production, we mean pair production in which the electron of the pair is created is a bound state of the atom. Ordinarily, discussions of pair production involve the situation where both the electron and positron are created in the continuum.] This insight settles a long standing quandary as it is known that the relativistic photoeffect cross section goes as 1/E at high energies, and as a consequence, the dispersion integral is not convergent if photoeffect were the only contribution to lowest order at high energy. As with photoeffect, bound-electron pair production goes as 1/E at high energies so the integration is naturally cut off at high energies. We utilize an estimate of the bound-electron pair production cross section for the non-relativistic K shell due to Costescu (see, for example, Bergstrom et al., 1997). In our evaluation of the anomalous scattering factors on 0-10 MeV, we carry out the integration for energies of 0-100 MeV.

The validity of this estimate for the total photon-atom cross section above 2mc^2 has not been examined in detail. It is possible that this approximation in our evaluation of the anomalous scattering factors is responsible for the decreasing validity of our MFASF predictions at non-forward angles for energies above 100 keV.

Previous evaluations of anomalous scattering factors have ignored the issue of energies above 2mc^2 by simply cutting off the integration at some suitably high energy. Besides being important for the evaluation of anomalous scattering factors in the 100-10,000 keV range, the behavior of the total cross section above 2mc^2 effects values at lower energies by an overall constant. This is related to the high-energy limit corrections to the anomalous scattering factors of Cromer and Liberman prepared by Kissel and Pratt (1990).

All our anomalous scattering factors satisfy the TRK sum rule to a fraction of 1% (typically at the 0.01% level). The sum rule check is included as part of the output in the ASF files in 'data_ASF'. As an example, the TRK sum rule check in 'data_ASF/006_asf0sl' is

GPRIME: SUM-RULE CHECK ON CROSS SECTIONS;
     COMPUTED=   1.25405E+00 FROM BOUND-BOUND TRANSITIONS
                 4.74291E+00 FROM BOUND-FREE  TRANSITIONS
                 2.07488E-03 FROM HIGH-ENERGY LIMIT
                 5.99903E+00 TOTAL
    PREDICTED=   6.00000E+00
   DIFFERENCE=  -9.69774E-04  (-1.62E-02%)
The resonance structure from inclusion of bound-bound transitions is evident as the narrow peaks in the differential cross sections shown in Figure 1. For use in other contexts, it would likely be useful to convolve our values with a smoothing function of width 1-2 eV, especially for comparisons with predictions of other authors or with experiment in the resonance region.


Table 5. Values of the anomalous scattering factors for selected atoms at 8.04778 keV from the RTAB database.
return to tables

 

atom

 

Filename

 

f′ = f’CL

 

g′ = f′ - f′(¥ )

 

f1 = N + f′

f″= g″

= -f″CL

= -f2

C (Z=6)

tables_ASF/006_asftab0sl

1.91297E-02

2.12046E-02

6.01913E+00

-9.55063E-03

Ne (Z=10)

tables_ASF/010_asftab0sl

1.04035E-01

1.11061E-01

1.01040E+01

-8.57632E-02

Al (Z=13)

tables_ASF/013_asftab0sl

2.16604E-01

2.29718E-01

1.32166E+01

-2.50211E-01

Zn (Z=30)

tables_ASF/030_asftab0sl

-1.55537E+00

-1.46193E+00

2.84446E+01

-6.94718E-01

Sn (Z=50)

tables_ASF/050_asftab0sl

3.73172E-02

3.46528E-01

5.00373E+01

-5.50974E+00

Pb (Z=82)

tables_ASF/082_asftab0sl

-4.08364E+00

-3.08775E+00

7.79164E+01

-8.68796E+00

PREVIOUS SECTION | CONTENTS | NEXT SECTION


LLNL home page and LLNL Disclaimers
UCRL-MI-137026

Information date: Sep. 2, 2000 lk