electromagnetismreflectionrefraction
I am trying to relate the equation for reflection coefficient in oblique mediums to the frequency but can't figure out how the frequency affects the reflection of light.
$n_1$= intrinsic impedance of medium 1
$n_2$= intrinsic impedance of medium 2
$\theta_1$ = incident angle
$\theta_2$= transmited angle
$$\Gamma= \frac{n_2 \cos(\theta_1) -n_1 \cos(\theta_2)}{n_2 \cos(\theta_1) + n_1 \cos(\theta_2)}$$
Can somebody explain to me how we can make this equation a function of frequency?
The phase-shift occurs due to the complex nature of the reflection coefficient. The so-called critical angle is given by:
$sin(\theta_M) = \frac{n_2}{n_1}$
As long as $\theta<\theta_M$ we have only partial reflection and a real valued reflection coefficient $R$.
As soon as the critical angle is exceeded $(\theta>\theta_M)$, we have $\mid R \mid=1$ and total reflection of the light occurs. $R$ is now complex and a phase shift is imposed on the reflected light, we write:
$R=e^{2j\phi}$
wehre $\phi$ is the phase shift between the incident and reflected wave.
Now, we know what the value of $R$ is, both for TE and for TM modes, i.e. the Fresnel coefficients for the two different states of polarization of the electric and magnetic fields. To get the formula for the phase-shift you only need to equal this complex expression for $R$ to the Fresnel formulas and solve for $\phi$. The formula you wrote is the one corresponding to TE modes, so it's the TE phase-shift (with electric fields perpendicular to the plane of incidence spanned by the wave normal and the normal to the interface). For TM modes you would get exactly the same result but multiplied by a factor of $n_{1}^{2}/n_{2}^{2}$.
Snell's law in its sine-over-sine form is not the best tool to handle situations that involve incidence angles above the critical angle, and the complex angles that you're getting will only lead you astray.
Instead, the way to see things is through an explicit evanescent wave, as in this recent answer of mine. When light refracts, you have an incoming wave of the form $$ \mathbf E_\mathrm{in}(\mathbf r,t) = \mathrm{Re}\mathopen{}\left( \mathbf E_{0,\mathrm{in}} e^{i(\mathbf k_\mathrm{in}\cdot \mathbf r- \omega t)}\right)\mathclose{}, $$ and on the other side you get a transmitted beam of the form $$ \mathbf E_\mathrm{tr}(\mathbf r,t) = \mathrm{Re}\mathopen{}\left( \mathbf E_{0,\mathrm{tr}} e^{i(\mathbf k_\mathrm{tr}\cdot \mathbf r- \omega t)}\right)\mathclose{}, $$ where Snell's law comes in through the matching of the two wavevectors, $\mathbf k_\mathrm{in}$ and $\mathbf k_\mathrm{tr}$. These are subject to two conditions:
The critical angle comes in when the right-hand side is too small for $k_{z,\mathrm{tr}}^2$ to be positive, i.e. when $n_1$ is so much bigger than $n_2$, and $k_x$ is so big, that the square root in $$k_{z,\mathrm{tr}} = \pm \frac{1}{n_1}\sqrt{ n_2^2k_{z,\mathrm{in}}^2 +(n_2^2-n_1^2)k_x^2 }$$ is no longer real.
If this is the case, it's not actually a problem: you just write $k_{z,\mathrm{tr}} = i\kappa_z$, and it gives you a perfectly workable 'transmitted' field, of the form $$ \mathbf E_\mathrm{tr}(\mathbf r,t) = \mathrm{Re}\mathopen{}\left( \mathbf E_{0,\mathrm{tr}} e^{-\kappa_{z} z}e^{i(k_x x- \omega t)}\right)\mathclose{}, $$ where what used to be an oscillatory $e^{ik_z z}$ has turned into the decaying exponential $e^{-\kappa_{z} z}$: this is the evanescent wave, which is the decaying trace of the light that's been unable to make it past the boundary, but which cannot drop discontinuously to zero.
Similarly, if you have a complex index of refraction, it means that your medium absorbs light, and there you have no alternative but to use explicit exponential decays in your description of the beams: if $n$ is complex, then the wavevector still needs to satisfy $$k^2 = n^2\omega^2,$$ which means that it too will be complex, and you'll have a decaying field in that medium.
In this situation, I would advise you to stay well away from a complex-angles formalism, and to keep things grounded in the cartesian components of the wavevector - or to be very, very careful with those complex angles if you insist on using them. The reason for this is that the cartesian components are related to the angle and the magnitude via \begin{align} k_{x,\mathrm{tr}} & = k_\mathrm{tr} \sin(\theta_\mathrm{tr}) \\ k_{z,\mathrm{tr}} & = k_\mathrm{tr} \cos(\theta_\mathrm{tr}), \end{align} and while it is tempting to make $\theta_\mathrm{tr}$ complex and leave it at that, it's important to note that $k_\mathrm{tr}$ also needs to be complex (because $k^2 = n^2\omega^2$), but that the two imaginary parts need to cancel out when you form $k_{x,\mathrm{tr}} = k_\mathrm{tr} \sin(\theta_\mathrm{tr})$, and that's not easy to visualize at all.
That said, you do have a good point that when the second medium has a nonzero absorbance, it is no longer easy to distinguish between the 'sub-critical' incidence angles and over-critical cases with total internal reflection - and indeed there is no longer any case with full reflection of the incident power, since the transmitted wave always has a nonzero $\mathrm{Re}(k_z)$ and therefore always carries up some amount of power away from the boundary.
The way to analyze this in a way that will bring out the relevant behaviour is to plot the wavevector component $k_z$ in the complex plane as the angle of incidence is shifted from $0$ to $\pi/2$. If $n_2$ is real, then $k_z$ will be real at normal incidence, it will approach the origin as $\theta_\mathrm{in}$ increases, and once it reaches the origin it will advance back up along the imaginary axis.
If $n_2$ has a nonzero imaginary part, though, then $k_z$ will follow a similar path but it will be detached: off with a nonzero imaginary part before reaching the origin, and then up the imaginary axis but with a positive real part. The image below shows this behaviour, tracking $k_z$ for $\theta_\mathrm{in}$ between $0$ and $\pi/2$, and with $n_2$ between $2$ (dark blue) to $2+i$ (light blue).
Mathematica source through Import["http://halirutan.github.io/Mathematica-SE-Tools/decode.m"]["http://i.stack.imgur.com/f5Zz6.png"]
As you can see, if $\mathrm{Im}(n_2)$ is very small (which it is in real materials), then you can define the critical angle as the 'kink' when $k_z$ is close to the origin. However, as that imaginary part increases, the concept of critical angle simply makes less and less sense: you just have a decaying wave away from the boundary, and the incidence angle has very little effect on how fast that wave decays. It mostly stops being a useful concept to begin with, so it doesn't make a lot of sense to dig deep into the mathematics to look for a precise definition of it in that regime.
Best Answer
The frequency (wavelength) dependence comes in through the refractive index, and it is specific to the medium. In the lab, the manufacturer of your optics will generally provide this dependence, as in e.g. this resource. Calculating this theoretically is generally a hard problem, though.