# Fresnel's equation with complex n

• Nick89
In summary, The Fresnel equations for reflectivity for s- and p-polarization can be simplified by substituting values for a vacuum to metal interface, resulting in simpler expressions for R_s and R_p. However, finding a clear derivation or statement of these formulas has been difficult due to conflicting references and definitions of polarization. One possible derivation can be found in Chew's "Waves and Fields in Inhomogeneous Media."
Nick89
Hi,

I posted a similar question a while back but never got a decent answer, most probably due to having too many details that obscured the real problem. This time around I'll try it more simply.

Fresnel's equation for reflectivity for s- and p-polarization are, according to Wikipedia, given by:
$$R_s = \left( \frac{n_1 \cos \theta - n_2 \sqrt{1 - \left( \frac{n_1}{n_2} \sin \theta \right)^2}}{n_1 \cos \theta + n_2 \sqrt{1 - \left( \frac{n_1}{n_2} \sin \theta \right)^2}} \right)^2$$
$$R_p = \left( \frac{ n_1 \sqrt{1 - \left( \frac{n_1}{n_2}\sin\theta \right)^2} - n_2 \cos\theta }{ n_1 \sqrt{1 - \left( \frac{n_1}{n_2}\sin\theta \right)^2} + n_2 \cos\theta } \right)^2$$

Now suppose we have a vacuum to metal interface so that the light traveling through the vacuum strikes and reflects off a metal surface. I think in this case the vacuum should have an index of refraction equal to unity while the metal gets a complex index of refraction: $n_1 = 1$ and $n_2 = n(1 + i \kappa) = n + ik$ where $k$ is the extinction coefficient (or is it $\kappa$?).

Apparently, by substituting these values and making one assumption, that $n^2+k^2 \gg 1$, one should be able to arrive at the following much simpler expressions:
$$R_s = \frac{\left( n - \cos \theta \right)^2 + k^2}{\left( n + \cos \theta \right)^2 + k^2}$$
$$R_p = \frac{\left( n - \sec \theta \right)^2 + k^2}{\left( n + \sec \theta \right)^2 + k^2}$$

Graphing these formulas tells me that $R_p$ has a minimum reflectivity at some high angle, similar to non-metals (where this minimum actually goes to 0 for the Brewster angle) while $R_s$ is a smoothly increasing function.The problem is as follows... I have been looking everywhere for a derivation, or even a simple statement of these formulas but I haven't found a decent reference anywhere. The only times I can find solid references to the reflectivity of a metal is for normal incidence, for $\theta = 0$ so that both $\cos \theta = \sec \theta = 1$, so that doesn't help.
I did find a few other references, but they either contradict themselves (describe in detail that Rp has the minimum in reflectivity, then show a graph and write Rs in it ...?!) or contradict other papers (some seem to mention Rp has the minimum, some mention it is Rs).

Some papers / books also don't mention s- or p-polarization at all but merely TE or TM modes. The same problem here: they contradict each other, some papers say Rs = TE, others seem to say Rs = TM...The questions are simple (but maybe not to answer):

1. What's the correspondance between TE/TM mode and s-/p-polarization? Or is there no one definition and does it depend on something else..? I would think TM = p polarization and TE = s polarization, or vice versa. Which is it?
Same goes for perpendicular and parallel polarization. Perpendicular = TE or TM, s or p?

My view is this:
s (senkrecht) = perpendicular = TE
p (parallel) = parallel = TM.

Correct, or not?2. Are the final formulas I've given for Rs and Rp correct, or did I switch them around (formula for Rs labeled Rp and vice versa)?

3. Finally: how does one get to that result...? I don't see how the assumption that n^2+k^2 is large is any help, I've been trying for days to derive it but I am getting nowhere.
I did find one other 'derivation' (not in detail, which is what I want to see) which arrived at
$$R_p = \frac{(n^2+k^2) \cos^2 \theta - 2n\cos\theta + 1}{(n^2+k^2) \cos^2 \theta + 2n\cos\theta + 1}$$
after dividing out $\cos^2 \theta$ and completing the square that does yield $R_p$. The problem is that another statement (without derivation) I've got shows this to be $R_s$ instead?!
Any help is MUCH appreciated! Thanks!

You could probably glean far more insight into this problem by just simply deriving these equations from the boundary conditions. One text that discusses the reflection coefficients is Chew's "Waves and Fields in Inhomogeneous Media." For the case of reflection off of a planar boundary lying in the x-y plane, we can derive the following ODE:
$$\left[ \mu \frac{d}{dz} \mu^{-1} \frac{d}{dz} + \omega^2\mu\epsilon - k_x^2 \right] e_y = 0$$
$$\left[ \epsilon \frac{d}{dz} \epsilon^{-1} \frac{d}{dz} + \omega^2\mu\epsilon - k_x^2 \right] h_y = 0$$

where e_y is the amplitude of the TE wave and h_y is the amplitude for the TM wave that depend upon $e^{\pm ik_xx}$. It is easy to see from this that for constant permittivity and permeability that we have plane wave solutions fo the form $e^{\pm ik_zz}$ where $k_z = \sqrt{\omega^2\mu\epsilon - k_x^2}$. In addition, by inspection of the ODE we can deduce that the following boundary conditions must be made
$$e_{1y} = e_{2y}, \qquad \mu_1^{-1}\frac{d}{dz}e_{1y} = \mu_2^{-1}\frac{d}{dz} e_{2y}$$
$$h_{1y} = e_{2y}, \qquad \epsilon _1^{-1}\frac{d}{dz}e_{1y} = \epsilon _2^{-1}\frac{d}{dz} h_{2y}$$
From this you can derive the Fresnel reflection coefficients for any form of the permittivity and permeability,
$$R^{TE} = \frac{\mu_2k_{1z} - \mu_1k_{2z}}{\mu_2k_{1z} + \mu_1k_{2z}}$$
$$R^{TM} = \frac{\epsilon_2k_{1z} - \epsilon_1k_{2z}}{\epsilon_2k_{1z} + \epsilon_1k_{2z}}$$
What is s and p? I have no idea but you should be able to easily verify for yourself from the above TE and TM equations. $n = \sqrt{\frac{\mu\epsilon}{\mu_0\epsilon_0}}= \frac{k}{\omega\sqrt{\mu_0\epsilon_0}}$ relates the index of refraction with the permittivity and permeabilities above. Now \theta is the angle between the normal and the direction of incidence so $k_x = k\sin\theta$ and $k_z = k\cos\theta$ I believe. Barring the mathematical exercise, my guess is that p is TM since the Brewster angle only occurs for the TM wave if the material is nonmagnetic (\mu_1 = \mu_2).

But to your more salient point, I think the reason why you cannot reduce your equations to those given in your paper is because the equation from Wikipedia is WRONG when using complex numbers. It is obvious that if you plug in a complex index of refraction into the first set of equations that the result is complex. But reflectivity is real, it is the magnitude squared of the reflective coefficients so of course we are going to see some problem in just blind simple substitution. Once you account for this it comes out trivially. Just take the first two terms of the Taylor's expansion of the square root. But instead of squaring the numerator and denominator, multiply them by the complex conjugate of themselves so that you get the magnitude squared. Then it will fall out in place (I've already checked this by hand myself).

Last edited:

Born2bwire said:
But to your more salient point, I think the reason why you cannot reduce your equations to those given in your paper is because the equation from Wikipedia is WRONG when using complex numbers. It is obvious that if you plug in a complex index of refraction into the first set of equations that the result is complex. But reflectivity is real, it is the magnitude squared of the reflective coefficients so of course we are going to see some problem in just blind simple substitution. Once you account for this it comes out trivially. Just take the first two terms of the Taylor's expansion of the square root. But instead of squaring the numerator and denominator, multiply them by the complex conjugate of themselves so that you get the magnitude squared. Then it will fall out in place (I've already checked this by hand myself).

I have tried this but I'm still getting nowhere.
If I understand you correctly you are saying that I should approximate the square root by
$$\sqrt{1-\left(\frac{\sin \theta}{n+ik}\right)^2} \approx 1-\frac{1}{2} \left( \frac{\sin \theta}{n+ik} \right)^2$$

Then use
$$r = \frac{1-\frac{1}{2} \left( \frac{\sin \theta}{n+ik} \right)^2 - (n+ik) \cos \theta}{1-\frac{1}{2} \left( \frac{\sin \theta}{n+ik} \right)^2 + (n+ik) \cos \theta}$$

And finally
$$R = r \bar{r}$$

I have been trying this but still getting nowhere... Problem is that I cannot take the complex conjugate unless I get r in the form a+ib which I cannot see how...

Then I read your post again and saw that you said I should do this to the numerator and denominator separately, does that mean this?
$$x = 1-\frac{1}{2} \left( \frac{\sin \theta}{n+ik} \right)^2 - (n+ik) \cos \theta$$
$$y = 1-\frac{1}{2} \left( \frac{\sin \theta}{n+ik} \right)^2 + (n+ik) \cos \theta$$
$$r = \frac{x}{y}$$
$$R = \frac{x \bar{x}}{y \bar{y}}$$
That doesn't seem right to me mathematically... But even that couldn't get me anywhere.

As for the rest of your post, I will look into deriving them manually.

Whoops, sorry, I should have said that you take only the first term of the Taylor's expansion. The fact that n^2+k^2 >> 1 it looks obvious that we drop the second term (and working it out with the second term we find that it can be dropped in the end). So anyway, let's just look at the numerator for the reflectivity of s.
$$\left| n_1 \cos\theta - n_2\sqrt{1 - \left(\frac{n_1}{n_2}\sin\theta\right)^2} \right|^2 = \left| \cos\theta - \left(n+ik\right) \sqrt{1 - \left(\frac{1}{n+ik}\sin\theta\right)^2} \right|^2 \approx \left| \cos\theta - (n+ik) \right|^2$$
$$\left| \cos\theta - (n+ik) \right|^2 = \left(\cos\theta-(n+ik)\right)\left(\cos\theta-(n-ik)\right) = \cos^2\theta-2n\cos\theta+n^2+k^2 = \left(n-\cos\theta\right)^2+k^2$$
And you can work out the rest of them from that. Sorry about that mistake earlier. If you had taken the second term included you would get something like
$$\left(n-\cos\theta\right)^2+k^2 + \frac{\sin^2\theta}{n^2+k^2} \left[ k^2-n^2+n\cos\theta+\frac{1}{2}\sin^2\theta\right]$$
By virture of the fact that n^2+k^2 >> 1 we can drop these additional terms.

EDIT:
As for your talk about the complex conjugates and the like. Yes, you can split them up accordingly. Take for example that any complex number can be written as follows:
$$a+ib = Ae^{i\theta_1}$$
Now let us just multiply this number to our heart's content:
$$X = \frac{(a+ib)(c+id)}{e+if} = Ae^{i\theta_1}Be^{i\theta_2}C^{i\theta_3} = ABC e^{i(\theta_1+\theta_2+\theta_3)}$$
If we take the complex conjugate of the number X,
$$X^* = ABCe^{-i(\theta_1+\theta_2+\theta_3)} = Ae^{-i\theta_1}Be^{-i\theta_2}C^{-i\theta_3} = \frac{(a-ib)(c-id)}{e-if}$$
So we note that if we have a product or quotient of complex numbers, then the complex conjugate of the final number is equal to the product and quotient of the conjugates themselves.

Last edited:
Thanks... I finally got it :) At least now I am sure that Rp is the reflectivity with the minimum. The paper I found these formulas in actually had them reversed, they stated Rs as the formula with sec theta and Rp with cos theta... I knew they were wrong, I just couldn't prove it to myself because I couldn't find any other decent references or derivations. But now I can finally see how you can get to those approximations and I get exactly the opposite (Rp is with sec theta while Rs is with cos theta). Thanks!

Why do you not have a look in the relevant textbooks? E.g. Born/Wolf, Principles of Optics,7.Edition, Cambridge University Press 1999Capter 14.2-14.4.
Good luck H. Weber

## 1. What is Fresnel's equation with complex n?

Fresnel's equation with complex n is a mathematical formula that describes the reflection and transmission of electromagnetic waves at the interface between two different media with complex refractive indices. It takes into account both the real and imaginary components of the refractive index, which can affect the behavior of light at the interface.

## 2. How is Fresnel's equation with complex n different from the regular Fresnel's equation?

The regular Fresnel's equation only considers the real part of the refractive index, which assumes that the medium is purely transparent. However, in real-world scenarios, the refractive index can also have an imaginary component, which accounts for absorption of light. Fresnel's equation with complex n takes this into account and provides a more accurate description of light behavior at the interface.

## 3. What is the significance of the complex refractive index in Fresnel's equation?

The complex refractive index takes into account the absorption of light at the interface between two media. This is particularly important for materials that are not perfectly transparent, such as metals. The imaginary component of the refractive index represents the amount of light that is absorbed by the material.

## 4. How is Fresnel's equation with complex n used in practical applications?

Fresnel's equation with complex n is used in a variety of fields, including optics, materials science, and engineering. It is particularly useful in designing and optimizing optical devices, such as lenses, mirrors, and filters. It is also used in the study of thin films and coatings, as well as in the development of anti-reflective coatings.

## 5. Are there any limitations to Fresnel's equation with complex n?

Like any mathematical model, Fresnel's equation with complex n has its limitations. It assumes that the interface between the two media is smooth and flat, and that the incident light is a plane wave. In reality, interfaces can have imperfections and the incident light can have a more complicated wave pattern. In addition, this equation does not take into account the effects of multiple reflections or scattering, which can occur in some situations.

Replies
17
Views
621
Replies
33
Views
2K
Replies
3
Views
263
Replies
0
Views
494
Replies
5
Views
1K
Replies
2
Views
1K
Replies
1
Views
649
Replies
3
Views
2K
Replies
5
Views
2K
Replies
16
Views
2K