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

Section Pawley method

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


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