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. 2.6, pp. 101-103
Section 2.6.1.6.3. Smoothing, desmearing, and Fourier transformation
O. Glattera
|
There are many methods published that offer a solution for this problem. Most are referenced and some are reviewed in the textbooks (Glatter, 1982a; Feigin & Svergun, 1987). The indirect transformation method in its original version (Glatter, 1977a,b, 1980a,b) or in modifications for special applications (Moore, 1980; Feigin & Svergun, 1987) is a well established method used in the majority of laboratories for different applications. This procedure solves the problems of smoothing, desmearing, and Fourier transformation [inversion of equations (2.6.1.9) or (2.6.1.54), (2.6.1.55)] in one step. A short description of this technique is given in the following.
Indirect transformation methods . The indirect transformation method combines the following demands: single-step procedure, optimized general-function system, weighted least-squares approximation, minimization of termination effect, error propagation, and consideration of the physical smoothing condition given by the maximum intraparticle distance. This smoothing condition requires an estimate Dmax as an upper limit for the largest particle dimension: For the following, it is not necessary for to be a perfect estimate, but it must not be smaller than D.
As for , we can use a function system for the representation of p(r) that is defined only in the subspace . A linear combination is used as an approximation to the PDDF. Let N be the number of functions and be the unknowns. The functions are chosen as cubic B splines (Greville, 1969; Schelten & Hossfeld, 1971) as they represent smooth curves with a minimum second derivative.
Now we take advantage of two facts. The first is that we know precisely how to calculate a smeared scattering function from I(h) [equation (2.6.1.56)] and how p(r) or D(R) is transformed into I(h) [equations (2.6.1.9) or (2.6.1.54), (2.6.1.55)], but we do not know the inverse transformations. The second fact is that all these transformations are linear, i.e. they can be applied to all terms in a sum like that in equation (2.6.1.58) separately. So it is easy to start with our approximation in real space [equation (2.6.1.58)] taking into account the a priori information Dmax. The approximation to the ideal (unsmeared) scattering function can be written as where the functions are calculated from by the transformations (2.6.1.9) or (2.6.1.54), (2.6.1.55), the coefficients remain unknown. The final fit in the smeared, experimental space is given by a similar series where the are functions calculated from by the transform (2.6.1.56). Equations (2.6.1.58), (2.6.1.59), and (2.6.1.60) are similar because of the linearity of the transforms. We see that the functions are calculated from in the same way as the data were produced by the experiment from p(r). Now we can minimize the expression where M is the number of experimental points. Such least-squares problems are in most cases ill conditioned, i.e. additional stabilization routines are necessary to find the best solution. This problem is far from being trivial, but it can be solved with standard routines (Glatter, 1977a,b; Tikhonov & Arsenin, 1977).
The whole process of data evaluation is shown in Fig. 2.6.1.14 . Similar routines cannot be used in crystallography (periodic structures) because there exists no estimate for Dmax [equation (2.6.1.57)].
Function systems φv(r); Ψv(h); and χv(h) used for the approximation of the scattering data in the indirect transformation method. |
Maximum particle dimension . The sampling theorem of Fourier transformation (Shannon & Weaver, 1949; Bracewell, 1986) gives a clear answer to the question of how the size of the particle D is related to the smallest scattering angle h1. If the scattering curve is observed at increments starting from a scattering angle h1, the scattering data contain, at least theoretically, the full information for all particles with maximum dimension D The first application of this theorem to the problem of data evaluation was given by Damaschun & Pürschel (1971a,b). In practice, one should always try to stay below this limit, i.e. taking into account the loss of information due to counting statistics and smearing effects. An optimum value for Δh = π/(6D) is claimed by Walter, Kranold & Becherer (1974).
Information content . The number of independent parameters contained in a small-angle scattering curve is given by with and being the lower and upper limits of h. In practice, this limit certainly depends on the statistical accuracy of the data. It should be noted that the number of functions N in equations (2.6.1.58) to (2.6.1.60) may be larger than because they are not independent. They are correlated by the stabilization routine. An example of this problem can be found in Glatter (1980a).
Resolution . There is no clear answer to the question concerning the smallest structural details, i.e. details in the function that can be recognized from an experimental scattering function. The limiting factors are the maximum scattering angle , the statistical error , and the weighting functions , and (Glatter, 1982a). The resolution of standard experiments is not better than approximately 10% of the maximum dimension of the particle for a monodisperse system. In the case of polydisperse systems, resolution can be defined as the minimum relative peak distance that can be resolved in a bimodal distribution. We know from simulations that this value is of the order of 25%.
Special transforms . The PDDF or the size distribution function D(R) is related to I(h) by equations (2.6.1.9) or (2.6.1.54), (2.6.1.55). In the special case of particles elongated in one direction (like cylinders), we can combine equations (2.6.1.41) and (2.6.1.43) and obtain This Hankel transform can be used in the indirect transformation method for the calculation of in (2.6.1.59). Doing this, we immediately obtain the PDDF of the cross section from the smeared experimental data. It is not necessary to know the length L of the particle if the results are not needed on an absolute scale. For this application, we only need the information that the scatterers are elongated in one direction with a constant cross section. This information can be found from the overall PDDF of the particle or can be a priori information from other experiments, like electron microscopy. The estimate for the maximum dimension (2.6.1.57) is related to the cross section in this application, i.e. the maximum dimension of the cross section must not be larger than .
The situation is quite similar for flat particles. If we combine (2.6.1.47) and (2.6.1.49), we obtain being the distance distribution function of the thickness. We have to check that the particles are flat with a constant thickness with maximum thickness . A is the area of the particles and would be needed only for experiments on an absolute scale.
References
Bracewell, R. (1986). Fourier transform and its applications. New York: McGraw-Hill.Google ScholarDamaschun, G. & Pürschel, H.-V. (1971a). Röntgen-Kleinwinkelstreuung von isotropen Proben ohne Fernordnung. I. Allgemeine Theorie. Acta Cryst. A27, 193–197.Google Scholar
Damaschun, G. & Pürschel, H. V. (1971b). Berechnung von Streumassenradien aus unverschmierten und spaltverschmierten Röntgen-Kleinwinkelstreukurven. Monatsh. Chem. 102, 1146–1155.Google Scholar
Feigin, L. A. & Svergun, D. I. (1987). Structure analysis by small-angle X-ray and neutron scattering. New York: Plenum.Google Scholar
Glatter, O. (1977a). Data evaluation in small-angle scattering: calculation of the radial electron density distribution by means of indirect Fourier transformation. Acta Phys. Austriaca, 47, 83–102.Google Scholar
Glatter, O. (1977b). A new method for the evaluation of small-angle scattering data. J. Appl. Cryst. 10, 415–421.Google Scholar
Glatter, O. (1980a). Evaluation of small-angle scattering data from lamellar and cylindrical particles by the indirect transformation method. J. Appl. Cryst. 13, 577–584.Google Scholar
Glatter, O. (1980b). Determination of particle-size distribution functions from small-angle scattering data by means of the indirect transformation method. J. Appl. Cryst. 13, 7–11.Google Scholar
Glatter, O. (1982a). In Small angle X-ray scattering, edited by O. Glatter & O. Kratky, Chap. 4. London: Academic Press.Google Scholar
Greville, T. N. E. (1969). Theory and applications of spline functions. New York: Academic Press.Google Scholar
Moore, P. B. (1980). Small-angle scattering. Information content and error analysis. J. Appl. Cryst. 13, 168–175.Google Scholar
Schelten, J. & Hossfeld, F. (1971). Application of spline functions to the correction of resolution errors in small-angle scattering. J. Appl. Cryst. 4, 210–223.Google Scholar
Shannon, C. E. & Weaver, W. (1949). The mathematical theory of communication. Urbana: University of Illinois Press.Google Scholar
Tikhonov, A. N. & Arsenin, V. Ya. (1977). Solution of ill-posed problems. New York: John Wiley.Google Scholar
Walter, G., Kranold, R. & Becherer, G. (1974). Zu Problemen der Ver- und Entschmierung von Röntgen-Kleinwinkel-Streukurven. Stud. Biophys. 47, 49–62.Google Scholar