International
Tables for Crystallography Volume D Physical properties of crystals Edited by A. Authier © International Union of Crystallography 2013 
International Tables for Crystallography (2013). Vol. D, ch. 1.6, pp. 155157
Section 1.6.3. Linear optics^{a}Department of Physics, University of Oxford, Parks Roads, Oxford OX1 3PU, England, and ^{b}Department of Earth Sciences, University of Oxford, Parks Roads, Oxford OX1 3PR, England 
It is necessary, in order to understand fully the propagation of light through a general anisotropic crystal, to address the question of the way in which an electromagnetic wave is affected by its passage through a regular array of atoms or molecules. A full analysis of this problem at a microscopical level is complicated and was treated, for example, by Ewald (1916), who showed through consideration of a `halfcrystal' how to link the electromagnetic field outside the crystal to that inside (a good description of Ewald's work on this can be read in the book P. P. Ewald and his Dynamical Theory of Xray Diffraction, published by the International Union of Crystallography, Oxford Science Publications, 1992). For the purposes needed here, it is sufficient to apply Maxwell's equations to a bulk anisotropic continuum crystal, thus taking a macroscopic approach. The treatment here follows that given by Nussbaum & Phillips (1976).
Consider the relationship between the dielectric displacement and an electric field which in tensor terms is given by where is the vacuum dielectric permittivity and is a secondrank tensor, the relative dielectric tensor. Correspondingly, there is an induced polarization related to via where is another secondrank tensor, called the dielectric susceptibility tensor. Note that the restriction to a linear relationship between and (or and ) confines the theory to the region of linear optics. Addition of higherorder terms (see above) gives nonlinear optics. (Nonlinear optics is discussed in Chapter 1.7 .)where and are the magnetic induction and magnetic field intensity, respectively. It is customary at this point to assume that the crystal is nonmagnetic, so that , where is the vacuum magnetic permeability. If planewave solutions of the form are substituted into equations (1.6.3.3) and (1.6.3.4), the following results are obtained: These equations taken together imply that , and are vectors that are mutually orthogonal to one another: note that in general and need not be parallel. Similarly (and hence ), and are mutually orthogonal. Now, on substituting (1.6.3.9) into (1.6.3.8), Defining the propagation vector (or wave normal) by where is the unit vector in the direction of and n is the refractive index for light propagating in this direction, equation (1.6.3.10) then becomes Via the vector identity , this result can be transformed to is equal to and can be expressed simply in tensor form as . Now, with equation (1.6.3.1), the fundamental equation of linear crystal optics is found: where I is the unit matrix.
Equation (1.6.3.14) is the relevant starting point for the derivation of the way in which light propagates in an anisotropic medium. To solve it in a particular case, treat it as an eigenvector–eigenvalue problem: the are the eigenvectors and the eigenvalues. For example, take the case of a uniaxial crystal. The dielectric tensor is then given by Assume that light propagates along a direction in the plane, at an angle to the axis. Then, using (1.6.3.11), it is seen that and Substituting into equation (1.6.3.14) yields Solving this for the eigenvalues n gives as one solution and as the other. This latter solution can be rewritten as showing how the observed refractive index varies between the limits set by and , called the ordinary and extraordinary refractive index, respectively (sometimes these are denoted by and , respectively). Equation (1.6.3.21) can be thought of as the equation of a uniaxial ellipsoid (circular cross section) with the lengths of the semiaxes given by and . This is illustrated in Fig. 1.6.3.1, where is the direction of propagation of the light ray at an angle to . Perpendicular to , an elliptical cross section is cut from the uniaxial ellipsoid with semiaxes equal to and given by (1.6.3.21): the directions and also correspond to the eigenvectors of equation (1.6.3.4).
Direction is therefore the direction of the electric polarization transverse to the propagation direction, so that the refractive index measured for light polarized along is given by the value . For light polarized along , the refractive index would be given by equation (1.6.3.7). When is aligned along , a circular cross section of radius is obtained, indicating that for light travelling along and with any polarization the crystal would appear to be optically isotropic. The ellipsoid described here is commonly known as the optical indicatrix, in this case a uniaxial indicatrix.
Two cases are recognized (Fig. 1.6.3.2). When , the indicatrix is a prolate ellipsoid and is defined to be positive; when , it is oblate and defined to be negative. Note that when the indicatrix is a sphere, indicating that the refractive index is the same for light travelling in any direction, i.e. the crystal is optically isotropic. The quantity is called the linear birefringence (often simply called birefringence). In general, then, for light travelling in any direction through a uniaxial crystal, there will be two rays, the ordinary and the extraordinary, polarized perpendicular to each other and travelling with different velocities. This splitting of a ray of light into two rays in the crystal is also known as double refraction.
The origin of the birefringence in terms of the underlying crystal structure has been the subject of many investigations. It is obvious that birefringence is a form of optical anisotropy (the indicatrix is not spherical) and so it must be linked to anisotropy in the crystal structure. Perhaps the most famous early study of this link, which is still worth reading, is that of Bragg (1924), who showed that it was possible to calculate rough values for the refractive indices, and hence birefringence, of calcite and aragonite. His theory relied upon the summation of polarizability contributions from the Ca^{2+} and O^{2−} ions.
Returning now to the theory of the indicatrix, more general solution of the fundamental equation (1.6.3.14) leads to a triaxial ellipsoid, i.e. one in which all three semiaxes are different from one another (Fig. 1.6.3.3).
It is conventional to label the three axes according to the size of the refractive index by (or simply ). In such an ellipsoid, there are always two special directions lying in the γ–α plane, known as the optic axial plane, and perpendicular to which there are circular cross sections (shown shaded) of radius β. Thus these two directions are optic axes down which the crystal appears to be optically isotropic, with a measured refractive index β for light of any polarization. For this reason, crystals with this type of indicatrix are known as biaxial. When the angle between the optic axes, denoted conventionally as 2V, is acute about the γ axis, the crystal is positive biaxial, and when it is acute about α the crystal is negative biaxial. Note that as 2V becomes smaller, the biaxial indicatrix becomes closer to a uniaxial indicatrix (positive or negative). In all general directions the crystal is optically anisotropic. Thus, for light along , the measured refractive indices will be α and β for light polarized along and , respectively; for light along , α and γ are measured for light polarized along and , respectively; and along , β and γ are measured for light polarized along and , respectively. There are therefore three different linear birefringences to measure: γ–β, β–α and γ–α.
The different indicatrices are oriented in the crystal according to symmetry considerations (Table 1.6.3.1), and so their observation can form valuable and reliable indicators of the crystal system.

It has been seen how the refractive indices can be described in a crystal in terms of an ellipsoid, known as the indicatrix. Thus for orthogonal axes chosen to coincide with the ellipsoid axes, one can write where , and . One can write this equation alternatively as where the are the relative dielectric impermeabilities. For the indicatrix in any general orientation with respect to the coordinate axesThus the dielectric impermeability tensor is described by a secondrank tensor, related inversely to the dielectric tensor.
References
Bragg, W. L. (1924). The refractive indices of calcite and aragonite. Proc. R. Soc. London Ser. A, 105, 370.Ewald, P. P. (1916). Zur Begründung der Kristalloptik. Ann. Phys. (Leipzig), 49, 1–38, 117–143.
Nussbaum, A. & Phillips, R. A. (1976). Contemporary optics for scientists and engineers. New Jersey: Prentice Hall.