52 Plane Waves in Anisotropic Media
Note: There is no text section in Jackson on this important topic. The student is assumed to figure it out by working through problems 7.16 and 7.17, and also analyzing the ionosphere as an anisotropic medium. The presentation below follows Landau, Continuous Media, p 333. There is a similar treatment in Zangwill.
We’ve already discussed one anisotropic medium, the ionosphere. A common anisotropic crystal found in elementary labs is calcite, CaCO3, whose strange optical properties are readily observed by putting a crystal on top of anythingthere is a double image! More spectacularly, a laser beam entering a calcite crystal will be split into two, and two parallel beams come out. In this picture, from Wikipedia, the blue laser light is only visible outside the crystal, from dust particle scattering. (The orange rays are fluoresence inside the crystal.)
Wave Vector and Energy Flow
Suppose we have a monochromatic plane wave in a transparent linear medium having the magnetic properties of the vacuum, but with a non-isotropic permittivity (a 3x3 tensor), that is:
This means that the vector is not, in general, parallel to Physically, it means the molecules polarize more easily in response to electric fields in particular directions.
Taking all quantities to be proportional to Maxwell's equations
for this plane wave are:
Staring at these equations, from the first is perpendicular to and and from the second is perpendicular to
That is, are all mutually perpendicular: see figure.
Next, is perpendicular to from the second equation, so must be in the plane.
This means that since is not parallel to is not perpendicular to .
It therefore follows that the direction of energy flow, , is not in the same direction as the wave vector !
The Wave Normal Vector
Following Landau, we define the dimensionless "wave normal" vector by
In vacuum, where would be a unit vector, and in an isotropic material, its length would be In an anisotropic medium, its length will vary with direction, as we’ll soon see.
The fundamental formulae now become (switching from to )
(Notice that no longer appears in the equations.)
From these equations: the energy flux Poynting vector can be written in terms of :
Exercise: Check that for an isotropic material, this gives the correct answer.
The displacement vector can be similarly expressed, using the fundamental formulae above:
Now combining this equation with we have:
This is a set of three homogeneous simultaneous equations for the components , of the form . The equations will be self-consistent provided that is,
In the coordinates where the dielectric tensor is diagonal, with elements it is straightforward to write this determinantal equation (in terms of )
For given anisotropic dielectric properties, this tells us what are allowed values of the wave normal vector. These values lie on a surface in the three-dimensional space.
Exercise: Check that this makes sense in the isotropic case, verify the expansion given just below.
Fresnel’s Equation and the Wave Vector Surface
Expanding out the determinant, the condition the equations be consistent becomes:
This is Fresnel's equation. In the space spanned by it describes a fourth-order surface, termed the wave vector surface. (There are sixth-order terms in the above expression, but they cancel out.)
Going back momentarily from to and writing this wave vector surface as (remembering the may vary with as well) makes clear it is the surface of all the vector points corresponding to waves having frequency
Group Velocity and Ray Vector
Suppose now a ray of light enters one of these crystals. How does its direction of propagation in the crystal relate to ? (Recall and the wave fields go as )
The direction of the ray is not because the energy, hence the physical ray, moves at the wave group velocity, and since the medium is anisotropic, is not in general parallel to
The wave vector surface is from which
so the ray direction is that of the normal to the wave vector surface at the given (You can check the sign by taking in 1D.) Landau calls the physical ray direction the ray vector, written and normalizes it by requiring
Two Refractive Indices
If we take a fixed direction, Fresnel's equation becomes a quadratic in meaning there are two physical refractive indices, two different wave speeds. How is that possible? It's because different polarizations propagate at different speeds: the electric oscillating field gets different response from the medium depending on the direction it's tugging the molecules.
It's all a lot simpler for uniaxial crystals, where (diagonalizing the dielectric matrix) Standard notation, in terms of the direction of the axis of symmetry of the crystal, is to write these diagonal elements as
In this case, the Fresnel equation factors to
This means the wave vector surface becomes two separate surfaces, a sphere and an ellipsoid:
(Shown here is a negative uniaxial crystal, a positive one also has the ellipsoid touching the sphere at the poles, but the ellipsoid now contains the sphere.)
The wave corresponding to has isotropic index of refraction, so is always parallel to This is just an ordinary ray: for a given vector pointing to a spot on the spherical surface, the normal to that surface is of course in the same direction, so the group velocity direction coincides with the ray direction.
On the other hand, for the allowed vector on the ellipsoidal surface, the normal to the surface will not be in the direction (except for the axes end points), and the group velocity, the direction of energy transfer, will not be : this is termed the extraordinary ray.
For this uniaxial crystal, the extraordinary ray will be in the plane of the optical axis and the ordinary ray (see figure). (Landau calls this plane the "principal section".) Differentiating the equation for the ellipse,
along the ellipse, so the direction of the normal is and the direction of the extraordinary ray is therefore given by
When a ray of light enters such a crystal, it separates into these two beams, ordinary and extraordinary, having different velocities, and orthogonal polarizations.
Simple Classical Model
To construct a simple classical model having this kind of behavior, assume that the simple harmonic oscillator restoring force potential in the Lorentz molecule model is replaced by not unreasonable for modeling the bonds between atoms holding the electrons in place in a non-isotropic lattice having an axis of symmetry.
Now a signal going through the crystal in the direction of the optic axis will have its electric field vector in the plane and will have the same speed irrespective of its polarizationall directions in that plane are equivalent. But what about a ray going in along the axis? If it's polarized along the axis, will govern the oscillators’ response. But for polarization along the axis, is what counts. So the two polarizations will travel at different speeds.
This is how a quarter-wave plate works: if linearly polarized light enters, the electric field in the plane at 45 degrees to the axis, and the speed difference and plate thickness is such that one component of the polarized ray is delayed a quarter cycle relative to the other on passing through the plate, circularly polarized light will emerge.