Skip to content

Volume Cosine Charge at the Center

Irodov 3.18DifficultyDifficulty 3 of 5UniquenessUniqueness 3 of 5

Problem Irodov 3.18

A solid ball of radius RR carries volume charge with density ρ=a⋅r\rho = \vect{a}\cdot\vect{r}, where a\vect{a} is a constant vector and r\vect{r} is the position vector measured from the ball’s center. Find the electric field vector at the center of the ball.

The ball of radius R in section along a, with ρ = a· r: positive on the side a points to, negative opposite, and denser farther from the plane through O normal to a.

The ball of radius RR in section along a\vect a, with ρ=a⋅r\rho = \vect a\cdot\vect r: positive on the side a\vect a points to, negative opposite, and denser farther from the plane through OO normal to a\vect a.

Given
  • ball of radius RR
  • ρ=a⋅r\rho = \vect a\cdot\vect r, a\vect a constant
Find
the field vector E\vect E at the center OO

Solution

Step 1The density in spherical coordinates

Take the xx axis along a\vect a and spherical coordinates rr, θ\theta, ϕ\phi about it, with θ\theta measured from a\vect a. The density is

ρ=arcos⁡θ(1) \rho = a r\cos\theta \tag{1}

and the volume element is dV=r2sin⁡θ dr dθ dϕdV = r^2\sin\theta\,dr\,d\theta\,d\phi.

Step 2The field of one volume element

The element dVdV at (r,θ,ϕ)(r, \theta, \phi) is the distance rr from OO, and the line from it to OO makes the angle θ\theta with the xx axis.

A volume element dV at distance r from the center and polar angle θ from a, its field d E at O, and the component dE_x along a.

A volume element dVdV at distance rr from the center and polar angle θ\theta from a\vect a, its field dEd\vect E at OO, and the component dExdE_x along a\vect a.

The components across the axis cancel between elements at ϕ\phi and ϕ+π\phi + \pi, so only the component along a\vect a survives. For ρ>0\rho > 0 the element pushes the field at OO away from itself, toward −x-x:

dEx=−14πε0 ρ dVr2cos⁡θdE_x = -\frac{1}{4\pi\varepsilon_0}\,\frac{\rho\,dV}{r^2}\cos\theta

Substituting (1) and dVdV, the factor r2r^2 cancels:

dEx=−a4πε0 rcos⁡2θsin⁡θ dr dθ dϕ(2) dE_x = -\frac{a}{4\pi\varepsilon_0}\,r\cos^2\theta\sin\theta\,dr\,d\theta\,d\phi \tag{2}
Step 3Integrate over the ball

The integrand of (2) separates into three one-variable integrals:

Ex=−a4πε0∫0Rr dr∫0πcos⁡2θsin⁡θ dθ∫02πdϕ=−a4πε0⋅R22⋅[−cos⁡3θ3]0π⋅2π=−a4πε0⋅R22⋅23⋅2π=−aR26ε0\begin{aligned} E_x &= -\frac{a}{4\pi\varepsilon_0}\int_0^{R}r\,dr\int_0^{\pi}\cos^2\theta\sin\theta\,d\theta\int_0^{2\pi}d\phi\\ &= -\frac{a}{4\pi\varepsilon_0}\cdot\frac{R^2}{2}\cdot\left[-\frac{\cos^3\theta}{3}\right]_0^{\pi}\cdot 2\pi\\ &= -\frac{a}{4\pi\varepsilon_0}\cdot\frac{R^2}{2}\cdot\frac{2}{3}\cdot 2\pi = -\frac{aR^2}{6\varepsilon_0} \end{aligned}

The field lies along a\vect a, so the vector is

E=−R2 a6ε0\vect E = -\frac{R^2\,\vect a}{6\varepsilon_0}
(3)

The field at the center has magnitude aR2/(6ε0)aR^2/(6\varepsilon_0) and points against a\vect a, from the positive half of the ball toward the negative half.

The total field E at the center of the ball, against a.

The total field E\vect E at the center of the ball, against a\vect a.

Check

By shells instead. The shell between rr and r+drr + dr is a sphere carrying the surface density ρ dr=(a dr) rcos⁡θ\rho\,dr = (a\,dr)\,r\cos\theta, a σ0cos⁡θ\sigma_0\cos\theta shell with σ0=ar dr\sigma_0 = a r\,dr; and such a shell is what two balls of densities ±ρ\pm\rho leave when shifted apart by l\vect l with ρl=σ0\rho l = \sigma_0. Inside, each ball gives ±ρr/(3ε0)\pm\rho\vect r/(3\varepsilon_0) from its own center by Gauss’s law, and the two differ by the uniform −ρl/(3ε0)-\rho\vect l/(3\varepsilon_0), so the shell’s field at its center is −σ0/(3ε0)-\sigma_0/(3\varepsilon_0) along a\vect a, whatever its radius. Summing the shells,

Ex=−∫0Rar dr3ε0=−aR26ε0,E_x = -\int_0^R \frac{a r\,dr}{3\varepsilon_0} = -\frac{aR^2}{6\varepsilon_0},

as found. Units: aa is a charge density per length, so aR2/ε0aR^2/\varepsilon_0 is (charge per area)/ε0/\varepsilon_0, a field.

Remark

The ball’s dipole moment is p=∫rρ dV=(4πR5/15) a\vect p = \int\vect r\rho\,dV = (4\pi R^5/15)\,\vect a, by the same three integrals with an extra factor rr. A uniformly polarized ball with that moment, whose charge sits on its surface, would give only −p/(4πε0R3)=−R2a/(15ε0)-\vect p/(4\pi\varepsilon_0R^3) = -R^2\vect a/(15\varepsilon_0) at its center; here the charge fills the interior and the inner layers, though weak, sit closer to OO, so the field at the center is two and a half times larger.