Multipole radiation is a theoretical framework for the description of electromagnetic or gravitational radiation from time-dependent distributions of distant sources. These tools are applied to physical phenomena which occur at a variety of length scales - from gravitational waves due to galaxy collisions to gamma radiation resulting from nuclear decay.[1][2][3] Multipole radiation is analyzed using similar multipole expansion techniques that describe fields from static sources, however there are important differences in the details of the analysis because multipole radiation fields behave quite differently from static fields. This article is primarily concerned with electromagnetic multipole radiation, although the treatment of gravitational waves is similar.
Electromagnetic radiation depends on structural details of the source system of electric charge and electric current. Direct analysis can be intractable if the structure is unknown or complicated. Multipole analysis offers a way to separate the radiation into moments of increasing complexity. Since the electromagnetic field depends more heavily on lower-order moments than on higher-order moments, the electromagnetic field can be approximated without knowing the structure in detail.
Properties of multipole radiation
Linearity of moments
Since Maxwell's equations are linear, the electric field and magnetic field depend linearly on source distributions. Linearity allows the fields from various multipole moments to be calculated independently and added together to give the total field of the system. This is the well-known principle of superposition.
Origin dependence of multipole moments
Multipole moments are calculated with respect to a fixed expansion point which is taken to be the origin of a given coordinate system. Translating the origin changes the multipole moments of the system with the exception of the first non-vanishing moment.[4][5] For example, the monopole moment of charge is simply the total charge in the system. Changing the origin will never change this moment. If the monopole moment is zero then the dipole moment of the system will be translation invariant. If both the monopole and dipole moments are zero then the quadrupole moment is translation invariant, and so forth. Because higher-order moments depend on the position of the origin, they cannot be regarded as invariant properties of the system.
Field dependence on distance
The field from a multipole moment depends on both the distance from the origin and the angular orientation of the evaluation point with respect to the coordinate system.[4] In particular, the radial dependence of the electromagnetic field from a stationary
-pole scales as
.[2] That is, the electric field from the electric monopole moment scales as inverse distance squared. Likewise, the electric dipole moment creates a field that scales as inverse distance cubed, and so on. As distance increases, the contribution of high-order moments becomes much smaller than the contribution from low-order moments, so high-order moments can be ignored to simplify calculations.
The radial dependence of radiation waves is different from static fields because these waves carry energy away from the system. Since energy must be conserved, simple geometric analysis shows that the energy density of spherical radiation, radius
, must scale as
. As a spherical wave expands, the fixed energy of the wave must spread out over an expanding sphere of surface area
. Accordingly, every time-dependent multipole moment must contribute radiant energy density that scales as
, regardless of the order of the moment. Hence, high-order moments cannot be discarded as easily as in static case. Even so, the multipole coefficients of a system generally diminish with increasing order, usually as
, so radiation fields can still be approximated by truncating high-order moments.[5]
Time-dependent electromagnetic fields
Sources
Time-dependent source distributions can be expressed using Fourier analysis. This allows separate frequencies to be analyzed independently. Charge density is given by
![{\displaystyle \rho (\mathbf {x} ,t)=\int _{-\infty }^{\infty }d\omega \,{\hat {\rho }}(\mathbf {x} ,\omega )e^{-i\omega t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/11f4ea8dde18692bb252c364188a6b1e03176c5e)
and current density by
[6]
![{\displaystyle \mathbf {J} (\mathbf {x} ,t)=\int _{-\infty }^{\infty }d\omega \,{\hat {\mathbf {J} }}(\mathbf {x} ,\omega )e^{-i\omega t}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5ed592169aa4b8e3c5f8538c6a56745578d13522)
For convenience, only a single angular frequency ω is considered from this point forward; thus
![{\displaystyle \rho (\mathbf {x} ,t)=\rho (\mathbf {x} )e^{-i\omega t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8374e092f8f4ba26a55968365bf20913b92ffe95)
![{\displaystyle \mathbf {J} (\mathbf {x} ,t)=\mathbf {J} (\mathbf {x} )e^{-i\omega t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fa6e79b7129925ccc266c0192a16603c6854cfea)
The
superposition principle may be applied to generalize results for multiple frequencies.
[5] Vector quantities appear in bold. The standard convention of taking the real part of complex quantities to represent physical quantities is used.
The intrinsic angular momentum of elementary particles (see Spin (physics)) may also affect electromagnetic radiation from some source materials. To account for these effects, the intrinsic magnetization of the system
would have to be taken into account. For simplicity however, these effects will be deferred to the discussion of generalized multipole radiation.
Potentials
The source distributions can be integrated to yield the time-dependent electric potential and magnetic potential φ and A respectively. Formulas are expressed in the Lorenz Gauge in SI units.[5][6]
![{\displaystyle \phi (\mathbf {x} ,t)={\frac {1}{4\pi \varepsilon _{0}}}\int d^{3}\mathbf {x'} \int dt'{\frac {\rho (\mathbf {x'} ,t')}{\|\mathbf {x} -\mathbf {x'} \|_{2}}}\delta \left(t'-\left(t-{\frac {\|\mathbf {x} -\mathbf {x'} \|_{2}}{c}}\right)\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2b7a0d66016fc65a99d1431e91eb3a2d6e91cba5)
![{\displaystyle \mathbf {A} (\mathbf {x} ,t)={\frac {\mu _{0}}{4\pi }}\int d^{3}\mathbf {x'} \int dt'{\frac {\mathbf {J} (\mathbf {x'} ,t')}{\|\mathbf {x} -\mathbf {x'} \|_{2}}}\delta \left(t'-\left(t-{\frac {\|\mathbf {x} -\mathbf {x'} \|_{2}}{c}}\right)\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/31c830ef807a15ed8db901a3eefdfd5398d1f64d)
In these formulas c is the speed of light in vacuum,
is the Dirac delta function, and
is the Euclidean distance from the source point x′ to the evaluation point x. Integrating the time-dependent source distributions above yields
![{\displaystyle \phi (\mathbf {x} ,t)={\frac {1}{4\pi \varepsilon _{0}}}e^{-i\omega t}\int d^{3}\mathbf {x'} \rho (\mathbf {x'} ){\frac {e^{ik\|\mathbf {x} -\mathbf {x'} \|_{2}}}{\|\mathbf {x} -\mathbf {x'} \|_{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/22c2760c6dc2c590a6270650965de31561c14b75)
![{\displaystyle \mathbf {A} (\mathbf {x} ,t)={\frac {\mu _{0}}{4\pi }}e^{-i\omega t}\int d^{3}\mathbf {x'} \mathbf {J} (\mathbf {x'} ){\frac {e^{ik\|\mathbf {x} -\mathbf {x'} \|_{2}}}{\|\mathbf {x} -\mathbf {x'} \|_{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d96a6b18ea4bf911a8ba62166246fda1e60e276a)
where
k = ω / c. These formulas provide the basis for analyzing multipole radiation.
Multipole expansion in near field
The near field is the region around a source where the electromagnetic field can be evaluated quasi-statically. If target distance from the multipole origin
is much smaller than the radiation wavelength
, then
. As a result, the exponential can be approximated in this region as:
![{\displaystyle e^{ik\|\mathbf {x} -\mathbf {x'} \|_{2}}=1+{\mathcal {O}}(kr)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b40d71895459de3f0d538e802a055eb86d217190)
See Taylor expansion. By using this approximation, the remaining x′ dependence is the same as it is for a static system, the same analysis applies.[4][5] Essentially, the potentials can be evaluated in the near field at a given instant by simply taking a snapshot of the system and treating it as though it were static - hence it is called quasi-static.[5] See near and far field and multipole expansion. In particular, the inverse distance
is expanded using spherical harmonics which are integrated separately to obtain spherical multipole coefficients.
Multipole expansion in far field: Multipole radiation
At large distances from a high frequency source,
, the following approximations hold:
![{\displaystyle {\frac {1}{\|\mathbf {x} -\mathbf {x'} \|_{2}}}={\frac {1}{r}}+{\mathcal {O}}(1/r^{2})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f6c9c90d7d5b890565fd49d550889f4bc45fb16f)
![{\displaystyle e^{ik\|\mathbf {x} -\mathbf {x'} \|_{2}}=e^{ik(r-\mathbf {n} \cdot \mathbf {x'} +{\mathcal {O}}(1/r))}=e^{ikr-ik\mathbf {n} \cdot \mathbf {x'} }(1+{\mathcal {O}}(1/r))}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7060fa5de8cd5bbad2cdb0fd122cce16c97062f1)
Since only the first-order term in
is significant at large distances, the expansions combine to give
![{\displaystyle {\frac {e^{ik\|\mathbf {x} -\mathbf {x'} \|_{2}}}{\|\mathbf {x} -\mathbf {x'} \|_{2}}}={\frac {e^{ikr}}{r}}\left(1-ik(\mathbf {n} \cdot \mathbf {x'} )+{\frac {(-ik)^{2}}{2}}(\mathbf {n} \cdot \mathbf {x'} )^{2}+\cdots \right)+{\mathcal {O}}(1/r^{2})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2e458864c9d434247d55f533977fbdcbfb75811f)
Each power of
corresponds to a different multipole moment. The first few moments are evaluated directly below.
Electric monopole radiation, nonexistence
The zeroth order term,
, applied to the scalar potential gives
![{\displaystyle \phi _{\text{Electric monopole}}(\mathbf {x} ,t)={\frac {1}{4\pi \varepsilon _{0}}}{\frac {e^{ikr-i\omega t}}{r}}\int d^{3}\mathbf {x'} \rho (\mathbf {x'} )={\frac {e^{ikr-i\omega t}}{4\pi \varepsilon _{0}r}}q}](https://wikimedia.org/api/rest_v1/media/math/render/svg/53f2ab3a769284b9e4446e2ead4e01e786abb94e)
where the total charge
![{\textstyle q=\int d^{3}\mathbf {x'} \rho (\mathbf {x'} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6471915e48fd3743da1015ecbdc879afecee065b)
is the electric monopole moment oscillating at frequency
ω.
Conservation of charge requires
q = 0 since
![{\displaystyle q(t)=\int d^{3}\mathbf {x'} \rho (\mathbf {x'} ,t)=\int d^{3}\mathbf {x'} \rho (\mathbf {x'} )e^{-i\omega t}=qe^{-i\omega t}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ac3352ba706c1a771aab169c9a16dcf2a1fda0f4)
If the system is closed then the total charge cannot fluctuate which means the oscillation amplitude q must be zero. Hence,
. The corresponding fields and radiant power must also be zero.[5]
Electric dipole radiation
Electric dipole potential
Electric dipole radiation can be derived by applying the zeroth-order term to the vector potential.[5]
![{\displaystyle \mathbf {A} _{\text{Electric dipole}}(\mathbf {x} ,t)={\frac {\mu _{0}}{4\pi }}{\frac {e^{ikr-i\omega t}}{r}}\int d^{3}\mathbf {x'} \mathbf {J} (\mathbf {x'} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4506dca01c675cac6695b34753ba62aa998c340e)
Integration by parts yields[7]
![{\displaystyle \int d^{3}\mathbf {x'} \mathbf {J} (\mathbf {x'} )=-\int d^{3}\mathbf {x'} \mathbf {x'} ({\boldsymbol {\nabla }}\cdot \mathbf {J} (\mathbf {x'} )).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/22b2177ccb27e509f58b77e292bfb4491ff7defb)
and the charge continuity equation shows
![{\displaystyle {\frac {\partial \rho (\mathbf {x} ,t)}{\partial t}}+{\boldsymbol {\nabla }}\cdot \mathbf {J} (\mathbf {x} ,t)=\left(-i\omega \rho (\mathbf {x} )+{\boldsymbol {\nabla }}\cdot \mathbf {J} (\mathbf {x} )\right)e^{-i\omega t}=0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/245ec1869fee20fe383573705a8be295eabc7669)
It follows that
![{\displaystyle \mathbf {A} _{\text{Electric dipole}}(\mathbf {x} ,t)={\frac {-i\omega \mu _{0}}{4\pi }}{\frac {e^{ikr-i\omega t}}{r}}\int d^{3}\mathbf {x'} \mathbf {x'} \rho (\mathbf {x'} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3c4010189687b3f87e78740771472371b67279b5)
Similar results can be obtained by applying the first-order term,
to the scalar potential. The amplitude of the electric dipole moment of the system is
, which allows the potentials to be expressed as
![{\displaystyle \phi _{\text{Electric dipole}}(\mathbf {x} ,t)={\frac {-ik}{4\pi \varepsilon _{0}}}{\frac {e^{ikr-i\omega t}}{r}}\mathbf {n} \cdot \mathbf {p} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/0c605b1ec6f62e13751e291c08a00f1a70ad8cc6)
![{\displaystyle \mathbf {A} _{\text{Electric dipole}}(\mathbf {x} ,t)={\frac {-i\omega \mu _{0}}{4\pi }}{\frac {e^{ikr-i\omega t}}{r}}\mathbf {p} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/227c6a6cc0b5df1dc3d5a75ea7dcd8bccd8e63b4)
Electric dipole fields
Once the time-dependent potentials are understood, the time-dependent electric field and magnetic field can be calculated in the usual way. Namely,
![{\displaystyle \mathbf {E} (\mathbf {x} ,t)=-{\boldsymbol {\nabla }}\phi (\mathbf {x} ,t)-{\frac {\partial \mathbf {A} (\mathbf {x} ,t)}{\partial t}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4481055fba7ee79781da8e999fbd498b57dd697b)
![{\displaystyle \mathbf {B} (\mathbf {x} ,t)={\boldsymbol {\nabla }}\times \mathbf {A} (\mathbf {x} ,t),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/83398850c31bc548ed21337b11c64f62e2570b1a)
or, in a source-free region of space, the relationship between the magnetic field and the electric field can be used to obtain
![{\displaystyle \mathbf {H} (\mathbf {x} ,t)={\frac {1}{\mu _{0}}}{\boldsymbol {\nabla }}\times \mathbf {A} (\mathbf {x} ,t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eb48c6b39ccd15a69cdda9fd5b64e157b44117d8)
![{\displaystyle \mathbf {E} (\mathbf {x} ,t)={\frac {iZ_{0}}{k}}{\boldsymbol {\nabla }}\times \mathbf {H} (\mathbf {x} ,t)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/49db9adb1a73cf8e4bff680e25e14b03c3f26c96)
where
is the impedance of free space. The electric and magnetic fields that correspond to the potentials above are
![{\displaystyle \mathbf {H} _{\text{Electric dipole}}(\mathbf {x} ,t)={\frac {ck^{2}}{4\pi }}(\mathbf {n} \times \mathbf {p} ){\frac {e^{ikr-i\omega t}}{r}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ded60286f3cdb92827cce65dda5fe815e3b9ed70)
![{\displaystyle \mathbf {E} _{\text{Electric dipole}}(\mathbf {x} ,t)=Z_{0}(\mathbf {H} _{\text{Electric dipole}}\times \mathbf {n} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0e307cec88ea45befeadf45d29360ce26a0d9f63)
which is consistent with spherical radiation waves.[5]
Pure electric dipole power
The power density, energy per unit area per unit time, is expressed by the Poynting vector
. It follows that the time averaged power density per unit solid angle is given by
![{\displaystyle {\frac {dP(\mathbf {x} )}{d\Omega }}={\frac {r^{2}}{2}}\Re (\mathbf {n} \cdot \mathbf {E} \times \mathbf {H} ).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/01b8f39013673aceff13fc95740d7a9c585537c5)
The dot product with
extracts the emission magnitude and the factor of 1/2 comes from averaging over time. As explained above, the
cancels the radial dependence of radiation energy density. Application to a pure electric dipole gives
![{\displaystyle {\frac {dP_{\text{Electric dipole}}(\mathbf {x} )}{d\Omega }}={\frac {c^{2}Z_{0}}{32\pi ^{2}}}k^{4}\|\mathbf {p} \|_{2}^{2}\sin ^{2}\theta }](https://wikimedia.org/api/rest_v1/media/math/render/svg/c6b3fcd6b37c0e3e6dbf3e06c23b66e81a6c2939)
where θ is measured with respect to
.[5] Integration over a sphere yields the total power radiated:
![{\displaystyle P_{\text{Electric dipole}}={\frac {c^{2}Z_{0}}{12\pi }}k^{4}\|\mathbf {p} \|_{2}^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/815fe49615e8f85a7a5b5bf7b3d9c195dddacd36)
Magnetic dipole radiation
Magnetic dipole potential
The first-order term,
, applied to the vector potential gives magnetic dipole radiation and electric quadrupole radiation.[5]
![{\displaystyle \mathbf {A} _{\text{Magnetic dipole / Electric quadrupole}}(\mathbf {x} ,t)={\frac {\mu _{0}}{4\pi }}{\frac {e^{ikr-i\omega t}}{r}}(-ik)\int d^{3}\mathbf {x'} (\mathbf {n} \cdot \mathbf {x'} )\mathbf {J} (\mathbf {x'} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c974064c94342cd567c8280cbbacc64b9a0e8515)
The integrand can be separated into symmetric and anti-symmetric parts in J and x′
![{\displaystyle (\mathbf {n} \cdot \mathbf {x'} )\mathbf {J} (\mathbf {x'} )={\frac {1}{2}}\left((\mathbf {n} \cdot \mathbf {x'} )\mathbf {J} (\mathbf {x'} )+(\mathbf {n} \cdot \mathbf {J} (\mathbf {x'} ))\mathbf {x'} \right)+{\frac {1}{2}}(\mathbf {x'} \times \mathbf {J} (\mathbf {x'} ))\times \mathbf {n} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/c0b8f815de0bded6ef275584b960740517390f08)
The second term contains the effective magnetization due to the current
and integration gives the magnetic dipole moment.
![{\displaystyle \int d^{3}\mathbf {x'} \mathbf {M} _{\text{effective}}(\mathbf {x'} )=\mathbf {m} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/ece7013c98d7c7209617060d2b9c74d8ea556c95)
![{\displaystyle \mathbf {A} _{\text{Magnetic dipole}}(\mathbf {x} ,t)={\frac {-ik\mu _{0}}{4\pi }}{\frac {e^{ikr-i\omega t}}{r}}\mathbf {m} \times \mathbf {n} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/f4102bb7b816cf3f5cd2bcfc226c15a76998e56a)
Notice that
has a similar form to
. That means the magnetic field from a magnetic dipole behaves similarly to the electric field from an electric dipole. Likewise, the electric field from a magnetic dipole behaves like the magnetic field from an electric dipole. Taking the transformations
![{\displaystyle \mathbf {E} _{\text{Electric dipole}}\rightarrow Z_{0}\mathbf {H} _{\text{Magnetic dipole}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f4a8e41ffefeb1ad768da810ad9afacc6ba83f08)
![{\displaystyle \mathbf {H} _{\text{Electric dipole}}\rightarrow {\frac {-1}{Z_{0}}}\mathbf {E} _{\text{Magnetic dipole}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/453878627760b6a99f0e785845c570eb374be1be)
![{\displaystyle \mathbf {p} \rightarrow \mathbf {m} /c}](https://wikimedia.org/api/rest_v1/media/math/render/svg/80bd044284711b6557c4239831f4a6db6ab110a2)
on previous results yields magnetic dipole results.[5]
Magnetic dipole fields
![{\displaystyle \mathbf {E} _{\text{Magnetic dipole}}(\mathbf {x} ,t)={\frac {-k^{2}Z_{0}}{4\pi }}(\mathbf {n} \times \mathbf {m} ){\frac {e^{ikr-i\omega t}}{r}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/82f9cd0e7a77357e9b504fb740205e138bba6534)
![{\displaystyle \mathbf {H} _{\text{Magnetic dipole}}(\mathbf {x} ,t)={\frac {-1}{Z_{0}}}(\mathbf {E} _{\text{Magnetic dipole}}\times \mathbf {n} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4eb34d9a5e563e0463551bd559eaf43651349bc4)
[5]
Pure magnetic dipole power
The average power radiated per unit solid angle by a magnetic dipole is
![{\displaystyle {\frac {dP_{\text{Magnetic dipole}}(\mathbf {x} )}{d\Omega }}={\frac {Z_{0}}{32\pi ^{2}}}k^{4}\|\mathbf {m} \|_{2}^{2}\sin ^{2}\theta }](https://wikimedia.org/api/rest_v1/media/math/render/svg/3cb154868553fea009343d9dfa6ed0be18cd5f39)
where θ is measured with respect to the magnetic dipole
. The total power radiated is:[5]
![{\displaystyle P_{\text{Magnetic dipole}}={\frac {Z_{0}}{12\pi }}k^{4}\|\mathbf {m} \|_{2}^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/04f4045ce38de1e4f18f7f6f8b7fdd1ccee0416b)
Electric quadrupole radiation
Electric quadrupole potential
The symmetric portion of the integrand from the previous section can be resolved by applying integration by parts and the charge continuity equation as was done for electric dipole radiation.
![{\displaystyle {\frac {1}{2}}\int d^{3}\mathbf {x} \left((\mathbf {n} \cdot \mathbf {x'} )\mathbf {J} (\mathbf {x'} )+(\mathbf {n} \cdot \mathbf {J} (\mathbf {x'} ))\mathbf {x'} \right)={\frac {-i\omega }{2}}\int d^{3}\mathbf {x'} \mathbf {x'} (\mathbf {n} \cdot \mathbf {x'} )\rho (\mathbf {x'} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8fd2d6b13b14e7adfb052d1cc0ae067127f206f7)
![{\displaystyle \mathbf {A} _{\text{Electric quadrupole}}(\mathbf {x} ,t)={\frac {-k\omega \mu _{0}}{8\pi }}{\frac {e^{ikr-i\omega t}}{r}}\int d^{3}\mathbf {x'} \mathbf {x'} (\mathbf {n} \cdot \mathbf {x'} )\rho (\mathbf {x'} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/383f6e650adbb3f2563f115bc977c659dd5058e5)
This corresponds to the traceless electric quadrupole moment tensor
. Contracting the second index with the normal vector
allows the vector potential to be expressed as[5]
![{\displaystyle \mathbf {A} _{\text{Electric quadrupole}}(\mathbf {x} ,t)={\frac {-k\omega \mu _{0}}{8\pi }}{\frac {e^{ikr-i\omega t}}{r}}{\frac {1}{3}}\mathbf {Q(n)} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/0aa481e7f97a2dd9e0c6d2c11705fbebf77172fd)
Electric quadrupole fields
The resulting magnetic and electric fields are:[5]
![{\displaystyle \mathbf {H} _{\text{Electric quadrupole}}(\mathbf {x} ,t)={\frac {-ick^{3}}{24\pi }}{\frac {e^{ikr-i\omega t}}{r}}\mathbf {n} \times \mathbf {Q(n)} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/731ec8ea6c033182aa721cdbe8e47b8e41eea72e)
![{\displaystyle \mathbf {E} _{\text{Electric quadrupole}}(\mathbf {x} ,t)=Z_{0}(\mathbf {H} _{\text{Electric quadrupole}}\times \mathbf {n} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8a16f350ed01744e6dbe376a96063add55cb51df)
Pure electric quadrupole power
The average power radiated per unit solid angle by an electric quadrupole is
![{\displaystyle {\frac {dP_{\text{Electric quadrupole}}(\mathbf {x} )}{d\Omega }}={\frac {c^{2}Z_{0}}{1152\pi ^{2}}}k^{6}\|(\mathbf {n} \times \mathbf {Q(n)} )\times \mathbf {n} \|_{2}^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c7af6fd6c0e9849a4fbdc2bafa31f4e6584ae80a)
where θ is measured with respect to the magnetic dipole
. The total power radiated is:[5]
![{\displaystyle P_{\text{Electric quadrupole}}={\frac {c^{2}Z_{0}}{1440\pi }}k^{6}\sum _{\alpha ,\beta }Q_{\alpha \beta }^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc81460f3d0a464f6dcb7dc048540ee8d56fbbb3)
Generalized multipole radiation
As the multipole moment of a source distribution increases, the direct calculations employed so far become too cumbersome to continue. Analysis of higher moments requires more general theoretical machinery. Just as before, a single source frequency
is considered. Hence the charge, current, and intrinsic magnetization densities are given by
![{\displaystyle \rho (\mathbf {x} ,t)=\rho (\mathbf {x} )e^{-i\omega t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8374e092f8f4ba26a55968365bf20913b92ffe95)
![{\displaystyle \mathbf {J} (\mathbf {x} ,t)=\mathbf {J} (\mathbf {x} )e^{-i\omega t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fa6e79b7129925ccc266c0192a16603c6854cfea)
![{\displaystyle \mathbf {M} (\mathbf {x} ,t)=\mathbf {M} (\mathbf {x} )e^{-i\omega t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dc37cec51e5e880ec9b07018ea3ba4ae9f438ea4)
respectively. The resulting electric and magnetic fields share the same time-dependence as the sources.
![{\displaystyle \mathbf {E} (\mathbf {x} ,t)=\mathbf {E} (\mathbf {x} )e^{-i\omega t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1f6ee86fe5f5fc9b2c2071d090953ef0ca0eb0af)
![{\displaystyle \mathbf {H} (\mathbf {x} ,t)=\mathbf {H} (\mathbf {x} )e^{-i\omega t}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6a95838c2b1914607d62c6072d9c469b7f321f7e)
Using these definitions and the continuity equation allows Maxwell's equations to be written as
![{\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {E} (\mathbf {x} )=-{\frac {iZ_{0}}{k}}{\boldsymbol {\nabla }}\cdot \mathbf {J} (\mathbf {x} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9a70e356a96810b3b297f12ba3c1ceadc1a4f75d)
![{\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {H} (\mathbf {x} )=-{\boldsymbol {\nabla }}\cdot \mathbf {M} (\mathbf {x} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4574b1ed791d1bc565d2de32b8d0014e5ddcb3c6)
![{\displaystyle {\boldsymbol {\nabla }}\times \mathbf {E} (\mathbf {x} )=ikZ_{0}\left(\mathbf {H} (\mathbf {x} )+\mathbf {M} (\mathbf {x} )\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/26632067ea9b7dcbbda1c574547f9070a86eaf20)
![{\displaystyle {\boldsymbol {\nabla }}\times \mathbf {H} (\mathbf {x} )=-{\frac {ik}{Z_{0}}}\mathbf {E} (\mathbf {x} )+\mathbf {J} (\mathbf {x} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e654e1f172955ec3d702270e8ebfdfc7c4fcddb)
These equations can be combined by taking the curl of the last equations and applying the identity
. This gives the vector forms of the non-homogeneous Helmholtz equation.
![{\displaystyle (\nabla ^{2}+k^{2})\mathbf {E} (\mathbf {x} )=-\left[ikZ_{0}\mathbf {J} (\mathbf {x} )+ikZ_{0}{\boldsymbol {\nabla }}\times \mathbf {M} (\mathbf {x} )+{\frac {iZ_{0}}{k}}{\boldsymbol {\nabla }}({\boldsymbol {\nabla }}\cdot \mathbf {J} (\mathbf {x} ))\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/887e287c3a6934872d333a21bf6dc8855301bd00)
![{\displaystyle (\nabla ^{2}+k^{2})\mathbf {H} (\mathbf {x} )=-\left[k^{2}\mathbf {M} (\mathbf {x} )+{\boldsymbol {\nabla }}\times \mathbf {J} (\mathbf {x} )+{\boldsymbol {\nabla }}({\boldsymbol {\nabla }}\cdot \mathbf {M} (\mathbf {x} ))\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5e361c12fa36f7edc47f64714cf5d7b5bb8ad8c7)
Solutions of the wave equation
The homogeneous wave equations that describes electromagnetic radiation with frequency
in a source-free region have the form.
![{\displaystyle ({\boldsymbol {\nabla }}^{2}+k^{2}){\boldsymbol {\Psi }}(\mathbf {x} )=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e671603e56188491699c04fbc748dc0012ede69a)
The wave function
can be expressed as a sum of vector spherical harmonics
![{\displaystyle {\boldsymbol {\Psi }}(\mathbf {x} )=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }f_{\ell m}(kr)\mathbf {X} _{\ell m}(\theta ,\phi )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8449f4cc8fd3cd10abfc1662d3aa8d7d8982c2f8)
![{\displaystyle f_{\ell m}(kr)=A_{\ell m}^{(1)}h_{\ell }^{(1)}(kr)+A_{\ell m}^{(2)}h_{\ell }^{(2)}(kr)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/33c886180081bd096bf6d06404f325c06d2ddcb1)
Where
are the normalized vector spherical harmonics and
and
are spherical Hankel functions. See spherical Bessel functions. The differential operator
is the angular momentum operator with the property
. The coefficients
and
correspond to expanding and contracting waves respectively. So
for radiation. To determine the other coefficients, the Green's function for the wave equation is applied. If the source equation is
![{\displaystyle ({\boldsymbol {\nabla }}^{2}+k^{2}){\boldsymbol {\Psi }}(\mathbf {x} )=-\mathbf {V} (\mathbf {x} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3dcb93af0e1f73058db8f1ac446d04985a2d67f4)
then the solution is:
![{\displaystyle \Psi _{\alpha }(\mathbf {x} )=\sum _{\beta }\int d^{3}\mathbf {x'} G_{\alpha \beta }(\mathbf {x} ,\mathbf {x'} )V_{\beta }(\mathbf {x'} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b3d47232fc97d61f6a9da8beb7a0d9b92974011e)
The Green function can be expressed in vector spherical harmonics.
![{\displaystyle G_{\alpha \beta }(\mathbf {x} ,\mathbf {x'} )=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }ikh_{\ell }^{(1)}(kr)j_{\ell }(kr')X_{\ell m\alpha }(\theta ,\phi )X_{\ell m\beta }^{*}(\theta ',\phi ')}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3cb86acb2f43e3fee4cffee1518626ca6c3935b7)
Note that
is a differential operator that acts on the source function
. Thus, the solution to the wave equation is:
![{\displaystyle {\boldsymbol {\Psi }}(\mathbf {x} )=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }{\frac {ik}{\sqrt {\ell (\ell +1)}}}h_{\ell }^{(1)}(kr)\mathbf {X} _{\ell m}(\theta ,\phi )\int d^{3}\mathbf {x'} j_{\ell }(kr')Y_{\ell m}^{*}(\theta ',\phi ')\mathbf {L'} \cdot \mathbf {V} (\mathbf {x'} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d3b01202380a3630ae2a453574e5b9c2355af5ce)
Electric multipole fields
Applying the above solution to the electric multipole wave equation
![{\displaystyle (\nabla ^{2}+k^{2})\mathbf {H} (\mathbf {x} )=-\left[k^{2}\mathbf {M} (\mathbf {x} )+{\boldsymbol {\nabla }}\times \mathbf {J} (\mathbf {x} )+{\boldsymbol {\nabla }}({\boldsymbol {\nabla }}\cdot \mathbf {M} (\mathbf {x} ))\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5e361c12fa36f7edc47f64714cf5d7b5bb8ad8c7)
gives the solution for the magnetic field:[5]
![{\displaystyle \mathbf {H} ^{(E)}(\mathbf {x} )=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }a_{\ell m}^{(E)}h_{\ell }^{(1)}(kr)\mathbf {X} _{\ell m}(\theta ,\phi )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/49da3099bf2bdc4ccdd9077ae811c7b5c281c5e9)
![{\displaystyle a_{\ell m}^{(E)}={\frac {ik}{\sqrt {\ell (\ell +1)}}}\int d^{3}\mathbf {x'} j_{\ell }(kr')Y_{\ell m}^{*}(\theta ',\phi ')\mathbf {L'} \cdot \left[k^{2}\mathbf {M} (\mathbf {x'} )+{\boldsymbol {\nabla '}}\times \mathbf {J} (\mathbf {x'} )+{\boldsymbol {\nabla '}}({\boldsymbol {\nabla '}}\cdot \mathbf {M} (\mathbf {x'} ))\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3a9cf50f26c24ea4020f50f8b8e648c1b2dbf658)
The electric field is:
![{\displaystyle \mathbf {E} ^{(E)}(\mathbf {x} )={\frac {iZ_{0}}{k}}{\boldsymbol {\nabla }}\times \mathbf {H} ^{(E)}(\mathbf {x} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0186508b070dec19b68b5373300895ab54321e38)
The formula can be simplified by applying the identities
![{\displaystyle \mathbf {L} \cdot \mathbf {V} (\mathbf {x} )=i{\boldsymbol {\nabla }}\cdot (\mathbf {x} \times \mathbf {V} (\mathbf {x} ))}](https://wikimedia.org/api/rest_v1/media/math/render/svg/474976dcd536c492fdd88eeb7b2e933905d224a3)
![{\displaystyle \mathbf {L} \cdot ({\boldsymbol {\nabla }}\times \mathbf {V} (\mathbf {x} ))=i\nabla ^{2}(\mathbf {x} \cdot \mathbf {V} (\mathbf {x} ))-{\frac {i\partial }{r\partial r}}(r^{2}{\boldsymbol {\nabla }}\cdot \mathbf {V} (\mathbf {x} ))}](https://wikimedia.org/api/rest_v1/media/math/render/svg/023ec13fa7f0278d0401236935a77ceef5ec53d1)
![{\displaystyle \mathbf {L} \cdot {\boldsymbol {\nabla }}s(\mathbf {x} )=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6a79decc1de10b1eb1247b8b69f8c452d279376a)
to the integrand, which results in[5]
![{\displaystyle a_{\ell m}^{(E)}={\frac {-ik^{2}}{\sqrt {\ell (\ell +1)}}}\int d^{3}\mathbf {x'} j_{\ell }(kr')Y_{\ell m}^{*}(\theta ',\phi ')\left[-ik{\boldsymbol {\nabla }}\cdot (\mathbf {x'} \times \mathbf {M} (\mathbf {x'} ))-{\frac {i}{k}}\nabla ^{2}(\mathbf {x'} \cdot \mathbf {J} (\mathbf {x'} ))-{\frac {c\partial }{r'\partial r'}}(r'^{2}\rho (\mathbf {x'} ))\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/42ed83e63d8325be3500428fdaaf6df95e5243e9)
Green's theorem and integration by parts manipulates the formula into
![{\displaystyle a_{\ell m}^{(E)}={\frac {-ik^{2}}{\sqrt {\ell (\ell +1)}}}\int d^{3}\mathbf {x'} j_{\ell }(kr')Y_{\ell m}^{*}(\theta ',\phi ')\left[-ik{\boldsymbol {\nabla }}\cdot (\mathbf {x'} \times \mathbf {M} (\mathbf {x'} ))+ik\mathbf {x'} \cdot \mathbf {J} (\mathbf {x'} )\right]+cY_{\ell m}^{*}(\theta ',\phi ')\rho (\mathbf {x'} ){\frac {\partial }{\partial r'}}(r'j_{\ell }(kr'))}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c665faa316752b3fbb85c5c5467941c1cd034002)
The spherical bessel function
can also be simplified by assuming that the radiation length scale is much larger than the source length scale, which is true for most antennas.
![{\displaystyle j_{\ell }(kr')={\frac {(kr')^{\ell }}{(2\ell +1)!!}}+{\mathcal {O}}((kr')^{\ell +2})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9aa0c438e78e72cc8c08adacbe58c65a8b1c8a75)
Retaining only the lowest order terms results in the simplified form for the electric multipole coefficients:[5]
![{\displaystyle a_{\ell m}^{(E)}={\frac {-ick^{\ell +2}}{(2\ell +1)!!}}\left({\frac {\ell +1}{\ell }}\right)^{1/2}[Q_{\ell m}+Q_{\ell m}']}](https://wikimedia.org/api/rest_v1/media/math/render/svg/efb0d33262fc8a8755d8be4324b064eccdfcea92)
![{\displaystyle Q_{\ell m}=\int d^{3}\mathbf {x'} r'^{\ell }Y_{\ell m}^{*}(\theta ',\phi ')\rho (\mathbf {x'} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/575e4fc5270d9f291dc2c970b91829c7a2efa5d4)
![{\displaystyle Q_{\ell m}'=-{\frac {ik}{c(\ell +1)}}\int d^{3}\mathbf {x'} r'^{\ell }Y_{\ell m}^{*}(\theta ',\phi '){\boldsymbol {\nabla }}\cdot (\mathbf {x'} \times \mathbf {M} (\mathbf {x'} ))}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9b650b529b4c3ee3bef5ed950914b4e031d43518)
is the same as the electric multipole moment in the static case if it were applied to the static charge distribution
whereas
corresponds to an induced electric multipole moment from the intrinsic magnetization of the source material.
Magnetic multipole fields
Applying the above solution to the magnetic multipole wave equation
![{\displaystyle (\nabla ^{2}+k^{2})\mathbf {E} (\mathbf {x} )=-\left[ikZ_{0}\mathbf {J} (\mathbf {x} )+ikZ_{0}{\boldsymbol {\nabla }}\times \mathbf {M} (\mathbf {x} ))+{\frac {iZ_{0}}{k}}{\boldsymbol {\nabla }}({\boldsymbol {\nabla }}\cdot \mathbf {J} (\mathbf {x} ))\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4992f4deb02e653f13aa005c5014a693176cc159)
gives the solution for the electric field:[5]
![{\displaystyle \mathbf {E} ^{(M)}(\mathbf {x} )=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }a_{\ell m}^{(M)}h_{\ell }^{(1)}(kr)\mathbf {X} _{\ell m}(\theta ,\phi )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9d8312ef8aaf9ba2d31e0539178345d6bd371973)
![{\displaystyle a_{\ell m}^{(M)}={\frac {ik}{\sqrt {\ell (\ell +1)}}}\int d^{3}\mathbf {x'} j_{\ell }(kr')Y_{\ell m}^{*}(\theta ',\phi ')\mathbf {L'} \cdot \left[ikZ_{0}\mathbf {J} (\mathbf {x} )+ikZ_{0}{\boldsymbol {\nabla }}\times \mathbf {M} (\mathbf {x} ))+{\frac {iZ_{0}}{k}}{\boldsymbol {\nabla }}({\boldsymbol {\nabla }}\cdot \mathbf {J} (\mathbf {x} ))\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/402a1495ac68b40a17a0623aae75852df24fd3bb)
The magnetic field is:
![{\displaystyle \mathbf {H} ^{(M)}(\mathbf {x} )=-{\frac {i}{kZ_{0}}}{\boldsymbol {\nabla }}\times \mathbf {E} ^{(M)}(\mathbf {x} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4dce67bdd10d3988edb211c7aa2f7129dee37151)
As before, the forumula simplifies to:
![{\displaystyle a_{\ell m}^{(M)}={\frac {-ik^{2}}{\sqrt {\ell (\ell +1)}}}\int d^{3}\mathbf {x'} j_{\ell }(kr')Y_{\ell m}^{*}(\theta ',\phi ')\left[{\boldsymbol {\nabla }}\cdot (\mathbf {x'} \times \mathbf {J} (\mathbf {x'} ))-k^{2}\mathbf {x'} \cdot \mathbf {M} (\mathbf {x'} )\right]+Y_{\ell m}^{*}(\theta ',\phi '){\boldsymbol {\nabla }}\cdot \mathbf {M} (\mathbf {x'} ){\frac {\partial }{\partial r'}}(r'j_{\ell }(kr'))}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c6479a704aaa68a696144d5a7e576e4f968758b8)
Retaining only the lowest order terms results in the simplified form for the magnetic multipole coefficients:[5]
![{\displaystyle a_{\ell m}^{(M)}={\frac {-ik^{\ell +2}}{(2\ell +1)!!}}\left({\frac {\ell +1}{\ell }}\right)^{1/2}[M_{\ell m}+M_{\ell m}']}](https://wikimedia.org/api/rest_v1/media/math/render/svg/518c7e76ab207a925394e2db86d7044e1d0ca13a)
![{\displaystyle M_{\ell m}={\frac {1}{\ell +1}}\int d^{3}\mathbf {x'} r'^{\ell }Y_{\ell m}^{*}(\theta ',\phi '){\boldsymbol {\nabla }}\cdot (\mathbf {x'} \times \mathbf {J} (\mathbf {x'} ))}](https://wikimedia.org/api/rest_v1/media/math/render/svg/df0533bd4ad38fab1f9131ea195d343ee983e108)
![{\displaystyle M_{\ell m}'=\int d^{3}\mathbf {x'} r'^{\ell }Y_{\ell m}^{*}(\theta ',\phi '){\boldsymbol {\nabla }}\cdot \mathbf {M} (\mathbf {x'} )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9a3b0e1c35cbe9005a69a01822e8b050118b0467)
is the magnetic multipole moment from the effective magnetization
while
corresponds to the intrinsic magnetization
.
General solution
The electric and magnetic multipole fields combine to give the total fields:[5]
![{\displaystyle \mathbf {E} (\mathbf {x} ,t)=\Re \left(\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }\left[a_{\ell m}^{(M)}h_{\ell }^{(1)}(kr)\mathbf {X} _{\ell m}(\theta ,\phi )+{\frac {iZ_{0}}{k}}a_{\ell m}^{(E)}{\boldsymbol {\nabla }}\times (h_{\ell }^{(1)}(kr)\mathbf {X} _{\ell m}(\theta ,\phi ))\right]e^{-i\omega t}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/14fdb7c26f5c1b653412668dc85c8d8113015e11)
![{\displaystyle \mathbf {H} (\mathbf {x} ,t)=\Re \left(\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }\left[a_{\ell m}^{(E)}h_{\ell }^{(1)}(kr)\mathbf {X} _{\ell m}(\theta ,\phi )-{\frac {i}{kZ_{0}}}a_{\ell m}^{(M)}{\boldsymbol {\nabla }}\times (h_{\ell }^{(1)}(kr)\mathbf {X} _{\ell m}(\theta ,\phi ))\right]e^{-i\omega t}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/23801ec759d86edf96be060dddf94bf9357a82c4)
Note that the radial function
can be simplified in the far field limit
.
![{\displaystyle h_{\ell }^{(1)}(kr)=(-i)^{\ell +1}{\frac {e^{ikr}}{kr}}+{\mathcal {O}}(1/r^{2})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/38aa253bc2b6426d188a895b1efb881732008f8f)
Thus the radial dependence of radiation is recovered.
See also
References
[1]
[2]
[3]
[4]
[5]
[6]