Field Strength Distribution of a Non-uniformly Charged Sphere

In introductory electromagnetism, students are often introduced to the idealized model of a uniformly charged sphere, where the charge density $\rho$ is constant throughout the volume. While this model is mathematically convenient, it rarely reflects the complexities of real-world physical systems. In advanced applications—such as modeling the charge distribution within an atomic nucleus, analyzing plasma physics, or understanding the polarization of dielectrics—the charge density is rarely constant. Instead, it varies as a function of the radial distance from the center.

A non-uniformly charged sphere is defined by a charge density $\rho(r)$ that depends on the radius $r$. Depending on the physical context, this distribution typically follows one of several mathematical patterns:

  • Linear Distribution: $\rho(r) = kr$, where the density increases or decreases linearly with the radius.
  • Power-Law Distribution: $\rho(r) = \rho_0 (r/R)^n$, a versatile model where the exponent $n$ dictates how sharply the charge concentrates toward the surface or the core.
  • Exponential Distribution: $\rho(r) = \rho_0 e^{-\alpha r}$, often used to describe systems where the charge density decays rapidly away from a central point.

To determine the electric field strength produced by such a distribution, we must move beyond simple algebraic formulas and employ integral calculus combined with fundamental electrostatic laws.

Theoretical Framework: Gauss's Law and Symmetry

The most efficient way to solve for the electric field of a non-uniform sphere is to exploit its radial symmetry. Because the charge distribution depends only on the distance $r$ from the origin, the resulting electric field $\mathbf{E}$ must also be purely radial. Its magnitude $E(r)$ is constant at any given radius, meaning the field lines point directly outward (or inward) and do not depend on the angular coordinates ($\theta, \phi$).

This symmetry allows us to apply Gauss's Law effectively. We choose a spherical Gaussian surface of radius $r$ centered at the origin. Gauss's Law states:

$$\oint_{S} \mathbf{E} \cdot d\mathbf{A} = \frac{Q_{\text{encl}}}{\epsilon_0}$$

For a spherically symmetric system, the surface integral simplifies significantly. Since $\mathbf{E}$ is parallel to the area vector $d\mathbf{A}$ and its magnitude is constant over the sphere, the integral becomes:

$$E(r) \cdot 4\pi r^2 = \frac{Q_{\text{encl}}(r)}{\epsilon_0}$$

The central challenge in these problems is not the application of Gauss's Law itself, but rather calculating the enclosed charge $Q_{\text{encl}}(r)$, which requires integrating the varying charge density over the volume.

Methodology for Calculating Enclosed Charge

To find the total charge contained within a Gaussian sphere of radius $r$, we perform a volume integral of the charge density function $\rho(r')$. In spherical coordinates, the differential volume element is $dV = 4\pi r'^2 dr'$. The general formula is:

$$Q_{\text{encl}}(r) = \int_{0}^{r} \rho(r') \cdot 4\pi r'^2 dr'$$

When solving these problems, it is critical to distinguish between two distinct spatial regions:

  1. Interior Region ($r < R$): The integration limit is the variable radius $r$. We are calculating the charge contained within the specific shell we are investigating.
  2. Exterior Region ($r \ge R$): The integration limit is the physical radius of the sphere $R$. In this case, the Gaussian surface encloses the total charge ($Q_{\text{total}}$) of the sphere, regardless of how large $r$ is.

Case Study: Linear Charge Density Distribution

To illustrate the mathematical rigor required, let us analyze a sphere of radius $R$ where the charge density increases linearly from the center to the surface.

Problem Definition:
Assume a sphere of radius $R$ has a charge density given by:
$$\rho(r) = \rho_0 \frac{r}{R} \quad \text{for } r \le R$$
And $\rho = 0$ for $r > R$. We seek to find the electric field strength $E(r)$ both inside and outside the sphere.

1. The Internal Electric Field ($r < R$)

First, we calculate the charge enclosed by a Gaussian surface at a radius $r$ within the sphere:

$$Q_{\text{encl}}(r) = \int_{0}^{r} \left( \rho_0 \frac{r'}{R} \right) 4\pi r'^2 dr'$$
$$Q_{\text{encl}}(r) = \frac{4\pi \rho_0}{R} \int_{0}^{r} r'^3 dr'$$
$$Q_{\text{encl}}(r) = \frac{4\pi \rho_0}{R} \left[ \frac{r'^4}{4} \right]_{0}^{r} = \frac{\pi \rho_0 r^4}{R}$$

Substituting this into the Gaussian equation:

$$E(r) \cdot 4\pi r^2 = \frac{\pi \rho_0 r^4}{R \epsilon_0}$$

Solving for $E(r)$ gives the internal field:
$$E_{\text{in}}(r) = \frac{\rho_0 r^2}{4 \epsilon_0 R}$$

2. The External Electric Field ($r \ge R$)

To find the field outside, we first determine the total charge $Q_{\text{total}}$ by integrating up to the boundary $R$:

$$Q_{\text{total}} = \frac{\pi \rho_0 R^4}{R} = \pi \rho_0 R^3$$

Applying Gauss's Law for a radius $r \ge R$:

$$E(r) \cdot 4\pi r^2 = \frac{\pi \rho_0 R^3}{\epsilon_0}$$

Solving for $E(r)$ gives the external field:
$$E_{\text{out}}(r) = \frac{\rho_0 R^3}{4 \epsilon_0 r^2}$$

Physical Interpretation and Analysis

The results of our derivation provide several profound insights into the nature of electrostatic fields:

  • Boundary Continuity: If we evaluate both expressions at the surface ($r = R$), we find:
    • $E_{\text{in}}(R) = \frac{\rho_0 R}{4 \epsilon_0}$
    • $E_{\text{out}}(R) = \frac{\rho_0 R}{4 \epsilon_0}$
      The field is continuous at the boundary. This is expected because there is no discrete surface charge density $\sigma$ at $r=R$; the charge is distributed through the volume.
  • Scaling Laws: In a uniform sphere, the internal field grows linearly with $r$ ($E \propto r$). However, in our linear non-uniform model, the field grows quadratically ($E \propto r^2$). This demonstrates that as charge becomes more concentrated toward the outer edges of the sphere, the electric field strength increases more aggressively as one moves outward from the center.
  • The Point Charge Equivalence: Notably, the external field $E_{\text{out}}$ follows the inverse-square law ($1/r^2$), identical to the field of a point charge. This confirms a fundamental principle of electrostatics: for any spherically symmetric distribution, the field outside the object behaves as if all the charge were concentrated at a single point at the center.

Summary of the Analytical Workflow

To master the study of non-uniform charge distributions, one should follow this systematic approach:

  1. Define the Density Function: Identify the mathematical form of $\rho(r)$.
  2. Integrate for Enclosed Charge: Use $Q_{\text{encl}}(r) = \int \rho(r') dV$ to find the charge within a radius $r$.
  3. Apply Gauss's Law: Use the symmetry-simplified form $E(r) = \frac{Q_{\text{encl}}(r)}{4\pi \epsilon_0 r^2}$.
  4. Evaluate Segmented Regions: Always perform separate calculations for the interior ($r < R$) and exterior ($r > R$) to account for the change in enclosed charge.