research papers
Computation of diffuse scattering arising from onephonon excitations in a neutron timeofflight singlecrystal Laue diffraction experiment
^{a}ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, UK, ^{b}Thomas Young Centre, London Centre for Nanotechnology and Department of Chemistry, University College London, London WC1E 6BT, UK, ^{c}Department of Physics, Royal Holloway, University of London, Egham TW20 0EX, UK, and ^{d}Deutsches Elektronen Synchrotron DESY, Notkestrasse 85, Hamburg 22603, Germany
^{*}Correspondence email: matthias.gutmann@stfc.ac.uk
Direct phonon excitation in a neutron timeofflight singlecrystal Laue diffraction experiment has been observed in a single crystal of NaCl. At room temperature both phonon emission and excitation leave characteristic features in the diffuse scattering and these are well reproduced using ab initio phonons from density functional theory (DFT). A measurement at 20 K illustrates the effect of thermal population of the phonons, leaving the features corresponding to phonon excitation and strongly suppressing the phonon annihilation. A recipe is given to compute these effects combining DFT results with the geometry of the neutron experiment.
1. Introduction
The most common type of diffuse scattering at roomtemperature in a singlecrystal diffraction experiment designed to measure it is caused by lattice vibrations and is known as thermal diffuse scattering (TDS) (Born & Lonsdale, 1942). This can be seen using laboratory or synchrotron Xray sources as well as neutrons at largescale facilities (Osborn & Welberry, 1990; Welberry & Goossens, 2008). In principle, such data can be used to reconstruct phonon dispersion curves at least for simple crystal structures and using model potentials (Holt et al., 1999). Diffuse scattering was among the early techniques used to obtain phonon dispersion curves along highsymmetry directions in elemental crystals before inelastic neutron and Xray experiments became routine (Olmer, 1948; Placzek & Van Hove, 1954, Squires, 1978; Burkel, 2000). Inelastic techniques using single crystals provide complete momentum and energyresolved dispersion curves, whilst diffraction only has momentum resolution. Interesting effects were noted by Willis et al. in the case of timeofflight (TOF) Laue neutron diffraction where the neutron can accidentally match the wavevector and energy of phonons leading to a splitting of the thermal diffuse scattering features corresponding to phonon emission and absorption (Willis et al., 1986; Schofield & Willis, 1987). This effect has not been observed with Xrays owing to the high energy of the order of several keV compared with up to hundreds of meV for lattice vibrations but in principle takes place as well. Willis derived a formalism to extract the sound velocity of a crystal using the splitting of the diffuse lines of the acoustic phonons which dominate the TDS pattern, at least close to the Bragg peaks (Willis, 1986). A few experiments were conducted to demonstrate the feasibility, and good agreement was obtained for the sound velocity obtained by other means (Willis, 1986; Carlile & Willis, 1989; Carlile et al., 1992). Inelastic effects have also been noticed in quartz but were not considered further (Tucker et al., 2001). It should be noted though that very few observations of direct phonon excitation in neutron diffraction experiments have been reported to date, limiting somewhat the use of this technique for molecular crystals in particular, with benzil being the only known example to the best of our knowledge (Welberry et al., 2003).
Nowadays, density functional theory (DFT) has established itself as a tool of choice to predict phonon dispersion curves from a et al., 2005). These computations also give the phonon eigenvectors for arbitrary q points which are not easily accessible experimentally though, in principle, they are available from the intensity of the TDS or inelastic peaks provided the inelastic differs from zero. Examples are now emerging where DFT computations are used to aid in the interpretation of diffuse scattering experiments (Bosak et al., 2009; Wehinger et al., 2014; Gutmann et al., 2013).
and is routinely used in the planning and interpretation of inelastic experiments (MitchellIn this paper, we combine DFT calculations with the geometry of the Laue neutron diffraction experiment to compute the `inelastic diffraction' associated with phonons. This is compared with diffuse scattering calculated in the instrument geometry using the quasistatic approximation, which is when the phonon energy is negligible compared with the energy of the scattering probe. Our motivation is to be able to identify and reproduce the diffuse scattering arising from onephonon scattering over a large region in ab initio calculations, in order to aid in the interpretation of diffuse scattering patterns beyond the Monte Carlo modelling approach (Welberry, 2004).
in combination with2. Experimental
2.1. Diffraction experiments
NaCl single crystals with a purity of 99.99% were purchased from Sigma–Aldrich; hereafter, they are referred to as synthetic. A large, irregularly shaped crystal was shaped into a sphere of 6 mm diameter. This sample is usually used as a calibration standard to ensure reproducibility between ISIS experiment cycles. Natural specimens were obtained in Hallstatt (Austria) from the local salt mine. They were of various colour, indicating impurities, but there were also colourless, optically clear crystals. An optically clear, colourless crystal of approximate dimensions 2 × 4 × 5 mm was extracted. However, the purity is not known explicitly and hence data from the synthetic crystal are used throughout the paper. The natural sample was only used to verify that effects similar to those described below occur and hence are sample independent, although they may vary in magnitude owing to impurity effects on the phonons (Caldwell, 1967). Neutron diffraction data were collected at room temperature using the SXD singlecrystal diffractometer at the ISIS spallation neutron source (Oxfordshire, UK; Keen et al., 2006) on various occasions on the synthetic crystal with exposures ranging between four hours and up to two days at single settings or a series of crystal settings about 20° apart. This crystal was also run at 20 K. The natural specimen was exposed for 12 h in various orientations but only at room temperature. Reciprocalspace volumes from these data were obtained using SXD2001 (Gutmann, 2005). A natural crystal was also used for a complementary highenergy Xray diffraction experiment on the beamline BW5 installed at the storage ring DORIS III at DESY in Hamburg, Germany. These data merely confirmed the absence of inelastic excitation effects seen in the neutron data.
2.2. Computation of phonon dispersions using DFT
Electronic structure calculations were performed using the planewave pseudopotential methods as implemented in the CASTEP code (Clark et al., 2005). The Perdew–Burke–Ernzerhof (Perdew et al., 1996) generalized gradient approximation functional was used for both the planewave calculation and the generation of pseudopotentials. These were of the optimized normconserving variety (Rappe et al., 1990). Energies and forces were well converged at a planewave cutoff of 600 eV. A sampling of 8 × 8 × 8 q points (60 points when symmetry reduced) was found to be sufficient to converge energy and atomic forces below 2.2 × 10^{−8} eV per ion and 1.0 × 10^{−3} eV Å^{−1}, respectively. Sufficient selfconsistent cycles were performed to achieve a convergence tolerance of 1.0 × 10^{−10} eV per atom. Geometry optimizations used the Broyden–Fletcher–Goldfarb–Shanno algorithm with a force tolerance of 1.0 × 10^{−3} eV Å^{−1}. Phonon dispersion calculations were performed on the resulting minimumenergy structures via diagonalization of dynamical matrices computed using density functional perturbation theory (DFPT) and linearresponse methods (Refson et al., 2006) on a total of 60 q points. Because our intention was to enable the program to use output from either CASTEP or VASP/PHONOPY, similar computations were also carried out with VASP (Kresse & Hafner, 1994), which equally well reproduced the diffraction features. Throughout the paper, the CASTEP calculations are used.
3. Computation of the thermal diffuse scattering
We first establish a coordinate system in the neutron diffractometer based on spherical coordinates. The location of each pixel on the detector is characterized by two angles, δ and ν, which correspond to the longitude and latitude, respectively, as illustrated in Fig. 1. The incident neutron beam travels along , which is along the positive y direction, and is scattered into a pixel located at the endpoint of . Expressed in the lefthanded orthonormal coordinate system (x, y, z) shown in Fig. 1, the incident and scattered wavevectors read
where λ_{i} and λ_{f} are the wavelengths (Å) of the incident and scattered neutrons, respectively. In the case of an event these wavelengths have the same magnitude, but they generally differ in the case of Table 1 lists the δ and ν angles for the 11 detector modules of the SXD instrument, and they are illustrated graphically in Fig. 2 for an NaCl data set at room temperature.

The angular ranges for the δ and ν angles in equation (1) for the various detectors of the SXD instrument are also given in Table 1. For a single crystal, the scattering vector is
where [ω], [χ] and [ϕ] denote rotation matrices corresponding to the angle settings of the goniometer stage, [UB] is the orientation matrix of the crystal, and h, k, l are which need not be integer valued (Busing & Levy, 1967). For phonons it is convenient to decompose h, k, l into an integer part, where Bragg reflection may occur, and a fractional part such that the latter reflects the wavevector of the phonon in the first Brillouin zone.
In TOF neutron scattering, every pixel has a timeofflight range associated with it, as a fixed histogram with a given number of time channels and a fixed minimum and maximum time. More recently eventmode has become available, where neutron events are time stamped as they occur without imposing a histogram. The SXD instrument at ISIS uses the fixed histogrammode. The timeofflight recorded corresponds to the total travel time of the neutron from the source to the detector pixel and is measured in units of microseconds. The total travel time t_{tot} can be decomposed into the time from the source to the sample, t_{i}, and the time from the sample to the pixel, t_{f}. These components can be related to the wavelengths of the incident and scattered neutron:
where h is the Planck constant, m_{n} is the neutron mass, L_{1} is the primary flight path (mm) from source to sample and L_{2} is the secondary flight path (mm) from the sample to the pixel. The energy (meV) of the neutron is related to the wavelength as follows:
For calculating the TDS in the instrument geometry, every pixel and TOF channel is converted to E_{f} for a given E_{i} and the scattering vector which is transformed back to h, k, l using equation (1) and decomposed into an integer and fractional h, k, l in the first in order to look up the eigenvectors and phonon frequencies at that point in the DFT output. We note that the DFT calculations sometimes use a compared with the conventional crystallographic cell and hence a transformation between these two cells should also be applied when converting to h, k, l of the DFT cell.
For calculating firstorder TDS it is convenient to define a onephonon ):
(Xu & Chiang, 2005The summation extends over all n_{atoms} atoms in the Further, m_{k} is the mass of atom k, b_{k} is its coherent scattering length, u_{k}^{2} is the meansquare atomic displacement, taken as isotropic for the case of NaCl, is the complexvalued eigenvector of atom k in the phonon mode j at wavevector in the first with being the nearest Bragg point, and is the position of atom k in the Here, refers to the primitive used by the DFT program. It should be noted that the form of equation (1) depends on whether the dynamical matrix in the DFT program is defined such as to give periodic eigenvectors, as is the case in CASTEP (Refson et al., 2006), or otherwise, as is the case in PHONOPY (Togo et al., 2008). The formula given here applies to CASTEP. For the case of nonperiodic eigenvectors see work by Xu & Chiang (2005).
The TDS intensity can be written as follows assuming energy and momentum conservation are fulfilled:
with C = N/2 and N being the number of atoms in the sample. This constant is treated as a scale factor and is set to 1 for convenience. ω_{j} denotes the frequency of the phonon mode j, and `+' and `−' refer to phonon creation and annihilation, respectively. The energy integration of the doubledifferential [equations (6) and (7)] is effectively carried out here numerically using many wavelengths and counting the number of phonons falling into a given Q bin and normalizing to this. The phonon occupancy at temperature T is given by
where ℏ is the reduced Planck constant and k_{B} is the Boltzmann constant.
In the quasistatic approximation, (6) and (7) are summed over all n_{modes} phonon modes to yield
The constant is treated again as a scale factor and set to one. Initially a master mapping from pixel coordinates and TOF to a reciprocalspace map is established, assuming only elastic events. This master mapping is essentially assigning a given pixel and TOF channel to a voxel in the reciprocalspace volume. The recipe to compute inelastic diffraction can be summarized as follows:
(1) Select a fixed wavelength λ_{i} and a fixed energy E_{i}, and .
(2) Compute for all pixels and TOF of the instrument λ_{f}, E_{f} and
(3) Compute as well as and use the experimental UB matrix for a given experimental run to derive h, k, l.
(4) Transform h, k, l from the conventional crystallographic cell to the DFT cell, if necessary, and decompose it into an integer part and a fractional part in the first Brillouin zone.
(5) Using the DFT phonons interpolated on a fine grid over the first h, k, l there is a matching phonon frequency corresponding to ω within a given tolerance in meV, and compute the TDS intensity using equations (6) or (7) as appropriate and equation (8). Add this to the reciprocalspace voxel using the master map.
check whether for the given(6) Repeat until all wavelengths are covered.
This recipe has been coded in OpenCL to enable parallel computation using multicore CPUs and GPUs.
4. Results
Sections through the (0kl) plane are shown in Fig. 3(a) employing data from the highangle, 90° and lowangle banks. With reference to Fig. 2, which illustrates the detector numbering used in this paper, the highangle banks are detectors 1, 6 and 10, the 90° banks are detectors 2, 5, 7, 9 and 11, and the lowangle banks are detectors 3, 4 and 8, respectively.
Diffuse scattering is concentrated close to the Bragg spots and takes the appearance of a butterfly shape. However, comparing the data from the various detectors, it becomes apparent that the diffuse scattering in the highangle detectors near the Bragg spots is geometrically distorted rather than symmetric as would be expected. Weak arcs join or can be thought of as joining various neighbouring reflections on either side of the Bragg peak and these arcs are absent in the data from the lowangle banks. It should be noted that the indexing is taken as determined by the software. The assignment of which h, k or l merely corresponds to a permutation of the symbols and is not relevant for comparing the various data sets. Sections were chosen to provide as large a coverage as possible. More complete volumetric data sets are available in the form of movies in the supplementary information.
isNaCl phonons have been determined experimentally (Raunio et al., 1969), albeit over a somewhat limited and can quite easily be computed using stateoftheart density functional codes such as CASTEP and VASP. NaCl is often included as an example tutorial for learning these programs, with the resulting phonon dispersion curves closely matching the experimental ones. Owing to the limited comparison and the fact that DFT phonons are calculated at 0 K with no adjustable parameters, only qualitative agreement would be expected in our calculations below. In the absence of direct phonon excitations in the quasistatic treatment, the thermal diffuse scattering can be calculated as the sum over all the phonon branches at a given q point using the formulas (5) and (9). It should be noted that in all calculations only the diffuse scattering is computed. Bragg peaks are not included, as is common practice when modelling diffuse scattering (Welberry, 2004).
The presence of direct phonon excitations in the data requires a different treatment. When creating the reciprocalspace maps, the underlying assumption is that the scattering events are purely elastic, i.e. no occurs and the initial and final momenta of the scattered neutron have the same magnitude. Neither of these conditions is fulfilled when takes place, which is governed by energy and momentum conservation.
The computations were carried out in the detector geometry of the instrument and then mapped to crystallographic . Of the order of 5600 incident wavelengths were used from 560 meV down to 1 meV covering the timeofflight range measured in the detector. The calculations were performed on Intel Xeon E52687 CPUs comprising effectively 32 cores. Despite parallelization of the codes, the calculations were rather formidable. Computing a volumetric data set of 300 × 300 × 300 points took 24 h per detector module and 5600 wavelengths. Three modules were computed in parallel and hence it would take of the order of four days to compute these volumetric data sets for all 11 detector modules and one orientation/temperature of the crystal. For comparison the TDS patterns took ∼100 s per detector module. At this point, the codes have not been fully optimized for speed or memory and further gains are possible. For each wavelength, the phonon spectrum was calculated for cases of both up and down scattering and all resulting maps were summed. Only inelastic events arising from exciting phonons from the ground state and vice versa were considered. Events such as between occupied phonon states, multiphonon scattering and higher harmonics, i.e. any processes involving more than one phonon, were not included and did not appear essential as their contribution is thought to diminish rapidly. The resulting reciprocalspace maps were not convoluted with the instrument resolution. The timeofflight peak shapes are typically asymmetric along the timeofflight direction, leading to a tail at the base of the peak that points towards the origin of giving the peaks a teardrop shape in the reciprocalspace maps. The footprint of the Bragg peaks is indicated by contours in Figs. 3, 4 and 5 and also the supplementary movies. These were derived from the threedimensional profile fitting routine in SXD2001 after transforming them to and taking a 1% contour level from the maximum of each peak. The contours were elongated in the radial direction to better capture the asymmetric tail. Missing contours on some Bragg peaks are due to the peaks being too close to the detector edge and the integration not succeeding on them. However, they can be easily inferred from the neighbouring Bragg peaks. Trial runs were initially carried out to obtain a good value for the energy tolerance of 1 meV. It can be seen that the results obtained using both the standard TDS formalism and our method of explicitly taking into account direct phonon excitations resemble the data recorded with the lowangle bank, corresponding to detector 3 in Figs. 3(a)–3(c). However, for the 90° [detectors 2 and 9 in Figs. 3(a)–3(c)] and backscattering data (detectors 1 and 10) using the latter formalism results in a qualitatively better agreement compared with the TDS formalism with the observed diffuse scattering. A more complete illustration of the inelastic effects as they occur throughout and away from the (0kl) layer is contained in the movies in the supplementary information. It was already noted by Willis (1986) that inelastic effects would best be seen in backscattering.
as a volumetric data set as outlined in §3To further ascertain the phonon origin of the scattering leading to the splitting of diffuse lines, two more experiments were conducted. In one, the same synthetic specimen was measured at room temperature in a different orientation by rotating the crystal around a vertical axis from ω = −150 to +150° and the diffuse scattering calculated accordingly using the full phonon dispersions and the quasistatic approach. The results are shown in Figs. 3(d)–3(f). An enlarged portion around the 040 reflection is shown in Fig. 4. Again, the diffuse scattering computed using the inelastic treatment shows better agreement. This illustrates that the diffuse features change appearance in the presence of inelastic effects as the crystal is rotated, whilst no changes are expected when using the quasistatic approximation. The former is expected, because upon rotating the crystal a given feature in is probed with neutrons of a different wavelength and hence energy. Intuitively, this can be understood from Bragg's law: given a constant d spacing, when changing the scattering angle θ, the wavelength has to change. A similar effect, also reported by Welberry et al. (2003), is seen by comparing Figs. 3(a) and 5(a) in their paper. The angular dependence of this effect close to the Bragg peaks is discussed rather more formally by Schofield & Willis (1987). If the diffuse scattering changes appearance and diffuse lines appear split when measured in different orientations, this observation can be taken as an indication that these diffuse features originate from excitations arising in the sample rather than static lattice distortions. The second experiment again used the synthetic crystal in a somewhat different orientation compared with the roomtemperature measurement, but this time the crystal was cooled to 20 K and the calculations carried out using this temperature and an energy tolerance of 1 meV. It is expected that this would change the population of the phonons leading to a more pronounced asymmetry in the intensity distribution of the features corresponding to the phonon excitation and annihilation. The results are shown in Fig. 5. The differences between the calculated patterns using the TDS and our approach in the (0kl) layer (Figs. 5a–5c) are a bit more subtle and from a comparison with the experimental data it is not clear which would be preferred. In comparison, the (4kl) layer contains stronger diffuse features and again the inelastic treatment is clearly favoured. The intensity scales are the same as in Fig. 3. This allows assignment of the arclike features on the higher and lowerQ side of a Bragg peak to phonon emission and absorption, respectively. The second supplementary movie shows a sidebyside comparison of the 20 K data and calculations for a volumetric data set. Finally, we note that the inelastic effects break the symmetry of the diffraction pattern, and hence averaging of the diffuse scattering using the Laue symmetry of the as is commonly done in singlecrystal diffuse scattering experiments, should be avoided in this case.
5. Summary and conclusions
In this paper, we have illustrated the effect of inelastic excitations that can occur in a timeofflight neutron singlecrystal Laue experiment. The signal arising from onephonon excitations takes the form of split arcs originating from Bragg reflections, and this is most pronounced in the data from the backscattering banks whilst being virtually absent in the lowangle detectors. The consequence is the breaking of the symmetry of the diffraction pattern, which means that such data should not be symmetrized using the Laue symmetry as is commonly done. Furthermore, recording of the diffuse scattering at a much lower temperature results in a pronounced asymmetry in the intensities of the diffuse split arcs on the high and lowQ side of Bragg reflections due to the change in the phonon population. The 1/ω factor in equations (6) and (7) makes the acoustic modes the dominant contribution close to the Bragg reflections. Further away there is a dependence on the details of the phonon dispersions of a particular material, e.g. whether branches cross and the relative alignment of the eigenvectors and Q vector. Close to the h00, 0k0 and 00l reflections these would be the longitudinal acoustic modes allowing the derivation of the sound velocity as reported by Willis et al. (1986). The effects have been computed combining the TDS formulas for the with phonon eigenvectors and frequencies as well as the temperature, and simulating the experiment in the actual instrument geometry using many wavelengths as a discretized representation of the incident white beam. It has been shown that in the case of neutron diffraction at a pulsed source the quasistatic approach is occasionally insufficient and the new approach presented in this paper provides a better way to simulate such features, despite the phonons being at the origin in both cases, i.e. the underlying physical cause of the diffuse scattering is the same. Whilst such effects may be considered a nuisance when modelling diffuse scattering using a Monte Carlo approach (Welberry, 2004), they can help in augmenting the information and classifying diffuse features with respect to their dynamic or static origin (Welberry et al., 2003). If ab initio phonons, or magnons in the case of magnetic diffuse scattering, are available these can be used to aid in the interpretation of such features and contrasting them with diffuse scattering arising from other types of disorder. In the case of magnons, polarized neutrons can further establish their magnetic origin (Brückel & Schweika, 2002). As more instruments similar to SXD are now operational at other spallation neutron sources such as TOPAZ at the Spallation Neutron Source (Oak Ridge National Laboratory, USA) or SENJU at JPARC (Ibaraki Prefecture, Japan), it is anticipated that such inelastic features will increasingly be observed (Schultz et al., 2014; Tamura et al., 2012). It can also be used complementarily to neutron powder inelastic experiments or optical spectroscopies such as IR, UV–Vis and Raman. The work presented in this paper can serve as a starting point to address a number of interesting topics arising from such effects, such as to what extent one can derive phonon dispersion curves, possibly taking into account several crystal orientations. This does not appear straightforward, as there is no simple geometric mapping between features before and after rotating the crystal. The conditions for phonon excitations are different, and the patterns are a superposition of phonon excitation curves measured with many different wavelengths and, in addition, are dependent on the instrument geometry. Another point concerns the TDS correction on Bragg intensities: Schofield & Willis (1987) noted that in the presence of inelastic excitation effects under favourable conditions phonon excitations may be forbidden around certain reflections, leading to Bragg intensities completely free of TDS, but details of the influence in the presence of such effects in general have not been worked out to date. Finally, as further extensions of this work one may envisage using such effects to fit either model potentials (Holt et al., 1999) or coupling this with DFT to derive better model potentials, possibly including other kinds of spectroscopic data as well. From a theory point of view, DFT has the advantage of being applicable to different kinds of materials by choice of the appropriate functional without the need to fit many parameters to a model potential. The computational effort in the current work for our method and the DFT calculations on top of this seem rather discouraging for the time being but this may change sometime in the future with growing computing power.
Supporting information
Quicktime movie scanning layerbylayer through https://doi.org/10.1107/S1600576715010912/ks5457sup2.mov
comparing the 20 K data measured on the SXD diffractometer on the left with the corresponding computed section using the inelastic excitation formalism on the right. DOI:Quicktime movie scanning layerbylayer through https://doi.org/10.1107/S1600576715010912/ks5457sup1.mov
comparing the roomtemperature data measured on the SXD diffractometer on the left with the corresponding computed section using the inelastic excitation formalism on the right. DOI:Acknowledgements
We would like to thank Drs T. G. Perring and R. Ewings for helpful discussions. We thank the Scientific Computing Department of the Science and Technology Facilities Council for continued access to the High Performance Computing facility SCARF cluster at the Rutherford Appleton Laboratory. Some of the computations have been performed on the UK Supercomputing Facility HECToR under the auspices of the Materials Chemistry Consortium (EPSRC grant EP/D504872). GG is grateful to the UK Science and Technology Facilities Council for financial support in the form of a NEXTDTC studentship under contract No. 4070006444. We are grateful to the London Centre for Nanotechnology, UCL Research Computing, and the eScience facility of the STFC for the computational resources.
References
Born, M. & Lonsdale, K. (1941). Nature, 150, 490. CrossRef Google Scholar
Bosak, A., Hoesch, M., Krisch, M., Chernyshov, D., Pattison, P., SchulzeBriese, C., Winkler, B., Milman, V., Refson, K., Antonangeli, D. & Farber, D. (2009). Phys. Rev. Lett. 103, 076403. Web of Science CrossRef PubMed Google Scholar
Brückel, T. & Schweika, W. (2002). Polarized Neutron Scattering. Forschungzentrum Jülich GmbH. Google Scholar
Burkel, E. (2000). Rep. Prog. Phys. 63, 171–232. Web of Science CrossRef CAS Google Scholar
Busing, W. R. & Levy, H. A. (1967). Acta Cryst. 22, 457–464. CrossRef IUCr Journals Web of Science Google Scholar
Caldwell, R. F. & Klein, M. V. (1967). Phys. Rev. 158, 851–875. CrossRef CAS Web of Science Google Scholar
Carlile, C. J., Keen, D. A., Wilson, C. C. & Willis, B. T. M. (1992). Acta Cryst. A48, 826–829. CrossRef CAS Web of Science IUCr Journals Google Scholar
Carlile, C. J. & Willis, B. T. M. (1989). Acta Cryst. A45, 708–715. CrossRef CAS Web of Science IUCr Journals Google Scholar
Clark, S. J., Segall, M. D., Pickard, C. J., Hasnip, P. J., Probert, M. I. J., Refson, K. & Payne, M. C. (2005). Z. Kristallogr. 220, 567–570. Web of Science CrossRef CAS Google Scholar
Gutmann, M. J. (2005). SXD2001. ISIS Facility, Rutherford Appleton Laboratory, Oxfordshire, UK. Google Scholar
Gutmann, M. J., Refson, K., von Zimmermann, M., Swainson, I. P., Dabkowski, A. & Dabkowska, H. (2013). J. Phys. Condens. Matter, 25, 315402. Web of Science CrossRef PubMed Google Scholar
Holt, M., Wu, Z., Hong, H., Zschack, P., Jemian, P., Tischler, J., Chen, H. & Chiang, T.C. (1999). Phys. Rev. Lett. 83, 3317–3319. Web of Science CrossRef CAS Google Scholar
Keen, D. A., Gutmann, M. J. & Wilson, C. C. (2006). J. Appl. Cryst. 39, 714–722. Web of Science CrossRef CAS IUCr Journals Google Scholar
Kresse, G. & Hafner, J. (1994). Phys. Rev. B, 49, 14251–14269. CrossRef CAS Web of Science Google Scholar
Mitchell, P. C. H., Parker, S. F., RamirezCuesta, A. J. & Tomkinson, J. (2005). Vibrational Spectroscopy with Neutrons. Singapore: World Scientific. Google Scholar
Olmer, P. (1948). Acta Cryst. 1, 57–63. CrossRef CAS IUCr Journals Web of Science Google Scholar
Osborn, J. C. & Welberry, T. R. (1990). J. Appl. Cryst. 23, 476–484. CrossRef Web of Science IUCr Journals Google Scholar
Perdew, J. P., Burke, K. & Ernzerhof, M. (1996). Phys. Rev. Lett. 77, 3865–3868. CrossRef PubMed CAS Web of Science Google Scholar
Placzek, G. & Van Hove, L. (1954). Phys. Rev. 93, 1207–1214. CrossRef CAS Web of Science Google Scholar
Rappe, A. M., Rabe, K. M., Kaxiras, E. & Joannopoulos, D. (1990). Phys. Rev. B, 41, 1227–1230. CrossRef Web of Science Google Scholar
Raunio, G., Almqvist, L. & Stedman, R. (1969). Phys. Rev. 178, 1496–1501. CrossRef CAS Web of Science Google Scholar
Refson, K., Tulip, P. R. & Clark, S. J. (2006). Phys. Rev. B, 73, 155114. Web of Science CrossRef Google Scholar
Schofield, P. & Willis, B. T. M. (1987). Acta Cryst. A43, 803–809. CrossRef CAS Web of Science IUCr Journals Google Scholar
Schultz, A. J., Jørgensen, M. R. V., Wang, X., Mikkelson, R. L., Mikkelson, D. J., Lynch, V. E., Peterson, P. F., Green, M. L. & Hoffmann, C. M. (2014). J. Appl. Cryst. 47, 915–921. Web of Science CrossRef CAS IUCr Journals Google Scholar
Squires, G. (1978). Introduction to the Theory of Thermal Neutron Scattering. Cambridge University Press. Google Scholar
Tamura, I., Oikawa, K., Kawasaki, T., Ohhara, T., Kaneko, K., Kiyanagi, R., Kimura, H., Takahashi, M., Kiyotani, T., Arai, M., Noda, Y. & Ohshima, K. (2012). J. Phys. Conf. Ser. 340, 012040. CrossRef Google Scholar
Togo, A., Oba, F. & Tanaka, I. (2008). Phys. Rev. B, 78, 134106. Web of Science CrossRef Google Scholar
Tucker, M. G., Keen, D. A. & Dove, M. T. (2001). Miner. Mag. 65, 489–507. Web of Science CrossRef CAS Google Scholar
Wehinger, B., Bosak, A., Piccolboni, G., Refson, K., Chernyshov, D., Ivanov, A., Rumiantsev, A. & Krisch, M. (2014). J. Phys. Condens. Matter, 26, 115401. Web of Science CrossRef PubMed Google Scholar
Welberry, T. R. (2004). Diffuse Xray Scattering and Models of Disorder. Oxford University Press. Google Scholar
Welberry, T. R. & Goossens, D. J. (2008). Acta Cryst. A64, 23–32. Web of Science CrossRef CAS IUCr Journals Google Scholar
Welberry, T. R., Goossens, D. J., David, W. I. F., Gutmann, M. J., Bull, M. J. & Heerdegen, A. P. (2003). J. Appl. Cryst. 36, 1440–1447. Web of Science CrossRef CAS IUCr Journals Google Scholar
Willis, B. T. M. (1986). Acta Cryst. A42, 514–525. CrossRef CAS Web of Science IUCr Journals Google Scholar
Willis, B. T. M., Carlile, C. J., Ward, R. C., David, W. I. F. & Johnson, M. W. (1986). Europhys. Lett. 2, 767–774. CrossRef CAS Web of Science Google Scholar
Xu, R. & Chiang, T.C. (2005). Z. Kristallogr. 220, 1009–1016. CrossRef CAS Google Scholar
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.