Rietveld refinement

Rietveld refinement is a technique devised 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 peaks in intensity at certain positions. The height, width and position of these peaks can be used to determine many aspects of the materials 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 reported for the diffraction of monochromatic neutrons where the peak 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

$Q\; =\; frac\; \{4\; pi\; sin\; left\; (\; heta\; ight\; )\}\{lambda\}$

Peak shape

The shape of a powder diffraction peak 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 peak almost exactly Gaussian in shape. If this distribution is assumed then the contribution of a given peak to the profile y_i at position 2"θ"_i is:

$y\_i\; =\; I\_k\; exp\; left\; [\; -4\; ln\; left\; (\; frac\{2\}\{H\_k^2\}\; ight\; )\; left\; (2\; heta\_i\; -\; 2\; heta\_k\; ight\; )^2\; ight\; ]$

where "H"_k is the full width at half peak height (full-width half-maximum), 2"θ"_k is the centreof the peak, and I_k is the calculated intensity of the peak (determined from the structure factor, the Lorentz factor, and multiplicity of the reflection)

At very low diffraction angles the peaks may acquire an asymmetry due to the vertical divergence of the beam. Rietveld used a semi-empirical correction factor, A_s to account for this asymmetry

$A\_s\; =\; 1\; -\; left\; [\; frac\; \{sP\; left\; (2\; heta\_i\; -\; 2\; heta\_k\; ight\; )^2\}\{tan\; heta\_k\}\; ight\; ]$

where P is the asymmetry factor and s is +1,0,-1 depending on the difference 2θ_i-2θ_kbeing 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 peaks contributing at the point 2θ_i.

Peak width

The width of the diffraction peaks are found to broaden at higher Bragg angles. This angular dependency was originallyrepresented by

$H\_k^2\; =\; U\; an^2\; heta\_k\; +\; V\; an\; heta\_k\; +\; W$

where U, V and W are the halfwidth 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 peak intensities will vary from that predicted for a completely random distribution. Rietveld allowed for moderate cases of the former by introducing a correction factor:

$I\_\{corr\}\; =\; I\_\{obs\}\; exp\; left\; (-Galpha^2\; ight\; )$

where I_obs is the intensity expected for a random sample, G is the preferred orientation parameter and α is the acute angle between the scattering vector and the normal of the crystallites.

Refinement

The principle of the Rietveld Method is to minimise a function M which represents the difference between a calculated profile y(calc) and the observed data y(obs). Rietveld defined such an equation as:

$M\; =\; sum\_\{i\}\; W\_i\; left\; \{\; y\_i^\{obs\}\; -\; frac\{1\}\{c\}\; y\_i^\{calc\}\; ight\; \}^2$

where W_i is the statistical weight and c is an overall scale factor such that $y^\{calc\}\; =\; c\; y^\{obs\}$

References

*cite journal
author=H. M. Rietveld
year=1969
title=A profile refinement method for nuclear and magnetic structures
journal=Journal of Applied Crystallography
volume=2
pages=65–71
doi=10.1107/S0021889869006558
*http://www.ccp14.ac.uk/ccp/web-mirrors/hugorietveld/xtal/paper2/paper2.html

Links

* [http://home.planet.nl/~rietv025/ The Rietveld Method] - the story.
* [http://www.ing.unitn.it/~maud/ The MAUD program] - software for Rietveld refinement.