research papers
Converting three-space matrices to equivalent six-space matrices for Delone scalars in S6
aRonin Institute, 9515 NE 137th Street, Kirkland, WA 98034-1820, USA, bRonin Institute, c/o NSLS-II, Brookhaven National Laboratory, Upton, NY 11973, USA, and cLawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
*Correspondence e-mail: lawrence.andrews@ronininstitute.org
The transformations from the primitive cells of the centered Bravais lattices to the corresponding centered cells have conventionally been listed as three-by-three matrices that transform three-space lattice vectors. Using those three-by-three matrices when working in the six-dimensional space of lattices represented as Selling scalars as used in Delone (Delaunay) reduction, one could transform to the three-space representation, apply the three-by-three matrices and then back-transform to the six-space representation, but it is much simpler to have the equivalent six-by-six matrices and apply them directly. The general form of the transformation from the three-space matrix to the corresponding matrix operating on Selling scalars (expressed in space S6) is derived, and the particular S6matrices for the centered Delone types are listed. (Note: in his later publications, Boris Delaunay used the Russian version of his surname, Delone.)
Keywords: Delaunay; Delone; centering transformations; centered lattices; reduced cells; lattice centering; Niggli; Selling; matrix transformations.
1. Introduction
The transformations from the primitive cells of the centered Bravais lattices to the centered cells and between alternative unit cells have conventionally been listed as matrices that are applied to three-space lattice vectors (Burzlaff & Zimmermann, 1985; Burzlaff et al., 1992). However, for both the major cell reductions [Niggli (1928) and Delone (1933)], it is convenient to work in a higher-dimension space than E3, as reported by Andrews & Bernstein (1988) for G6 reduction and Andrews et al. (2019) for S6. Therefore, as we did for G6 (Andrews & Bernstein, 1988), we need to provide the mathematically equivalent six-by-six matrices for centering in S6. This reduces the need to convert repeatedly from S6 into three-space vectors, transform the three-space vectors, and then transform back into S6. We derive the general form and list the particular matrices for converting from the 24 canonical Delone types to centered lattices in S6.
2. Background and notation
2.1. The space
Andrews et al. (2019) introduced the space S6 as an alternative representation of crystallographic lattices. The space is defined in terms of the `Selling scalars' used in Selling reduction (Selling, 1874) and by Delone (1933) for the classification of lattices. A point s in S6 is defined by
where d = −a − b − c.
2.2. The space
A crystallographic a, b, c, α, β, γ], but, when presenting common operations on unit cells, it is convenient to express each of the cell edges as a vector in the three-dimensional space of real numbers E3 (also written as R3) with each cell expressed as a 3 × 3 matrix of real numbers, i.e. as an element of E3 × E3 (see Burzlaff et al., 1992). An issue with this approach is that we should get the same crystallographic after any proper rotation of E3, i.e. by any unitary matrix of determinant +1. Such proper rotation matrices form the Lie group SO(3). Therefore, formally we should treat any matrix representation of a cell c ∈ E3 × E3 as equivalent to rc, for all r ∈ SO(3), and work in the space of (E3 × E3)/SO(3). We call this space of equivalence classes E3×3.
is commonly represented as three cell edge lengths and three angles, [Because the matrices that multiply cells represented in E3×3 are indistinguishable from ordinary 3 × 3 matrices, we will designate them as ME3, with the understanding that they may be applied in either space E3 or space E3×3.
The convention in E3 is to use the cell edges as the basis vectors of the space. There are infinitely many choices of the orientation to form the basis. Currently, it is the common convention to orient one edge vector along the x axis etc. in a right-handed setting. The convention in S6 is to use unit vectors [100000], [010000], ….
2.3. The method for deriving a transformation in one space from a transformation in another
Consider two spaces X and Y with one invertible conversion MYX mapping
and a not necessarily invertible mapping
and a transformation of X into X
We compose the mappings and the transformation to define a new transformation U of Y into Y
If Y is a finite-dimensional linear and U is linear, then we can represent U as a matrix [https://en.wikipedia.org/wiki/Linear_map] by choosing an appropriate basis. Section A in the supporting information considers the linearity in more detail.
If X is in E3×3 and Y is in S6, we can map E3×3 to S6,
Because MXY is invariant under rotation, there are infinitely many choices for the inverse. We can choose, for example,
where
which is applicable for a reasonable set of valid S6 cells. This would then allow a similar demonstration to that given in Section A of the supporting information with the mapping from Y to X being the simple square root that a linear T generates a linear U in this more complex but similar case. The details are left as an exercise for the reader.
The point is that, because U is linear, the components of its representation as a matrix can be determined by applying it to basis vectors each with only one non-zero component, letting X be in E3×3, Y be in S6 and MXY be E3toS6 (see Section 3.1).
3. Converting an matrix to an matrix
If we represent the cell as an S6 vector (Andrews et al., 2019), we can define an operator E3toS6 where E3toS6(a, b, c) = [b · c, a · c, a · b, a · d, b · d, c · d], where d = −a − b − c.
We form a matrix operating in E3×3, ME3 = [[m1,1, m1,2, m1,3], [m2,1, m2,2, m2,3], [m3,1, m3,2, m3,3]]. We need to compute a 6 × 6 matrix, MS6, to operate on S6 vectors.
The obvious basis vectors for S6, [1, 0, 0, 0, 0, 0], [0, 1, 0, 0, 0, 0], [0, 0, 1, 0, 0, 0], [0, 0, 0, 1, 0, 0], [0, 0, 0, 0, 1, 0] and [0, 0, 0, 0, 0, 1], do not correspond to real vectors in E3, since a dot product of 1 for real non-zero unit basis vectors would imply an angle between them of zero, i.e. that they are identical, but if two unit basis E3 vectors, say a and b, are identical, and one, say c, is perpendicular to both a and b, then the d = −a − b − c vector cannot be perpendicular to a or b, because a · d = b · d = −a · a − a · b − a · c = −2a · a, which cannot be zero. Therefore we use the negatives of those S6 basis vectors.
The E3×3 basis vectors we choose are shown in Table 1 with the corresponding S6 vectors.
|
3.1. Relationship to
We use the operator E3toS6 that converts a vector in E3×3 to one in S6 (see above). In our case, we are starting from reduced unit cells, which means that in S6 all six scalars are zero or negative. We choose the S6 basis vectors to have zero scalars except for a single −1 in each. In E3×3 we choose an orthogonal set (see above), where for each E3×3 vector E3toS6 produces only the corresponding S6 basis vector.
As an illustrative example, we choose the first basis vector [above, and matrix E in step (b) in Fig. 1],
We apply ME3 to that vector [step (a) in Fig. 1] and the corresponding MS6 to the corresponding S6 vector. Because only one element of the S6 basis vector is non-zero, the result of multiplying by MS6 produces only the elements of the corresponding column of MS6, with the other elements being zero [step (c) in Fig. 1]. When we multiply that E3×3 [step (b) in Fig. 1] basis vector by ME3 and then convert to S6 using E3toS6 [step (d) in Fig. 1], the resulting elements of S6 are the same S6 column elements expressed in terms of the elements of ME3.
In each case, only one of the Selling scalars will be −1 and the others will be 0 [step (c) in Fig. 1]. Because S6 is invariant under rotations of E3, we could have used any unit vector on E3 in place of [1, 0, 0], and we would have obtained the same set of S6 basis vectors.
The computer algebra system Maxima (Version 5.36.1; Chou & Schelter, 1986; https://maxima.sourceforge.net) was used to generate the following equations.
For simplicity, we show the definition of the first row of the S6 matrix (the complete matrix is listed in Section B of the supporting information). For any three-by-three matrix, ME3, the equation below computes the first column of the negative of the complete matrix (the supporting information has all the columns):
4. Conversion of reduced primitive lattices to centered lattices
Tables 2 and 3 list the matrices (Burzlaff & Zimmermann, 1985) for converting from primitive to standard centered lattices, computed from the above derivations. The designations of the 24 Delone types are slight modifications of the symbols of Delone (1933) to more modern forms. His cubic lattices are changed from `K' to `C', tetragonal from `Q' to `T' and triclinic from `T' to `A'. For each lattice, the E3×3 matrix is listed, followed on the next line by the corresponding S6 matrix.
|
|
Burzlaff & Zimmermann (1985) renumbered the lattice types of Delone (1933). For example, the cubic lattices in Delone (1933) are K1, K3 and K5. In the reports by Burzlaff & Zimmermann (1985) and Burzlaff et al. (1992), they are listed as K1, K2 and K3. Here they are listed as C1, C3 and C5. The full enumeration of the types is shown in Fig. 2. It is important to note that Burzlaff et al. (1992) showed the matrices as the transposes of the corresponding matrices of Burzlaff & Zimmermann (1985). We have chosen to start from the earlier paper. The MS6 matrices produced are then applied to the left of the S6 vectors. The International Tables for Crystallography (Burzlaff et al., 2016) use the same convention as Burzlaff & Zimmermann (1985).
5. Summary
This paper is a reference for researchers who need a method that applies to the space S6, but for which only the matrices applicable to the edge vectors of the are available. In addition, we have provided a list of the matrices required for conversion of primitive cells in S6 to the more standard centered presentations.
6. Availability of code
The C++ code for distance calculations in S6 is available at github.com, both at https://github.com/duck10/LatticeRepLib.git and https://github.com/yayahjb/ncdist.git For E3toS6 see LatticeRepLib/MatS6.cpp.
Supporting information
Discussion of linearity and complete definition of the transformation matrix. DOI: https://doi.org/10.1107/S2053273319014542/ae5074sup1.pdf
Acknowledgements
Careful copy editing and corrections by Frances C. Bernstein are gratefully acknowledged. Our thanks to Jean Jakoncic and Alexei Soares for helpful conversations and access to data and facilities at Brookhaven National Laboratory.
Funding information
Funding for this research was provided in part by: Dectris Ltd, US Department of Energy Offices of Biological and Environmental Research and of Basic Energy Sciences (grant No. DE-AC02-98CH10886; grant No. E-SC0012704); US National Institutes of Health (grant No. P41RR012408; grant No. P41GM103473; grant No. P41GM111244; grant No. R01GM117126; grant No. P30GM133893; grant No. R21GM129570).
References
Andrews, L. C. & Bernstein, H. J. (1988). Acta Cryst. A44, 1009–1018. CrossRef CAS Web of Science IUCr Journals Google Scholar
Andrews, L. C., Bernstein, H. J. & Sauter, N. K. (2019). Acta Cryst. A75, 115–120. Web of Science CrossRef IUCr Journals Google Scholar
Burzlaff, H., Grimmer, H., Gruber, B., deWolff, P. & Zimmermann, H. (2016). Crystal Lattices. International Tables for Crystallography, Vol. A, ch. 3.1, pp. 698–718. Chester: International Union of Crystallography. Google Scholar
Burzlaff, H. & Zimmermann, H. (1985). Z. Kristallogr. 170, 247–262. CrossRef CAS Web of Science Google Scholar
Burzlaff, H., Zimmermann, H. & de Wolff, P. M. (1992). Crystal Lattices. International Tables for Crystallography, Vol. A, ch. 9, pp. 734–744. Dordrecht: Kluwer Academic Publishers. Google Scholar
Chou, S.-C. & Schelter, W. F. (1986). J. Autom. Reasoning, 2, 253–273. CrossRef Google Scholar
Delone, B. N. (1933). Z. Krystallogr. 84, 109–149. Google Scholar
Niggli, P. (1928). Krystallographische und Strukturtheoretische Grundbegriffe. Handbuch der Experimentalphysik, Vol. 7, Part 1. Leipzig: Akademische Verlagsgesellschaft. Google Scholar
Selling, E. (1874). J. Reine Angew. Math. (Crelle's J.), 1874(77), 143–229. Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.