In this problem we first investigate the radiation properties of a very short and thin filament of current. Although this filament of current as an antenna is not physically realizable, as a mathematical model it will help us to understand the general behavior of antennas. Here we also introduce some of the fundamental antenna parameters.

## Problem Statement

A very thin and short antenna which works at frequency \(\omega\) has the the following current distribution.

\[\mathbf{J}\left(\mathbf{r},\omega\right)=\mathbf{J_0} \, \delta (\mathbf{r})=J_0 \, \delta (\mathbf{r}) \mathbf{\hat{z}}\]

where \(J_0\) is the amplitude of the current distribution, \(\delta (\mathbf{r})\) is delta function, and \(\mathbf{\hat{z}}\) is the unit vector in z direction of the Cartesian coordinate . Assume that the antenna (current distribution) is placed in a homogeneous environment with permitivity and permeability of \(\varepsilon\) and \(\mu\).

1. Find the generated electric and magnetic fields \(\big(\mathbf{E}\left(\mathbf{r},\omega\right), \mathbf{H}\left(\mathbf{r},\omega\right)\big)\) by the antenna.

2. Sketch the far filed radiation pattern of the antenna in transverse and vertical planes (e.g. x-y and x-z planes).

3. Calculate the half power beamwidth (HW) of the antenna in x-z plane and its directivity.

## Hints

1. First find vector potential \(\mathbf{A}\left(\mathbf{r},\omega\right)\) by solving the equation below.

\[\nabla^2\mathbf{A}\left(\mathbf{r},\omega\right)+\omega^2\mu\,\varepsilon\,\mathbf{A}\left(\mathbf{r},\omega\right)=-\mu\,\mathbf{J}\left(\mathbf{r},\omega\right)\]

Try to use the symmetries in the problem to simplify this equation. Then, find \(\mathbf{E}\left(\mathbf{r},\omega\right)\) and \(\mathbf{H}\left(\mathbf{r},\omega\right)\) using the following equations.

\[\mathbf{E}\left(\mathbf{r},\omega\right)=i \omega\mathbf{A}\left(\mathbf{r},\omega\right)-\frac{\nabla\big(\nabla\cdot\mathbf{A}\left(\mathbf{r},\omega\right)\big)}{i\omega\mu\,\varepsilon}\]

\[\mathbf{H}\left(\mathbf{r},\omega\right)=\frac{1}{\mu} \nabla\times\mathbf{A}\left(\mathbf{r},\omega\right)\]

2. Approximate \(\mathbf{E}\left(\mathbf{r},\omega\right)\) and \(\mathbf{H}\left(\mathbf{r},\omega\right)\) for \(r \gg \lambda \) to find far fields (e.g. \(\mathbf{E_{ff}}\left(\mathbf{r},\omega\right)\) and \(\mathbf{H_{ff}}\left(\mathbf{r},\omega\right)\)). Then, using far fields calculate the Poynting vector. Finally, sketch the radial component of the Poynting vector versus \(\varphi\) (azimuth angle) in x-y plane, and versus \(\theta\) (polar angle) in x-z plane.

3. Half power beamwidth (HW), in x-z plane, is defined as the width of the angular interval in which the radiation intensity is more than half of the maximum radiation intensity.

Directivity is defined as below.

\[D(\theta,\varphi)=\frac{\text{radiation intensity}}{\text{average of radiation intensity over all directions}}\]

## Show Solution

1. The first step to solve this problem is to find the vector potential. Therefore, we start with the equation below.

\begin{equation}

\nabla^2\mathbf{A}\left(\mathbf{r},\omega\right)+\omega^2\mu\,\varepsilon\,\mathbf{A}\left(\mathbf{r},\omega\right)=-\mu\,J_0 \, \delta (\mathbf{r}) \mathbf{\hat{z}} \tag{1}\label{AJ}

\end{equation}

We can split this equation into three equations, each in one direction of the Cartesian coordinate (e.g. \(\mathbf{\hat{x}, \hat{y}, \text{and}\, \hat{z}}\)). If you do that, you will see that for \(A_x\left(\mathbf{r},\omega\right)\) and \(A_y\left(\mathbf{r},\omega\right)\) there is no source term (right hand side of the their governing equations is zero). Therefore, \(A_x\left(\mathbf{r},\omega\right)\) and \(A_y\left(\mathbf{r},\omega\right)\) are zero.

\begin{equation}

\Rightarrow \mathbf{A}\left(\mathbf{r},\omega\right)=A_z\left(\mathbf{r},\omega\right)\mathbf{\hat{z}} \tag{2}

\end{equation}

Using this fact, equation [\ref{AJ}] can be simplified as below.

\begin{equation}

\nabla^2 A_z\left(\mathbf{r},\omega\right)+\omega^2\mu\,\varepsilon\,A_z\left(\mathbf{r},\omega\right)=-\mu\,J_0 \, \delta (\mathbf{r}) \tag{3}\label{AzJ}

\end{equation}

Since the source term (right hand side of equation [\ref{AzJ}]) has a spherical symmetry (since it does not depend on \(\theta\) and \(\varphi\)), \(A_z\left(\mathbf{r},\omega\right)\) is also spherically symmetric (there is not any preferred direction for \(A_z\left(\mathbf{r},\omega\right)\) ). Thus, \(A_z\left(\mathbf{r},\omega\right)\) is only a function of radial distance from the source point, i.e. \(A_z\left(\mathbf{r},\omega\right)=A_z\left(r,\omega\right)\). Due to this spherical symmetry we use spherical coordinates to expand equation [\ref{AzJ}] as shown below.

\begin{equation}

\frac{1}{r^2}\frac{\partial}{\partial r}\left(r^2\,\frac{\partial A_z\left(r,\omega\right)}{\partial r} \right)+\omega^2\mu\,\varepsilon\,A_z\left(r,\omega\right)=-\mu\,J_0 \, \delta (\mathbf{r}) \tag{4}\label{ArJ}

\end{equation}

The right hand side of equation [\ref{ArJ}] is only non-zero in the center of coordinate. So, our strategy to tackle equation [\ref{ArJ}] is to solve it for all space excluding the center, and then satisfying the source condition at center (we will see the second part later). Therefore, we solve the following equation first.

\begin{equation}

\frac{1}{r^2}\frac{\partial}{\partial r}\left(r^2\,\frac{\partial A_z\left(r,\omega\right)}{\partial r} \right)+\omega^2\mu\,\varepsilon\,A_z\left(r,\omega\right)=0\,, \quad (r\neq 0) \tag{5}\label{Ar0}

\end{equation}

To make equation [\ref{Ar0}] a little bit simpler, without loss of generality we define \(A_z\left(r,\omega\right)\) as below.

\begin{equation}

A_z\left(r,\omega\right)= \frac{\phi(r,\omega)}{r}\tag{6}\label{APhi}

\end{equation}

By substituting equation [\ref{APhi}] in equation [\ref{Ar0}] we obtain

\[\frac{\partial^2 \phi(r,\omega)}{\partial r^2}+\omega^2\mu\,\varepsilon\,\phi(r,\omega)=0\,, \quad (r\neq 0)\]

\[\Rightarrow \phi(r,\omega)=C_1\,e^{i k r}+C_2\,e^{-i k r}\]

which results in,

\[A_z(r,\omega)=C_1\,\frac{ e^{i k r}}{r}+C_2\,\frac{ e^{-i k r}}{r}\]

where, \(k=\omega\sqrt{\mu\,\varepsilon}\), and \(C_1\) and \(C_2\) are two arbitrary constants.

Now, we need to determine \(C_1\) and \(C_2\). Firstly, as we have used \(e^{-i \omega t}\) as the time dependency of harmonic fields (see Basics), \(C_1\,e^{i k r}\) corresponds to an outgoing wave while \(C_2\,e^{-i k r}\) is an incoming wave. Since we assumed that the antenna is in a homogeneous space, electromagnetic waves go away from the antenna. Therefore, there is not any incoming wave. So, \(C_2\) is zero.

\begin{equation}

\Rightarrow A_z(r,\omega)=C_1\,\frac{ e^{i k r}}{r} \tag{7}\label{AC1}

\end{equation}

Note that equation [\ref{AC1}] is valid for \(r\neq 0\) and we need to determine \(C_1\) in a way that equation [\ref{AC1}] satisfies equation [\ref{AzJ}] at \(r=0\) as well (source condition). This step is a little bit tricky. Since \(r=0\) is a singular point of \(A_z(r,\omega)\), it is not possible to discuss about the value of \(A_z(r,\omega)\) exactly at \(r=0\) and relate it to the source. However, \(A_z(r,\omega)\) is integrable over a volume around \(r=0\), let say a small sphere with radius \(a\). Therefore, we can relate the integral value of \(\mathbf{A}(\mathbf{r},\omega)\) to the integral value of source \(\big(\mathbf{J}\left(\mathbf{r},\omega\right)\big)\) over this volume to determine \(C_1\). To do that, we integrate both sides of the equation [\ref{AzJ}] over a very small sphere around \(r=0\) with radius \(a\) as it approaches zero. For every term in equation [\ref{AzJ}] we calculate the integral as shown below.

\begin{eqnarray*}

\int_{v_a}\nabla^2 A_z\left(\mathbf{r},\omega\right)dv &=&\int_{v_a}\nabla \cdot \big(\nabla A_z\left(\mathbf{r},\omega\right)\big)dv\\&=& \int_{s_a}\nabla A_z\left(\mathbf{r},\omega\right).d\mathbf{s}\\&=&\int_{s_a}C_1\,e^{ikr}\left(\frac{ik}{r}-\frac{1}{r^2}\right)r^2 \sin(\theta)\,d\theta\,d\varphi\\&=& 4 \pi C_1\,e^{ika}\left(ika-1\right)\\

\end{eqnarray*}

\begin{equation}

\Rightarrow\lim_{a \to 0} \int_{v_a}\nabla^2 A_z\left(\mathbf{r},\omega\right)dv=-4 \pi C_1 \tag{8}\label{a1}

\end{equation}

\begin{equation}

\lim_{a \to 0}\int_{v_a}\omega^2\mu\,\varepsilon\,A_z\left(\mathbf{r},\omega\right)dv =\omega^2\mu\,\varepsilon\lim_{a \to 0}\int_{v_a}C_1\,\frac{e^{ikr}}{r} r^2 \sin(\theta)\,d\theta\,d\varphi\,dr=0 \tag{9}\label{a2}

\end{equation}

\begin{equation}

\lim_{a \to 0}\int_{v_a}-\mu\,J_0 \, \delta (\mathbf{r})dv =-\mu\,J_0 \lim_{a \to 0}\int_{v_a}\delta (\mathbf{r})dv=-\mu\,J_0 \tag{10}\label{a3}

\end{equation}

By putting the results of equation [\ref{a1}], [\ref{a2}], and [\ref{a3}] together, \(C_1\) can be found

\[C_1=\frac{\mu\,J_0}{4\pi}\]

and therefore,

\[A_z(r,\omega)=\frac{\mu\,J_0}{4\pi}\,\frac{ e^{i k r}}{r}\]

\begin{equation}

\Rightarrow \mathbf{A}(\mathbf{r},\omega)=\frac{\mu\,\mathbf{J_0}}{4\pi}\,\frac{ e^{i k r}}{r} \tag{11}

\end{equation}

Now, one can calculate the magnetic field easily as follows.

\begin{eqnarray*}

\mathbf{H}\left(\mathbf{r},\omega\right)&=&\frac{1}{\mu} \nabla\times\mathbf{A}\left(\mathbf{r},\omega\right)\\&=&\frac{1}{\mu} \nabla\times\left(\frac{\mu\,J_0}{4\pi}\,\frac{ e^{i k r}}{r}\mathbf{\hat{z}}\right)\\&=&\frac{J_0}{4\pi}\,\nabla\left(\frac{e^{i k r}}{r}\right)\times \mathbf{\hat{z}}\\&=&\frac{J_0}{4\pi}e^{i k r}\left(\frac{ik}{r}-\frac{1}{r^2}\right)\mathbf{\hat{r}}\times \mathbf{\hat{z}}

\end{eqnarray*}

\begin{equation}

\Rightarrow\mathbf{H}\left(\mathbf{r},\omega\right)=-\frac{ikJ_0}{4\pi}\frac{e^{i k r}}{r} \left(1+\frac{i}{kr}\right) \sin(\theta)\,\hat{\varphi} \tag{12}\label{H}

\end{equation}

to calculate the electric field, we do the following.

\begin{eqnarray*}\nabla\big(\nabla\cdot\mathbf{A}\left(\mathbf{r},\omega\right)\big)&=&\nabla\left[\frac{\mu J_0}{4\pi}e^{i k r}\left(\frac{i k}{r}-\frac{1}{r^2}\right)\mathbf{\hat{r}}\cdot \mathbf{\hat{z}}\right]\\&=& \frac{\mu J_0}{4\pi}\nabla\left[e^{i k r}\left(\frac{i k}{r}-\frac{1}{r^2}\right)\\cos(\theta)\right]\\&=& \frac{\mu J_0}{4\pi}e^{i k r}\left[\\cos(\theta)\left(-\frac{k^2}{r}-\frac{2ik}{r^2}+\frac{2}{r^3}\right)\mathbf{\hat{r}}+\sin(\theta)\left(-\frac{ik}{r^2}+\frac{1}{r^3}\right)\mathbf{\hat{\theta}}\right]\end{eqnarray*} \begin{eqnarray}\Rightarrow\mathbf{E}\left(\mathbf{r},\omega\right)&=&i \omega\mathbf{A}\left(\mathbf{r},\omega\right)-\frac{\nabla\big(\nabla\cdot\mathbf{A}\left(\mathbf{r},\omega\right)\big)}{i\omega\mu\,\varepsilon}\nonumber\\&=&\frac{i\omega\mu J_0}{4\pi}\frac{e^{i k r}}{r}\left[\\cos(\theta)\left(-\frac{2i}{kr}+\frac{2}{(kr)^2}\right)\mathbf{\hat{r}}+\sin(\theta)\left(-1-\frac{i}{kr}+\frac{1}{(kr)^2}\right)\mathbf{\hat{\theta}}\right] \tag{13}\label{E}\end{eqnarray}

2. Lets first clarify what we mean by radiation pattern of an antenna. Generally, the radiation pattern is defined as the radiation intensity (radial component of the Poynting vector) as a function of direction in space (\(\varphi\) and \(\theta\)). Since in reality the receiver and transmitter antennas are located very far from each other, the calculation and measurement of radiation patterns are usually done using the fields in far distances form antennas. Hence, the fields are called far fields and the resultant pattern is called far filed radiation pattern or briefly far field pattern (we will see in Problem One how to precisely define the far field zone of an antenna). It should be noted that the radiation pattern is usually normalized to its maximum.

Based on what is mentioned, firstly, we need to calculate far fields. According to equation [\ref{H}] and [\ref{E}], for \(kr\gg 1\), the electric and magnetic fields can be approximated as below.

\begin{equation}

\mathbf{H_{ff}}\left(\mathbf{r},\omega\right)\cong \frac{-ik J_0}{4\pi}\frac{e^{i k r}}{r}\sin(\theta)\,\hat{\varphi} \tag{14}\label{Hff}

\end{equation}

\begin{equation}

\mathbf{E_{ff}}\left(\mathbf{r},\omega\right)\cong \frac{-i\omega\mu J_0}{4\pi}\frac{e^{i k r}}{r}\sin(\theta)\mathbf{\hat{\theta}} \tag{15}\label{Eff}

\end{equation}

This approximation is only valid for fields at large distance form the antenna (\(kr\gg 1\)).This condition can also be rewritten in the following way.

\begin{equation}

kr\gg 1 \Rightarrow r\gg \frac{1}{k}=\frac{1}{\omega \sqrt(\mu\varepsilon)}=\frac{c}{\omega}=\frac{\lambda}{2 \pi}\Rightarrow r\gg \lambda \tag{16}

\end{equation}

where \(\lambda\) is the wavelength of the wave in the medium. As can be seen, far field zone for an infinitesimal antenna is all points of space for which the distance to the antenna is much larger than wavelength.

Now, Poynting vector can be calculated as below.

\begin{eqnarray}

\mathbf{S_{ff}}\left(\mathbf{r},\omega\right)&=&\frac{1}{2}\,\text{Real}\Big\{\mathbf{E_{ff}}\left(\mathbf{r},\omega\right)\times \mathbf{H_{ff}^*}\left(\mathbf{r},\omega\right)\Big\} \nonumber\\&=&\frac{ \mu^\frac{3}{2} \varepsilon^\frac{1}{2}\omega^2 J_0^2}{32 \pi ^2}\frac{1}{r^2}\,sin^2(\theta)\mathbf{\hat{r}}\tag{17}\label{Sff}

\end{eqnarray}

As can be seen the intensity drops as \(\frac{1}{r^2}\) which is expected from the law of conservation of energy.

The far field radiation pattern can be calculated easily by normalizing equation [\ref{Sff}].

\begin{equation}

S^n_{f\!f}\left(\mathbf{\theta,\varphi}\right)=sin^2(\theta) \tag{18}

\end{equation}

Here, the superscript \textit{'n'} denotes the normalization. The far filed radiation pattern is depicted in figure (1) in three dimensions.

Conventionally, radiation patterns are drawn in some cut-planes (there are some standard planes which we will define in Problem One), since it is easier to see the details of the pattern which are not visible in a three dimensional figure. Here, due to rotational symmetry, it is sufficient to draw the pattern only in x-z plane. However, to show the rotational symmetry it is better to draw the pattern in x-y plane as well. The far field radiation pattern in x-y and x-z planes are shown in figure (2) and (3).

3. In this last step we just want to introduce some simple antenna parameters and the way to calculate them, such as half power beamwidth (HW) and "directivity". All of these parameters are calculated using the Poynting vector, therefore, they just provide different pictures of the radiation characteristics of the antenna.

HW, in a determined plane, is defined as the width of the angular interval in which the radiation intensity is more than half of the maximum radiation intensity.

According to equation [\ref{Sff}], one can calculate the HW of the antenna in x-z plane as below.

\[\frac{ \mu^\frac{3}{2} \varepsilon^\frac{1}{2}\omega^2 J_0^2}{32 \pi ^2}\frac{1}{r^2}\,sin^2(\theta) \frac{1}{2}\times\frac{ \mu^\frac{3}{2} \varepsilon^\frac{1}{2}\omega^2 J_0^2}{32 \pi ^2}\frac{1}{r^2}\]

\[\Rightarrow sin^2(\theta)\frac{1}{2}\Rightarrow \theta\in (-45^\circ,45^\circ)\]

\[\Rightarrow HW=90^\circ\]

This a measure which means that the antenna radiates effectively over 90 degree in the x-z plane which is a quite large angle. Although HW is an important parameter, it does not indicate in which directions the antenna has higher radiation. Therefore, we define "directivity" to address this question.

Since the radiated intensity of an antenna is a function of the direction of radiation, a criterion is needed to characterize and quantify this dependency. This criterion is the so-called "directivity function" and defined as below.

\begin{equation}

D(\theta,\varphi)=\frac{\text{radiation intensity}}{\text{average of radiation intensity over all directions}}\tag{19}\label{Dirc}

\end{equation}

Physically, this function describes the ability of an antenna to send power in a specific direction and compares it to the average power radiated in all directions.

Using equation [\ref{Dirc}], the directivity of our infinitesimal antenna can be obtained as below.

\begin{eqnarray}

D(\theta,\varphi)&=&\frac{S^n_{f\!f}\left(\mathbf{\theta,\varphi}\right)}{\frac{1}{4\pi}\int_0^{2\pi}\int_0^{\pi}S^n_{f\!f}\left(\mathbf{\theta,\varphi}\right)\sin(\theta) d\theta d\varphi}\nonumber\\&=&\frac{sin^2(\theta)}{\frac{1}{4\pi}\int_0^{2\pi}\int_0^{\pi}sin^3(\theta) d\theta d\varphi}\nonumber\\&=&\frac{3}{2}\,sin^2(\theta) \tag{20}\label{D}

\end{eqnarray}

For example the directivity of the antenna in the direction of \(\theta=60^\circ, \varphi=30^\circ\) is equal to 1.125, which mean that the antenna radiates 1.125 times larger than the average radiation in this direction( note that the maximum directivity for an infinitesimal antenna is 1.5). As can be seen in equation [\ref{D}], the diectivity is not a function of azimuth angle. Therefore, it is said that the antenna is omni-directional, which means that there is one plane in which the directitity is constant (here x-y plane).