short communications\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoFOUNDATIONS
ADVANCES
ISSN: 2053-2733

On the allowed values for the triclinic unit-cell angles

CROSSMARK_Color_square_no_text.svg

aMPL, Imperial College, London SW7 2AZ, England, and bDiamond Light Source Ltd, Harwell Science and Innovation Campus, Didcot, Oxfordshire OX11 0DE, England
*Correspondence e-mail: j.foadi@imperial.ac.uk

(Received 20 September 2010; accepted 29 October 2010; online 25 November 2010)

This short article questions and investigates the possible range of values for the three angles of a unit cell in the triclinic system. Although no constraints are reported in manuals and tables for crystallography, the three angles are not really independent; the range of allowed values is calculated and presented in this paper.

1. Introduction

There is a broad, if not universal, acceptance by the community of crystallographers that cell parameters in the triclinic system can assume any arbitrary value. This is reported everywhere in textbooks, tables and papers on crystallography. It is certainly true that the a, b and c lengths of the cell can be given any value, and that they are quantities independent of each other. The same, though, cannot be stated for the angles; although each individual angle can assume any value between 0 and 180°, they are not independent of each other. This feature of the triclinic system is strongly related to the volume of the unit cell; wrong combinations of angles can result in the volume being zero or, even, an imaginary number.

In this short paper we will examine the whole issue quantitatively. It will be readily shown that the determinant appearing in the expression for the unit-cell volume is at the heart of the problem. A factorized version of that expression will enable us to single out the set of allowed values for the cell's angles. Finally, a simple geometric construction will help in clarifying the reason why the volume expression can have a null or imaginary value.

2. The expression for the volume of a unit cell

Throughout the text we will adopt a right-handed coordinate system; the main conclusions are identical if a left-handed system is chosen instead. The easiest way to derive the equation for the volume of a unit cell is through the triple product of its axes. This, in turn, is equivalent to the determinant of the matrix whose rows are the cell axis components (see, for instance, Giacovazzo et al., 2002[Giacovazzo, C., Monaco, H. L., Artioli, G., Viterbo, D., Ferraris, G., Gilli, G., Zanotti, G. & Catti, M. (2002). Fundamentals of Crystallography, 2nd ed. Oxford University Press.], pp. 68–69, or Woolfson, 1997[Woolfson, M. M. (1997). An Introduction to X-ray Crystallography, 2nd ed. Cambridge University Press.], pp. 61–62),

[V = {\bf a} \cdot [{\bf b} \times {\bf c}] = \det (A) \equiv \Biggl| \matrix { a_{x} & a_{y} & a_{z} \cr b_{x} & b_{y} & b_{z} \cr c_{x} & c_{y} & c_{z}} \Biggr|.]

When cell lengths and cell angles are preferred to Cartesian components, it is customary to square the volume. As the determinant of a transpose matrix is equivalent to the determinant of the matrix itself, we have

[\eqalign{V^2 &= [\det (A)]^2 = \det (A) \det (A^t) = \det (AA^t) \cr &= \Biggl| \matrix { {\bf a} \cdot {\bf a} & {\bf a} \cdot {\bf b} & {\bf a} \cdot {\bf c} \cr {\bf b} \cdot {\bf a} & {\bf b} \cdot {\bf b} & {\bf b} \cdot {\bf c} \cr {\bf c} \cdot {\bf a} & {\bf c} \cdot {\bf b} & {\bf c} \cdot {\bf c}}\Biggr|\cr}]

or

[V^2 = a^2b^2c^2 \Biggl| \matrix { 1 & \cos\gamma & \cos\beta \cr \cos\gamma & 1 & \cos\alpha \cr \cos\beta & \cos\alpha & 1} \Biggr| \equiv a^2b^2c^2P .\eqno (1)]

Equation (1)[link] for the square of the volume is interesting as it makes it immediately evident that cell lengths are independent of cell angles. From the same expression it also appears that determinant P is of key importance in the calculation of the unit-cell volume. In the mathematical context of trihedral angles, its square root has been known as the generalized sine of the trihedral angle (Allendoerfer, 1965[Allendoerfer, C. B. (1965). Math. Mag. 38, 253-259.]; Nahir, 1998[Nahir, Y. (1998). Int. J. Math. Educ. Sci. Technol. 29, 65-82.]). It is supposed to be a positive quantity, otherwise the volume will be zero or an imaginary number. Unfortunately, P could assume null or negative values for certain combinations of angles. In fact, P can assume any real value between −4 and 1. Therefore, we will have to find the conditions on [\alpha], [\beta] and [\gamma] that make P positive. The analytic expression for P is

[P = 1 - \cos ^2 \alpha - \cos ^2 \beta - \cos ^2 \gamma + 2 \cos\alpha\cos\beta\cos\gamma .\eqno (2)]

Finding the conditions that make equation (2)[link] positive can be quite complicated. It is preferable to use an equivalent factorized form (Allendoerfer, 1965[Allendoerfer, C. B. (1965). Math. Mag. 38, 253-259.]; Sivardière, 1996[Sivardière, J. (1996). Acta Cryst. A52, 955-956.]),

[\eqalignno{P &= 4\sin \left ({{\alpha+\beta+\gamma} \over {2}}\right) \sin \left ({{\alpha+\beta-\gamma} \over {2}}\right) \sin \left ({{\alpha+\gamma-\beta} \over {2}}\right)&\cr &\quad\times\sin \left ({{\beta+\gamma-\alpha} \over {2}}\right). &(3)\cr}]

This equation has also been reported in International Tables for Crystallography (Donnay & Donnay, 1959[Donnay, J. D. H. & Donnay, G. (1959). International Tables for X-ray Crystallography, Vol. II. Birmingham: Kynoch Press.]; Koch, 2004[Koch, E. (2004). International Tables for Crystallography, Vol. C. Dordrecht: Kluwer Academic Publishers.]), and it is frequently used in computational programs for crystallography. A demonstration of how equation (3)[link] can be derived from equation (2)[link] is shown in Appendix A[link].

3. Determination of the allowed range for the unit-cell angles

It is straightforward to find out when P is positive if equation (3)[link] is used. Simply, the product of its factors will have to be positive. This can happen if all factors have a positive sign, or if all have a negative sign, or if two factors are positive and two are negative. We can limit [\alpha], [\beta] and [\gamma] to values in the range [0,180], as numbers higher than these correspond to a specular image of the same unit cell, better described with a left-handed coordinate system. We will, therefore, be looking to a region where P is positive, which is contained in the cube with a vertex at (0,0,0) and with side equal to 180. Each sine factor in the expression will be positive or negative on either of the sides of a specific plane. For example, [\sin[(\alpha+\beta+\gamma)/2]] will be positive above the plane of equation

[\alpha + \beta + \gamma = 0]

and below the plane of equation

[\alpha + \beta + \gamma - 360 = 0.]

This intersects the cube in a closed, convex region; in this part of the cube all combinations of the three angles will give a positive value for the sine, in the complementary region of the cube to this the sine will be negative. By considering all possible combinations of signs for the sines and by inspecting, for each combination, which region of the cube is positive and which is negative, we are able to determine a convex region within which P is positive. This is the interior of the tetrahedron depicted in Fig. 1[link], whose vertices have coord­inates [A \equiv (0,0,0)], [B \equiv (180,0,180)], [C \equiv (0,180,180)] and [D \equiv (180,180,0)].

[Figure 1]
Figure 1
Each point inside this tetrahedron has coordinates [(\alpha,\beta,\gamma)] which constitute allowed values for the unit cell of a triclinic system. The extra lines and points depicted inside the tetrahedron correspond to the angles of crystal systems different from the triclinic (see text).

As the volume of the tetrahedron is one third of the volume of the cube, there are, in fact, going to be two times more combinations of ([\alpha,\beta,\gamma]) giving a null or imaginary volume than those giving a positive, admissible volume. From a purely quantitative point of view, the three angles [\alpha], [\beta] and [\gamma] will be valid unit-cell angles if they satisfy at the same time the following four inequalities:

[\left \{ \matrix { 0 & \lt & \alpha + \beta + \gamma & \lt & 360, \cr 0 & \lt & \alpha + \beta - \gamma & \lt & 360, \cr 0 & \lt & \alpha - \beta + \gamma & \lt & 360, \cr 0 & \lt & - \alpha + \beta + \gamma & \lt & 360.}\right . \eqno (4)]

For example, a unit cell with angles 60, 60 and 130 cannot exist, despite the fact that they numerically appear perfectly sensible values for a triclinic cell. Fig. 1[link] also shows angle combinations corresponding to those systems different from the triclinic. Point [K \equiv (90,90,90)] corresponds to the orthorhombic, tetragonal and cubic systems, point [H \equiv (90,90,120)] to the hexagonal one. The straight line going from [L \equiv (90,0,90)] to [L^{\prime} \equiv (90,180,90)] corresponds to the monoclinic system, where the [\beta] angle can assume any value in the range (0,180). Finally, the straight line going from [A \equiv (0,0,0)] to [A^{\prime} \equiv (120,120,120)] corresponds to the rhombohedral system, where all three cell angles are equal to a single value between 0 and 120. It is immediately obvious that all these points lie inside the tetrahedron, as it should be.

4. A geometrically intuitive explanation for null and imaginary volumes

In the previous section, quantitative arguments for unit-cell angle constraints have been presented. It would be desirable to find a qualitative or intuitive explanation of why there are constraints in the first place. The following simple geometric construction should provide the reader with such an explanation (Fig. 2[link]).

[Figure 2]
Figure 2
In this geometric construction, two of the three unit-cell axes, [{\bf a}] and [{\bf b}], are the central axes for two cones of half-angles [\beta] and [\alpha], respectively. When the cones intersect in two straight lines, then two specularly equivalent unit cells are constructed, and the third unit-cell axis, [{\bf c}], coincides with any of the two intersections. If the cones' apertures are progressively made to shrink, the two intersections will get closer and closer, until they coincide with a single straight line lying on the ab plane; at this point the unit cell will have collapsed to a plane, and its volume is null. Further shrinking will have the effect of making even that single intersection disappear; in this case we have no unit cell and, consequently, its volume is an imaginary number.

Let [\gamma] be the angle between axes [{\bf a}] and [{\bf b}] of the cell. The third axis, [{\bf c}], will form an angle [\alpha] with the [{\bf b}] axis and an angle [\beta] with the [{\bf a}] axis. This means that [{\bf c}] must lie at the same time on the surface of two cones with vertices at the origin of the unit cell. The first has its axis coincident with [{\bf a}] and half-angle equal to [\beta], while the second has its axis coincident with [{\bf b}] and half-angle equal to [\alpha]. For certain combinations of the three angles the two cones will intersect in two straight lines starting from the origin and diverging specularly from the plane formed by [{\bf a}] and [{\bf b}]; in this case [{\bf c}] coincides with any of the two straight lines, and two specular unit cells will be constructed. When the cones' half-angles shrink beyond a certain value, though, the two straight lines will converge towards a single line on the ab plane, and further shrinking will result in no intersection at all. In the first case the unit-cell volume will become zero, and the cell will degenerate into a parallelogram; in the second case the volume will be imaginary and no unit cell can be built.

5. Discussion

A unit cell describing a crystallographic structure must, quite obviously, have a finite, not imaginary, volume. In this sense, Nature is never going to be fooled by wrong combinations of angles, as they only occur for null or imaginary volumes. Constraints [equation (4)[link]] are only useful in the context of a mathematical theory or a computer algorithm. The authors were led to the investigations underlying this paper when they tested their cRy package for statistical applications in crystallography (Waterman et al., 2011[Waterman, D., Evans, G. & Foadi, J. (2011). In preparation.]) with unit-cell angles 60, 20, 80 (corresponding to a null volume) and found inconsistent, infinite results. Their software now includes conditions [equation (4)[link]] to check against wrongly designed cells.

APPENDIX A

A factorized expression for the determinant P

Using the following trigonometric identity

[\cos^2 \alpha = [1+\cos(2\alpha)]/2, \eqno (5)]

let us replace all squared cosines in equation (2)[link]. This yields

[P = - 1/2 - [\cos(2\alpha) + \cos(2\beta) + \cos(2\gamma)]/2 + 2\cos\alpha\cos\beta\cos\gamma. \eqno (6)]

With another trigonometric identity, turning a sum into a product,

[\cos\alpha + \cos\beta = 2\cos [(\alpha + \beta)/2] \cos [(\alpha-\beta)/2],]

equation (6)[link] is readily transformed into

[P = -[1 - \cos (2\gamma)]/2 - \cos (\alpha + \beta) \cos (\alpha - \beta) + 2\cos\alpha\cos\beta\cos\gamma]

and, using again the identity [equation (5)[link]],

[\eqalign{P &= - \cos^2\gamma - \cos (\alpha + \beta)\cos (\alpha-\beta) + \cos(\alpha + \beta)\cos\gamma \cr &\quad+ \cos(\alpha-\beta)\cos\gamma.\cr}]

This last expression can easily be factorized as follows

[P = [\cos\gamma - \cos(\alpha-\beta)][\cos(\alpha + \beta) - \cos\gamma]. \eqno (7)]

Equation (7)[link] can be turned into the final factorized formula [equation (3)[link]] once the trigonometric identity

[\cos\alpha - \cos\beta = - 2\sin[(\alpha + \beta)/2] \sin [(\alpha-\beta)/2]]

is applied to each of its factors.

Acknowledgements

The authors would like to thank one of the referees for the suggestion to stress the right-handedness of the system of axis chosen in the text. This work was supported by grant No. WT089809MA from the Wellcome Trust.

References

First citationAllendoerfer, C. B. (1965). Math. Mag. 38, 253–259.  CrossRef Google Scholar
First citationDonnay, J. D. H. & Donnay, G. (1959). International Tables for X-ray Crystallography, Vol. II. Birmingham: Kynoch Press.  Google Scholar
First citationGiacovazzo, C., Monaco, H. L., Artioli, G., Viterbo, D., Ferraris, G., Gilli, G., Zanotti, G. & Catti, M. (2002). Fundamentals of Crystallography, 2nd ed. Oxford University Press.  Google Scholar
First citationKoch, E. (2004). International Tables for Crystallography, Vol. C. Dordrecht: Kluwer Academic Publishers.  Google Scholar
First citationNahir, Y. (1998). Int. J. Math. Educ. Sci. Technol. 29, 65–82.  CrossRef Google Scholar
First citationSivardière, J. (1996). Acta Cryst. A52, 955–956.  CrossRef Web of Science IUCr Journals Google Scholar
First citationWaterman, D., Evans, G. & Foadi, J. (2011). In preparation.  Google Scholar
First citationWoolfson, M. M. (1997). An Introduction to X-ray Crystallography, 2nd ed. Cambridge University Press.  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.

Journal logoFOUNDATIONS
ADVANCES
ISSN: 2053-2733
Follow Acta Cryst. A
Sign up for e-alerts
Follow Acta Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds