- 1. Some history: the idea behind the polar lattice
- 2. Construction of the polar lattice
- 3. From the polar lattice to the reciprocal lattice
- 4. The polar lattice as a special case of the pseudo-reciprocal lattice
- 5. The antireciprocal lattice
- 6. Deriving dual lattices from dual space
- 7. Conclusions
- References

- 1. Some history: the idea behind the polar lattice
- 2. Construction of the polar lattice
- 3. From the polar lattice to the reciprocal lattice
- 4. The polar lattice as a special case of the pseudo-reciprocal lattice
- 5. The antireciprocal lattice
- 6. Deriving dual lattices from dual space
- 7. Conclusions
- References

## teaching and education

## The Bravais

as a didactic tool for diffraction beginners^{a}Cristallographie, Résonance Magnétique et Modélizations (CRM2), UMR–CNRS 7036, Faculté des Sciences et Technologies, Institut Jean Barriol, Nancy Université, Boulevard des Aiguillettes, BP 70239, F54506 Vandœuvre-lès-Nancy Cedex, France, and ^{b}Institute for Mathematics, Astrophysics and Particle Physics, Faculty of Science, Mathematics and Computing Science, Radboud University Nijmegen, Postbus 9010, 6500 GL Nijmegen, The Netherlands^{*}Correspondence e-mail: massimo.nespolo@crm2.uhp-nancy.fr

When undergraduate students discover crystallography for the first time, they are usually already familiar with the phenomenon of diffraction as the `bending' of waves around small obstacles. The special (periodic) nature of crystals acting as `diffraction gratings' that produce interference of diffracted waves is typically rationalized in terms of the

of the crystal. The concept of the however, remains somewhat abstract for beginners, until they perform a diffraction experiment. It can be made more easily understandable through an intermediate step, namely its ancestor, the Bravais By means of a short historical trip through pre-X-ray crystallography, a generalized introduction to the notion of the dual lattice is given, of which the is the most common but by no means the only example, and it is shown how the use of the Bravais can ease the introduction of the reciprocal lattice.Keywords: crystallographic education; Bravais polar lattice; crystallographic history; dual lattice.

### 1. Some history: the idea behind the polar lattice

The atomic positions in a crystalline solid are usually described by fractional coordinates inside the ), who obtained a number of 15. It was Bravais (1850; English translation by Shaler, 1949) who found the correct number of 14, and the three-dimensional lattices are now universally known as Bravais lattices.] The points of the are those points to which a fixed corner of the can be moved by a translational symmetry of the crystal. The Bravais lattices express the periodicity in the space that physically contains the crystals and where we observe them. This space is usually called in crystallography, or, in particular when decorated with atoms, The concept of a in is very well established and normally its introduction in undergraduate courses does not present difficulties.

and the periodicity of the atomic distribution is expressed through the [The derivation of the independent types of lattices in three-dimensional space was first obtained by Frankenheim (1842Until the beginning of the 20th century, crystals were studied morphologically, optically, and by physical and chemical tests, but there were no direct means to determine their atomic structure. With the discovery of X-rays, and the pioneer experiments on their diffraction by crystals (Friedrich *et al.*, 1912), scientists working on condensed matter had a new world accessible to them. To interpret the diffraction pattern produced when a crystal is subject to an X-ray beam, the concept of the was introduced (Ewald, 1913, 1921), and the lattice in was then renamed the the two being naturally dual to each other.

Actually, Ewald did not conjure the concept of the ; reprinted by Yale University, 1947; https://www.archive.org/details/117714283 ) had introduced a *reciprocal system of vectors*, for which** a**_{i} · **a**_{j}* = δ_{ij}, *i.e.* precisely the same mathematical definition that is often used to introduce the However, well before Gibbs, Bravais (1850) had introduced what is the true ancestor of the a lattice dual to the which is nevertheless not reciprocal to it because its parameters are measured in ångströms, and not in ångströms^{−1}. It is called the and played an important role in the study of crystal morphology.

As is well known, the law of the constancy of interfacial angles (Stenon's law) states that interfacial angles between corresponding faces are constant. The relative development (size) of a face, and even its presence, depends on the growth conditions of the crystal and represents an accessory feature. The most evident consequence is that the habitus – the external shape of a crystal – may not reveal much about its symmetry. This is a common trap for students, because the eye is caught by the relative *sizes* of the faces, whereas it is their *orientation* that reflects the symmetry of the crystal. When studying wooden crystal models, for example, one must always be careful to stress that those models (usually) represent an idealized crystal, and that the relative sizes of the faces may be substantially different, although their relative orientation remains unchanged. In order to eliminate the effect of the accessory character of the crystal and to retain the essential features, the study of crystal morphology is performed by taking, from the centre of the crystal, the directions normal to each face, extended to intersect a sphere circumscribed around the crystal. These intersections, called spherical poles, are either projected onto a plane passing through the centre of the sphere, resulting in the or are extended to reach a plane tangential to the sphere, giving the (see, for example, Barker, 1922; Terpstra & Codd, 1961). The difference between these projections is the position of the point of sight (at the centre of the sphere or on the surface of it, respectively), while the common feature is precisely the elimination of the accessory character of the crystal – the relative size of the faces. The angle between the normals of two dihedral faces is taken as the angle between the faces, although it is actually the supplement of it: this comes from the way these angles are actually measured, namely by finding the orientation of the crystal, mounted on an optical goniometer, corresponding to the directions producing a reflection of the light incident on the faces.

A crystal face is simply the last plane of a family (*hkl*). Although the number of families is infinite, clearly only a limited number of faces can develop in the crystal morphology: the development of a face is in fact related to the growth conditions and the density of atoms (Bravais–Donnay–Harker law; Donnay & Harker, 1937). When the construction of a set of normals is performed on all the families of planes (*hkl*) and not only on the faces developed in the crystal habitus, the set of poles obtained when the normals are extended to meet the sphere represents the families of planes. However, the metric part of the information, namely the interplanar distance *d*_{(hkl)}, is lost. This metrical information was added by Bravais (1850), who transformed the set of *directions* normal to the families (*hkl*) to a system of *vectors*, each vector having a norm directly related to *d*_{(hkl)}. The periodic repetition of these vectors builds up a lattice that is dual to the original (direct) lattice, because the same construction applied twice returns to the original lattice. Because the intersections of the normals to the families (*hkl*) with the sphere circumscribed around the crystal are called poles, Bravais' construction was called the (`*réseau polaire*' in the French original).

### 2. Construction of the polar lattice

The volume of the three-dimensional **a**, **b** and **c** is

where **G** is the containing the dot products of the basis vectors. The norm |**a** × **b**| is the surface of the two-dimensional cell defined by the vectors **a** and **b**, which Bravais indicated as *S*(001). It corresponds to a vector normal to (001), the norm of which is |**a**||**b**|sinγ. The scalar product of **c** and a unit vector perpendicular to (001), *i.e.* **c** · (**a** × **b**)/|**a** × **b**|, gives the perpendicular distance between two planes of the (001) family, *i.e.* *d*_{(001)} (Fig. 1). The volume of the can therefore be expressed as *V* = *S*(001)*d*_{(001)}. This result is immediately generalized to any family (*hkl*) because all cells with the same multiplicity have the same volume. Therefore, the volume of the is *V* = *S*(*hkl*)*d*_{(hkl)}, provided that the multiplicity is kept unchanged when changing the family of planes.

Bravais took a coordinate system *O***abc**, defining a primitive of volume *V*, and drew from the origin *O* the normals to three conjugated planes of the lattice. On each of these normals he marked a series of equidistant points, *S*(*hkl*)/*V*^{1/3} apart. The factor *S*(*hkl*) makes the distance between these points proportional to the size of the two-dimensional cell of the (*hkl*) family of lattice planes. The cube root of *V* is a normalization factor corresponding to the mean distance between the nodes of a lattice, *i.e.* the side of a cube having the same volume as the primitive of the lattice. As will be demonstrated, this normalization is required to obtain a of the same volume as the of the original lattice. Bravais' construction produces a new axial setting, *O***a**^{p}**b**^{p}**c**^{p}, the basis vectors of which have parameters *S*(100)/*V*^{1/3}, *S*(010)/*V*^{1/3} and *S*(001)/*V*^{1/3} (Fig. 2).

The lattice built on *O***a**^{p}**b**^{p}**c**^{p} is the (indicated by the superscript p) and has, evidently, the same dimensions as the original lattice, namely ångströms. A row [*hkl*]^{p} of the corresponds to each family (*hkl*) of the original lattice. The period of [*hkl*]^{p} is the norm of the vector from the origin to the first node along this row, *i.e.* by construction

In particular, the cell parameters of the

arewhere **u**_{v} is a unit vector in the direction of **v**^{p}.

Let **B** be the 3 × 3 matrix with **a**, **b** and **c** as its columns, and let **B**^{p} be the matrix with **a**^{p}, **b**^{p} and **c**^{p} as its columns. By the definition of the basis of the we have

from which one sees that **B**^{T}**B**^{p} = *V*^{2/3}**I**_{3}, where **I**_{3} is the 3 × 3 identity matrix. This means that **B**^{p} = *V*^{2/3}**B**^{−T} and (**B**^{p})^{T} = *V*^{2/3}**B**^{−1}.

For the metric tensors, this gives **G** = **B**^{T}**B** and **G**^{p} = (**B**^{p})^{T}**B**^{p} = *V*^{2/3}**B**^{−1}*V*^{2/3}**B**^{−T} = V^{4/3}**G**^{−1}. In particular, one has det(**G**^{p}) = (*V*^{4/3})^{3}/det(**G**) = *V*^{4}/*V*^{2} = *V*^{2}, and thus *V*^{p} = *V*, *i.e.*

Bravais' normalization factor for the period along the lattice row [*hkl*]^{p} is thus chosen such that the of the has the same size as that of the original lattice.

The result that taking the **B**^{p} = *V*^{2/3}**B**^{−T}, since this implies (**B**^{p})^{p} = (*V*^{p})^{2/3}(**B**^{p})^{−T} = *V*^{2/3}(*V*^{2/3}**B**^{−T})^{−T} = *V*^{2/3}(*V*^{−2/3}**B**) = **B**, but it also follows directly from the properties of the vectors. By definition, **b**^{p} is perpendicular to **a** and **c**, and **c**^{p} is perpendicular to **a** and **b**, hence **a** is perpendicular to **b**^{p} and **c**^{p}. Since also (**a**^{p})^{p} is perpendicular to **b**^{p} and **c**^{p}, **a** and (**a**^{p})^{p} lie on the same line. Finally, since **a** · **a**^{p} = **a**^{p} · (**a**^{p})^{p} = *V*^{2/3}, **a** and (**a**^{p})^{p} actually have to be equal.

Determining the cell parameters is now straightforward. The norms of the basis vectors are

taking *V*^{p} = *V*. The interaxial angles are the supplement of the dihedral angles between corresponding planes,

and from equation (3) one obtains immediately

The *F* and *I* types, which are interchanged. Therefore, we find all the main features of the later but expressed in the same metric as the The has the merit of being easily understandable by any student who has handled a couple of wooden models and has realized that the orientation of the faces, not their size, is the external character reflecting the true symmetry of the crystal. The dihedral angle (actually, the supplement of it, as we have seen) is easily computed as the angle φ between two vectors of the

where the *bra* 〈**hkl**| simply denotes the row vector (*h*,*k*,*l*) and the *ket* |**hkl**〉 the column vector (*h*,*k*,*l*)^{T}. The *bra* and *ket* (together forming a *bracket*) thus provide a convenient notation for the dot product calculation **v**^{T} · **G** · **w** involving the **G**.

### 3. From the to the reciprocal lattice

The discovery of the wave nature of the X-rays through their diffraction by a crystal (Friedrich *et al.*, 1912) led Ewald (1913, 1921) to introduce the notion of the which rapidly replaced the to the extent that the latter was almost completely forgotten, with the exception of a few citations in rare texts like Terpstra & Codd (1961), Brasseur (1967) (where it was, however, erroneously mentioned simply as a synonym of the reciprocal lattice) and Rigault (1999), while it is not even mentioned in Volume B of *International Tables for Crystallography* (2008). Nowadays, it is common to introduce the in a rather axiomatic way, by imposing *k* = 1 in the analogue of equation (4) when **v**^{p} is replaced by **v*** (see *e.g.* Giacovazzo *et al.*, 2002). Actually, the is more easily understood through a more physical introduction. The amplitude of a diffracted wave is proportional to the Fourier transform of the electron density of the scattering body, but the scaling properties of the Fourier transform lead to a metric of the space that is reciprocal to the metric of the By repeating the construction of the but taking the period along the row [*hkl*]* as 1/*d*_{(hkl)} instead of *V*^{2/3}/*d*_{(hkl)} = *S*(*hkl*)/*V*^{1/3} (thus obtaining vectors measured in Å^{−1} instead of Å), all the well known properties of the are immediately obtained, in particular

Since **a*** = **a**^{p}/*V*^{2/3}, **b*** = **b**^{p}/*V*^{2/3} and **c*** = **c**^{p}/*V*^{2/3}, the interaxial angles for the are the same as those for the *i.e.* equation (8) remains unchanged.

Furthermore, for the **G*** we have **G*** = *V*^{−4/3}**G**^{p} = **G**^{−1}. From this, one sees that the dihedral angles for the are the same as for the since the scaling factor *V*^{−4/3} cancels out:

The introduction of the *V*^{2/3}.

### 4. The as a special case of the pseudo-reciprocal lattice

The polar and reciprocal lattices can actually both be seen as special cases of a generalized dual lattice, defined by taking a parametric value for *k* in equation (4) written for two dual bases (**v**_{1},…, **v**_{n}) and (**w**_{1},…, **w**_{n}). Grimmer (2003) introduced the term pseudo-reciprocal lattice for such a lattice, with the constraint, however, of having *k* of dimension length squared. If we remove this restriction, the pseudo-reciprocal lattice becomes the parent lattice of any lattice dual to the The is then clearly pseudo-reciprocal (with *k* = 1). Bravais' is pseudo-reciprocal with *k* = *V*^{2/3}, complying also with Grimmer's original definition.

A particularly pleasing feature of the pseudo-reciprocal lattice is that it can actually be defined without recourse to the basis. We demonstrate this for the *L* be a lattice with primitive lattice basis (**v**_{1},…, **v**_{n}). Because the chosen basis is primitive, the vectors in *L* are all linear combinations of **v**_{i} (*i* = 1,…, *n*) with integral coefficients. Let us then introduce a lattice as the set of all vectors **r** in having integral dot products with all the vectors in *L*, *i.e.*

If we now consider the *L** with basis (,…, \vf *v*_{n}^{*}) as defined above, it is clear that *L** , since are contained in (because all inner products are integral). Conversely, any vector **r** ∈ can be expressed as a linear combination of the vectors **v*** ∈ *L**: ** r** = with λ

_{i}∈ . But

**r**·

**v**

_{i}= λ

_{i}, and since

**r**∈ this gives λ

_{i}∈ . Thus,

**r**is an integral linear combination of the values of and hence contained in

*L**. For an arbitrary pseudo-reciprocal lattice with constant

*k*, the above argument is easily modified and shows that the pseudo-reciprocal lattice consists of all vectors

**r**in having dot products in with all vectors in

*L*.

### 5. The antireciprocal lattice

Bravais' choice to take *V*^{2/3}/*d*_{(hkl)} as the period along the rows of the dual lattice was quite evidently dictated by the desire that the length of the vector representing the (*hkl*) family of lattice planes should be proportional to the size of the two-dimensional cell in a plane of that family, and that the primitive of the new lattice should have the same volume as that of the A different natural choice is to take *d*_{(hkl)} as the period along the lattice row perpendicular to the (*hkl*) family. Such an approach comes out quite naturally from a morphological analysis and is guided by the fundamental features of the (*hkl*) families of lattice planes, namely the orientation of the planes, as expressed by their normals, and the interplanar distance. For the vector **r**^{a} representing the (*hkl*) family one immediately obtains

and this shows that some care must be taken with this approach. Since the vectors are unrestricted in their lengths, looking at all families (*hkl*) of lattice planes would produce vectors **r**^{a} of arbitrarily small length. Thus, running over all families (*hkl*), the corresponding vectors **r**^{a} do not form a lattice. In order to arrive at a lattice, we only construct the vectors **a**^{a}, **b**^{a} and **c**^{a} for a (**a**, **b**, **c**) and take (**a**^{a}, **b**^{a}, **c**^{a}) as the for a lattice that we call antireciprocal (indicated by the superscript a). Since the basis of the antireciprocal basis is defined exactly the same as that of the but with inverse parameters, we easily derive the analogous relations

While taking *d*_{(hkl)} as the period for the basis vectors of the antireciprocal lattice is very intuitive, the price to be paid for this is that the construction is basis-dependent. Equation (18) shows that the volume of the may change upon transition to a different basis of the since the interplanar distances may alter. Furthermore, the antireciprocal lattice is only pseudo-reciprocal if the interplanar distances *d*_{(100)}, *d*_{(010)} and *d*_{(001)} are all equal for the chosen basis (**a**, **b**, **c**).

#### 5.1. Example

Let (**e**_{1}, **e**_{2}, **e**_{3}) be a Cartesian (orthonormal) basis of (**e**_{i} · **e**_{j} = δ_{ij}) and *L* the standard lattice . For this basis, the polar, reciprocal and antireciprocal lattices coincide. However, in crystallography a conventional basis is usually chosen to comply with the symmetry directions of the lattice. Let such a basis be (**a**, **b**, **c**), with **a** = (100)^{T}, **b** = (110)^{T} and **c** = (001)^{T} with respect to (**e**_{1}, **e**_{2}, **e**_{3}). Then *V* = 1 and, for the reciprocal basis (**a***, **b***, **c***), one has **a*** = (10)^{T}, **b*** = (010)^{T} and **c*** = (001)^{T}, and *V** = 1/*V* = 1. For the antireciprocal basis (**a**^{a}, **b**^{a}, **c**^{a}), however, one obtains **a**^{a} = (10)^{T}, **b**^{a} = (010)^{T} and **c**^{a} = (001)^{T}, giving volume *V*^{a} = .

### 6. Deriving dual lattices from dual space

Underlying the concept of dual lattice is the much more abstract notion of dual space. Given a **V**, its dual space **V*** is defined to be the space of linear functions (functionals) from **V** to . Considering an element **v** of **V** as a column vector (*x*_{1},…, *x*_{n})^{T}, a function φ is no more than a 1 × *n* matrix (*u*_{1},…, *u*_{n}) (*i.e.* a row vector), and the application of φ to **v** is simply φ(**v**) = *u*_{1}*x*_{1} + ⋯ + *u*_{n}*x*_{n}. However, this is just the usual dot product of the column vector **u** = (*u*_{1},…, *u*_{n})^{T} with **v**. Thus, identifying row vectors with column vectors, **V*** is identified with **V** such that **u** ∈ **V** gives rise to the function φ_{u}: **v** **u** · **v**.

The vectors mapped to 0 by the function φ_{u} are precisely those vectors lying in the (*n* − 1)-dimensional subspace (called a hyperplane) of **V** that is perpendicular to **u**, which represents the kernel of φ_{u}. *Vice versa*, for every (*n* − 1)-dimensional subspace of **V** there is a function in **V*** (unique up to scalings) mapping this subspace to 0, namely φ_{u} such that **u** is perpendicular to the subspace. This makes it natural to define, for a given basis (**v**_{1},…, **v**_{n}) of **V**, a basis (φ_{1},…, φ_{n}) of **V*** such that φ_{i}(**v**_{j}) = δ_{ij}. With the identification of **V*** with **V**, the row vector φ_{i} read as a column vector is a vector , such that is perpendicular to **v**_{j} for *j* ≠ *i* and = 1. The basis (,…, ) is then called the of (**v**_{1},…, **v**_{n}).

The different types of dual lattices share the property that they have a λ_{1}**v**_{1}*, λ_{2}**v**_{2}*,… λ_{n}**v**_{n}*), *i.e.* such that the *i*th basis vector lies on the line perpendicular to all **v**_{j} with *j* ≠ *i*.

For the special case of **V** = , the vector product is a convenient means of finding the dual lattices, since for (**a**, **b**, **c**) the basis (**b** × **c**, **c** × **a**, **a** × **b**) fulfils the required orthogonality conditions and the vectors only have to be scaled as desired.

We note that the two dual spaces occur naturally in crystallography. The *E*^{n} is an affine space with underlying **V** = , *i.e.* every point of *E*^{n} is obtained as the translation of a single (but arbitrary) point of *E*^{n} by all vectors of **V**. The dual space **V*** of **V** (identified with **V**) is the in which the face normals and vectors reside.

### 7. Conclusions

The Bravais

is a natural intermediate between morphological and diffraction studies of crystals, both considering vectors that are normal to certain planes in The construction of the Bravais is intuitive and facilitates a good understanding of the before any diffraction experiment is performed.### References

Barker, T. V. (1922). *Graphical and Tabular Methods in Crystallography as the Foundation of a New System of Practice*. London: Thomas Murby and Co. Google Scholar

Brasseur, H. (1967). *Cristallographie*. Liège: F. Gothier. Google Scholar

Bravais, M. A. (1850). *J. École Polytech.* Cahier 33, **29**, 1–128. Google Scholar

Donnay, J. D. H. & Harker, D. (1937). *Am. Mineral.* **22**, 446–467. CAS Google Scholar

Ewald, P. P. (1913). *Phys. Z.* **14**, 465–472. CAS Google Scholar

Ewald, P. P. (1921). *Z. Kristallogr.* **56**, 129–156. Google Scholar

Frankenheim, M. L. (1842). *System der Krystalle*. Breslau: Von Grass, Barth und Comp. Google Scholar

Friedrich, W., Knipping, P. & von Laue, M. T. F. (1912). *Sitzungsber. Math. Phys. Kl. K. Bayer. Akad. Wiss. München*, pp. 303–322. Google Scholar

Giacovazzo, C., Monaco, H. L., Artioli, G., Viterbo, D., Ferraris, G., Gilli, G., Zanotti, G. & Catti, M. (2002). *Fundamentals of Crystallography*, 2nd ed. Chester, Oxford: IUCr/Oxford University Press. Google Scholar

Gibbs, J. W. (1881). *Vector Analysis*. New Haven: privately printed. Google Scholar

Grimmer, H. (2003). *Acta Cryst.* A**59**, 287–296. Web of Science CrossRef CAS IUCr Journals Google Scholar

*International Tables for Crystallography* (2008). Vol. B, 3rd ed. Heidelberg: Springer. Google Scholar

Rigault, G. (1999). *Atti Accad. Sci. Torino Cl. Sci. Fis. Mat. Nat.* **133**, 1–8. Google Scholar

Shaler, A. J. (1949). *Crystallogr. Soc. Am. Mem.* **1**. Google Scholar

Terpstra, P. & Codd, L. W. (1961). *Crystallometry.* London: Longmans. Google Scholar

© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.