When I consult most resources, I get a set of equations for Fresnel reflection and transmission which look like this:
$$ r_\text{s} = \frac{n_1 \cos \theta_\text{i} - n_2 \cos \theta_\text{t}} {n_1 \cos \theta_\text{i} + n_2 \cos \theta_\text{t}} \\ t_\text{s} = \frac{2 n_1 \cos \theta_\text{i}} {n_1 \cos \theta_\text{i} + n_2 \cos \theta_\text{t}} \\ r_\text{p} = \frac{n_2 \cos \theta_\text{i} - n_1 \cos \theta_\text{t}} {n_2 \cos \theta_\text{i} + n_1 \cos \theta_\text{t}} \\ t_\text{p} = \frac{2 n_1 \cos \theta_\text{i}} {n_2 \cos \theta_\text{i} + n_1 \cos \theta_\text{t}} $$ Where,
- $r_\text{s/p}$ are the reflected intensities of perpendicular/parallel polarised light
- $t_\text{s/p}$ are the reflected transmissions of perpendicular/parallel polarised light
- $n_1$ and $n_2$ are the (potentially complex) indices of refraction
- $\theta_\text{i}$ is the angle of incidence
- $\theta_\text{t}$ is the angle of transmission
My aim is to convert these to usable Mueller matrices so that I can use them to transform Stokes vectors, which would determine the amount of light of each polarisation type which is refracted.
A ray tracer project I found had done this for reflection only, but their equations look considerably different:
$$ 2a^2 = \sqrt{(n^2 - \kappa^2 - \sin^2 \theta)^2 + 4 n^2 \kappa^2} + n^2 - \kappa^2 - \sin^2 \theta \\ 2b^2 = \sqrt{(n^2 - \kappa^2 - \sin^2 \theta)^2 + 4 n^2 \kappa^2} - n^2 + \kappa^2 + \sin^2 \theta \\ \\ F_\text{s} = \frac{a^2 + b^2 - 2 a \cos \theta + \cos^2 \theta} {a^2 + b^2 + 2 a \cos \theta + \cos^2 \theta} \\ F_\text{p} = \frac{a^2 + b^2 - 2 a \sin \theta \tan \theta + \sin^2 \theta \tan^2 \theta} {a^2 + b^2 + 2 a \sin \theta \tan \theta + \sin^2 \theta \tan^2 \theta} \\ δ_\text{s} = \arctan{\frac{2b \cos \theta} {\cos^2 \theta - a^2 - b^2}} \\ δ_\text{p} = \arctan{\frac{2 \cos \theta((n^2 - \kappa^2) b - (2 n \kappa) a} {(n^2 + \kappa^2)^2 \cos^2 \theta - a^2 - b^2}} \\ \\ A = \frac {F_\text{s} + F_\text{p}}{2}\\ B = \frac {F_\text{s} - F_\text{p}}{2}\\ C = \cos(\delta_\text{s} - \delta_\text{p})\sqrt{F_\text{s} F_\text{p}}\\ S = \sin(\delta_\text{s} - \delta_\text{p})\sqrt{F_\text{s} F_\text{p}}\\ \\ M = \begin{bmatrix} A & B & 0 & 0\\ B & A & 0 & 0\\ 0 & 0 & C & S\\ 0 & 0 & -S & S\\ \end{bmatrix} $$ Where,
- $\tilde{n} = n + i \kappa$ is the relative index of refraction ($= \frac{n_2}{n_1}$)
- $\theta$ is the angle of incidence
I guess the real question is,
How do I compute a Mueller matrix for the transmitted (refracted) ray?
Subtracting the reflected result from the incoming Stokes vector doesn't appear to be correct unless $\kappa = 0$.
Alternatively, how do I get from the two Fresnel equations for reflectance $r_\text{s}$ and $r_\text{p}$, to the four equations for reflectance $F_\text{s}$, $F_\text{p}$, $\delta_\text{s}$ and $\delta_\text{p}$? If I figure out how to do that, then I think I could apply the same logic to the equations for transmission and get a result from there.
Parts of this which I don't have issue with:
- The second set of equations only use the angle of incidence so they are clearly using Snell's Law and trigonometric identities to remove the angle of transmission from the equations. This does mean that the angles it's dealing with are complex angles, but I guess that's fine.
- The second set of equations splits the complex index of refraction into two parts so it has probably done that and then expanded the result. But there is no $i$ anywhere in the result so I assume they're doing something to the complex values when they get them back. Is it something like, converting to polar coordinates? 🤔
I attempted to follow the trail of sources to see if some paper spelled it out, but it really seems like it's papers all the way down. (The best lead I have is an expensive textbook which I can't confirm actually contains the answer.)