Laplacian

Spherical gradient, divergence, curl and Laplacian

November 9, 2016 math and physics play , , , , , , , , , ,

[Click here for a PDF of this post with nicer formatting]

Unit vectors

Two of the spherical unit vectors we can immediately write by inspection.

\begin{equation}\label{eqn:sphericalLaplacian:20}
\begin{aligned}
\rcap &= \Be_1 \sin\theta \cos\phi + \Be_2 \sin\theta \sin\phi + \Be_3 \cos\theta \\
\phicap &= -\Be_1 \sin\theta + \Be_2 \cos\phi
\end{aligned}
\end{equation}

We can compute \( \thetacap \) by utilizing the right hand triplet property

\begin{equation}\label{eqn:sphericalLaplacian:40}
\begin{aligned}
\thetacap
&=
\phicap \cross \rcap \\
&=
\begin{vmatrix}
\Be_1 & \Be_2 & \Be_3 \\
-S_\phi & C_\phi & 0 \\
S_\theta C_\phi & S_\theta S_\phi & C_\theta \\
\end{vmatrix} \\
&=
\Be_1 \lr{ C_\theta C_\phi }
+\Be_2 \lr{ C_\theta S_\phi }
+\Be_3 \lr{ -S_\theta \lr{ S_\phi^2 + C_\phi^2 } } \\
&=
\Be_1 \cos\theta \cos\phi
+\Be_2 \cos\theta \sin\phi
-\Be_3 \sin\theta.
\end{aligned}
\end{equation}

Here I’ve used \( C_\theta = \cos\theta, S_\phi = \sin\phi, \cdots \) as a convenient shorthand. Observe that with \( i = \Be_1 \Be_2 \), these unit vectors admit a small factorization that makes further manipulation easier

\begin{equation}\label{eqn:sphericalLaplacian:80}
\boxed{
\begin{aligned}
\rcap &= \Be_1 e^{i\phi} \sin\theta + \Be_3 \cos\theta \\
\thetacap &= \cos\theta \Be_1 e^{i\phi} – \sin\theta \Be_3 \\
\phicap &= \Be_2 e^{i\phi}
\end{aligned}
}
\end{equation}

It should also be the case that \( \rcap \thetacap \phicap = I \), where \( I = \Be_1 \Be_2 \Be_3 = \Be_{123}\) is the \R{3} pseudoscalar, which is straightforward to check

\begin{equation}\label{eqn:sphericalLaplacian:60}
\begin{aligned}
\rcap \thetacap \phicap
&=
\lr{ \Be_1 e^{i\phi} \sin\theta + \Be_3 \cos\theta }
\lr{ \cos\theta \Be_1 e^{i\phi} – \sin\theta \Be_3 }
\Be_2 e^{i\phi} \\
&=
\lr{ \sin\theta \cos\theta – \cos\theta \sin\theta + \Be_{31} e^{i\phi} \lr{ \cos^2\theta + \sin^2\theta } }
\Be_2 e^{i\phi} \\
&=
\Be_{31} \Be_2 e^{-i\phi} e^{i\phi} \\
&=
\Be_{123}.
\end{aligned}
\end{equation}

This property could also have been used to compute \(\thetacap\).

Gradient

To compute the gradient, note that the coordinate vectors for the spherical parameterization are
\begin{equation}\label{eqn:sphericalLaplacian:120}
\begin{aligned}
\Bx_r
&= \PD{r}{\Br} \\
&= \PD{r}{\lr{r \rcap}} \\
&= \rcap + r \PD{r}{\rcap} \\
&= \rcap,
\end{aligned}
\end{equation}
\begin{equation}\label{eqn:sphericalLaplacian:140}
\begin{aligned}
\Bx_\theta
&= \PD{\theta}{\lr{r \rcap} } \\
&= r \PD{\theta}{} \lr{ S_\theta \Be_1 e^{i\phi} + C_\theta \Be_3 } \\
&= r \PD{\theta}{} \lr{ C_\theta \Be_1 e^{i\phi} – S_\theta \Be_3 } \\
&= r \thetacap,
\end{aligned}
\end{equation}
\begin{equation}\label{eqn:sphericalLaplacian:160}
\begin{aligned}
\Bx_\phi
&= \PD{\phi}{\lr{r \rcap} } \\
&= r \PD{\phi}{} \lr{ S_\theta \Be_1 e^{i\phi} + C_\theta \Be_3 } \\
&= r S_\theta \Be_2 e^{i\phi} \\
&= r \sin\theta \phicap.
\end{aligned}
\end{equation}

Since these are all normal, the dual vectors defined by \( \Bx^j \cdot \Bx_k = \delta^j_k \), can be obtained by inspection
\begin{equation}\label{eqn:sphericalLaplacian:180}
\begin{aligned}
\Bx^r &= \rcap \\
\Bx^\theta &= \inv{r} \thetacap \\
\Bx^\phi &= \inv{r \sin\theta} \phicap.
\end{aligned}
\end{equation}

The gradient follows immediately
\begin{equation}\label{eqn:sphericalLaplacian:200}
\spacegrad =
\Bx^r \PD{r}{} +
\Bx^\theta \PD{\theta}{} +
\Bx^\phi \PD{\phicap}{},
\end{equation}

or
\begin{equation}\label{eqn:sphericalLaplacian:240}
\boxed{
\spacegrad
=
\rcap \PD{r}{} +
\frac{\thetacap}{r} \PD{\theta}{} +
\frac{\phicap}{r\sin\theta} \PD{\phicap}{}.
}
\end{equation}

More information on this general dual-vector technique of computing the gradient in curvilinear coordinate systems can be found in
[2].

Partials

To compute the divergence, curl and Laplacian, we’ll need the partials of each of the unit vectors \( \PDi{\theta}{\rcap}, \PDi{\phi}{\rcap}, \PDi{\theta}{\thetacap}, \PDi{\phi}{\thetacap}, \PDi{\phi}{\phicap} \).

The \( \thetacap \) partials are

\begin{equation}\label{eqn:sphericalLaplacian:260}
\begin{aligned}
\PD{\theta}{\thetacap}
&=
\PD{\theta}{} \lr{
C_\theta \Be_1 e^{i\phi} – S_\theta \Be_3
} \\
&=
-S_\theta \Be_1 e^{i\phi} – C_\theta \Be_3 \\
&=
-\rcap,
\end{aligned}
\end{equation}
\begin{equation}\label{eqn:sphericalLaplacian:280}
\begin{aligned}
\PD{\phi}{\thetacap}
&=
\PD{\phi}{} \lr{
C_\theta \Be_1 e^{i\phi} – S_\theta \Be_3
} \\
&=
C_\theta \Be_2 e^{i\phi} \\
&=
C_\theta \phicap.
\end{aligned}
\end{equation}

The \( \phicap \) partials are

\begin{equation}\label{eqn:sphericalLaplacian:300}
\begin{aligned}
\PD{\theta}{\phicap}
&=
\PD{\theta}{} \Be_2 e^{i\phi} \\
&=
0.
\end{aligned}
\end{equation}
\begin{equation}\label{eqn:sphericalLaplacian:320}
\begin{aligned}
\PD{\phi}{\phicap}
&=
\PD{\phi}{} \Be_2 e^{i \phi} \\
&=
-\Be_1 e^{i \phi} \\
&=
-\rcap \gpgradezero{ \rcap \Be_1 e^{i \phi} }
– \thetacap \gpgradezero{ \thetacap \Be_1 e^{i \phi} }
– \phicap \gpgradezero{ \phicap \Be_1 e^{i \phi} } \\
&=
-\rcap \gpgradezero{ \lr{
\Be_1 e^{i\phi} S_\theta + \Be_3 C_\theta
} \Be_1 e^{i \phi} }
– \thetacap \gpgradezero{ \lr{
C_\theta \Be_1 e^{i\phi} – S_\theta \Be_3
} \Be_1 e^{i \phi} } \\
&=
-\rcap \gpgradezero{ e^{-i\phi} S_\theta e^{i \phi} }
– \thetacap \gpgradezero{ C_\theta e^{-i\phi} e^{i \phi} } \\
&=
-\rcap S_\theta
– \thetacap C_\theta.
\end{aligned}
\end{equation}

The \( \rcap \) partials are were computed as a side effect of evaluating \( \Bx_\theta \), and \( \Bx_\phi \), and are

\begin{equation}\label{eqn:sphericalLaplacian:340}
\PD{\theta}{\rcap}
=
\thetacap,
\end{equation}
\begin{equation}\label{eqn:sphericalLaplacian:360}
\PD{\phi}{\rcap}
=
S_\theta \phicap.
\end{equation}

In summary
\begin{equation}\label{eqn:sphericalLaplacian:380}
\boxed{
\begin{aligned}
\partial_{\theta}{\rcap} &= \thetacap \\
\partial_{\phi}{\rcap} &= S_\theta \phicap \\
\partial_{\theta}{\thetacap} &= -\rcap \\
\partial_{\phi}{\thetacap} &= C_\theta \phicap \\
\partial_{\theta}{\phicap} &= 0 \\
\partial_{\phi}{\phicap} &= -\rcap S_\theta – \thetacap C_\theta.
\end{aligned}
}
\end{equation}

Divergence and curl.

The divergence and curl can be computed from the vector product of the spherical coordinate gradient and the spherical representation of a vector. That is

\begin{equation}\label{eqn:sphericalLaplacian:400}
\spacegrad \BA
= \spacegrad \cdot \BA + \spacegrad \wedge \BA
= \spacegrad \cdot \BA + I \spacegrad \cross \BA.
\end{equation}

That gradient vector product is

\begin{equation}\label{eqn:sphericalLaplacian:420}
\begin{aligned}
\spacegrad \BA
&=
\lr{
\rcap \partial_{r}
+ \frac{\thetacap}{r} \partial_{\theta}
+ \frac{\phicap}{rS_\theta} \partial_{\phi}
}
\lr{ \rcap A_r + \thetacap A_\theta + \phicap A_\phi} \\
&=
\rcap \partial_{r}
\lr{ \rcap A_r + \thetacap A_\theta + \phicap A_\phi} \\
&+ \frac{\thetacap}{r} \partial_{\theta}
\lr{ \rcap A_r + \thetacap A_\theta + \phicap A_\phi} \\
&+ \frac{\phicap}{rS_\theta} \partial_{\phicap}
\lr{ \rcap A_r + \thetacap A_\theta + \phicap A_\phi} \\
&=
\lr{ \partial_r A_r + \rcap \thetacap \partial_r A_\theta + \rcap \phicap \partial_r A_\phi} \\
&+ \frac{1}{r}
\lr{
\thetacap (\partial_\theta \rcap) A_r + \thetacap (\partial_\theta \thetacap) A_\theta + \thetacap (\partial_\theta \phicap) A_\phi
+\thetacap \rcap \partial_\theta A_r + \partial_\theta A_\theta + \thetacap \phicap \partial_\theta A_\phi
} \\
&+ \frac{1}{rS_\theta}
\lr{
\phicap (\partial_\phi \rcap) A_r + \phicap (\partial_\phi \thetacap) A_\theta + \phicap (\partial_\phi \phicap) A_\phi
+\phicap \rcap \partial_\phi A_r + \phicap \thetacap \partial_\phi A_\theta + \partial_\phi A_\phi
} \\
&=
\lr{ \partial_r A_r + \rcap \thetacap \partial_r A_\theta + \rcap \phicap \partial_r A_\phi} \\
&+ \frac{1}{r}
\lr{
\thetacap (\thetacap) A_r + \thetacap (-\rcap) A_\theta + \thetacap (0) A_\phi
+\thetacap \rcap \partial_\theta A_r + \partial_\theta A_\theta + \thetacap \phicap \partial_\theta A_\phi
} \\
&+ \frac{1}{r S_\theta}
\lr{
\phicap (S_\theta \phicap) A_r + \phicap (C_\theta \phicap) A_\theta – \phicap (\rcap S_\theta + \thetacap C_\theta) A_\phi
+\phicap \rcap \partial_\phi A_r + \phicap \thetacap \partial_\phi A_\theta + \partial_\phi A_\phi
}.
\end{aligned}
\end{equation}

The scalar component of this is the divergence
\begin{equation}\label{eqn:sphericalLaplacian:440}
\begin{aligned}
\spacegrad \cdot \BA
&=
\partial_r A_r
+ \frac{A_r}{r}
+ \inv{r} \partial_\theta A_\theta
+ \frac{1}{r S_\theta}
\lr{ S_\theta A_r + C_\theta A_\theta + \partial_\phi A_\phi
} \\
&=
\partial_r A_r
+ 2 \frac{A_r}{r}
+ \inv{r} \partial_\theta A_\theta
+ \frac{1}{r S_\theta}
C_\theta A_\theta
+ \frac{1}{r S_\theta} \partial_\phi A_\phi \\
&=
\partial_r A_r
+ 2 \frac{A_r}{r}
+ \inv{r} \partial_\theta A_\theta
+ \frac{1}{r S_\theta}
C_\theta A_\theta
+ \frac{1}{r S_\theta} \partial_\phi A_\phi,
\end{aligned}
\end{equation}

which can be factored as
\begin{equation}\label{eqn:sphericalLaplacian:460}
\boxed{
\spacegrad \cdot \BA
=
\inv{r^2} \partial_r (r^2 A_r)
+ \inv{r S_\theta} \partial_\theta (S_\theta A_\theta)
+ \frac{1}{r S_\theta} \partial_\phi A_\phi.
}
\end{equation}

The bivector grade of \( \spacegrad \BA \) is the bivector curl
\begin{equation}\label{eqn:sphericalLaplacian:480}
\begin{aligned}
\spacegrad \wedge \BA
&=
\lr{
\rcap \thetacap \partial_r A_\theta + \rcap \phicap \partial_r A_\phi
} \\
&\quad + \frac{1}{r}
\lr{
\thetacap (-\rcap) A_\theta
+\thetacap \rcap \partial_\theta A_r + \thetacap \phicap \partial_\theta A_\phi
} \\
&\quad +
\frac{1}{r S_\theta}
\lr{
-\phicap (\rcap S_\theta + \thetacap C_\theta) A_\phi
+\phicap \rcap \partial_\phi A_r + \phicap \thetacap \partial_\phi A_\theta
} \\
&=
\lr{
\rcap \thetacap \partial_r A_\theta – \phicap \rcap \partial_r A_\phi
} \\
&\quad + \frac{1}{r}
\lr{
\rcap \thetacap A_\theta
-\rcap \thetacap \partial_\theta A_r + \thetacap \phicap \partial_\theta A_\phi
} \\
&\quad +
\frac{1}{r S_\theta}
\lr{
-\phicap \rcap S_\theta A_\phi + \thetacap \phicap C_\theta A_\phi
+\phicap \rcap \partial_\phi A_r – \thetacap \phicap \partial_\phi A_\theta
} \\
&=
\thetacap \phicap \lr{
\inv{r S_\theta} C_\theta A_\phi
+\frac{1}{r} \partial_\theta A_\phi
-\frac{1}{r S_\theta} \partial_\phi A_\theta
} \\
&\quad +\phicap \rcap \lr{
-\partial_r A_\phi
+
\frac{1}{r S_\theta}
\lr{
-S_\theta A_\phi
+ \partial_\phi A_r
}
} \\
&\quad +\rcap \thetacap \lr{
\partial_r A_\theta
+ \frac{1}{r} A_\theta
– \inv{r} \partial_\theta A_r
} \\
&=
I
\rcap \lr{
\inv{r S_\theta} \partial_\theta (S_\theta A_\phi)
-\frac{1}{r S_\theta} \partial_\phi A_\theta
}
+ I \thetacap \lr{
\frac{1}{r S_\theta} \partial_\phi A_r
-\inv{r} \partial_r (r A_\phi)
}
+ I \phicap \lr{
\inv{r} \partial_r (r A_\theta)
– \inv{r} \partial_\theta A_r
}
\end{aligned}
\end{equation}

This gives
\begin{equation}\label{eqn:sphericalLaplacian:500}
\boxed{
\spacegrad \cross \BA
=
\rcap \lr{
\inv{r S_\theta} \partial_\theta (S_\theta A_\phi)
-\frac{1}{r S_\theta} \partial_\phi A_\theta
}
+ \thetacap \lr{
\frac{1}{r S_\theta} \partial_\phi A_r
-\inv{r} \partial_r (r A_\phi)
}
+ \phicap \lr{
\inv{r} \partial_r (r A_\theta)
– \inv{r} \partial_\theta A_r
}.
}
\end{equation}

This and the divergence result above both check against the back cover of [1].

Laplacian

Using the divergence and curl it’s possible to compute the Laplacian from those, but we saw in cylindrical coordinates that it was much harder to do it that way than to do it directly.

\begin{equation}\label{eqn:sphericalLaplacian:540}
\begin{aligned}
\spacegrad^2 \psi
&=
\lr{
\rcap \partial_{r} +
\frac{\thetacap}{r} \partial_{\theta} +
\frac{\phicap}{r S_\theta} \partial_{\phi}
}
\lr{
\rcap \partial_{r} \psi
+ \frac{\thetacap}{r} \partial_{\theta} \psi
+ \frac{\phicap}{r S_\theta} \partial_{\phi} \psi
} \\
&=
\partial_{rr} \psi
+ \rcap \thetacap \partial_r \lr{ \inv{r} \partial_\theta \psi}
+ \rcap \phicap \inv{S_\theta} \partial_r \lr{ \inv{r} \partial_\phi \psi } \\
&
\quad + \frac{\thetacap}{r} \partial_{\theta} \lr{ \rcap \partial_{r} \psi }
+ \frac{\thetacap}{r^2} \partial_{\theta} \lr{ \thetacap \partial_{\theta} \psi }
+ \frac{\thetacap}{r^2} \partial_{\theta} \lr{ \frac{\phicap}{S_\theta} \partial_{\phi} \psi } \\
&
\quad + \frac{\phicap}{r S_\theta} \partial_{\phi} \lr{ \rcap \partial_{r} \psi }
+ \frac{\phicap}{r^2 S_\theta} \partial_{\phi} \lr{ \thetacap \partial_{\theta} \psi }
+ \frac{\phicap}{r^2 S_\theta^2} \partial_{\phi} \lr{ \phicap \partial_{\phi} \psi } \\
&=
\partial_{rr} \psi
+ \rcap \thetacap \partial_r \lr{ \inv{r} \partial_\theta \psi}
+ \rcap \phicap \inv{S_\theta} \partial_r \lr{ \inv{r} \partial_\phi \psi } \\
&
\quad + \frac{\thetacap\rcap}{r} \partial_{\theta} \lr{ \partial_{r} \psi }
+ \frac{1}{r^2} \partial_{\theta \theta} \psi
+ \frac{\thetacap \phicap}{r^2} \partial_{\theta} \lr{ \frac{1}{S_\theta} \partial_{\phi} \psi } \\
&
\quad + \frac{\phicap \rcap}{r S_\theta} \partial_{\phi r} \psi
+ \frac{\phicap\thetacap}{r^2 S_\theta} \partial_{\phi\theta} \psi
+ \frac{1}{r^2 S_\theta^2} \partial_{\phi \phi} \psi \\
&
\quad + \frac{\thetacap}{r} (\partial_\theta \rcap) \partial_{r} \psi
+ \frac{\thetacap}{r^2} (\partial_\theta \thetacap) \partial_{\theta} \psi
+ \frac{\thetacap}{r^2} (\partial_\theta \phicap) \frac{\phicap}{S_\theta} \partial_{\phi} \psi \\
&
\quad + \frac{\phicap}{r S_\theta} (\partial_\phi \rcap) \partial_{r} \psi
+ \frac{\phicap}{r^2 S_\theta} (\partial_\phi \thetacap) \partial_{\theta} \psi
+ \frac{\phicap}{r^2 S_\theta^2} (\partial_\phi \phicap) \partial_{\phi} \psi \\
&=
\partial_{rr} \psi
+ \rcap \thetacap \partial_r \lr{ \inv{r} \partial_\theta \psi}
+ \rcap \phicap \inv{S_\theta} \partial_r \lr{ \inv{r} \partial_\phi \psi } \\
&
\quad + \frac{\thetacap\rcap}{r} \partial_{\theta} \lr{ \partial_{r} \psi }
+ \frac{1}{r^2} \partial_{\theta \theta} \psi
+ \frac{\thetacap \phicap}{r^2} \partial_{\theta} \lr{ \frac{1}{S_\theta} \partial_{\phi} \psi } \\
&
\quad + \frac{\phicap \rcap}{r S_\theta} \partial_{\phi r} \psi
+ \frac{\phicap\thetacap}{r^2 S_\theta} \partial_{\phi\theta} \psi
+ \frac{1}{r^2 S_\theta^2} \partial_{\phi \phi} \psi \\
&
\quad + \frac{\thetacap}{r} (\thetacap) \partial_{r} \psi
+ \frac{\thetacap}{r^2} (-\rcap) \partial_{\theta} \psi
+ \frac{\thetacap}{r^2} (0) \frac{\phicap}{S_\theta} \partial_{\phi} \psi \\
&
\quad + \frac{\phicap}{r S_\theta} (S_\theta \phicap) \partial_{r} \psi
+ \frac{\phicap}{r^2 S_\theta} (C_\theta \phicap) \partial_{\theta} \psi
+ \frac{\phicap}{r^2 S_\theta^2} (-\rcap S_\theta – \thetacap C_\theta) \partial_{\phi} \psi
\end{aligned}
\end{equation}

All the bivector factors are expected to cancel out, but this should be checked. Those with an \( \rcap \thetacap \) factor are

\begin{equation}\label{eqn:sphericalLaplacian:560}
\partial_r \lr{ \inv{r} \partial_\theta \psi}
– \frac{1}{r} \partial_{\theta r} \psi
+ \frac{1}{r^2} \partial_{\theta} \psi
=
-\inv{r^2} \partial_\theta \psi
+\inv{r} \partial_{r \theta} \psi
– \frac{1}{r} \partial_{\theta r} \psi
+ \frac{1}{r^2} \partial_{\theta} \psi
= 0,
\end{equation}

and those with a \( \thetacap \phicap \) factor are
\begin{equation}\label{eqn:sphericalLaplacian:580}
\frac{1}{r^2} \partial_{\theta} \lr{ \frac{1}{S_\theta} \partial_{\phi} \psi }
– \frac{1}{r^2 S_\theta} \partial_{\phi\theta} \psi
+ \frac{1}{r^2 S_\theta^2} C_\theta \partial_{\phi} \psi
=
– \frac{1}{r^2} \frac{C_\theta}{S_\theta^2} \partial_{\phi} \psi
+ \frac{1}{r^2 S_\theta} \partial_{\theta \phi} \psi
– \frac{1}{r^2 S_\theta} \partial_{\phi\theta} \psi
+ \frac{1}{r^2 S_\theta^2} C_\theta \partial_{\phi} \psi
= 0,
\end{equation}

and those with a \( \phicap \rcap \) factor are
\begin{equation}\label{eqn:sphericalLaplacian:600}
– \inv{S_\theta} \partial_r \lr{ \inv{r} \partial_\phi \psi }
+ \frac{1}{r S_\theta} \partial_{\phi r} \psi
– \frac{1}{r^2 S_\theta^2} S_\theta \partial_{\phi} \psi
=
\inv{S_\theta} \frac{1}{r^2} \partial_\phi \psi
– \inv{r S_\theta} \partial_{r \phi} \psi
+ \frac{1}{r S_\theta} \partial_{\phi r} \psi
– \frac{1}{r^2 S_\theta} \partial_{\phi} \psi
= 0.
\end{equation}

This leaves
\begin{equation}\label{eqn:sphericalLaplacian:620}
\spacegrad^2 \psi
=
\partial_{rr} \psi
+ \frac{2}{r} \partial_{r} \psi
+ \frac{1}{r^2} \partial_{\theta \theta} \psi
+ \frac{1}{r^2 S_\theta} C_\theta \partial_{\theta} \psi
+ \frac{1}{r^2 S_\theta^2} \partial_{\phi \phi} \psi.
\end{equation}

This factors nicely as

\begin{equation}\label{eqn:sphericalLaplacian:640}
\boxed{
\spacegrad^2 \psi
=
\inv{r^2} \PD{r}{} \lr{ r^2 \PD{r}{ \psi} }
+ \frac{1}{r^2 \sin\theta} \PD{\theta}{} \lr{ \sin\theta \PD{\theta}{ \psi } }
+ \frac{1}{r^2 \sin\theta^2} \PDSq{\phi}{ \psi}
,
}
\end{equation}

which checks against the back cover of Jackson. Here it has been demonstrated explicitly that this operator expression is valid for multivector fields \( \psi \) as well as scalar fields \( \psi \).

References

[1] JD Jackson. Classical Electrodynamics. John Wiley and Sons, 2nd edition, 1975.

[2] A. Macdonald. Vector and Geometric Calculus. CreateSpace Independent Publishing Platform, 2012.

Gradient, divergence, curl and Laplacian in cylindrical coordinates

November 6, 2016 math and physics play , , , , , , , , , , , ,

[Click here for a PDF of this post with nicer formatting]

In class it was suggested that the identity

\begin{equation}\label{eqn:laplacianCylindrical:20}
\spacegrad^2 \BA =
\spacegrad \lr{ \spacegrad \cdot \BA }
-\spacegrad \cross \lr{ \spacegrad \cross \BA },
\end{equation}

can be used to compute the Laplacian in non-rectangular coordinates. Is that the easiest way to do this?

How about just sequential applications of the gradient on the vector? Let’s start with the vector product of the gradient and the vector. First recall that the cylindrical representation of the gradient is

\begin{equation}\label{eqn:laplacianCylindrical:80}
\spacegrad = \rhocap \partial_\rho + \frac{\phicap}{\rho} \partial_\phi + \zcap \partial_z,
\end{equation}

where
\begin{equation}\label{eqn:laplacianCylindrical:100}
\begin{aligned}
\rhocap &= \Be_1 e^{\Be_1 \Be_2 \phi} \\
\phicap &= \Be_2 e^{\Be_1 \Be_2 \phi} \\
\end{aligned}
\end{equation}

Taking \( \phi \) derivatives of \ref{eqn:laplacianCylindrical:100}, we have

\begin{equation}\label{eqn:laplacianCylindrical:120}
\begin{aligned}
\partial_\phi \rhocap &= \Be_1 \Be_1 \Be_2 e^{\Be_1 \Be_2 \phi} = \Be_2 e^{\Be_1 \Be_2 \phi} = \phicap \\
\partial_\phi \phicap &= \Be_2 \Be_1 \Be_2 e^{\Be_1 \Be_2 \phi} = -\Be_1 e^{\Be_1 \Be_2 \phi} = -\rhocap.
\end{aligned}
\end{equation}

The gradient of a vector \( \BA = \rhocap A_\rho + \phicap A_\phi + \zcap A_z \) is

\begin{equation}\label{eqn:laplacianCylindrical:60}
\begin{aligned}
\spacegrad \BA
&=
\lr{ \rhocap \partial_\rho + \frac{\phicap}{\rho} \partial_\phi + \zcap \partial_z }
\lr{ \rhocap A_\rho + \phicap A_\phi + \zcap A_z } \\
&=
\quad \rhocap \partial_\rho \lr{ \rhocap A_\rho + \phicap A_\phi + \zcap A_z } \\
&\quad + \frac{\phicap}{\rho} \partial_\phi \lr{ \rhocap A_\rho + \phicap A_\phi + \zcap A_z } \\
&\quad + \zcap \partial_z \lr{ \rhocap A_\rho + \phicap A_\phi + \zcap A_z } \\
&=
\quad \rhocap \lr{ \rhocap \partial_\rho A_\rho + \phicap \partial_\rho A_\phi + \zcap \partial_\rho A_z } \\
&\quad + \frac{\phicap}{\rho} \lr{ \partial_\phi(\rhocap A_\rho) + \partial_\phi(\phicap A_\phi) + \zcap \partial_\phi A_z } \\
&\quad + \zcap \lr{ \rhocap \partial_z A_\rho + \phicap \partial_z A_\phi + \zcap \partial_z A_z } \\
&=
\quad \partial_\rho A_\rho + \rhocap \phicap \partial_\rho A_\phi + \rhocap \zcap \partial_\rho A_z \\
&\quad +\frac{1}{\rho} \lr{ A_\rho + \phicap \rhocap \partial_\phi A_\rho – \phicap \rhocap A_\phi + \partial_\phi A_\phi + \phicap \zcap \partial_\phi A_z } \\
&\quad + \zcap \rhocap \partial_z A_\rho + \zcap \phicap \partial_z A_\phi + \partial_z A_z \\
&=
\quad \partial_\rho A_\rho + \frac{1}{\rho} \lr{ A_\rho + \partial_\phi A_\phi } + \partial_z A_z \\
&\quad +
\zcap \rhocap \lr{
\partial_z A_\rho
-\partial_\rho A_z
} \\
&\quad +
\phicap \zcap \lr{
\inv{\rho} \partial_\phi A_z
– \partial_z A_\phi
} \\
&\quad +
\rhocap \phicap \lr{
\partial_\rho A_\phi
– \inv{\rho} \lr{ \partial_\phi A_\rho – A_\phi }
},
\end{aligned}
\end{equation}

As expected, we see that the gradient splits nicely into a dot and curl

\begin{equation}\label{eqn:laplacianCylindrical:160}
\begin{aligned}
\spacegrad \BA
&= \spacegrad \cdot \BA + \spacegrad \wedge \BA \\
&= \spacegrad \cdot \BA + \rhocap \phicap \zcap (\spacegrad \cross \BA ),
\end{aligned}
\end{equation}

where the cylindrical representation of the divergence is seen to be

\begin{equation}\label{eqn:laplacianCylindrical:140}
\spacegrad \cdot \BA
=
\inv{\rho} \partial_\rho (\rho A_\rho) + \frac{1}{\rho} \partial_\phi A_\phi + \partial_z A_z,
\end{equation}

and the cylindrical representation of the curl is

\begin{equation}\label{eqn:laplacianCylindrical:180}
\spacegrad \cross \BA
=
\rhocap
\lr{
\inv{\rho} \partial_\phi A_z
– \partial_z A_\phi
}
+
\phicap
\lr{
\partial_z A_\rho
-\partial_\rho A_z
}
+
\inv{\rho} \zcap \lr{
\partial_\rho ( \rho A_\phi )
– \partial_\phi A_\rho
}.
\end{equation}

Should we want to, it is now possible to evaluate the Laplacian of \( \BA \) using
\ref{eqn:laplacianCylindrical:20}
, which will have the following components

\begin{equation}\label{eqn:laplacianCylindrical:220}
\begin{aligned}
\rhocap \cdot \lr{ \spacegrad^2 \BA }
&=
\partial_\rho
\lr{
\inv{\rho} \partial_\rho (\rho A_\rho) + \frac{1}{\rho} \partial_\phi A_\phi + \partial_z A_z
}
–
\lr{
\inv{\rho} \partial_\phi \lr{
\inv{\rho} \lr{
\partial_\rho ( \rho A_\phi ) – \partial_\phi A_\rho
}
}
– \partial_z \lr{
\partial_z A_\rho -\partial_\rho A_z
}
} \\
&=
\partial_\rho \lr{ \inv{\rho} \partial_\rho (\rho A_\rho)}
+ \partial_\rho \lr{ \frac{1}{\rho} \partial_\phi A_\phi}
+ \partial_{\rho z} A_z
– \inv{\rho^2}\partial_{\phi \rho} ( \rho A_\phi )
+ \inv{\rho^2}\partial_{\phi\phi} A_\rho
+ \partial_{zz} A_\rho
– \partial_{z\rho} A_z \\
&=
\partial_\rho \lr{ \inv{\rho} \partial_\rho (\rho A_\rho)}
+ \inv{\rho^2}\partial_{\phi\phi} A_\rho
+ \partial_{zz} A_\rho
– \frac{1}{\rho^2} \partial_\phi A_\phi
+ \frac{1}{\rho} \partial_{\rho\phi} A_\phi
– \inv{\rho^2}\partial_{\phi} A_\phi
– \inv{\rho}\partial_{\phi\rho} A_\phi \\
&=
\partial_\rho \lr{ \inv{\rho} \partial_\rho (\rho A_\rho)}
+ \inv{\rho^2}\partial_{\phi\phi} A_\rho
+ \partial_{zz} A_\rho
– \frac{2}{\rho^2} \partial_\phi A_\phi \\
&=
\inv{\rho} \partial_\rho \lr{ \rho \partial_\rho A_\rho}
+ \inv{\rho^2}\partial_{\phi\phi} A_\rho
+ \partial_{zz} A_\rho
– \frac{A_\rho}{\rho^2}
– \frac{2}{\rho^2} \partial_\phi A_\phi,
\end{aligned}
\end{equation}

\begin{equation}\label{eqn:laplacianCylindrical:240}
\begin{aligned}
\phicap \cdot \lr{ \spacegrad^2 \BA }
&=
\inv{\rho} \partial_\phi
\lr{
\inv{\rho} \partial_\rho (\rho A_\rho) + \frac{1}{\rho} \partial_\phi A_\phi + \partial_z A_z
}
–
\lr{
\lr{
\partial_z \lr{
\inv{\rho} \partial_\phi A_z – \partial_z A_\phi
}
-\partial_\rho \lr{
\inv{\rho} \lr{ \partial_\rho ( \rho A_\phi ) – \partial_\phi A_\rho}
}
}
} \\
&=
\inv{\rho^2} \partial_{\phi\rho} (\rho A_\rho)
+ \frac{1}{\rho^2} \partial_{\phi\phi} A_\phi
+ \inv{\rho}\partial_{\phi z} A_z
– \inv{\rho} \partial_{z\phi} A_z
+ \partial_{z z} A_\phi
+\partial_\rho \lr{ \inv{\rho} \partial_\rho ( \rho A_\phi ) }
– \partial_\rho \lr{ \inv{\rho} \partial_\phi A_\rho} \\
&=
\partial_\rho \lr{ \inv{\rho} \partial_\rho ( \rho A_\phi ) }
+ \frac{1}{\rho^2} \partial_{\phi\phi} A_\phi
+ \partial_{z z} A_\phi
+ \inv{\rho^2} \partial_{\phi\rho} (\rho A_\rho)
+ \inv{\rho}\partial_{\phi z} A_z
– \inv{\rho} \partial_{z\phi} A_z
– \partial_\rho \lr{ \inv{\rho} \partial_\phi A_\rho} \\
&=
\partial_\rho \lr{ \inv{\rho} \partial_\rho ( \rho A_\phi ) }
+ \frac{1}{\rho^2} \partial_{\phi\phi} A_\phi
+ \partial_{z z} A_\phi
+ \inv{\rho^2} \partial_{\phi} A_\rho
+ \inv{\rho} \partial_{\phi\rho} A_\rho
+ \inv{\rho^2} \partial_\phi A_\rho
– \inv{\rho} \partial_{\rho\phi} A_\rho \\
&=
\partial_\rho \lr{ \inv{\rho} \partial_\rho ( \rho A_\phi ) }
+ \frac{1}{\rho^2} \partial_{\phi\phi} A_\phi
+ \partial_{z z} A_\phi
+ \frac{2}{\rho^2} \partial_{\phi} A_\rho \\
&=
\inv{\rho} \partial_\rho \lr{ \rho \partial_\rho A_\phi }
+ \frac{1}{\rho^2} \partial_{\phi\phi} A_\phi
+ \partial_{z z} A_\phi
+ \frac{2}{\rho^2} \partial_{\phi} A_\rho
– \frac{A_\phi}{\rho^2},
\end{aligned}
\end{equation}

\begin{equation}\label{eqn:laplacianCylindrical:260}
\begin{aligned}
\zcap \cdot \lr{ \spacegrad^2 \BA }
&=
\partial_z
\lr{
\inv{\rho} \partial_\rho (\rho A_\rho) + \frac{1}{\rho} \partial_\phi A_\phi + \partial_z A_z
}
–
\inv{\rho} \lr{
\partial_\rho \lr{ \rho \lr{
\partial_z A_\rho -\partial_\rho A_z
}
}
– \partial_\phi \lr{
\inv{\rho} \partial_\phi A_z – \partial_z A_\phi
}
} \\
&=
\inv{\rho} \partial_{z\rho} (\rho A_\rho)
+ \frac{1}{\rho} \partial_{z\phi} A_\phi
+ \partial_{zz} A_z
– \inv{\rho}\partial_\rho \lr{ \rho \partial_z A_\rho }
+ \inv{\rho}\partial_\rho \lr{ \rho \partial_\rho A_z }
+ \inv{\rho^2} \partial_{\phi\phi} A_z
– \inv{\rho} \partial_{\phi z} A_\phi \\
&=
\inv{\rho}\partial_\rho \lr{ \rho \partial_\rho A_z }
+ \inv{\rho^2} \partial_{\phi\phi} A_z
+ \partial_{zz} A_z
+ \inv{\rho} \partial_{z} A_\rho
+\partial_{z\rho} A_\rho
+ \frac{1}{\rho} \partial_{z\phi} A_\phi
– \inv{\rho}\partial_z A_\rho
– \partial_{\rho z} A_\rho
– \inv{\rho} \partial_{\phi z} A_\phi \\
&=
\inv{\rho}\partial_\rho \lr{ \rho \partial_\rho A_z }
+ \inv{\rho^2} \partial_{\phi\phi} A_z
+ \partial_{zz} A_z
\end{aligned}
\end{equation}

Evaluating these was a fairly tedious and mechanical job, and would have been better suited to a computer algebra system than by hand as done here.

Explicit cylindrical Laplacian

Let’s try this a different way. The most obvious potential strategy is to just apply the Laplacian to the vector itself, but we need to include the unit vectors in such an operation

\begin{equation}\label{eqn:laplacianCylindrical:280}
\spacegrad^2 \BA =
\spacegrad^2 \lr{ \rhocap A_\rho + \phicap A_\phi + \zcap A_z }.
\end{equation}

First we need to know the explicit form of the cylindrical Laplacian. From the painful expansion, we can guess that it is

\begin{equation}\label{eqn:laplacianCylindrical:300}
\spacegrad^2 \psi
=
\inv{\rho}\partial_\rho \lr{ \rho \partial_\rho \psi }
+ \inv{\rho^2} \partial_{\phi\phi} \psi
+ \partial_{zz} \psi.
\end{equation}

Let’s check that explicitly. Here I use the vector product where \( \rhocap^2 = \phicap^2 = \zcap^2 = 1 \), and these vectors anticommute when different

\begin{equation}\label{eqn:laplacianCylindrical:320}
\begin{aligned}
\spacegrad^2 \psi
&=
\lr{ \rhocap \partial_\rho + \frac{\phicap}{\rho} \partial_\phi + \zcap \partial_z }
\lr{ \rhocap \partial_\rho \psi + \frac{\phicap}{\rho} \partial_\phi \psi + \zcap \partial_z \psi } \\
&=
\rhocap \partial_\rho
\lr{ \rhocap \partial_\rho \psi + \frac{\phicap}{\rho} \partial_\phi \psi + \zcap \partial_z \psi }
+ \frac{\phicap}{\rho} \partial_\phi
\lr{ \rhocap \partial_\rho \psi + \frac{\phicap}{\rho} \partial_\phi \psi + \zcap \partial_z \psi }
+ \zcap \partial_z
\lr{ \rhocap \partial_\rho \psi + \frac{\phicap}{\rho} \partial_\phi \psi + \zcap \partial_z \psi } \\
&=
\partial_{\rho\rho} \psi
+ \rhocap \phicap \partial_\rho \lr{ \frac{1}{\rho} \partial_\phi \psi}
+ \rhocap \zcap \partial_{\rho z} \psi
+ \frac{\phicap}{\rho} \partial_\phi \lr{ \rhocap \partial_\rho \psi }
+ \frac{\phicap}{\rho} \partial_\phi \lr{ \frac{\phicap}{\rho} \partial_\phi \psi }
+ \frac{\phicap \zcap }{\rho} \partial_{\phi z} \psi
+ \zcap \rhocap \partial_{z\rho} \psi
+ \frac{\zcap \phicap}{\rho} \partial_{z\phi} \psi
+ \partial_{zz} \psi \\
&=
\partial_{\rho\rho} \psi
+ \inv{\rho} \partial_\rho \psi
+ \frac{1}{\rho^2} \partial_{\phi \phi} \psi
+ \partial_{zz} \psi
+ \rhocap \phicap
\lr{
-\frac{1}{\rho^2} \partial_\phi \psi
+\frac{1}{\rho} \partial_{\rho \phi} \psi
-\inv{\rho} \partial_{\phi \rho} \psi
+ \frac{1}{\rho^2} \partial_\phi \psi
}
+ \zcap \rhocap \lr{
-\partial_{\rho z} \psi
+ \partial_{z\rho} \psi
}
+ \phicap \zcap \lr{
\inv{\rho} \partial_{\phi z} \psi
– \inv{\rho} \partial_{z\phi} \psi
} \\
&=
\partial_{\rho\rho} \psi
+ \inv{\rho} \partial_\rho \psi
+ \frac{1}{\rho^2} \partial_{\phi \phi} \psi
+ \partial_{zz} \psi,
\end{aligned}
\end{equation}

so the Laplacian operator is

\begin{equation}\label{eqn:laplacianCylindrical:340}
\boxed{
\spacegrad^2
=
\inv{\rho} \PD{\rho}{} \lr{ \rho \PD{\rho}{} }
+ \frac{1}{\rho^2} \PDSq{\phi}{}
+ \PDSq{z}{}.
}
\end{equation}

All the bivector grades of the Laplacian operator are seen to explicitly cancel, regardless of the grade of \( \psi \), just as if we had expanded the scalar Laplacian as a dot product
\( \spacegrad^2 \psi = \spacegrad \cdot \lr{ \spacegrad \psi} \).
Unlike such a scalar expansion, this derivation is seen to be valid for any grade \( \psi \). We know now that we can trust this result when \( \psi \) is a scalar, a vector, a bivector, a trivector, or even a multivector.

Vector Laplacian

Now that we trust that the typical scalar form of the Laplacian applies equally well to multivectors as it does to scalars, that cylindrical coordinate operator can now be applied to a
vector. Consider the projections onto each of the directions in turn

\begin{equation}\label{eqn:laplacianCylindrical:360}
\spacegrad^2 \lr{ \rhocap A_\rho }
=
\rhocap \inv{\rho} \partial_\rho \lr{ \rho \partial_\rho A_\rho }
+ \frac{1}{\rho^2} \partial_{\phi\phi} \lr{\rhocap A_\rho}
+ \rhocap \partial_{zz} A_\rho
\end{equation}

\begin{equation}\label{eqn:laplacianCylindrical:380}
\begin{aligned}
\partial_{\phi\phi} \lr{\rhocap A_\rho}
&=
\partial_\phi \lr{ \phicap A_\rho + \rhocap \partial_\phi A_\rho } \\
&=
-\rhocap A_\rho
+\phicap \partial_\phi A_\rho
+ \phicap \partial_\phi A_\rho
+ \rhocap \partial_{\phi\phi} A_\rho \\
&=
\rhocap \lr{ \partial_{\phi\phi} A_\rho -A_\rho }
+ 2 \phicap \partial_\phi A_\rho
\end{aligned}
\end{equation}

so this component of the vector Laplacian is

\begin{equation}\label{eqn:laplacianCylindrical:400}
\begin{aligned}
\spacegrad^2 \lr{ \rhocap A_\rho }
&=
\rhocap
\lr{
\inv{\rho} \partial_\rho \lr{ \rho \partial_\rho A_\rho }
+ \inv{\rho^2} \partial_{\phi\phi} A_\rho
– \inv{\rho^2} A_\rho
+ \partial_{zz} A_\rho
}
+
\phicap
\lr{
2 \inv{\rho^2} \partial_\phi A_\rho
} \\
&=
\rhocap \lr{
\spacegrad^2 A_\rho
– \inv{\rho^2} A_\rho
}
+
\phicap
\frac{2}{\rho^2} \partial_\phi A_\rho
.
\end{aligned}
\end{equation}

The Laplacian for the projection of the vector onto the \( \phicap \) direction is

\begin{equation}\label{eqn:laplacianCylindrical:420}
\spacegrad^2 \lr{ \phicap A_\phi }
=
\phicap \inv{\rho} \partial_\rho \lr{ \rho \partial_\rho A_\phi }
+ \frac{1}{\rho^2} \partial_{\phi\phi} \lr{\phicap A_\phi}
+ \phicap \partial_{zz} A_\phi,
\end{equation}

Again, since the unit vectors are \( \phi \) dependent, the \( \phi \) derivatives have to be treated carefully

\begin{equation}\label{eqn:laplacianCylindrical:440}
\begin{aligned}
\partial_{\phi\phi} \lr{\phicap A_\phi}
&=
\partial_{\phi} \lr{-\rhocap A_\phi + \phicap \partial_\phi A_\phi} \\
&=
-\phicap A_\phi
-\rhocap \partial_\phi A_\phi
– \rhocap \partial_\phi A_\phi
+ \phicap \partial_{\phi \phi} A_\phi \\
&=
– 2 \rhocap \partial_\phi A_\phi
+
\phicap
\lr{
\partial_{\phi \phi} A_\phi
– A_\phi
},
\end{aligned}
\end{equation}

so the Laplacian of this projection is
\begin{equation}\label{eqn:laplacianCylindrical:460}
\begin{aligned}
\spacegrad^2 \lr{ \phicap A_\phi }
&=
\phicap
\lr{
\inv{\rho} \partial_\rho \lr{ \rho \partial_\rho A_\phi }
+ \phicap \partial_{zz} A_\phi,
\inv{\rho^2} \partial_{\phi \phi} A_\phi
– \frac{A_\phi }{\rho^2}
}
– \rhocap \frac{2}{\rho^2} \partial_\phi A_\phi \\
&=
\phicap \lr{
\spacegrad^2 A_\phi
– \frac{A_\phi}{\rho^2}
}
– \rhocap \frac{2}{\rho^2} \partial_\phi A_\phi.
\end{aligned}
\end{equation}

Since \( \zcap \) is fixed we have

\begin{equation}\label{eqn:laplacianCylindrical:480}
\spacegrad^2 \zcap A_z
=
\zcap \spacegrad^2 A_z.
\end{equation}

Putting all the pieces together we have
\begin{equation}\label{eqn:laplacianCylindrical:500}
\boxed{
\spacegrad^2 \BA
=
\rhocap \lr{
\spacegrad^2 A_\rho
– \inv{\rho^2} A_\rho
– \frac{2}{\rho^2} \partial_\phi A_\phi
}
+\phicap \lr{
\spacegrad^2 A_\phi
– \frac{A_\phi}{\rho^2}
+ \frac{2}{\rho^2} \partial_\phi A_\rho
}
+
\zcap \spacegrad^2 A_z.
}
\end{equation}

This matches the results of \ref{eqn:laplacianCylindrical:220}, …, from the painful expansion of
\( \spacegrad \lr{ \spacegrad \cdot \BA } – \spacegrad \cross \lr{ \spacegrad \cross \BA } \).

Line charge field and potential.

October 26, 2016 math and physics play , , , , , , , , , ,

[Click here for a PDF of this post with nicer formatting]

When computing the most general solution of the electrostatic potential in a plane, Jackson [1] mentions that \( -2 \lambda_0 \ln \rho \) is the well known potential for an infinite line charge (up to the unit specific factor). Checking that statement, since I didn’t recall what that potential was offhand, I encountered some inconsistencies and non-convergent integrals, and thought it was worthwhile to explore those a bit more carefully. This will be done here.

Using Gauss’s law.

For an infinite length line charge, we can find the radial field contribution using Gauss’s law, imagining a cylinder of length \( \Delta l \) of radius \( \rho \) surrounding this charge with the midpoint at the origin. Ignoring any non-radial field contribution, we have

\begin{equation}\label{eqn:lineCharge:20}
\int_{-\Delta l/2}^{\Delta l/2} \ncap \cdot \BE (2 \pi \rho) dl = \frac{\lambda_0}{\epsilon_0} \Delta l,
\end{equation}

or

\begin{equation}\label{eqn:lineCharge:40}
\BE = \frac{\lambda_0}{2 \pi \epsilon_0} \frac{\rhocap}{\rho}.
\end{equation}

Since

\begin{equation}\label{eqn:lineCharge:60}
\frac{\rhocap}{\rho} = \spacegrad \ln \rho,
\end{equation}

this means that the potential is

\begin{equation}\label{eqn:lineCharge:80}
\phi = -\frac{2 \lambda_0}{4 \pi \epsilon_0} \ln \rho.
\end{equation}

Finite line charge potential.

Let’s try both these calculations for a finite charge distribution. Gauss’s law looses its usefulness, but we can evaluate the integrals directly. For the electric field

\begin{equation}\label{eqn:lineCharge:100}
\BE
= \frac{\lambda_0}{4 \pi \epsilon_0} \int \frac{(\Bx – \Bx’)}{\Abs{\Bx – \Bx’}^3} dl’.
\end{equation}

Using cylindrical coordinates with the field point \( \Bx = \rho \rhocap \) for convience, the charge point \( \Bx’ = z’ \zcap \), and a the charge distributed over \( [a,b] \) this is

\begin{equation}\label{eqn:lineCharge:120}
\BE
= \frac{\lambda_0}{4 \pi \epsilon_0} \int_a^b \frac{(\rho \rhocap – z’ \zcap)}{\lr{\rho^2 + (z’)^2}^{3/2}} dz’.
\end{equation}

When the charge is uniformly distributed around the origin \( [a,b] = b[-1,1] \) the \( \zcap \) component of this field is killed because the integrand is odd. This justifies ignoring such contributions in the Gaussing cylinder analysis above. The general solution to this integral is found to be

\begin{equation}\label{eqn:lineCharge:140}
\BE
=
\frac{\lambda_0}{4 \pi \epsilon_0}
\evalrange{
\lr{
\frac{z’ \rhocap }{\rho \sqrt{ \rho^2 + (z’)^2 } }
+\frac{\zcap}{ \sqrt{ \rho^2 + (z’)^2 } }
}
}{a}{b},
\end{equation}

or
\begin{equation}\label{eqn:lineCharge:240}
\boxed{
\BE
=
\frac{\lambda_0}{4 \pi \epsilon_0}
\lr{
\frac{\rhocap }{\rho}
\lr{
\frac{b}{\sqrt{ \rho^2 + b^2 } }
-\frac{a}{\sqrt{ \rho^2 + a^2 } }
}
+ \zcap
\lr{
\frac{1}{ \sqrt{ \rho^2 + b^2 } }
-\frac{1}{ \sqrt{ \rho^2 + a^2 } }
}
}.
}
\end{equation}

When \( b = -a = \Delta l/2 \), this reduces to

\begin{equation}\label{eqn:lineCharge:160}
\BE
=
\frac{\lambda_0}{4 \pi \epsilon_0}
\frac{\rhocap }{\rho}
\frac{\Delta l}{\sqrt{ \rho^2 + (\Delta l/2)^2 } },
\end{equation}

which further reduces to \ref{eqn:lineCharge:40} when \( \Delta l \gg \rho \).

Finite line charge potential. Wrong but illuminating.

Again, putting the field point at \( z’ = 0 \), we have

\begin{equation}\label{eqn:lineCharge:180}
\phi(\rho)
= \frac{\lambda_0}{4 \pi \epsilon_0} \int_a^b \frac{dz’}{\lr{\rho^2 + (z’)^2}^{1/2}},
\end{equation}

which integrates to
\begin{equation}\label{eqn:lineCharge:260}
\phi(\rho)
= \frac{\lambda_0}{4 \pi \epsilon_0 }
\ln \frac{ b + \sqrt{ \rho^2 + b^2 }}{ a + \sqrt{\rho^2 + a^2}}.
\end{equation}

With \( b = -a = \Delta l/2 \), this approaches

\begin{equation}\label{eqn:lineCharge:200}
\phi
\approx
\frac{\lambda_0}{4 \pi \epsilon_0 }
\ln \frac{ (\Delta l/2) }{ \rho^2/2\Abs{\Delta l/2}}
=
\frac{-2 \lambda_0}{4 \pi \epsilon_0 } \ln \rho
+
\frac{\lambda_0}{4 \pi \epsilon_0 }
\ln \lr{ (\Delta l)^2/2 }.
\end{equation}

Before \( \Delta l \) is allowed to tend to infinity, this is identical (up to a difference in the reference potential) to \ref{eqn:lineCharge:80} found using Gauss’s law. It is, strictly speaking, singular when \( \Delta l \rightarrow \infty \), so it does not seem right to infinity as a reference point for the potential.

There’s another weird thing about this result. Since this has no \( z \) dependence, it is not obvious how we would recover the non-radial portion of the electric field from this potential using \( \BE = -\spacegrad \phi \)? Let’s calculate the elecric field from \ref{eqn:lineCharge:180} explicitly

\begin{equation}\label{eqn:lineCharge:220}
\begin{aligned}
\BE
&=
-\frac{\lambda_0}{4 \pi \epsilon_0}
\spacegrad
\ln \frac{ b + \sqrt{ \rho^2 + b^2 }}{ a + \sqrt{\rho^2 + a^2}} \\
&=
-\frac{\lambda_0 \rhocap}{4 \pi \epsilon_0 }
\PD{\rho}{}
\ln \frac{ b + \sqrt{ \rho^2 + b^2 }}{ a + \sqrt{\rho^2 + a^2}} \\
&=
-\frac{\lambda_0 \rhocap}{4 \pi \epsilon_0}
\lr{
\inv{ b + \sqrt{ \rho^2 + b^2 }} \frac{ \rho }{\sqrt{ \rho^2 + b^2 }}
-\inv{ a + \sqrt{ \rho^2 + a^2 }} \frac{ \rho }{\sqrt{ \rho^2 + a^2 }}
} \\
&=
-\frac{\lambda_0 \rhocap}{4 \pi \epsilon_0 \rho}
\lr{
\frac{ -b + \sqrt{ \rho^2 + b^2 }}{\sqrt{ \rho^2 + b^2 }}
-\frac{ -a + \sqrt{ \rho^2 + a^2 }}{\sqrt{ \rho^2 + a^2 }}
} \\
&=
\frac{\lambda_0 \rhocap}{4 \pi \epsilon_0 \rho}
\lr{
\frac{ b }{\sqrt{ \rho^2 + b^2 }}
-\frac{ a }{\sqrt{ \rho^2 + a^2 }}
}.
\end{aligned}
\end{equation}

This recovers the radial component of the field from \ref{eqn:lineCharge:240}, but where did the \( \zcap \) component go? The required potential appears to be

\begin{equation}\label{eqn:lineCharge:340}
\phi(\rho, z)
=
\frac{\lambda_0}{4 \pi \epsilon_0 }
\ln \frac{ b + \sqrt{ \rho^2 + b^2 }}{ a + \sqrt{\rho^2 + a^2}}
–
\frac{z \lambda_0}{4 \pi \epsilon_0 }
\lr{ \frac{1}{\sqrt{\rho^2 + b^2}}
-\frac{1}{\sqrt{\rho^2 + a^2}}
}.
\end{equation}

When computing the electric field \( \BE(\rho, \theta, z) \), it was convienent to pick the coordinate system so that \( z = 0 \). Doing this with the potential gives the wrong answers. The reason for this appears to be that this kills the potential term that is linear in \( z \) before taking its gradient, and we need that term to have the \( \zcap \) field component that is expected for a charge distribution that is non-symmetric about the origin on the z-axis!

Finite line charge potential. Take II.

Let the point at which the potential is evaluated be

\begin{equation}\label{eqn:lineCharge:360}
\Bx = \rho \rhocap + z \zcap,
\end{equation}

and the charge point be
\begin{equation}\label{eqn:lineCharge:380}
\Bx’ = z’ \zcap.
\end{equation}

This gives

\begin{equation}\label{eqn:lineCharge:400}
\begin{aligned}
\phi(\rho, z)
&= \frac{\lambda_0}{4\pi \epsilon_0} \int_a^b \frac{dz’}{\Abs{\rho^2 + (z – z’)^2 }} \\
&= \frac{\lambda_0}{4\pi \epsilon_0} \int_{a-z}^{b-z} \frac{du}{ \Abs{\rho^2 + u^2} } \\
&= \frac{\lambda_0}{4\pi \epsilon_0}
\evalrange{\ln \lr{ u + \sqrt{ \rho^2 + u^2 }}}{b-z}{a-z} \\
&=
\frac{\lambda_0}{4\pi \epsilon_0}
\ln \frac
{ b-z + \sqrt{ \rho^2 + (b-z)^2 }}
{ a-z + \sqrt{ \rho^2 + (a-z)^2 }}.
\end{aligned}
\end{equation}

The limit of this potential \( a = -\Delta/2 \rightarrow -\infty, b = \Delta/2 \rightarrow \infty \) doesn’t exist in any strict sense. If we are cavilier about the limits, as in \ref{eqn:lineCharge:200}, this can be evaluated as

\begin{equation}\label{eqn:lineCharge:n}
\phi \approx
\frac{\lambda_0}{4\pi \epsilon_0} \lr{ -2 \ln \rho + \textrm{constant} }.
\end{equation}

however, the constant (\( \ln \Delta^2/2 \)) is infinite, so there isn’t really a good justification for using that constant as the potential reference point directly.

It seems that the “right” way to calculate the potential for the infinite distribution, is to

  • Calculate the field from the potential.
  • Take the PV limit of that field with the charge distribution extending to infinity.
  • Compute the corresponding potential from this limiting value of the field.

Doing that doesn’t blow up. That field calculation, for the finite case, should include a \( \zcap \) component. To verify, let’s take the respective derivatives

\begin{equation}\label{eqn:lineCharge:420}
\begin{aligned}
-\PD{z}{} \phi
&=
-\frac{\lambda_0}{4\pi \epsilon_0}
\lr{
\frac{ -1 + \frac{z – b}{\sqrt{ \rho^2 + (b-z)^2 }} }{
b-z + \sqrt{ \rho^2 + (b-z)^2 }
}
–
\frac{ -1 + \frac{z – a}{\sqrt{ \rho^2 + (a-z)^2 }} }{
a-z + \sqrt{ \rho^2 + (a-z)^2 }
}
} \\
&=
\frac{\lambda_0}{4\pi \epsilon_0}
\lr{
\frac{ 1 + \frac{b – z}{\sqrt{ \rho^2 + (b-z)^2 }} }{
b-z + \sqrt{ \rho^2 + (b-z)^2 }
}
–
\frac{ 1 + \frac{a – z}{\sqrt{ \rho^2 + (a-z)^2 }} }{
a-z + \sqrt{ \rho^2 + (a-z)^2 }
}
} \\
&=
\frac{\lambda_0}{4\pi \epsilon_0}
\lr{
\inv{\sqrt{ \rho^2 + (b-z)^2 }}
-\inv{\sqrt{ \rho^2 + (a-z)^2 }}
},
\end{aligned}
\end{equation}

and

\begin{equation}\label{eqn:lineCharge:440}
\begin{aligned}
-\PD{\rho}{} \phi
&=
-\frac{\lambda_0}{4\pi \epsilon_0}
\lr{
\frac{ \frac{\rho}{\sqrt{ \rho^2 + (b-z)^2 }} }{
b-z + \sqrt{ \rho^2 + (b-z)^2 }
}
–
\frac{ \frac{\rho}{\sqrt{ \rho^2 + (a-z)^2 }} }{
a-z + \sqrt{ \rho^2 + (a-z)^2 }
}
} \\
&=
-\frac{\lambda_0}{4\pi \epsilon_0}
\lr{
\frac{\rho \lr{
-(b-z) + \sqrt{ \rho^2 + (b-z)^2 }
}}{ \rho^2 \sqrt{ \rho^2 + (b-z)^2 } }
–
\frac{\rho \lr{
-(a-z) + \sqrt{ \rho^2 + (a-z)^2 }
}}{ \rho^2 \sqrt{ \rho^2 + (a-z)^2 } }
} \\
&=
\frac{\lambda_0}{4\pi \epsilon_0 \rho}
\lr{
\frac{b-z}{\sqrt{ \rho^2 + (b-z)^2 }}
-\frac{a-z}{\sqrt{ \rho^2 + (a-z)^2 }}
}
.
\end{aligned}
\end{equation}

Putting the pieces together, the electric field is
\begin{equation}\label{eqn:lineCharge:460}
\BE =
\frac{\lambda_0}{4\pi \epsilon_0}
\lr{
\frac{\rhocap}{\rho} \lr{
\frac{b-z}{\sqrt{ \rho^2 + (b-z)^2 }}
-\frac{a-z}{\sqrt{ \rho^2 + (a-z)^2 }}
}
+
\zcap \lr{
\inv{\sqrt{ \rho^2 + (b-z)^2 }}
-\inv{\sqrt{ \rho^2 + (a-z)^2 }}
}
}.
\end{equation}

For has a PV limit of \ref{eqn:lineCharge:40} at \( z = 0 \), and also for the finite case, has the \( \zcap \) field component that was obtained when the field was obtained by direct integration.

Conclusions

  • We have to evaluate the potential at all points in space, not just on the axis that we evaluate the field on (should we choose to do so).
  • In this case, we found that it was not directly meaningful to take the limit of a potential distribution. We can, however, compute the field from a potential for a finite charge distribution,
    take the limit of that field, and then calculate the corresponding potential for the infinite distribution.

Is there a more robust mechanism that can be used to directly calculate the potential for an infinite charge distribution, instead of calculating the potential from the field of such an infinite distribution?

I think that were things go wrong is that the integral of \ref{eqn:lineCharge:180} does not apply to charge distributions that are not finite on the infinite range \( z \in [-\infty, \infty] \). That solution was obtained by utilizing an all-space Green’s function, and the boundary term in that Green’s analysis was assumed to tend to zero. That isn’t the case when the charge distribution is \( \lambda_0 \delta( z ) \).

References

[1] JD Jackson. Classical Electrodynamics. John Wiley and Sons, 2nd edition, 1975.

Jackson’s electrostatic self energy analysis

October 10, 2016 math and physics play , , , , , , , , , ,

[Click here for a PDF of this post with nicer formatting]

Motivation

I was reading my Jackson [1], which characteristically had the statement “the […] integral can easily be shown to have the value \( 4 \pi \)”, in a discussion of electrostatic energy and self energy. After a few attempts and a couple of pages of calculations, I figured out how this can be easily shown.

Context

Let me walk through the context that leads to the “easy” integral, and then the evaluation of that integral. Unlike my older copy of Jackson, I’ll do this in SI units.

The starting point is a statement that the work done (potential energy) of one charge \( q_i \) in a set of \( n \) charges, where that charge is brought to its position \( \Bx_i \) from infinity, is

\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:20}
W_i = q_i \Phi(\Bx_i),
\end{equation}

where the potential energy due to the rest of the charge configuration is

\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:40}
\Phi(\Bx_i) = \inv{4 \pi \epsilon} \sum_{i \ne j} \frac{q_j}{\Abs{\Bx_i – \Bx_j}}.
\end{equation}

This means that the total potential energy, making sure not to double count, to move all the charges in from infinity is

\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:60}
W = \inv{4 \pi \epsilon} \sum_{1 \le i < j \le n} \frac{q_i q_j}{\Abs{\Bx_i - \Bx_j}}. \end{equation} This sum over all unique pairs is somewhat unwieldy, so it can be adjusted by explicitly double counting with a corresponding divide by two \begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:80} W = \inv{2} \inv{4 \pi \epsilon} \sum_{1 \le i \ne j \le n} \frac{q_i q_j}{\Abs{\Bx_i - \Bx_j}}. \end{equation} The point that causes the trouble later is the continuum equivalent to this relationship, which is \begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:100} W = \inv{8 \pi \epsilon} \int \frac{\rho(\Bx) \rho(\Bx')}{\Abs{\Bx - \Bx'}} d^3 \Bx d^3 \Bx', \end{equation} or \begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:120} W = \inv{2} \int \rho(\Bx) \Phi(\Bx) d^3 \Bx. \end{equation} There's a subtlety here that is often passed over. When the charge densities represent point charges \( \rho(\Bx) = q \delta^3(\Bx - \Bx') \) are located at, notice that this integral equivalent is evaluated over all space, including the spaces that the charges that the charges are located at. Ignoring that subtlety, this potential energy can be expressed in terms of the electric field, and then integrated by parts \begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:140} \begin{aligned} W &= \inv{2 } \int (\spacegrad \cdot (\epsilon \BE)) \Phi(\Bx) d^3 \Bx \\ &= \frac{\epsilon}{2 } \int \lr{ \spacegrad \cdot (\BE \Phi) - (\spacegrad \Phi) \cdot \BE } d^3 \Bx \\ &= \frac{\epsilon}{2 } \oint dA \ncap \cdot (\BE \Phi) + \frac{\epsilon}{2 } \int \BE \cdot \BE d^3 \Bx. \end{aligned} \end{equation} The presumption is that \( \BE \Phi \) falls off as the bounds of the integration volume tends to infinity. That leaves us with an energy density proportional to the square of the field \begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:160} w = \frac{\epsilon}{2 } \BE^2. \end{equation}

Inconsistency

It’s here that Jackson points out the inconsistency between \ref{eqn:electrostaticJacksonSelfEnergy:160} and the original
discrete analogue \ref{eqn:electrostaticJacksonSelfEnergy:80} that this was based on. The energy density is positive definite, whereas the discrete potential energy can be negative if there is a difference in the sign of the charges.

Here Jackson uses a two particle charge distribution to help resolve this conundrum. For a superposition \( \BE = \BE_1 + \BE_2 \), we have

\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:180}
\BE
=
\inv{4 \pi \epsilon} \frac{q_1 (\Bx – \Bx_1)}{\Abs{\Bx – \Bx_1}^3}
+ \inv{4 \pi \epsilon} \frac{q_2 (\Bx – \Bx_2)}{\Abs{\Bx – \Bx_2}^3},
\end{equation}

so the energy density is
\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:200}
w =
\frac{1}{32 \pi^2 \epsilon} \frac{q_1^2}{\Abs{\Bx – \Bx_1}^4 }
+
\frac{1}{32 \pi^2 \epsilon} \frac{q_2^2}{\Abs{\Bx – \Bx_2}^4 }
+
2 \frac{q_1 q_2}{32 \pi^2 \epsilon}
\frac{(\Bx – \Bx_1)}{\Abs{\Bx – \Bx_1}^3} \cdot
\frac{(\Bx – \Bx_2)}{\Abs{\Bx – \Bx_2}^3}.
\end{equation}

The discrete potential had only an interaction energy, whereas the potential from this squared field has an interaction energy plus two self energy terms. Those two strictly positive self energy terms are what forces this field energy positive, independent of the sign of the interaction energy density. Jackson makes a change of variables of the form

\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:220}
\begin{aligned}
\Brho &= (\Bx – \Bx_1)/R \\
R &= \Abs{\Bx_1 – \Bx_2} \\
\ncap &= (\Bx_1 – \Bx_2)/R,
\end{aligned}
\end{equation}

for which we find

\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:240}
\Bx = \Bx_1 + R \Brho,
\end{equation}

so
\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:260}
\Bx – \Bx_2 =
\Bx_1 – \Bx_2 + R \Brho
R (\ncap + \Brho),
\end{equation}

and
\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:280}
d^3 \Bx = R^3 d^3 \Brho,
\end{equation}

so the total interaction energy is
\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:300}
\begin{aligned}
W_{\textrm{int}}
&=
\frac{q_1 q_2}{16 \pi^2 \epsilon}
\int d^3 \Bx
\frac{(\Bx – \Bx_1)}{\Abs{\Bx – \Bx_1}^3} \cdot
\frac{(\Bx – \Bx_2)}{\Abs{\Bx – \Bx_2}^3} \\
&=
\frac{q_1 q_2}{16 \pi^2 \epsilon}
\int R^3 d^3 \Brho
\frac{ R \Brho }{ R^3 \Abs{\Brho}^3 } \cdot
\frac{R (\ncap + \Brho)}{R^3 \Abs{\ncap + \Brho}^3} \\
&=
\frac{q_1 q_2}{16 \pi^2 \epsilon R}
\int d^3 \Brho
\frac{ \Brho }{ \Abs{\Brho}^3 } \cdot
\frac{(\ncap + \Brho)}{ \Abs{\ncap + \Brho}^3}.
\end{aligned}
\end{equation}

Evaluating this integral is what Jackson calls easy. The technique required is to express the integrand in terms of gradients in the \( \Brho \) coordinate system

\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:320}
\begin{aligned}
\int d^3 \Brho
\frac{ \Brho }{ \Abs{\Brho}^3 } \cdot
\frac{(\ncap + \Brho)}{ \Abs{\ncap + \Brho}^3}
&=
\int d^3 \Brho
\lr{ – \spacegrad_\Brho \inv{\Abs{\Brho}} }
\cdot
\lr{ – \spacegrad_\Brho \inv{\Abs{\ncap + \Brho}} } \\
&=
\int d^3 \Brho
\lr{ \spacegrad_\Brho \inv{\Abs{\Brho}} }
\cdot
\lr{ \spacegrad_\Brho \inv{\Abs{\ncap + \Brho}} }.
\end{aligned}
\end{equation}

I found it somewhat non-trivial to find the exact form of the chain rule that is required to simplify this integral, but after some trial and error, figured it out by working backwards from
\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:340}
\spacegrad_\Brho^2 \inv{ \Abs{\Brho} \Abs{\ncap + \Brho}}
=
\spacegrad_\Brho \cdot \lr{ \inv{\Abs{\Brho}} \spacegrad_\Brho \inv{ \Abs{\ncap + \Brho} } }
+
\spacegrad_\Brho \cdot \lr{ \inv{\Abs{\ncap + \Brho}} \spacegrad_\Brho \inv{ \Abs{\Brho} } }.
\end{equation}

In integral form this is
\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:360}
\begin{aligned}
\oint dA’ \ncap’ \cdot \spacegrad_\Brho \inv{ \Abs{\Brho} \Abs{\ncap + \Brho}}
&=
\int d^3 \Brho’
\spacegrad_{\Brho’} \cdot \lr{ \inv{\Abs{\Brho’ – \ncap}} \spacegrad_{\Brho’} \inv{ \Abs{\Brho’} } }
+
\int d^3 \Brho
\spacegrad_\Brho \cdot \lr{ \inv{\Abs{\ncap + \Brho}} \spacegrad_\Brho \inv{ \Abs{\Brho} } } \\
&=
\int d^3 \Brho’
\lr{ \spacegrad_{\Brho’} \inv{\Abs{\Brho’ – \ncap} } \cdot \spacegrad_{\Brho’} \inv{ \Abs{\Brho’} } }
+
\int d^3 \Brho’
\inv{\Abs{\Brho’ – \ncap}} \spacegrad_{\Brho’}^2 \inv{ \Abs{\Brho’} } \\
&+
\int d^3 \Brho
\lr{ \spacegrad_\Brho \inv{\Abs{\ncap + \Brho}}} \cdot \spacegrad_\Brho \inv{ \Abs{\Brho} }
+
\int d^3 \Brho
\inv{\Abs{\ncap + \Brho}} \spacegrad_\Brho^2 \inv{ \Abs{\Brho} } \\
&=
2 \int d^3 \Brho
\lr{ \spacegrad_\Brho \inv{\Abs{\ncap + \Brho}}} \cdot \spacegrad_\Brho \inv{ \Abs{\Brho} } \\
&- 4 \pi
\int d^3 \Brho’
\inv{\Abs{\Brho’ – \ncap}} \delta^3(\Brho’)
– 4 \pi
\int d^3 \Brho
\inv{\Abs{\Brho + \ncap}} \delta^3(\Brho) \\
&=
2 \int d^3 \Brho
\lr{ \spacegrad_\Brho \inv{\Abs{\ncap + \Brho}}} \cdot \spacegrad_\Brho \inv{ \Abs{\Brho} }
– 8 \pi.
\end{aligned}
\end{equation}

This used the Laplacian representation of the delta function \( \delta^3(\Bx) = -(1/4\pi) \spacegrad^2 (1/\Abs{\Bx}) \). Back-substitution gives

\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:380}
\int d^3 \Brho
\frac{ \Brho }{ \Abs{\Brho}^3 } \cdot
\frac{(\ncap + \Brho)}{ \Abs{\ncap + \Brho}^3}
=
4 \pi
+
\oint dA’ \ncap’ \cdot \spacegrad_\Brho \inv{ \Abs{\Brho} \Abs{\ncap + \Brho}}.
\end{equation}

We can argue that this last integral tends to zero, since

\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:400}
\begin{aligned}
\oint dA’ \ncap’ \cdot \spacegrad_\Brho \inv{ \Abs{\Brho} \Abs{\ncap + \Brho}}
&=
\oint dA’ \ncap’ \cdot \lr{
\lr{ \spacegrad_\Brho \inv{ \Abs{\Brho}} } \inv{\Abs{\ncap + \Brho}}
+
\inv{ \Abs{\Brho}} \lr{ \spacegrad_\Brho \inv{\Abs{\ncap + \Brho}} }
} \\
&=
-\oint dA’ \ncap’ \cdot \lr{
\frac{ \Brho } {\inv{ \Abs{\Brho}}^3 } \inv{\Abs{\ncap + \Brho}}
+
\inv{ \Abs{\Brho}} \frac{ (\Brho + \ncap) }{ \Abs{\ncap + \Brho}^3 }
} \\
&=
-\oint dA’ \inv{\Abs{\Brho} \Abs{\Brho + \ncap}}
\lr{
\frac{ \ncap’ \cdot \Brho }{
{\Abs{\Brho}}^2 }
+\frac{ \ncap’ \cdot (\Brho + \ncap) }{
{\Abs{\Brho + \ncap}}^2 }
}.
\end{aligned}
\end{equation}

The integrand in this surface integral is of \( O(1/\rho^3) \) so tends to zero on an infinite surface in the \( \Brho \) coordinate system. This completes the “easy” integral, leaving

\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:420}
\int d^3 \Brho
\frac{ \Brho }{ \Abs{\Brho}^3 } \cdot
\frac{(\ncap + \Brho)}{ \Abs{\ncap + \Brho}^3}
=
4 \pi.
\end{equation}

The total field energy can now be expressed as a sum of the self energies and the interaction energy
\begin{equation}\label{eqn:electrostaticJacksonSelfEnergy:440}
W =
\frac{1}{32 \pi^2 \epsilon} \int d^3 \Bx \frac{q_1^2}{\Abs{\Bx – \Bx_1}^4 }
+
\frac{1}{32 \pi^2 \epsilon} \int d^3 \Bx \frac{q_2^2}{\Abs{\Bx – \Bx_2}^4 }
+ \inv{ 4 \pi \epsilon}
\frac{q_1 q_2}{\Abs{\Bx_1 – \Bx_2} }.
\end{equation}

The interaction energy is exactly the potential energies for the two particles, the this total energy in the field is biased in the positive direction by the pair of self energies. It is interesting that the energy obtained from integrating the field energy density contains such self energy terms, but I don’t know exactly what to make of them at this point in time.

References

[1] JD Jackson. Classical Electrodynamics. John Wiley and Sons, 2nd edition, 1975.

Helmholtz theorem

October 1, 2016 math and physics play , , , , , , , , , , , , , ,

[Click here for a PDF of this post with nicer formatting]

This is a problem from ece1228. I attempted solutions in a number of ways. One using Geometric Algebra, one devoid of that algebra, and then this method, which combined aspects of both. Of the three methods I tried to obtain this result, this is the most compact and elegant. It does however, require a fair bit of Geometric Algebra knowledge, including the Fundamental Theorem of Geometric Calculus, as detailed in [1], [3] and [2].

Question: Helmholtz theorem

Prove the first Helmholtz’s theorem, i.e. if vector \(\BM\) is defined by its divergence

\begin{equation}\label{eqn:helmholtzDerviationMultivector:20}
\spacegrad \cdot \BM = s
\end{equation}

and its curl
\begin{equation}\label{eqn:helmholtzDerviationMultivector:40}
\spacegrad \cross \BM = \BC
\end{equation}

within a region and its normal component \( \BM_{\textrm{n}} \) over the boundary, then \( \BM \) is
uniquely specified.

Answer

The gradient of the vector \( \BM \) can be written as a single even grade multivector

\begin{equation}\label{eqn:helmholtzDerviationMultivector:60}
\spacegrad \BM
= \spacegrad \cdot \BM + I \spacegrad \cross \BM
= s + I \BC.
\end{equation}

We will use this to attempt to discover the relation between the vector \( \BM \) and its divergence and curl. We can express \( \BM \) at the point of interest as a convolution with the delta function at all other points in space

\begin{equation}\label{eqn:helmholtzDerviationMultivector:80}
\BM(\Bx) = \int_V dV’ \delta(\Bx – \Bx’) \BM(\Bx’).
\end{equation}

The Laplacian representation of the delta function in \R{3} is

\begin{equation}\label{eqn:helmholtzDerviationMultivector:100}
\delta(\Bx – \Bx’) = -\inv{4\pi} \spacegrad^2 \inv{\Abs{\Bx – \Bx’}},
\end{equation}

so \( \BM \) can be represented as the following convolution

\begin{equation}\label{eqn:helmholtzDerviationMultivector:120}
\BM(\Bx) = -\inv{4\pi} \int_V dV’ \spacegrad^2 \inv{\Abs{\Bx – \Bx’}} \BM(\Bx’).
\end{equation}

Using this relation and proceeding with a few applications of the chain rule, plus the fact that \( \spacegrad 1/\Abs{\Bx – \Bx’} = -\spacegrad’ 1/\Abs{\Bx – \Bx’} \), we find

\begin{equation}\label{eqn:helmholtzDerviationMultivector:720}
\begin{aligned}
-4 \pi \BM(\Bx)
&= \int_V dV’ \spacegrad^2 \inv{\Abs{\Bx – \Bx’}} \BM(\Bx’) \\
&= \gpgradeone{\int_V dV’ \spacegrad^2 \inv{\Abs{\Bx – \Bx’}} \BM(\Bx’)} \\
&= -\gpgradeone{\int_V dV’ \spacegrad \lr{ \spacegrad’ \inv{\Abs{\Bx – \Bx’}}} \BM(\Bx’)} \\
&= -\gpgradeone{\spacegrad \int_V dV’ \lr{
\spacegrad’ \frac{\BM(\Bx’)}{\Abs{\Bx – \Bx’}}
-\frac{\spacegrad’ \BM(\Bx’)}{\Abs{\Bx – \Bx’}}
} } \\
&=
-\gpgradeone{\spacegrad \int_{\partial V} dA’
\ncap \frac{\BM(\Bx’)}{\Abs{\Bx – \Bx’}}
}
+\gpgradeone{\spacegrad \int_V dV’
\frac{s(\Bx’) + I\BC(\Bx’)}{\Abs{\Bx – \Bx’}}
} \\
&=
-\gpgradeone{\spacegrad \int_{\partial V} dA’
\ncap \frac{\BM(\Bx’)}{\Abs{\Bx – \Bx’}}
}
+\spacegrad \int_V dV’
\frac{s(\Bx’)}{\Abs{\Bx – \Bx’}}
+\spacegrad \cdot \int_V dV’
\frac{I\BC(\Bx’)}{\Abs{\Bx – \Bx’}}.
\end{aligned}
\end{equation}

By inserting a no-op grade selection operation in the second step, the trivector terms that would show up in subsequent steps are automatically filtered out. This leaves us with a boundary term dependent on the surface and the normal and tangential components of \( \BM \). Added to that is a pair of volume integrals that provide the unique dependence of \( \BM \) on its divergence and curl. When the surface is taken to infinity, which requires \( \Abs{\BM}/\Abs{\Bx – \Bx’} \rightarrow 0 \), then the dependence of \( \BM \) on its divergence and curl is unique.

In order to express final result in traditional vector algebra form, a couple transformations are required. The first is that

\begin{equation}\label{eqn:helmholtzDerviationMultivector:800}
\gpgradeone{ \Ba I \Bb } = I^2 \Ba \cross \Bb = -\Ba \cross \Bb.
\end{equation}

For the grade selection in the boundary integral, note that

\begin{equation}\label{eqn:helmholtzDerviationMultivector:740}
\begin{aligned}
\gpgradeone{ \spacegrad \ncap \BX }
&=
\gpgradeone{ \spacegrad (\ncap \cdot \BX) }
+
\gpgradeone{ \spacegrad (\ncap \wedge \BX) } \\
&=
\spacegrad (\ncap \cdot \BX)
+
\gpgradeone{ \spacegrad I (\ncap \cross \BX) } \\
&=
\spacegrad (\ncap \cdot \BX)
–
\spacegrad \cross (\ncap \cross \BX).
\end{aligned}
\end{equation}

These give

\begin{equation}\label{eqn:helmholtzDerviationMultivector:721}
\boxed{
\begin{aligned}
\BM(\Bx)
&=
\spacegrad \inv{4\pi} \int_{\partial V} dA’ \ncap \cdot \frac{\BM(\Bx’)}{\Abs{\Bx – \Bx’}}
–
\spacegrad \cross \inv{4\pi} \int_{\partial V} dA’ \ncap \cross \frac{\BM(\Bx’)}{\Abs{\Bx – \Bx’}} \\
&-\spacegrad \inv{4\pi} \int_V dV’
\frac{s(\Bx’)}{\Abs{\Bx – \Bx’}}
+\spacegrad \cross \inv{4\pi} \int_V dV’
\frac{\BC(\Bx’)}{\Abs{\Bx – \Bx’}}.
\end{aligned}
}
\end{equation}

References

[1] C. Doran and A.N. Lasenby. Geometric algebra for physicists. Cambridge University Press New York, Cambridge, UK, 1st edition, 2003.

[2] A. Macdonald. Vector and Geometric Calculus. CreateSpace Independent Publishing Platform, 2012.

[3] Garret Sobczyk and Omar Le’on S’anchez. Fundamental theorem of calculus. Advances in Applied Clifford Algebras, 21:221–231, 2011. URL https://arxiv.org/abs/0809.4526.