 1. Introduction
 2. Ellipse geometry – the intersection of a cone's axis and its elliptical conic section
 3. Experiment and method
 4. Beam energy determination and calibration of a twodimensional detector
 5. Discussion
 6. Conclusion
 A1. Error due to thermal expansion of the calibration standard
 A2. Error due to detector stage displacement
 A3. Error due to Xray attenuation
 A4. Other sources of error
 References
 1. Introduction
 2. Ellipse geometry – the intersection of a cone's axis and its elliptical conic section
 3. Experiment and method
 4. Beam energy determination and calibration of a twodimensional detector
 5. Discussion
 6. Conclusion
 A1. Error due to thermal expansion of the calibration standard
 A2. Error due to detector stage displacement
 A3. Error due to Xray attenuation
 A4. Other sources of error
 References
research papers
Complete elliptical ring geometry provides energy and instrument calibration for synchrotronbased twodimensional Xray diffraction
^{a}Diamond Light Source Ltd, Diamond House, Harwell Science and Innovation Campus, Didcot, Oxfordshire OX11 0DE, UK
^{*}Correspondence email: michael.hart@diamond.ac.uk
A complete calibration method to characterize a static planar twodimensional detector for use in Xray diffraction at an arbitrary wavelength is described. This method is based upon geometry describing the point of intersection between a cone's axis and its elliptical conic section. This point of intersection is neither the ellipse centre nor one of the ellipse focal points, but some other point which lies in between. The presented solution is closed form, algebraic and noniterative in its application, and gives values for the Xray beam energy, the sampletodetector distance, the location of the beam centre on the detector surface and the detector tilt relative to the incident beam. Previous techniques have tended to require prior knowledge of either the Xray beam energy or the sampletodetector distance, whilst other techniques have been iterative. The new calibration procedure is performed by collecting diffraction data, in the form of diffraction rings from a powder standard, at known displacements of the detector along the beam path.
Keywords: instrument calibration; synchrotron radiation; Xray diffraction.
1. Introduction
The use of twodimensional detectors is widespread in the collection of Xray diffraction data. Their high resolution and suitability for use at high Xray energies allows a variety of experiments to be performed, such as protein crystallography (Helliwell, 1982), singlecrystal diffraction (Riekel, 2000), the study of grain orientation and strain in polycrystals (Hofmann et al., 2012), the evolution of in response to strain (Evans et al., 2012), and the study of thick samples in complex environments (Sui et al., 2010).
We present a noniterative method for the calibration of a planar twodimensional detector used for synchrotronbased Xray diffraction at an arbitrary wavelength. The calibrated parameters include the sampletodetector distance, the position of the Xray beam trajectory on the detector, the detector tilt and the beam energy. Our approach is made possible through the presentation of a single coherent geometrical description and uses diffraction ring data collected from a number of detector positions along the beam path. Our method accounts for the random variation in detector tilt and lateral position as the detector is moved during data collection. This calibration technique is useful in a facility where the energy can be freely chosen but difficult to evaluate accurately via other methods.
There are many diffraction geometries in use for the collection of Xray crystallographic data, and these have previously been described (Helliwell, 2004). Here we are concerned with the geometry of the detector relative to the Xray beam. Where possible, we have adopted a notation used elsewhere (Hinrichsen et al., 2008). A description for the geometry of the intersection of a cone's axis and its elliptical conic section has been attempted elsewhere (Hinrichsen et al., 2008) and will be discussed in more detail later. We do not consider parabolas or hyperbolas, and thus our method applies to cases where both the energy is high enough and the detector's tilt is set appropriately such that calibration data from a powder standard can be recorded as ellipses upon the detector surface.
At low energies, the diffraction pattern of a calibration standard can be recorded and the topology of that pattern used to estimate the diffraction geometry and incident wavelength (Norby, 1997). This method can be sufficiently accurate at large Xray wavelengths ( Å), when large diffraction angles are present. With twodimensional detectors, this indirect calibration method is common practice, owing to its practical implementation in various software packages (Heiney, 2005; Hammersley et al., 1996; Ilavsky, 2012) which are currently used in many facilities around the world.
With a twodimensional detector it is possible and often desirable to conduct the diffraction experiment at short wavelengths – or high photon energies – in order to collect a large section of et al., 2012).
At high Xray energies, the small diffraction angles (°) lead to an interaction between the sampletodetector distance and the wavelength that can be hard to separate mathematically when solving the relevant nonlinear equations (Bragg's law). This interaction makes their simultaneous difficult (HongA diffraction ring resulting from an interplanar spacing given by d, recorded on an orthogonal detector at a distance of D from the sample, using highenergy Xrays of wavelength λ, will have a radius approximately given by . It is not possible to solve this equation for both λ and D using diffraction data taken from a single sampletodetector distance. At best the product of λ and D can be found. Thus, at high energies, it is often the case that either the sampletodetector distance or the beam wavelength needs to be accurately determined prior to the of the other calibration parameters. One possible method to indirectly determine the incident wavelength is to configure the monochromator to produce a wavelength that matches the of a known sample with a well defined K edge; however, this will only cover a discrete set of energies and not always the desired one. An alternative method for calibration at high energies, which uses the extra information obtained by varying the position of the detector, exists (Hong et al., 2012) but is iterative in its application and thus requires the design of an exit strategy.
Another difficulty with calibrating a twodimensional detector using diffraction ring data is that determination of the detector tilt relative to the beam direction is nontrivial. The eccentricity of a diffraction ring,^{1} seen as an ellipse, can be expressed as a function of the detector tilt and the [see equation (11) derived later]; as such it is impossible to calculate the detector tilt using a single diffraction ring without knowledge of the (which in itself is a function of the interplanar spacing and the beam wavelength). Fig. 1(a) demonstrates this ambiguity, where it can be seen that cones of various opening angle and orientation can intersect a plane to produce an ellipse of identical eccentricity.
For each diffraction cone, as seen in Fig. 1(a), that projects to give an identically shaped diffraction ring, the location of the beam centre within each ellipse will differ as per equation (5). The beam centre is not located at either the ellipse centre or the ellipse focal point of a diffraction ring [see equation (10) derived later]. As a consequence the beam centre cannot be defined from one diffraction ring. For a set of concentric cones the positions of the ellipse centres and the ellipse focal points differ for rings originating from different (hkl) planes. This viewpoint is in contrast to that presented elsewhere (Norby, 1997; Rajiv et al., 2007; Hinrichsen et al., 2008), where it was incorrectly stated that the cone's axis intersects the elliptical conic section at its focal point. Subsequent geometrical expressions involving the focal point were thus incorrectly derived (Rajiv et al., 2007; Hinrichsen et al., 2008).
By using more than one diffraction ring it is possible to derive values for the detector tilt and the Xray energy using equation (11), if the interplanar spacings for each diffraction ring are known. However, for small Bragg angles or detector tilt . Equation (11) becomes insensitive to the precise value of the and so cannot be used to attain a good value for the beam energy. To overcome this problem, we use diffraction ring data collected from multiple detector distances to first calculate the beam energy and then give a value for the detector tilt.
By using data from multiple detector positions of known separation we can solve for both λ and D. Our method considers geometries where the detector normal may exhibit an arbitrary angle relative to the incoming beam, regardless of whether the tilt is applied intentionally or not. It is important to allow for a variation in tilt, due to mechanical runout, as the detector stage is translated whilst collecting diffraction data for calibration. For this reason simple ray tracing techniques may be inaccurate.
The closed form solution we present here yields results for all calibration parameters in a noniterative manner and at the same time adds understanding to the underlying geometrical problem. Iterative techniques (Hong et al., 2012) lacking a complete model may yield sufficiently accurate results, if the exit strategy is suitably designed.
The quality of our method derives from a complete analytical threedimensional model which utilizes powder diffraction patterns from a standard, recorded at different sampletodetector distances (Fig. 2). This allows us to remove the ambiguity in sampletodetector distance and wavelength (Bragg's law) as well as the ambiguity in detector tilt and opening angle of a diffraction cone (Fig. 1a).
Our method can be broken down into a number of distinct steps (Fig. 3). Each diffraction ring recorded on a twodimensional detector is an ellipse. The location of the ellipse centre and the length of the semimajor axis in the coordinate system of the detector need to be accurately determined. Diffraction data are collected from a number of different nominal sampletodetector distances. At each detector position, the ratio of the sampletodetector distance to the sine of the tilt and the location of the beam centre on the detector are found by solving a nonlinear equation. Data from all detector positions are then brought together in another nonlinear equation which provides values for the sampletodetector distance and the beam wavelength. Finally, the detector tilt and direction of tilt can be calculated from the sampletodetector distance and the beam centre, respectively. By using data from all available diffraction rings we can improve the accuracy of these calculated parameters. Solving an overdetermined system reduces the sensitivity of the calculated parameters to noise.
Having fully calibrated the detector in terms of its five spatial parameters and the incident beam wavelength, Q, can be assigned to each detector pixel. This can be achieved by using equations (19) and (20).
values,The method presented in this article requires the parameterization of the elliptical diffraction ring formed at the intersection of the diffraction cone and the plane of the twodimensional detector. This procedure is beyond the scope of this article; however, a number of techniques are known (Hanan et al., 2004; Hart & Drakopoulos, 2013). The technique used in this article directly fits an intensityweighted ellipse to peak centres determined at equal intervals in azimuthal angle. This technique works well for both complete and spotty rings (Hart & Drakopoulos, 2013).
To aid our attempt in providing an accurate method for detector calibration, we first present a complete geometrical description of the intersection between a cone's axis and its elliptical conic section.
2. Ellipse geometry – the intersection of a cone's axis and its elliptical conic section
If a twodimensional detector is positioned perfectly orthogonal to the Xray beam, then a diffraction ring will appear as a circle. A deviation from orthogonality will result in an elliptical ring. The righthand diagram of Fig. 2 demonstrates the convention used in this article in the case where a detector exhibits a tilt. In the plane of an orthogonal detector we consider a single tilt, t, of the detector, whilst the azimuthal direction of the tilt relative to the x axis, β, is perpendicular to the axis of tilt. We use a subset of the Euler angles.
Fig. 1(b) shows a view (Kumar, 2006) along the tilt axis of the detector, depicting the Xray beam (propagating along line S–H), a diffraction cone of halfangle , where θ is the and a representation of the detector surface (bold offhorizontal line) intercepting the cone. The intersection of a cone and the twodimensional detector (placed at an angle t) is an ellipse. This ellipse represents the recorded diffraction ring on the detector. The ellipse is shown projected from the detector surface, and both the semimajor, a, and semiminor, b, axes are indicated. The distance from the centre of the ellipse to the point where the Xray beam strikes the detector surface, c, is also indicated. D is the distance from the sample, S, to the detector at H. Labelled on both the plan view and the ellipse, u is the radius of a ring as it would be recorded on an orthogonal, t = 0, detector at position H.
Other authors have stated that the direct Xray beam or axis of a cone intersects a detector or plane at the focal point of the diffraction ring or ellipse (Norby, 1997; Rajiv et al., 2007; Hinrichsen et al., 2008) or at a symmetry centre defined by the diffraction ellipses (Cervellino et al., 2006, 2008). We postulate that, should a cone intersect a plane to produce an elliptical conic section, then the cone axis intersects the plane at a point that lies between the ellipse centre and closest focal point of the ellipse. This can be visualized by considering the threedimensional geometrical construct known as Dandelin spheres (De Villiers, 2007; Dandelin, 1822).
A number of expressions describing the intersection of a plane and a cone, in the context of twodimensional Xray diffraction, can now be derived.
The similarity of the triangle spanned by u and D, the triangle spanned by Gcos(t) and D+Gsin(t) (Kumar, 2006), and the triangle spanned by gcos(t) and Dgsin(t) gives the following geometrical relationship:
Using the substitutions g = ac and G = a+c, equation (1) can be rewritten to give the semimajor axis, a, in terms of the sampletodetector distance, the detector tilt and the distance from the ellipse centre to the beam centre:
Expressions for a and c can also be found in terms of the sampletodetector distance, the detector tilt and the Equation (1) gives
and
and thus
Equations (3) and (4) are equivalent to expressions derived elsewhere (Hinrichsen et al., 2008).
In polar coordinates, the distance from the ellipse centre to any point on the ellipse as a function of the angular coordinate, γ, measured from the major axis is
Use of the substitutions and gives an expression for the semiminor axis, b, in terms of the sampletodetector distance, the detector tilt and the Bragg angle:
and thus
The distance from the centre of the ellipse to the focus, f, can now be given:
and
The eccentricity of the ellipse, e, can also be described:
Our expressions for the semiminor axis and the eccentricity [equations (7) and (11), respectively] disagree with those derived elsewhere (Hinrichsen et al., 2008). We do, however, agree with Hinrichsen in that the eccentricity is a function of the tilt, t, and the θ, which is in contrast to an alternative view that a single value for eccentricity can be used to describe the detector tilt (Cervellino et al., 2006, 2008). Furthermore, it is possible to conclude from equation (10) that, if the detector tilt is such that the detector intersects the diffraction cone to form an ellipse, then the point where the cone's axis intercepts the ellipse is closer than either focal point is to the ellipse centre. In other words, if , then .
Directly derived from the Cartesian equation of an ellipse, the following relationship linking the semimajor axis, a, the semiminor axis, b, the distance from the ellipse centre to the beam centre, c, and the radius of the ring on an orthogonal detector, u, holds:
The expressions derived are summarized in Table 1.

This geometrical description provides a framework upon which our calibration procedure can be based.
3. Experiment and method
The calibration of the diffraction geometry for a twodimensional flatpanel detector was carried out at the Joint Engineering, Environmental and Processing beamline (I12) at Diamond Light Source Limited (UK). The detector used was a flatpanel Pixium RF4343 (Thales), with a pixel size of 148 × 148 µm. The Pixium family of detectors is widely used for highenergy Xray diffraction (Daniels & Drakopoulos, 2009; Sui et al., 2010; Evans et al., 2012; Hofmann et al., 2012). Note that this article does not address any spatial distortion which may result from the curvature of a generic twodimensional detector. Images used in our calibration procedure should already be corrected for such distortions.
A monochromatic beam with a nominal energy of 88.012 keV was generated using a doublecrystal Laue Xray monochromator. The energy was set by scanning the monochromator through the Kshell of 125 µmthick lead foil. The beam size was defined to be 200 × 200 µm, using a set of tungsten slits positioned before the sample stage.
Calibration data were collected, in transmission mode, from a fine powder cerium dioxide (CeO_{2}) standard (NIST Standard Reference Material 674b). The powder standard was mounted in a 0.5 mmthick planar arrangement so as to minimize any azimuthal geometrical sample effects typically produced by a sample in a capillary arrangement (Norby, 1997). The detector was fixed to a highprecision translation stage (Advanced Design Consulting USA Inc.), allowing the detector to be positioned at different sampletodetector distances. All data were collected on the Pixium detector in highresolution mode (2880 × 2881 pixels) with an exposure of 4.0 s.
Data were collected for two commonly occurring situations: the detector was positioned to be orthogonal (as determined by the back reflection from a guidance laser), and the detector was positioned with a nominal horizontal tilt of 6.4°. In both cases data were collected at eight nominal sampletodetector distances, from 580 to 1980 mm in steps of 200 mm.
All data analysis was performed using customwritten software.
4. Beam energy determination and calibration of a twodimensional detector
Calibration of the spatial position of the twodimensional detector and the determination of the beam energy requires the collection of twodimensional diffraction data from a known standard at a variety of different sampleto detector distances.
Determination of each beam centre is achieved through a purely geometrical calculation. For this part of the calibration procedure the sample only needs to produce diffraction rings that can be accurately described. However, the part of our procedure that determines the sampletodetector distance and the beam energy uses knowledge of the sample's lattice spacing. The measurement error in the CeO_{2} standard has been selected for use.
as measured from a fitted ellipse is greater than that of the lattice spacing from a known standard. Use of the sample's lattice spacing allows for a more accurate calibration procedure, whilst at the same time introducing the beam energy as an important parameter. This is why a fine powderDuring data collection, the sample maintains a constant position whilst the detector is moved to different fixed positions along the beam path. It is important to know accurately the relative change in detector position along the beam path. At each detector position, diffraction data are recorded for calibration. The mechanical runout displayed by the translation stage is taken into account by our calibration procedure.
For each diffraction ring at each detector position, the ellipse centre, (o_{x},o_{y}), and the length of the semimajor axis, a, in the coordinate system of the detector are required. Also note that only equations (2) and (3) are used in our calibration procedure.
Table 2 shows a summary of the parameters that are to be calculated in order to calibrate a twodimensional detector for use in Xray diffraction. The required observables are also summarized. To aid understanding we segregate fitted parameters, with the superscripts fitA, fitB and fitC indicating which part of the procedure they are calculated from (Fig. 3).

4.1. Beam centre calculation
To calculate the beam centre, the parameters describing the diffraction ellipses from different (hkl) are required.
The location of the beam centre, at a given detector position, can be calculated from the data provided in a single diffraction image containing multiple rings. We calculate the location of the beam centre for multiple detector positions.
As a result of translation stage runout, the beam centre on a diffraction image is likely to be nonstationary as the detector is translated along the beam's path. For a given image, the recorded diffraction rings result from diffraction cones with a range of opening angles . These diffraction cones share a common axis, and thus a common beam centre exists within the diffraction image for all diffraction rings. The ellipses that describe the diffraction rings do not, however, share a single common ellipse centre; instead the ellipse centres are positioned along a line coincident with the major axes of the ellipses. The beam centre is also located along this line, and the distance from the beam centre to each ellipse centre can be given by equation (4).
For each detector position, i, employing the method of linear least squares, the ellipse centres are used to give a line of best fit with intercept, , and gradient, . Each line of best fit can be considered to be coincident with the major axis of the ellipse, and thus to relate the x and y coordinates of the beam centre, (h_{xi},h_{yi}):
For diffraction ellipse α, where α denotes different { hkl}, at detector position i we can write a general expression for the semimajor axis, . From equation (2),
where
and can be written as a function of the coordinates of the ellipse centre, , and the x coordinate of the beam centre, h_{xi}, of the diffraction image:
D_{i} is the sampletodetector distance of a detector at position i, and is the difference between the sampletodetector distance at a chosen detector position, D_{0}, and D_{i}. For the purpose of this calculation we have arbitrarily chosen D_{0} to be the distance from the sample to the detector at the detector's furthest position from the sample as used during data collection. Calibration parameters can be found for all detector positions.
Equations (14) and (16) are used for nonlinear leastsquares optimization of the x coordinate for the beam centre, h_{xi}^{fitA}, and the ratio k_{i}^{fitA} at each detector position i. Only data from two diffraction rings are required to solve for h_{xi}^{fitA} and k_{i}^{fitA}; however, the use of all available rings will minimize fitting noise. The y coordinates of the beam centres, h_{yi}^{fitA}, for all i are subsequently found using equation (13). The optimization is performed using previously determined ellipse parameters, and (Hart & Drakopoulos, 2013), from a number of diffraction rings that are recorded on the twodimensional detector. For this nonlinear leastsquares optimization, a good initial value for the ratio k_{i}^{fitA} is the length of the major axis squared from the largest ellipse used, and the initial value for the x coordinate of the beam centre can be taken as the x coordinate of the ellipse centre from the smallest ellipse. The fitted value for the x coordinate of the beam centre is also constrained to lie to one side of all ellipse centres (Fig. 1b).
For two cases – `orthogonal' (the detector is positioned orthogonal to the Xray beam) and `tilted' (the detector is positioned with a tilt measured to be approximately 6.4° about the vertical axis) – Fig. 4 shows the coordinates of the ellipse centres (dots) for eight different sampletodetector distances, D_{i}. Successive images are separated by detector movements of 200.0 mm along the beam path. At each detector position, data from the (hkl) = (111), (002), (022), (311), (004), (422) and (333)/(511) diffraction rings were used where available, and best fit ellipses were found (Hart & Drakopoulos, 2013). For each set of ellipse centres, a line of best fit is shown together with the calculated beam centre (circles). For ideal detector displacements, the beam centres show a linear relationship in h_{x} and h_{y}. Deviations from this are caused by mechanical runout of the translation stage (Fig. 4). For the `tilted' case, we can see that the ellipse centres for each image lie along straight lines; these lines are conceptually coincident with the major axes of the ellipses. In the `orthogonal' case we can see that the ellipse centres also exhibit a linear relationship, which is contrary to what may be expected for an orthogonal detector (a perfectly orthogonal detector will exhibit concentric circular diffraction rings). This is due to the difficulty in aligning the detector to be perfectly orthogonal, and thus a small but quantifiable tilt remains. The small tilt coupled with the variation in the angular position of the detector results in a larger variation in the azimuthal direction of the detector tilt for the `orthogonal' case as compared to the `tilted' case.
4.2. Beam energy, sampletodetector distance and detector tilt calculations
This part of the calibration procedure makes use of data collected from multiple detector positions. This information allows for the simultaneous determination of the beam energy and the sampletodetector distance. Knowledge of the detector's movement along the beam path is required.
In cases where the path of the translation stage is parallel to the beam path, the differences in distance of the detector position along the beam path, , can be taken as the reported movement of the translation stage itself. If the translation stage axis is not parallel to the beam path, then the best fit values for the beam centre, (h_{xi}^{fitA},h_{yi}^{fitA}), can be used to calculate . In the nonparallel case, if denotes the changes in reported position of the translation stage along its own axis, then , where m_{det} is the gradient of the best fit line relating the values [(h_{xi}^{fitA}h_{x0}^{fitA})^{2}+(h_{yi}^{fitA}h_{y0}^{fitA})^{2}]^{1 / 2} to .
As a function of detector position the azimuthal direction of detector tilt may be nonstationary. This does not, however, affect the accuracy of our method, as our calculation makes use of the length of the semimajor axis, a, and the ellipse major axes should all lie within the same plane for a given detector position.
Rewriting equation (3) as a general expression using equation (15) gives
where
from Bragg's law. By substituting the Bragg angles for the interplanar spacing of the crystal lattice, we not only introduce the beam wavelength as an important parameter but also improve the accuracy of our calibration procedure, as the interplanar spacings for a diffraction standard are known to a high degree of accuracy.
Equations (17) and (18) are used for nonlinear leastsquares optimization of the beam wavelength, , and the sampletodetector distance, D_{0}^{fitB}. Only data from two diffraction rings at different detector positions are required to reasonably solve for and D_{0}^{fitB}; however, use of all available data across all rings from all detector positions will improve the accuracy of the fitted parameters. Once D_{0}^{fitB} has been found, the sampletodetector distances for all detector positions can easily be calculated from .
The observables used when solving equation (17) are the previously determined values for the semimajor axes, , and the ratios k_{i}^{fitA}, obtained from solving equation (14). As already stated, the d spacings, , used in equation (18) are known to a high level of accuracy as a diffraction standard is used for detector calibration. The wavelength and the sampletodetector distance are often approximately known, and these can be used as initial values to improve the speed of this nonlinear optimization. In this case the wavelength is approximately set to be at the K edge of Pb. The sampletodetector distance at the detector's furthest position for the `orthogonal' case has been measured, using a tape measure, to be approximately 1977 mm.
For each detector position i, once the sampletodetector distance has been calculated, the detector tilt, t_{i}^{fitB}, can be determined using equation (15), and the azimuthal direction of the tilt, , can be determined from the gradient, , of the line of best fit through the fitted ellipse centres (Hart & Drakopoulos, 2013).
The complete characterization of the detector geometry can be used to calculate theoretical diffraction rings for our CeO_{2} standard. These theoretical rings can then be compared with the observed rings used during calibration (Table 3). The difference between the values for the observed, , and theoretical semimajor axis lengths, , are shown in Fig. 5. The small values of these residuals together with a lack of any apparent trend give an indication of quality for both the model used in this article and the performance of the nonlinear leastsquares optimization. We cannot show residuals for the beam centre as the direct Xray beam is blocked from reaching the detector during data collection.

In practice, accurate calibration values can be attained from one set of diffraction data. Table 4 demonstrates the repeatability of our method and shows average fitted values for the beamline energy and the calibration parameters of the detector at its furthest position from the sample. Mean values are calculated from 29 separate sets of data, whilst the calculated errors are taken as the standard deviation. Also shown is the mean coefficient of determination, R^{2}, for the ellipse semimajor axes, associated with the final fit that determines the beam energy and the sampletodetector distance. Note that the error reported for the measured beam energy reflects our uncertainty in configuring the monochromator to output an energy corresponding to the K edge of lead. This error is based upon our knowledge of the monochromator and a Monte Carlo method (Lenz & Ayres, 1992) which is used to estimate the error in fitting a Gaussian to the first derivative of the absorption edge.

Given accurate wavelength and detector parameters for a fixed experimental detector position (i.e. the subscript i can now be dropped), or Q values can be assigned to each pixel of the detector by using equations (19) (Hammersley et al., 1996) and (20).
where and are the x and y coordinates of the pixels relative to the beam centre and are given by and .
Now each detector pixel has associated with it not just an intensity value but also a or a Q value. Each pixel can then be assigned to an appropriate bin of a given interval in or Q, and a onedimensional diffraction pattern can be directly formed (Fig. 6).
5. Discussion
Twodimensional Xray diffraction is an established scientific technique. Owing to the relationship between beam energy and sampletodetector distance, other calibration techniques have required the accurate determination of one of these values during the course of an experiment prior to data reduction or a process of iteration. For twodimensional detectors, accurate independent measurement of the sampletodetector distance is difficult because the sensor plane is normally within the detector behind a protective layer. Instead the beam's energy can be set to match the absorption K edge of an element such as those provided by the National Institute of Standards and Technology. This procedure can, however, be both time consuming and of limited accuracy. The measured width of the K edge is often appreciably wider than the theoretical width of the absorption K edge, as a result of several instrumental broadening effects. In addition, such an approach restricts an experiment to be performed at an energy matching the K edge of an available element.
With the aid of a highprecision translation stage used to alter the position of the detector along the beam axis, this article has presented an accurate method to simultaneously determine the beam energy and spatial position of the detector relative to both the beam and the sample. The entire calibration procedure, from data collection to data processing, can be accomplished within 10 min, limited only by motor speed, exposure time and computation power. The method does not require prior knowledge of any calibration parameter, other than an estimate of the sampletodetector distance. Furthermore, each best fit of the calibration parameters is set up such that at most only two parameters are determined at a time; this increases both the accuracy and the speed at which a final solution can be found.
To help implement our method of calibration, the geometry of a cone and an intersecting plane was revisited. We have found that the point of intersection between the cone axis and its elliptical conic section has been neglected in the literature and thus have not been able to find any nomenclature for this point. It is the only common point for a set of concentric elliptical conic sections generated from a single intersecting plane. In addition, it appears that the diffraction community considers this point of intersection to be one of the focal points of the ellipse formed by the intersecting cone and plane (Norby, 1997; Hinrichsen et al., 2008), in other words, that the beam centre as recorded on a detector coincides with the focal point of the surrounding diffraction ellipses. In our opinion this assumption is incorrect. By visualizing the construct of Dandelin spheres (De Villiers, 2007; Dandelin, 1822) it is seen that the cone axis does not intersect the plane at a focal point of the ellipse. Our subsequent derivation for the focal point, [equation (9)], and the distance from an ellipse centre to this point of intersection, [equation (4)], supports this view.
As all diffraction ellipses recorded on a twodimensional detector result from a single incident Xray beam striking a sample, all diffraction ellipses share a common beam centre. The distance from the ellipse centre to the point on an ellipse that is coincident with the beam centre, , is a function of θ. This implies that the ellipse centres for the diffraction ellipses do not occupy the same point in space. The same can be said for the focal points, because in general . It is also the case that diffraction ellipses recorded on a single image do not have the same eccentricity: they are dissimilar in shape. The eccentricity of an ellipse is a function of θ [equation (11)].
Our new understanding of the geometry allowed us to create a method for determining the beam centre in relation to the diffraction rings recorded on a detector [equation (14)]. Our method utilizes information from multiple diffraction ellipses but does not require accurate prior knowledge of the sampletodetector distance. To guarantee a highquality fit it is important to ensure that the diffraction data are of high resolution and that a decent method for the determination of ellipse parameters is used. Azimuthal peak broadening effects should also be minimized by using a radially symmetric calibration sample.
Having found the beam centre for each diffraction image, which is later used in assigning a Q value to each detector pixel, it is actually the ratio of the sampletodetector distance to the sine of the detector tilt [equation (15)] that is used in a subsequent calculation to determine the sampletodetector distance itself and the beam wavelength [equation (17)].
The discrepancy between the calculated beam energy values and the theoretical value (Table 4) is considered to be within the limits of uncertainty to which the monochromator can be set to the K edge of Pb. We discuss errors and potential systematic corrections in Appendix A.
The sampletodetector distances for the two cases differ, but this is expected as the vertical rotation axis of the stage on which our detector lies does not pass through the plane of the detecting surface.
In the `orthogonal' case it was found that a small detector tilt exists. This is expected and results from the difficulties in ensuring true orthogonality of the detector to the incoming Xray beam during experimental setup. This small tilt can explain the relatively large error in the reported value for the azimuthal direction of the tilt. The diffraction rings are very close to circular. In the `tilted' case the tilt is relatively large, allowing the azimuthal tilt direction to be more easily determined. In both cases the coefficient of determination, R^{2}, for the ellipse semimajor axes, as determined when solving equation (17), is very close to unity. Such high values of R^{2} provide confidence in the quality of both the model used and the data provided.
With the help of Nika (Ilavsky, 2012), a software package that can calibrate and reduce twodimensional diffraction data, we can check whether our calculated calibration parameters are selfconsistent. For accurate calibration, Nika requires knowledge of either the beam energy or the sampletodetector distance, but by using our calculated beam energy, consistency between the other calibrated parameters and the beam energy can be confirmed. Using an image from the `tilted' case with a calculated energy of 88.0452 keV as an input (note that this is different from the reported average energy shown in Table 4, as this comparison is performed using a single set of calibration data) provides a very good agreement between our calibration parameters and those calculated using the Nika software package. The differences in value are sampletodetector distance 0.055 mm, beam centre 0.009 mm and detector tilt 0.025°. Reasonable agreement was also found when using the Fit2D package (Hammersley et al., 1996). To perform the comparison in Fit2D we first used the `tilt' command to find the detector tilt and the beam centre, and then used the `calibrant' command to find the sampletodetector distance. For Fit2D, the differences in values are sampletodetector distance 0.23 mm, beam centre 0.11 mm and detector tilt 0.052°.
Within a single procedure we have demonstrated that it is possible to accurately determine not just the beam energy but also the spatial position and orientation of the twodimensional detector relative to the sample and the beam, and thus the complete diffraction geometry.
The performance of our approach in assigning or Q values to each detector pixel can be further examined by looking at the peak locations for CeO_{2} as measured from a onedimensional diffraction pattern. Each peak position is found by locating the maxima of a fitted parabola. Table 3 compares theoretical values (generated using a lattice parameter of 5.41165 Å) for a selection of CeO_{2} diffraction peaks with values generated from a single set of data from the `tilted' case at the furthest detector position. The coefficient of determination, R^{2}, for the peak locations is 1.00–1.07 × 10^{−9}. Across 29 data sets at the furthest detector position using ten diffraction rings, the average absolute percentage error in calculated peak position was 0.00238%, the standard deviation of the absolute error was 0.00186% and the maximum absolute error was 0.0124%. In other words, the majority of peak positions were calculated to be within 1 in 25 000 of the theoretical value, whilst in the worst case an error of 1 in 8000 occurred.
6. Conclusion
Our method of calibrating a twodimensional detector for use in Xray diffraction has a number of benefits. Provided multiple rings can be recorded from various positions of the detector along the beam path, the experimenter is free to choose whatever energy they desire. This energy is easily determined as part of the calibration procedure. There is no longer a requirement to know either the sampletodetector distance or the beam wavelength, prior to calibration, as is the case with other calibration procedures such as Fit2D (Hammersley et al., 1996) or Nika (Ilavsky, 2012). By collecting data from a range of detector positions and accounting for arbitrary detector tilt, we have provided a method for accurately determining all calibration parameters.
This calibration procedure is based upon our geometrical description for the intersection of a cone's axis and its elliptical conic section.
APPENDIX A
Sources of error
There are a number of factors that can affect the accuracy of the calibration procedure, some random and some systematic. Note, however, that we are working in the superresolution paradigm where we are determining parameters at a scale that is up to two orders of magnitude better than the instrument resolution (given by the energy bandwidth of the monochromator, the sample thickness and the angular resolution of the detector).
The thermal stability of the sample and detector environment can be controlled; however, it is normal for some thermal variation to occur within an experimental hutch. This thermal variation affects not just the expansion of the calibration standard but also the displacement of the stage upon which the detector rests. The entire instrumental setup including the steel support structures is subject to
The translation stage also exhibits a systematic error due to a variation in flatness of the floormounted guide rails, which thus creates a variation in tilt and displacement at the level of the detector. An additional random error can be attributed to mechanical runout.Geometrical attenuation effects can also be considered, in particular the attenuation of the Xrays passing through the standard and the attenuation that occurs within the scintillating plane of the flatpanel detector. Both effects lead to a systematic shift in the recorded position for the centre of mass of a diffraction peak as recorded by the detector. Such attenuation affects can be quantified and applied to the calibration procedure. These corrections should improve the accuracy of the calculated beam energy; however, it could be argued that, unless equivalent transformations to correct for attenuation can be applied to experimental data for all Q, these corrections should not be applied to the calibration procedure. For the purpose of facilitating highquality data reduction we feel it is better to allow the calibration parameters to take values that compensate for the attenuated diffraction geometry, assuming that the experimental data cannot be mapped into this attenuationfree `perfect' space.
A1. Error due to of the calibration standard
At room temperature the coefficient of linear ^{−5} K^{−1}. This translates to a thermal dependence in the calculated energy of approximately 0.8 eV K^{−1}. The thermal stability of the experimental hutch at I12 is approximately K. This thermal dependence can be considered as an uncertainty in the calibration procedure, or if the temperature of the standard is known at the time of data collection, then the lattice parameter of the standard can be suitably scaled.
for cerium dioxide is approximately 10A2. Error due to detector stage displacement
The detector stage displacement is an important parameter for the calibration procedure, as the diffraction patterns recorded at the various detector positions define the geometry of the diffraction cones. By using a Renishaw XL80 interferometer we are able to determine the displacement of the translation stage more accurately. Compared to the nominal displacements used in the main body of this article, we found a systematic error in displacement of up to 125 µm in size dependent upon translation stage position. By performing repeated movements of the translation stage we were able to determine that runout and and Fig. 5] by approximately 30%. However, by applying this correction for the pseudorandom detector position, underlying drifts in the residuals (ordered by increasing semimajor axis length) grouped by diffraction image become apparent.
of the translation assembly due to the motor contribute to an uncertainty in the stage position of up to µm, which results in a eV uncertainty in the calculated energy. Correcting for the nominal sampletodetector distances in a representative calibration calculation reduced the spread in residuals for the ellipse semimajor axis [see equation (17)These drifts can be removed by considering Xray attenuation of the sample and the detector's scintillating surface, as presented in the next section.
A3. Error due to Xray attenuation
Xray attenuation at the sample and the detector's scintillator shift the effective position of a diffraction ring. In the case of sample attenuation (in a planar arrangement), an integral over all possible diffraction paths results in a shift in the centre of mass of the diffraction peak towards lower Q. A similar, but smaller, centre of mass shift occurs in the scintillating volume, due to the absorption of photons arriving at an angle to the detector surface. For an orthogonal detector with scintillator thickness w_{scin}, the corresponding change in spatial peak position, , at θ can be approximated by
where
and is the
of the scintillator.For an orthogonal detector and a planar sample of thickness w_{sam}, the change in peak position due to attenuation at the sample, , at θ can be approximated by
where
and is the
of the sample.At 88 keV correcting for attenuation at the sample ( cm^{−1}) and the detector ( cm^{−1}) reduces the calculated beam energy by 1.3 eV. Note that the effect of the increasing shift in diffraction peak location towards lower Q, with increasing diffraction angle, can be considered equivalent to diffraction from a point source onto a curved virtual detector surface/scintillator of infinitesimal thickness.
A4. Other sources of error
This list of errors is by no means exhaustive but is representative of the types of corrections that can be implemented to improve the calibration procedure. Other sources of error include the flatness of the scintillator surface, the
of the pixels, Xray attenuation at the scintillator for a tilted detector and energy drift of the monochromator. Errors can be segregated in terms of whether they are directly associated with our technique or with the process of attaining calibration data.In our measurements we experienced experimental errors of 3 eV, which we believe can be attributed to monochromator instability. The monochromator used does not have active thermal stabilization and the `orthogonal' and `tilted' data sets were collected two hours apart. We believe that energy drift from the monochromator can account for the majority of the energy difference. The monochromator was configured to use the Si 111 reflection. Using Bragg's law and , with eV and E = 88 keV, gives 7.7 ×10^{1} µrad.
Footnotes
^{1}For the purpose of this article a diffraction ring can be circular or elliptical. A diffraction ring recorded upon an orthogonal detector will be circular in shape; however, if the detector is tilted then the diffraction ring will be elliptical in shape.
References
Cervellino, A., Giannini, C., Guagliardi, A. & Ladisa, M. (2006). J. Appl. Cryst. 39, 745–748. Web of Science CrossRef CAS IUCr Journals
Cervellino, A., Giannini, C., Guagliardi, A. & Ladisa, M. (2008). J. Appl. Cryst. 41, 701–704. Web of Science CrossRef CAS IUCr Journals
Dandelin, P. G. (1822). Mémoire sur Quelques Propriétés Remarquables de la Focale Parabolique, p. 4. Brussels: Nabu Press.
Daniels, J. E. & Drakopoulos, M. (2009). J. Synchrotron Rad. 16, 463–468. Web of Science CrossRef CAS IUCr Journals
De Villiers, M. (2007). Teach. Learn. Math. 4, 25–27.
Evans, A., Jones, N., Rugg, D., Lindley, T. & Dye, D. (2012). J. Nucl. Mater. 424, 123–131. Web of Science CrossRef CAS
Hammersley, A., Svensson, S. O., Hanfland, M., Fitch, A. N. & Hausermann, D. (1996). High Pressure Res. 14, 235–248. CrossRef Web of Science
Hanan, J., Ustundag, E. & Almer, J. D. (2004). Adv. Xray Anal. 47, 174–180. CAS
Hart, M. L. & Drakopoulos, D. (2013). J. Synchrotron Rad. In preparation.
Heiney, P. A. (2005). IUCr Commission on Powder Diffraction Newsletter, No. 32, pp. 9–11.
Helliwell, J. R. (1982). Nucl. Instrum. Methods, 201, 153–174. CrossRef CAS
Helliwell, J. R. (2004). International Tables for Crystallography, Vol. C, edited by E. Prince, ch. 2.2, pp. 26–41. Dordrecht: Kluwer Academic Publishers.
Hinrichsen, B., Dinnebier, R. E. & Jansen, M. (2008). Powder Diffraction. Theory and Practice, edited by R. E. Dinnebier & S. J. L. Billinge, ch. 14, pp. 414–439. Cambridge: RSC Publishing.
Hofmann, F., Abbey, B., Connor, L., Baimpas, N., Song, X., Keegan, S. & Korsunsky, A. M. (2012). J. Mater. Res. 2, 192–199.
Hong, X., Chen, Z. & Duffy, T. S. (2012). Rev. Sci. Instrum. 83, 063901. Web of Science CrossRef PubMed
Ilavsky, J. (2012). J. Appl. Cryst. 45, 324–328. Web of Science CrossRef CAS IUCr Journals
Kumar, A. (2006). PhD thesis, Austin College, Sherman, Texas, USA.
Lenz, D. D. & Ayres, T. R. (1992). Publ. Astron. Soc. Pac. 104, 1104–1106. CrossRef Web of Science
Norby, P. (1997). J. Appl. Cryst. 30, 21–30. CrossRef CAS Web of Science IUCr Journals
Rajiv, P., Hinrichsen, B., Dinnebier, R. & Jansen, M. (2007). Powder Diffr. 22, 3–19. Web of Science CrossRef CAS
Riekel, C. (2000). Rep. Prog. Phys. 63, 233–262. Web of Science CrossRef CAS
Sui, T., Landini, G. & Korsunsky, A. (2010). AIP Conf. Proc. 1394, 113–124.
This is an openaccess article distributed under the terms of the Creative Commons Attribution (CCBY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.