Next Article in Journal
Sustainable Housing Supporting Health and Well-Being
Previous Article in Journal
Systematic Model-Based Steady State and Dynamic Optimization of Combined Cooling and Antisolvent Multistage Continuous Crystallization Processes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Mathematical Tools that Connect Different Indexing Analyses †

by
Ryoko Oishi-Tomiyasu
Institute of Mathematics for Industry, Kyushu University, 819-0395 Fukuoka, Japan
Presented at the 2nd International Electronic Conference on Crystals, 10–20 November 2020; Available online: https://iocc_2020.sciforum.net/.
Proceedings 2020, 62(1), 8; https://doi.org/10.3390/proceedings2020062008
Published: 31 December 2020
(This article belongs to the Proceedings of The 2nd International Online Conference on Crystals)

Abstract

:
As mathematical tools that can be commonly used for indexing analyses from different types of experimental patterns, we have recently developed (i) rules on forbidden hkl’s that can be used even when the space group and setting are unknown, (ii) an algorithm for error-stable Bravais lattice determination, (iii) generalization of the de Wolff figure of merit for powder diffraction (1D data) to data in higher-dimensions such as Kikuchi patterns (2D data) by electron backscatter diffraction (EBSD). In particular, (ii) could be used in a variety of situations, not just for indexing. It is explained how (i)–(iii) are used in the mathematical framework of our indexing algorithms.

1. Introduction

Mathematical tools that can be commonly used in ab-initio indexing analyses are introduced herein. They were originally invented for powder diffraction [1], and subsequently applied to the indexing of Kikuchi bands in electron backscatter diffraction (EBSD) patterns [2]. “Ab-initio” means that the indexing is carried out without any prior information on the parameters and Bravais type of the unit cell.
In the case of powder diffraction, the values of d-spacings (hence, lengths of reciprocal-lattice vectors) are obtained from positions of diffraction peaks. In the EBSD case, the orientations of reciprocal-lattice vectors are provided from the positions of Kikuchi bands. Our indexing algorithms for them use a common mathematical framework shown in Figure 1. First, the parameters of the primitive cell are determined, because (i) simple rules of systematic absence are available, if only basis vectors of a primitive lattice are considered. Subsequently, (ii) Bravais-type (and centering) determination is carried out. This process can be error-stable enough to deal with unit-cell parameters containing large errors due to zero-point shifts (powder [3]) or projection-center shifts (EBSD [4]). We also (iii) generalized the idea of the de Wolff figure of merit Mn [5], which has been the most efficient indicator in powder indexing. The generalized one presents similar properties for EBSD patterns [2]. The developed software is now available on the web: Powder auto-indexing: https://z-code.kek.jp/zrg/ (CONOGRAPH). EBSD ab-initio indexing: https://osdn.net/projects/ebsd-conograph/
In what follows, mainly the mathematics used for (ii) is discussed. Due to the limitation of the space, (i) and (iii) are only mentioned referring to published papers. The author believes that these theoretical results will be also useful in different analyses of crystallography.

1.1. Notation

We summarize the notation and symbols used in the article. The inner product of the Euclidean space N is denoted by u v , and the Euclidean norm u u is denoted by | u | 2 . Any basis v 1 , , v N  of an N-dimensional(N-D) lattice L is associated to a quadratic form:
f ( x 1 , , x N ) = | x 1 v 1 + + x N v N | 2 = x T S x ,
where x = ( x 1 , , x N ) T is a vector, and S is the symmetric matrix with v i v j in the (i,j)-entry. S is also the Gramian (or metric tensor [6]) of L. The stabilizer of S is defined as the following subgroup of G L N ( )  (=the group of integral matrices with the determinant ±1):
Stab ( S ) = { g G L N ( ) : g S g T = S } .
The Gramians S1, S2 belong to the same Bravais type, if Stab ( S 1 ) , Stab ( S 2 ) are conjugate in G L N ( ) (i.e., there exists σ G L N ( ) such that σ Stab ( S 1 ) σ 1 = Stab ( S 2 ) ) [7].
On the linear space S N consisting of N-by-N symmetric matrices, an inner product is defined by S T : = Trace ( S T ) , which makes S N the metric space (the distance between S and T equals ( S T ) ( S T ) ). The subset of S N consisting of all the positive definite symmetric matrices is denoted by S N + . The action of G L N ( ) on S N + is given by S g S g T .
The following is an overview of the lattice-basis reduction theory that discusses methods to provide the representatives for the orbits G L N ( ) \ S N + . Namely, D S N +  is the subset that fulfills the following (i), (ii):
(i)
S N + = g G L N ( ) D [ g ] ,
(ii)
D [ g 1 ] D [ g 2 ] = for any g 1 ± g 2 G L N ( ) , where D [ g ] = d e f { g S g T : S D } .
As the boundaries of D are prone to complications, overlaps of the boundary D : = D D i n ( D i n : set of interior points of D ) are frequently allowed. In such a case, D should satisfies (i) and the following (ii)’ and (iii)’:
(ii)’
D i n [ g 1 ] D i n [ g 2 ] = for any g 1 ± g 2 G L N ( ) .
(iii)’
D D [ g ] for only finitely many g G L N ( ) .
It is straightforward to see that any S in D i n satisfies Stab ( S ) = { ± 1 } . Thus, all the S with non-triclinic Bravais types belong to the boundary of D . The following are the definitions of Venkov [8] and Delaunay reductions used in Section 3; for any fixed S 0 S N + , define D S 0  by:
D S 0 : = { S S N + : S S 0 ( g S g T ) S 0   for   any   g G L N ( ) } .
From the definition, D S 0 [ g ] = D S 0   g T Stab ( S 0 ) holds. If S belongs to D S 0 , S is Venkov-reduced with regard to S 0 . In particular, S is Selling-reduced, if S belongs to D A N , where A N is the symmetric matrix with 2 in the diagonal entries and 1 in the other entries.
A N ( i , j ) = { 2 i = j , 1 i j .

2. Determination of the Primitive Lattice

For some types of systematic absence (SA), the ratio of forbidden reflections is not small, but considerably high rate. The rules of SA stated in International Tables depend on the space group and setting of atomic positions. Simple rules of SA are useful for developing algorithms that generally work, if they are available.
In order to obtain such simple rules, only basis vectors of the primitive lattice are considered herein. L * is the reciprocal lattice of the crystal lattice L. { l 1 * , l 2 * } is a primitive set, if it is a subset of some basis l 1 * , l 2 * , l 3 * of L * .
Theorem 1.
(Theorem 2, [9]) Regardless of the type of SA, there are infinitely many primitive sets { l 1 * , l 2 * } of L * such that none of l 1 * , l 2 * , l 1 * + 2 l 2 * , 2 l 1 * + l 2 * correspond to an extinct reflection due to the SA. Furthermore, there exist infinitely many 2D sublattices L 2 * of L * such that L 2 * is expanded by such l 1 * , l 2 * .
Theorem 1 is not true, if l 1 * , l 2 * , l 1 * + 2 l 2 * , 2 l 1 * + l 2 * are replaced e.g., by l 1 * , l 2 * , l 1 * + l 2 * , l 1 * l 2 * (vectors in Ito’s formula [10]: 2 ( | l 1 * | 2 + | l 2 * | 2 ) = | l 1 * + l 2 * | 2 + | l 1 * l 2 * | 2 ). The theorem assures us that some combinations of observed reflections correspond to l 1 * , l 2 * , l 1 * + 2 l 2 * , 2 l 1 * + l 2 * , for some two vectors l 1 * , l 2 * contained in a basis of L*. In the powder case, the inner product l 1 * l 2 * is computed by
l 1 * l 2 * = ( | l 1 * + 2 l 2 * | 2 | l 1 * | 2 4 | l 2 * | 2 ) / 4 = ( | 2 l 1 * + l 2 * | 2 4 | l 1 * | 2 | l 2 * | 2 ) / 4 .
Similarly, in the EBSD case, the direction l * / | l * | of the reciprocal-lattice vector l are obtained from the coordinates of Kikuchi bands. The vector-length ratio | l 1 * | : | l 2 * | : | l 1 * + 2 l 2 * | can be calculated from the directions of l 1 * / | l 1 * | , l 2 * / | l 2 * | , ( l 1 * + 2 l 2 * ) / | l 1 * + 2 l 2 * | by solving the linear equation.
( l 1 * / | l 1 * | 2 l 2 * / | l 2 * | ( l 1 * + 2 l 2 * ) / | l 1 * + 2 l 2 * | ) x = 0
In both Equations (5) and (6), the lengths (or directions) of l 1 * , l 2 * , l 1 * + 2 l 2 * are sufficient to obtain the matrix (or the ratio of its components) in Equation (7). The remaining length (or direction) of 2 l 1 * + l 2 * can be used to remove unlikely solutions quickly.
l 1 * · l 1 * l 1 * · l 2 * l 1 * · l 2 * l 2 * · l 2 *
Theorem 2 is a 3D version of Theorem 1.
Theorem 2.
(Theorem 4 in [9]) Regardless of the type of SA, there are infinitely many bases l 1 * , l 2 * , l 3 * of L * such that the following hold:
(a)
the reflections of ± l 1 * + l 2 * + l 3 * are not forbidden.
(b)
For both i = 2 ,   3 , (i) none of the reflections of m l 1 * + ( m 1 ) ( l 1 * + l i * ) are forbidden for any integer m, or (ii) none of the reflections of m l i * + ( m 1 ) ( l 1 * l i * ) are forbidden for any integer m 0 .
As a result, CONOGRAPH assigns l 1 * ± l 2 * , l 1 * ± l 3 * , l 1 * + l 2 * + l 3 * and either of l 1 * or { l 2 * , l 3 * } to various combinations of observed reflections. See [2] for the EBSD case.

3. Bravais-Lattice Determination from Unit-Cell Parameters Containing Large Observation Errors

3.1. Theoretical Background

After the parameters of the primitive cell are obtained in the indexing process, it is necessary to convert them into parameters of the conventional cell. For a Gramian matrix S o b s extracted from observed data, how can one estimate the Bravais type of the unknown true value S ^ of S o b s ? The error can be observational errors or rounding errors of floating-point numbers [11].
If S o b s is exact (i.e., S o b s = S ^ ), the symmetry group of n-by-n  S o b s can be determined e.g., by the method of [12]. However, if S o b s S ^ , no matter how close S o b s is to S ^ , g Stab ( S o b s ) is not generally true even if 1 g Stab ( S ^ ) . As a result, it is only possible to estimate likely ones as Stab ( S ^ ) . For this reason, error-stable methods have been investigated in mathematical crystallography.
This determination can be done by step 1 and 2 in Table 1 by using a finite set H0 with the following property, where D is a domain that fulfills (i), (ii)’, (iii)’ in Section 1.1.
H0: if S obs D , then S ^ g H 0 D [ g ] .
Namely, H0 is a finite set containing all g G L N ( ) such that g 1 S o b s ( g 1 ) T is nearly reduced (i.e., close to D ) for some S o b s that belongs to D ).
If D is the Niggli-reduced domain (Chap.9.2.2, [6]), G0 in Table 1 consists of 168 elements. The number m of lattice characters Lk is 42, after two triclinic cases are excluded (Table 9.2.5.1 [6]). From the definition, H0 must contain G0, hence, the computation time of the method of Table 1 is roughly estimated as |H0| × m ≥ 168 × 42 = 7056. This is a little time consuming if it is applied to multiple primitive cells generated in the indexing process.
The methods of Andrews and Bernstein [14,15] are basically same as this Niggli-reduced case, although it is not assured that their heuristics can always generate all the necessary operations (in their method, 25 operations in [16] are used to generate the elements of H0).
The use of the Delaunay reduced domain was proposed Burzlaff and Zimmermann [17,18]. This reduces the number of lattice characters from 44 to 30. However, H0 is set to {1} in their method, so it can basically handle only the exact case.
Thus, the following are the problems, in order to develop a faster and more reliable Bravais-lattice determination method.
Q1: Which reduction method minimizes the computation time for Table 1?
Q2: Under which assumption on the error size of S o b s is it possible to output all the S with S o b s S and Stab ( S ) = Stab ( S ^ ) ?
Our idea for Q1 was to use the following Venkov-reduced domain D S 0 as D in Table 1.
  • S0 = I3 (3 × 3 identity matrix).
  • S0 = A3 in Equation (4).
By the choice of S0, D S 0 include non-triclinic Lk in the interior distant from its boundary D S 0 . If S ^ is in Lk and D S 0 , S o b s S ^ is also in the interior of D S 0 . As a result, S ^ and S o b s are both reduced with regard to the same basis. In this case, it is not necessary to consider the nearly-reduced case (namely, H0 may be set to {1}).
Based on this idea, the author proved that error-stable determination is possible under the following condition C on the error size of S o b s [19] (This is an answer to Q2):
C: for any 3-by-3 symmetric matrix T and 0 v n , if S ^ T v S ^ v T / 2 , S o b s · T > 0 also holds.
Namely, C excludes only the case: v T S v / 2 S ^ · T S o b s · T 0 (If L is the crystal lattice with the Gramian S ^ , the half of the squared-length of some non-zero vector in L is observed as a non-positive value). Hence, C assumes that the error of S o b s is not extraordinary large. Under this condition, the following is proved:
Theorem 3.
(Theorems 1–4 in [19]) For a given S o b s D S 0 , assume that S ^ belongs to the Bravais type B, in addition to . In this case, S ^ belongs to the VB, a union of finitely many linear subspaces in Table 2.
After the centering types of the unit cell are classified to those in Table 2, sit is straightforward to classify them into higher-symmetric Bravais types.
Therefore, contrary to our intuition, regardless of the magnitude of the error in S o b s , it is possible to output S with S o b s S and Stab ( S ) = Stab ( S ^ ) generally, without increasing the computation time at all. However, the error of S o b s affects the distance between the output S and its true value S ^ .

3.2. Computation Results

The implemented program is used in our indexing software [1,2]. In [2], indexing analysis was carried out for EBSD patterns with large projection-center shifts as follows (z: the camera length).
Δ x z , Δ y z , Δ z z = 0 , ± 0.005 , ± 0.01 , ± 0.02
Nevertheless, Bravais-lattice determination failed only for a small number of cases among them (see Tables and Figure in [2] for more details).

4. Discussion

The theorems presented in Section 2 and Section 3 hold true for any symmetry types the crystal structures can have. Our error-stable Bravais-type determination is mathematically guaranteed, even for parameters containing large error. Although the number of operations |H0| × m ≥ 168 × 42 = 7056 cannot be decreased as long as the Niggli reduction is used, it can be reduced from 7056 to 154 (58, conditionally) by using the Venkov reduction for I3 and A3. However, other reduction methods (or S0) might be able to provide a faster method. No studies have been reported for lattices of dimensions more than 3.
As ab-initio indexing software for powder diffraction patterns, ITO [21,22], TREOR [23], and DICVOL [24] are well known. EBSD ab-initio indexing have been also studied in [25], although more accurate methods for band extraction and projection center identification are also needed for this indexing analysis.
From a theoretical point of view, the two indexing analyses have much in common. This suggests that updating the mathematical crystallography is effective in obtaining reliable and efficient analytical methods in a short time.

Funding

This study was financially supported by the PREST (JPMJPR14E6) and JST Mirai (JPMJMI18GD).

Acknowledgments

We would like to extend our gratitude to R. Taniguchi and S. E. Graiff-Zurita of Kyushu University, who helped us in coding the software, and performing the computation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Oishi-Tomiyasu, R. Robust powder auto-indexing using many peaks. J. Appl. Crystallogr. 2014, 47, 593–598. [Google Scholar] [CrossRef]
  2. Oishi-Tomiyasu, R.; Tanaka, T.; Nakagawa, J. Distribution Rules of Systematic Absence and Generalized de Wolff Figure of Merit Applied to EBSD Ab-Initio Indexing. Available online: https://arxiv.org/abs/2003.13403 (accessed on 31 March 2020).
  3. Dong, C.; Wu, F.; Chen, H. Correction of zero shift in powder diffraction patterns using the reflection-pair method. J. Appl. Crystallogr. 1999, 32, 850–853. [Google Scholar] [CrossRef]
  4. Ram, F.; Zaefferer, S.; Jäpel, T.; Raabe, D. Error analysis of the crystal orientations and disorientations obtained by the classical electron backscatter diffraction technique. J. Appl. Crystallogr. 2015, 48, 797–813. [Google Scholar] [CrossRef]
  5. de Wolff, P.M. A Simplified Criterion for the Reliability of a Powder Pattern Indexing. J. Appl. Crystallogr. 1968, 1, 108–113. [Google Scholar] [CrossRef]
  6. Hahn, T. International Tables for Crystallography Vol. A, 5th ed.; Springer: Berlin/Heidelberg, Germany, 2005. [Google Scholar]
  7. Michel, L. Symmetry, invariants, topology. iv. fundamental concepts for the study of crystal symmetry. Phys. Rep. 2001, 341, 265–336. [Google Scholar] [CrossRef]
  8. Venkov, B.A. On the reduction of positive quadratic forms. Izv. Akad. Nauk SSSR Ser. Mat. 1940, 134, 37–52. [Google Scholar]
  9. Oishi-Tomiyasu, R. Distribution rules of systematic absences on the Conway topograph and their application to powder auto-indexing. Acta Crystallogr. A 2013, 69, 603–610. [Google Scholar] [CrossRef]
  10. Ito, T. A general powder X-ray photography. Nature 1949, 164, 755–756. [Google Scholar] [CrossRef]
  11. Grosse-Kunstleve, R.W.; Sauter, N.K.; Adams, P.D. Numerically stable algorithms for the computation of reduced unit cells. Acta Crystallogr. A 2004, 60, 1–6. [Google Scholar] [CrossRef]
  12. Plesken, W.; Souvignier, B.B. Computing isometries of lattices. J. Symb. Comput. 1997, 24, 327–334. [Google Scholar] [CrossRef]
  13. Andrews, L.C.; Bernstein, H.J. Lattices and reduced cells as points in 6-space and selection of bravais lattice type by projections. Acta Crystallogr. A 1988, 44, 1009–1018. [Google Scholar] [CrossRef]
  14. Andrews, L.C.; Bernstein, H.J. The geometry of Niggli reduction: BGAOL–embedding Niggli reduction and analysis of boundaries. J. Appl. Crystallogr. 2014, 47, 346–359. [Google Scholar] [CrossRef]
  15. Gruber, B. The relationship between reduced cells in a general Bravais lattice. Acta Crystallogr. A 1973, 29, 433–440. [Google Scholar] [CrossRef]
  16. Burzlaff, H.; Zimmermann, H. On the metrical properties of lattices. Z. Krist. Cryst. Mater. 1985, 170, 247–262. [Google Scholar] [CrossRef]
  17. Zimmermann, H.; Burzlaff, H. Delos—A computer program for the determination of a unique conventional cell. Z. Krist. Cryst. Mater. 1985, 170, 241–246. [Google Scholar] [CrossRef]
  18. Oishi-Tomiyasu, R. Rapid bravais-lattice determination algorithm for lattice parameters containing large observation errors. Acta Crystallogr. Sect. A 2012, 68, 525–535. [Google Scholar] [CrossRef] [PubMed]
  19. de Wolff, P.M. Detection of simultaneous zone relations among powder diffraction lines. Acta Crystallogr. 1958, 11, 664–665. [Google Scholar] [CrossRef]
  20. Visser, J.W. A fully automatic program for finding the unit cell from powder data. J. Appl. Crystallogr. 1969, 2, 89–95. [Google Scholar] [CrossRef]
  21. Kohlbeck, F.; Horl, E.M. Trial and error indexing program for powder patterns of monoclinic substances. J. Appl. Crystallogr. 1978, 11, 60–61. [Google Scholar] [CrossRef]
  22. Boultif, A.; Louër, D. Powder pattern indexing with the dichotomy method. J. Appl. Crystallogr. 2004, 37, 724–731. [Google Scholar] [CrossRef]
  23. Dingley, D.J.; Wright, S.I. Determination of crystal phase from an electron backscatter diffraction pattern. J. Appl. Crystallogr. 2009, 42, 234–241. [Google Scholar] [CrossRef]
  24. Li, L.; Ouyang, S.; Yang, Y.; Han, M. EBSDL: A computer program for determining an unknown Bravais lattice using a single electron backscatter diffraction pattern. J. Appl. Crystallogr. 2014, 47, 1466–1468. [Google Scholar] [CrossRef]
  25. Li, L.; Han, M. Determining the Bravais lattice using a single electron backscatter diffraction pattern. J. Appl. Crystallogr. 2015, 48, 107–115. [Google Scholar] [CrossRef]
Figure 1. Common mathematical framework of our indexing algorithms.
Figure 1. Common mathematical framework of our indexing algorithms.
Proceedings 62 00008 g001
Table 1. Outline of error-stable Bravais lattice determination methods 1.
Table 1. Outline of error-stable Bravais lattice determination methods 1.
Prepared sets in codesFor a domain D that fulfills (i), (ii)’, (iii)’ in Section 1.1, and its topological closure D ¯ , let G 0 be the finite set consisting of all g G L N ( ) with D ¯ D ¯ [ g ] . For each finite group G k ( k = 1 , , m ) contained in G 0 , prepare the set of linear subspace L k consisting of all S S n with Stab ( S ) G k . (Namely, L 1 , , L m are lattice characters [6].)Finite set H0 consisting of operations g for which D [ g ] may contain S ^ when S o b s is in D .
Input parametersGramian S o b s (assume S o b s D by exchanging the basis)
Step 1For any g H 0 , if S 2 o b s = g 1 S o b s ( g 1 ) T is close to the domain D ¯ (i.e., nearly reduced), do the following; for each L k ( k = 1 , , m ), calculate S L k closest to S 2 o b s by projecting S 2 o b s on Lk. If S 2 o b s and S are close to each other, store g, S in the array for the Bravais type of L k .
Step 2Output the stored g, S after removing duplicates.
1 The same calculation can be done, even if D is replaced by a union of finitely many D [ g ] such as the Venkov reduced domain D S 0 . The only difference is that Lk may not be in the boundary of D ¯ .
Table 2. B, S0, VB in Theorem 1.
Table 2. B, S0, VB in Theorem 1.
Bravais Type BS0H0 = {1}?The Number of Linear Subspaces VB
(The Number When H0 = {1})
Primitive monoclinic I 3 Yes3 (3) … Table 3 in [19]
Face-centered orthorhombic A 3 Yes3 (3) … Table 4 in [19]
Body-centered orthorhombic 1 A 3 Yes … Table 5 in [19]
Rhombohedral A 3 conditionally yes 264 (16) … Table 6 in [19]
Base-centered monoclinic A 3 conditionally yes 269 (21) … Table 8 in [19]
1 As for the face-centered case, our method simply uses the fact that S o b s has the face-centered symmetry if and only if the inverse of S o b s has the body-centered symmetry. 2 If H0 = {1}, the number of steps can be reduced to the same number as the case when S o b s is exact. H0 = {1} also holds even for rhombohedral and base-centered case, by adding another condition to those of Theorem 3 [19].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Oishi-Tomiyasu, R. Mathematical Tools that Connect Different Indexing Analyses. Proceedings 2020, 62, 8. https://doi.org/10.3390/proceedings2020062008

AMA Style

Oishi-Tomiyasu R. Mathematical Tools that Connect Different Indexing Analyses. Proceedings. 2020; 62(1):8. https://doi.org/10.3390/proceedings2020062008

Chicago/Turabian Style

Oishi-Tomiyasu, Ryoko. 2020. "Mathematical Tools that Connect Different Indexing Analyses" Proceedings 62, no. 1: 8. https://doi.org/10.3390/proceedings2020062008

APA Style

Oishi-Tomiyasu, R. (2020). Mathematical Tools that Connect Different Indexing Analyses. Proceedings, 62(1), 8. https://doi.org/10.3390/proceedings2020062008

Article Metrics

Back to TopTop