Multipole Expansion for Radiation;Vector Spherical Harmonics Michael Dine Department of Physics University of California, Santa Cruz February 2013
We seek a more systematic treatment of the multipole expansion for radiation. The strategy will be to consider three regions: 1 Radiation zone: r λ d. This is the region we have already considered for the dipole radiation. But we will see that there is a deeper connection between the usual multipole moments and the radiation at large distances (which in all cases falls as 1/r). 2 Intermediate zone (static zone): λ r d. Note that time derivatives are of order 1/λ (c = 1), while derivatives with respect to r are of order 1/r, so in this region time derivatives are negligible, and the fields appear static. Here we can do a conventional multipole expansion. 3 Near zone: d r. Here is is more difficult to find simple approximations for the fields.
Our goal is to match the solutions in the intermediate and radiation zones. We will see that in the intermediate zone, because of the static nature of the field, there is a multipole expansion identical to that of electrostatics (where moments are evaluated at each instant). This solution will match onto outgoing spherical waves, all falling as eikr iωt r, but with a sequence of terms suppressed by powers of d/λ. So we have two kinds of expansion going on, a different one in each region.
Let s first review some of the special functions we will need for this analysis. The basic equation which will interest is the Helmholtz equation (similar to the Schrodinger equation): ( 2 + ω 2 )ψ( x, ω) = 0. (1) One can expand the solutions for fixed r, in spherical coordinates: ψ( x, ω) = l,m f l,m (r)y l,m (θ, φ). (2) f l,m obeys: [ d 2 dr 2 + 2 r d dr + k 2 ] l(l + 1) r 2 f l (r) = 0. (3)
The f l (r) s are spherical Bessel functions," f l (r) = g l (kr). The g l s are rather nice functions (you will encounter them in quantum mechanics in studying scattering theory). Typically they are presented as four types, the j l s, n l s, and h l s, spherical Bessel, Neumann, and Hankel functions. They are distinguished by their behaviors at the origin and at. The j l s are regular at the origin, n l s are singular; both behave like sines at. There are actually two types of h s, h (1) l = j l + in l ; h (2) l = j l in l. (4) These behave as e ±ix at, corresponding to outgoing and incoming spherical waves. They are both singular at the origin (a suitable linear combination is not).
j 0 (x) = sin x n 0 (x) = cos(x) x x with relatively simple asymptotics: h (1) eix 0 (x) = ix (5) x 0 : j l (x) x l (2l + 1)!! (2l 1)!! n l (x) x l+1. (6) x : j l (x) 1 x sin(x lπ 2 ) n l(x) 1 x cos(x lπ ). (7) 2 h (1) l l+1 eix ( i) x h (2) l (i) l+1 e ix x (8) With e iωt time dependence, h (1) l e iωt+ikr is an outgoing spherical wave.
Green s function for the Helmholtz equation: ( 2 + k 2 )G( x, x ) = δ( x x ) (9) G( x, x ) = eik x x 4π x x (10) = l,m g l (r, r )Y l,m (θ φ )Y lm (θφ) where g l (r, r ) = ikj l (kr < )h (1) l (kr > ).
Jackson proves this by looking at the equations for the coefficients of each Y lm. But we can see why this must be correct rather simply. First, for r < r, say, for each l the equation for g l is the same equation we studied above, in both r and r. As r 0, the result should be non-singular, so we must take the j l s. For r, as r, we would like outgoing spherical wave behavior, which fixes h (1) l. The argument is symmetric in r and r, which explains the r <, r >.
Now we can fix the coefficients by requiring that as k 0, this should go over to 1 1 2l + 1 l,m r > ( r< r > ) l Y lm (θ, φ )Y lm (θ, φ). (11) Examining the expansions of h,j for small argument, yields the coefficient above: ik (kr <) l ( ) (2l 1)!! (2l + 1)!! (kr > ) l+1 (12) = 1 1 2l + 1 r > ( ) l r< r >
So we are ready to attack the problem. We seek a generalization of the static multipole expansion. We will follow a slight modification of Jackson s strategy. We start by considering the intermediate zone. It turns out to be simpler to work with E and H (I will follow Jackson in using H rather than B). E and H are vector quantities. It is helpful to start by considering the scalar quantities r E and r H. The first, in electrostatics, can be completely determined from the charge distribution; this is the same problem we will encounter in the intermediate zone. r H is similar. We will focus mostly on the first.
Since 2 r E = r 2 E + 2 E = r 2 E (13) and similarly for H, it follows that: { ( 2 + ω 2 r ) H } r E = 0. (14) These are scalar quantities, so they can readily be expanded in spherical harmonics: r E = l,m k l(l + 1)a E (l, m)g l,m (kr)y l,m (Ω). (15) and similarly for r H.
Now for the problems which typically interest us, in which we have a localized source of radiation, of size much smaller than a wavelength, we can rather easily figure out the coefficients a E. Consider, first, the intermediate zone", λ r d. In this regime, the fields are essentially static. You can see this by comparing derivatives: r 1 r 1 c t 1 λ. (16) So we can neglect time derivatives in this regime. But this means that E φ (17)
In this static regime, φ d 3 x ρ( x ) x x. (18) (this also follows from the form of the Green s function of the Helmholtz equation in the limit that we can neglect k). Then r E = r φ = φ r (19)
But we know the form of φ far from a localized charge distribution: φ( x) = 4π q l,m Y l,m (Ω) 2l + 1 r l+1. (20) l,m We can construct φ r, so for small r, a l,m l(l + 1)gl (kr) 4π l + 1 1 2l + 1 r l+1 q l,m. (21)
Now we want to match to the general solution. Away from r = 0, we have g l (kr) = h (1) l (kr) (22) due to the requirement of outgoing spherical wave boundary conditions. Now we match, noting that (eqns. 9.85,9.88 in Jackson) h (1) (2l 1)!! l i x l+1 (23) so a E (l, m) = ck l+2 i(2l + 1)!! ( ) l + 1 1/2 q lm. (24) l
Knowledge of r E determines H and E, through the usual equations for plane waves in free space: H = i k E E = i k H (25) So, e.g., for the TM" ( r H = 0) case, let s solve the full set of equations. Take (a guess) H lm = A lm (r) LY l,m (26) Note first that LA(r) = 0. Using the language of quantum mechanics, this is because A(r) is rotationally invariant, and L is the generator of infinitesimal rotations. But it follows from: iɛ ijk x j k A(r) = iɛ ijk x j x k r A(r) r = 0. (27)
Then it is easy to see that: r E = i k r ( H) = i k r iɛ ijk A lm j L k Y lm (28) 1 k L 2 A lm (r)y l,m = 1 k A lml(l + 1)Y lm. so that the connection of A(r) and the expansion coefficients of r E can be read off immediately (Jackson s 9.122). You can check that this configuration (with the appropriate H, as above) solves the full set of Maxwell s equations; similarly for the "TE" modes.
It is natural to define vector spherical harmonics": X lm (θ, φ) = 1 l(l + 1) LYlm (θ, φ). (29) From the properties of the spherical harmonics, it is easy to see that: X lm X l m dω = δ l,l δ m,m (30) X lm ( r X l m )dω = 0 (31) (the last follows from an integration by parts and use of the commutation relations of L and r). H = lm A lm (r) l(l + 1) X lm. (32)
Following your text, we can derive this more directly. Rather than work with A and φ, it is convenient to work with r E and r H directly. We focus on the latter. Defining a field, E which is divergence free, E = E + i J (33) ωɛ 0 we have the Maxwell equations (in Jacksons 9.160, set Z 0 = 1, and M = 0) H = 0 E = 0 (34) E ik H = i ωɛ 0 J H + ik E = 0.
Now take the curl of the third equation, substitute for H from the fourth, to give: ( 2 + k 2 ) E = i/k ( J). (35) (There is a similar equation for H). So ( 2 + k 2 ) r E = i/k ( L J). (36)
So now we can solve for r E using the Green s function for the Helmholtz equation: r E (x) = 1 e ik x x 4πk x x L J( x )d 3 x. (37) Now L J = i( r ) ( J) (38) can be simplified using our ɛ identity trick.
We need to be careful about the ordering of the derivatives. The expression above is: iɛ ijk ɛ klm (r i j l J m ) (39) = i(r i j i J j r i j j J i ) = ir i 2 J i i r ( J) = i 2 (r i J i ) 2i J i r ( J)
This can be rewritten: L J = i 2 ( r J) i r r (r 2 A) (40)
So using this, we have (using the expansion of the Green s function in spherical Bessel functions, eqns. 9.95-9.97 in Jackson: a E (l, m) = k 2 l(l + 1) j l (kr)y l,m [ 1 k 2 2 ( r J) ic k 1 r ] r (r 2 ρ) d 3 x. (41) Integrating by parts on the 2 term twice, and using the wave equation (without sources) gives a factor of k 2 ; but we can actually drop this term in the case of a localized source, leaving our earlier expression for a E (using the small argument limit of the Bessel functions). Exercise: Verify the matching for the electric dipole and quadrupole moments, i.e. check that these results reproduce our calculation of the dipole and quadrupole radiation, checking the connection between the large and small distance expressions.
Applications Jackson considers two applications of this method. First are cases where low moments vanish. This happens in some atoms and nuclei. Alternatively, it is interesting, for example, for antennae when the wavelength is not too different than the size of the antenna. (Think of FM, 100 MHz). Jackson also (in 9.4) discusses the problem of determining the current in an antenna, and shows that this is actually a boundary value problem, where one must solve simultaneously for the current and the field.