[Click here for a PDF of this post with nicer formatting]
Problem 1. Energy, momentum, etc., of EM waves.
\item Calculate the energy density, energy flux, and momentum density of a plane monochromatic linearly polarized electromagnetic wave.
\item Calculate the values of these quantities averaged over a period.
\item Imagine that a plane monochromatic linearly polarized wave incident on a surface (let the angle between the wave vector and the normal to the surface be ) is completely reflected. Find the pressure that the EM wave exerts on the surface.
\item To plug in some numbers, note that the intensity of sunlight hitting the Earth is about ( the intensity is the average power per unit area transported by the wave). If sunlight strikes a perfect absorber, what is the pressure exerted? What if it strikes a perfect reflector? What fraction of the atmospheric pressure does this amount to?
Part 1. Energy and momentum density.
Because it doesn’t add too much complexity, I’m going to calculate these using the more general elliptically polarized wave solutions. Our vector potential (in the Coulomb gauge , ) has the form
The elliptical polarization case only differs from the linear by allowing to be complex, rather than purely real or purely imaginary. Observe that the Coulomb gauge condition implies
a fact that will kill of terms in a number of places in the following manipulations.
Also observe that for this to be a solution to the wave equation operator
the frequency and wave vector must be related by the condition
For the time and spatial phase let’s write
In the Coulomb gauge, our electric and magnetic fields are
Similar to \S 48 of the text , let’s split into a phase and perpendicular vector components so that
where has a real square
This allows a split into two perpendicular real vectors
where since is real.
Our electric and magnetic fields are now reduced to
or explicitly in terms of and
The special case of interest for this problem, since it only strictly asked for linear polarization, is where and one of or is zero (i.e. is strictly real or strictly imaginary). The case with strictly real, as done in class, is
Now lets calculate the energy density and Poynting vectors. We’ll need a few intermediate results.
Let’s use arrowed vectors for the phasor parts
where we can recover our vector quantities by taking real parts , . Our energy density in terms of these phasors is then
Note that , and (since ). Also , so we have
Now, for the Poynting vector. We have
Reducing the terms we get , and , leaving
Now, the text in \S 47 defines the energy flux as the Poynting vector, and the momentum density as , so we just divide 1.22 by for the momentum density and we are done. For the linearly polarized case (all that was actually asked for, but less cool to calculate), where is real, we have
Part 2. Averaged.
We want to average over one period, the time such that , so the average is
It is clear that this will just kill off the sinusoidal terms, leaving
Part 3. Pressure.
The magnitude of the momentum of light is related to its energy by
and can thus loosely identify the magnitude of the force as
With pressure as the force per area, we could identify
as the instantaneous (directed) pressure on a surface. What is that for linearly polarized light? We have from above for the linear polarized case (where )
If we look at the magnitude of the average pressure from the radiation, we have
Part 4. Sunlight.
With atmospheric pressure at , and the pressure from the light at , we have roughly of pressure from the sunlight being only of the total atmospheric pressure. Wow. Very tiny!
Would it make any difference if the surface is a perfect absorber or a reflector? Consider a ball hitting a wall. If it manages to embed itself in the wall, the wall will have to move a bit to conserve momentum. However, if the ball bounces off twice the momentum has been transferred to the wall. The numbers above would be for perfect absorbtion, so double them for a perfect reflector.
Problem 2. Spherical EM waves.
Suppose you are given:
where and is the unit vector in the -direction. This is a simple example of a spherical wave.
\item Show that obeys all four Maxwell equations in vacuum and find the associated magnetic field.
\item Calculate the Poynting vector. Average over a full cycle to get the intensity vector . Where does it point to? How does it depend on ?
\item Integrate the intensity vector flux through a spherical surface centered at the origin to find the total power radiated.
Part 1. Maxwell equation verification and magnetic field.
Our vacuum Maxwell equations to verify are
We’ll also need the spherical polar forms of the divergence and curl operators, as found in \S 1.4 of 
We can start by verifying the divergence equation for the electric field. Observe that our electric field has only an component, so our divergence is
We have a zero divergence since the component has no dependence (whereas itself does since the unit vector ).
All of the rest of Maxwell’s equations require so we’ll have to first calculate that before progressing further.
A aside on approaches attempted to find
I tried two approaches without success to calculate . First I hoped that I could just integrate to obtain and then take the curl. Doing so gave me a result that had . I hunted for an algebraic error that would account for this, but could not find one.
The second approach that I tried, also without success, was to simply take the cross product . This worked in the monochromatic plane wave case where we had
since one can easily show that . Again, I ended up with a result for that did not have a zero divergence.
Finding with a more systematic approach.
Following  \S 16.2, let’s try a phasor approach, assuming that all the solutions, whatever they are, have all the time dependence in a term.
Let’s write our fields as
Substitution back into Maxwell’s equations thus requires equality in the real parts of
With we can now directly compute the magnetic field phasor
The electric field of this problem can be put into phasor form by noting
which allows for reading off the phasor part directly
Now we can compute the magnetic field phasor . Since we have only a component in our field, the curl will have just and components. This is reasonable since we expect it to be perpendicular to .
Chugging through all the algebra we have
so our magnetic phasor is
Multiplying by and taking real parts gives us the messy magnetic field expression
Since this was constructed directly from , this implicitly verifies one more of Maxwell’s equations, leaving only , and . Neither of these looks particularly fun to verify, however, we can take a small shortcut and use the phasors to verify without the explicit time dependence.
From 2.54 we have for the divergence
Let’s also verify the last of Maxwell’s equations in phasor form. The time dependence is knocked out, and we want to see that taking the curl of the magnetic phasor returns us (scaled) the electric phasor. That is
With only and components in the magnetic phasor we have
Immediately, we see that with no explicit dependence in the coordinates, we have no nor terms in the curl, which is good. Our curl is now just
What we expect is which is
FIXME: Somewhere I must have made a sign error, because these aren’t matching! Have an extra term and the wrong sign on the term.
Part 2. Poynting and intensity.
Our Poynting vector is
which we could calculate from 2.34, and 2.55. However, that looks like it’s going to be a mess to multiply out. Let’s use instead the trick from \S 48 of the course text , and work with the complex quantities directly, noting that we have
Now we can do the Poynting calculation using the simpler relations 2.52, 2.54.
Let’s also write
So our Poynting vector is
Note that our unit vector basis was rotated from , so we have
and plug this into our Poynting expression
Now we have to multiply out our terms. We have
Since this has no real part, there is no average contribution to in the direction. What do we have for the time dependent part
This is non zero, so we have a time dependent contribution that averages out. Moving on
This is non-zero, so the steady state Poynting vector is in the outwards radial direction. The last piece is
Assembling all the results we have
We can read off the intensity directly
Part 3. Find the power.
Through a surface of radius , integration of the intensity vector 2.68 is
Our average power through the surface is therefore
Notes on grading of my solution.
Problem 2 above was the graded portion.
FIXME1: I lost a mark in the spot I expected, where I failed to verify one of the Maxwell equations. I’ll still need to figure out what got messed up there.
What occured to me later, also mentioned in the grading of the solution was that Maxwell’s equations in the space-time domain could have been used to solve for instead of all the momentum space logic (which simplified some things, but probably complicated others).
FIXME2: I lost a mark on 2.68 with a big beside it. I’ll have to read the graded solution to see why.
FIXME3: Lost a mark for the final average power result 2.69. Again, I’ll have to go back and figure out why.
 L.D. Landau and E.M. Lifshitz. The classical theory of fields. Butterworth-Heinemann, 1980.
 D.J. Griffith. Introduction to Electrodynamics. Prentice-Hall, 1981.
 JD Jackson. Classical Electrodynamics Wiley. John Wiley and Sons, 2nd edition, 1975.