research papers
Multigrain indexing of unknown multiphase materials
^{a}Department of Physics, Technical University of Denmark, Building 307, Fysikvej, Kongens Lyngby, 2800, Denmark
^{*}Correspondence email: hfpo@fysik.dtu.dk
A multigrain indexing algorithm for use with samples comprising an arbitrary number of known or unknown phases is presented. No a priori crystallographic knowledge is required. The algorithm applies to data acquired with a monochromatic beam and a conventional twodimensional detector for diffraction. Initially, candidate grains are found by searching for crystallographic planes, using a Dirac comb convoluted with a box function as a filter. Next, candidate grains are validated and the is optimized. The algorithm is validated by simulations. Simulations of 500 cementite grains and ∼100 reflections per grain resulted in 99.2% of all grains being indexed correctly and 99.5% of the reflections becoming associated with the right grain. Simulations with 200 grains associated with four mineral phases and 50–700 reflections per grain resulted in 99.9% of all grains being indexed correctly and 99.9% of the reflections becoming associated with the right grain. The main limitation is in terms of overlap of diffraction spots and computing time. Potential areas of use include threedimensional grain mapping, structural solution and studies of complex samples, and studies of dilute phases.
Keywords: indexing; multigrain crystallography; threedimensional Xray diffraction microscopy; 3DXRD.
1. Introduction
Multigrain crystallography is a relatively new approach (Lauridsen et al., 2001; Poulsen, 2004; Sørensen, Schmidt et al., 2012) whose aim is to provide a crystallographic description of each grain within a polycrystal or a powder specimen. The technique is complementary to traditional crystallographic analysis based on either single crystals or averaging over an ensemble of grains. The experimental setup is in the simplest case identical to that typically used in singlecrystal Xray crystallography, with a monochromatic beam, a fully illuminated sample in transmission geometry on a rotary table and a twodimensional detector. The images acquired during a rotation of the sample may comprise up to a million diffraction spots from the grains simultaneously illuminated. A key step in the analysis of such data is a multigrain indexing program. Provided spot overlap is not excessive, programs have been developed that can index up to 3000 grains simultaneously (Lauridsen et al., 2001; Wright, 2005; Ludwig et al., 2009; Moscicki et al., 2009; Bernier et al., 2011; Sharma et al., 2012; Schmidt, 2014). Once grains have been indexed, all the tools of singlecrystal crystallography can be exploited for analysis of each of the grains. As examples, we mention mapping (Jakobsen et al., 2007; Wejdemann et al., 2013), defect studies (Ungár et al., 2010), and in particular the solution and of each grain (Schmidt et al., 2003; Sørensen, Schmidt et al., 2012). As illustrated, for example, in work on the compound [Cu(C_{2}O_{2}H_{3})_{2}]·H_{2}O, the can be on a par with singlecrystal results and clearly superior to results from stateoftheart powder diffraction (Vaughan et al., 2004).
Provided diffraction data are acquired with a high spatial resolution twodimensional camera close to the sample, the multigrain indexing routines can also be used to generate threedimensional maps of grains, their orientations and their stresses. This is enabled by tomographic type reconstruction algorithms, similar to classic computed tomography scanning but with diffraction contrast replacing absorption contrast. This has led to the establishment of the techniques known as threedimensional Xray diffraction microscopy (3DXRD; Poulsen, 2004; Hefferan et al., 2009) and diffraction contrast tomography (Ludwig et al., 2009; Reischig et al., 2013). Typically the sample dimensions are ∼1 mm, while the spatial resolution is 1–3 µm. By stitching subvolumes together, maps with up to 20 000 grains have been assembled. As examples of applications, see Offerman et al. (2002), Schmidt et al. (2008), King et al. (2008), Aydıner et al. (2009), Oddershede et al. (2012) and Hefferan et al. (2012).
A main limitation of the previous work has been that it applies almost exclusively to monophase materials. Furthermore, the indexing programs above all assume the et al., 2011; Sørensen, Hakim et al., 2012), but this still requires the phases to be known. To our knowledge there have only been two proposals for dealing with unknown phases, both summarized by Sørensen, Schmidt et al. (2012). In the first study the diffraction data from five grains with different unitcell parameters and orthorhombic or monoclinic symmetry were superposed, and a fastFouriertransformbased approach was used to index them. In the second study, 12 crystals of an `unknown' monophase compound with a unitcell volume of 2942 Å^{3} were successfully indexed by rotating two copies of the same data sets with respect to each other and searching for resonances (see also Schmidt, 2014).
and of the material to be known. One straightforward way to generalize the previous work is to apply the multigrain indexing and/or grain mapping algorithms repeatedly, once for each phase (JimenezMeleroGenerally speaking, the
of any structure is defined by three lattice vectors. Hence, a search and optimization procedure in the ninedimensional space spanned by these three lattice vectors will provide an indexing of all grains. However, to our knowledge, such an approach is computationally not feasible.In this article, we report a multigrain indexing routine that involves searching and optimization in a threedimensional space. This is computationally feasible, and the work presented is in fact performed with MATLAB (The MathWorks Inc., Natick, MA, USA) code on a singlecore PC. The algorithm can be applied to an arbitrary ensemble of grains belonging to an arbitrary number of known or unknown phases. Notably, the efficiency does not depend on the number of phases. The only limitations are in terms of overlap of diffraction spots and the signaltonoise ratio of the reflections. We demonstrate the method using fullscale simulations.
2. Algorithm
The indexing procedure below is based only on the position of the reflections, not their intensity. There are no a priori assumptions except for a minimum and maximum length of diffraction vectors given by parameters q_{min} and q_{max}, respectively. For ease of presentation, in the following we assume an Xray diffraction setup with a monochromatic beam, a rotary table and a farfield twodimensional detector (implying that the sample dimensions are negligible or comparable to the detector pixel size). Furthermore, we assume the diffraction spots are already harvested and represented as diffraction vectors (reflections) in a coordinate system fixed to the sample.
The algorithm that we have developed essentially indexes grains independently. It is inspired by the DIRAX algorithm (Duisenberg, 1992) for indexing of single crystals in the presence of outliers. The concept is first to search for sets of equidistant lattice planes in the full set of reflections. Such sets are represented by the direction of the plane normal, , a unit vector, and by the distance between adjacent planes, d* (see Fig. 1). The number of experimentally determined reflections `lying on such planes' can be counted. Candidate grains are defined by local maxima in the number of reflections on such planes as determined by a threedimensional search in and d*. Candidate grains are associated with a subset of the reflections – those `lying on the planes'. This subset is refined by defining a new direction of the lattice plane normal and applying the above procedure to the subset already found.
In a second step a candidate grain is validated or rejected. If validated, a set of basis vectors in
is determined, and the reflections corresponding to the grain are identified within the full set of reflections. In a third step, the is optimized by techniques similar to those well known from conventional crystallography.Once a grain has been identified the associated reflections are removed from the pool of all reflections. The entire procedure is then repeated. The indexing algorithm stops if the number of remaining (nonindexed) reflections becomes very low, or if no progress in terms of number of grains indexed has been reached within a certain number of iterations.
Once the grains have been indexed, multigrain studies can exploit the tools of singlecrystal crystallography. This is outside the scope of this article; we refer the reader to the review by Sørensen, Schmidt et al. (2012).
The main steps of the algorithm are now described in detail.
2.1. Identifying candidate grains
A number N_{u1} of directions are randomly chosen in sample space (see Fig. 2). These are candidate lattice plane normals. For each the reflections in the entire data set are projected onto the line defined by . We then search for the onedimensional lattice that fits most of these projected reflections. This is done by introducing a filter comprising a regular array of box functions (mathematically speaking a Dirac comb – also known as an impulse train function – convoluted with a box function; see Fig. 3). Let d* be the distance between the centers of the boxes, and 2∊ be the width of each box, with d* > ∊. ∊ is fixed and will typically be slightly larger than the experimental centerofmass errors of the position of the reflections in sample space. At the expense of data analysis speed, this filter is used as it is more selective than a classical Fourier transform.
The parameter d*, the lattice plane spacing, is increased from a value d*_{min} (chosen to be smaller than q_{min}) to d*_{max} = q_{max}/2 in increasingly larger steps in such a way that the onedimensional lattice point furthest from the origin (and closer than q_{max}) only moves a distance of ∊ in each step, ensuring that no lattice points are missed by the counting. For each d*, we now count the number of reflections within the boxes and subtract a similar count rate for the case of the reflections being randomly placed in sample space. The optimal value of d* is defined as the value resulting in the highest number of counts.
Among the N_{u1} candidate lattice plane normals with corresponding optimal d* values, the ten candidates with the highest counts are kept for further investigation. For each candidate, the direction and the d* value are then further optimized by a local grid search, and these optimized values are used in the next step.
The step described above tends to produce subsets of the reflections that contain (almost) all of the reflections from one particular grain but also a significant number of reflections from other grains that by chance happen to be projected such that they fall within the boxes of the comb. In order to clean up these subsets, the above procedure is repeated on the subsets corresponding to the N_{j} best candidates: N_{u2} directions (not parallel to the original direction for that subset) are chosen randomly, and the best value of d* is found for each direction. Because of the smaller number of reflections the second search can be performed with a larger value of ∊ in order to speed up the algorithm. Out of these N_{u2} searches the ten best combinations (resulting in most counts) are then chosen, and the direction and value of d* are again optimized by a local grid search. After this the best candidate is saved. Since this is done for each of the N_{j} best candidates from the first step, this results in N_{j} subsets of reflections, each predominately originating from a single grain, and the task is now to identify the lattice basis vectors for each of these candidate grains.
2.2. Indexing of candidate grains
For each of the N_{j} candidate grains found in step 1, a search is now performed for a (reciprocal) lattice basis. This is done by first listing all directions perpendicular to a plane spanned by three of the reflections (using only the 50 reflections closest to the origin, and not reflections with a cross product that is almost zero). The chosen subset of reflections is then as before projected onto each of these directions, and a search for the best onedimensional lattice is performed. In this case the possible d* values are selected from the list of projected reflections (again between d*_{min} and d*_{max}). The best value of d* is chosen as the largest d* value among those d* values corresponding to counts within 80% of the maximum count.
This results in a best d* value and the corresponding count from each of the searched directions, and from this set all d* values lower than 80% of the maximum count are removed. The remaining d* values and directions correspond to potential directlattice vectors pointing along the given direction and having length 1/d* (e.g. Giacovazzo et al., 2011). From these potential directlattice vectors the three shortest, linearly independent candidates are found.
At this point the algorithm has produced a candidate N_{j} potential grains in the sample. From each of these the vectors are calculated, and for each point in this potential the nearest reflection (closer than 5∊) is found. This search for reflections in the potential lattice is performed in the full set of reflections. The number of reflections matching the potential lattice found in this way is then used to determine which of the N_{j} candidate unit cells is best by selecting the one corresponding to the most reflections, and if the number of reflections is high enough this candidate is then chosen for further optimization; if not, a new search is performed.
(three realspace vectors) for2.3. Optimization of the unit cell
It is possible that the steps described above result in a potential ; Křivý & Gruber, 1976; GrosseKunstleve et al., 2004). As cell reduction is not a main focus point here, we chose a simple heuristic solution.
that either is smaller or larger than the real of the grain or is given by a nonstandard set of basis vectors. There are well known methods to determine a standard reduced (Gruber, 1973By trying simple linear combinations of the found basis vectors, it is checked if a N_{t} of consecutive tries where no potential with enough reflections has been found.
with a larger volume results in significantly (more than 20%) more reflections corresponding to the lattice. If this is the case this new is chosen over the original guess. Next, the orientation and length of the three potential basis vectors are optimized, and the reflections corresponding to the found (candidate) grain are removed from the full set of reflections to obtain a new set of reflections which is then used as the full set of reflections, and the steps describes above are repeated. This continues until there have been a number3. Simulations
The aim of the simulation is to prove the concept and to provide an understanding of the accuracy and robustness of the algorithm.
The diffraction data that served as input to the simulations were generated by the program PolyXSim (Sørensen, 2006). The Xray energy was defined to be 50 keV. Reflections were harvested and a normally distributed noise was added to each component of the position of each reflection.
Two studies were performed (each comprising six identical simulations): one with 500 grains of cementite, and one with 200 grains of four minerals typically found in granite (50 grains of each). The materials chosen reflect potential use for studies of minority phases in steel and for geological specimens, respectively. Furthermore, the minerals exhibit different crystal symmetries ranging from trigonal to monoclinic. The unitcell parameters for the five materials are given in Table 1. All grains were assumed to be fully illuminated and of the same size, sufficiently large that signaltonoise issues can be neglected. The parameters for the simulations are shown in Table 2. σ is the standard deviation of the noise added to each of the components of the position in of the various reflections.


As a measure of the success of the simulations a grain was defined as successfully identified if the unitcell volume was within 1% of the nominal value and if the fraction of correctly identified reflections was above 0.9.
The results are shown in Table 3. For each material class (cementite and granite) the average results from ten simulations are presented. The table shows that almost all grains are successfully identified, with an accurate determination of the unitcell volume, a very high level of average completeness and a very low level of falsely attributed reflections.

Notably, the figures of merit for the multiphase granite simulation in Table 3 are even better than for the singlephase cementite simulation. The data analysis speed was also much faster. This illustrates the fact that grains are indexed independently, and hence the figure of merit and running times are assumed to be approximately independent of the number of phases, but strongly dependent on the total number of grains.
It is also noteworthy that the program identified the grains of the four different phases in the granite simulation in the order of decreasing number of reflections per grain, that is plagioclase, orthoclase, biotite, quartz. This happens to be in order of decreasing unitcell volume and increasing crystal symmetry.
4. Discussion
4.1. Potential applications
The output of the above program – the list of grains and associated reflections – can immediately be used in connection with a variety of existing singlecrystal and multicrystal analysis programs. An overview of the status of (singlephase) multigrain crystallography was provided recently by Sørensen, Schmidt et al. (2012). Here examples are provided of resulting threedimensional maps of grains, orientations and stresses in materials science and geoscience, as well as the use of multigrain structure solution and in smallmolecule drug design and studies. Included is also preliminary work on proteins. The fact that the algorithm presented above requires neither a priori information nor that a certain number of grains is present for each phase generalizes all of these applications to studies of multiphase materials with unknown space groups and an arbitrary number of grains of each phase.
In addition we point to the following prospective uses:
(a) Minority phases. Using powder diffraction it is often difficult to observe phases with a smaller than 1%. In contrast, with multigrain methods, it is realistic to detect volume fractions of 10^{−6}, provided the grains of the minority phase are sufficiently large to give rise to detectable diffraction spots (Poulsen et al., 2001). Implemented with a suitable sample translation, it becomes possible to screen for diffraction signals from particles associated with parts per million concentrations.
(b) Local diffraction. With the upcoming synchrotron nanobeam beamlines it becomes possible to use scanning procedures to make threedimensional maps of a set of nanoscale grains. As scanning methods tend to be slow, for larger samples it will be natural to map only the grains within an intrinsic gauge volume, which may be a small fraction of the total sample volume. The diffraction signal will in such cases be dominated by spurious diffraction spots moving in and out of the beam during the rotation. Analogous to the case of the DIRAX algorithm for use in singlecrystal diffractometry with an obstinate list of reflections, the current algorithm is seen as being robust towards such outliers.
(c) Total crystallography. Sørensen, Schmidt et al. (2012) defined this concept as the simultaneous characterization of the threedimensional atomic and the threedimensional grainscale structure of polycrystalline samples with unknown phase(s), as well as the temporal characteristics of such samples. In other words total crystallography is a (hypothetical) method that allows one to study an arbitrary polycrystal and for each grain characterize both its atomic structure and its mesocale structure (the threedimensional shape, orientation and stress state). In our view the current work demonstrates the feasibility of total crystallography. Specifically, in a 3DXRD setup, farfield and nearfield twodimensional detectors can be used simultaneously (Poulsen, 2012; Hefferan et al., 2012). It is therefore possible first to identify and index the grains using the farfield data and the method described here and then for each phase to use the nearfield data and existing programs to map grains in three dimensions.
The generalization of the above approach to monochromatic or timeofflight neutron multicrystal diffraction is seen as straightforward.
4.2. Limitations
The indexing algorithm described here does not use any a priori information, and as such we anticipate that known and unknown phases will be handled equally well. Likewise, crystal symmetry is not used directly, and as such we anticipate that, everything else being equal, figures of merits will be independent of crystal symmetry. These hypotheses are corroborated by the simulations performed. Instead, applications will be restricted by the following inherent limitations (which apply to all multigrain indexing approaches):
(a) Spot overlap. The probability of spot overlap on the detector is determined by the number of grains illuminated, the texture of the sample, the size of the and the orientation spread of each grain. Experience from singlephase indexing methods shows that for typical samples of relevance to hard materials science one can aim at indexing up to a few thousand grains simultaneously, while for mediumsized crystal structures such as crambin, indexing of up to 100 grains is feasible (Sørensen, Schmidt et al., 2012).
(b) Grain sizes. Typically grain volumes vary by several orders of magnitude. As a consequence, all reflections associated with large grains may be above the intensity threshold, all reflections for small grains may be below, while a fraction of the diffraction spots are identified for grains of intermediate size. For the latter grains an iterative approach may be needed.
An additional limitation of the current implementation is the speed of data processing. Indexing of the 500 cementite grains above took 5 d, while the 200 mineral grains took 1.5 d (using MATLAB and a singlecore PC). We anticipate that the speed could be greatly increased by the use of a different programming language, GPUs and parallel computing, but still the approach may be too slow, for example, for online analysis of more challenging samples. There are a variety of ways to increase the speed, involving, for example, a priori information (some phases are known), sorting the data (searching first for the larger grains using only the more intense spots) or gaining additional information on grain position by acquiring data on an additional semitransparent nearfield detector. However, this is outside the scope of this article, where the aim has been to demonstrate the feasibility of multigrain indexing of multiphase samples.
5. Conclusion
Conventionally, Xray crystallography is based on two extreme sample morphologies: perfect single crystals and homogeneous powders. The multigrain methods developed over the past decade enable us to treat polycrystalline specimens as an ensemble of individual crystals, creating the possibility to rigorously characterize such samples in terms not just of average properties but of the distributions of those properties. The work in this article generalizes the previous work in the direction of multiphase materials. It also points to the feasibility of total crystallography: the synthesis of methodologies for threedimensional grain mapping (mesoscale structure) and structure solution and
(atomic scale). More specifically, for a fixed number of grains and reflections the resulting figure of merit and efficiency of the algorithm presented do not depend on the number of phases. Likewise, the figure of merit and data processing speed are not strongly dependent on the group symmetry.Acknowledgements
The authors are grateful to Gavin Vaughan, Jon Wright and Henning Osholm Sørensen for advice. We acknowledge funding from the ERC grant `Transmission Xray microscopy' and the Danish instrument center for users of large scale neutron and Xray facilities, DANSCATT. We thank the anonymous referees for constructive advice.
References
Aydıner, C. C., Bernier, J. V., Clausen, B., Lienert, U., Tomé, C. & Brown, D. W. (2009). Phys. Rev. B, 80, 024113. Google Scholar
Bernier, J. V., Barton, N. R., Lienert, U. & Miller, M. P. (2011). J. Strain Anal. Eng. Des. 46, 527–547. Web of Science CrossRef Google Scholar
Duisenberg, A. J. M. (1992). J. Appl. Cryst. 25, 92–96. CrossRef CAS Web of Science IUCr Journals Google Scholar
Giacovazzo, C., Monaco, H. L., Artioli, G., Viterbo, D., Milanesio, M. & Ferraris, G. (2011). Fundamentals of Crystallography, IUCr Texts on Crystallography 15. Oxford University Press. Google Scholar
GrosseKunstleve, R. W., Sauter, N. K. & Adams, P. D. (2004). Acta Cryst. A60, 1–6. Web of Science CrossRef CAS IUCr Journals Google Scholar
Gruber, B. (1973). Acta Cryst. A29, 433–440. CrossRef IUCr Journals Web of Science Google Scholar
Hefferan, C. M., Li, S. F., Lind, J., Lienert, U., Rollett, A. D., Wynblatt, P. & Suter, R. M. (2009). Comput. Mater. Continua, 14, 209–219. Google Scholar
Hefferan, C. M., Lind, J., Li, S. F., Lienert, U., Rollett, A. D. & Suter, R. M. (2012). Acta Mater. 60, 4311–4318. Web of Science CrossRef CAS Google Scholar
Jakobsen, B., Poulsen, H. F., Lienert, U. & Pantleon, W. (2007). Acta Mater. 55, 3421–3430. Web of Science CrossRef CAS Google Scholar
JimenezMelero, E., van Dijk, N. H., Zhao, L., Sietsma, J., Wright, J. P. & van der Zwaag, S. (2011). Mater. Sci. Eng. A, 528, 6407–6416. CAS Google Scholar
King, A., Johnson, G., Engelberg, D., Ludwig, W. & Marrow, J. (2008). Science, 321, 382–385. Web of Science CrossRef PubMed CAS Google Scholar
Křivý, I. & Gruber, B. (1976). Acta Cryst. A32, 297–298. CrossRef IUCr Journals Web of Science Google Scholar
Lauridsen, E. M., Schmidt, S., Suter, R. M. & Poulsen, H. F. (2001). J. Appl. Cryst. 34, 744–750. Web of Science CrossRef CAS IUCr Journals Google Scholar
Ludwig, W., Reischig, P., King, A., Herbig, M., Lauridsen, E. M., Johnson, G., Marrow, T. J. & Buffière, J. Y. (2009). Rev. Sci. Instrum. 80, 033905. Web of Science CrossRef PubMed Google Scholar
Moscicki, M., Kenesei, P., Wright, J., Pinto, H., Lippmann, T., Borbély, A. & Pyzalla, A. R. (2009). Mater. Sci. Eng. A, 524, 64–68. Web of Science CrossRef Google Scholar
Oddershede, J., Camin, B., Schmidt, S., Mikkelsen, L. P., Sørensen, H. O., Lienert, U., Poulsen, H. F. & Reimers, W. (2012). Acta Mater. 60, 3570–3580. Web of Science CrossRef CAS Google Scholar
Offerman, S. E., van Dijk, N. H., Sietsma, J., Grigull, S., Lauridsen, E. M., Margulies, L., Poulsen, H. F., Rekveldt, M. T. & van der Zwaag, S. (2002). Science, 298, 1003–1005. Web of Science CrossRef PubMed CAS Google Scholar
Poulsen, H. F. (2004). ThreeDimensional Xray Diffraction Microscopy. Berlin: Springer. Google Scholar
Poulsen, H. F. (2012). J. Appl. Cryst. 45, 1084–1097. Web of Science CrossRef CAS IUCr Journals Google Scholar
Poulsen, H. F., Juul Jensen, D., Tschentscher, T., Wcislak, L., Lauridsen, E. M., Margulies, L. & Schmidt, S. (2001). Textures Microstruct. 35, 39–54. CrossRef Google Scholar
Reischig, P., King, A., Nervo, L., Viganó, N., Guilhem, Y., Palenstijn, W. J., Batenburg, K. J., Preuss, M. & Ludwig, W. (2013). J. Appl. Cryst. 46, 297–311. Web of Science CrossRef CAS IUCr Journals Google Scholar
Schmidt, S. (2014). J. Appl. Cryst. 47, 276–284. Web of Science CrossRef CAS IUCr Journals Google Scholar
Schmidt, S., Olsen, U. L., Poulsen, H. F., Sørensen, H. O., Lauridsen, E. M., Margulies, L., Maurice, C. & Juul Jensen, D. (2008). Scr. Mater. 59, 491–494. Web of Science CrossRef CAS Google Scholar
Schmidt, S., Poulsen, H. F. & Vaughan, G. B. M. (2003). J. Appl. Cryst. 36, 326–332. Web of Science CrossRef CAS IUCr Journals Google Scholar
Sharma, H., Huizenga, R. M. & Offerman, S. E. (2012). J. Appl. Cryst. 45, 705–718. Web of Science CrossRef CAS IUCr Journals Google Scholar
Sørensen, H. O. (2006). PolyXSim, http://sourceforge.net/apps/trac/fable/wiki/PolyXSim. Google Scholar
Sørensen, H. O., Hakim, S., Pedersen, S., Christiansen, B. C., Balogh, Z. I., Hem, C. P., Pasarin, I. S., Schmidt, S., Olsen, U. L., Oddershede, J., Frandsen, C., Feidenhans'l, R. K. & Stipp, S. L. S. (2012). Can. Mineral. 50, 501–509. Google Scholar
Sørensen, H. O., Schmidt, S., Wright, J. P., Vaughan, G. B. M., Techert, S., Garman, E. F., Oddershede, J., Davaasambuu, J., Paithankar, K. S., Gundlach, C. & Poulsen, H. F. (2012). Z. Kristallogr. 227, 63–78. Google Scholar
Ungár, T., Ribárik, G., Balogh, L., Salem, A. A., Semiatin, S. L. & Vaughan, G. B. M. (2010). Scr. Mater. 63, 69–72. Google Scholar
Vaughan, G. B. M., Schmidt, S. & Poulsen, H. F. (2004). Z. Kristallogr. 219, 813–825. Web of Science CSD CrossRef CAS Google Scholar
Wejdemann, C., Poulsen, H. F., Lienert, U. & Pantleon, W. (2013). JOM, 65, 35–43. CrossRef CAS Google Scholar
Wright, J. P. (2005). ImageD11, http://sourceforge.net/projects/fable/files/ImageD11/. 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.