Tables for
Volume H
Powder diffraction
Edited by C. J. Gilmore, J. A. Kaduk and H. Schenk

International Tables for Crystallography (2018). Vol. H, ch. 3.6, pp. 292-293

Section Domain size and shape

M. Leonia*

aDepartment of Civil, Environmental and Mechanical Engineering, University of Trento, via Mesiano 77, 38123 Trento, Italy
Correspondence e-mail: Domain size and shape

| top | pdf |

In nanostructured materials, the finite size of the scattering domains is usually the dominant source of line-profile broadening. Actually, when dealing with size, we should consider a size and a shape distribution of the domains. In most cases, one or more distributions of similar objects are considered. For an up-to-date description of issues related to size broadening, see Chapter 5.1[link] . The domain shape is not a property of the material and therefore the use of symmetry constraints [e.g. spherical harmonics to describe the shape of the scattering object as in the model of Popa (1998[link]) or as a size extension of Stephens' (1999[link]) work] is not justified in the general case (Nye, 1987[link]). Exceptions, however, exist.

The size-broadening contribution in WPPM follows the ideas of Bertaut (1949a[link],b[link], 1950[link]) and of Stokes & Wilson (1942[link]). Bertaut proposed the division of the domains into columns and the analysis of the independent scattering of these columns. The column-length distribution can always be extracted from the data: more complex models involving given shapes or distributions simply modify the way in which the columns are rearranged. Stokes and Wilson introduced the concept of a ghost to calculate the Fourier transform for a given shape: the volume common to a domain of shape c and its ghost, i.e. a copy of the same domain displaced by a quantity L along the scattering direction hkl, is proportional to the (size) Fourier transform [A_{c,hkl}^S(L,D)] for the given domain. The calculation has been already carried out analytically for several simple shapes characterized by a single length parameter (Stokes & Wilson, 1942[link]; Lele & Anantharaman, 1966[link]; Wilson, 1969[link]; Langford & Louër, 1982[link]; Vargas et al., 1983[link]; Grebille & Bérar, 1985[link]; Scardi & Leoni, 2001[link]), and can be performed numerically in the general case (Leonardi et al., 2012[link]). It is possible to relate the Fourier coefficients to the size values obtained from traditional methods. In particular, the area-weighted average size <L>S (Warren–Averbach method) and the volume-weighted average size <L>V (Williamson–Hall method) are obtained as[\eqalignno{\langle L \rangle_S &= - \left [\left. dA_c^S (L,D)/dL \right |_{L = 0}\right]^{ - 1} = D/K_k = - D/H_1, &(3.6.22)\cr \langle L \rangle_V &= [\beta(s)]^{-1} = 2\textstyle \int\limits_0^{D/K} A_c^S(L,D)\, {\rm d}L = D/K_\beta, &(3.6.23)}]where β(s) is the integral breadth and Kk and Kβ are the initial slope and integral breadth Scherrer constants, respectively (Langford & Wilson, 1978[link]; Scardi & Leoni, 2001[link]).

The average size might have little physical significance in real cases: the size (and shape) distribution can in fact play a key role in determining both the properties and the diffraction line-profile shapes of the powder under analysis. Fortunately, the Fourier coefficients for the polydisperse case can be easily calculated for any given distribution: the log-normal and the gamma distributions are the most common (and the most flexible). The equations and the corresponding moments are[\eqalignno{g_l(D) &= {1 \over {D\sigma_l \sqrt {2\pi} }} \exp \left[- {{(\ln D - \mu_l)^2} \over {2\sigma_l^2}} \right], \cr M_{l,n} &= \exp\left(n\mu_l + {{n^2} \over 2}\sigma_l^2\right), & (3.6.24)}][\eqalignno{g_g(D) &= {{\sigma_p} \over {\mu_g \Gamma (\sigma_g)}} \left({{\sigma_g D} \over {\mu_g}} \right)^{{\sigma_g}-1} \exp \left(-{{\sigma_g D} \over {\mu_g}}\right), &\cr M_{g,n} &= \left({{\mu_g} \over {\sigma_g}} \right)^n\,{{\Gamma (n + \sigma_g)} \over {\Gamma (\sigma_g)}}. & (3.6.25)}]The scattered intensity for the given distribution gi, and the given shape c, is[\eqalignno {I_{c,hkl} (s) &\propto {{\textstyle \int\limits _0^\infty \left [\int\limits_{L = 0}^{D/K_{hkl}} A_{c,hkl}^S(L,D)\exp (2\pi iLs)\, {\rm d} L \right] w(D)\,{\rm d}D} \over {\textstyle \int\limits_0^\infty w(D)\,{\rm d}D}} \cr & \propto \textstyle \int \limits_{L = 0}^\infty \left [{{\textstyle \int\limits _{D = LK_{hkl}}^\infty A_{c,hkl}^S(L,D)w(D)\, {\rm d}D} \over {\textstyle \int\limits _0^\infty w(D)\, {\rm d}D}} \right] \exp (2\pi iLs)\,{\rm d}L \cr & \propto \textstyle\int \limits_{L = 0}^\infty A_{hkl}^S(L)\exp(2\pi iLs)\,{\rm d}L, & (3.6.26)}]where w(D) = gi(D)Vc(D) and where[A_{hkl}^S(L) = {{\textstyle\int \limits_{D = LK_{hkl}}^\infty A_{c,hkl}^S(L,D)w(D)\,{\rm d} D} \over {\textstyle\int \limits_0^\infty w(D)\, {\rm d}D}}. \eqno (3.6.27)]

With a suitable definition of the Fourier coefficients, the polydisperse case therefore becomes analogous to the monodisperse case. Analytic expressions can be obtained in particular cases. For instance, the Fourier coefficients for the log-normal and gamma distributions (Scardi & Leoni, 2001[link]) are [A_l^S(L) = { \sum\limits_{n = 0}^3} {\rm erfc} \left [{{\ln (L K_c) - \mu_l - (3 - n)\sigma_l^2} \over {\sigma_l \sqrt 2 }} \right] {{M_{l,3 - n}} \over {2M_{l,3}}} H_n^c L^n \eqno (3.6.28)] and [A_g^S(L) = {\sum\limits_{n = 0}^3} \left({{\sigma _g} \over {M_{g,1}}} \right)^n \,{{\Gamma [\sigma_g + (3 - n), K_c L\sigma_g/M_{g,1}]} \over {\Gamma (\sigma_g + 3)}}H_n^c L^n, \eqno (3.6.29)]respectively, where[\eqalign {\Gamma (x,a) &= \textstyle \int \limits_a^\infty y^{x - 1} \exp (-y)\, {\rm d}y, \cr \Gamma (x) &= \Gamma (x,0), \cr {\rm erfc}(x) &= 2\pi^{-1/2}\textstyle\int \limits_x^\infty \exp (-y^2)\,{\rm d}y,}]with the definitions already given in equations (3.6.24)[link] and (3.6.25)[link].

The functional forms of equations (3.6.28)[link] and (3.6.29)[link] clearly suggest that the profile for a log-normal distribution of domains (which is frequently encountered in practice) is not Voigtian: all traditional line-profile analysis methods (based on Voigt or pseudo-Voigt functions) are therefore unable to correctly deal with a log-normally dispersed powder.

By analogy to the monodisperse case, it is possible to relate the parameters of the polydisperse system to the size obtained with traditional methods (Warren–Averbach and Williamson–Hall, respectively). The following holds (Krill & Birringer, 1998[link]; Scardi & Leoni, 2001[link]):[\langle L\rangle_S = {1 \over {K_k}} {{M_{i,3}} \over {M_{i,2}}}, \quad \langle L \rangle_V = {1 \over {K_\beta}} {{M_{i,4}} \over {M_{i,3}}}. \eqno (3.6.30)]

Here, it is clear that diffraction does not provide the first moment of the distribution directly: ratios between high-order moments are involved.

Using an analytical expression for the description of a size distribution can help in stabilizing the results (as the size distribution curve is forced to be zero at very small and very large size values). Some doubts can, however, arise as to the physical validity of this forcing. An example is the case of a multimodal system. The traditional LPA techniques are unable to directly deal with multimodal size distributions. In cases where the multimodal character is clear and the various distribution are well behaved (i.e. when they can be modelled with analytical functions), the pattern can be usually modelled by considering the material as made of different fractions, each of them characterized by a different size distribution.

A possible alternative has been proposed in the literature: replacing the analytical distribution with a histogram. The ability of this model to fit the experimental data has been demonstrated (Leoni & Scardi, 2004[link]; Matěj et al., 2011[link]); a regularization might be necessary to stabilize the shape and/or smoothness of the size distribution. The quality of the measurement and the availability of models describing all contributions to the peak broadening are in most cases the limiting factors for extensive use of the histogram model: correlations of small sizes with the background and with features such as thermal diffuse scattering (Beyerlein et al., 2012[link]) can in fact occur. So far, this is the only available method for exploring cases where the analytical models are unable to correctly describe the observed broadening.


Bertaut, E. F. (1949a). Etude aux rayons X de la répartition des dimensions des cristallites dans une poudre cristalline. C. R. Acad. Sci. 228, 492–494.Google Scholar
Bertaut, E. F. (1949b). Signification de la dimension cristalline mesurée d'apres la largeur de raie Debye-Scherrer. C. R. Acad. Sci. 228, 187–189.Google Scholar
Bertaut, E. F. (1950). Raies de Debye–Scherrer et repartition des dimensions des domaines de Bragg dans les poudres polycristallines. Acta Cryst. 3, 14–18.Google Scholar
Beyerlein, K. R., Leoni, M. & Scardi, P. (2012). Temperature diffuse scattering of nanocrystals. Acta Cryst. A68, 382–392.Google Scholar
Grebille, D. & Bérar, J.-F. (1985). Calculation of diffraction line profiles in the case of a major size effect: application to boehmite AlOOH. J. Appl. Cryst. 18, 301–307.Google Scholar
Krill, C. E. & Birringer, R. (1998). Estimating grain-size distributions in nanocrystalline materials from X-ray diffraction profile analysis. Philos. Mag. A, 77, 621–640.Google Scholar
Langford, J. I. & Louër, D. (1982). Diffraction line profiles and Scherrer constants for materials with cylindrical crystallites. J. Appl. Cryst. 15, 20–26.Google Scholar
Langford, J. I. & Wilson, A. J. C. (1978). Scherrer after sixty years: a survey and some new results in the determination of crystallite size. J. Appl. Cryst. 11, 102–113.Google Scholar
Lele, S. & Anantharaman, T. R. (1966). Influence of crystallite shape on particle size broadening of Debye–Scherrer reflections. Proc. Indian Acad. Sci. A, 64, 261–274.Google Scholar
Leonardi, A., Leoni, M., Siboni, S. & Scardi, P. (2012). Common volume functions and diffraction line profiles of polyhedral domains. J. Appl. Cryst. 45, 1162–1172.Google Scholar
Leoni, M. & Scardi, P. (2004). Nanocrystalline domain size distributions from powder diffraction data. J. Appl. Cryst. 37, 629–634.Google Scholar
Matěj, Z., Matějová, L., Novotný, F., Drahokoupil, J. & Kuzel, R. (2011). Determination of crystallite size distribution histogram in nanocrystalline anatase powders by XRD. Z. Kristallogr. Proc. 1, 87–92.Google Scholar
Nye, J. F. (1987). Physical Properties of Crystals: Their Representation by Tensors and Matrices, reprint edition. Oxford University Press.Google Scholar
Popa, N. C. (1998). The (hkl) dependence of diffraction-line broadening caused by strain and size for all Laue groups in Rietveld refinement. J. Appl. Cryst. 31, 176–180.Google Scholar
Scardi, P. & Leoni, M. (2001). Diffraction line profiles from polydisperse crystalline systems. Acta Cryst. A57, 604–613.Google Scholar
Stephens, P. W. (1999). Phenomenological model of anisotropic peak broadening in powder diffraction. J. Appl. Cryst. 32, 281–289.Google Scholar
Stokes, A. R. & Wilson, A. J. C. (1942). A method of calculating the integral breadths of Debye–Scherrer lines. Math. Proc. Cambridge Philos. Soc. 38, 313–322.Google Scholar
Vargas, R., Louër, D. & Langford, J. I. (1983). Diffraction line profiles and Scherrer constants for materials with hexagonal crystallites. J. Appl. Cryst. 16, 512–518.Google Scholar
Wilson, A. J. C. (1969). Variance apparent particle sizes for cylinders, prisms and hemispheres. J. Appl. Cryst. 2, 181–183.Google Scholar

to end of page
to top of page