Radiation from Linear Centre-Fed Antenna

As an illustration of the use of a multipole expansion for a source whose dimensions are comparable to a wavelength, consider the radiation from a linear centre-fed antenna. We assume that the antenna runs along the $z$ -axis, and extends from $z=-d/2$ to $z=d/2$ . The current flowing along the antenna vanishes at the end points, and is an even function of $z$ . Thus, we can write

$$I(z, t) = I(\vert z\vert)\, {\rm e}^{-{\rm i}\,\omega\,t},$$ (1547)

where $I(d/2)=0$ . Because the current flows radially, ${\bf r}\times{\bf j} = 0$ . Furthermore, there is no intrinsic magnetization. Thus, according to Equations (1547)-(1548), all of the magnetic multipole coefficients, $a_M(l,m)$ , vanish. In order to calculate the electric multipole coefficients, $a_E(l,m)$ , we need expressions for the charge and current densities. In spherical polar coordinates, the current density ${\bf j}$ can be written in the form

$${\bf j}({\bf r}) = \frac{I(r)}{2\pi \,r^{\,2}} \,\left[\delta(\cos\theta-1)-\delta(\cos\theta+1)\right]{\bf e}_r,$$ (1548)

for $r<d/2$ , where the delta functions cause the current to flow only upwards and downwards along the $z$ -axis. From the continuity equation (1522), the charge density is given by

$$\rho({\bf r}) = \frac{1}{{\rm i}\,\omega} \frac{d I(r)}{dr} \left[\frac{\delta(\cos\theta-1) -\delta(\cos\theta+1)}{2\pi\, r^{\,2}}\right],$$ (1549)

for $r<d/2$ .

The above expressions for ${\bf j}$ and $\rho$ can be substituted into Equation (1543) to give

$$a_E(l,m) =\frac{\mu_0\, c\,k^{\,2}}{2\pi\sqrt{l\,(l+1)}} \int_0^{d/2}\left... ...(k\,r)\,I(r) -\frac{1}{k} \frac{dI(r)}{dr} \frac{d[r\,j_l(k\,r)]}{dr}\right\}dr$$

$$\oint Y_{lm}^{\,\ast}\left[\delta(\cos\theta-1) -\delta(\cos\theta+1)\right]d{\mit\Omega}.$$ (1550)

The angular integral has the value

$$\oint Y_{lm}^{\,\ast}\left[\delta(\cos\theta-1) -\delta(\cos\thet... ...right]d{\mit\Omega} = 2\pi\,\delta_{m0}\, \left[Y_{l0}(0) - Y_{l0}(\pi)\right],$$ (1551)

showing that only $m=0$ multipoles are generated. This is hardly surprising, given the cylindrical symmetry of the antenna. The $m=0$ spherical harmonics are even (odd) about $\theta=\pi/2$ for $l$ even (odd). Hence, the only nonvanishing multipoles have $l$ odd, and the angular integral reduces to

$$\oint Y_{lm}^{\,\ast}\left[\delta(\cos\theta-1) -\delta(\cos\theta+1)\right] d\Omega = \sqrt{4\pi\,(2\,l+1)}.$$ (1552)

After some slight rearrangement, Equation (1552) can be written

$$a_E(l,0) =\frac{\mu_0 \,c\,k}{2\pi} \left[\frac{4\pi\, (2\,l+1)} ... ... +r j_l(k\,r) \left(\frac{d^{\,2} I}{dr^{\,2}} + k^{\,2}\,I\right)\right\}\,dr.$$ (1553)

In order to evaluate the above integral, we need to specify the current $I(z)$ along the antenna. In the absence of radiation, the sinusoidal time variation at frequency $\omega$ implies a sinusoidal space variation with wavenumber $k=\omega/c$ . However, the emission of radiation generally modifies the current distribution. The correct current $I(z)$ can only be found be solving a complicated boundary value problem. For the sake of simplicity, we assume that $I(z)$ is a known function: specifically,

$$I(z) = I\,\sin(k\,d/2 - k\,\vert z\vert),$$ (1554)

for $z<d/2$ , where $I$ is the peak current. With a sinusoidal current, the second term in curly brackets in Equation (1555) vanishes. The first term is a perfect differential. Consequently, Equations (1555) and (1556) yield

$$a_E(l, 0) = \frac{\mu_0\, c\, I}{\pi \,d}\left[\frac{4\pi\, (2\,l+1)} {l\,(l+1)}\right]^{1/2} \left(\frac{k\,d}{2}\right)^2\,j_l(k\,d/2),$$ (1555)

for $l$ odd.

Let us consider the special cases of a half-wave antenna (i.e., $k\, d=\pi$ , so that the length of the antenna is half a wavelength) and a full-wave antenna (i.e., $k\,d=2\pi$ ). For these two values of $k\,d$ , the $l=1$ electric multipole coefficient is tabulated in Table 3, along with the relative values for $l= 3$ and $5$ . It is clear, from the table, that the coefficients decrease rapidly in magnitude as $l$ increases, and that higher $l$ coefficients are more important the larger the source dimensions. However, even for a full-wave antenna, it is generally sufficient to retain only the $l=1$ and $l= 3$ coefficients when calculating the angular distribution of the radiation. It is certainly adequate to keep only these two harmonics when calculating the total radiated power (which depends on the sum of the squares of the coefficients).

Table 3: The first few electric multipole coefficients for a half-wave and a full-wave antenna.

$k\,d$	$a_E(1,0)$	$a_E(3,0)/a_E(1,0)$	$a_E(5,0)/a_E(1,0)$
$\pi$	$4\sqrt{6\pi}\,(\mu_0 \,c\,I/4\pi \,d)$	$4.95\times 10^{-2}$	$1.02\times 10^{-3}$
$2\pi$	$4\pi\sqrt{6\pi}\,(\mu_0\, c\,I/4\pi\, d)$	$3.25\times 10^{-1}$	$3.09\times 10^{-2}$

In the radiation zone, the multipole fields (1479)-(1480) reduce to

$$c\,{\bf B} \simeq \frac{{\rm e}^{\,{\rm i}\,(k\,r -\omega \,t)}}{k\,r} \sum_... ...+1} \left[ a_E(l,m)\,{\bf X}_{lm}+ a_M(l,m) \,{\bf n}\times{\bf X}_{lm}\right],$$ (1556)

$${\bf E} \simeq c\,{\bf B}\times{\bf n},$$ (1557)

where use has been made of the asymptotic form (1432). The time-averaged power radiated per unit solid angle is given by

$$\frac{dP}{d{\mit \Omega}} = \frac{{\rm Re}\,( {\bf n}\cdot{\bf E}\times{\bf B}^{\,\ast})\,r^{\,2}} {2\,\mu_0},$$ (1558)

or

$$\frac{dP}{d{\mit \Omega}} = \frac{1}{2\,\mu_0\, c \, k^{\,2}}\lef... ...,m)\,{\bf X}_{lm} +a_M(l,m)\,{\bf n}\times{\bf X}_{lm}\right]\right\vert^{\,2}.$$ (1559)

Retaining only the $l=1$ and $l= 3$ electric multipole coefficients, the angular distribution of the radiation from the antenna takes the form

$$\frac{dP}{d{\mit \Omega}} = \frac{\vert a_E(l,0)\vert^{\,2}}{4\,\... ..._{1,0} - \frac{a_E(3,0)}{\sqrt{6}\,a_E(1,0)}\,{\bf L} Y_{3,0}\right\vert^{\,2},$$ (1560)

where use has been made of Equation (1476). The various factors in the absolute square are

$$\vert{\bf L} Y_{1,0}\vert^{\,2} = \frac{3}{4\pi}\,\sin^2\theta,$$ (1561)

$$\vert{\bf L} Y_{3,0}\vert^{\,2} = \frac{63}{16\pi}\,\sin^2\theta \,(5\cos^2\theta -1)^2,$$ (1562)

$$({\bf L} Y_{1,0})^{\,\ast}\cdot ({\bf L} Y_{3,0}) =\frac{3\sqrt{21}} {8\pi}\,\sin^2\theta\,(5\cos^2\theta-1).$$ (1563)

With these angular factors, Equation (1561) becomes

$$\frac{dP}{d{\mit \Omega}} = \lambda\,\frac{3\,\mu_0 \,c\, I^{\,2}... ...qrt{\frac{7}{8}}\frac{a_E(3,0)}{a_E(1,0)} \,(5\cos^2\theta-1)\right\vert^{\,2},$$ (1564)

where $\lambda$ equals 1 for a half-wave antenna, and $\pi^2/4$ for a full-wave antenna. The coefficient in front of $(5\,\cos^2\theta -1)$ is $0.0463$ and $0.304$ for the half-wave and full-wave antenna, respectively. It turns out that the radiation pattern obtained from the two-term multipole expansion specified in the previous equation is almost indistinguishable from the exact result for the case of a half-wave antenna. For the case of a full-wave antenna, the two-term expansion yields a radiation pattern that differs from the exact result by less than 5 percent.

The total power radiated by the antenna is

$$P = \frac{1}{2\,\mu_0\, c\,k^{\,2}} \sum_{l\,{\rm odd}} \vert a_E(l,0)\vert^{\,2},$$ (1565)

where use has been made of Equation (1500). It is evident from Table 3 that a two-term multipole expansion gives an accurate expression for the power radiated by both a half-wave and a full-wave antenna. In fact, a one-term multipole expansion gives a fairly accurate result for the case of a half-wave antenna.

It is clear, from the previous analysis, that the multipole expansion converges rapidly when the source dimensions are of order the wavelength of the radiation. It is also clear that if the source dimensions are much less than the wavelength then the multipole expansion is likely to be completely dominated by the term corresponding to the lowest value of $l$ .