Here we systematically derive the above relations from electromagnetic premises.
Material parameters In order to compute meaningful Fresnel coefficients, we must assume that the medium is (approximately)
linear and
homogeneous. If the medium is also
isotropic, the four field vectors are
related by \begin{align} \mathbf{D} &= \epsilon \mathbf{E} \\ \mathbf{B} &= \mu \mathbf{H}\,, \end{align} where and are scalars, known respectively as the (electric)
permittivity and the (magnetic)
permeability of the medium. For vacuum, these have the values and , respectively. Hence we define the
relative permittivity (or
dielectric constant) , and the
relative permeability . In optics it is common to assume that the medium is non-magnetic, so that . For
ferromagnetic materials at radio/microwave frequencies, larger values of must be taken into account. But, for optically transparent media, and for all other materials at optical frequencies (except possible
metamaterials), is indeed very close to 1; that is, . In optics, one usually knows the
refractive index of the medium, which is the ratio of the
speed of light in vacuum () to the speed of light in the medium. In the analysis of partial reflection and transmission, one is also interested in the electromagnetic
wave impedance , which is the ratio of the amplitude of to the amplitude of . It is therefore desirable to express and in terms of and , and thence to relate to . The last-mentioned relation, however, will make it convenient to derive the reflection coefficients in terms of the wave
admittance , which is the reciprocal of the wave impedance . In the case of
uniform plane sinusoidal waves, the wave impedance or admittance is known as the
intrinsic impedance or admittance of the medium. This case is the one for which the Fresnel coefficients are to be derived.
Electromagnetic plane waves In a uniform plane sinusoidal
electromagnetic wave, the
electric field has the form {{NumBlk|:|\mathbf{E_k}e^{i(\mathbf{k\cdot r}-\omega t)},|}} where is the (constant) complex amplitude vector, is the
imaginary unit, is the
wave vector (whose magnitude is the angular
wavenumber), is the
position vector, is the
angular frequency, is time, and it is understood that the
real part of the expression is the physical field. The value of the expression is unchanged if the position varies in a direction normal to ; hence
is normal to the wavefronts. To advance the
phase by the angle
ϕ, we replace by (that is, we replace by ), with the result that the (complex) field is multiplied by . So a phase
advance is equivalent to multiplication by a complex constant with a
negative argument. This becomes more obvious when the field () is factored as , where the last factor contains the time-dependence. That factor also implies that differentiation w.r.t. time corresponds to multiplication by . If
ℓ is the component of in the direction of , the field () can be written . If the argument of is to be constant,
ℓ must increase at the velocity \omega/k\,,\, known as the
phase velocity . This in turn is equal to Solving for gives As usual, we drop the time-dependent factor , which is understood to multiply every complex field quantity. The electric field for a uniform plane sine wave will then be represented by the location-dependent
phasor {{NumBlk|:|\mathbf{E_k}e^{i\mathbf{k\cdot r}}.|}} For fields of that form,
Faraday's law and the
Maxwell-Ampère law respectively reduce to \begin{align} \omega\mathbf{B} &= \mathbf{k}\times\mathbf{E}\\ \omega\mathbf{D} &= -\mathbf{k}\times\mathbf{H}\,. \end{align} Putting and , as above, we can eliminate and to obtain equations in only and : \begin{align} \omega\mu\mathbf{H} &= \mathbf{k}\times\mathbf{E}\\ \omega\epsilon\mathbf{E} &= -\mathbf{k}\times\mathbf{H}\,. \end{align} If the material parameters and are real (as in a lossless dielectric), these equations show that form a
right-handed orthogonal triad, so that the same equations apply to the magnitudes of the respective vectors. Taking the magnitude equations and substituting from (), we obtain \begin{align} \mu cH &= nE\\ \epsilon cE &= nH\,, \end{align} where and are the magnitudes of and . Multiplying the last two equations gives {{NumBlk|:|n = c\,\sqrt{\mu\epsilon}\,.|}} Dividing (or cross-multiplying) the same two equations gives , where {{NumBlk|:|Y = \sqrt{\epsilon/\mu}\,.|}} This is the
intrinsic admittance. From () we obtain the phase velocity {{nowrap|1=c/n=1\big/\!\sqrt{\mu\epsilon\,}.}} For vacuum this reduces to {{nowrap|1=c=1\big/\!\sqrt{\mu_0\epsilon_0}.}} Dividing the second result by the first gives n=\sqrt{\mu_{\text{rel}}\epsilon_{\text{rel}}}\,. For a
non-magnetic medium (the usual case), this becomes {{tmath|1= n=\sqrt{\epsilon_{\text{rel} } } }}. Taking the reciprocal of (), we find that the intrinsic
impedance is {{nowrap|1=Z=\sqrt{\mu/\epsilon}.}} In vacuum this takes the value Z_0=\sqrt{\mu_0/\epsilon_0}\,\approx 377\,\Omega\,, known as the
impedance of free space. By division, {{nowrap|1=Z/Z_0=\sqrt{\mu_{\text{rel}}/\epsilon_{\text{rel}}}.}} For a
non-magnetic medium, this becomes Z=Z_0\big/\!\sqrt{\epsilon_{\text{rel}}}=Z_0/n.
Wave vectors In Cartesian coordinates , let the region have refractive index , intrinsic admittance , etc., and let the region have refractive index , intrinsic admittance , etc. Then the plane is the interface, and the axis is normal to the interface (see diagram). Let and (in bold
roman type) be the unit vectors in the and directions, respectively. Let the plane of incidence be the plane (the plane of the page), with the angle of incidence measured from towards . Let the angle of refraction, measured in the same sense, be , where the subscript stands for
transmitted (reserving for
reflected). In the absence of
Doppler shifts,
ω does not change on reflection or refraction. Hence, by (), the magnitude of the wave vector is proportional to the refractive index. So, for a given , if we
redefine as the magnitude of the wave vector in the
reference medium (for which ), then the wave vector has magnitude in the first medium (region in the diagram) and magnitude in the second medium. From the magnitudes and the geometry, we find that the wave vectors are \begin{align} \mathbf{k}_\text{i} &= n_1 k(\mathbf{i}\sin\theta_\text{i} + \mathbf{j}\cos\theta_\text{i})\\[.5ex] \mathbf{k}_\text{r} &= n_1 k(\mathbf{i}\sin\theta_\text{i} - \mathbf{j}\cos\theta_\text{i})\\[.5ex] \mathbf{k}_\text{t} &= n_2 k(\mathbf{i}\sin\theta_\text{t} + \mathbf{j}\cos\theta_\text{t})\\ &= k(\mathbf{i}\,n_1\sin\theta_\text{i} + \mathbf{j}\,n_2\cos\theta_\text{t})\,, \end{align} where the last step uses Snell's law. The corresponding
dot products in the phasor form () are {{NumBlk|:|\begin{align} \mathbf{k}_\text{i}\mathbf{\cdot r} &= n_1 k(x\sin\theta_\text{i} + y\cos\theta_\text{i})\\ \mathbf{k}_\text{r}\mathbf{\cdot r} &= n_1 k(x\sin\theta_\text{i} - y\cos\theta_\text{i})\\ \mathbf{k}_\text{t}\mathbf{\cdot r} &= k(n_1 x\sin\theta_\text{i} + n_2 y\cos\theta_\text{t})\,. \end{align}|}} Hence: {{NumBlk|:|At y=0\,,~~~\mathbf{k}_\text{i}\mathbf{\cdot r}=\mathbf{k}_\text{r}\mathbf{\cdot r}=\mathbf{k}_\text{t}\mathbf{\cdot r}=n_1 kx\sin\theta_\text{i}\,.|}}
s components For the
s polarization, the field is parallel to the axis and may therefore be described by its component in the direction. Let the reflection and transmission coefficients be and , respectively. Then, if the incident field is taken to have unit amplitude, the phasor form () of its -component is {{NumBlk|:|E_\text{i}=e^{i\mathbf{k}_\text{i}\mathbf{\cdot r}},|}} and the reflected and transmitted fields, in the same form, are {{NumBlk|:|\begin{align} E_\text{r} &= r_{s\,} e^{i\mathbf{k}_\text{r}\mathbf{\cdot r}}\\ E_\text{t} &= t_{s\,} e^{i\mathbf{k}_\text{t}\mathbf{\cdot r}}. \end{align}|}} Under the sign convention used in this article, a positive reflection or
transmission coefficient is one that preserves the direction of the
transverse field, meaning (in this context) the field normal to the plane of incidence. For the
s polarization, that means the field. If the incident, reflected, and transmitted fields (in the above equations) are in the -direction ("out of the page"), then the respective fields are in the directions of the red arrows, since form a right-handed orthogonal triad. The fields may therefore be described by their components in the directions of those arrows, denoted by . Then, since , {{NumBlk|:|\begin{align} H_\text{i} &=\, Y_1 e^{i\mathbf{k}_\text{i}\mathbf{\cdot r}}\\ H_\text{r} &=\, Y_1 r_{s\,} e^{i\mathbf{k}_\text{r}\mathbf{\cdot r}}\\ H_\text{t} &=\, Y_2 t_{s\,} e^{i\mathbf{k}_\text{t}\mathbf{\cdot r}}. \end{align}|}} At the interface, by the usual
interface conditions for electromagnetic fields, the tangential components of the and fields must be continuous; that is, {{NumBlk|:|\left.\begin{align} E_\text{i} + E_\text{r} &= E_\text{t}\\ H_\text{i}\cos\theta_\text{i} - H_\text{r}\cos\theta_\text{i} &= H_\text{t}\cos\theta_\text{t} \end{align}~~\right\}~~~\text{at}~~ y=0\,.|}} When we substitute from equations () to () and then from (), the exponential factors cancel out, so that the interface conditions reduce to the simultaneous equations {{NumBlk|:|\begin{align} 1 + r_\text{s} &=\, t_\text{s}\\ Y_1\cos\theta_\text{i} - Y_1 r_\text{s}\cos\theta_\text{i} &=\, Y_2 t_\text{s}\cos\theta_\text{t} \,, \end{align}|}} which are easily solved for and , yielding {{NumBlk|:|r_\text{s}=\frac{Y_1\cos\theta_\text{i}-Y_2\cos\theta_\text{t}}{Y_1\cos\theta_\text{i}+Y_2\cos\theta_\text{t}}|}} and {{NumBlk|:|t_\text{s}=\frac{2Y_1\cos\theta_\text{i}}{Y_1\cos\theta_\text{i}+Y_2\cos\theta_\text{t}}\,.|}} At
normal incidence , indicated by an additional subscript 0, these results become {{NumBlk|:|r_\text{s0}=\frac{Y_1-Y_2}{Y_1+Y_2}|}} and {{NumBlk|:|t_\text{s0}=\frac{2Y_1}{Y_1+Y_2}\,.|}} At
grazing incidence , we have , hence and .
p components For the
p polarization, the incident, reflected, and transmitted fields are parallel to the red arrows and may therefore be described by their components in the directions of those arrows. Let those components be (redefining the symbols for the new context). Let the reflection and transmission coefficients be and . Then, if the incident field is taken to have unit amplitude, we have {{NumBlk|:|\begin{align} E_\text{i} &= e^{i\mathbf{k}_\text{i}\mathbf{\cdot r}}\\ E_\text{r} &= r_{p\,} e^{i\mathbf{k}_\text{r}\mathbf{\cdot r}}\\ E_\text{t} &= t_{p\,} e^{i\mathbf{k}_\text{t}\mathbf{\cdot r}}. \end{align}|}} If the fields are in the directions of the red arrows, then, in order for to form a right-handed orthogonal triad, the respective fields must be in the -direction ("into the page") and may therefore be described by their components in that direction. This is consistent with the adopted sign convention, namely that a positive reflection or transmission coefficient is one that preserves the direction of the transverse field the field in the case of the
p polarization. The agreement of the
other field with the red arrows reveals an alternative definition of the sign convention: that a positive reflection or transmission coefficient is one for which the field vector in the plane of incidence points towards the same medium before and after reflection or transmission. So, for the incident, reflected, and transmitted fields, let the respective components in the -direction be . Then, since , {{NumBlk|:|\begin{align} H_\text{i} &=\, Y_1 e^{i\mathbf{k}_\text{i}\mathbf{\cdot r}}\\ H_\text{r} &=\, Y_1 r_{p\,} e^{i\mathbf{k}_\text{r}\mathbf{\cdot r}}\\ H_\text{t} &=\, Y_2 t_{p\,} e^{i\mathbf{k}_\text{t}\mathbf{\cdot r}}. \end{align}|}} At the interface, the tangential components of the and fields must be continuous; that is, {{NumBlk|:|\left.\begin{align} E_\text{i}\cos\theta_\text{i} - E_\text{r}\cos\theta_\text{i} &= E_\text{t}\cos\theta_\text{t}\\ H_\text{i} + H_\text{r} &= H_\text{t} \end{align}~~\right\}~~~\text{at}~~ y=0\,.|}} When we substitute from equations () and () and then from (), the exponential factors again cancel out, so that the interface conditions reduce to {{NumBlk|:|\begin{align} \cos\theta_\text{i} - r_\text{p}\cos\theta_\text{i} &=\, t_\text{p}\cos\theta_\text{t}\\ Y_1 + Y_1 r_\text{p} &=\, Y_2 t_\text{p} \,. \end{align}|}} Solving for and , we find {{NumBlk|:|r_\text{p}=\frac{Y_2\cos\theta_\text{i}-Y_1\cos\theta_\text{t}}{Y_2\cos\theta_\text{i}+Y_1\cos\theta_\text{t}}|}} and {{NumBlk|:|t_\text{p}=\frac{2Y_1\cos\theta_\text{i}}{Y_2\cos\theta_\text{i}+Y_1\cos\theta_\text{t}}\,.|}} At
normal incidence indicated by an additional subscript 0, these results become {{NumBlk|:|r_\text{p0}=\frac{Y_2-Y_1}{Y_2+Y_1}|}} and {{NumBlk|:|t_\text{p0}=\frac{2Y_1}{Y_2+Y_1}\,.|}} At , we again have , hence and . Comparing () and () with () and (), we see that at
normal incidence, under the adopted sign convention, the transmission coefficients for the two polarizations are equal, whereas the reflection coefficients have equal magnitudes but opposite signs. While this clash of signs is a disadvantage of the convention, the attendant advantage is that the signs agree at
grazing incidence.
Power ratios (reflectivity and transmissivity) The
Poynting vector for a wave is a vector whose component in any direction is the
irradiance (power per unit area) of that wave on a surface perpendicular to that direction. For a plane sinusoidal wave the Poynting vector is , where and are due
only to the wave in question, and the asterisk denotes complex conjugation. Inside a lossless dielectric (the usual case), and are in phase, and at right angles to each other and to the wave vector ; so, for s polarization, using the and components of and respectively (or for p polarization, using the and components of and ), the
irradiance in the direction of is given simply by , which is in a medium of intrinsic impedance . To compute the irradiance in the direction normal to the interface, as we shall require in the definition of the power transmission coefficient, we could use only the component (rather than the full component) of or or, equivalently, simply multiply by the proper geometric factor, obtaining . From equations () and (), taking squared magnitudes, we find that the
reflectivity (ratio of reflected power to incident power) is {{NumBlk|:|R_\text{s}=\left|\frac{Y_1\cos\theta_\text{i}-Y_2\cos\theta_\text{t}}{Y_1\cos\theta_\text{i}+Y_2\cos\theta_\text{t}}\right|^2|}} for the s polarization, and {{NumBlk|:|R_\text{p}=\left|\frac{Y_2\cos\theta_\text{i}-Y_1\cos\theta_\text{t}}{Y_2\cos\theta_\text{i}+Y_1\cos\theta_\text{t}}\right|^2|}} for the p polarization. Note that when comparing the powers of two such waves in the same medium and with the same cos
θ, the impedance and geometric factors mentioned above are identical and cancel out. But in computing the power
transmission (below), these factors must be taken into account. The simplest way to obtain the power transmission coefficient (
transmissivity, the ratio of transmitted power to incident power
in the direction normal to the interface, i.e. the direction) is to use (conservation of energy). In this way we find {{NumBlk|:|T_\text{s} =1-R_\text{s} =\,\frac{4\,\text{Re}\{ Y_1 Y_2\cos\theta_\text{i}\cos\theta_\text{t}\}}{\left|Y_1\cos\theta_\text{i}+Y_2\cos\theta_\text{t}\right|^2}|}} for the s polarization, and {{NumBlk|:|T_\text{p} =1-R_\text{p} =\,\frac{4\,\text{Re}\{Y_1 Y_2\cos\theta_\text{i}\cos\theta_\text{t}\}}{\left|Y_2\cos\theta_\text{i}+Y_1\cos\theta_\text{t}\right|^2}|}} for the p polarization. In the case of an interface between two lossless media (for which ϵ and μ are
real and positive), one can obtain these results directly using the squared magnitudes of the amplitude transmission coefficients that we found earlier in equations () and (). But, for given amplitude (as noted above), the component of the Poynting vector in the direction is proportional to the geometric factor and inversely proportional to the wave impedance . Applying these corrections to each wave, we obtain two ratios multiplying the square of the amplitude transmission coefficient: {{NumBlk|:|T_\text{s} =\left(\frac{2Y_1\cos\theta_\text{i}}{Y_1\cos\theta_\text{i}+Y_2\cos\theta_\text{t}}\right)^2\frac{\,Y_2\,}{Y_1}\,\frac{\cos\theta_\text{t}}{\cos\theta_\text{i}} =\frac{4Y_1 Y_2\cos\theta_\text{i}\cos\theta_\text{t}}{\left(Y_1\cos\theta_\text{i}+Y_2\cos\theta_\text{t}\right)^2}|}} for the s polarization, and {{NumBlk|:|T_\text{p} =\left(\frac{2Y_1\cos\theta_\text{i}}{Y_2\cos\theta_\text{i}+Y_1\cos\theta_\text{t}}\right)^2\frac{\,Y_2\,}{Y_1}\,\frac{\cos\theta_\text{t}}{\cos\theta_\text{i}} =\frac{4Y_1 Y_2\cos\theta_\text{i}\cos\theta_\text{t}}{\left(Y_2\cos\theta_\text{i}+Y_1\cos\theta_\text{t}\right)^2}|}} for the p polarization. The last two equations apply only to lossless dielectrics, and only at incidence angles smaller than the critical angle (beyond which, of course, ). For unpolarized light: T={1 \over 2}(T_s+T_p) R={1 \over 2}(R_s+R_p) where R+T=1.
Equal refractive indices From equations () and (), we see that two dissimilar media will have the same refractive index, but different admittances, if the ratio of their permeabilities is the inverse of the ratio of their permittivities. In that unusual situation we have (that is, the transmitted ray is undeviated), so that the cosines in equations (), (), (), (), and () to () cancel out, and all the reflection and transmission ratios become independent of the angle of incidence; in other words, the ratios for normal incidence become applicable to all angles of incidence. When extended to spherical reflection or scattering, this results in the Kerker effect for
Mie scattering.
Non-magnetic media Since the Fresnel equations were developed for optics, they are usually given for non-magnetic materials. Dividing () by ()) yields Y=\frac{n}{\,c\mu\,}\,. For non-magnetic media we can substitute the
vacuum permeability for , so that Y_1=\frac{n_1}{\,c\mu_0} ~~;~~~ Y_2=\frac{n_2}{\,c\mu_0}\,; that is, the admittances are simply proportional to the corresponding refractive indices. When we make these substitutions in equations () to () and equations () to (), the factor
cμ0 cancels out. For the amplitude coefficients we obtain: Similarly, Fresnel's sine law will apply to the p polarization instead of the s polarization, and his tangent law to the s polarization instead of the p polarization. This switch of polarizations has an analog in the old mechanical theory of light waves (see
§History, above). One could predict reflection coefficients that agreed with observation by supposing (like Fresnel) that different refractive indices were due to different
densities and that the vibrations were
normal to what was then called the
plane of polarization, or by supposing (like
MacCullagh and
Neumann) that different refractive indices were due to different
elasticities and that the vibrations were
parallel to that plane. Thus the condition of equal permittivities and unequal permeabilities, although not realistic, is of some historical interest. == See also ==