# User:Ddcampayo/Birefrigence

## Theory

More generally, birefringence can be defined by considering a dielectric permittivity and a refractive index that are tensors. Consider a plane wave propagating in an anisotropic medium, with a relative permittivity tensor ε, where the refractive index n, is defined by ${\displaystyle n\cdot n=\epsilon }$. If the wave has an electric vector of the form:

${\displaystyle \mathbf {E=E_{0}} \exp \left[i(\mathbf {k\cdot r} -\omega t)\right]\,}$ (2)

where r is the position vector and t is time, then the wave vector k and the angular frequency ω must satisfy Maxwell's equations in the medium, leading to the equations:

${\displaystyle -\nabla \times \nabla \times \mathbf {E} ={\frac {1}{c^{2}}}(\mathbf {\epsilon } \cdot {\frac {\partial ^{2}\mathbf {E} }{\partial t^{2}}})}$ (3a)

${\displaystyle \nabla \cdot (\mathbf {\epsilon } \cdot \mathbf {E} )=0}$ (3b)

where c is the speed of light in a vacuum. Substituting eqn. 2 in eqns. 3a-b leads to the conditions:

${\displaystyle |\mathbf {k} |^{2}\mathbf {E_{0}} -\mathbf {(k\cdot E_{0})k} ={\frac {\omega ^{2}}{c^{2}}}(\mathbf {\epsilon } \cdot \mathbf {E_{0}} )}$ (4a)

${\displaystyle \mathbf {k} \cdot (\mathbf {\epsilon } \cdot \mathbf {E_{0}} )=0}$ (4b)

For the matrix product ${\displaystyle (\epsilon \cdot \mathbf {E} )}$ often a separate name is used, the dielectric displacement vector ${\displaystyle \mathbf {D} }$. So essentially birefringence concerns the general theory of linear relationships between these two vectors in anisotropic media.

To find the allowed values of k, E0 can be eliminated from eq 4a. One way to do this is to write eqn 4a in Cartesian coordinates, where the x, y and z axes are chosen in the directions of the eigenvectors of ε, so that

${\displaystyle \mathbf {\epsilon } ={\begin{bmatrix}n_{x}^{2}&0&0\\0&n_{y}^{2}&0\\0&0&n_{z}^{2}\end{bmatrix}}\,}$ (4c)

Hence eqn 4a becomes

${\displaystyle (-k_{y}^{2}-k_{z}^{2}+{\frac {\omega ^{2}n_{x}^{2}}{c^{2}}})E_{x}+k_{x}k_{y}E_{y}+k_{x}k_{z}E_{z}=0}$ (5a)

${\displaystyle k_{x}k_{y}E_{x}+(-k_{x}^{2}-k_{z}^{2}+{\frac {\omega ^{2}n_{y}^{2}}{c^{2}}})E_{y}+k_{y}k_{z}E_{z}=0}$ (5b)

${\displaystyle k_{x}k_{z}E_{x}+k_{y}k_{z}E_{y}+(-k_{x}^{2}-k_{y}^{2}+{\frac {\omega ^{2}n_{z}^{2}}{c^{2}}})E_{z}=0}$ (5c)

where Ex, Ey, Ez, kx, ky and kz are the components of E0 and k. This is a set of linear equations in Ex, Ey, Ez, and they have a non-trivial solution if their determinant is zero:

${\displaystyle \det {\begin{bmatrix}(-k_{y}^{2}-k_{z}^{2}+{\frac {\omega ^{2}n_{x}^{2}}{c^{2}}})&k_{x}k_{y}&k_{x}k_{z}\\k_{x}k_{y}&(-k_{x}^{2}-k_{z}^{2}+{\frac {\omega ^{2}n_{y}^{2}}{c^{2}}})&k_{y}k_{z}\\k_{x}k_{z}&k_{y}k_{z}&(-k_{x}^{2}-k_{y}^{2}+{\frac {\omega ^{2}n_{z}^{2}}{c^{2}}})\end{bmatrix}}=0\,}$ (6)

Multiplying out eqn (6), and rearranging the terms, we obtain

${\displaystyle {\frac {\omega ^{4}}{c^{4}}}-{\frac {\omega ^{2}}{c^{2}}}\left({\frac {k_{x}^{2}+k_{y}^{2}}{n_{z}^{2}}}+{\frac {k_{x}^{2}+k_{z}^{2}}{n_{y}^{2}}}+{\frac {k_{y}^{2}+k_{z}^{2}}{n_{x}^{2}}}\right)+\left({\frac {k_{x}^{2}}{n_{y}^{2}n_{z}^{2}}}+{\frac {k_{y}^{2}}{n_{x}^{2}n_{z}^{2}}}+{\frac {k_{z}^{2}}{n_{x}^{2}n_{y}^{2}}}\right)(k_{x}^{2}+k_{y}^{2}+k_{z}^{2})=0\,}$ (7)

In the case of a uniaxial material, where nx=ny=no and nz=ne say, eqn 7 can be factorised into

${\displaystyle \left({\frac {k_{x}^{2}}{n_{o}^{2}}}+{\frac {k_{y}^{2}}{n_{o}^{2}}}+{\frac {k_{z}^{2}}{n_{o}^{2}}}-{\frac {\omega ^{2}}{c^{2}}}\right)\left({\frac {k_{x}^{2}}{n_{e}^{2}}}+{\frac {k_{y}^{2}}{n_{e}^{2}}}+{\frac {k_{z}^{2}}{n_{o}^{2}}}-{\frac {\omega ^{2}}{c^{2}}}\right)=0\,.}$ (8)

Each of the factors in eqn 8 defines a surface in the space of vectors k — the surface of wave normals. The first factor defines a sphere and the second defines an ellipsoid. Therefore, for each direction of the wave normal, two wavevectors k are allowed. Values of k on the sphere correspond to the ordinary rays while values on the ellipsoid correspond to the extraordinary rays.

For a biaxial material, eqn (7) cannot be factorized in the same way, and describes a more complicated pair of wave-normal surfaces.[1]

Birefringence is often measured for rays propagating along one of the optical axes (or measured in a two-dimensional material). In this case, n has two eigenvalues that can be labeled n1 and n2. n can be diagonalized by:

${\displaystyle \mathbf {n} =\mathbf {R(\chi )} \cdot {\begin{bmatrix}n_{1}&0\\0&n_{2}\end{bmatrix}}\cdot \mathbf {R(\chi )} ^{\textrm {T}}}$ (9)

where R(χ) is the rotation matrix through an angle χ. Rather than specifying the complete tensor n, we may now simply specify the magnitude of the birefringence Δn, and extinction angle χ, where Δn = n1 − n2.

A crystal of rutile exhibiting birefringence for rotating polarizer
1. ^ Born M, and Wolf E, Principles of Optics, 7th Ed. 1999 (Cambridge University Press), §15.3.3