Tables for
Volume D
Physical properties of crystals
Edited by A. Authier

International Tables for Crystallography (2013). Vol. D, ch. 1.1, pp. 14-16

Section Symmetry of tensors imposed by the crystalline medium

A. Authiera*

aInstitut de Minéralogie et de Physique des Milieux Condensés, 4 Place Jussieu, 75005 Paris, France
Correspondence e-mail: Symmetry of tensors imposed by the crystalline medium

| top | pdf |

Many papers have been devoted to the derivation of the invariant components of physical property tensors under the influence of the symmetry elements of the crystallographic point groups: see, for instance, Fumi (1951[link], 1952a[link],b[link],c[link], 1987[link]), Fumi & Ripamonti (1980a[link],b[link]), Nowick (1995)[link], Nye (1957[link], 1985[link]), Sands (1995)[link], Sirotin & Shaskol'skaya (1982)[link], and Wooster (1973)[link]. There are three main methods for this derivation: the matrix method (described in Section[link]), the direct inspection method (described in Section[link]) and the group-theoretical method (described in Section 1.2.4[link] and used in the accompanying software, see Section[link] ). Matrix method – application of Neumann's principle

| top | pdf |

An operation of symmetry turns back the crystalline edifice on itself; it allows the physical properties of the crystal and the tensors representing them to be invariant. An operation of symmetry is equivalent to a change of coordinate system. In a change of system, a tensor becomes [t'^{\alpha \beta}_{\gamma \delta} = t^{ij}_{kl}A^{\alpha}_{i} A^{\beta}_{j}B^{k}_{\gamma}B^{l}_{\delta}.]If A represents a symmetry operation, it is a unitary matrix: [A = B^{T} = B^{-1}.]

Since the tensor is invariant under the action of the symmetry operator A, one has, according to Neumann's principle, [t'^{\alpha \beta}_{\gamma \delta } = t^{\alpha \beta}_{\gamma \delta}]and, therefore, [t^{\alpha \beta}_{\gamma \delta} = t^{ij}_{kl} A^{\alpha}_{i} A^{\beta}_{j}B^{k}_{\gamma}B^{l}_{\delta}. \eqno(]

There are therefore a certain number of linear relations between the components of the tensor and the number of independent components is reduced. If there are p components and q relations between the components, there are [p-q] independent components. This number is independent of the system of axes. When applied to each of the 32 point groups, this reduction enables one to find the form of the tensor in each case. It depends on the rank of the tensor. In the present chapter, the reduction will be derived for tensors up to the fourth rank and for all crystallographic groups as well as for the isotropic groups. An orthonormal frame will be assumed in all cases, so that co- and contravariance will not be apparent and the positions of indices as subscripts or superscripts will not be meaningful. The [Ox_{3}] axis will be chosen parallel to the threefold, fourfold or sixfold axis in the trigonal, tetragonal and hexagonal systems. The accompanying software to the present volume enables the reduction for tensors of any rank to be derived. The operator A is in diagonal form

| top | pdf | Introduction

| top | pdf |

If one takes as the system of axes the eigenvectors of the operator A, the matrix is written in the form [\pmatrix { \exp {\rm i} \theta & 0 & 0\cr 0 & \exp -{\rm i} \theta & 0 \cr 0 & 0 & \pm 1 \cr },]where θ is the rotation angle, [Ox_{3}] is taken parallel to the rotation axis and coefficient [A_{3}] is equal to +1 or −1 depending on whether the rotation axis is direct or inverse (proper or improper operator).

The equations ([link] can then be simplified and reduce to [t^{ij}_{kl} = t^{ij}_{kl} A^i_{i}A^j_{j}B^k_{k}B^l_{l } \eqno{(}](without any summation).

If the product [A^{i}_{i}A\hskip1pt^{j}_{j}B^{k}_{k}B^{l}_{l}] (without summation) is equal to unity, equation ([link] is trivial and there is significance in the component [t_{kl}]. On the contrary, if it is different from 1, the only solution for ([link] is that [t^{ij}_{kl } = 0]. One then finds immediately that certain components of the tensor are zero and that others are unchanged. Case of a centre of symmetry

| top | pdf |

All the diagonal components are in this case equal to −1. One thus has:

  • (i) Tensors of even rank, [t\hskip1pt^{ij\ldots}= (-1)^{2p}t\hskip1pt^{ij\ldots}]. The components are not affected by the presence of the centre of symmetry. The reduction of tensors of even rank is therefore the same in a centred group and in its noncentred subgroups, that is in any of the 11 Laue classes: [\matrix{&\bar{{\bf 1}} \quad 1\hfill\cr &{\bf 2/{\bi m}} \quad 2,\ m \hfill\cr &{\bi mmm} \quad 222,\ 2mm \hfill\cr &\bar{{\bf 3}} \quad 3\hfill\cr &{\bar{\bf 3}}{\bi m} \quad 32, \ 3m\hfill\cr &{\bf 4/{\bi m}} \quad {\bar 4}, \ 4 \hfill\cr &{\bf 4/{\bi m m}} \quad {\bar 4}2m, \ 422, \ 4mm\hfill\cr &{\bf 6/{\bi m}} \quad {\bar 6}, \ 6 \hfill\cr &{\bf 6/{\bi m m}} \quad {\bar 6}2m, \ 622, \ 6mm\hfill\cr &{\bi m}{\bar{\bf 3}} \quad 23\hfill\cr &{\bi m}{\bar{\bf 3}}{\bi m} \quad 432, \ {\bar 4}32.\hfill\cr}]If a tensor is invariant with respect to two elements of symmetry, it is invariant with respect to their product. It is then sufficient to make the reduction for the generating elements of the group and (since this concerns a tensor of even rank) for the 11 Laue classes.

  • (ii) Tensors of odd rank, [t\hskip1pt^{ij\ldots}= (-1)^{2p+1}t\hskip1pt^{ij\ldots}]. All the components are equal to zero. The physical properties represented by tensors of rank 3, such as piezoelectricity, piezomagnetism, nonlinear optics, for instance, will therefore not be present in a centrosymmetric crystal. General case

| top | pdf |

By replacing the matrix coefficients [A^i_{i}] by their expression, ([link] becomes, for a proper rotation, [t\hskip1pt^{jk\ldots }= t\hskip1pt^{jk\ldots}\exp(ir\theta)\exp (-is\theta) (1)^{t}= t\hskip1pt^{jk\ldots} \exp i(r-s)\theta,]where r is the number of indices equal to 1, s is the number of indices equal to 2, t is the number of indices equal to 3 and [r + s + t = p] is the rank of the tensor. The component [t\hskip1pt^{jk\ldots}] is not affected by the symmetry operation if [(r - s)\theta = 2K\pi,]where K is an integer, and is equal to zero if [(r - s)\theta \neq 2K\pi.]

The angle of rotation θ can be put into the form [2\pi /q], where q is the order of the axis. The condition for the component not to be zero is then [r - s = Kq.]

The condition is fulfilled differently depending on the rank of the tensor, p, and the order of the axis, q. Indeed, we have [r - s \leq p] and

  • [p = 2], [r - s \leq 2]: the result of the reduction will be the same for any [q \geq 3];

  • [p = 3], [r - s \leq 3]: the result of the reduction will be the same for any [q \geq 4];

  • [p = 4], [r - s \leq 4]: the result of the reduction will be the same for any [q \geq 5].

It follows that:

  • (i) for tensors of rank 2, the reduction will be the same for trigonal (threefold axis), tetragonal (fourfold axis) and hexagonal (sixfold axis) groups;

  • (ii) for tensors of rank 3, the reduction will be the same for tetragonal and hexagonal groups;

  • (iii) for tensors of rank 4, the reduction will be different for trigonal, tetragonal and hexagonal groups.

The inconvenience of the diagonalization method is that the vectors and eigenvalues are, in general, complex, so in practice one uses another method. For instance, we may note that equation ([link] can be written in the case of [p = 2] by associating with the tensor a [3 \times 3] matrix T: [T = B T B^{T},]where B is the symmetry operation. Through identification of homologous coefficients in matrices T and [B T B^{T}], one obtains relations between components [t_{ij}] that enable the determination of the independent components. The method of direct inspection

| top | pdf |

The method of `direct inspection', due to Fumi (1952a[link],b[link], 1987[link]), is very simple. It is based on the fundamental properties of tensors; the components transform under a change of basis like a product of vector components (Section[link]).


  • (1) Let us consider a tensor of rank 3 invariant with respect to a twofold axis parallel to [Ox_{3}]. The matrix representing this operator is [\pmatrix{{\bar 1} &0 &0\cr 0 &{\bar 1} &0\cr 0 &0 &1\cr}.]

    The component [t_{ijk}] behaves under a change of axes like the product of the components [x_{i}, x_{j}, x_{k}]. The components [x_{1}, x_{2}, x_{3}] of a vector become, respectively, [-x_{1}], [ -x_{2}], [x_{3}]. To simplify the notation, we shall denote the components of the tensor simply by ijk. If, amongst the indices i, j and k, there is an even number (including the number zero) of indices that are equal to 3, the product [x_{i}x_{j}x_{k}] will become [-x_{i}x_{j}x_{k}] under the rotation. As the component `ijk' remains invariant and is also equal to its opposite, it must be zero. 14 components will thus be equal to zero: [111, 122, 133, 211, 222, 133, 112, 121, 212, 221, 323, 331, 332, 313.]

  • (2) Let us now consider that the same tensor of rank 3 is invariant with respect to a fourfold axis parallel to [Ox_{3}]. The matrix representing this operator and its action on a vector of coordinates [x_{1}, x_{2}, x_{3}] is given by [\pmatrix{\phantom{-}x_{2}\cr -x_{1}\cr \phantom{-}x_{3}\cr} = \pmatrix{0 &1 &0\cr {\bar 1} &0 &0\cr 0 &0 &1\cr} \pmatrix{x_{1}\cr x_{2}\cr x_{3}\cr}. \eqno(]Coordinate [x_{1}] becomes [x_{2}], [x_{2}] becomes [-x_{1}] and [x_{3}] becomes [x_{3}]. Component ijk transforms like product [x^{i}x\hskip1pt^{j}x^{k}] according to the rule given above. Since the twofold axis parallel to [Ox_{3}] is a subgroup of the fourfold axis, we can start from the corresponding reduction. We find [\matrix{311 &\Longleftrightarrow &\phantom{-}322 &: &t_{311} = &\phantom{-}t_{322}\phantom{.}\cr 123 &\Longleftrightarrow &-(213) &: &t_{123} = &-t_{213}\phantom{.}\cr 113 &\Longleftrightarrow &\phantom{-}223 &: &t_{113} = &\phantom{-}t_{223}\phantom{.}\cr 333 &\Longleftrightarrow &\phantom{-}333 &: &t_{333} = &\phantom{-}t_{333}\phantom{.}\cr 132 &\Longleftrightarrow &-(231) &: &t_{132} = &-t_{231}\phantom{.}\cr 131 &\Longleftrightarrow &\phantom{-}232 &: &t_{131} = &\phantom{-}t_{232}\phantom{.}\cr 312 &\Longleftrightarrow &-(321) &: &t_{312} = &-t_{321}.\cr}]All the other components are equal to zero.

It is not possible to apply the method of direct inspection for point group 3. One must in this case use the matrix method described in Section[link]; once this result is assumed, the method can be applied to all other point groups.


Fumi, F. G. (1951). Third-order elastic coefficients of crystals. Phys. Rev. 83, 1274–1275.
Fumi, F. G. (1952a). Physical properties of crystals: the direct inspection method. Acta Cryst. 5, 44–48.
Fumi, F. G. (1952b). The direct-inspection method in systems with a principal axis of symmetry. Acta Cryst. 5, 691–694.
Fumi, F. G. (1952c). Third-order elastic coefficients in trigonal and hexagonal crystals. Phys. Rev. 86, 561.
Fumi, F. G. (1987). Tables for the third-order elastic tensors in crystals. Acta Cryst. A43, 587–588.
Fumi, F. G. & Ripamonti, C. (1980a). Tensor properties and rotational symmetry of crystals. I. A new method for group 3(3z) and its application to general tensors up to rank 8. Acta Cryst. A36, 535–551.
Fumi, F. G. & Ripamonti, C. (1980b). Tensor properties and rotational symmetry of crystals. II. Groups with 1-, 2- and 4-fold principal symmetry and trigonal and hexagonal groups different from group 3. Acta Cryst. A36, 551–558.
Nowick, A. S. (1995). Crystal Properties via Group Theory. Cambridge University Press.
Nye, J. F. (1957). Physical Properties of Crystals, 1st ed. Oxford: Clarendon Press.
Nye, J. F. (1985). Physical Properties of Crystals, revised ed. Oxford University Press.
Sands, D. E. (1995). Vectors and Tensors in Crystallography. New York: Dover.
Sirotin, Y. I. & Shaskol'skaya, M. P. (1982). Fundamentals of Crystal Physics. Moscow: Mir.
Wooster, W. A. (1973). Tensors and Group Theory for the Physical Properties of Crystals. Oxford: Clarendon Press.

to end of page
to top of page