Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

In this chapter, we return to (9.46)-(9.48) and examine the consequences of Maxwell’s equations in a homogeneous material for a general traveling electromagnetic plane wave. The extra complication is polarization.

12.1: The String in Three Dimensions

In most of our discussions of wave phenomena so far, we have assumed that the motion is taking place in a plane, so that we can draw pictures of the system on a sheet of paper. We have implicitly been restricting ourselves to two-dimensional waves. This is all right for longitudinal oscillations in three dimensions, because all the action is taking place along a single line. However, for transverse oscillations, going from two dimensions to three dimensions makes an enormous difference because there are two transverse directions in which the system can oscillate.

For example, consider a string in three dimensions, stretched in the zz direction. Each point on the string can oscillate in both the xx direction and the yy direction. If the system were not approximately linear, this could be a horrendous problem. Linearity allows us to solve the problem of oscillation in the xx-zz plane separately from the problem of oscillation in the yy-zz plane. We have already solved these two-dimensional problems in chapter 5. Then we can simply put the results together to get the most general motion of the three-dimensional system. In other words, we can treat the xx component of the transverse oscillation and the yy component as completely independent.

Suppose that there is a harmonic traveling wave in the +zz direction in the string. The displacement of the string at zz from its equilibrium position, (0,0,z)(0, 0, z), can be written as

Ψ(z,t)=Re[(ψ1x^+ψ2y^)ei(kzωt)](12.1)\vec{\Psi}(z, t)=\operatorname{Re}\left[\left(\psi_{1} \hat{x}+\psi_{2} \hat{y}\right) e^{i(k z-\omega t)}\right] \tag{12.1}

where x^\hat{x} and y^\hat{y} are unit vectors in the xx and yy direction and ψ1\psi_{1} and ψ2\psi_{2} are complex parameters describing the amplitude and phase of the oscillations in the xx-zz plane and the yy-zz plane,

ψj=Ajeiϕj for j=1 to 2.(12.2)\psi_{j}=A_{j} e^{i \phi_{j}} \text { for } j=1 \text { to } 2. \tag{12.2}

It is convenient to arrange these parameters into a complex vector

Z=(ψ1ψ2),(12.3)Z=\left(\begin{array}{l} \psi_{1} \\ \psi_{2} \end{array}\right) , \tag{12.3}

which gives a complete description of the motion of the string.

Polarization

Figure

“Polarization” refers to the nature of the motion of a point on the string (or other transverse oscillation). This motion is animated in program 12-1. You may want to read the discussion below with this program running.

If ϕ1=ϕ2\phi_{1}=\phi_{2}, or A1A_{1} or A2A_{2} is zero, then 12.3 represent a linearly polarized string. Linear polarization is easy to understand. It means that each point on the string is oscillating back and forth in a fixed plane. For example,

u1=(10)(12.4)u_{1}=\left(\begin{array}{l} 1 \\ 0 \end{array}\right) \tag{12.4}
u2=(01)(12.5)u_{2}=\left(\begin{array}{l} 0 \\ 1 \end{array}\right) \tag{12.5}

represent strings oscillating in the xx-zz plane and the yy-zz plane respectively. A string oscillating in a plane an angle θ\theta from the positive xx axis (towards the positive yy axis) is represented by

uθ=(cosθsinθ).(12.6)u_{\theta}=\left(\begin{array}{l} \cos \theta \\ \sin \theta \end{array}\right) \tag{12.6} .

This is shown in the xx-yy plane in Figure 12.1. The polarization vectors 12.4-12.6 can be multiplied by a phase factor, eiϕe^{i \phi}, without affecting the polarization state in any important way. This just corresponds to an overall resetting of the clock.

u_{1}, u_{2} and u_{\theta}.

Figure 12.1:u1u_{1}, u2u_{2} and uθu_{\theta}.

More interesting is circular polarization. A circularly polarized wave in a string is represented by either

(1i)(12.7)\left(\begin{array}{l} 1 \\ i \end{array}\right) \tag{12.7}

or

(1i).(12.8)\left(\begin{array}{c} 1 \\ -i \end{array}\right) . \tag{12.8}

In 12.7, the yy component lags behind the xx component by π/2(=ϕ2)\pi / 2\left(=\phi_{2}\right). Thus, at any fixed point in space, the field rotates from xx to yy, or in the counterclockwise direction viewed from the positive zz axis (with the wave coming at you), as shown in Figure 12.2. This is called “left-circular polarization” because the string resembles a left-handed screw. Likewise,

(1i).(12.9)\left(\begin{array}{c} 1 \\ -i \end{array}\right) . \tag{12.9}

represents clockwise rotation of the string. This is called “right-circular polarization.”

Circular polarization.

Figure 12.2:Circular polarization.

The vector

(AiB)(12.10)\left(\begin{array}{c} A \\ i B \end{array}\right) \tag{12.10}

with A>B>0A > B > 0 represents elliptical polarization. A point on the string traces out an ellipse with semi-major axis AA along the 1 axis and semi-minor axis BB along the 2 axis, with counterclockwise rotation, as shown in Figure 12.3

Elliptical polarization with long axis in the x direction.

Figure 12.3:Elliptical polarization with long axis in the xx direction.

General elliptical polarization.

Figure 12.4:General elliptical polarization.

A completely general vector can be written in the following form:

(ψ1ψ2)=eiϕ(AcosθiBsinθAsinθ+iBcosθ)(12.11)\left(\begin{array}{l} \psi_{1} \\ \psi_{2} \end{array}\right)=e^{i \phi}\left(\begin{array}{c} A \cos \theta-i B \sin \theta \\ A \sin \theta+i B \cos \theta \end{array}\right) \tag{12.11}

with ABA \geq|B| and 0θ<π0 \leq \theta<\pi and ϕ\phi is real phase (which is not very relevant relevant to the physics but can be there to make the math look uglier). This represents elliptical polarization with semi-major axis AA at an angle θ\theta with the 1 axis, as in

uθ=(cosθsinθ).(12.12)u_{\theta}=\left(\begin{array}{c} \cos \theta \\ \sin \theta \end{array}\right) . \tag{12.12}

and semi-minor axis BB as shown in Figure 12.4. If BB is positive (negative), the rotation is counterclockwise (clockwise). The physically interesting parameters AA, BB and θ\theta can be found from ψ1\psi_{1} and ψ2\psi_{2} as follows:

A2+B2=ψ12+ψ22,AB=Im(ψ1ψ2).(12.13)\begin{gathered} A^{2}+B^{2}=\left|\psi_{1}\right|^{2}+\left|\psi_{2}\right|^{2} , \\ A B=-\operatorname{Im}\left(\psi_{1} \psi_{2}^{*}\right) . \tag{12.13} \end{gathered}

Thus,

A±B=ψ12+ψ222Im(ψ1ψ2),(12.14)A \pm B=\sqrt{\left|\psi_{1}\right|^{2}+\left|\psi_{2}\right|^{2} \mp 2 \operatorname{Im}\left(\psi_{1} \psi_{2}^{*}\right)} , \tag{12.14}

gives AA and BB. Then θ\theta satisfies

(A2B2)cos2θ=ψ12ψ22,(A2B2)sin2θ=2Re(ψ1ψ2).(12.15)\begin{aligned} &\left(A^{2}-B^{2}\right) \cos 2 \theta=\left|\psi_{1}\right|^{2}-\left|\psi_{2}\right|^{2} , \\ &\left(A^{2}-B^{2}\right) \sin 2 \theta=2 \operatorname{Re}\left(\psi_{1} \psi_{2}^{*}\right) . \tag{12.15} \end{aligned}

Notice that the overall phase factor eiϕe^{i \phi} cancels out in 12.11-12.13.

12.2: Electromagnetic Waves

General Electromagnetic Plane Waves

Figure12-1

We saw in chapters 8 and 9 that an electromagnetic plane wave traveling in the +z+z direction looks like this,

Ex(z,t)=εxei(kzωt),Ey(z,t)=εyei(kzωt),(12.16)E_{x}(z, t)=\varepsilon_{x} e^{i(k z-\omega t)}, \quad E_{y}(z, t)=\varepsilon_{y} e^{i(k z-\omega t)}, \tag{12.16}
Bx(z,t)=βxei(kzωt),By(z,t)=βyei(kzωt),(12.17)B_{x}(z, t)=\beta_{x} e^{i(k z-\omega t)}, \quad B_{y}(z, t)=\beta_{y} e^{i(k z-\omega t)}, \tag{12.17}
Ez(z,t)=Bz(z,t)=0,(12.18)E_{z}(z, t)=B_{z}(z, t)=0, \tag{12.18}

where the β\betas are determined by Maxwell’s equations to be

βy=ncεx,βx=ncεy.(12.19)\beta_{y}=\frac{n}{c} \varepsilon_{x}, \quad \beta_{x}=-\frac{n}{c} \varepsilon_{y} . \tag{12.19}

As usual, we have written the wave with the irreducible time dependence, eiωte^{-i \omega t}. To get the real electric and magnetic fields, we take the real part of 12.14-12.15. Note, in particular, that the constants εj\varepsilon_{j} and βj\beta_{j} for j=xj = x and yy may be complex.

The restriction to motion in the zz direction is not important. Because the physics of Maxwell’s equations is invariant under rotations in three-dimensional space, we can write down the form of a plane wave moving with an arbitrary k\vec{k} vector by extracting the features of 12.14-12.17 that do not depend on the direction. These are:

  1. k\vec{k}, E\vec{E} and B\vec{B} are mutually orthogonal vectors,

    kE=kB=EB=0;(12.20)\vec{k} \cdot \vec{E}=\vec{k} \cdot \vec{B}=\vec{E} \cdot \vec{B}=0 ; \tag{12.20}
  2. B\vec{B} is determined by the cross product

    B=nck^×E=1ωk×E,(12.21)\vec{B}=\frac{n}{c} \hat{k} \times \vec{E}=\frac{1}{\omega} \vec{k} \times \vec{E} , \tag{12.21}

    where k^\hat{k} is a unit vector in the direction of the k\vec{k} vector, the direction of propagation of the wave.

These two conditions imply that a general real electromagnetic plane wave can be written as

E=Re(e(k)eikriωt)B=Re(b(k)eikriωt)(12.22)\begin{aligned} &\vec{E}=\operatorname{Re}\left(\vec{e}(\vec{k}) e^{i \vec{k} \cdot \vec{r}-i \omega t}\right) \\ &\vec{B}=\operatorname{Re}\left(\vec{b}(\vec{k}) e^{i \vec{k} \cdot \vec{r}-i \omega t}\right) \tag{12.22} \end{aligned}

where the vectors, e\vec{e} and b\vec{b}, are complex, in general, satisfying

b(k)=1ωk×e(k)=nck^×e(k) and k^e(k)=0.(12.23)\vec{b}(\vec{k})=\frac{1}{\omega} \vec{k} \times \vec{e}(\vec{k})=\frac{n}{c} \hat{k} \times \vec{e}(\vec{k}) \quad \text { and } \quad \hat{k} \cdot \vec{e}(\vec{k})=0 . \tag{12.23}

There are two things to note about the relations, 12.21.

  1. It is enough to specify the direction of the electric field, e(k)\vec{e}(\vec{k}). The magnetic field is then determined by 12.21. The vector, e\vec{e} is called the “polarization” of the electromagnetic wave.

  2. Because of 12.21, the polarization is perpendicular to k\vec{k}, and thus lives in a two-dimensional vector space.

In the two-dimensional space perpendicular to k\vec{k}, we can choose a basis of real vectors, e^1\hat{e}_{1} and e^2\hat{e}_{2}, where

e^1k^=e^2k^=e^1e^2=0,e^1×e^2=k^.(12.24)\hat{e}_{1} \cdot \hat{k}=\hat{e}_{2} \cdot \hat{k}=\hat{e}_{1} \cdot \hat{e}_{2}=0, \quad \hat{e}_{1} \times \hat{e}_{2}=\hat{k}. \tag{12.24}

For example, for a plane wave traveling in the zz direction, k^=z^\hat{k}=\hat{z}, we could take e1=x^e_{1}=\hat{x} and e2=y^e_{2}=\hat{y}. Then we can write

e(k)=ψ1e^1+ψ2e^2.(12.25)\vec{e}(\vec{k})=\psi_{1} \hat{e}_{1}+\psi_{2} \hat{e}_{2}. \tag{12.25}

The components, ψ1\psi_{1} and ψ2\psi_{2} go into the the two-dimensional vector, 12.3, that describes the polarization state of the electromagnetic wave, just as it describes the polarization state of the string.[1] We can always go back to the components of the electric field using 12.23 and 12.20 and then find the magnetic field using 12.21.

Now the entire discussion of transverse waves on a string from 12.4 to 12.13 can be taken over to describe polarized light. The direction of displacement of the string goes over directly into the direction of the electric field. Thus the animation in program 12-1 applies just as well to the electric field in a polarized wave as to polarization in a string.

Energy and Intensity

The energy density in an electromagnetic field is

E=12(ϵE2+1μB2).(12.26)\mathcal{E}=\frac{1}{2}\left(\epsilon \vec{E}^{2}+\frac{1}{\mu} \vec{B}^{2}\right) . \tag{12.26}

Because the energy density is a nonlinear function of the field strengths, we must use real fields in 12.24. The momentum density is

P=ϵE×B.(12.27)\overrightarrow{\mathcal{P}}=\epsilon \vec{E} \times \vec{B} . \tag{12.27}

The Poynting vector, a measure of energy flow, is

S=c2P.(12.28)\vec{S}=c^{2} \overrightarrow{\mathcal{P}} . \tag{12.28}

These quantities satisfy

tE+S=0.(12.29)\frac{\partial}{\partial t} \mathcal{E}+\vec{\nabla} \cdot \vec{S}=0 . \tag{12.29}

The Poynting vector is useful because it measures the intensity of the wave, the energy per unit time per unit area carried by the electromagnetic wave. The relation, 12.27, then expresses conservation of energy. The sum of the change in the energy density at any point plus the energy flowing away from it is zero.

To see what these quantities look like in terms of the vector, ZZ, let us compute the electric and magnetic fields explicitly using 12.20 and 12.21:

E=Re(e(k)eikriωt)B=Re(b(k)eikriωt)(12.30)\begin{aligned} &\vec{E}=\operatorname{Re}\left(\vec{e}(\vec{k}) e^{i \vec{k} \cdot \vec{r}-i \omega t}\right) \\ &\vec{B}=\operatorname{Re}\left(\vec{b}(\vec{k}) e^{i \vec{k} \cdot \vec{r}-i \omega t}\right) \tag{12.30} \end{aligned}
b(k)=1ωk×e(k)=nck^×e(k) and k^e(k)=0.(12.31)\vec{b}(\vec{k})=\frac{1}{\omega} \vec{k} \times \vec{e}(\vec{k})=\frac{n}{c} \hat{k} \times \vec{e}(\vec{k}) \quad \text { and } \quad \hat{k} \cdot \vec{e}(\vec{k})=0 . \tag{12.31}

The result is

E=A1e^1cos(krωt+ϕ1)+A2e^2cos(krωt+ϕ2),B=μϵ(A1e^2cos(krωt+ϕ1)A2e^1cos(krωt+ϕ2)).(12.32)\begin{gathered} \vec{E}=A_{1} \hat{e}_{1} \cos \left(\vec{k} \cdot \vec{r}-\omega t+\phi_{1}\right)+A_{2} \hat{e}_{2} \cos \left(\vec{k} \cdot \vec{r}-\omega t+\phi_{2}\right) , \\ \vec{B}=\sqrt{\mu \epsilon}\left(A_{1} \hat{e}_{2} \cos \left(\vec{k} \cdot \vec{r}-\omega t+\phi_{1}\right)-A_{2} \hat{e}_{1} \cos \left(\vec{k} \cdot \vec{r}-\omega t+\phi_{2}\right)\right) . \tag{12.32} \end{gathered}

Putting this into 12.24 and 12.26 gives

E=ϵ4π(A12cos2(krωt+ϕ1)+A22cos2(krωt+ϕ2)),(12.33)\mathcal{E}=\frac{\epsilon}{4 \pi}\left(A_{1}^{2} \cos ^{2}\left(\vec{k} \cdot \vec{r}-\omega t+\phi_{1}\right)+A_{2}^{2} \cos ^{2}\left(\vec{k} \cdot \vec{r}-\omega t+\phi_{2}\right)\right), \tag{12.33}
S=k^ϵμc4π(A12cos2(krωt+ϕ1)+A22cos2(krωt+ϕ2)).(12.34)\vec{S}=\hat{k} \sqrt{\frac{\epsilon}{\mu}} \frac{c}{4 \pi}\left(A_{1}^{2} \cos ^{2}\left(\vec{k} \cdot \vec{r}-\omega t+\phi_{1}\right)+A_{2}^{2} \cos ^{2}\left(\vec{k} \cdot \vec{r}-\omega t+\phi_{2}\right)\right) . \tag{12.34}

You can check explicitly that 12.27 is satisfied. Because ω\omega is very large for interesting electromagnetic waves, we are almost always interested in only the time averaged values of E\mathcal{E} and S\vec{S}. These are

E=ϵ8π(A12+A22),(12.35)\langle\mathcal{E}\rangle=\frac{\epsilon}{8 \pi}\left(A_{1}^{2}+A_{2}^{2}\right), \tag{12.35}
S=k^ϵμc8π(A12+A22).(12.36)\langle\vec{S}\rangle=\hat{k} \sqrt{\frac{\epsilon}{\mu}} \frac{c}{8 \pi}\left(A_{1}^{2}+A_{2}^{2}\right) . \tag{12.36}

Note that the time averaged values depend only on the quantity

Z2ψ12+ψ22=A12+A22.(12.37)|Z|^{2} \equiv\left|\psi_{1}\right|^{2}+\left|\psi_{2}\right|^{2}=A_{1}^{2}+A_{2}^{2} . \tag{12.37}

The intensity of the light is proportional to Z2|Z|^{2}.

Circular Polarization and Spin

Although linear polarization is more familiar and perhaps easier to understand, there is a sense in which circular polarization is the more fundamental. The plane electromagnetic wave in the k^\hat{k} direction can be rotated around the k^\hat{k} axis without changing anything but its polarization state. The rotation symmetry of the physics suggests that we ought to be able to find states that behave simply under such a rotation, and just get multiplied by a phase factor. These states are, in fact, the circular polarization states. The action of a rotation by an angle θ\theta about the k^\hat{k} axis on the polarization vector, ZZ, is represented by the matrix

Rθ=(cosθsinθsinθcosθ).(12.38)R_{\theta}=\left(\begin{array}{cc} \cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{array}\right) . \tag{12.38}

For example, RθR_{\theta} acting on u1u_{1}, 12.4, gives uθu_{\theta}, 12.6:

Rθ(10)=(cosθsinθ).(12.39)R_{\theta}\left(\begin{array}{l} 1 \\ 0 \end{array}\right)=\left(\begin{array}{l} \cos \theta \\ \sin \theta \end{array}\right) . \tag{12.39}

But on the left- and right-circularly polarized states,

Rθ(1i)=eiθ(1i),Rθ(1i)=eiθ(1i).(12.40)R_{\theta}\left(\begin{array}{l} 1 \\ i \end{array}\right)=e^{-i \theta}\left(\begin{array}{l} 1 \\ i \end{array}\right), \quad R_{\theta}\left(\begin{array}{c} 1 \\ -i \end{array}\right)=e^{i \theta}\left(\begin{array}{c} 1 \\ -i \end{array}\right) . \tag{12.40}

This is related to the fact that the circularly polarized states carry the maximum angular momentum possible, which in turn is related to the quantum mechanical property of the spin of the photon.


12.3: Wave Plates and Polarizers

One reason that polarization is important is that the polarization state of an electromagnetic wave can be easily manipulated. Two of the most important devices for such manipulation are polarizers and wave plates.

Unpolarized Light

Figure12-2

In any beam of light, at any given point and time, the electric field points in a particular direction. Likewise, because any plane electromagnetic wave with a definite angular frequency can be described by 12.20 and 12.21,

E=Re(e(k)eikriωt)B=Re(b(k)eikriωt)(12.41)\begin{aligned} \vec{E} &=\operatorname{Re}\left(\vec{e}(\vec{k}) e^{i \vec{k} \cdot \vec{r}-i \omega t}\right) \\ \vec{B} &=\operatorname{Re}\left(\vec{b}(\vec{k}) e^{i \vec{k} \cdot \vec{r}-i \omega t}\right) \tag{12.41} \end{aligned}
b(k)=1ωk×e(k)=nck^×e(k) and k^e(k)=0.(12.42)\vec{b}(\vec{k})=\frac{1}{\omega} \vec{k} \times \vec{e}(\vec{k})=\frac{n}{c} \hat{k} \times \vec{e}(\vec{k}) \quad \text { and } \quad \hat{k} \cdot \vec{e}(\vec{k})=0. \tag{12.42}

every plane wave is polarized. However, in an “unpolarized” beam, the light wave consists of a range of angular frequencies with different polarizations. As a result of the interference of the different harmonic components of the wave, the polarization wanders more or less randomly as a function of time and space, and on the average, no particular polarization is picked out. A simple example of what this looks like is animated in program 12-2, where we plot an electric field of the form

Ex(t)=cos(ω1t+ϕ1)+cos(ω2t+ϕ2),Ey(t)=cos(ω3t+ϕ3)+cos(ω4t+ϕ4),(12.43)\begin{aligned} &E_{x}(t)=\cos \left(\omega_{1} t+\phi_{1}\right)+\cos \left(\omega_{2} t+\phi_{2}\right), \\ &E_{y}(t)=\cos \left(\omega_{3} t+\phi_{3}\right)+\cos \left(\omega_{4} t+\phi_{4}\right), \tag{12.43} \end{aligned}

where the phases are random and the frequencies are chosen at random in a small range around a central frequency. You can watch the E\vec{E} field wandering in the xx-yy plane, eventually filling it up. The narrower the range of frequencies in the wave, the more slowly the polarization wanders. In the example in program 12-2, the range of frequencies is of the order of 10% of the central frequency, so the polarization wanders rapidly. But for a beam with a fairly well-defined frequency, the polarization will be nearly constant over many cycles of the wave. The time over which the polarization is approximately constant is called the coherence time of the wave. For a plane wave of definite frequency, the coherence time is infinite.

Polarizers

A “polarizer” is a device that allows light polarized in a particular direction (the “easy transmission axis” of the polarizer) to pass through with very little absorption, but absorbs most of the light polarized in the perpendicular direction. Thus an unpolarized light beam, passing though the polarizer, emerges polarized along the easy axis.

For the transverse oscillations of a string, a polarizer is simply a slit that allows the string to oscillate in one transverse direction but not in the perpendicular direction.

For electromagnetic waves, the most familiar example of a polarizer, Polaroid, was invented by Edwin Land over 50 years ago, partly in experiments done in the attic of the Jefferson Physical Laboratory, where he worked as an undergraduate at Harvard. The idea of polaroid is to make material that conducts electricity (poorly) in one direction, but not in the other. Then the electric field in the conducting direction will be absorbed (the energy going to resistive loss), while the electric field in the nonconductive direction will be unaffected. One way of doing this is to make sheets of polymer (polyvinyl alcohol) stretched (to align the polymer molecules along a preferred axis) and doped with iodine (to allow conduction along the polymer molecules).[2]

Wave Plates

“Wave plates” are optical elements that change the relative phase of the two components of ZZ. Wave plates are possible because there are materials in which the index of refraction depends on the polarization. This property is called “birefringence.” It can happen in various ways.

For example, the transparent polymer material cellophane is made into thin sheets by stretching. Because of the stretching, the polymer strands tend to be oriented along the stretch direction. The dielectric constant in this material depends on the direction of the electric field. It is easier for charges to move along the polymer strands than across them. Thus the dielectric constant is larger for electric fields in the stretch direction.

The same effect may arise because of the inherent structure of a transparent crystal. An example is the naturally occurring mineral, calcite, a crystalline form of calcium carbonate, CaCO3CaCO_{3}. Crystals of calcite have the fascinating property of splitting a beam of unpolarized light into its two polarization states. Birefringence can even be produced mechanically, by stressing a transparent material, squeezing the electronic structure in one direction.

However the birefringence is produced, we can make a wave plate by orienting the material so that the xx and yy directions correspond to different indices of refraction, nxn_{x} and nyn_{y}, and then making a slice of the material in the form of a plate in the xx-yy plane, with some thickness \ell in the zz direction. Now an electromagnetic wave traveling in the zz direction through the plate has different kk values depending on its polarization:

k={nxcω for polarization in the x direction nycω for polarization in the y direction (12.44)k=\left\{\begin{array}{l} \frac{n_{x}}{c} \omega \text { for polarization in the } x \text { direction } \\ \frac{n_{y}}{c} \omega \text { for polarization in the } y \text { direction } \end{array}\right. \tag{12.44}

In particular, the phase difference, between xx and yy polarized light in going through the plate is

Δϕ=nxnycω.(12.45)\Delta \phi=\frac{n_{x}-n_{y}}{c} \omega \ell . \tag{12.45}

Note that in general the phase difference, Δϕ\Delta \phi, depends on the frequency of the light. Even if nxn_{x} and nyn_{y} depend on frequency, it would be a bizarre accident if that dependence canceled the ω\omega dependence from the explicit factor of ω\omega in 12.39.

Initially unpolarized light passing through a pair of crossed polarizers with a wave plate in between.

Figure 12.5:Initially unpolarized light passing through a pair of crossed polarizers with a wave plate in between.

Consider, now, putting such a wave plate between two crossed polarizers, oriented at ±45\pm 45^{\circ}, as shown in Figure 12.5. Without the wave plate, no light would get through because the first polarizer transmits only light polarized at 4545^{\circ}, described by the ZZ vector

Z=(1/21/2)(12.46)Z=\left(\begin{array}{l} 1 / \sqrt{2} \\ 1 / \sqrt{2} \end{array}\right) \tag{12.46}

and the second polarizer absorbs it.

Coming out of the first polarizer, the vector, ZZ, looks like 12.40 for all the frequency components in the white light. But when the wave plate is inserted in between, a frequency dependent phase difference is added, so that the ZZ vector coming out of the wave plate (up to an irrelevant overall phase) looks like

Z=(1/2eiΔϕ/2).(12.47)Z=\left(\begin{array}{c} 1 / \sqrt{2} \\ e^{-i \Delta \phi} / \sqrt{2} \end{array}\right) . \tag{12.47}

For frequencies such that eiΔϕe^{-i \Delta \phi} is −1, the light is polarized in the −4545^{\circ} direction, and gets eiΔϕe^{-i \Delta \phi} through the second polarizer without further attenuation. But for frequencies such that e is 1, the light is still absorbed by the second polarizer. Intermediate frequencies are partially absorbed.

It is this frequency dependence that produces the interesting patterns of color that you see when you put cellophane or a stressed piece of plastic between polarizers.

Matrices

The effects of wave plates and polarizers and the like can be summarized by multiplication of the ZZ vector by 2×2 matrices. For example, a perfect polarizer with an axis at an angle θ\theta from the 1 axis can be represented by

Pθ=(cos2θcosθsinθcosθsinθsin2θ).(12.48)P_{\theta}=\left(\begin{array}{cc} \cos ^{2} \theta & \cos \theta \sin \theta \\ \cos \theta \sin \theta & \sin ^{2} \theta \end{array}\right) . \tag{12.48}

The object PθP_{\theta} is called a “projection operator,” because it projects the vector onto the direction parallel to uθu_{\theta}. It satisfies

PθPθ=Pθ,(12.49)P_{\theta} P_{\theta}=P_{\theta}, \tag{12.49}

as it must, since the first polarizer produces polarized light and the second one transmits it perfectly. PθP_{\theta} acting on a vector transmits the component in the θ\theta direction. This is easiest to visualize if θ=0\theta = 0 or π/2\pi / 2. The matrices

P0=(1000),Pπ/2=(0001),(12.50)P_{0}=\left(\begin{array}{ll} 1 & 0 \\ 0 & 0 \end{array}\right), \quad P_{\pi / 2}=\left(\begin{array}{ll} 0 & 0 \\ 0 & 1 \end{array}\right) , \tag{12.50}

represent polarizers along the 1 and 2 axes respectively.

A wave plate in which the phase difference is π/2\pi / 2 is called a “quarter wave plate.” For a wave plate in which the phase difference is between 0 and π\pi, it is conventional to call the axis with the smaller phase the “fast axis.” A quarter wave plate with fast axis along the 1 axis is represented by

Q0=(100i).(12.51)Q_{0}=\left(\begin{array}{ll} 1 & 0 \\ 0 & i \end{array}\right) . \tag{12.51}

Notice that we can write

Q0=P0+iPπ/2.(12.52)Q_{0}=P_{0}+i P_{\pi / 2} . \tag{12.52}

This should convince you that in general if the fast axis is in the θ\theta direction, the quarter wave plate looks like

Qθ=Pθ+iPθ+π/2.(12.53)Q_{\theta}=P_{\theta}+i P_{\theta+\pi / 2} . \tag{12.53}

The discussion of 12.39 shows that in general, a wave plate will only be a quarter wave plate for light of a definite frequency.

A wave plate in which the phase difference is π\pi is called a “half wave plate.” A half wave plate is obtained by replacing the ii in 12.45-12.47 by -1. Thus,

Hθ=PθPθ+π/2.(12.54)H_{\theta}=P_{\theta}-P_{\theta+\pi / 2} . \tag{12.54}

Notice that

Hθ=QθQθ;(12.55)H_{\theta}=Q_{\theta} Q_{\theta}; \tag{12.55}

two quarter wave plates make a half wave plate.

Producing circularly polarized light.

Figure 12.6:Producing circularly polarized light.

Here are two amusing devices that you can make with these optical elements (or matrices). Consider the combination of first a polarizer at 4545^{\circ} and then a quarter wave plate, as shown in Figure 12.6. By forming the matrix product, Q0Pπ/4Q_{0} P_{\pi / 4}, you can see that this produces counterclockwise circularly polarized light from anything with a component of polarization in the π/4\pi / 4 direction. The argument goes like this. The product is

Q0Pπ/4=(100i)(1/21/21/21/2)=(1/21/2i/2i/2).(12.56)Q_{0} P_{\pi / 4}=\left(\begin{array}{cc} 1 & 0 \\ 0 & i \end{array}\right)\left(\begin{array}{cc} 1 / 2 & 1 / 2 \\ 1 / 2 & 1 / 2 \end{array}\right)=\left(\begin{array}{cc} 1 / 2 & 1 / 2 \\ i / 2 & i / 2 \end{array}\right) . \tag{12.56}

When this acts on an arbitrary vector you get circularly polarization unless the vector is annihilated by Pπ/4P_{\pi / 4}.

Q0Pπ/4(ψ1ψ2)=ψ1+ψ22(1i).(12.57)Q_{0} P_{\pi / 4}\left(\begin{array}{l} \psi_{1} \\ \psi_{2} \end{array}\right)=\frac{\psi_{1}+\psi_{2}}{2}\left(\begin{array}{l} 1 \\ i \end{array}\right) . \tag{12.57}

In the opposite order, Pπ/4Q0P_{\pi / 4} Q_{0} is an analyzer for circularly polarized light. It annihilates counterclockwise light and converts clockwise polarized light to light linearly polarized in the π/4\pi / 4 direction.

Optical Activity

“Optical activity” is a property of many organic and some inorganic compounds. An optically active material rotates the polarization of light without absorbing either component of the polarization. A familiar example of such a material is corn syrup, a thick aqueous solution of sugar that you probably have in your kitchen. If you put a rectangular container of corn syrup between polarizers, as shown in Figure 12.7, and rotate the second polarizer until the intensity of the light getting through is a maximum, you will find that direction of the second polarizer is not the same as that of the first. The plane of the polarization has been rotated by some angle θ\theta. The rotation angle, θ\theta, is proportional to the thickness of the container, the length of the region of syrup that the light goes through.

A rectangular container of corn syrup between polarizers.

Figure 12.7:A rectangular container of corn syrup between polarizers.

Clearly, the optical activity of corn syrup cannot depend on crystal structure, because the stuff is a perfectly uniform liquid, completely invariant under rotations in three-dimensional space. It can have no special axes, or any such thing. Optical activity must work very differently from birefringence.

You can find a clue to the nature of optical activity by considering what it looks like if you look at it in a mirror. If you reflect the system illustrated in Figure 12.7 in the xx-zz plane, by changing the sign of all the yy coordinates, the angle θ\theta changes to −θ\theta. Thus the corn syrup that you see in a mirror must be fundamentally different from the corn syrup in your kitchen. This is not so strange. After all, your right hand looks like a left hand when you look at it in a mirror. The corn syrup must have the same property and have a definite “handedness.” In fact, because of the tetrahedral bonding of the carbon atoms of which they are built, the sugar molecules in the corn syrup can and do have such a handedness.

Because of the handedness of the sugar molecules, the index of refraction of the corn syrup actually depends on the handedness of the light. It is slightly different for left- and right-circularly polarized light. This happens because the E\vec{E} field of a circularly polarized beam twists slightly as it traverses each sugar molecule and sees a slightly different electronic structure depending on the direction of the twist. Then, because the indices of refraction are slightly different, the left- and right-circularly polarized components of the light get different phase factors (kk \ell) in passing through a thickness, \ell, of the syrup.

We can now use our matrix language to see how this leads to optical activity. Up to an irrelevant overall phase, we can choose the phase produced on the left-circularly polarized light to be −θ\theta and that on the right-circularly polarized light to be θ\theta. Then we can represent the action of the syrup on an arbitrary wave by the matrix

eiθP++eiθP,(12.58)e^{-i \theta} P_{+}+e^{i \theta} P_{-} , \tag{12.58}

where P±P_{\pm} are matrices that pick out the left- and right-circularly polarized components, respectively. They satisfy

P±(1±i)=(1±i),P±(1i)=0.(12.59)P_{\pm}\left(\begin{array}{c} 1 \\ \pm i \end{array}\right)=\left(\begin{array}{c} 1 \\ \pm i \end{array}\right), \quad P_{\pm}\left(\begin{array}{c} 1 \\ \mp i \end{array}\right)=0 . \tag{12.59}

You can check that the matrices are

P±=12(1i±i1).(12.60)P_{\pm}=\frac{1}{2}\left(\begin{array}{cc} 1 & \mp i \\ \pm i & 1 \end{array}\right) . \tag{12.60}

Then 12.52 becomes

eiθ12(1ii1)+eiθ12(1ii1)=(cosθsinθsinθcosθ).(12.61)e^{-i \theta} \frac{1}{2}\left(\begin{array}{cc} 1 & -i \\ i & 1 \end{array}\right)+e^{i \theta} \frac{1}{2}\left(\begin{array}{cc} 1 & i \\ -i & 1 \end{array}\right)=\left(\begin{array}{cc} \cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{array}\right) . \tag{12.61}

This is just the rotation matrix RθR_{\theta}, of 12.34! RθR_{\theta} rotates both components of any light by an angle θ\theta.

One might wonder about the reason for the handedness of the sugar molecules. In fact, there are physical processes, the weak interactions that give rise to β\beta-radioactivity, that look different when reflected in a mirror[3] and thus in principle could distinguish between left-handed and right-handed molecules. However, these interactions are most likely irrelevant to the handedness of corn syrup. Probably, the reason is biology rather than physics. Long ago, when the beginnings of life emerged from the primordial soup, purely by accident, the right-handed sugars were used. From then on, the handedness was maintained by the processes of reproduction.

Initially unpolarized light passing through a pair of crossed polarizers.

Figure 12.8:Initially unpolarized light passing through a pair of crossed polarizers.

Crossed Polarizers and Quantum Mechanics

Polarization offers many opportunities to get confused when you think of the light wave in terms of photons. Let us imagine turning down the intensity of the light to the point where one photon at a time is going through the polarizers and consider first the deceptively simple situation of light moving in the zz direction through crossed polarizers in the xx-yy plane. Suppose that the first polarizer transmits light polarized in the xx direction, and the second transmits light polarized in the yy direction. This is deceptively simple because it seems that we can interpret what is going on simply in terms of photons. The situation is depicted in Figure 12.8. This seems simple enough to interpret in terms of photons. The unpolarized light in region II is composed equally of photons polarized in the xx direction and in the yy direction (goes the wrong “classical” argument). Those polarized in the xx direction get through the first polarizer, so half the photons are still around in region IIII, where the intensity is reduced by half. Then none of these get through the second polarizer, so that the intensity in region IIIIII is zero.

But compare this with the apparently similar situation in which the second polarizer transmits light polarized at 4545^{\circ} in the xx-yy plane, as shown in Figure 12.9. Now the wave description tells us that the intensity in region IIIIII is reduced by another factor of 2 from that in region IIII. This is impossible to interpret in terms of classical particles. To see this, it is only necessary to turn down the intensity so that only one photon comes through at a time. Then the first polarizer is OK. As before, if the photon is polarized in the xx direction, it get through. But now what happens at the second polarizer. The photon cannot split up. Either it gets through or it doesn’t. To be consistent with the wave description, in which the intensity is reduced by another factor of two, the transmission at the second polarizer must be a probabilistic event. Half the time the photon gets through. Half the time it is absorbed. There is no way for the

Initially unpolarized light passing through a pair of polarizers at with axes at 45^{\circ}.

Figure 12.9:Initially unpolarized light passing through a pair of polarizers at with axes at 4545^{\circ}.

photon in region IIII to tell whether it is going to make it! It is random. God plays dice.


12.4: Boundary between Dielectrics

Let us return to the infinite plane boundary between two dielectrics that we discussed in chapter 9, but now consider an electromagnetic wave coming in at an arbitrary angle. As in chapter 5, we will assume that the boundary is the z=0z = 0 plane, and that for z<0z < 0 we have dielectric constant ϵ\epsilon, while for z>0z > 0, dielectric constant ϵ\epsilon^{\prime}. We assume μ=1\mu=1 everywhere.

On the general grounds of translation invariance and local interactions discussed in the previous chapter, all the components of the electric and magnetic fields will have the general form

ψ(r,t)eikr+Reikr for z0ψ(r,t)τeikr for z0(12.62)\begin{array}{cl} \psi(r, t) \propto e^{i \vec{k} \cdot \vec{r}}+R e^{i \vec{k} \cdot \vec{r}} & \text { for } z \leq 0 \\ \psi(r, t) \propto \tau e^{i \vec{k}^{\prime} \cdot \vec{r}} & \text { for } z \geq 0 \end{array} \tag{12.62}

where

k~x=kx,kx=kx,(12.63)\tilde{k}_{x}=k_{x}, \quad k_{x}^{\prime}=k_{x} , \tag{12.63}

and

k~z=ω2/v2kx2=kzkz=ω2/v2kx2.(12.64)\begin{gathered} \tilde{k}_{z}=-\sqrt{\omega^{2} / v^{2}-k_{x}^{2}}=-k_{z} \\ k_{z}^{\prime}=\sqrt{\omega^{2} / v^{\prime 2}-k_{x}^{2}} \tag{12.64} \end{gathered} .

Thus Snell’s law is satisfied, with θ\theta and θ\theta^{\prime} defined as shown in Figure 12.10.

ksinθ=ksinθ.k=μϵωc=nωcnsinθ=nsinθ.(12.65)\begin{aligned} &k \cdot \sin \theta=k^{\prime} \sin \theta^{\prime} . \\ &|\vec{k}|=\sqrt{\mu \epsilon} \frac{\omega}{c}=n \frac{\omega}{c} \\ &n \sin \theta=n^{\prime} \sin \theta^{\prime} . \tag{12.65} \end{aligned}
Scattering of plane waves from a plane boundary.

Figure 12.10:Scattering of plane waves from a plane boundary.

The details of the scattering will depend on the polarization. It is clear (by symmetry as usual) that the two cases will be polarization in the xx-zz plane and polarization perpendicular to the xx-zz plane. Of course, we lose nothing by considering these two separately, because of linearity. Any polarization for the incoming wave can be dealt with by forming a linear combination of the parallel and perpendicular solutions.

Polarization Perpendicular to the Scattering Plane

Let us first consider perpendicular polarization. This means that the electric field is in the yy direction (out of the plane of the paper), while the magnetic field is the xx-zz plane:4

Ey(r,t)=Aei(krωt)+RAei(k~rωt) for z0Ey(r,t)=τAei(krωt) for z0Ez=Ex=0(12.66)\begin{gathered} E_{y}(r, t)=A e^{i(\vec{k} \vec{r}-\omega t)}+R_{\perp} A e^{i(\tilde{k} \cdot \vec{r}-\omega t)} \quad \text { for } z \leq 0 \\ E_{y}(r, t)=\tau_{\perp} A e^{i\left(\vec{k}^{\prime} \cdot \vec{r}-\omega t\right)} \quad \text { for } z \geq 0 \\ E_{z}=E_{x}=0 \tag{12.66} \end{gathered}

Using 12.19)

B=nck^×E=1ωk×E,(12.67)\vec{B}=\frac{n}{c} \hat{k} \times \vec{E}=\frac{1}{\omega} \vec{k} \times \vec{E} , \tag{12.67}

we can write

Bx(r,t)=ncAcosθei(krωt)+nccosθRAei(krωt) for z0Bx(r,t)nccosθτAei(krωt) for z0Bz(r,t)=ncsinθAei(krωt)+ncsinθRAei(krωt) for z0Bz(r,t)=ncsinθτAei(krωt) for z0By=0(12.68)\begin{array}{rlr} B_{x}(r, t)=-\frac{n}{c} A \cos \theta e^{i(\vec{k} \cdot \vec{r}-\omega t)}+\frac{n}{c} \cos \theta R_{\perp} A e^{i(\vec{k} \cdot \vec{r}-\omega t)} & \text { for } z \leq 0 \\ B_{x}(r, t)--\frac{n^{\prime}}{c} \cos \theta^{\prime} \tau_{\perp} A e^{i\left(\vec{k}^{\prime} \cdot \vec{r}-\omega t\right)} & & \text { for } z \geq 0 \\ B_{z}(r, t)=\frac{n}{c} \sin \theta A e^{i(\vec{k} \cdot \vec{r}-\omega t)}+\frac{n}{c} \sin \theta R_{\perp} A e^{i(\vec{k} \cdot \vec{r}-\omega t)} & & \text { for } z \leq 0 \\ B_{z}(r, t)=\frac{n^{\prime}}{c} \sin \theta^{\prime} \tau_{\perp} A e^{i\left(\vec{k}^{\prime} \cdot \vec{r}-\omega t\right)} & & \text { for } z \geq 0 \\ B_{y}=0 & \end{array} \tag{12.68}

The system is shown in Figure 12.11. This figure shows the directions of the magnetic fields of the incoming (Bi\vec{B}_{i}), reflected (Br\vec{B}_{r}), and transmitted (Bt\vec{B}_{t}) component waves in scattering of an electromagnetic plane wave polarized parallel to a plane dielectric boundary. The k\vec{k} vectors are shown directly beneath the magnetic fields. The nontrivial boundary conditions are that EyE_{y} and BxB_{x} are continuous (the latter because we have assumed μ=1\mu=1 so there is no sheet of bound current on the boundary). BzB_{z} is also continuous, but that provides no new information. Thus

1+R=τ(12.69)1+R_{\perp}=\tau_{\perp} \tag{12.69}

where

ξ=kzkz.(12.70)\xi_{\perp}=\frac{k_{z}^{\prime}}{k_{z}} . \tag{12.70}

Polarization in the Scattering Plane

Polarization in the xx-zz plane looks like

By(r,t)=Aei(krωt)+RAei(krωt) for z0By(r,t)=τAei(krωt) for z0Bz=Bx=0,(12.71)\begin{array}{cl} B_{y}(r, t)=A e^{i(\vec{k} \cdot \vec{r}-\omega t)}+R_{\|} A e^{i(\vec{k} \cdot \vec{r}-\omega t)} & \text { for } z \leq 0 \\ B_{y}(r, t)=\tau_{\|} A e^{i\left(\vec{k}^{\prime} \cdot \vec{r}-\omega t\right)} & \text { for } z \geq 0 \\ B_{z}=B_{x}=0, & \end{array} \tag{12.71}

where, for convenience, we have defined the reflection and transmission coefficients in terms of the magnetic fields, and

Ex(r,t)=cncosθAei(krωt)cncosθRAei(krωt) for z0Ex(r,t)=cncosθτAei(krωt) for z0Ez(r,t)=cnsinθAei(krωt)cnsinθRAei(k~rωt) for z0Ez(r,t)=cnsinθτAei(krωt) for z0Ez=0.(12.72)\begin{array}{rlr} E_{x}(r, t)=\frac{c}{n} \cos \theta A e^{i(\vec{k} \cdot \vec{r}-\omega t)}-\frac{c}{n} \cos \theta R_{\|} A e^{i(-\overrightarrow{\vec{k}} \cdot \vec{r}-\omega t)} & \text { for } z \leq 0 \\ E_{x}(r, t)=\frac{c}{n^{\prime}} \cos \theta^{\prime} \tau_{\|} A e^{i\left(\vec{k}^{\prime} \cdot \vec{r}-\omega t\right)} & & \text { for } z \geq 0 \\ E_{z}(r, t)=-\frac{c}{n} \sin \theta A e^{i(\vec{k} \cdot \vec{r}-\omega t)}-\frac{c}{n} \sin \theta R_{\|} A e^{i(-\tilde{\vec{k}} \cdot \vec{r}-\omega t)} & & \text { for } z \leq 0 \\ E_{z}(r, t)=-\frac{c}{n^{\prime}} \sin \theta^{\prime} \tau_{\|} A e^{i\left(\vec{k}^{\prime} \cdot \vec{r}-\omega t\right)} & & \text { for } z \geq 0 \\ E_{z}=0 . & \end{array} \tag{12.72}

Now the nontrivial boundary conditions are the continuity of ByB_{y} and ExE_{x}. EzE_{z} is not continuous because a surface bound charge density builds up on the dielectric boundary. The boundary conditions yield

1+R=τ(12.73)1+R_{\|}=\tau_{\|} \tag{12.73}
cosθn(1R)=cosθnτ(12.74)\frac{\cos \theta}{n}\left(1-R_{\|}\right)=\frac{\cos \theta^{\prime}}{n^{\prime}} \tau_{\|} \tag{12.74}

or

τ=21+ξ(12.75)\tau_{\|}=\frac{2}{1+\xi_{\|}} \tag{12.75}
R=1ξ1+ξ(12.76)R_{\|}=\frac{1-\xi \|}{1+\xi_{\|}} \tag{12.76}

where

ξ=cosθ/ncosθ/n=n2kzn2kz.(12.77)\xi_{\|}=\frac{\cos \theta^{\prime} / n^{\prime}}{\cos \theta / n}=\frac{n^{2} k_{z}^{\prime}}{n^{\prime 2} k_{z}} . \tag{12.77}

One of the interesting things about 12.74 is that when

n2kzn2kz1(12.78)\frac{n^{2} k_{z}^{\prime}}{n^{\prime 2} k_{z}}-1 \tag{12.78}

there is no reflection. This condition is satisfied for a special angle of incidence called Brewster’s angle. We can understand the significance of Brewster’s angle as follows:

 from Snell’s law, n2n2sin2θsin2θ(12.79)\text { from Snell's law, } \frac{n^{2}}{n^{\prime 2}}-\frac{\sin ^{2} \theta^{\prime}}{\sin ^{2} \theta} \tag{12.79}
kzkz=kx/kzkx/kz=tanθtanθ(12.80)\frac{k_{z}^{\prime}}{k_{z}}=\frac{k_{x} / k_{z}}{k_{x}^{\prime} / k_{z}^{\prime}}=\frac{\tan \theta}{\tan \theta^{\prime}} \tag{12.80}
n2kzn2kz=sinθcosθsinθcosθ=1.(12.81)\frac{n^{2} k_{z}^{\prime}}{n^{\prime 2} k_{z}}=\frac{\sin \theta^{\prime} \cos \theta^{\prime}}{\sin \theta \cos \theta}=1 . \tag{12.81}

Thus sin2θ=sin2θ\sin 2 \theta=\sin 2 \theta^{\prime}. Because θθ\theta \neq \theta^{\prime} (that would be the trivial situation with no boundary), this means that

θ=π/2θ.(12.82)\theta=\pi / 2-\theta^{\prime} . \tag{12.82}

In other words, Brewster’s angle is defined by the condition that the reflected and transmitted plane waves are perpendicular, as shown in the diagram in Figure 12.12. The relevance of this condition is that the reflected wave can be thought of as being produced by the motion of the charges on the boundary. But if these are moving in a direction perpendicular to the electric field in the would-be reflected wave, then the wave cannot be produced.

Brewster's angel.

Figure 12.12:Brewster’s angel.


4 The quantities, RR_{\perp} and τ\tau_{\perp} in this section and RR_{\|} and τ\tau_{\|} in the next are conventionally called “Fresnel coefficients.” See Hecht, page 97.

12.5: Radiation

In this section, we write down the electric and magnetic fields associated with changing charge and current densities.

Fields of moving charges

Because Maxwell’s equations are partial differential equations, lots of initial conditions or boundary conditions must be specified to determine the solutions. For example, a constant electric field everywhere is a solution to the free-space Maxwell’s equations, and therefore you can add a constant field to any solution and it will still be a solution. Such things must be determined by physical initial conditions or boundary conditions. One set of conditions that is frequently interesting is an analog of the boundary condition at infinity that we discussed for one-dimensional waves. Suppose that you have a universe which is initially stationary, with no electric currents, no magnetic fields, and only electric fields due to stationary charges (which you know how to compute from Physics 15b). At some time, you begin to move charges around in some finite region of space. What are the electric and magnetic fields produced in this way? This question has a relatively simple answer that is a nice intuitive generalization of the relations you learned in 15b for the electric and vector potentials from stationary charge and current distributions. These relations were

ϕ(r)=d3rρ(r)rr(12.83)\phi(\vec{r})=\int d^{3} r^{\prime} \frac{\rho\left(\vec{r}^{\prime}\right)}{\left|\vec{r}-\overrightarrow{r^{\prime}}\right|} \tag{12.83}
A(r)=1cd3rJ(r)rr(12.84)\vec{A}(\vec{r})=\frac{1}{c} \int d^{3} r^{\prime} \frac{\overrightarrow{\mathcal{J}}\left(\vec{r}^{\prime}\right)}{\left|\vec{r}-\overrightarrow{r^{\prime}}\right|} \tag{12.84}

The generalizations are

ϕ(r,t)=d3rρ(r,trr/c)rr(12.85)\phi(\vec{r}, t)=\int d^{3} r^{\prime} \frac{\rho\left(\overrightarrow{r^{\prime}}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right)}{\left|\vec{r}-\vec{r}^{\prime}\right|} \tag{12.85}
A(r,t)=1cd3rJ(r,trr/c)rr(12.86)\vec{A}(\vec{r}, t)=\frac{1}{c} \int d^{3} r^{\prime} \frac{\overrightarrow{\mathcal{J}}\left(\overrightarrow{r^{\prime}}, t-|\vec{r}-\vec{r}| / c\right)}{\left|\vec{r}-\overrightarrow{r^{\prime}}\right|} \tag{12.86}

It is a straightforward, but tedious, exercise in vector calculus to show these satisfy Maxwell’s equations. I am not going to talk about this (I’ll write down the derivation in an appendix for those of you who are interested), but it is worth trying to understand what these relations mean physically. The important physical point that these relations imply is that if the charge and current distributions depend on time, and if they are producing the fields, then what determines what the field is at some point r\vec{r} is the values of the charge and current distributions at earlier times. The farther away the charge is, the earlier the time has to be. That is what the factor of trr/ct-\left|\vec{r}-\vec{r}^{\prime}\right| / c is telling us. The appearance of this factor is a kind of boundary condition at infinity. It is consistent with the relativistic version of the principle of causality. Because information cannot be transferred faster than light, a charge distribution at a space-time point (r,t)\left(\vec{r}^{\prime}, t^{\prime}\right) can effect the fields at the space-time point (r,t)(\vec{r}, t), only if ttt \geq t^{\prime} and

rrttc(12.87)\frac{\left|\vec{r}-\vec{r}^{\prime}\right|}{t-t^{\prime}} \leq c \tag{12.87}

In these relations, 12.82 and 12.83, however, the condition is even stronger — a charge distribution at a space time point (r,t)\left(\vec{r}^{\prime}, t^{\prime}\right) can effect the fields at the space time point (r,t)(\vec{r}, t) only if light can travel directly from (r,t)\left(\vec{r}^{\prime}, t^{\prime}\right) to (r,t)(\vec{r}, t) — that is ifttt \geq t^{\prime} and

 ck (12.88)\text { ck } \tag{12.88}

or

tt=rr/c(12.89)t-t^{\prime}=\left|\vec{r}-\overrightarrow{r^{\prime}}\right| / c \tag{12.89}

or

t=trr/c(12.90)t^{\prime}=t-\left|\vec{r}-\vec{r}^{\prime}\right| / c \tag{12.90}

These are just words. We have not derived this! The real justification of this discussion comes when you check that the relations actually satisfy Maxwell’s equations. That can wait for Physics 153 or 232 (or the appendix if you are in a hurry). However, I hope that this discussion at least makes the result reasonable. In fact you have already seen the result in action in 15b in Purcell’s discussion of the electric field from a charge that starts and stops. Look at the ANIMATIONS - PURCELL - the field from a charge that suddenly accelerates. This is an animation of a famous figure in Purcell’s book. The interesting thing about the animation is the kink in the electric field that propagates out from the acceleration event at the velocity of light — because it is light. Inside the kink, the fields are those of the moving charge. Outside the kind, the fields are those of the stationary charge. The kink — the electromagnetic way — is what connects the two asymptotic regions together. It is also fun to compare with PURCELL2 which illustrates what happens if an initially moving charge stops suddenly.

Now let’s see at what the electric and magnetic fields look like in an important limit. The connection between the potentials and the fields is the following:

E=ϕ1ctA(12.91)\vec{E}=-\vec{\nabla} \phi-\frac{1}{c} \frac{\partial}{\partial t} \vec{A} \tag{12.91}
B=×A(12.92)\vec{B}=\vec{\nabla} \times \vec{A} \tag{12.92}

These relations are completely general. The special limit that I want to consider is one in which the charges and currents are confined to a small region around r=0\vec{r} = 0. Then we will look at the electric and magnetic fields produced by the moving charges far away, for large r|\vec{r}|. It is actually easiest to look at the magnetic field:

B=×A=×1cd3rJ(r,trr/c)rr(12.93)\vec{B}=\vec{\nabla} \times \vec{A}=\vec{\nabla} \times \frac{1}{c} \int d^{3} r^{\prime} \frac{\overrightarrow{\mathcal{J}}\left(\overrightarrow{r^{\prime}}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right)}{\left|\vec{r}-\vec{r}^{\prime}\right|} \tag{12.93}

The point is that the curl (×)(\vec{\nabla} \times) can operate in two different places, either on the 1/rr1 /\left|\vec{r}-\overrightarrow{r^{\prime}}\right| or on the rr/c-\left|\vec{r}-\overrightarrow{r^{\prime}}\right| / c in the time dependence of J\mathcal{J}. The first gives a contribution that drops of like 1/r21 / r^{2} for large rr, just like the magnetic field from a time-independent distribution of currents. But the second gives a contribution that only falls off like 1/r1/r. Thus this contribution dominates for large rr. Explicitly (using the chain rule), it is

B=1c2d3rrrrr2×ddtJ(r,trr/c)(12.94)\vec{B}=-\frac{1}{c^{2}} \int d^{3} r^{\prime} \frac{\vec{r}-\vec{r}^{\prime}}{\left|\vec{r}-\vec{r}^{\prime}\right|^{2}} \times \frac{d}{d t} \overrightarrow{\mathcal{J}}\left(\vec{r}^{\prime}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right) \tag{12.94}
1c21rd3rr^×ddtJ(r,trr/c)(12.95)\rightarrow-\frac{1}{c^{2}} \frac{1}{r} \int d^{3} r^{\prime} \hat{r} \times \frac{d}{d t} \overrightarrow{\mathcal{J}}\left(\overrightarrow{r^{\prime}}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right) \tag{12.95}

where in 12.92 we have dropped a r\overrightarrow{r^{\prime}} in the numerator because this term falls like 1/r21 / r^{2} for large rr.

This is the magnetic field of an electromagnetic wave. Notice that it is perpendicular to the direction of motion (r^)(\hat{r}). The 1/r1/r fall-off is what we expect for an electromagnetic wave, because the energy density goes like the square of the field, and falls off like 1/r21 / r^{2} as the wave spreads.

The electric field can be computed in a similar way, although you also need to use the conservation of electric charge.

tρ+J=0(12.96)\frac{\partial}{\partial t} \rho+\vec{\nabla} \cdot \overrightarrow{\mathcal{J}}=0 \tag{12.96}

As you would expect, the result is that the electric field has the same magnitude as the magnetic field and is perpendicular to both direction of motion and to the magnetic field. The piece that corresponds to a traveling electromagnetic wave can be written as

E1c2d3rrrrr×(rrrr2×ddtJ(r,trr/c))(12.97)\vec{E} \rightarrow-\frac{1}{c^{2}} \int d^{3} r^{\prime} \frac{\vec{r}-\overrightarrow{r^{\prime}}}{\left|\vec{r}-\overrightarrow{r^{\prime}}\right|} \times\left(\frac{\vec{r}-\overrightarrow{r^{\prime}}}{\left|\vec{r}-\overrightarrow{r^{\prime}}\right|^{2}} \times \frac{d}{d t} \overrightarrow{\mathcal{J}}\left(\overrightarrow{r^{\prime}}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right)\right) \tag{12.97}
1c21rd3rr^×(r^×ddtJ(r,trr/c))(12.98)\rightarrow \frac{1}{c^{2}} \frac{1}{r} \int d^{3} r^{\prime} \hat{r} \times\left(\hat{r} \times \frac{d}{d t} \overrightarrow{\mathcal{J}}\left(\overrightarrow{r^{\prime}}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right)\right) \tag{12.98}

To lowest order in 1/r1/r for charges moving with velocities much smaller than cc, we can simplify the electric field in 12.95 by substituting

rrr(12.99)|\vec{r}-\vec{r}| \rightarrow r \tag{12.99}

and write the result as

E(r,t)1c21rd3rr^×(r^×ddtJ(r,tr/c))(12.100)\vec{E}(\vec{r}, t) \approx \frac{1}{c^{2}} \frac{1}{r} \int d^{3} r^{\prime} \hat{r} \times\left(\hat{r} \times \frac{d}{d t} \overrightarrow{\mathcal{J}}\left(\overrightarrow{r^{\prime}}, t-r / c\right)\right) \tag{12.100}

The reason for the restriction to nonrelativistic motion of charges is that if a charged particle is moving at a speed close to the velocity of light, then we cannot neglect its position, r\overrightarrow{r^{\prime}}, when it is moving towards r\vec{r}. To see this, consider the impossible limit in which the charge is moving towards the point r\vec{r} at the speed of light. Then if the charge contributes to the electric field at r\vec{r} at one time, then it also contributes at later times because the particle keeps up with the moving light wave. While v=cv = c is impossible, for vcv \approx c, the r\overrightarrow{r^{\prime}} dependence cannot be ignored because it leads to very rapid time dependence of the potentials, and hence to large fields. What happens is that the contribution of charges moving relativistically to the electric fields in front of them get enhanced by factors of ccv\frac{c}{c-v}. This effect is widely used today to produce intense “light” from particle accelerators — so called synchrotron radiation. You can see this effect in the ANIMATIONS if you make vv close to 1.

A particularly important and instructive case of 12.97 is the nonrelativistic motion of a single charge, QQ, moving along a trajectory R(t)\vec{R}(t). For this system,5

J(r,t)=Qv(t)δ3(rR(t))=QdR(t)dtδ3(rR(t))(12.101)\overrightarrow{\mathcal{J}}(\vec{r}, t)=Q \vec{v}(t) \delta^{3}(\vec{r}-\vec{R}(t))=Q \frac{d \vec{R}(t)}{d t} \delta^{3}(\vec{r}-\vec{R}(t)) \tag{12.101}

Then the integration over d3rd^{3} r^{\prime} in (a)(a) eliminates the δ\delta-function, and the electric field of the outgoing electromagnetic wave is proportional to the acceleration,

E(r,t)1c21rQr^×(r^×a(tr/c))(12.102)\vec{E}(\vec{r}, t) \approx \frac{1}{c^{2}} \frac{1}{r} Q \hat{r} \times(\hat{r} \times \vec{a}(t-r / c)) \tag{12.102}

where

a(t)=d2R(l)dt2(12.103)\vec{a}(t)=\frac{d^{2} \vec{R}(l)}{d t^{2}} \tag{12.103}

All that the cross products with r^\hat{r} do is to pick out minus the component of a(lr/c)\vec{a}(l-r / c) that is perpendicular to r\vec{r}. It follows from the famous “bac-cab” identity,

a×(b×c)=b(ac)c(ab),(12.104)\vec{a} \times(\vec{b} \times \vec{c})=\vec{b}(\vec{a} \cdot \vec{c})-\vec{c}(\vec{a} \cdot \vec{b}), \tag{12.104}

that

E(r,t)1c21rQ(a(tr/c)r^(r^a(tr/c))).(12.105)\vec{E}(\vec{r}, t) \approx-\frac{1}{c^{2}} \frac{1}{r} Q(\vec{a}(t-r / c)-\hat{r}(\hat{r} \cdot \vec{a}(t-r / c))) . \tag{12.105}

This had to happen because the electric field of an electromagnetic wave is perpendicular to its direction of motion. In this case, for large rr, the wave is nearly a plane wave moving in direction r\vec{r}.

Antenna Pattern

Let us do an even more explicit example by considering a charge that oscillates harmonically along the zz axis,

R(t)=z^cosωt.(12.106)\vec{R}(t)=\ell \hat{z} \cos \omega t . \tag{12.106}

so that

a(t)=ω2z^cosωt.(12.107)\vec{a}(t)=-\ell \omega^{2} \hat{z} \cos \omega t . \tag{12.107}
E(r,t)ω2c21rQ(z^r^(r^z^))cos[ω(tr/c)].(12.108)\vec{E}(\vec{r}, t) \approx \frac{\ell \omega^{2}}{c^{2}} \frac{1}{r} Q(\hat{z}-\hat{r}(\hat{r} \cdot \hat{z})) \cos [\omega(t-r / c)] . \tag{12.108}

The vector z^r^(r^z^)\hat{z}-\hat{r}(\hat{r} \cdot \hat{z}) is the component of z^\hat{z} perpendicular to r\vec{r}, as illustrated in Figure 12.13.

Evidently, the magnitude of \hat{z}-\hat{r}(\hat{r} \cdot \hat{z}) is \sin \theta. This means that the intensity of the electromagnetic wave at an angle \theta from the z axis is proportional to \sin ^{2} \theta. The intensity pattern can be conveniently represented in polar coordinates, where we…

Figure 12.13:Evidently, the magnitude of z^r^(r^z^)\hat{z}-\hat{r}(\hat{r} \cdot \hat{z}) is sinθ\sin \theta. This means that the intensity of the electromagnetic wave at an angle θ\theta from the zz axis is proportional to sin2θ\sin ^{2} \theta. The intensity pattern can be conveniently represented in polar coordinates, where we plot the intensity as a function of θ\theta. The result is shown below. This is the “antenna pattern” for the oscillating

dipole in the z direction. It is shown in [Figure 12.14](#fig-12-14). The two lobes of the pattern arise because the field is highest in the x-y plane, for \theta=\pi / 2, and drops to zero as we approach the z axis, \theta = 0 or \theta = \pi.

Figure 12.14:dipole in the zz direction. It is shown in Figure 12.14. The two lobes of the pattern arise because the field is highest in the xx-yy plane, for θ=π/2\theta=\pi / 2, and drops to zero as we approach the zz axis, θ=0\theta = 0 or θ=π\theta = \pi.

Checking Maxwell’s equations

These things are called retarded potentials. This is a confusing name, since there is really nothing special about the potentials themselves. What is special is the assumption of a particular relation between the potentials and the charges and currents — that the fields are being produced entirely by the charges and currents. Here I show that they satisfy Maxwell’s equations. I call this an appendix because you are NOT responsible for knowing the details. I include it for your general education.

Some mathematical things to notice about the solution:

tρ+J=0(12.109)\frac{\partial}{\partial t} \rho+\vec{\nabla} \cdot \overrightarrow{\mathcal{J}}=0 \tag{12.109}

implies

tρ+J=0(12.110)\frac{\partial}{\partial t} \rho+\vec{\nabla} \cdot \overrightarrow{\mathcal{J}}=0 \tag{12.110}

This is called the Lorentz gauge condition.

E=2ϕ1ctA(12.111)\vec{\nabla} \cdot \vec{E}=-\nabla^{2} \phi-\frac{1}{c} \frac{\partial}{\partial t} \vec{\nabla} \cdot \vec{A} \tag{12.111}
=(2+1c22t2)ϕ(12.112)=\left(-\nabla^{2}+\frac{1}{c^{2}} \frac{\partial^{2}}{\partial t^{2}}\right) \phi \tag{12.112}
=d3r(ρ(r,trr/c)(2)1rr(12.113)=\int d^{3} r^{\prime}\left(\rho\left(\overrightarrow{r^{\prime}}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right)\left(-\nabla^{2}\right) \frac{1}{\left|\vec{r}-\overrightarrow{r^{\prime}}\right|}\right. \tag{12.113}
+1rr(2+1c22t2)ρ(r,trr/c)(12.114)+\frac{1}{\left|\vec{r}-\vec{r}^{\prime}\right|}\left(-\nabla^{2}+\frac{1}{c^{2}} \frac{\partial^{2}}{\partial t^{2}}\right) \rho\left(\overrightarrow{r^{\prime}}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right) \tag{12.114}
2(ρ(r,trr/c))(1rr))(12.115)\left.-2\left(\vec{\nabla} \rho\left(\overrightarrow{r^{\prime}}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right)\right) \cdot\left(\vec{\nabla} \frac{1}{\left|\vec{r}-\overrightarrow{r^{\prime}}\right|}\right)\right) \tag{12.115}

The first term is the one we want. It is

=d3rρ(r,trr/c)4πδ3(rr)(12.116)=\int d^{3} r^{\prime} \rho\left(\vec{r}^{\prime}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right) 4 \pi \delta^{3}\left(\vec{r}-\vec{r}^{\prime}\right) \tag{12.116}
=4πρ(r,t)(12.117)=4 \pi \rho(\vec{r}, t) \tag{12.117}

The other two terms cancel because of the special form of the variable trr/ct-\left|\vec{r}-\overrightarrow{r^{\prime}}\right| / c.

1rr(2+1c22t2)ρ(r,trr/c)(12.118)\frac{1}{\left|\vec{r}-\overrightarrow{r^{\prime}}\right|}\left(-\nabla^{2}+\frac{1}{c^{2}} \frac{\partial^{2}}{\partial t^{2}}\right) \rho\left(\vec{r}^{\prime}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right) \tag{12.118}
2(ρ(r,trr/c))(1rr)(12.119)-2\left(\vec{\nabla} \rho\left(\overrightarrow{r^{\prime}}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right)\right) \cdot\left(\vec{\nabla} \frac{1}{\left|\vec{r}-\vec{r}^{\prime}\right|}\right) \tag{12.119}
=1rr1c2ρ¨(r,trr/c)(12.120)=\frac{1}{\left|\vec{r}-\vec{r}^{\prime}\right|} \frac{1}{c^{2}} \ddot{\rho}\left(\vec{r}^{\prime}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right) \tag{12.120}
+1rr1crrrrρ˙(r,trr/c)(12.121)+\frac{1}{\left|\vec{r}-\vec{r}^{\prime}\right|} \frac{1}{c} \vec{\nabla} \cdot \frac{\vec{r}-\vec{r}^{\prime}}{\left|\vec{r}-\vec{r}^{\prime}\right|} \dot{\rho}\left(\overrightarrow{r^{\prime}}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right) \tag{12.121}
21c1rrρ˙(r,trr/c)(12.122)-2 \frac{1}{c} \frac{1}{\left|\vec{r}-\vec{r}^{\prime}\right|} \dot{\rho}\left(\vec{r}^{\prime}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right) \tag{12.122}

where * means differentiation with respect to the time variable:

ρ˙(r,trr/c)tρ(r,t)t=trr/c(12.123)\left.\dot{\rho}\left(\overrightarrow{r^{\prime}}, t-\left|\vec{r}-\vec{r}^{\prime}\right| / c\right) \equiv \frac{\partial}{\partial t^{\prime}} \rho\left(\overrightarrow{r^{\prime}}, t^{\prime}\right)\right|_{t^{\prime}=t-\left|\vec{r}-\vec{r}^{\prime}\right| / c} \tag{12.123}

ii and iiii come from aa (from the 2t2\frac{\partial^{2}}{\partial t^{2}} and 2−\nabla^{2} terms respectively) and iiiiii comes from bb. Now the \vec{\nabla} in iiii gives two terms — acting on ρ˙\dot{\rho} cancels ii and acting on cancels iiiiii. Thus

E=4πρ.(12.124)\vec{\nabla} \cdot \vec{E}=4 \pi \rho . \tag{12.124}

Likewise,

×B1ctE(12.125)\vec{\nabla} \times \vec{B}-\frac{1}{c} \frac{\partial}{\partial t} \vec{E} \tag{12.125}
=×(×A)1ct(ϕ1ctA)(12.126)=\vec{\nabla} \times(\vec{\nabla} \times \vec{A})-\frac{1}{c} \frac{\partial}{\partial t}\left(-\vec{\nabla} \phi-\frac{1}{c} \frac{\partial}{\partial t} \vec{A}\right) \tag{12.126}
=(A)2A(1ctϕ)+1c22t2A=\vec{\nabla}(\vec{\nabla} \cdot \vec{A})-\nabla^{2} \vec{A}-\vec{\nabla}\left(\frac{1}{c} \frac{\partial}{\partial t} \phi\right)+\frac{1}{c^{2}} \frac{\partial^{2}}{\partial t^{2}} \vec{A}

or using the Lorentz gauge condition,

=(2+1c22t2)A=\left(-\nabla^{2}+\frac{1}{c^{2}} \frac{\partial^{2}}{\partial t^{2}}\right) \vec{A}

Form here on, the derivation is the same as for E\vec{\nabla} \cdot \vec{E}, and we find

×B1ctE=4πcJ\vec{\nabla} \times \vec{B}-\frac{1}{c} \frac{\partial}{\partial t} \vec{E}=\frac{4 \pi}{c} \overrightarrow{\mathcal{J}}

QED.


5 This equation makes use of δ-function notation. To a physicist, a δ-function δ(x) is just a function that has area 1 and is so sharply peaked around x = 0 that we don’t care exactly what it looks like. All that matters is the area and where the peak is. The δ3 ~ ¡ ~r − R(t) ¢ in the equation is actually the product of three delta functions, for the x, y and z components, and just tells you that ~r = (x, y, z) = R~(t) = ¡ X(t), Y (t), Z(t) ¢ — that is that the particle is moving along the trajectory R~(t). For a mathematical discussion of the δ-function you can look at http://mathworld.wolfram.com/DeltaFunction.html. But don’t be frightened. It is just a simple device for ignoring small details that we don’t care about. If you translate the integral into words or pictures, it may help.

Problems

Footnotes
  1. This is sometimes called the Jones vector. See Hecht, page 323.

  2. See Sears, Zemansky and Young, page 813.

  3. They violate what is called “parity” symmetry.