Transfer-matrix method (optics)

From Wikipedia, the free encyclopedia

For the transfer-matrix method in geometric optics, see Ray transfer matrix analysis. For the transfer-matrix method in statistical physics, see Transfer-matrix method.
Propagation of a ray through a layer
Propagation of a ray through a layer

The transfer-matrix method is a method used in optics and acoustics to analyze the propagation of electromagnetic or acoustic waves through a stratified (layered) medium.[1] This is for example relevant for the design of anti-reflective coatings and dielectric mirrors.

The reflection of light from a single interface between two media is described by the Fresnel equations. However, when there are multiple interfaces, such as in the figure, the reflections themselves are also partially reflected. Depending on the exact path length, these reflections can interfere destructively or constructively. The overall reflection of a layer structure is the sum of an infinite number of reflections, which is cumbersome to calculate.

The transfer-matrix method is based on the fact that, according to Maxwell's equations, there are simple continuity conditions for the electric field across boundaries from one medium to the next. If the field is known at the beginning of a layer, the field at the end of the layer can be derived from a simple matrix operation. A stack of layers can then be represented as a system matrix, which is the product of the individual layer matrices. The final step of the method involves converting the system matrix back into reflection and transmission coefficients.

Contents

[edit] Formalism for electromagnetic waves

Below is described how the transfer matrix is applied to electromagnetic waves (for example light) of a given frequency propagating through a stack of layers at normal incidence. It can be generalized to deal with incidence at an angle, absorbing media, and media with magnetic properties. We assume that the stack layers are normal to the z axis and that the field within one layer can be represented as the superposition of a left- and right-traveling wave with wave number k,

E(z) = Ereikz + Ele − ikz.

Because it follows from Maxwell's equation that E and F = dE / dz must be continuous across a boundary, it is convenient to represent the field as the vector (E(z),F(z)), where

F(z) = ikEreikz − ikEle − ikz.

Since there are two equations relating E and F to Er and El, these two representations are equivalent. In the new representation, propagation over a distance L into the positive z direction is described by the matrix

M = \left( \begin{array}{cc} \cos kL & \frac{1}{k} \sin kL \\ -k \sin kL & \cos kL \end{array} \right),

and

\left(\begin{array}{c} E(z+L) \\ F(z+L) \end{array} \right) = M\cdot \left(\begin{array}{c} E(z) \\ F(z) \end{array} \right)

Such a matrix can represent propagation through a layer if k is the wave number in the medium and L the thickness of the layer: For a system with N layers, each layer j has a transfer matrix Mj, where j increases towards higher z values. The system transfer matrix is then

M_s = M_N \cdot \ldots \cdot M_2 \cdot M_1.

Typically, one would like to know the reflectance and transmittance of the layer structure. If the layer stack starts at z = 0, then for negative z, the field is described as

E_L(z) = E_0 e^{ik_Lz} + r E_0 e^{-ik_Lz},

where E0 is the amplitude of the incoming wave, kL the wave number in the left medium, and r is the amplitude (not intensity!) reflectance coefficient of the layer structure. On the other side of the layer structure, the field consists of a right-propagating transmitted field

E_R(z) = t E_0 e^{ik_R z},

where t is the amplitude transmittance and kR is the wave number in the rightmost medium. If FL = dEL / dz and FR = dER / dz, then we can solve

\left(\begin{array}{c} E(z_R) \\ F(z_R) \end{array} \right) = M\cdot \left(\begin{array}{c} E(0) \\ F(0) \end{array} \right)

in terms of the matrix elements Mmn of the system matrix Ms and obtain

r = \frac{-i k_R M_{11} + k_L k_R M_{12} + M_{21} + ik_L M_{22}}{-M_{21} + ik_L M_{22} + i k_R M_{11} + k_L k_R M_{12}},

and

t = \frac{(-M_{21} + i k_L M_{22})(M_{11} + i k_L M_{12}) + (M_{11} - i k_L M_{12})(M_{21} + i k_L M_{22})}{-M_{21} + ik_L M_{22} + i k_R M_{11} + k_L k_R M_{12}}.

The intensity transmittance and reflectance, which are often of more practical use, are T = | t | 2 and R = | r | 2, respectively.

[edit] Example

As an illustration, consider a single layer of glass with a refractive index n and thickness d suspended in air at a wave number k (in air). In glass, the wave number is k' = nk. The transfer matrix is

M=\left(\begin{array}{cc}\cos k'd & \sin(k'd)/k' \\ -k' \sin k'd & \cos kd \end{array}\right).

The amplitude reflection coefficient can be simplified to

r = \frac{(1/n - n) \sin k'd}{(n+1/n)\sin k'd + 2 i \cos(k'd)}.

This configuration effectively describes a Fabry-Pérot interferometer or etalon: for k'd=0, \pi, 2\pi, \cdots, the reflection vanishes.

[edit] Acoustic waves

It is possible to apply the transfer-matrix method to sound waves. Instead of the electric field E and its derivative F, the displacement u and the stress σ = Cdu / dz, where C is Young's modulus should be used.

[edit] References

  1. ^ Born, M.; Wolf, E., Principles of optics: electromagnetic theory of propagation, interference and diffraction of light. Oxford, Pergamon Press, 1964.