Jashcraf / poke

Poke (pronounced poh-keh) is a Polarization Ray Tracing and Gaussian Beamlet module for Python
BSD 3-Clause "New" or "Revised" License
32 stars 6 forks source link

1.2 Assess need for Berreman calculus #8

Open Jashcraf opened 1 year ago

Jashcraf commented 1 year ago

Birefringent Thin Films and Polarizing Elements

Section 3.3: Propagation in Layered Biaxial Media

A plane wave incident on a parallel layer of biaxial medium will excite four plane waves in the medium. These waves are linearly polarized in directions specified by the D fields and share the same value $\beta$ for snells law

$$\beta = nsin(\theta)$$

Where the angles change based on the effective refractive index experienced. We also know that there are four parameters $\alpha$ defined by

$$\alpha = ncos\theta$$

amounts to knowledge of the four indices and angles, because

$$ n = \sqrt{\alpha^2 + \beta^2} $$ $$ \theta = arcsin(\beta/n) $$

These are used to solve Fresnel's quartic equation, which is kind of a beast

image

But there are simplifications for Uniaxial and Isotropic media later in the chapter. To determine the basis vectors you need to solve for the Auxilliary Matrix L

image

Which I think is just a function of the layer thickness and electric permitivity. Computing the field matrix is done by solving for the eigenvectors of the auxilliary matrix.

image

Jashcraf commented 1 year ago

Berreman Calculus

Section 4.4: Berreman Calculus

Armed with the field matrix you have the basis vectors for your four principal directions. Note that the eigenequation can't be used for the isotropic case (but we would just call the characteristic matrix in this case). Beginning with the field matrix and the field coeficients, we can compute the total field vector.

image

The four waves also change phase linearly with displacement but at different rates. We account for this with a diagonal phase matrix.

image

Computing the characteristic matrix is done like a rotation of the phase matrix by the field matrix

image

Through which the m vector can be propagated

image

For an arbitrary set of thin films this is done with

image

Repeating this procedure for the full characteristic matrix with the field matrix for the cover medium and the substrate medium to derive the system matrix A.

image