Rod of Charge
 

As an exercise in the "brute force" integration of COULOMB'S LAW (unavoidable in most cases), here is one way to find the electric field due to a uniformly charged, skinny rod of finite length $L$:

If the total charge $Q$ is uniformly distributed along the rod, then the charge per unit length is

\begin{displaymath}
\lambda = {Q \over L} \; .
\end{displaymath} (1)

\epsfig{figure=PS/Rod_of_Charge.eps,angle=0,width=0.95\textwidth}


We want to evaluate the electric field $\Vec{E}$ at an arbitrary "test point" in space. Such a point can be characterized completely by its perpendicular distance $r$ from the rod and its distance $h$ up from the bottom end of the rod, measured parallel to the rod, as shown. We choose to look at the rod and the point in their common plane.

Then we pick an arbitrarily position a distance $z$ up from the bottom end of the rod, as shown. A small slice of the rod (width $dz$) at that position contains a "charge element" $dq = \lambda \, dz$ which contributes $d\Vec{E}$ to the electric field vector $\Vec{E}$ at the test point. Coulomb's Law says that $d\Vec{E}$ points away from the charge element (assuming positive charge) and has a magnitude

\begin{displaymath}
dE = {k_E \lambda \, dz \over R^2}
\end{displaymath} (2)

where
\begin{displaymath}
R = \sqrt{ r^2 + (z-h)^2 }
\end{displaymath} (3)

is the distance from the charge element to the test point. In general $d\Vec{E}$ makes an angle $\theta$ with the direction $\Hat{r}$ perpendicular to the rod. If we define the $\Hat{z}$ direction to be "up" parallel to the rod, then we can separate $d\Vec{E}$ until its $r$ and $z$ components:
\begin{displaymath}
dE_r = dE \; \cos \theta
\end{displaymath} (4)

\begin{displaymath}
dE_z = - dE \; \sin \theta
\end{displaymath} (5)

To integrate these equations we need to convert all variables to match the differential (over which we integrate). We could use Eq. (3) to express $R$ in terms of $z$ (where $r$ and $h$ are constants) and use
\begin{displaymath}
\cos \theta = {r \over R}
\end{displaymath} (6)

\begin{displaymath}
\sin \theta = {(z-h) \over R}
\end{displaymath} (7)

but this would leave us with integrals that cannot be solved by inspection.

If we want to solve this problem without reference to external aids (like tables of integrals), it is better to convert into angles and trigonometric functions as follows:

Equation (6) can be rewritten

\begin{displaymath}
{1 \over R^2} = {\cos^2 \theta \over r^2}
\end{displaymath} (8)

and since
\begin{displaymath}
z - h = r \tan \theta \; ,
\end{displaymath} (9)

giving
\begin{displaymath}
dz = r \sec^2 \theta \; d\theta = {r \, d\theta \over \cos^2 \theta} \; ,
\end{displaymath} (10)

we can write Eq. (2) as
\begin{displaymath}
dE = k_E \lambda \, \left(\cos^2 \theta \over r^2 \right)
 . . . 
 . . . ver \cos^2 \theta \right)
= {k_E \lambda \over r} \, d\theta
\end{displaymath} (11)

and from that, Eqs. (4) and (5), respectively, as
\begin{displaymath}
dE_r = {k_E \lambda \over r} \, \cos \theta \; d\theta
= {k_E \lambda \over r} \, du
\end{displaymath} (12)

where $u \equiv \sin \theta$, and
\begin{displaymath}
dE_z = - {k_E \lambda \over r} \, \sin \theta \; d\theta
= {k_E \lambda \over r} \, dv
\end{displaymath} (13)

where $v \equiv \cos \theta$.

Integrating these differentials is trivial; we are left with just the differences between $u$ (or $v$) at the limits of integration (the top and bottom of the rod):

\begin{displaymath}
E_r = {k_E \lambda \over r} \, \left[
{ (L-h) \over \sqrt{r^2 + (L-h)^2} }
+ { h \over \sqrt{r^2 + h^2} }
\right]
\end{displaymath} (14)

(note that $u$ is negative at the bottom) and
\begin{displaymath}
E_z = k_E \lambda \, \left[
{ 1 \over \sqrt{r^2 + (L-h)^2} }
- { 1 \over \sqrt{r^2 + h^2} }
\right]
\end{displaymath} (15)

These equations express a completely general solution to this problem.

Let's check to see what these give for the field directly out from the midpoint of the rod - i.e. for $h = L/2$:

$\displaystyle E_r$ $\textstyle =$ $\displaystyle {k_E \lambda \over r} \, \left[
{ L/2 \over \sqrt{r^2 + L^2/4} }
+ { L/2 \over \sqrt{r^2 + L^2/4} }
\right]
\cr$ (16)

and
$\displaystyle E_z$ $\textstyle =$ $\displaystyle k_E \lambda \, \left[
{ 1 \over \sqrt{r^2 + L^2/4} }
- { 1 \over \sqrt{r^2 + L^2/4} }
\right]
\cr$ (17)


Let's also check to see what we get for $E_r$ (at the midpoint) very far from the rod ($r \gg L$):

\begin{displaymath}
E_r \goestoas{r \to \infty} {k_E \lambda L \over r^2}
= {k_E Q \over r^2}
\end{displaymath} (18)

(i.e. Coulomb's Law) $^\surd$

and very close to the rod ($r \ll L$):

\begin{displaymath}
E_r \goestoas{r \to 0} {2 k_E \lambda \over r} \; .
\end{displaymath} (19)

The last result can be used as the field due to an infinitely long uniform line of charge. But there is a much easier way to obtain it . . . .


Jess H. Brewer - Last modified: Mon Nov 16 17:22:34 PST 2015