# 9.3: Solution of the Wave Equation for Magnetic Vector Potential

- Page ID
- 24834

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

The magnetic vector potential \(\widetilde{\bf A}\) due to a current density \(\widetilde{\bf J}\) is given by the following wave equation:

\[\nabla^2 \widetilde{\bf A} -\gamma^2\widetilde{\bf A} = -\mu\widetilde{\bf J} \label{m0196_ePDEA} \]

where \(\gamma\) is the *propagation constant*, defined in the usual manner^{1}

\[\gamma^2\triangleq -\omega^2\mu\epsilon \nonumber \]

Equation \ref{m0196_ePDEA} is a partial differential equation which is inhomogeneous (in the mathematical sense) and can be solved given appropriate boundary conditions. In this section, we present the solution for arbitrary distributions of current in free space.

The strategy is to first identify a solution for a distribution of current that exists at a single point. This distribution is the *current moment*. An example of a current moment is shown in Figure \(\PageIndex{1}\).

A general expression for a current moment located at the origin is:

\[\widetilde{\bf J}({\bf r}) = \hat{\bf l}~\widetilde{I}~\Delta l~\delta({\bf r}) \label{m0196_eCM} \]

where \(\widetilde{I}\) has units of current (SI base units of A), \(\Delta l\) has units of length (SI base units of m), \(\hat{\bf l}\) is the direction of current flow, and \(\delta({\bf r})\) is the volumetric sampling function^{2} defined as follows:

\begin{align} ~ &\delta({\bf r}) \triangleq 0~~~ \mbox{for}~~~ {\bf r}\neq 0;~~~ \mbox{and} \label{m0196_eDelta1} \\ ~ &\int_{\mathcal{V}}\delta({\bf r})~dv \triangleq 1 \label{m0196_eDelta2}\end{align}

where \(\mathcal{V}\) is any volume which includes the origin (\({\bf r}=0\)). It is evident from Equation \ref{m0196_eDelta2} that \(\delta({\bf r})\) has SI base units of m\(^{-3}\); i.e., inverse volume. Subsequently \({\bf J}({\bf r})\) has SI base units of A/m\(^2\), indicating that it is a volume current density.^{3} – as indicated by Equation \ref{m0196_eDelta1} – it exists only at the origin and nowhere else.

Substituting Equation \ref{m0196_eCM} into Equation \ref{m0196_ePDEA}, we obtain:

\[\nabla^2 \widetilde{\bf A} -\gamma^2\widetilde{\bf A} = -\mu\hat{\bf l}~\widetilde{I}~\Delta l~\delta({\bf r}) \label{m0196_ePDEAC} \]

The general solution to this equation presuming homogeneous and time-invariant media (i.e., \(\mu\) and \(\epsilon\) constant with respect to space and time) is:

\[\widetilde{\bf A}({\bf r}) = \hat{\bf l}~\mu~\widetilde{I}~\Delta l~\frac{e^{\pm \gamma r}}{4\pi r} \label{m0196_eA} \]

To confirm that this is a solution to Equation \ref{m0196_ePDEAC}, substitute this expression into Equation \ref{m0196_ePDEA} and observe that the equality holds.

Note that Equation \ref{m0196_eA} indicates two solutions, corresponding to the signs of the exponent in the factor \(e^{\pm \gamma r}\). We can actually be a little more specific by applying some common-sense physics. Recall that \(\gamma\) in general may be expressed in terms of real and imaginary components as follows:

\[\gamma = \alpha + j\beta \nonumber \]

where \(\alpha\) is a real-valued positive constant known as the *attenuation constant* and \(\beta\) is a real-valued positive constant known as the *phase propagation constant*. Thus, we may rewrite Equation \ref{m0196_eA} as follows:

\[\widetilde{\bf A}({\bf r}) = \hat{\bf l}~\mu~\widetilde{I}~\Delta l~\frac{ e^{\pm \alpha r} e^{\pm j\beta r} }{4\pi r} \label{m0196_eA2} \]

Now consider the factor \(e^{\pm \alpha r}/r\), which determines the dependence of magnitude on distance \(r\). If we choose the negative sign in the exponent, this factor decays exponentially with increasing distance from the origin, ultimately reaching zero at \(r\rightarrow\infty\). This is precisely the expected behavior, since we expect the magnitude of a radiated field to diminish with increasing distance from the source. If on the other hand we choose the positive sign in the exponent, this factor increases to infinity as \(r\rightarrow\infty\). That outcome can be ruled out on physical grounds, since it implies the introduction of energy independently from the source. The requirement that field magnitude diminishes to zero as distance from the source increases to infinity is known as the *radiation condition*, and is essentially a boundary condition that applies at \(r\rightarrow\infty\). Invoking the radiation condition, Equation \ref{m0196_eA} becomes:

\[\widetilde{\bf A}({\bf r}) = \hat{\bf l}~\mu~\widetilde{I}~\Delta l~\frac{e^{-\gamma r}}{4\pi r} \label{m0196_eA3} \]

In the loss-free (\(\alpha=0\)) case, we cannot rely on the radiation condition to constrain the sign of \(\gamma\). However, in practical engineering work there is always some media loss; i.e., \(\alpha\) might be negligible but is not quite zero. Thus, we normally assume that the solution for lossless (including free space) conditions is given by Equation \ref{m0196_eA3} with \(\alpha=0\):

\[\widetilde{\bf A}({\bf r}) = \hat{\bf l}~\mu~\widetilde{I}~\Delta l~\frac{e^{-j \beta r}}{4\pi r} \label{m0196_eA4} \]

Now consider the factor \(e^{-j \beta r}\) in Equation \ref{m0196_eA4}. This factor exclusively determines the dependence of the phase of \(\widetilde{\bf A}({\bf r})\) with increasing distance \(r\) from the source. In this case, we observe that surfaces of constant phase correspond to spherical shells which are concentric with the source. Thus, \(\widetilde{\bf A}({\bf r})\) is a spherical wave.

Let us now consider a slightly more complicated version of the problem in which the current moment is no longer located at the origin, but rather is located at \({\bf r}'\). This is illustrated in Figure \(\PageIndex{2}\).

Figure \(\PageIndex{2}\): Current moment displaced from the origin. ( CC BY-SA 4.0; C. Wang)This current distribution can be expressed as follows:

\[\widetilde{\bf J}({\bf r}) = \hat{\bf l}~\widetilde{I}~\Delta l~\delta({\bf r}-{\bf r}') %\label{m0196_eCM2} \]

The solution in this case amounts to a straightforward modification of the existing solution. To see this, note that \(\widetilde{\bf A}\) depends only on \(r\), and not at all on \(\theta\) or \(\phi\). In other words, \(\widetilde{\bf A}\) depends only on the distance between the “field point” \({\bf r}\) at which we observe \(\widetilde{\bf A}\), and the “source point” \({\bf r}'\) at which the current moment lies. Thus we replace \(r\) with this distance, which is \(\left|{\bf r}-{\bf r}'\right|\). The solution becomes:

\[\widetilde{\bf A}({\bf r}) = \hat{\bf l}~\mu~\widetilde{I}~\Delta l~\frac{e^{-\gamma \left|{\bf r}-{\bf r}'\right|}}{4\pi \left|{\bf r}-{\bf r}'\right|} \nonumber \]

Continuing to generalize the solution, let us now consider a scenario consisting of a filament of current following a path \(\mathcal{C}\) through space. This is illustrated in Figure \(\PageIndex{3}\).

Figure \(\PageIndex{3}\): A filament of current lying along the path \(\mathcal{C}\). This current distribution may be interpreted as a collection of current moments lying along \(\mathcal{C}\). ( CC BY-SA 4.0; C. Wang)Such a filament may be viewed as a collection of a large number \(N\) of discrete current moments distributed along the path. The contribution of the \(n^{\mbox{th}}\) current moment, located at \({\bf r}_n\), to the total magnetic vector potential is:

\[\Delta \widetilde{\bf A}({\bf r};{\bf r}_n) = \hat{\bf l}({\bf r}_n)~\mu~\widetilde{I}({\bf r}_n)~\Delta l~\frac{e^{-\gamma \left|{\bf r}-{\bf r}_n\right|}}{4\pi \left|{\bf r}-{\bf r}_n\right|} \nonumber \]

Note that we are allowing both the current \(\widetilde{I}\) and current direction \(\hat{\bf l}\) to vary with position along \(\mathcal{C}\). Assuming the medium is linear, superposition applies. So we add these contributions to obtain the total magnetic vector potential:

\begin{align} \widetilde{\bf A}({\bf r}) &\approx \sum_{n=1}^{N} \Delta \widetilde{\bf A}({\bf r};{\bf r}_n) \\ &\approx \frac{\mu}{4\pi} ~\sum_{n=1}^N \hat{\bf l}({\bf r}_n)~\widetilde{I}({\bf r}_n) ~\frac{e^{-\gamma \left|{\bf r}-{\bf r}_n\right|}}{\left|{\bf r}-{\bf r}_n\right|}~\Delta l\end{align}

Now letting \(\Delta l\rightarrow 0\) so that we may replace \(\Delta l\) with the differential length \(dl\):

\[\boxed{ \widetilde{\bf A}({\bf r}) = \frac{\mu}{4\pi}~\int_{\mathcal{C}} \hat{\bf l}({\bf r}')~\widetilde{I}({\bf r}')~\frac{e^{-\gamma \left|{\bf r}-{\bf r}'\right|}}{\left|{\bf r}-{\bf r}'\right|} ~dl } \label{m0196_eAL} \]

where we have also replaced \({\bf r}_n\) with the original notation \({\bf r}'\) since we once again have a continuum (as opposed to a discrete set) of source locations.

The magnetic vector potential corresponding to radiation from a line distribution of current is given by Equation \ref{m0196_eAL}.

Given \(\widetilde{\bf A}({\bf r})\), the magnetic and electric fields may be determined using the procedure developed in Section 9.2.

**Additional Reading:**

- “Magnetic potential” on Wikipedia.
- “Dirac delta function” on Wikipedia.

- Alternatively, \(\gamma^2\triangleq -\omega^2\mu\epsilon_c\) accounting for the possibility of lossy media.↩
- Also a form of the
*Dirac delta function*; see “Additional Reading” at the end of this section.↩ - Remember, the density of a volume current is with respect to the
*area*through which it flows, therefore the units are A/m\(^2\).↩