Rietveld refinement
Rietveld refinement is a technique described by Hugo Rietveld for use in the characterisation of crystalline materials. The neutron and X-ray diffraction of powder samples results in a pattern characterised by reflections (peaks in intensity) at certain positions. The height, width and position of these reflections can be used to determine many aspects of the material's structure. The Rietveld method uses a least squares approach to refine a theoretical line profile until it matches the measured profile. The introduction of this technique was a significant step forward in the diffraction analysis of powder samples as, unlike  other techniques at that time, it was able to deal reliably with strongly overlapping reflections. The method was first implemented in 1967, and reported in 1969 for the diffraction of monochromatic neutrons where the reflection-position is reported in terms of the Bragg angle, 2"θ". This terminology will be used here although the technique is equally applicable to alternative scales such as x-ray energy or neutron time-of-flight. The only wavelength and technique independent scale is in reciprocal space units or momentum transfer "Q", which is historically rarely used in powder diffraction but very common in all other  diffraction and optics techniques. The relation is Introduction. The most common powder X-ray diffraction (XRD) refinement technique used today is based on the method proposed in the 1960s by Hugo Rietveld. The Rietveld method fits a calculated profile (including all structural and instrumental parameters) to experimental data. It employs the non-linear least squares method, and requires the reasonable initial approximation of many free parameters, including peak shape, unit cell dimensions and coordinates of all atoms in the crystal structure. Other parameters can be guessed while still being reasonably refined. In this way one can refine the crystal structure  of a powder material from PXRD data. The successful outcome of the refinement is directly related to the quality of the data, the quality of the model (including initial approximations), and the experience of the user. The Rietveld method is an incredibly powerful technique which began a remarkable era for powder XRD and materials science in general. Powder XRD is at heart a very basic experimental technique with diverse applications and experimental options. Despite being slightly limited by the one-dimensionality of PXRD data and limited resolution, powder XRD's power is astonishing. It is possible to determine the accuracy  of a crystal structure model by fitting a profile to a 1D plot of observed intensity vs angle. It is important to remember that Rietveld refinement requires a crystal structure model and offers no way to come up with such a model on its own. However, it can be used to find structural details missing from a partial or complete ab initio structure solution, such as unit cell dimensions, phase quantities, crystallite sizes/shapes, atomic coordinates/bond lengths, micro strain in crystal lattice, texture, and vacancies. Powder diffraction profiles: peak positions and shapes. Before exploring Rietveld refinement, it is necessary  to establish a greater understanding of powder diffraction data and what information is encoded therein in order to establish a notion of how to create a model of a diffraction pattern, which is of course necessary in Rietveld refinement. A typical diffraction pattern can be described by the positions, shapes, and intensities of multiple Bragg reflections. Each of the three mentioned properties encodes some information relating to the crystal structure, the properties of the sample, and the properties of the instrumentation. Some of these contributions are shown in Table 1, below. The structure of a powder pattern is essentially defined  by instrumental parameters and two crystallographic parameters: unit cell dimensions, and atomic content and coordination. So, a powder pattern model can be constructed as follows: It is easy to model a powder pattern given the crystal structure of a material. The opposite, determining the crystal structure from a powder pattern, is much more complicated. A brief explanation of the process follows, though it is not the focus of this article. To determine structure from a powder diffraction pattern the following steps should be taken. First, Bragg peak positions and intensities should be found by fitting to a peak shape function  including background. Next, peak positions should be indexed and used to determine unit cell parameters, symmetry, and content. Third, peak intensities determine space group symmetry and atomic coordination. Finally, the model is used to refine all crystallographic and peak shape function parameters. To do this successfully, there is a requirement for excellent data which means good resolution, low background, and a large angular range. Peak shape functions. For general application of the Rietveld method, irrespective of the software used, the observed Bragg peaks in a powder diffraction pattern are best described by the so-called peak shape function (PSF). The  PSF is a convolution of three functions: the instrumental broadening formula_2, wavelength dispersion formula_3, and the specimen function formula_4, with the addition of a background function, formula_5. It is represented as follows: where formula_7 denotes convolution, which is defined for two functions formula_8 and formula_9 as an integral: The instrumental function depends on the location and geometry of the source, monochromator, and sample. Wavelength function accounts for the distribution of the wavelengths in the source, and varies with the nature of the source and monochromatizing technique. The specimen function depends on several things. First is dynamic scattering, and secondly the  physical properties of the sample such as crystallite size, and microstrain. A short aside: unlike the other contributions, those from the specimen function can be interesting in materials characterization. As such, average crystallite size, formula_11, and microstrain, formula_12, effects on Bragg peak broadening, formula_13 (in radians), can be described as follows, where formula_14 is a constant: Returning to the peak shape function, the goal is to correctly model the Bragg peaks which exist in the observed powder diffraction data. In the most general form, the intensity, formula_17, of the formula_18 point (formula_19, where formula_20 is the number of measured points ) is the sum of the contributions formula_21 from the m overlapped Bragg peaks (formula_22), and the background, formula_23, and can be described as follows: where formula_25 is the intensity of the formula_26 Bragg peak, and formula_27. Since formula_25 is a multiplier, it is possible to analyze the behaviour of different normalized peak functions formula_29 independently of peak intensity, under the condition that the integral over infinity of the PSF is unity. There are various functions that can be chosen to do this with varying degrees of complexity. The most basic functions used in this way to represent Bragg reflections are  the Gauss, and Lorentzian functions. Most commonly though, is the pseudo-Voigt function, a weighted sum of the former two (the full Voigt profile is a convolution of the two, but is computationally more demanding). The pseudo-Voigt profile is the most common and is the basis for most other PSF's. The pseudo-Voigt function can be represented as: where and are the Gaussian and Lorentzian contributions, respectively. Thus, where: The pseudo-Voigt function, like the Gaussian and Lorentz functions, is a centrosymmetric function, and as such does not model asymmetry. This can be problematic for non-ideal powder  XRD data, such as those collected at synchrotron radiation sources, which generally exhibit asymmetry due to the use of multiple focusing optics. The Finger–Cox–Jephcoat function is similar to the pseudo-Voigt, but has better handling of asymmetry 12, which is treated in terms of axial divergence. The function is a convolution of pseudo-Voigt with the intersection of the diffraction cone and a finite receiving slit length using two geometrical parameters, formula_57, and formula_58, where formula_59 and formula_60 are the sample and the detector slit dimensions in the direction parallel to the goniometer axis, and formula_61 is the  goniometer radius 12. Peak shape as described in Rietveld's paper. The shape of a powder diffraction reflection is influenced by the characteristics of the beam, the experimental arrangement, and the sample size and shape. In the case of monochromatic neutron sources the convolution of the various effects has been found to result in a reflex almost exactly Gaussian in shape. If this distribution is assumed then the contribution of a given reflection to the profile y at position 2"θ" is: where formula_63 is the full width at half peak height (full-width half-maximum), formula_64 is the center of  the reflex, and formula_65 is the calculated intensity of the reflex (determined from the structure factor, the Lorentz factor, and multiplicity of the reflection). At very low diffraction angles the reflections may acquire an asymmetry due to the vertical divergence of the beam. Rietveld used a semi-empirical correction factor, formula_66 to account for this asymmetry: where formula_68 is the asymmetry factor and formula_69 is +1,0, or –1 depending on the difference formula_70 being positive, zero or negative respectively. At a given position more than one diffraction peak may contribute to the profile. The intensity is simply the sum of  all reflections contributing at the point formula_71. Integrated intensity. For a Bragg peak formula_72, the observed integrated intensity, formula_73, as determined from numerical integration is where formula_75 is the total number of data points in the range of the Bragg peak. The integrated intensity depends on multiple factors, and can be expressed as the following product: where: Peak width as described in Rietveld's paper. The width of the diffraction peaks are found to broaden at higher Bragg angles. This angular dependency was originally represented by where formula_47, formula_48, and formula_49 are the half-width parameters and may be refined  during the fit. Preferred orientation. In powder samples there is a tendency for plate- or rod-like crystallites to align themselves along the axis of a cylindrical sample holder. In solid polycrystalline samples the production of the material may result in greater volume fraction of certain crystal orientations (commonly referred to as texture). In such cases the reflex intensities will vary from that predicted for a completely random distribution. Rietveld allowed for moderate cases of the former by introducing a correction factor: where formula_90 is the intensity expected for a random sample, formula_91 is the preferred orientation parameter and formula_92  is the acute angle between the scattering vector and the normal of the crystallites. Refinement. The principle of the Rietveld method is to minimize a function formula_93 which analyzes the difference between a calculated profile formula_94 and the observed data formula_95. Rietveld defined such an equation as: where formula_97 is the statistical weight and formula_98 is an overall scale factor such that formula_99. Least squares method. The fitting method used in Rietveld refinement is the non-linear least squares approach. A detailed derivation of non-linear least squares fitting will not be given here. Further detail can be found in  Chapter 6 of Pecharsky and Zavalij's text 12 . There are a few things to note however. First, non-linear least squares fitting has an iterative nature for which convergence may be difficult to achieve if the initial approximation is too far from correct, or when the minimized function is poorly defined. The latter occurs when correlated parameters are being refined at the same time, which may result in divergence and instability of the minimization. This iterative nature also means that convergence to a solution does not occur immediately for the method is not exact. Each iteration depends on the  results of the last which dictate the new set of parameters used for refinement. Thus, multiple refinement iterations are required to eventually converge to a possible solution. Rietveld method basics. Using non-linear least squares minimization, the following system is solved: where formula_101is the calculated intensity and formula_102is the observed intensity of a point formula_103in the powder pattern, formula_14, is a scale factor, and formula_20is the number of measured data points. The minimized function is given by: where formula_107is the weight, and formula_14from the previous equation is unity (since formula_14 is usually absorbed in the phase scale factor). The  summation extends to all formula_20 data points. Considering the peak shape functions and accounting for the overlapping of Bragg peaks because of the one-dimensionality of XRD data, the expanded form of the above equation for the case of a single phase measured with a single wavelength becomes: where: For a material that contains several phases (formula_120), the contribution from each is accounted for by modifying the above equation as follows: It can easily be seen from the above equations that experimentally minimizing the background, which holds no useful structural information, is paramount for a successful profile fitting. For a  low background, the functions are defined by contributions from the integrated intensities and peak shape parameters. But with a high background, the function being minimized depends on the adequacy of the background and not integrated intensities or peak shapes. Thus, a structure refinement cannot adequately yield structural information in the presence of a large background. It is also worth noting the increased complexity brought forth by the presence of multiple phases. Each additional phase adds to the fitting, more Bragg peaks, and another scale factor tied to corresponding structural parameters, and peak shape. Mathematically they are easily accounted for , but practically, due to the finite accuracy and limited resolution of experimental data, each new phase can lower the quality and stability of the refinement. It is advantageous to use single phase materials when interested in finding precise structural parameters of a material. However, since the scale factors of each phase are determined independently, Rietveld refinement of multi phase materials can quantitatively examine the mixing ratio of each phase in the material. Refinement parameters. Background. Generally, the background is calculated as a Chebyshev polynomial. In GSAS and GSAS-II they appear as follows. Again, background is treated as a Chebyshev  polynomial of the first kind ("Handbook of Mathematical Functions", M. Abramowitz and IA. Stegun, Ch. 22), with intensity given by: where formula_123are the coefficients of the Chebyshev polynomial taken from Table 22.3, pg. 795 of the Handbook. The coefficients have the form: and the values for formula_125are found in the Handbook. The angular range (formula_126) is converted to formula_52 to make the Chebyshev polynomial orthogonal by And, the orthogonal range for this function is –1 to +1. Other parameters. Now, given the considerations of background, peak shape functions, integrated intensity, and non-linear least squares minimization, the parameters used in  the Rietveld refinement which put these things together can be introduced. Below are the groups of independent least squares parameters generally refined in a Rietveld refinement. Each Rietveld refinement is unique and there is no prescribed sequence of parameters to include in a refinement. It is up to the user to determine and find the best sequence of parameters for refinement. It is worth noting that it is rarely possible to refine all relevant variables simultaneously from the beginning of a refinement, nor near the end since the least squares fitting will be destabilized or lead to a false minimum . It is important for the user to determine a stopping point for a given refinement. Given the complexity of Rietveld refinement it is important to have a clear grasp of the system being studied (sample, and instrumentation) to ensure that results are accurate, realistic, and meaningful. High data quality, a large enough formula_51 range, and a good model – to serve as the initial approximation in the least squares fitting – are necessary for a successful, reliable, and meaningful Rietveld refinement. Figures of merit. Since refinement depends on finding the best fit between a calculated and experimental pattern, it is important to  have a numerical figure of merit quantifying the quality of the fit. Below are the figures of merit generally used to characterize the quality of a refinement. They provide insight to how well the model fits the observed data. Profile residual (reliability factor): Weighted profile residual: Bragg residual: Expected profile residual: Goodness of fit: It is worth mentioning that all but one (formula_135) figure of merit include a contribution from the background. There are some concerns about the reliability of these figures, as well there is no threshold or accepted value which dictates what represents a good fit. The most  popular and conventional figure of merit used is the goodness of fit which should approach unity given a perfect fit, though this is rarely the case. In practice, the best way to assess quality is a visual analysis of the fit by plotting the difference between the observed and calculated data plotted on the same scale.