5.6: Calculating Electric Fields of Charge Distributions
By the end of this section, you will be able to:
- Explain what a continuous source charge distribution is and how it is related to the concept of quantization of charge
- Describe line charges, surface charges, and volume charges
- Calculate the field of a continuous source charge distribution of either sign
The charge distributions we have seen so far have been discrete: made up of individual point particles. This is in contrast with a continuous charge distribution , which has at least one nonzero dimension. If a charge distribution is continuous rather than discrete, we can generalize the definition of the electric field. We simply divide the charge into infinitesimal pieces and treat each piece as a point charge.
Note that because charge is quantized, there is no such thing as a “truly” continuous charge distribution. However, in most practical cases, the total charge creating the field involves such a huge number of discrete charges that we can safely ignore the discrete nature of the charge and consider it to be continuous. This is exactly the kind of approximation we make when we deal with a bucket of water as a continuous fluid, rather than a collection of \(\textrm{H}_2\textrm{O}\) molecules.
Our first step is to define a charge density for a charge distribution along a line, across a surface, or within a volume, as shown in Figure \(\PageIndex{1}\).
Definitions of charge density:
- linear charge density: \(\lambda \equiv \) charge per unit length (Figure \(\PageIndex{1a}\)); units are coulombs per meter (\(\textrm{C}/\textrm{m}\))
- surface charge density: \(\sigma \equiv \) charge per unit area (Figure \(\PageIndex{1b}\)); units are coulombs per square meter \((\textrm{C}/\textrm{m}^2)\)
- volume charge density: \(\rho \equiv \) charge per unit volume (Figure \(\PageIndex{1c}\)); units are coulombs per square meter \((\textrm{C}/\textrm{m}^3)\)
For a line charge, a surface charge, and a volume charge, the summation in the definition of an Electric field discussed previously becomes an integral and \(q_i\) is replaced by \(dq = \lambda\,dl\), \(\sigma\,dA\), or \(\rho\,dV\), respectively:
\[ \begin{align} \vec{E}(P) &= \underbrace{\dfrac{1}{4\pi \epsilon_0} \sum_{i=1}^N \left(\dfrac{q_i}{r^2}\right)\hat{r}}_{\text{Point charges}} \label{eq1} \\[4pt] \vec{E}(P) &= \underbrace{\dfrac{1}{4\pi \epsilon_0} \int_{\textrm{line}} \left(\dfrac{\lambda \, dl}{r^2}\right) \hat{r}}_{\text{Line charge}} \label{eq2} \\[4pt] \vec{E}(P) &= \underbrace{\dfrac{1}{4\pi \epsilon_0} \int_{\textrm{surface}} \left(\dfrac{\sigma \,dA}{r^2}\right) \hat{r} }_{\text{Surface charge}}\label{eq3} \\[4pt] \vec{E}(P) &= \underbrace{\dfrac{1}{4\pi \epsilon_0} \int_{\textrm{volume}} \left(\dfrac{\rho \,dV}{r^2}\right) \hat{r}}_{\text{Volume charge}} \label{eq4} \end{align} \]
The integrals in Equations \ref{eq1}-\ref{eq4} are generalizations of the expression for the field of a point charge. They implicitly include and assume the principle of superposition. The “trick” to using them is almost always in coming up with correct expressions for \(dl\), \(dA\), or \(dV\), as the case may be, expressed in terms of r , and also expressing the charge density function appropriately. It may be constant; it might be dependent on location.
Note carefully the meaning of \(r\) in these equations: It is the distance from the charge element (\(q_i, \, \lambda \, dl, \, \sigma \, dA, \, \rho \, dV\)) to the location of interest, \(P(x, y, z)\) (the point in space where you want to determine the field). However, don’t confuse this with the meaning of \(\hat{r}\); we are using it and the vector notation \(\vec{E}\) to write three integrals at once. That is, Equation \ref{eq2} is actually
\[ \begin{align} E_x (P) &= \dfrac{1}{4\pi \epsilon_0} \int_{\textrm{line}} \left(\dfrac{\lambda \, dl}{r^2}\right)_x, \\[4pt] E_y(P) &= \dfrac{1}{4\pi \epsilon_0} \int_{\textrm{line}} \left(\dfrac{\lambda \, dl}{r^2}\right)_y, \\[4pt] E_z(P) &= \dfrac{1}{4\pi \epsilon_0} \int_{\textrm{line}} \left(\dfrac{\lambda \, dl}{r^2}\right)_z \end{align} \nonumber \]
Find the electric field a distance \(z\) above the midpoint of a straight line segment of length \(L\) that carries a uniform line charge density \(\lambda\).
Strategy
Since this is a continuous charge distribution, we conceptually break the wire segment into differential pieces of length \(dl\), each of which carries a differential amount of charge
\[dq = \lambda \, dl. \nonumber \]
Then, we calculate the differential field created by two symmetrically placed pieces of the wire, using the symmetry of the setup to simplify the calculation (Figure \(\PageIndex{2}\)). Finally, we integrate this differential field expression over the length of the wire (half of it, actually, as we explain below) to obtain the complete electric field expression.
Solution
Before we jump into it, what do we expect the field to “look like” from far away? Since it is a finite line segment, from far away, it should look like a point charge. We will check the expression we get to see if it meets this expectation.
The electric field for a line charge is given by the general expression
\[\vec{E}(P) = \dfrac{1}{4\pi \epsilon_0} \int_{\textrm{line}} \dfrac{\lambda dl}{r^2}\hat{r}. \nonumber \]
The symmetry of the situation (our choice of the two identical differential pieces of charge) implies the horizontal (\(x\))-components of the field cancel, so that the net field points in the \(z\)-direction. Let’s check this formally.
The total field \(\vec{E}(P)\) is the vector sum of the fields from each of the two charge elements (call them \(\vec{E}_1\) and \(\vec{E}_2\), for now):
\[ \begin{align*} \vec{E}(P) &= \vec{E}_1 + \vec{E}_2 \\[4pt] &= E_{1x}\hat{i} + E_{1z}\hat{k} + E_{2x} (-\hat{i}) + E_{2z}\hat{k}. \end{align*} \]
Because the two charge elements are identical and are the same distance away from the point \(P\) where we want to calculate the field, \(E_{1x} = E_{2x}\), so those components cancel. This leaves
\[ \begin{align*} \vec{E}(P) &= E_{1z}\hat{k} + E_{2z}\hat{k} \\[4pt] &= E_1 \, \cos\theta\,\hat{k} + E_2 \, \cos\theta\,\hat{k}. \end{align*} \]
These components are also equal, so we have
\[ \begin{align*} \vec{E}(P) &= \dfrac{1}{4 \pi \epsilon_0}\int \dfrac{\lambda dl}{r^2} \, \cos\theta\, \hat{k} + \dfrac{1}{4 \pi \epsilon_0}\int \dfrac{\lambda dl}{r^2} \, \cos\theta\, \hat{k} \\[4pt] &= \dfrac{1}{4 \pi \epsilon_0}\int_0^{L/2} \dfrac{2\lambda dx}{r^2} \, \cos\theta\, \hat{k} \end{align*} \]
where our differential line element \(dl\) is \(dx\), in this example, since we are integrating along a line of charge that lies on the \(x\)-axis. (The limits of integration are 0 to \(\frac{L}{2}\), not \(-\frac{L}{2}\) to \(+\frac{L}{2}\), because we have constructed the net field from two differential pieces of charge \(dq\). If we integrated along the entire length, we would pick up an erroneous factor of 2.)
In principle, this is complete. However, to actually calculate this integral, we need to eliminate all the variables that are not given. In this case, both \(r\) and \(\theta\) change as we integrate outward to the end of the line charge, so those are the variables to get rid of. We can do that the same way we did for the two point charges: by noticing that
\[r = (z^2 + x^2)^{1/2} \nonumber \]
and
\[\cos\theta = \dfrac{z}{r} = \dfrac{z}{(z^2 + x^2)^{1/2}}. \nonumber \]
Substituting, we obtain
\[ \begin{align*} \vec{E}(P) &= \dfrac{1}{4 \pi \epsilon_0}\int_0^{L/2} \dfrac{2\lambda dx}{(z^2 + x^2)} \dfrac{z}{(z^2 + x^2)^{1/2}} \hat{k} \\[4pt] &= \dfrac{1}{4 \pi \epsilon_0}\int_0^{L/2} \dfrac{2\lambda z}{(z^2 + x^2)^{3/2}} dx\,\hat{k} \\[4pt] &= \dfrac{2 \lambda z}{4 \pi \epsilon_0} \left[\dfrac{x}{z^2\sqrt{z^2 + x^2}}\right]_0^{L/2} \hat{k}. \end{align*} \]
which simplifies to
\[\vec{E}(z) = \dfrac{1}{4 \pi \epsilon_0} \dfrac{\lambda L}{z\sqrt{z^2 + \dfrac{L^2}{4}}} \, \hat{k}. \label{5.12} \]
Significance
Notice, once again, the use of symmetry to simplify the problem. This is a very common strategy for calculating electric fields. The fields of nonsymmetrical charge distributions have to be handled with multiple integrals and may need to be calculated numerically by a computer.
How would the strategy used above change to calculate the electric field at a point a distance \(z\) above one end of the finite line segment?
- Answer
-
We will no longer be able to take advantage of symmetry. Instead, we will need to calculate each of the two components of the electric field with their own integral.
Find the electric field a distance \(z\) above the midpoint of an infinite line of charge that carries a uniform line charge density \(\lambda\).
Strategy
This is exactly like the preceding example, except the limits of integration will be \(-\infty\) to \(+\infty\).
Solution
Again, the horizontal components cancel out, so we wind up with
\[\vec{E}(P) = \dfrac{1}{4 \pi \epsilon_0} \int_{-\infty}^{\infty} \dfrac{\lambda dx}{r^2} \, \cos \, \theta \hat{k} \nonumber \]
where our differential line element \(dl\) is \(dx\), in this example, since we are integrating along a line of charge that lies on the \(x\)-axis. Again,
\[ \begin{align*} \cos\theta &= \dfrac{z}{r} \\[4pt] &= \dfrac{z}{(z^2 + x^2)^{1/2}}. \end{align*} \nonumber \]
Substituting, we obtain
\[ \begin{align*} \vec{E}(P) &= \dfrac{1}{4 \pi \epsilon_0} \int_{-\infty}^{\infty} \dfrac{\lambda dx}{(z^2 + x^2)} \, \dfrac{z}{(z^2 + x^2)^{1/2}} \hat{k} \\[4pt] &= \dfrac{1}{4 \pi \epsilon_0} \int_{-\infty}^{\infty} \dfrac{\lambda z}{(z^2 + x^2)^{3/2}}dx\,\hat{k} \\[4pt] &= \dfrac{1}{4 \pi \epsilon_0} \left[ \dfrac{x}{z^2\sqrt{z^2 + x^2}}\right]_{-\infty}^{\infty} \, \hat{k} \end{align*} \]
which simplifies to
\[\vec{E}(z) = \dfrac{1}{4 \pi \epsilon_0} \dfrac{2\lambda}{z}\hat{k}. \nonumber \]
Significance
Our strategy for working with continuous charge distributions also gives useful results for charges with infinite dimension.
In the case of a finite line of charge, note that for \(z \gg L\), \(z^2\) dominates the \(L\) in the denominator, so that Equation \ref{5.12} simplifies to
\[\vec{E} \approx \dfrac{1}{4\pi \epsilon_0} \dfrac{\lambda L}{z^2} \hat{k}. \nonumber \]
If you recall that \(\lambda L = q\) the total charge on the wire, we have retrieved the expression for the field of a point charge, as expected.
In the limit \(L \rightarrow \infty\) on the other hand, we get the field of an infinite straight wire , which is a straight wire whose length is much, much greater than either of its other dimensions, and also much, much greater than the distance at which the field is to be calculated:
\[\vec{E}(z) = \dfrac{1}{4 \pi \epsilon_0} \dfrac{2\lambda}{z}\hat{k}. \label{infinite straight wire} \]
An interesting artifact of this infinite limit is that we have lost the usual \(1/r^2\) dependence that we are used to. This will become even more intriguing in the case of an infinite plane.
A ring has a uniform charge density \(\lambda\), with units of coulomb per unit meter of arc. Find the electric field at a point on the axis passing through the center of the ring.
Strategy
We use the same procedure as for the charged wire. The difference here is that the charge is distributed on a circle. We divide the circle into infinitesimal elements shaped as arcs on the circle and use polar coordinates shown in Figure \(\PageIndex{3}\).
Solution
The electric field for a line charge is given by the general expression
\[\vec{E}(P) = \dfrac{1}{4\pi \epsilon_0} \int_{\textrm{line}} \dfrac{\lambda dl}{r^2} \hat{r}. \nonumber \]
A general element of the arc between \(\theta\) and \(\theta + d\theta\) is of length \(R\,d\theta\) and therefore contains a charge equal to \(\lambda R \,d\theta\). The element is at a distance of \(r = \sqrt{z^2 + R^2}\) from \(P\), the angle is \(\cos\phi = \dfrac{z}{\sqrt{z^2+R^2}}\) and therefore the electric field is
\[ \begin{align*} \vec{E}(P) &= \dfrac{1}{4\pi \epsilon_0} \int_{\textrm{line}} \dfrac{\lambda dl}{r^2} \hat{r} = \dfrac{1}{4\pi \epsilon_0} \int_0^{2\pi} \dfrac{\lambda Rd\theta}{z^2 + R^2} \dfrac{z}{\sqrt{z^2 + R^2}} \hat{z} \\[4pt] &= \dfrac{1}{4\pi \epsilon_0} \dfrac{\lambda Rz}{(z^2 + R^2)^{3/2}} \hat{z} \int_0^{2\pi} d\theta \\[4pt] &= \dfrac{1}{4\pi \epsilon_0} \dfrac{2\pi \lambda Rz}{(z^2 + R^2)^{3/2}} \hat{z} \\[4pt] &= \dfrac{1}{4\pi \epsilon_0} \dfrac{q_{tot}z}{(z^2 + R^2)^{3/2}} \hat{z}. \end{align*} \]
Significance
As usual, symmetry simplified this problem, in this particular case resulting in a trivial integral. Also, when we take the limit of \(z \gg R\), we find that
\[\vec{E} \approx \dfrac{1}{4\pi \epsilon_0} \dfrac{q_{tot}}{z^2} \hat{z}, \nonumber \]
as we expect.
Find the electric field of a circular thin disk of radius \(R\) and uniform charge density at a distance \(z\) above the center of the disk (Figure \(\PageIndex{4}\))
Strategy
The electric field for a surface charge is given by
\[\vec{E}(P) = \dfrac{1}{4 \pi \epsilon_0} \int_{\textrm{surface}} \dfrac{\sigma\,dA}{r^2} \hat{r}. \nonumber \]
To solve surface charge problems, we break the surface into symmetrical differential “stripes” that match the shape of the surface; here, we’ll use rings, as shown in the figure. Again, by symmetry, the horizontal components cancel and the field is entirely in the vertical \((\hat{k})\) direction. The vertical component of the electric field is extracted by multiplying by \(\theta\), so
\[\vec{E}(P) = \dfrac{1}{4\pi \epsilon_0} \int_{surface} \dfrac{\sigma\,dA}{r^2} \, \cos\theta \, \hat{k}. \nonumber \]
As before, we need to rewrite the unknown factors in the integrand in terms of the given quantities. In this case,
\[dA = 2\pi r'dr' \nonumber \]
\[r^2 = r'^2 + z^2 \nonumber \]
\[\cos\theta = \dfrac{z}{(r'^2 + z^2)^{1/2}}. \nonumber \]
(Please take note of the two different “\(r\)’s” here; \(r\) is the distance from the differential ring of charge to the point \(P\) where we wish to determine the field, whereas \(r'\) is the distance from the center of the disk to the differential ring of charge.) Also, we already performed the polar angle integral in writing down \(dA\).
Solution
Substituting all this in, we get
\[ \begin{align*} \vec{E}(P) &= \vec{E}(z) \\[4pt] &= \dfrac{1}{4 \pi \epsilon_0} \int_0^R \dfrac{\sigma (2\pi r' dr')z}{(r'^2 + z^2)^{3/2}} \hat{k} \\[4pt] &= \dfrac{1}{4 \pi \epsilon_0} (2\pi \sigma z)\left(\dfrac{1}{z} - \dfrac{1}{\sqrt{R^2 + z^2}}\right) \hat{k} \end{align*} \]
or, more simply,
\[\vec{E}(z) = \dfrac{1}{4 \pi \epsilon_0} \left( 2 \pi \sigma - \dfrac{2 \pi \sigma z}{\sqrt{R^2 + z^2}}\right)\hat{k}. \label{5.14} \]
Significance
Again, it can be shown (via a Taylor expansion) that when \(z \gg R\), this reduces to
\[\vec{E}(z) \approx \dfrac{1}{4 \pi \epsilon_0} \dfrac{\sigma \pi R^2}{z^2} \hat{k},\nonumber \]
which is the expression for a point charge \(Q = \sigma \pi R^2\).
How would the above limit change with a uniformly charged rectangle instead of a disk?
- Answer
-
The point charge would be \(Q = \sigma ab\) where \(a\) and \(b\) are the sides of the rectangle but otherwise identical.
As \(R \rightarrow \infty\), Equation \ref{5.14} reduces to the field of an infinite plane , which is a flat sheet whose area is much, much greater than its thickness, and also much, much greater than the distance at which the field is to be calculated:
\[ \begin{align} \vec{E} &= \lim_{R \rightarrow \infty} \dfrac{1}{4 \pi \epsilon_0} \left( 2 \pi \sigma - \dfrac{2 \pi \sigma z}{\sqrt{R^2 + z^2}}\right)\hat{k} \\[4pt] &= \dfrac{\sigma}{2 \epsilon_0} \hat{k}. \label{5.15} \end{align} \]
Note that this field is constant. This surprising result is, again, an artifact of our limit, although one that we will make use of repeatedly in the future. To understand why this happens, imagine being placed above an infinite plane of constant charge. Does the plane look any different if you vary your altitude? No—you still see the plane going off to infinity, no matter how far you are from it. It is important to note that Equation \ref{5.15} is because we are above the plane. If we were below, the field would point in the \(- \hat{k}\) direction.
Find the electric field everywhere resulting from two infinite planes with equal but opposite charge densities (Figure \(\PageIndex{5}\)).
Strategy
We already know the electric field resulting from a single infinite plane, so we may use the principle of superposition to find the field from two.
Solution
The electric field points away from the positively charged plane and toward the negatively charged plane. Since the \(\sigma\) are equal and opposite, this means that in the region outside of the two planes, the electric fields cancel each other out to zero. However, in the region between the planes, the electric fields add, and we get
\[\vec{E} = \dfrac{\sigma}{\epsilon_0}\hat{i} \nonumber \]
for the electric field. The \(\hat{i}\) is because in the figure, the field is pointing in the \(+x\)-direction.
Significance
Systems that may be approximated as two infinite planes of this sort provide a useful means of creating uniform electric fields.
What would the electric field look like in a system with two parallel positively charged planes with equal charge densities?
- Answer
-
The electric field would be zero in between, and have magnitude \(\dfrac{\sigma}{\epsilon_0}\) everywhere else.