Ewald–Oseen extinction theorem explained

In optics, the Ewald–Oseen extinction theorem, sometimes referred to as just the extinction theorem, is a theorem that underlies the common understanding of scattering (as well as refraction, reflection, and diffraction). It is named after Paul Peter Ewald and Carl Wilhelm Oseen, who proved the theorem in crystalline and isotropic media, respectively, in 1916 and 1915. Originally, the theorem applied to scattering by an isotropic dielectric objects in free space. The scope of the theorem was greatly extended to encompass a wide variety of bianisotropic media.

Overview

An important part of optical physics theory is starting with microscopic physics—the behavior of atoms and electrons—and using it to derive the familiar, macroscopic, laws of optics. In particular, there is a derivation of how the refractive index works and where it comes from, starting from microscopic physics. The Ewald–Oseen extinction theorem is one part of that derivation (as is the Lorentz–Lorenz equation etc.).

When light traveling in vacuum enters a transparent medium like glass, the light slows down, as described by the index of refraction. Although this fact is famous and familiar, it is actually quite strange and surprising when you think about it microscopically. After all, according to the superposition principle, the light in the glass is a superposition of:

• The original light wave, and
• The light waves emitted by oscillating electrons in the glass.

(Light is an oscillating electromagnetic field that pushes electrons back and forth, emitting dipole radiation.)

Individually, each of these waves travels at the speed of light in vacuum, not at the (slower) speed of light in glass. Yet when the waves are added up, they surprisingly create only a wave that travels at the slower speed.

The Ewald–Oseen extinction theorem says that the light emitted by the atoms has a component traveling at the speed of light in vacuum, which exactly cancels out ("extinguishes") the original light wave. Additionally, the light emitted by the atoms has a component which looks like a wave traveling at the slower speed of light in glass. Altogether, the only wave in the glass is the slow wave, consistent with what we expect from basic optics.

A more complete description can be found in Classical Optics and its Applications, by Masud Mansuripur. A proof of the classical theorem can be found in Principles of Optics, by Born and Wolf., and that of its extension has been presented by Akhlesh Lakhtakia.

Derivation from Maxwell's equations

Introduction

When an electromagnetic wave enters a dielectric medium, it excites (resonates) the material's electrons whether they are free or bound, setting them into a vibratory state with the same frequency as the wave. These electrons will in turn radiate their own electromagnetic fields as a result of their oscillation (EM fields of oscillating charges). Due to the linearity of Maxwell equations, one expects the total field at any point in space to be the sum of the original field and the field produced by oscillating electrons. This result is, however, counterintuitive to the practical wave one observes in the dielectric moving at a speed of c/n, where n is the medium index of refraction. The Ewald–Oseen extinction theorem seek to address the disconnect by demonstrating how the superposition of these two waves reproduces the familiar result of a wave that moves at a speed of c/n.

Derivation

The following is a derivation based on a work by Ballenegger and Weber.^[1] Let's consider a simplified situation in which a monochromatic electromagnetic wave is normally incident on a medium filling half the space in the region z>0 as shown in Figure 1.The electric field at a point in space is the sum of the electric fields due to all the various sources. In our case, we separate the fields in two categories based on their generating sources. We denote the incident field$$\backslash mathbf\_$$and the sum of the fields generated by the oscillating electrons in the medium$$\backslash mathbf\_(z,\; t).$$

The total field at any point z in space is then given by the superposition of the two contributions,$$\backslash mathbf(z,\; t)=\; \backslash mathbf\_(z,\; t)\; +\; \backslash mathbf\_(z,\; t).$$

To match what we already observe,

E _{vac} has this form. However, we already know that inside the medium, z>0, we will only observe what we call the transmitted E-field E _{T} which travels through the material at speed c/n.

Therefore in this formalism,$$\backslash mathbf\_(z,\; t)\; =\; -\backslash mathbf\_(z,\; t)\; +\; \backslash mathbf\_(z,\; t)$$

This to say that the radiated field cancels out the incident field and creates a transmitted field traveling within the medium at speed c/n. Using the same logic, outside the medium the radiated field produces the effect of a reflected field

E _{R} traveling at speed c in the opposite direction to the incident field.$$\backslash mathbf\_(z,\; t)\; =\; -\backslash mathbf\_(z,\; t)\; -\backslash mathbf\_(z,\; t)$$assume that the wavelength is much larger than the average separation of atoms so that the medium can be considered continuous. We use the usual macroscopic E and B fields and take the medium to be nonmagnetic and neutral so that Maxwell's equations readboth the total electric and magnetic fields$$\backslash mathbf=\backslash mathbf\_+\backslash mathbf\_,\; \backslash quad\; \backslash mathbf=\backslash mathbf\_+\backslash mathbf\_$$the set of Maxwell equations inside the dielectric$$\backslash begin\; \backslash \backslash \; \backslash \backslash \; \backslash \backslash \; \backslash end$$where J includes the true and polarization current induced in the material by the outside electric field. We assume a linear relationship between the current and the electric field, hence$$\backslash mathbf=\; \backslash left(\backslash mathbf\_\; +\; \backslash mathbf\_\backslash right)$$

The set of Maxwell equations outside the dielectric has no current density term$$\backslash begin\; \backslash \backslash \; \backslash \backslash \; \backslash \backslash \backslash end$$

The two sets of Maxwell equations are coupled since the vacuum electric field appears in the current density term.

For a monochromatic wave at normal incidence, the vacuum electric field has the form$$\backslash mathbf\_(z,\; t)=\backslash mathbf\_\; \backslash exp\; [i(k\; z-\; \backslash omega\; t)],$$with

k = \omega /{c} .

Now to solve for

E _{rad} , we take the curl of the third equation in the first set of Maxwell equation and combine it with the fourth.

We simplify the double curl in a couple of steps using Einstein summation.

Hence we obtain,$$\backslash nabla\backslash times\backslash nabla\; \backslash times\; \backslash mathbf\_=\backslash nabla(\backslash nabla\backslash cdot\backslash mathbf\_)-\backslash nabla^2\backslash mathbf\_$$

Then substituting

J by {\sigma} \left( E _{vac} + E _{rad} \right) , using the fact that \nabla ⋅ E _{rad} = 0 we obtain,$$\backslash nabla^2\backslash mathbf\_\backslash mathrm\; =\; \backslash frac\; \backslash partial\; (\backslash mu\_\backslash mathbf\_\; +\; \backslash mu\_0\; \backslash mathbf\_\; +\; \backslash epsilon\_\; \backslash mu\_\; \backslash partial\; \backslash mathbf\_\; /\; \backslash partial\; t)$$

Realizing that all the fields have the same time dependence

\exp ( - i \omega t) , the time derivatives are straightforward and we obtain the following inhomogeneous wave equation$$\backslash nabla^\; \backslash mathbf\_+\backslash mu\_\; \backslash omega^\backslash left(\backslash epsilon\_+i\; \backslash sigma\; /\; \backslash omega\backslash right)\; \backslash mathbf\_=\; -i\; \backslash mu\_\; \backslash omega\; \backslash sigma\; \backslash mathbf\_(z)$$with particular solution$$\backslash mathbf\_^\; =\; -\backslash mathbf\_(z)$$

For the complete solution, we add to the particular solution the general solution of the homogeneous equation which is a superposition of plane waves traveling in arbitrary directions$$\backslash left(\backslash mathbf\_^\backslash right)\_\; =\; \backslash int\; g\_(\backslash boldsymbol,\; \backslash boldsymbol)\; \backslash exp\; \backslash left(i\; \backslash mathbf\text{'}\; \backslash cdot\; \backslash mathbf\backslash right)\; d\; \backslash Omega$$where

k' is found from the homogeneous equation to be$$k^=\backslash mu\_\; \backslash epsilon\_\; \backslash omega^2\; \backslash left(1+i\; \backslash frac\backslash right)$$

Note that we have taken the solution as a coherent superposition of plane waves. Because of symmetry, we expect the fields to be the same in a plane perpendicular to the

z axis. Hence k ' ⋅ a = 0 , where a is a displacement perpendicular to z .

Since there are no boundaries in the region

z> 0 , we expect a wave traveling to the right. The solution to the homogeneous equation becomes,$$\backslash mathbf\_^\; =\; \backslash mathbf\_\; \backslash exp\; \backslash left(i\; k\text{'}\; z\backslash right)$$

Adding this to the particular solution, we get the radiated wave inside the medium (

z > 0 )$$\backslash mathbf\_\; =\; -\backslash mathbf\_(z)\; +\; \backslash mathbf\_\; \backslash exp\; \backslash left(i\; k\text{'}\; z\backslash right)$$

The total field at any position

z is the sum of the incident and radiated fields at that position. Adding the two components inside the medium, we get the total field$$\backslash mathrm(z)\; =\; \backslash mathrm\_\; \backslash exp\; \backslash left(i\; k\text{'}\; z\backslash right),\; \backslash qquad\; z>0$$

This wave travels inside the dielectric at speed

c/n, $$n\; =\; c\; k\text{'}\; /\; \backslash omega\; =\; \backslash sqrt$$

We can simplify the above

n to a familiar form of the index of refraction of a linear isotropic dielectric. To do so, we remember that in a linear dielectric an applied electric field E induces a polarization P proportional to the electric field P = \epsilon _{0} \chi _{e} E . When the electric field changes, the induced charges move and produces a current density given by \partial P / \partial t . Since the time dependence of the electric field is \exp( - i\omega t) , we get$$\backslash mathbf=-i\; \backslash epsilon\_\; \backslash omega\; \backslash chi\_\; \backslash mathbf,$$which implies that the conductivity$$\backslash sigma=-i\; \backslash epsilon\_\; \backslash omega\; \backslash chi\_.$$

Then substituting the conductivity in the equation of

n , gives$$n=\backslash sqrt$$which is a more familiar form. For the region z< 0 , one imposes the condition of a wave traveling to the left. By setting the conductivity in this region \sigma =0 , we obtain the reflected wave$$\backslash mathrm(z)\; =\; \backslash mathrm\_\; \backslash exp\; \backslash left(-i\; k\; z\backslash right),$$traveling at the speed of light.

Note that the coefficients nomenclature,

E _{T} and E _{R} , are only adopted to match what we already expect.

Hertz vector approach

The following is a derivation based on a work by Wangsness ^[2] and a similar derivation found in chapter 20 of Zangwill's text, Modern Electrodynamics.^[3] The setup is as follows, let the infinite half-space

z< 0 be vacuum and the infinite half-space z> 0 be a uniform, isotropic, dielectric material with electric susceptibility, \chi.

The inhomogeneous electromagnetic wave equation for the electric field can be written in terms of the electric Hertz Potential,

\boldsymbol{\pi} _{e} , in the Lorenz gauge as$$\backslash nabla^\; \backslash boldsymbol\_-\backslash frac\; \backslash frac=-\backslash frac.$$

The electric field in terms of the Hertz vectors is given as$$\backslash mathbf=\backslash nabla\; \backslash times\; \backslash nabla\; \backslash times\; \backslash boldsymbol\backslash pi\_-\backslash frac-\backslash frac\backslash left(\backslash nabla\; \backslash times\; \backslash boldsymbol\_\backslash right),$$but the magnetic Hertz vector

\boldsymbol{\pi} _{m} is 0 since the material is assumed to be non-magnetizable and there is no external magnetic field. Therefore the electric field simplifies to$$\backslash mathbf\; =\; \backslash nabla\; \backslash times\; \backslash nabla\; \backslash times\; \backslash boldsymbol\backslash pi\_\; -\; \backslash frac.$$

In order to calculate the electric field we must first solve the inhomogeneous wave equation for

\boldsymbol{\pi} _{e} . To do this, split \boldsymbol{\pi} _{e} in the homogeneous and particular solutions$$\backslash boldsymbol\_(\backslash mathbf,\; t)\; =\; \backslash boldsymbol\_(\backslash mathbf,\; t)\; +\; \backslash boldsymbol\_(\backslash mathbf,\; t).$$

Linearity then allows us to write$$\backslash mathbf(\backslash mathbf,\; t)=\backslash mathbf\_(\backslash mathbf,\; t)+\backslash mathbf\_(\backslash mathbf,\; t).$$

The homogeneous solution,

E _{h} ( r , t) , is the initial plane wave traveling with wave vector k _{0} = \omega/c in the positive z direction$$\backslash mathbf\_(\backslash mathbf,\; t)\; =\; \backslash mathbf\_\; e^\; .$$

We do not need to explicitly find

\boldsymbol{\pi} _{e,h} ( r , t) since we are only interested in finding the field.

The particular solution,

\boldsymbol{\pi} _{e,p} ( r , t) and therefore, E _{p} ( r , t) , is found using a time dependent Green's function method on the inhomogeneous wave equation for \boldsymbol{\pi} _{e,p} which produces the retarded integral$$\backslash boldsymbol\_(\backslash mathbf,\; t)\; =\; \backslash frac\; \backslash int\; d^3\; r\text{'}\; \backslash frac.$$

Since the initial electric field is polarizing the material, the polarization vector must have the same space and time dependence

P ( r , t) = P _{0} e ^{i(k} . More detail about this assumption is discussed by Wangsness. Plugging this into the integral and expressing in terms of Cartesian coordinates produces$$\backslash boldsymbol\_(\backslash mathbf,\; t)=\; \backslash frac\; \backslash int\_^\; dz\text{'}\; e^\; \backslash int\_^\; dx\text{'}\; \backslash int\_^\; dy\text{'}\; \backslash frac.$$

First, consider only the integration over

x' and y' and convert this to cylindrical coordinates (x,y,z) → (\rho,\varphi,z) and call \left| r - r '\right| = R $$I:=\backslash int\_^\; d\; x^\; \backslash int\_^\; d\; y^\; \backslash frac=\backslash int\_^\; d\; \backslash varphi\text{'}\; \backslash int\_^\; d\; \backslash rho\text{'}\; \backslash frac=2\backslash pi\; \backslash int\_^\; d\; \backslash rho\text{'}\; \backslash frac.$$

Then using the substitution$$R^=\backslash rho^+\backslash left|z\text{'}-z\backslash right|^\; \backslash Rightarrow\; \backslash rho^\; =\; R^-\backslash left|z\text{'}-z\backslash right|^\; \backslash Rightarrow\; \backslash rho\; \backslash ,\; d\; \backslash rho\; =\; R\; \backslash ,\; dR$$and$$\backslash rho=\backslash sqrt$$so the limits become$$\backslash rho=0=\backslash sqrt\; \backslash Rightarrow\; R=\backslash left|z\text{'}-z\backslash right|$$and$$\backslash rho\; =\; \backslash infty=\backslash sqrt\; \backslash Rightarrow\; R=\backslash infty.$$

Then introduce a convergence factor

e ^{-\epsilon} with \epsilon \in \R into the integrand since it does not change the value of the integral,

Then

\epsilon \in \R implies \lim _{\epsilon} e ^{-\epsiloninfty} =0 , hence \lim _{\epsilon} e ^{(i k _{0} - \epsilon) inf ty} = \lim _{\epsilon} e ^{i k _{0} inf ty} e ^{-\epsiloninfty} =0 . Therefore,

\end

Now, plugging this result back into the z-integral yields$$\backslash boldsymbol\_(z,\; t)=\; \backslash frac\; \backslash int\_^\; d\; z\text{'}$$

e^

e^

Notice that

\boldsymbol{\pi} _{e,p} is now only a function of z and not r , which was expected for the given symmetry.

This integration must be split into two due to the absolute value

\left|z - z' \right| inside the integrand. The regions are z< 0 and z> 0 . Again, a convergence factor must be introduced to evaluate both integrals and the result is$$\backslash boldsymbol\_(z,\; t)=-\backslash frac$$

\begin & \\ & \end

Instead of plugging

\boldsymbol{\pi} _{e,p} directly into the expression for the electric field, several simplifications can be made. Begin with the curl of the curl vector identity,$$\backslash nabla\; \backslash times(\backslash nabla\; \backslash times\; \backslash mathbf)=\backslash nabla(\backslash nabla\; \backslash cdot\; \backslash mathbf)-\backslash nabla^\; \backslash mathbf,$$therefore,

Notice that

\nabla ⋅ \boldsymbol{\pi} _{e,p} = 0 because P has no { z } dependence and is always perpendicular to \hat{ z }. Also, notice that the second and third terms are equivalent to the inhomogeneous wave equation, therefore,

Therefore, the total field is$$\backslash mathbf(z,\; t)=\backslash mathbf\_\; e^+k\_^\; \backslash boldsymbol\backslash pi\_(z,\; t)$$which becomes,$$\backslash mathbf(z,\; t)=\backslash left\backslash$$

Notes and References

1. Ballenegger. Vincent C.. Weber. T. A.. 1999-07-01. The Ewald–Oseen extinction theorem and extinction lengths. American Journal of Physics . 67. 7. 599–605 . 10.1119/1.19330. 0002-9505.
2. Wangsness. Roald K.. 1981-10-01. Effect of matter on the phase velocity of an electromagnetic wave. American Journal of Physics. 49. 10. 950–953. 10.1119/1.12596. 0002-9505. 1981AmJPh..49..950W.
3. Book: Zangwill, Andrew. Modern Electrodynamics. Cambridge University Press. 2013. 9780521896979.