Fundamental Parameters Line Profile Fitting in Laboratory Diffractometers

The fundamental parameters approach to line profile fitting uses physically based models to generate the line profile shapes. Fundamental parameters profile fitting (FPPF) has been used to synthesize and fit data from both parallel beam and divergent beam diffractometers. The refined parameters are determined by the diffractometer configuration. In a divergent beam diffractometer these include the angular aperture of the divergence slit, the width and axial length of the receiving slit, the angular apertures of the axial Soller slits, the length and projected width of the x-ray source, the absorption coefficient and axial length of the sample. In a parallel beam system the principal parameters are the angular aperture of the equatorial analyser/Soller slits and the angular apertures of the axial Soller slits. The presence of a monochromator in the beam path is normally accommodated by modifying the wavelength spectrum and/or by changing one or more of the axial divergence parameters. Flat analyzer crystals have been incorporated into FPPF as a Lorentzian shaped angular acceptance function. One of the intrinsic benefits of the fundamental parameters approach is its adaptability any laboratory diffractometer. Good fits can normally be obtained over the whole 20 range without refinement using the known properties of the diffractometer, such as the slit sizes and diffractometer radius, and emission profile.

The fundamental parameters approach to line profile fitting uses physically based models to generate the line profile shapes. Fundamental parameters profile fitting (FPPF) has been used to synthesize and fit data from both parallel beam and divergent beam diffractometers. The refined parameters are determined by the diffractometer configuration. In a divergent beam diffractometer these include the angular aperture of the divergence slit, the width and axial length of the receiving slit, the angular apertures of the axial Soller slits, the length and projected width of the x-ray source, the absorption coefficient and axial length of the sample. In a parallel beam system the principal parameters are the angular aperture of the equatorial analyser/Soller slits and the angular apertures of the axial Soller slits. The presence of a monochromator in the beam path is normally accommodated by modifying the wavelength spectrum and/or by changing one or more of the axial divergence parameters. Flat analyzer crystals have been incorporated into FPPF as a Lorentzian shaped angular acceptance function. One of the intrinsic benefits of the fundamental parameters approach is its adaptability any laboratory diffractometer. Good fits can normally be obtained over the whole 20 range without refinement using the known properties of the diffractometer, such as the slit sizes and diffractometer radius, and emission profile. cial software packages [5,6,7] have become available for fundamental parameters profile fitting (FPPF) either for use as single line profile fitting, lattice parameter refinement or for Rietveld analysis.
FPPF has been used to synthesise and fit data from both parallel beam and divergent beam diffractometers. The refined parameters are determined by the diffractometer configuration. In a divergent beam diffractometer these include the angular aperture of the divergence slit, the width and axial length of the receiving slit, the angular apertures of the axial Soller slits, the length and projected width of the x-ray source, the absorption coefficient and axial length of the sample. In a parallel beam system the principal parameters are the angular aperture of the equatorial analyser/Soller slits and the angular apertures of the axial Soller slits. The presence of a monochromator in the beam path is normally accommodated by modifying the wavelength spectrum and/or by changing one or more of the axial divergence parameters. Flat analyser crystals have been incorporated into FPPF as a Lorentzian shaped angular acceptance function.
One of the intrinsic benefits of the fundamental parameters approach is its adaptability to any laboratory diffractometer. Good fits can normally be obtained over the whole 2θ range without refinement using the known properties of the diffractometer, such as the slit sizes and diffractometer radius, and the emission profile. Fine tuning is sometimes necessary to accommodate a monochromator or to compensate for the fact that certain aberrations are not completely independent [8]. Under these conditions some of the instrument parameters need to be refined, but the refined values normally are within ±10 % of the actual values. Correlation between refined instrument parameters can occur when fitting to data over a restricted 2θ range. Such correlation occurs between the axial divergence parameters and absorption as both of these aberrations can produce similar forms of asymmetric profiles between 2θ = 50° and 100° in diverging beam diffractometers. Correlation is minimised by using data with a large 2θ range so that the unique angular dependence of individual aberrations becomes evident. When a set of instrument profiles cannot be fitted by FPPF, this is usually an indication that either the model used is invalid (eg. incorrectly chosen slit value), the instrument is mis-aligned, there is an overlapping impurity line or, the specimen is generating crystallite size broadening or is inhomogeneously strained.
FPPF was designed originally as a tool for analysing diffraction line broadening. Fitting is done by convolution and corrections for instrument broadening and peak shift are intrinsic to the refinement. When an instrument is well characterised, line broadening can be analysed without a reference specimen. Moreover, when a reference standard is used, which has different properties from the specimen with line broadening, some compensation can be made for these differences. For example, when LaB 6 SRM 660a (µ powder ≈ 500 cm -1 ) is used as a reference, compensation can be made for differences in the absorption of the sample and LaB 6 . In the latest version of the commercial software package (TOPAS) 2 , the concept of fundamental parameters has been extended so that any user defined profile that accurately describes the physical broadening can be readily convoluted into the refinement.
In this paper we will discuss the physical origin of the instrumental profile shapes for various laboratory diffractometer configurations including both divergent beam and parallel beam instruments. This will include a description of the geometrical aberrations as well as discussion on the nature of the wavelength distribution and the influence of monochromators on this distribution. Some discussion is also presented to demonstrate how the FPPF may be fitted to experimental data from a material with a low attenuation coefficient.

Basic Objectives of the FPPF Technique
One of the basic objectives of the FPPF technique is to be able to fit any powder diffraction profile using a physically based model to describe both the instrument profile and any diffraction broadening generated by the specimen. In principle, therefore, the technique should be adaptable to any powder diffractometer and fit profiles of widely differing shapes, such as those in Fig. 2, by simply modifying the physical parameters of the diffractometer used to describe the profile.
Although most applications of FPPF have focussed on the conventional diffractometer it has also been utilised for analysing neutron diffraction data [9] and synchrotron data [10,11]. In the TOPAS implementation of FPPF there are a wide variety of possible aberration functions available within the package and these can put together to suit a particular diffractometer design and in terms of parameters that are relevant to the instrument. One of the most important achievements of the FPPF technique for practical users is speed of calculation. Accurate multiple convolution calculations over large 2θ ranges can be very time intensive and it is of central importance to minimise this time to enable Rietveld refinement to be completed in "seconds" without loss of accuracy within the profile function synthesis procedures. Various procedures have been implemented, some of which are described by Cheary and Coelho [2,3,4], but one of the most important has been to code the time intensive calculations at an assembler code level taking steps to optimise the use of the various registers within the PC chip.
years the number of diffractometer options available from manufacturers has increased and users operate with a wider range of x-ray optical designs. The types of geometrical aberrations encountered is somewhat broader than those discussed in the classic work by Wilson [12] as will be discussed below.

Divergent Beam Diffractometers-Symmetric Diffraction
The most widely used laboratory diffractometer in use today is still the divergent beam diffractometer with either a bent graphite monochromator in the diffracted beam (see Fig. 3a) or, a ground and bent asymmetrically cut germanium monochromator in the incident beam (see Fig. 3b). Both of these configurations possess a similar array of geometrical aberrations. The major difference between them is the wavelength distribution which normally consists of both the Kα 1 and Kα 2 components of the K spectrum in the graphite monochromator case and only the Kα 1 with the Ge monochromator [13]. Further discussion of the wavelength distribution and the effects of monochromators is given later.
The principal geometric aberrations contributing to profiles from the above diffractometers are, (i) the finite width of the x-ray source, (ii) a divergent incident beam on to a flat specimen (flat specimen error), (iii) the finite width of the receiving slit, (iv) the beam penetration into the specimen (specimen transparency), (v) the deviation of the beam from the equatorial plane (axial divergence). These aberrations all produce some degree of line broadening and, in the case of flat specimen error, specimen transparency and axial divergence, some asymmetry is also introduced. Zero 2θ and specimen surface displacement errors may also be present in a diffractometer, but these only affect the 2θ position of a profile and not its shape. In both configurations the monochromators not only determine the wavelength distribution, but they also act to reduce the axial divergence and it is often considered unnecessary to include Soller slits between the sample and the monochromator. necessary to operate under asymmetric conditions as illustrated in Fig. 4a. The problem with operating in this mode on most commercial diffractometers is that the receiving slit is no longer at the focus of the diffracted beam and profiles are broadened by "defocussing". The amount of defocussing is determined by the angle of divergence of the incident beam and distance of the focus from the receiving slit.
Defocussing also occurs in diffractometer configurations where, • the sample is oscillated ±δω about the diffractometer axis so that the angle of incidence on to the specimen varies between θ + δω and θδω. This oscillation moves the focus of the diffracted beam continuously back and forth in front of and behind the receiving slit, • the receiving slit is replaced by a position sensitive detector (PSD). The only position on the detector that is normally in focus is its centre (see Fig. 4b); all diffracted beams entering the detector at off-centre positions are defocussed. In PSD systems, the aberrations contributing to a profile include all the standard diffractometer aberrations, except the aberrations formerly due to the receiving slit are replaced by defocussing, the discharge resolution of the detector, and parallax error [14,15]. In a scanning PSD diffractometer, the recorded profile shape is an average of all the profile shapes across the active window length.

Parallel Beam Diffractometers
There are two common forms of the parallel beam powder diffractometer which are illustrated in Fig. 5. These are based on using either analyser slits (otherwise referred to as equatorial Soller slits) or a flat Ge/Si analyser crystal as the angular discriminator of the diffracted beam. Amongst laboratory diffractometers the analyser slit set-up is the most widely used form as it offers more intensity but poorer resolution than the analyser crystal set-up.
Parallel beam diffractometers have emerged as one of the most popular forms of laboratory diffractometer over the past ten years and now constitute more than 30 % of the new diffractometer purchases. There are fewer geometric aberrations contributing to the profile and systematic errors arising from specimen displacement, specimen transparency and surface roughness are not significant. There are two geometric aberrations contributing to a parallel diffractometer, (i) the angular acceptance function of the analyser foils or analysing crystal, (ii) deviation of the beam from the equatorial plane (ie. axial divergence). In most laboratory diffractometers, the parallel beam is produced by using a parabolic graded multilayer mirror with the line x-ray source positioned at the focus of the mirror [16]. Although the beam may be parallel in the equatorial plane, it will not be parallel in axial plane and axial divergence can be expected in both the incident and diffracted beams. Low angle profiles will therefore be asymmetric although not to the same extent as diverging beam instruments.

The Instrument Aberrations
The geometric instrument aberrations tend to determine the shape of a diffractometer profile at low 2θ angles (ie., 2θ < 50°). At high 2θ angles (2θ > 100°), the profile conforms primarily to the shape of the wavelength distribution in the beam. With the exception of the aberrations associated with "receiving system" of the diffractometer and the x-ray source, all of the geometric instrument aberration profiles vary with 2θ. In the following sections the shapes of the major instrument aberrations used in FPPF analysis to describe the various laboratory diffractometer configurations are discussed for conditions that are typical of those encountered in practice. The aberration functions generated by mis-setting a diverging beam diffractometer or using it under asymmetric conditions are also discussed. Most of the results quoted here are for a diffractometer radius R = 215 mm.
The convention adopted here for describing the angular variables is that 2φ refers to the continuously variable angle measured on the diffractometer whereas 2θ or 2θ B refers to the Bragg angle of the diffraction line. The angle ε refers to the difference between the measured angle 2φ and 2θ, ε = 2φ -2θ. (4)

Finite X-Ray Source Width
The profile shape of this aberration is generally expressed as an impulse function of width ∆2θ x as shown Fig. 6a. Although the choice of an impulse function may be not be strictly valid for describing the x-ray  source aberration, the exact shape used is not critical when a long fine focus tube (target width ≈0.4 mm) is installed on the diffractometer. At a take-off angle of 6°t his appears as a projected width w x ≈ 0.04 mm and the aberration profile has a width ∆2θ x ≈ 0.01° and does not contribute significantly to the overall width of the instrument profile.
In broad focus tubes, the target width ≈2 mm and w x ≈ 0.2 mm at 6° take-off so that the aberration profile width ∆2θ x ≈ 0.056°. At this level the source width makes a much bigger contribution to the overall width of the instrument profile and a more accurate form for the aberration profile shape is necessary. A good approximation under these circumstances is a Gaussian shape rather than an impulse function. In diffractometers with curved crystal incident beam monochromators the source width can also have a greater contribution because of the magnification effect introduced by the monochromator. This occurs with asymmetrically-cut Johansson incident beam monochromators where the source-crystal and crystal-focal point distances are typically ≈120 mm and ≈230 mm, respectively. A fine focus tube with a projected width of 0.04 mm is then effectively magnified to ≈0.08 mm. Under these conditions the effective source width can be trimmed down by reducing the width of the focal line slit.
For accurate line profile analysis it is also necessary to modify the simple impulse model even with long fine focus tubes. Bergmann [17] has shown that most of the anode surface in an x-ray tube produces x rays albeit at a much lower intensity than the focal line on the anode. This is illustrated in Fig. 7 which shows the intensity recorded by scanning with a 50 µm slit across the image of a x-ray source formed through a 10 µm pinhole in platinum. A better approximation to the aberration function is a sharp impulse function superimposed on a broad impulse function to represent the so called "tube tails". This is illustrated in Fig. 6b. The parameters introduced to describe the "tube tails" are the extents of the high and low angle tails, Z 1 and Z 2 , and the intensity of the tail f relative to the intensity at the tube focus. In most instances the intensity of the tails is ≈0.1 % of the peak intensity and is only significant when analysing intense lines. The tails themselves are not necessarily symmetric with respect to the tube focus and can extend over a 2θ range up to 0.6°.

X-Ray Receiving System Models
In diverging beam diffractometers the receiving slit is placed at the focus of the diffracted beam and for per-fect focussing should have an infinitely small width. Owing to the many aberrations present, focussing is never perfect and the count rate incident on the receiving slit tends to increase with increasing slit width, but at the expense of resolution. In parallel beam diffractometers the receiving system is based on using either the Hart-Parrish system of analyser slits [18] or, a flat analyser system as illustrated in Fig. 5 earlier. In many glancing incidence diffractometers the receiving system consists of analyser slits and an analyser crystal in the diffracted beam. Although the aberration functions associated with the various receiving systems for parallel beam and divergent beam diffractometers all possess different shapes, they all possess the common property of being independent of 2θ.

Receiving Slit in a Diverging Beam Diffractometer
Most commercial diffractometers have a selection of receiving slits ranging in width from 0.05 mm up to 0.3 mm although occasionally larger slit sizes up to 0.6 mm are used to measure integrated intensity rapidly. The aberration function for a perfectly aligned receiving slit is an impulse function of width ∆2θ r given by, (5) where w r is the width of the receiving slit. The angular width ∆2θ r subtended by the receiving slit relative to the diffractometer axis is therefore normally between 0.013° (0.05 mm) and 0.08° (0.3 mm). When the slit size is larger than 0.15 mm, the receiving slit aberration is often the dominant aberration in a diffractometer over the angular range 2θ = 15° to 60°.

Parrish-Hart Analyser Slits
Analyser slits act as an angular filter in the diffracted beam. The aberration function or transmission function for these slits is a triangle function, as shown in Fig. 8, in which the base width ∆2θ r is given by the angular aperture ∆ of the slits.
In the original Parrish-Hart diffractometer on Station 2.3 at the Daresbury synchrotron, the analyser slits were 360 mm with a spacing of 0.2 mm between adjacent foils giving an angular 2θ aperture ∆ ≈ 0.06°. In laboratory diffractometers the angular aperture ∆ is typically ≈0.1°. A problem often encountered with analyser slits is specular x-ray reflection from the analyser foils [19]. Weak satellite peaks appear on both the high angle and low angle profile tails but not necessarily of the same intensity as shown in Fig.9a. This effect has been incorporated into the aberration profile by adding two Voigt functions of unequal intensity ,one on each side of the triangular aberration function, to represent the satellite reflections [11]. The parameters of the satellite peaks can then be determined by fitting profiles from a reference material such as the NIST reference standard LaB 6 , SRM 660a. An alternative approach to determining the aberration profile of an analyser slit system, without the effects of the wavelength profile distorting the result, is to simply carry out a 2θ scan across the incident beam. Provided the axial divergence of the incident beam is kept small, by including axial Soller slits, and the equatorial divergence is negligible then the incident beam scan will have exactly the same shape as the aberration profile of the analyser slits.

Analyser Crystals
The inclusion of an analyser crystal in the diffracted beam of a parallel beam diffractometer gives high res-  olution diffraction patterns with a low background, but the intensity is invariably less than the Parrish-Hart configuration. The aberration profile introduced by the analyser crystal is generally very narrow and can be determined by measuring the rocking curve of the crystal. For a perfect analyser crystal the aberration profile will be determined by the Darwin profile of the analyser crystal. In practice, however, the aberration profile will be broadened by the mosaic structure of the crystal, any stresses in the crystal and any waviness or curvature of the crystal surface [20]. As a consequence the aberration profile can be dependent on the size of the beam incident on the crystal. A first approximation to the shape of the aberration profile of an in-situ analyser can be obtained from a 2θ scan of the analyser/detector using a very fine incident beam as shown in Fig. 9b. Although the profile recorded in this way also has the wavelength distribution folded into it, the result does at least give an indication of the shape and upper limit of the FWHM of the aberration profile.

Flat Specimen Error
The basic optics of the focussing powder diffractometer set up for symmetric diffraction is illustrated in Fig. 10. The x rays are incident at an angle θ on an ideal polycrystalline specimen with a surface radius of curvature ρ. For diffraction from a particular hkl plane the common property of all the diffracted rays from the specimen is that they all deviate through the same angle 2θ. By simple geometry it can be shown that all the diffracted rays converge to a focus on a circle which has the same curvature as the specimen surface. The focus of the diffracted rays defines the position of the receiv-ing slit. In commercial diffractometers the specimen is invariably flat and the diffracted beam no longer focusses perfectly. Good focussing characteristics, however, can be maintained with appropriately chosen slits to limit the equatorial divergence and a diffractometer radius R sufficiently large to reduce defocussing errors to an acceptable level without losing too much diffracted intensity.
For an incident beam, with an equatorial divergence α centred on the diffractometer axis, the aberration profile J FS (2θ) is asymmetric and exists only for the region ε = 0. X rays diffracted from the centre of the specimen are detected at 2φ = 2θ where as x rays diffracted offcentre are detected at 2φ < 2θ as shown in Fig. 11.
The relationship between the difference ε = 2φ -2θ and the distance q from the diffractometer axis at which the ray is diffracted is, (6) assuming small angles of divergence for the incident beam (eg. α < 2°). Thus, when a specimen is illuminated over a length Q then the x rays diffracted from either extremity of the beam (ie., at q = ±Q/2) defines the limiting value of ε M of the aberration function J FS (ε) which exists within the range 0 ≥ ε ≥ ε M where (7) This relation can also be expressed in terms of the equatorial divergence of the beam α starting from the relation,  . Diffraction from a flat plate showing the relationship between the measured angle 2φ on the diffractometer and the diffraction angle 2θ for a ray at the outer limit of a beam of divergence α.
When 2θ > 10°, Eq. (8) can be approximated as and ε M becomes (9) When the incident beam is centred on the diffractometer axis the normalised equation for the aberration function J FS (ε) for flat specimen error is, Modern commercial diffractometers operate with either a fixed divergence α or a fixed illumination length Q.
In the fixed α mode the aberration function is broad at low 2θ and ε M has a cot θ dependence whereas in fixed Q mode the breadth rises from zero at 2θ = 0 up to a maximum at 2θ = 90°. The extent of the changes in the aberration function J SF (ε) for each mode of operation using typical operating values for α and Q are shown in Fig. 12.
With a fixed angle of divergence α = 1°, the effects of flat specimen error in commercial diffractometers are discernable as an increase in both the asymmetry and breadth below 2θ ≈ 40° [21]. Under conditions of constant α the beam size increases with decreasing 2θ and eventually the beam will cover the whole specimen. The angle 2θ lim at which this occurs is given by, (11) where L sp is the length of the specimen. Under these circumstances the value of ε M is given by Eq. (7) with Q = L sp . Below 2θ lim the aberration function remains the same as the beam extends beyond the specimen.
In diffractometers with a fixed illuminated length, the effects of flat specimen error are generally smaller at low 2θ values and increase with increasing 2θ. With an illuminated length of 20 mm on the specimen, flat specimen error is clearly discernable at 2θ > 20°. It is not always possible to collect data from fixed Q mode diffractometer over a large 2θ range (eg., up to 150° 2θ as the α angle required at large 2θ is larger than the diffractometer can accommodate. For example, to maintain a fixed beam length of 20 mm on a specimen over the range 2θ = 0° to 90°, the angle of divergence α will need to increase from 0° up to ≈4°. A divergence angle α = 4° is close to the maximum value at which most diffractometers can operate when a pyrolytic graphite monochromator is installed in the diffracted beam. At angles of α greater than 4° the diffracted beam may extend beyond the graphite crystal and not be diffracted into the detector. In practice, this is not usually a problem as fixed Q mode diffractometry is normally used for analysing materials such as clays with very low angle diffraction lines where the lines of interest start at 2θ ≈ 3°.

Specimen Transparency
Specimen transparency produces asymmetry and broadening of the instrumental profile function. For perfect focussing all the diffraction should occur on the focussing circle, but when the beam penetrates the surface, diffraction will occur over a range of depths within the specimen. The aberration function J µ (ε) for an infinitely thick specimen is given by, Comparison of flat specimen aberration functions at 2θ = 20° and 2θ =60° for a diffractometer with either a fixed divergence angle α = 1° or fixed illumination length Q = 20 mm (R = 215 mm). All the aberration functions shown were convoluted with a very narrow Lorentzian emission profile with a FWHM = 0.001 mÅ to overcome infinities at ε = 0. Each plot is normalised to the same I max and both plots cover a range of 0.2° 2θ.
where δ = (2/µR) sin2θ rad and µ is the linear attenuation coefficient. In low absorbing specimens which cannot be considered to be infinitely thick, the angular variable ε has a lower limit ε min so that in Eq. (12), where T is the specimen thickness. As a consequence the unit area normalising constant in (12) also changes and J µ (ε) becomes, The asymmetry and broadening from specimen transparency is greatest for low absorption materials and is clearly evident when the linear attenuation coefficient µ < 50 cm -1 . The contribution of specimen transparency is greatest at 2θ ≈ 90° and at this angle the aberration profile has a FWHM ≈ 0.03° 2θ when µ ≈ 50 cm -1 , but this drops to ~0.005° 2θ when µ ≈ 200 cm -1 . Fig. 13 shows the shapes of aberration profiles for an infinitely thick specimen and how they are affected by both µ and 2θ. Specimen transparency effects are quite strong in polymeric materials with very low attenuation coefficients (µ ≈ 30 cm -1 or less). They can also show up in loosely bound powders of low atomic number materials (e.g., MgO or Si) where the porosity ≈50 % or less so that a material with a µ = 100 cm -1 is reduced to a powder with a µ = 50 cm -1 .

Diffractometer Defocussing
Defocussing results in broadened diffraction lines and occurs when the receiving slit is not positioned at the focus of the diffracted beam. The most common causes of defocussing are, • mis-setting the incident beam angle ω so that it is no longer at the symmetric condition ω = θ as in asymmetric diffraction or, • wrongly positioning the receiving slit so that the distance of the slit from the sample is larger or smaller than the nominal radius R of the diffractometer. For both of these conditions the focus of the diffracted beam will be either in front of or behind the receiving slit as illustrated in Fig. 14 for the asymmetric diffraction case.
The aberration profile J DF (ε) for this condition is a impulse function of angular width ∆ DF = D/R (rad) where D is the width of the defocussed beam at the receiving slit. Assuming the equatorial divergence α is small (ie., 2° or smaller) the angular width ∆ DF of the aberration profile is given by (13) Under asymmetric diffraction conditions R 2 and R are related by the equation, (14) so that the width ∆ DF of the aberration function is, (15) In asymmetric diffraction, the breadth of diffraction lines increases as the deviation from the symmetric condition, |ωθ |, increases. Defocussing is larger at low 2θ angles and varies more rapidly with ωθ at low 2θ values. Conversely, at high 2θ values, diffractometers will tolerate quite large errors in ωθ (ie., up to ±5° at 2θ ≈ 150°) without the effects of defocussing being detectable in the line breadth. Also, by reducing the angle of divergence α the defocussing ∆ DF can also be reduced, but at the expense of the diffracted intensity. The plot of ∆ DF vs (ωθ) for a range of 2θ angles from 30° to 150° is shown in Fig. 15. and cos rad 13. Specimen transparency aberration profiles at 2θ = 30° and 90° for linear attenuation coefficients of 50, 100, and 200 cm

Axial Divergence
In a laboratory diffractometer only a small fraction of the photons that form the incident beam emerge from an x-ray source parallel to the equatorial plane. Likewise most of the x rays that reach the detector slit after diffracting from the sample are angled to the equatorial plane. Under these circumstances, a diffractometer will record x-ray counts over a range of angles 2φ other than the true diffraction angle 2θ. The only rays for which 2φ = 2θ will be those propagating parallel to the equatorial plane and incident on the diffractometer axis. In practice, axial divergence is most readily recognised by the asymmetry it introduces into low angle diffraction lines (2θ < 30°) where the low angle tails extend further than the high angles tails.
For a particular ray path the measured diffraction angle 2φ for a true diffraction angle 2θ depends on the axial divergence β and γ in the incident and diffracted rays (see Fig.16). Assuming small angles for β and γ, the difference ε = 2φ -2θ is given by,  From this equation it is evident that the effect of axial divergence on ε will not only be large when 2θ is small, but also when 2θ is large (ie., 2θ > 150°) where the asymmetry is opposite to that at low 2θ. Although axial divergence is as strong in very high angle lines as it is in low angle lines, the positive asymmetry developed tends to be less noticeable as it is overshadowed by the dispersion of the emission profile. When axial divergence in the incident beam is small (β ≈ 0) as in many synchrotron systems or laboratory diffractometers with very narrow incident beam Soller slits, all the axial divergence arises from the diffracted beam and Eq. (16) reduces to, For this condition, axial divergence is absent at 2θ = 90°, ε ≤ 0 when 2θ < 90° and ε ≥ 0 when 2θ > 90°. The aberration function J AX (ε) for this condition can be derived analytically [22,23,3,4] by considering an axially parallel beam from a line source incident on a narrow capillary specimen of randomly oriented crystallites as shown below in Fig. 17.
The emergent diffracted beam from each individual ray in the incident beam in Fig.18 is a collection of radiating cones with a semi-angle 2θ and J AX (ε) can be derived by determining the intensity profile when the receiving slit is scanned across these cones. The parameters of the diffractometer that define J AX (ε) in this instance are L s , the axial length of the specimen bathed in x rays and the length L r of the receiving slit and the function J AX (ε) is given by, where Examples of this aberration function at 2θ = 10° and 50° are shown in Fig. 18 for conditions that resemble a laboratory diffractometer with very narrow incident beam Soller slits (ie., ≤1°), R = 215 mm, an illuminated specimen length Ls = 12 mm and a receiving slit of length Lr = 16 mm.
When allowance is made for axial divergence in both the incident and diffracted beams and for the presence of Soller slits in each of these beams, the calculation of the aberration function J AX (ε) can no longer be done analytically. Eastabrook [24] first demonstrated how to calculate the axial divergence aberration profile for a conventional diffractometer, but restricted his discussion to instrumental conditions that could be solved analytically rather than for conditions that were widely used in practice. Pike [25] generalised the application conditions and included the effect of Soller slits, but restricted his calculations to the determination of the centre of gravity and variance of profiles rather than the aberration profile itself. More recently Cheary and Coelho [3,4] have developed a semi-analytical approach to the calculation and their results have been incorporated into a profile refinement procedure. In short, their procedure consists of, (i) calculating the analytical aberration function J AX (β,ε) arising from incident rays all with the same axial divergence β. The instrument parameters required to define this function are the axial lengths L x , L s , and L r of the x-ray source, the sample and the receiving slit, Arcs formed in the vicinity of the receiving slit by the diffraction cones radiating from all points across the sample axis. The axial divergence aberration function is the profile shape detected by an infinitely narrow receiving slit as it scans through the arcs. 2 1 (ii) incorporating the Soller slits into the calculation as angular intensity filters on the axial divergence β and γ in the incident and diffracted beams, respectively. The transmission functions, S I (β) and S D (γ), for the incident and diffracted beam Soller slits, respectively, are each triangle functions with 100 % transmission at β = 0 and γ = 0, and 0 % transmission at β = ± ∆ I /2 and γ = ± ∆ D /2 where ∆ I and ∆ D are the angular apertures of the slits (see Fig. 9 for definition of ∆), (iii) calculate the full aberration profile J AX (ε) by integrating the aberration functions determined at each β across all the allowed β values, In practice this integration is converted to a summation, but great care is needed carrying out the summation because of singularities in the J AX (β,ε) functions. Considerable savings in computing time without loss of accuracy can be achieved by using a smoothing procedure in which the function J AX (ε) is convoluted in ε space with an impulse of width = step size of the data being fitted. The saving in computing time comes from a reduction in the number of summation points at discrete values of β required to accurately reproduce Eq. (18) using a summation. In the absence of Soller slits the maximum axial divergence in the incident and diffracted beams, , and respectively, is determined by the length L x , L s , and L r and given by, Under these circumstances the axial divergence can be as large as 10° in the incident beam and 10° in the diffracted beam for a diffractometer with a 12 mm long xray source, a 12 mm receiving slit and a 20 mm wide specimen. However, most commercial diffractometers are supplied with Soller slits having angular apertures ∆ between 2° and 5°, and in such cases the breadth of the aberration profile is reduced quite dramatically. This is illustrated in Fig. 19 which shows aberration profile shapes calculated using the method of Cheary and Coelho for diffractometer configurations that include Soller slits as well as a configuration with no Soller slits. When narrow Soller slits are included in the beam path (ie., ∆ I and ∆ D ≈ 2°), the dimensions L x , L r , and L s tend to be redundant parameters in the calculation of J AX (ε) as the Soller slits control the maximum axial divergence in the incident and diffracted beams.
The major change brought on by axial divergence is extension of the low angle tails of profiles below 2θ ≈ 50°. At very high angles, 2θ ≥ 150°, the asymmetry reverses and the extension of the high angle tails increases with increasing 2θ. Broadening from axial divergence is evident at all 2θ values but passes through a minimum in the region 2θ ≈ 110°. However, the shift in peak angle 2θ max at I max relative to the true 2θ B is small at all 2θ values unlike the shift in centre of gravity 2θ cg -2θ B which varies quite considerably from very large negative values, near 2θ = 0°, to very large positive values as 2θ approaches180° [12]. In practice it is not always possible to calculate the exact form of the axial divergence aberration function for a particular specimen/diffractometer configuration. The two main reasons for this are, (i) in specimens with strong preferred orientation, such as thin films and rolled or extruded metals, the diffraction cones are no longer of uniform intensity along the arcs of the diffraction cones. (ii) the inclusion of a monochromator in the beam path reduces the axial divergence. When either diffracted beam or an incident beam monochromator is present the optical path length of the beam is extended. With a diffracted beam monochromator the optical path length of the diffracted beam is extended by the optical path between the receiving slit and the detector slit. In a diffractometer with a graphite monochromator tuned to Cu Kα radiation the optical path length in the diffracted beam is increased by ≈100 mm giving an effective radius of the detector arm of 100 + R mm. An example of the extent to which axial divergence is reduced by introducing a graphite diffracted beam monochromator is illustrated in Fig. 20. Incident beam monochromators also increase the optical path length of the incident beam. For the most common type of monochromator in use, the asymmetrically Ge ground and bent monochromator, the path is increased by over 300 mm. Monochromators also act as angular intensity filters and their effect on a profile is similar to the addition of Soller slits in the beam path [26]. The effect of a monochromator can therefore be represented as a Soller slit in a profile fitting model.

Linear Position Sensitive Detector (LPSD) Aberrations
LPS detectors with angular windows up to 10° (at 200 mm) are used in reflection mode on commercial diffractometers to increase the data collection rate particularly for kinetic based studies. In most cases LPS detectors are used in stationary mode and only a fixed angular region of a pattern is recorded. A number of manufacturers offer LPSD systems where the detector can be scanned in θ -2θ mode with its centre maintained at the focus condition. As the detector is scanned the total diffraction pattern is formed by adding together and averaging the patterns recorded at each step. In  this mode a diffraction pattern can be accumulated very rapidly with excellent counting statistics because individual lines are within the detector window for a considerable time. For example a LPSD with an angular window of 10° moving at 5° per min corresponds to a diffraction peak being detected for 120 s. Consequently, in this mode, a 100° 2θ diffraction pattern may only take approximately 20 min to collect. To obtain the same level of counting precision in a conventional single slit diffractometer with 0.02° 2θ steps would take approximately one week. In practice, this gain in time is never fully realised because of the longer dead times of LPSDs and the damage that can be caused on the anode wire by very high localised count rates. Also the diffraction peaks recorded at off-centre positions along the detector window are broadened and asymmetric. Some manufacturers make the angular window a software adjustable parameter to maintain the resolution although this means that only a fraction of the detector is being used.
In LPSD systems, the receiving slit aberration is no longer relevant and the flat specimen aberration is replaced by an aberration function that embodies three effects which are folded together in the final function; • flat specimen error including defocussing, • parallax error, • thermal noise. A full treatment of these aberrations and how they are modified by a scanning LPSD system is given in Cheary and Coelho [14]. For a stationary LPSD the aberrations are discussed in Secs. 4.7.1-4.7.3.

LPSD Flat Specimen Error J PSD (ε ε) Including Defocussing
In a conventional diffractometer defocussing and flat specimen error can be convoluted together independ-ently, but this is no longer valid in a LPSD system. The effect of both the specimen and the detector being flat is that the defocussing is no longer symmetric about the centre of the LPSD. In a stationary LPSD the aberration function J PSD (ε) depends on the 2θ value of a peak and on the offset angle β of the recorded peak from the centre of the LPSD. The angular variables used to define the LPSD system are given in Fig. 21 where the ζ is limited by the angle of divergence α of the incident beam (ie., -α/2 ≤ ζ ≤ +α/2).
When the diffractometer is operated symmetrically and the incident beam is centred on the diffractometer axis then assuming small angles for ε and ζ, the difference ε = 2φ -2θ is related parabolically to ζ, (20) where The equation for the aberration function J PSD (ε) is obtained by transforming the intensity across the incident beam IB(ζ) from ζ space into ε space, ie., (21) As the incident beam intensity IB(ζ) is reasonably constant then the un-normalised form of J PSD (ε) = |dζ /dε |. This is readily calculated by differentiating the transformation Eq. (20) and has the form,  θ β π θ β β θ β θ β ε ζ where ε ≤ ε 0 . When calculating this function allowance has to be made for the fact that Eq. (20) can be double valued. Physically this means that more than one part of the incident beam contributes at a particular ε value and J PSD (ε) possesses a discontinuity at the boundary between two rays and one ray contributing to the aberration function. This is illustrated in Fig. 22 for profile arising from a central ray entering the LPSD at 0.5° offcentre when set at a low 2θ (ie., β = 0.5°, 2θ = 20°). When ε 0 is outside the range of ε values dictated by the limiting values of ζ in the incident beam (typically ±0.5°) then there are no discontinuities within J PSD (ε), no infinities and J PSD (ε) is finite at all ε as illustrated in Fig. 23 for the profile with β = -2°. It should be noted that the aberration functions for a particular 2θ are not symmetric with respect to β. In Fig. 22 the profiles for β = +2° and β = -2° at 2θ = 20° are distinctly different and have different limits in ε.

Parallax Error
When a diffracted x-ray photon enters a LPSD its path in the detector gas is not perpendicular to the anode wire, except at the centre of the detector, and additional broadening, known as parallax broadening, is introduced into off-centre diffraction lines. This arises because the ionisation caused by incoming photons is likely to occur at any point along its path, but the subsequent avalanche is perpendicular to the anode wire (see Fig. 24). For a detector gas with low absorption, the profile shape recorded by the LPSD at an angle ±β from the centre of the LPSD is a unit area impulse function with a width ∆ PX (°2θ) given by, (23) where D is the depth of the detector and R is the radius of the diffractometer. Parallax error can be reduced considerably by increasing the detector gas pressure. In this way the x rays are absorbed by a thin layer of gas beneath the detector window, in which case the diffractometer radius R is defined by the front window of the detector rather than its anode wire.
Most detectors are designed with a quantum efficiency of at least 80 % and have depths somewhere between 5 mm and 10 mm. The aberration profiles for these conditions resemble those shown in Fig. 23 where the absorption of the detector µ gas > 0. The shape of the aberration profiles for parallax broadening J µ PSD (ε) at an angle β to the centre of the LPSD is given by, where µ g is the linear attenuation coefficient of the detector gas [27]. Parallax broadening can be quite large for long detectors and is usually the dominant aberration when LPSDs are used for peaks above 2θ ≈ 40°. For example, when β = ±5° and the detector has a depth D ≈ 8 mm, the breadth of the parallax function ∆ p can be as large as 0.2°.

Thermal Noise
The ultimate resolution of a LPSD is the spatial uncertainty of the position measurement of an individual x-ray photon incident normal to the anode wire of the detector. This is controlled by the spatial broadening of the discharge produced by incoming photons, by the thermal noise generated in the anode wire and by the accuracy of the position-encoding electronics [28]. In most instances the thermal noise profile J NPSD (ε) is modelled as a Gaussian function with a FWHM σ, (25) This function represents the distribution of "measured positions" of a discharge caused by an incoming x-ray photon when the photon is incident normally at exactly the same position along the wire. In most LPSDs, the angular width σ is equivalent to a positional uncertainty along the anode wire between ∆x ≈ 0.04 mm and 0.20 mm where σ = ∆x/R radians. This effect contributes to the broadening much like that from the typical receiving slit width in a conventional diffractometer. Smaller ∆x values down to 0.040 mm are possible when the gas in the LPSD is under pressure. Most LPSDs have a positional resolution ∆x around 0.1 mm although those that utilise delay lines can be as large 0.2 mm.

The Wavelength Distribution in Laboratory Diffractometers
The natural shape of the energy distribution W(E) of a single characteristic x-ray emission line well above the threshold energy is Lorentzian and given by (26) where E 0 is the peak emission energy and Γ is the lifetime broadening given by the sum of the widths of the two relevant atomic levels involved in the transition.
When a Lorentzian energy spectrum is transformed into 2θ space it remains Lorentzian, provided dE/d2θ is relatively constant across the profile, with a FWHM Γ 2θ in 2θ space, given by (27) where Γ and E 0 are in eV. The form of the Kα 1 α 2 emission profile from Cu is shown in Fig.24. For all of the transition element anodes used in x-ray diffraction neither the Kα 1 line nor the Kα 2 line has a Lorentzian shape and both the Kα 1 and Kα 2 lines are asymmetric with extended high angle tails. Moreover, the asymmetries and FWHM values of the Kα 1 and Kα 2 peaks are different. This is particularly evident with the first series of transition element and is related to the anomaly in the atomic number dependence of the atomic level widths of the L II and L III levels obtained from xray photoelectron spectroscopy data and x-ray emission spectroscopy [29]. The line widths Γ, asymmetry indices κ and energies E 0 of the four transition elements commonly used for x-ray targets are given in Table 1.
In FPPF an accurate model of the emission profile is essential particularly for the analysis of high angle profiles. The FWHM Γ 2θ of the Kα 1 line from Cu and other commonly used transition metal anodes is shown in Fig. 25. At 2θ = 40° the FWHM Γ 2θ is less than or equal to 0.01° (2θ) is considerably smaller than the FWHM values of actual diffraction lines from commercial diffractometers which are typically in the range 0.07° to 0.10° 2θ [21] depending on the choice of slits. In this region the contribution of the emission profile is swamped by the geometrical aberrations. When 2θ > 60°, however, the emission profile tends to dominate over the geometrical aberrations. Once 2θ > 100°, the total breadth of the geometrical aberrations is relatively minor and the profile shape conforms closely to the emission profile.
Line widths Γ, asymmetry indices κ and emission energies E 0 for the Kα 1 and Kα 2 emissions of selected transition elements (data from Salem and Lee [30]) The natural asymmetry of the emission lines arises from the multiplet structure of the transitions. In addition to the main transitions involving the change in vacancy state 1s → 2p, it has been recognised that 3d spectator transitions also contribute up to 30 % of the Kα 1 α 2 emissions [31]. In these transitions the atom is doubly ionised and the actual vacancy transition is still 1s → 2p, but the second vacancy in the 3d level is not directly involved in the transition. The notation for this transition is 1s3d → 2p3d. A phenomenological representation for accurately describing the asymmetric Cu emission profile was first used by Berger [32]. In this model the Kα 1 and Kα 2 lines were each represented as two Lorentzian profiles as shown in Fig. 26. This representation was also used successfully on accurate spectroscopic data from Cu obtained by Härtwig et al. [33]. A systematic study of the Kα and Kβ emission profiles from Cr, Mn, Fe, Co, Ni, and Cu [34] has shown that the phenomenological representation can be used to accurately represent these elements down to an R factor of 1 %, although in some cases it is necessary to use up to seven Lorentzians.
Another feature of Kα emission lines which needs to be included in an accurate line profile fitting model is the satellite multiplet structure in the high energy tail as shown in Fig. 27. As a group these have an intensity of ≈0.6 % of the Kα 1 emission line in the case of Cu rising uniformly with decreasing atomic number up to ≈1.4 % for Cr [35]. Satellites lines are also evident in Kβ 1 β 3 spectra and appear on both the low energy and high energy tails as the Kβ' and Kβ" lines [36]. The Kα 3 α 4 non-diagram lines arise from the transition 1s2p → 2p 2 in which the actual vacancy transition is 1s → 2p, but in the presence of a 2p spectator hole [37,38]. For most of the first transition series the Kα satellite structure can be fitted accurately with four or five Lorentzians [35,38]. In most x-ray diffraction studies this level of precision in fitting the satellites is unnecessary. For Cu it is sufficient to represent the Kα satellite group as a single broad Lorentzian so that the total Kα spectrum can be represented by five Lorentzians as given in Table 2.    Table 2. Relative intensities I 0 (areas), wavelengths λ and lifetimes widths Γ λ (in λ units) for representing the Cu Kα spectrum by five Lorentzians. The Kα 1 α 2 data were taken from Höltzer et al. [34] and the satellite data from Cheary and Coelho [2]. In the transformation of the distribution W(λ) from wavelength space to 2θ space at 2θ < 130° the Lorentzian shape in λ maps directly into a Lorentzian in 2θ with minimal error. At higher angles the effects of dispersion between λ and 2θ become increasingly evident [12] and the transformation needs to be carried more accurately. The full expression for the transformed wavelength distribution in 2θ space for a sum of Lorentzians in λ space is (28) where I 0i and Γ λi are the relative areas and wavelength FWHM of the ith Lorentzian. The d spacing is defined in relation to the reference wavelength λ ref of the emission line which in the TOPAS implementation of this procedure is the wavelength of the highest intensity Lorentzian. In the case of the Cu Kα spectrum given in Table 2, the reference wavelength is that of the Kα 1a component (λ = 1.540591 Å) so that λ ref = λ Kα1a = 2d sinθ Kα1a . At very high 2θ values (ie., 2θ ≥ 150°) the cosθ terms arising from d2θ /dλ can change quite significantly over a profile and elevate the intensity in the high angle tail. In addition to dispersion, the distortional effects of the Lorentz factor on high profiles should also be incorporated when analysing profile shapes at very high angles [39,40,41].
As the natural shape of an emission profile is Lorentzian, the tails can extend a considerable distance from the central peak as shown in Fig. 27. Most diffractometers, however, operate with Kβ filtering or some form of monochromatisation and as such the tails are attenuated to varying degrees. When a Ni Kβ is included in the beamline of Cu Kα instrument, the attenuation appears to be more or less uniform across the profile except below the Ni K absorption edge. Careful analysis however, reveals a small variation in the attenuation across the profile owing to the increase in linear attenuation coefficient with increasing wavelength (ie., µ α λ 3 ), but this is generally a small effect and is not expected to affect the profile shape significantly [12]. The inclusion of a curved graphite monochromator in the diffracted beam of a diffractometer greatly reduces the range of wavelengths entering the detector resulting in profile tails that diminish more rapidly than the natural emission profile as shown in Fig. 28. These monochromators can also affect the relative intensity I Kα2 /I Kα1 ratio by up to ±10 % depending on the alignment and setting of the crystal. For example, in the monochromated spectrum shown in Fig. 28 the relative intensity of the CuKα 1 /Kα 2 is approximately 0.46 rather than 0.50 as in the unfiltered spectrum.
Commercial pyrolytic graphite monochromators are strongly oriented polycrystals bent by a hot pressing operation. The focussing is imperfect and the resolution is relatively poor. Better resolution and, in some cases, high intensities can be obtained from "ground and bent" (Johannson) single crystals. These crystals achieve perfect focussing when correctly aligned and are able to select a very narrow wavelength band. The most common materials used for Johannson monochromators on conventional diffractometers are quartz, germanium and silicon. The example shown in Fig. 29 shows the wavelength spectrum from an asymmetrically cut ground and bent Ge crystal used as an incident beam monochromator. The wavelength passband is narrow enough to remove 99.98 % of the Kα 2 component and 100 % Kα satellites from the CuKα spectrum and, almost completely eradicate the Lorentzian tails of the emission profile.
In the presence of monochromators, even low resolution graphite monochromators, it is no longer possible to accurately model the wavelength distribution from an x-ray tube using the tabulated unfiltered spectra such as the one in Table 2. Although first principles calculations of the wavelength transmission function through ideal monochromator are possible, it is currently more practical to determine experimentally a "learned" spectrum for the curved graphite, germanium, quartz and lithium fluoride monochromators used in the majority of monochromated laboratory powder diffractometers. This can be done by modifying the "sum of Lorentzians" representation in energy space or λ space to fit the spectrum entering the detector. In broad terms, monochromators reduce the width of the wavelength distribution and tend to truncate the tails of spectra. A number of approaches can be used to accommodate these changes, • represent the components of the wavelength distribution as Voigt or pseudo-Voigt functions rather than Lorentzians, to limit the extension of the profile tails, and modify the relative intensities of the component. • represent the effect of the monochromator as a wavelength filter with a transmission function T(λ) represented by a simple function, such as a split pseudo-Voigt or split Pearson VII function with up to four refineable parameters, which operates on the tabulated Lorentzian emission profile data. The split functions are used to incorporate asymmetry in the T(λ). In both cases the parameters of each representation are obtained by analysing and fitting the high angle profiles from either a reference line profile standard, such as LaB 6 SRM 660a, or single crystal disc such as a 111 wafer of silicon.
The inclusion of a parabolic multilayer mirror in the incident beam of a diffractometer can also introduce a distortion into the wavelength spectrum [42]. This happens because the Kα 1 and Kα 2 components of the spectrum reflect off the mirror in slightly different directions as shown in Fig. 31. The separation of the Kα 1 and Kα 2 peak maxima, ∆2θ Kα21 = 2θ (Kα 2 ) -2θ (Kα 1 ), in a line profile is then either larger or smaller than the same profile from a diffractometer with no mirror. When a diffractometer is set up with its mirror in the orientation shown in Fig. 30, the separation ∆2θ Kα21 is smaller than the value expected from the known Kα 1 and Kα 2 wavelengths by an amount corresponding to the difference ∆ψ in the directions of the two incident beams on to the specimen. According to Toraya and Hibina [42] the difference δψ for a conventional long fine focus x-ray with a projected target width of 0.04 mm can rise to 0.0017° for a high resolution mirror, but decreases to a negligible level for low resolution mir-Volume 109, Number 1, January-February 2004 Journal of Research of the National Institute of Standards and Technology 21 Fig. 29. Wavelength spectrum emerging from an asymmetrically cut Ge111 ground and bent incident beam monochromator presented linearly and logarithmically. The Kα satellites are completely removed, but the Kα 2 is still present at ≈0.02 % of Kα 1 even in a well aligned system. rors. In any high accuracy lattice parameter determination using this type of diffractometer it is necessary to incorporate the wavelength dependent zero error into the analysis.

Fundamental Parameters Profile Fitting (FPPF) in Practice
In practice FPPF requires accurate numerical procedures for carrying out multiple convolution integrals. This can be done by representing the calculated profiles as a histogram and reducing the convolution integral to a summation. To avoid systematic errors with this approach the angular step size between calculated intensities needs to be very small but in so doing the operation becomes very time consuming [2]. In the XFIT and TOPAS implementations of FPPF a semianalytical procedure has been developed for convoluting the aberration functions. In this procedure the two aberration functions being folded together are calculated at the same 2θ values as the measured data points and then a continuous function is formed by interpolating between the calculated points. As the calculated functions are then each a series of linear sections it is possible to calculate the convolution integral analytically [2]. Some difficulties are experienced with aberration functions that possess singularities, but in all cases the functions are integrable and the effect of the singularity can be overcome either by a convolution process or by a smoothing operation [4].
Although FPPF is a powerful method of profile analysis it can fail when used without recognising its limitations. The physical parameters describing a diffractometer are not statistically independent within the least squares refinement procedure and strong correlations exist between many of the refineable parameters. For example, when the angular acceptance angles, ∆ I and ∆ D , of the incident beam and diffracted beam Soller slits are refined independently to a set of profiles, various combinations of refined values ∆ I and ∆ D can be obtained that give equally good fits to the data. The reason for this is that these two parameters are very strongly correlated and the values can in fact be interchanged without changing either the shape of a synthesised profile significantly or the quality of the fit. As a consequence the best fit with an unconstrained refinement often corresponds to refined parameters that do not make physical sense. When first using the FPPF approach it is advisable to investigate the validity of the refined parameters for a particular diffractometer using a well crystallised reference specimen and establish the parameters of the instrument that allow profiles over the whole 2θ range to be defined with one set of values. It would be unrealistic to expect the refined parameters for a diffractometer to match the directly measured val-Volume 109, Number 1, January-February 2004   Table 3. ues exactly as there are a number of second order effects in diffractometer profiles that are not incorporated in the fitting model. Moreover, not all of the instrumental aberrations are independent and as such the convolution model is not strictly valid with certain combinations of aberrations [8]. Nevertheless, experience with most of the commercially available diffractometers has shown that refined values reasonably close to the true instrumental values can be expected. In the original investigation of fundamental parameters fitting [2] it was shown that when deliberate changes were made to the diffractometer set-up, such as changing the receiving slit length, receiving slit width and diffractometer radii, the change in the refined instrument values corresponded well with the actual changes despite the limitations in the axial divergence model used at that time. In general, therefore, if physically unrealistic instrumental parameters are required to describe the diffraction pattern of a reference material then there is either a deficiency in the model used to describe the diffractometer, some sort of mis-setting of the diffractometer or the refinement is trapped in a false minimum.
Once the refined instrument parameters have been established (ie., the dimensions and apertures of the various slits, the source size and the wavelength distribution), and these are in reasonable agreement with the actual values, the only instrument parameter that may need refinement from specimen to specimen is the linear attenuation coefficient µ. Absorption effects in a line profile are only evident when µ < 100 cm -1 and under these circumstances there is some justification in making µ a refineable parameter. In specimens with µ ≥ 200 cm -1 refinement of µ normally has little effect on the profile shape because the specimen transparency profile is relatively narrow at all 2θ angles (ie., ≤0.01°2θ). Consequently, µ is fixed at a representative value and not refined.
Even in well crystallised powders with crystallite sizes up to 2 µm crystallite size broadening is detectable in high angle lines. Conversely, in many powders it is not uncommon to have crystallite sizes down to 0.5 µm and crystallite size broadening is evident even in low angle lines. Also, in many crushed powders the action of crushing even in a standard pestle and mortar can induce microstrain broadening. In TOPAS and XFIT therefore, the apparent crystallite size T app and the percent microstrain ε rms are frequently included as refineable parameters although their contribution is generally small. Various profile functions can be adopted for crystallite size and microstrain and commonly take the form of a Lorentzian function or a Gaussian function. The most common unit area functions used for crystallite size and microstrain, B cryst (2θ) and B µ (2θ), respectively, are where the FWHM of the Lorentzian = H cryst = 180 λ/(πT app cosθ o ) °2θ and the FWHM of the Gaussian = H µ = (18ε rms ) tan θ 0°2 θ. In some crushed powders it has been necessary to use a Lorentzian function to represent the microstrain but with a FWHM = ηtanθ where η is a constant related to the microstrain [4].
The quality of the fits obtainable from the FPPF approach and the degree to which the refined instrument parameters agree with the actual values is illustrated below for data collected from a polycrystalline MgO reference specimen over the range 2θ = 36° to 150° using CuKα radiation. The specimen used in this instance was a 20 mm disc prepared by sintering at 1400°C for 24 h in air. When a small segment of the disc was examined in a scanning electron microscope the vast majority of crystallites were between 1 µm and 2 µm in diameter. Table 3 shows the refined values obtained for the MgO data from TOPAS whilst Fig. 30 shows the fits obtained to two of the profiles that are sensitive to the instrument parameters and the specimen absorption. Only four profile shape parameters were refined and these were the absorption coefficient µ, the receiving slit width w r , the angular aperture of the incident beam Soller slits ∆ I and the apparent crystallite size T app . All the other instrument parameters were fixed at their actual values. Some parameters, such as the equatorial divergence α, the width of the x-ray source and the sample width L s , have very little effect on the profile shapes and a wide range of values can be accommodated for each of these parameters without affecting the quality of the fit. Although the receiving slit length L r and x-ray source length L x can have a significant effect on the profile shape, they were not refined because of their strong correlation with the Soller slit aperture ∆ I and the fact that the improvement in fit with their inclusion was minimal. Equally good fits are obtained by fixing L x and ∆ I and allowing L r to refine. The equatorial divergence angle α and other axial divergence parameters are normally included in a refinement when low angle diffraction lines, preferably below 2θ = 25°, exist in the data set as this region of the diffraction pattern that is sensitive to equatorial and axial divergence. The inclusion of crystallite size broadening in most refinements is very important. In the present MgO refinement its removal from the refinement resulted in the other refined parameters taking on physically unrealistic values and a rise in the R wp value from 3.5 % to 6.2 %. It was also clearly evident that the calculated profiles no longer fitted the tails of the observed profiles.

Concluding Remarks
Fundamental parameters profile fitting offers a number of benefits as a method of profile analysis. It is based on a physical model of the diffractometer and its refined parameters should be self consistent with physical dimensions of the diffractometer and the physical properties of the sample. On this basis it can therefore identify whether or not a diffractometer is operating at optimum resolution for the conditions used and provide a means for assessing the performance of a diffractometer in a particular application. For example, a knowledge of the performance of a diffractometer operating under asymmetric conditions can be important in choosing the best 2θ range and slits that will maintain sufficient resolution to monitor line positions in stress analysis. As the profile shape is known, the FPPF technique also provides greater certainty in the identification of weak peaks or impurity lines embedded in the tails of stronger line. In Rietveld analysis or quantitative analysis, FPPF allows the profile shapes across the whole 2θ range to be fitted without any instrument based parameters in the refinement. The focus of the refinement is therefore on the diffraction effects of the specimen and not the instrument. Moreover, line broadening analysis is an integral feature of FPPF for both individual lines or as part of a Rietveld analysis and, correction for instrument broadening is an intrinsic part of the analysis making reference specimens unnecessary in many circumstances. FPPF has its weaknesses. Up until now it has only had limited success as a method of accurate lattice parameter determination and in some diffractometer set-ups the FPPF instrument parameters can differ significantly from the actual values. Although it corrects the 2θ positions of lines for instrument line shift, including zero shift and specimen displacement, it has not been possible to obtain a set of lattice parameters for a particular specimen which are the same, within an uncertainty < ±0.0005 Å, for every hkl line in the pattern [43]. In addition, the FPPF approach does not have a physically based model for incorporating mirrors and monochromators into wavelength distribution and the axial divergence profile although recent work by Masson et al. [44] suggested an analytical approach in which devices within the optical path of a diffractometer can be readily incorporated into the theoretical formalism for describing the profile shape. The current physical modelling of laboratory based diffractometers and some of the diffraction processes within FPPF is adequate for many applications, but some degree of refinement of the theoretical procedures underpinning the technique is still necessary before the technique can claim to accurately describe both the shape and position of powder diffraction lines. Table 3. Comparison of actual diffractometer and sample parameters with values obtained for MgO data by FPPF refinement. In this refinement the radius of the diffractometer was fixed at 200 mm and the CuKα spectrum adopted for fitting was the one given earlier in Table 2