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.5, pp. 282-287

Chapter 3.5. Data reduction to |Fhkl| values

A. Le Baila*

aUniversité du Maine, Institut des Molécules et Matériaux du Mans, UMR CNRS 6283, Avenue Olivier Messiaen, 72085 Le Mans Cedex 9, France
Correspondence e-mail:

Methods for the fast extraction of all the peak intensities from a complete powder-diffraction pattern are discussed. The story of the modern whole-powder-pattern decomposition methods (the so-called Pawley and Le Bail methods) is provided together with a list of their main application areas. The conclusion is that the elaboration of these methods represented an essential step towards the possibility to solve more easily, if not routinely, a crystal structure from the diffraction pattern of a powder sample, i.e. from one-dimensional data instead of single-crystal three-dimensional data.

3.5.1. Introduction

| top | pdf |

Collecting the structure-factor amplitudes |Fhkl| is the key step leading to structure solution from powder-diffraction data, just as from single-crystal data. However, there are specific difficulties and pitfalls associated with powder data, mainly because of diffraction-peak overlap. Once indexing is realized, data reduction to |Fhkl| is a fast process, using whole-powder-pattern decomposition (WPPD) methods. This comfortable situation was not attained without past efforts, which are reviewed in this chapter. The introduction of modern WPPD methods occurred slowly and progressively over the past 30 years, thanks to increases in computer power, improvements in graphical user interfaces, diffractometer data digitalization, the availability of synchrotron and neutron radiation, and last but not least, the proposition of new algorithms. Innovations were not instantly accepted, this also being true for the Rietveld (1969[link]) method, or could not be applied immediately to every type of powder data. Predecessors of the current WPPD methods extracted peak intensities without restraining the cell, so that each peak position was a parameter to be refined (as well as the peak intensity, and the peak shape and width). This is still useful if the aim is to obtain peak positions for indexing, although simple derivative methods can make searching for peak positions faster. Taking advantage of the indexing (Bergmann et al., 2004[link]), new WPPD methods that applied cell restraints to the peak positions opened the door to a long list of new possibilities and applications (including first indexing confirmation and manual or automatic space-group estimation) which are detailed in this chapter. A partial review of the applications realized in thousands of published papers is given, and the evolution of the methods will be discussed. Additional information on the topic of reduction to |Fhkl| values can be found in the books by Young (1993[link]), Giacovazzo (1998[link]), David et al. (2002[link]), Pecharsky & Zavalij (2003[link]), Clearfield et al. (2008[link]) and Dinnebier & Billinge (2008[link]) or in selected reviews (Toraya, 1994[link]; Langford & Louër, 1996[link]; Le Bail, 2005[link]).

3.5.2. Algorithms

| top | pdf |

Whole-powder-pattern fitting (WPPF) is a general definition including WPPD as well as the Rietveld method (Rietveld, 1969[link]). In the latter method, the atomic coordinates are required for the intensity calculations, and the sum of all the peak contributions produces a calculated powder pattern which is compared to the observed one, allowing the least-squares refinement of profile and structural parameters. The Rietveld method historically preceded modern WPPD methods, though the latter are applicable without atomic coordinates. Of course, one may use WPPD methods if the structure is known, but in some cases one does not want to use that knowledge (not wanting to restrain the peak intensities by the structural model, for instance, nevertheless believing in the indexing or wanting to confirm it, using the restraint of the cell parameters). Any WPPF approach should be able to model the peak shape and width variations with diffraction angle (complications not considered here may occur in the case of anisotropic broadening intrinsic to the sample). This can be done by fitting some analytical profile shape and width parameters in a semi-empirical approach. The angular variation of these parameters is generally controlled by refining the U, V and W terms in the Caglioti et al. (1958[link]) expression [(FWHM)2 = U tan2 θ + V tan θ + W] or a variation (where FWHM = full width at half maximum). The alternative is to use the fundamental parameter approach (FPA) (Cheary & Coelho, 1992[link]). However, some of the original computer programs did not apply any cell restraint or even any restraint at all. Unrestrained cell

| top | pdf |

Without a cell hypothesis, no |Fhkl| values can be extracted; the intensity values collected will be noted by I(i) until Miller indices are attributed, enabling the multiplicity correction. Obtaining all the peak positions, areas, breadths and shape parameters as independent values for a whole powder pattern is limited to simple cases where there is not too much peak overlap. With such an approach (both cell and space group unknown or unused) one has to estimate the number of peaks to be fitted, so that the fit of a complex group of peaks will lead to large uncertainties. However, knowing the cell and space group provides at least the correct number of peaks and an estimate of their starting positions. Such calculations were made as an alternative to the Rietveld method, during the first stage of the so-called two-stage method for refinement of crystal structures (Cooper et al., 1981[link]). In the case of X-ray data, the profile shapes applied in the Rietveld method (Gaussian at the beginning for neutron data) evolved a great deal (Wiles & Young, 1981[link]), and on the WPPD side happened to be described in these two-stage approaches by a sum of Lorentzian curves, or double Gaussians (Will et al., 1983[link], 1987[link]). The computer program PROFIT (Scott, 1987[link]), derived from software for individual profile fitting (Sonneveld & Visser, 1975[link]) and extended to the whole pattern, was applied to the study of crystallite size and strain in zinc oxide (Langford et al., 1986[link]) and for the characterization of line broadening in copper oxide (Langford & Louër, 1991[link]). Studying a whole pattern can also be done in simple cases by using software designed for the characterization of single or small groups of peaks; an example is a ZnO study (Langford et al., 1993[link]) using the computer program FIT (Socabim/Bruker). WPPD on complex cases is mostly realized today by using peak positions controlled by the cell parameters, with the benefit of stronger accuracy of the |Fhkl| values, even if the lost degrees of freedom may lead to slightly worse fits, increasing the profile R factors. Before 1987, close to thirty structure determinations by powder diffractometry (SDPDs) were achieved using intensities extracted by using these old WPPD methods without cell constraints (see the SDPD database; Le Bail, 2007[link]). It can be argued that freeing the peak positions allows one to take into account subtle effects in position displacement (in stressed samples, for example). But systematic discrepancy of observed peak positions with regard to the theor­etical position, as expected from the cell parameters, can be modelled as well in modern WPPD methods or in the Rietveld method. Restrained cell

| top | pdf |

WPPD methods that strictly applied the peak positions calculated from a cell (hypothesized from indexing results) marked a great step forward in the quest for improving the SDPD success rate. This is essentially because the quality of the estimated intensities increased, and even if the main handicap of powder diffraction (peak overlap) could not be completely circumvented, it was at least more clearly delimited. Today two generic names are retained for such cell-restrained WPPD methods that can produce a set of extracted |Fhkl| values suitable for attempting a structure solution: the Pawley and Le Bail methods. Both were derived from the Rietveld (1969[link]) method, so they share with it many equations which will not be restated here (see Chapter 4.7[link] ). Pawley method

| top | pdf |

The idea of removing the crystal structure refinement part in a Rietveld program and adding the potential to refine an individual intensity for every expected Bragg peak produced a new software package (named ALLHKL) allowing refinement of the cell parameters very precisely and extraction of a set of structure-factor amplitudes (Pawley, 1981[link]). The process was much later called the `Pawley method'. Overcoming the least-squares ill conditioning due to peak overlap was achieved by using slack constraints (Waser, 1963[link]). Pawley clearly insisted on the usefulness of the procedure for the confirmation of the indexing of a powder pattern of an unknown. Nevertheless, the structure of the C6F10 (at 4.2 K) test case selected for demonstration purposes remained unsolved (but see Section[link] below). No SDPD of an unknown was realized using the Pawley method for several years (although successful tests were published corresponding to redeterminations of previously known structures). The first real SDPD of an unknown using the Pawley method seems to be that of I2O4 (Lehmann et al., 1987[link]); its powder pattern had been previously indexed, but the structure not determined because of the lack of a suitable single crystal. During these pioneering years, ALLHKL could not extract the intensities for more than 300 peaks, so that, in more complex cases, it was necessary to subdivide the pattern into several parts. Moreover, it was rather difficult to avoid completely the ill conditioning due to overlapping peaks. Successful fits yielded equipartitioned intensities (i.e., equal structure factors for those Bragg peaks with exact overlap). Unsuccessful fits could easily produce negative intensities which, combined with positive ones for other peak(s) at the same angle, reproduced the global positive value. Moreover, the first version to apply Gaussian peak shapes could not easily produce any SDPD because of the relatively poor resolution of constant-wavelength neutron data, so that it needed to be adapted to X-ray data, with the implementation of more complex peak shapes. Several programs were subsequently developed, based on the same principles as the original Pawley method. The first of them, by Toraya (1986[link]), extended the use to X-ray data with non-Gaussian profile shapes, and introduced two narrow band matrices instead of a large triangular matrix, saving both computation time and memory space in a program named WPPF. Some programs were used to produce intensities in order to apply the so-called two-stage method (Cooper et al., 1981[link]) for structure refinement, such as PROFIT (Scott, 1987[link]) and PROFIN (Will, 1988[link]) (no slack constraints, but equal division of the intensity between expected peaks when the overlap was severe). There was intense continuing activity on Pawley-like software with other programs such as FULFIT (Jansen et al., 1988[link]), LSQPROF (Jansen et al., 1992[link]) and POLISH (Byrom & Lucas, 1993[link]). Estimation of intensities of overlapping reflections was improved in LSQPROF by applying relations between structure-factor amplitudes derived from direct methods, and the Patterson function was considered in the satellite program DOREES (Jansen et al., 1992[link]). The question of how to determine the intensities of completely (or largely) overlapping reflections (either systematic overlap due to symmetry or fortuitous overlap) from a single powder pattern cannot have a definite simple answer, but continues to be discussed, since it is essential for improving our ability to solve structures. An early view with a probabilistic approach was given by David (1987[link]), later introducing Bayesian statistics (Sivia & David, 1994[link]) into the Pawley method. Early detection of preferred orientation on the basis of analysis of the E-value distribution was another way (Peschar et al., 1995[link]) to improve the structure-factor-amplitude estimate. New computer programs based on the Pawley method continue to be written even today. Le Bail method

| top | pdf |

In order to be able to estimate R factors related to integrated intensities, Rietveld (1969[link]) stated [see also the book The Rietveld Method edited by Young (1993)]: `a fair approximation to the observed integrated intensities can be made by separating the peaks according to the calculated values of the integrated intensities,' i.e.[I_{hkl}({\rm obs}) = \textstyle\sum\limits_j [ w_{j,hkl}S^2_{hkl}({\rm calc})y_j({\rm obs})/y_j({\rm calc}) ], \eqno(3.5.1)]where wj,hkl is a measure of the contribution of the Bragg peak at position 2θhkl to the diffraction profile yj at position 2θj [corresponding to equation 7 in Rietveld (1969[link])]. [S^2_{hkl}] is the sum of the nuclear and magnetic contributions for neutron diffraction, or is more simply [F^2_{hkl}] for X-rays. The sum is over all yj(obs) that can theoretically contribute to the integrated intensity Ihkl(obs). Bias is introduced here by apportioning the intensities according to the calculated intensities; this is why the observed intensities are said to be `observed', in quotation marks, in the Rietveld method. These `observed' intensities are used in the RB and RF calculations (residuals on intensities and structure-factor amplitudes, respectively). They are also required for Fourier-map estimations, which, as a consequence, are less reliable than those from single-crystal data.

A process using the Rietveld decomposition formula iteratively for WPPD purposes was first applied in 1988 (Le Bail et al., 1988[link]) and much later was called the `Le Bail method' or `Le Bail fit', or `pattern matching' as well as `profile matching' in the FULLPROF Rietveld program (Rodriguez-Carvajal, 1990[link]). In the original computer program (named ARITB), arbitrarily all equal [S^2_{hkl}({\rm calc})] values are first entered in the above equation, instead of using structure factors calculated from the atomic coordinates, resulting in `Ihkl(obs)' which are then re-entered as new [S^2_{hkl}({\rm calc})] values at the next iteration, while the usual profile and cell parameters (but not the scale) are refined by least squares (ARITB used profile shapes represented by Fourier series, either analytical or learned from experimental data, providing an easy way to realize convolution by broadening functions modelling size–strain sample effects, possibly anisotropic). Equipartition of exactly overlapping reflections comes from the strictly equal result from equation (3.5.1)[link] for Bragg peaks at the same angles which would have equal starting calculated intensities. Not starting from a set of all equal [S^2_{hkl}({\rm calc})] values avoids equipartition for the exactly overlapping reflections but produces Ihkl(obs) keeping the same original ratio as the [S^2_{hkl}({\rm calc})] ones. It is understandable that such an iterative process requires starting cell and profile parameters as good as the Rietveld method itself. The process is easier to incorporate within an existing Rietveld code than the Pawley method, so that most Rietveld codes now include structure-factor amplitudes extraction as an option (generally multiphase), with the possibility of combining Rietveld refinement(s) together with Le Bail fit(s).

A non-exhaustive list of programs applying this method (either exclusively or added within a Rietveld code) includes MPROF (Jouanneaux et al., 1990[link]), later renamed WinMPROF; FULLPROF (Rodriguez-Carvajal, 1990[link]); EXTRACT (Baer­locher, 1990[link]); EXTRA (Altomare et al., 1995[link]); EXPO (Altomare et al., 1999[link]), which is the integration of EXTRA and SIRPOW.92 for solution and refinement of crystal structures; and RIETAN (Izumi & Ikeda, 2000[link]). Then followed most well known Rietveld codes (BGMN, GSAS, MAUD, TOPAS etc.) or standalone programs (AJUST by Rius et al., 1996). In the work of the Giacovazzo group, many modifications of the [|F_{hkl}|] values for SDPD purposes were applied before or after the extraction and were integrated in EXPO2011 (Altomare et al., 2011[link]): obtaining information about the possible presence of preferred orientation by statistical analysis of the normalized structure-factor moduli; using the positivity of the Patterson function in the decomposition process, this having been considered previously (David, 1987[link]; Estermann & Gramlich, 1993[link]); characterization of pseudotranslational symmetry used as prior information in the pattern-decomposition process; multiple Le Bail fits with random attribution of intensity to the overlapping reflections, instead of equipartition, followed by application of direct methods to large numbers of such data sets; use of a located structure fragment for improving the pattern-decomposition process; and use of probability (triplet-invariant distribution functions) integrated with the Le Bail algorithm. Another approach for solving the overlapping problem was proposed by using maximum-entropy coupled with likelihood evaluation (Dong & Gilmore, 1998[link]). The list of structure solutions made from intensities extracted by using the Pawley and Le Bail methods is too long to be given here; a partial list (>1000 first cases, including those using |Fhkl| values extracted by other methods) can be found on the web (Le Bail, 2007[link]). The first application of the Le Bail method was to the structure solution of LiSbWO6 (Le Bail et al., 1988[link]) using the ARITB software.

3.5.3. Pitfalls in the extraction of accurate |Fhkl| values using the Pawley and Le Bail methods

| top | pdf |

In the Rietveld refinement guidelines published by the IUCr Commission on Powder Diffraction (McCusker et al., 1999[link]), it is said that the Rietveld Rwp value should approach the value obtained in a structure-free refinement (i.e. using WPPD methods). Such a refinement is recommended for the estimation of initial values of the Rietveld profile parameters. Consequently, |Fhkl| values extracted by WPPD can be used to make a good reproduction of the experimental powder pattern if the cell is correct (which is ultimately only proven if the structure is solved and refined). Pitfalls can occur during post-treatment and application of the |Fhkl| data if one neglects the possible errors that are inherently present due to exact or accidental overlap, preferred orientation effects or wrong background estimations, citing only three of the main possible causes of errors in these |Fhkl| values. Consequences of (exact or accidental) overlap

| top | pdf |

The uncertainties of the |Fhkl| values of overlapped reflections cannot be overcome in a single powder-diffraction experiment. This problem has led to various approaches, all being more or less inefficient: equipartition, non-equipartition by random distribution etc. If direct methods are applied, the trend is to multiply the number of solution attempts, trying to identify the most convinc­ing one by using structural arguments (such as atoms in chemically reasonable positions). When applying real-space methods (which require chemical knowledge, such as the three-dimensional molecular structure or the presence of definite polyhedra) one generally chooses to work either directly on the raw powder pattern or on a pseudo pattern built from the extracted |Fhkl| values, so that wrong individual values are less of a problem, since only the sums of the contributions in overlapping regions are checked during the search for the molecule, polyhedra or atom positions. Indeed, working on the raw powder pattern does not need reduction to |Fhkl| values in theory, but in practice either the Pawley or Le Bail methods are applied first in order to fix the zero point, background, cell and profile parameters which will then be applied during the structure model checking, and to speed the calculations. The extracted |Fhkl| values can be used in mathematical expressions defining correlations induced by the overlap. These equations were developed by David et al. (1998[link]) for the Pawley method in the real-space structure solution program DASH and by Pagola et al. (2000[link]) for the Le Bail method in PSSP. Regenerating a powder pattern from the extracted |Fhkl| values was carried out in the ESPOIR real-space computer program (Le Bail, 2001[link]) using a simple Gaussian peak shape whose width follows the Caglioti relation established from the raw pattern. With such a pseudo powder pattern (without profile asymmetry, background etc.), the calculations are much faster than if the raw pattern is used. When using direct methods instead of real-space methods, the approaches are different, because direct methods require a more complete data set (up to d = 1 Å) of accurate |Fhkl| values. However, removing up to half of them (those with too much overlap, i.e. where the overlap is greater than half the FWHM, for instance) can lead to success with direct methods. One can even remove up to 70–80% of the data if the Patterson method is applied and if only a small number of heavy atoms are to be located. Preferred-orientation effects

| top | pdf |

One has to ensure that preferred orientation is minimized during the recording of the powder pattern if the extracted |Fhkl| values are to be used for structure solution. In transmission geometry with a capillary specimen (provided that it is not composed of long needle-shaped particles that are all aligned), there is generally no problem. But in reflection geometry with samples pressed on the holder, preferred orientation is not rare, even if it is not obvious in the data. Collecting a second pattern from a sample dusted onto the holder through a fine sieve can be informative. However, some WPPD applications may not be sensitive to such a problem. If only the cell parameters have to be accurately estimated for thermal-expansion studies (Evans et al., 1996[link]), it can be much faster to use WPPD rather than the Rietveld method. However, it is not recommended to do this systematically, especially if the structure is complex and the resolution is low [see the warnings in Peterson (2005[link])]. Background-estimation effects

| top | pdf |

Logically, a background value will never be higher than the observed intensity at the diffraction angle where it is visually estimated. In the first refinement cycles by the Pawley or Le Bail methods, it is preferable to keep the background fixed as well as the cell parameters, which assumes that the starting values have been carefully estimated or even (in the cell-parameter case) have already been refined from the peak positions that were used for indexing. This is because all refinement processes need to start from parameters close to the final values. Selecting background values leading to negative intensities after background removal could result in negative |Fhkl| values if the software does not account for this.

3.5.4. Applications and by-products

| top | pdf |

The first modern WPPD method with cell restraints was developed for neutron data by Pawley (1981[link]), 12 years after the publication of the paper that described the Rietveld (1969[link]) method. In Le Bail et al. (1988[link]) a new WPPD approach was used to extract intensities, making use of iterations of the Rietveld decomposition formula. It is clear that both these WPPD methods are children of the Rietveld method. Today most users of the Rietveld method do not cite the original Rietveld papers, but only refer to the computer program that they used. This is also now increasingly the case for the WPPD methods.

From the Thomson Reuters ISI citation index consulted in May 2015, the papers for the Pawley and Le Bail methods scored 892 and 1425 direct citations, respectively. There are several highly cited papers that then cite these two papers. The most highly cited paper (>5100 times) that cites both WPPD methods concerns use of the Le Bail intensity-extraction method by FULLPROF for solving magnetic structures (Rodriguez-Carvajal, 1993[link]). This paper is also given as a reference for FULLPROF used in more standard Rietveld refinements. This suggests that the impact of WPPD methods is higher than commonly believed. The list of possible WPPD applications includes phase identification, quantitative phase analysis, measurement of crystallite sizes and strains, creation of Fourier maps for partially solved structures, structure refinement by the two-step method, studies of electron-density distribution, and characterization of pole figures, using either the Pawley or Le Bail methods. All routes to SDPD use at least one of them. WPPD has even entered into the indexing step, with Kariuki et al. (1999[link]) using the Le Bail fit for testing cell hypotheses (for which it is faster than the Pawley method) in a new computer program that uses a genetic algorithm. But the main applications of the WPPD methods are to provide support for cell-parameter refinement and the determination of the space group, as a prelude to the full use of the extracted |Fhkl| values for ab initio structure solution. Supporting indexing and space-group determination

| top | pdf |

As they yield the smallest possible profile Rp and Rwp factors (smaller than from the Rietveld method, which is limited by the crystal-structure refinement), the Pawley and Le Bail methods provide strong support for both the proposed indexing and the determination of the space group. Some computer programs provide an automated suggestion for the latter. This support is needed to show that it is worth attempting to solve the structure. Once the structure is solved, the structure constraint will remove the ambiguity between intensities of close Bragg peaks and necessarily improve the quality of the cell parameters. If the structure is already known, the best approach is the Rietveld method. There is a progression in the precision of the refined cell parameters from the lowest level (least squares from individually extracted peak positions) to a medium level (WPPD with cell restraint) to the highest possible level (Rietveld, adding the structure constraint). With both Pawley and Le Bail methods, the fit quality is checked using agreement factors which are the same as with the Rietveld method: Rp and Rwp (moreover, a careful visual check is recommended). The reliabilities relative to the structure (RB and RF), which can still be calculated, are meaningless (WPPD programs tending to obtain a value close to zero for both of them). Structure solution

| top | pdf |

SDPD can be undertaken by various approaches, depending on the chemical knowledge of the sample (formula, molecular formula, presence of defined polyhedra…), either directly using the |Fhkl| values for structure solution by direct or Patterson methods, or by rebuilding a pseudo powder pattern from them, or by applying fixed profile parameters from the Pawley or Le Bail fits during whole-powder-pattern fitting wherein the structure solution is attempted by real-space methods. In order to illustrate the power of WPPD methods and to show the progress realized over the last 30 years, the decafluorocyclohexene structure that was unsolved in the Pawley method paper of 1981[link] is reconsidered. As stated by Pawley, from plausible extinctions the space group of the C6F10 crystal structure at 4.2 K could well be P21/n. The |Fhkl| values were extracted from the rebuilt neutron powder pattern by applying the Le Bail method and used for attempting the structure solution by real-space methods. The neutron powder pattern was rebuilt from the 109 intensities extracted up to 54° 2θ, in space group P2/m, given in Table 2 of the original paper. The fit (using FULLPROF) in P21/n of the data rebuilt in P2/m is satisfactory (Fig. 3.5.1[link]). The three-dimensional C6F10 molecule was rotated and translated (six degrees of freedom) in the cell using the ESPOIR (Le Bail, 2001[link]) Monte Carlo program, leading to a plausible starting model (Rp = 13.6%) ready for Rietveld refinement. This program builds a pseudo powder pattern from the extracted |Fhkl| values, which is then compared to the data calculated from the model (Fig. 3.5.2[link]). Unrefined atomic coordinates are available from the Crystallography Open Database (COD, CIF No. 3500009) (Grazulis et al., 2009[link]); a projection of the corresponding structure is shown in Fig. 3.5.3[link]. The true crystal structure is apparently more complex (Solovyov et al., 2014[link]). Final resolution of the structure will require collection of a better experimental powder pattern. However, the coordinates have been refined by energy minimization in the solid state (Smrčok et al., 2013[link]).

[Figure 3.5.1]

Figure 3.5.1 | top | pdf |

Data reduction to |Fhkl| values for the C6F10 Pawley (1981[link]) test case by the Le Bail method using FULLPROF. The neutron powder pattern (4.2 K) was rebuilt (λ = 1.909 Å) from the intensities given in the original paper (P2/m). The extraction of |Fhkl| values was carried out in the space group P21/n.

[Figure 3.5.2]

Figure 3.5.2 | top | pdf |

The C6F10 Monte Carlo molecule positioning by the real-space ESPOIR program produces that best fit (Rp = 13.6%) of the pseudo powder pattern built from the previously extracted |Fhkl| values (Fig. 3.5.1[link]), overcoming the equipartition problem at the reduction stage. Compared to Fig. 3.5.1[link], which shows intensities, the multiplicity and geometrical factors are removed, leading to structure-factor amplitudes.

[Figure 3.5.3]

Figure 3.5.3 | top | pdf |

Projection along the b axis of the C6F10 structure model in P21/n before Rietveld refinement.

3.5.5. Conclusion

| top | pdf |

`Which is best: the Pawley or the Le Bail method?' is not a question with a simple conclusive answer. The fact is that both methods are able to estimate structure-factor amplitudes, which can lead to structure solution from powder-diffraction data in a more efficient way than was previously possible, even if the problem of peak overlap precludes attaining single-crystal quality data from only one powder pattern. The advantage of the Le Bail method over the Pawley method is its speed, which becomes apparent when several thousands of |Fhkl| values have to be extracted, for instance from high-resolution synchrotron data of a complex compound, since no more than about ten parameters have to be refined instead of thousands. Chemical knowledge may reduce the number of |Fhkl| values necessary for solving a structure to the first 100 Bragg peaks at low diffraction angles, as shown in the above example. Moreover, the small number of successful participants in the three SDPD round robins held in 1998, 2002 and 2008 (Le Bail et al., 2009[link]) did not allow us to conclude whether one approach is really better than the other, or even to be sure if all the further modifications of the extracted |Fhkl| values in particular computer programs are really decisive improvements (the conclusion was that SDPD `on demand' was still not an easy task). WPPD is not the only reef on the SDPD journey; indexing remains a considerable bottleneck for complex materials and low-resolution data.


Altomare, A., Burla, M. C., Camalli, M., Carrozzini, B., Cascarano, G. L., Giacovazzo, C., Guagliardi, A., Moliterni, A. G. G., Polidori, G. & Rizzi, R. (1999). EXPO: a program for full powder pattern decomposition and crystal structure solution. J. Appl. Cryst. 32, 339–340.Google Scholar
Altomare, A., Burla, M. C., Cascarano, G., Giacovazzo, C., Guagliardi, A., Moliterni, A. G. G. & Polidori, G. (1995). EXTRA: a program for extracting structure-factor amplitudes from powder diffraction data. J. Appl. Cryst. 28, 842–846.Google Scholar
Altomare, A., Cuocci, C., Giacovazzo, C., Moliterni, A. G. G. & Rizzi, R. (2011). Obtaining models suitable for Rietveld refinement: the EXPO2011 techniques. Z. Kristallogr. 226, 869–881.Google Scholar
Baerlocher, Ch. (1990). EXTRACT, a Fortran program for the extraction of integrated intensities from a powder pattern. Institut für Kristallographie, ETH, Zürich, Switzerland.Google Scholar
Bergmann, J., Le Bail, A., Shirley, R. & Zlokazov, V. (2004). Renewed interest in powder diffraction data indexing. Z. Kristallogr. 219, 783–790.Google Scholar
Byrom, P. G. & Lucas, B. W. (1993). POLISH: computer program for improving the accuracy of structure-factor magnitudes obtained from powder data. J. Appl. Cryst. 26, 137–139.Google Scholar
Caglioti, G., Paoletti, A. & Ricci, F. P. (1958). Choice of collimators for a crystal spectrometer for neutron diffraction. Nucl. Instrum. 3, 223–228.Google Scholar
Cheary, R. W. & Coelho, A. (1992). A fundamental parameters approach to X-ray line-profile fitting. J. Appl. Cryst. 25, 109–121.Google Scholar
Clearfield, C., Reibenspies, J. & Bhuvanesh, N. (2008). Principles and applications of powder diffraction, pp. 261–309. Oxford: Wiley.Google Scholar
Cooper, M. J., Rouse, K. D. & Sakata, M. (1981). An alternative to the Rietveld profile refinement method. Z. Kristallogr., 157, 101–117.Google Scholar
David, W. I. F. (1987). The probabilistic determination of intensities of completely overlapping reflections in powder diffraction patterns. J. Appl. Cryst. 20, 316–319.Google Scholar
David, W. I. F., Shankland, K., McCusker, L. B. & Baerlocher, Ch. (2002). Structure Determination from Powder Diffraction Data, IUCr Monographs on Crystallography, Vol 13. New York: Oxford University Press.Google Scholar
David, W. I. F., Shankland, K. & Shankland, N. (1998). Routine determination of molecular crystal structures from powder diffraction data. Chem. Commun. pp. 931–932.Google Scholar
Dinnebier, R. E. & Billinge, S. J. L. (2008). Powder diffraction: theory and practice, pp. 134–165. Cambridge: RSC Publishing.Google Scholar
Dong, W. & Gilmore, C. J. (1998). The ab initio solution of structures from powder diffraction data: the use of maximum entropy and likelihood to determine the relative amplitudes of overlapped reflections using the pseudophase concept. Acta Cryst. A54, 438–446.Google Scholar
Estermann, M. A. & Gramlich, V. (1993). Improved treatment of severely or exactly overlapping Bragg reflections for the application of direct methods to powder data. J. Appl. Cryst. 26, 396–404.Google Scholar
Evans, J. S. O., Mary, T. A., Vogt, T., Subramanian, M. A. & Sleight, A. W. (1996). Negative thermal expansion in ZrW2O8 and HfW2O8. Chem. Mater. 8, 2809–2823.Google Scholar
Giacovazzo, C. (1998). Direct Phasing in Crystallography: Fundamentals and Applications, IUCr Monographs on Crystallography, Vol. 8, pp. 410–444. New York: Oxford University Press.Google Scholar
Grazulis, S., Chateigner, D., Downs, R. T., Yokochi, A. F. T., Quirós, M., Lutterotti, L., Manakova, E., Butkus, J., Moeck, P. & Le Bail, A. (2009). Crystallography Open Database - an open-access collection of crystal structures. J. Appl. Cryst. 42, 726–729.Google Scholar
Izumi, F. & Ikeda, T. (2000). A Rietveld-analysis program RIETAN-98 and its applications to zeolites. Mater. Sci. Forum, 321–324, 198–205.Google Scholar
Jansen, E., Schäfer, W. & Will, G. (1988). Profile fitting and the two-stage method in neutron powder diffractometry for structure and texture analysis. J. Appl. Cryst. 21, 228–239.Google Scholar
Jansen, J., Peschar, R. & Schenk, H. (1992). On the determination of accurate intensities from powder diffraction data. I. Whole-pattern fitting with a least-squares procedure. J. Appl. Cryst. 25, 231–236.Google Scholar
Jouanneaux, A., Murray, A. D. & Fitch, A. N. (1990). MPROF: a multipattern Rietveld refinement program for neutron, X-ray and synchrotron radiation. SERC Daresbury Laboratory, Warrington, England.Google Scholar
Kariuki, B. M., Belmonte, S. A., McMahon, M. I., Johnston, R. L., Harris, K. D. M. & Nelmes, R. J. (1999). A new approach for indexing powder diffraction data based on whole-profile fitting and global optimization using a genetic algorithm. J. Synchrotron Rad. 6, 87–92.Google Scholar
Langford, J. I., Boultif, A., Auffrédic, J. P. & Louër, D. (1993). The use of pattern decomposition to study the combined X-ray diffraction effects of crystallite size and stacking faults in ex-oxalate zinc oxide. J. Appl. Cryst. 26, 22–33.Google Scholar
Langford, J. I. & Louër, D. (1991). High-resolution powder diffraction studies of copper(II) oxide. J. Appl. Cryst. 24, 149–155.Google Scholar
Langford, J. I. & Louër, D. (1996). Powder diffraction. Rep. Prog. Phys. 59, 131–234.Google Scholar
Langford, J. I., Louër, D., Sonneveld, E. J. & Visser, J. W. (1986). Applications of total pattern fitting to a study of crystallite size and strain in zinc oxide powder. Powder Diffr. 1, 211–221.Google Scholar
Le Bail, A. (2001). ESPOIR: A program for solving structures by Monte Carlo analysis of powder diffraction data. Mater. Sci. Forum, 378, 65–70.Google Scholar
Le Bail, A. (2005). Whole powder pattern decomposition methods and applications: A retrospection. Powder Diffr. 20, 316–326.Google Scholar
Le Bail, A. (2007). Structure Determination from Powder Diffraction Database – SDPD. .Google Scholar
Le Bail, A., Cranswick, L. M. D., Adil, K., Altomare, A., Avdeev, M., Cerny, R., Cuocci, C., Giacovazzo, C., Halasz, I., Lapidus, S. H., Louwen, J. N., Moliterni, A., Palatinus, L., Rizzi, R., Schilder, E. C., Stephens, P. W., Stone, K. H. & van Mechelen, J. (2009). Third structure determination by powder diffractometry round robin (SDPDRR-3). Powder Diffr. 24, 254–262.Google Scholar
Le Bail, A., Duroy, H. & Fourquet, J. L. (1988). Ab-initio structure determination of LiSbWO6 by X-ray powder diffraction. Mater. Res. Bull. 23, 447–452.Google Scholar
Lehmann, M. S., Christensen, A. N., Fjellvåg, H., Feidenhans'l, R. & Nielsen, M. (1987). Structure determination by use of pattern decomposition and the Rietveld method on synchrotron X-ray and neutron powder data; the structures of Al2Y4O9 and I2O4. J. Appl. Cryst. 20, 123–129.Google Scholar
McCusker, L. B., Von Dreele, R. B., Cox, D. E., Louër, D. & Scardi, P. (1999). Rietveld refinement guidelines. J. Appl. Cryst. 32, 36–50.Google Scholar
Pagola, S., Stephens, P. W., Bohle, D. S., Kosar, A. D. & Madsen, S. K. (2000). The structure of malaria pigment beta-haematin. Nature, 404, 307–310.Google Scholar
Pawley, G. S. (1981). Unit-cell refinement from powder diffraction scans. J. Appl. Cryst. 14, 357–361.Google Scholar
Pecharsky, V. K. & Zavalij, P. Y. (2003). Fundamentals of powder diffraction and structural characterization of materials. New York: Springer.Google Scholar
Peschar, R., Schenk, H. & Capkova, P. (1995). Preferred-orientation correction and normalization procedure for ab initio structure determination from powder data. J. Appl. Cryst. 28, 127–140.Google Scholar
Peterson, V. K. (2005). Lattice parameter measurement using Le Bail versus structural (Rietveld) refinement: a caution for complex, low symmetry systems. Powder Diffr. 20, 14–17.Google Scholar
Rietveld, H. M. (1969). A profile refinement method for nuclear and magnetic structures. J. Appl. Cryst. 2, 65–71.Google Scholar
Rius, J., Sane, J., Miravitlles, C., Amigo, J. M., Reventos, M. M. & Louër, D. (1996). Determination of crystal structures from powder diffraction data by direct methods: extraction of integrated intensities from partially overlapping Bragg reflections. Anal. Quim. 92, 223–227.Google Scholar
Rodriguez-Carvajal, J. (1990). FULLPROF: a program for Rietveld refinement and pattern-matching analysis. Abstracts of the meeting Powder Diffraction, Toulouse, France, pp. 127–128.Google Scholar
Rodriguez-Carvajal, J. (1993). Recent advances in magnetic structure detemination by neutron powder diffraction. Physica B, 192, 55–69.Google Scholar
Scott, H. G. (1987). PROFIT – a peak-fitting program for powder diffraction profile. IUCr Satellite Meeting on X-ray Powder Diffractometry, 20–23 August 1987, Fremantle, Western Australia. Abstract P28.Google Scholar
Sivia, D. S. & David, W. I. F. (1994). A Bayesian approach to extracting structure-factor amplitudes from powder diffraction data. Acta Cryst. A50, 703–714.Google Scholar
Smrčok, Ľ., Mach, P. & Le Bail, A. (2013). Decafluorocyclohex-1-ene at 4.2 K – crystal structure and theoretical analysis of weak interactions. Acta Cryst. B69, 395–404.Google Scholar
Solovyov, L. A., Fedorov, A. S. & Kuzubov, A. A. (2014). Complete crystal structure of decafluorocyclohex-1-ene at 4.2 K from original neutron diffraction data. Acta Cryst. B70, 395–397.Google Scholar
Sonneveld, E. J. & Visser, J. W. (1975). Automatic collection of powder data from photographs. J. Appl. Cryst. 8, 1–7.Google Scholar
Toraya, H. (1986). Whole-powder-pattern fitting without reference to a structural model: application to X-ray powder diffraction data. J. Appl. Cryst. 19, 440–447.Google Scholar
Toraya, H. (1994). Applications of whole-powder-pattern fitting technique in materials characterization. Mater. Charact. Adv. X-ray Anal. 37, 37–47.Google Scholar
Waser, J. (1963). Least-squares refinement with subsidiary conditions. Acta Cryst. 16, 1091–1094.Google Scholar
Wiles, D. B. & Young, R. A. (1981). A new computer program for Rietveld analysis of X-ray powder diffraction patterns. J. Appl. Cryst. 14, 149–151.Google Scholar
Will, G. (1988). Crystal structure analysis and refinement using integrated intensities from accurate profile fits. Aust. J. Phys. 41, 283–296.Google Scholar
Will, G., Masciocchi, N., Parrish, W. & Hart, M. (1987). Refinement of simple crystal structures from synchrotron radiation powder diffraction data. J. Appl. Cryst. 20, 394–401.Google Scholar
Will, G., Parrish, W. & Huang, T. C. (1983). Crystal-structure refinement by profile fitting and least-squares analysis of powder diffractometer data. J. Appl. Cryst. 16, 611–622.Google Scholar
Young, R. A. (1993). The Rietveld Method. IUCr Monographs on Crystallography, Vol. 5. New York: Oxford University Press.Google Scholar

to end of page
to top of page