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, p. 282

Section Restrained cell

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: 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.


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
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
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
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
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
Le Bail, A. (2007). Structure Determination from Powder Diffraction Database – SDPD. .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
Pawley, G. S. (1981). Unit-cell refinement from powder diffraction scans. J. Appl. Cryst. 14, 357–361.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
Rietveld, H. M. (1969). A profile refinement method for nuclear and magnetic structures. J. Appl. Cryst. 2, 65–71.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
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
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
Waser, J. (1963). Least-squares refinement with subsidiary conditions. Acta Cryst. 16, 1091–1094.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

to end of page
to top of page