Electric Field Intensity on the Axis of a Uniformly Charged Disk

In the study of electrostatics, moving from discrete point charges to continuous charge distributions is a pivotal step. One of the most fundamental problems involves determining the electric field produced by a uniformly charged disk. This problem serves as an excellent bridge between basic Coulombic principles and the application of integral calculus in physics.

In this tutorial, we will systematically derive the expression for the electric field intensity along the central axis of a disk and analyze its behavior under different physical limits.

Problem Definition and Parameters

Consider a thin, non-conducting disk with the following characteristics:

Radius of the disk: $R$
Surface charge density: $\sigma$ (assumed to be uniform, where $\sigma = Q/A$)
Observation point: A point $P$ located on the central axis of the disk at a distance $z$ from the center.

We define the $z$-axis as the axis of symmetry, where $z > 0$ represents points above the disk and $z < 0$ represents points below it.

1. Leveraging Symmetry

Before diving into the calculus, we must apply symmetry arguments to simplify the problem. Because the disk is rotationally symmetric about the $z$-axis, for every charge element $dq$ at a position $(r, \phi)$, there is a corresponding element on the opposite side of the axis.

The radial components (perpendicular to the $z$-axis) of the electric field produced by these elements will point in opposite directions and cancel each other out. Consequently, the net electric field $\mathbf{E}$ must point strictly along the $z$-axis. This allows us to focus our derivation solely on the $z$-component, $E_z$.

2. Mathematical Derivation

To find the total field, we divide the disk into an infinite number of infinitesimal concentric rings.

2.1 Defining the Differential Element

Consider a thin ring with radius $r$ and an infinitesimal width $dr$.

The area of this ring is $dA = 2\pi r , dr$.
The charge contained within this ring is $dq = \sigma , dA = \sigma (2\pi r , dr)$.

2.2 The Contribution of a Single Ring

Let $R'$ be the distance from any point on the ring to the observation point $P$ on the axis. Using the Pythagorean theorem:
[ R' = \sqrt{r^2 + z^2} ]

According to Coulomb's Law, the magnitude of the differential electric field $dE$ produced by the charge $dq$ is:
[ dE = \frac{1}{4\pi\varepsilon_{0}} \frac{dq}{(R')^2} ]

As established by our symmetry argument, we only need the $z$-component. We can find this by multiplying $dE$ by $\cos\theta$, where $\theta$ is the angle between the field vector and the $z$-axis. From the geometry of the setup, $\cos\theta = \frac{z}{R'}$.

Thus, the differential $z$-component is:
[ dE_z = dE \cos\theta = \frac{1}{4\pi\varepsilon_{0}} \frac{dq}{(r^2 + z^2)} \cdot \frac{z}{\sqrt{r^2 + z^2}} ]
[ dE_z = \frac{1}{4\pi\varepsilon_{0}} \frac{\sigma (2\pi r , dr) z}{(r^2 + z^2)^{3/2}} ]

Simplifying the constants, we get:
[ dE_z = \frac{\sigma z}{2\varepsilon_{0}} \frac{r , dr}{(r^2 + z^2)^{3/2}} ]

2.3 Integration

To find the total electric field $E_z$, we integrate this expression from the center of the disk ($r = 0$) to its outer edge ($r = R$):
[ E_z(z) = \frac{\sigma z}{2\varepsilon_{0}} \int_{0}^{R} \frac{r}{(r^2 + z^2)^{3/2}} , dr ]

To solve this integral, we use the substitution $u = r^2 + z^2$, which implies $du = 2r , dr$ (or $r , dr = \frac{1}{2} du$). Changing the limits of integration:

When $r = 0$, $u = z^2$.
When $r = R$, $u = R^2 + z^2$.

The integral becomes:
[ E_z(z) = \frac{\sigma z}{2\varepsilon_{0}} \cdot \frac{1}{2} \int_{z^2}^{R^2 + z^2} u^{-3/2} , du ]
[ E_z(z) = \frac{\sigma z}{4\varepsilon_{0}} \left[ \frac{u^{-1/2}}{-1/2} \right]{z^2}^{R^2 + z^2} = \frac{\sigma z}{4\varepsilon{0}} \left[ -2u^{-1/2} \right]{z^2}^{R^2 + z^2} ]
[ E_z(z) = \frac{\sigma z}{2\varepsilon{0}} \left( \frac{1}{z} - \frac{1}{\sqrt{R^2 + z^2}} \right) ]

Distributing the $z$ term, we arrive at the final analytical expression:
[ \boxed{E_z(z) = \frac{\sigma}{2\varepsilon_{0}} \left( 1 - \frac{z}{\sqrt{R^2 + z^2}} \right)} ]

(Note: For $z > 0$, the field points away from the disk if $\sigma > 0$.)

3. Physical Interpretation and Limiting Cases

A powerful way to verify a derivation is to check if it behaves logically in extreme scenarios.

Scenario	Mathematical Limit	Simplified Result	Physical Meaning
Far-field	$z \gg R$	$E_z \approx \frac{\sigma R^2}{4\varepsilon_0 z^2} = \frac{Q}{4\pi\varepsilon_0 z^2}$	At great distances, the disk looks like a single point charge $Q$.
Near-field	$z \to 0$	$E_z \to \frac{\sigma}{2\varepsilon_0}$	Very close to the surface, the disk behaves like an infinite plane.
Infinite Disk	$R \to \infty$	$E_z = \frac{\sigma}{2\varepsilon_0}$	For an infinite sheet, the field is uniform and independent of distance.

4. Numerical Example

Problem: A disk has a radius $R = 5\text{ cm}$ and a uniform surface charge density $\sigma = 2 \times 10^{-6}\text{ C/m}^2$. Calculate the electric field strength at a point on the axis $3\text{ cm}$ above the center.

Given:

$R = 0.05\text{ m}$
$z = 0.03\text{ m}$
$\sigma = 2 \times 10^{-6}\text{ C/m}^2$
$\varepsilon_0 \approx 8.854 \times 10^{-12}\text{ F/m}$

Step 1: Calculate the geometric term
[ \sqrt{R^2 + z^2} = \sqrt{0.05^2 + 0.03^2} = \sqrt{0.0025 + 0.0009} = \sqrt{0.0034} \approx 0.0583\text{ m} ]

Step 2: Apply the formula
[ E_z = \frac{2 \times 10^{-6}}{2(8.854 \times 10^{-12})} \left( 1 - \frac{0.03}{0.0583} \right) ]
[ E_z \approx (1.129 \times 10^5) \times (1 - 0.5146) ]
[ E_z \approx 1.129 \times 10^5 \times 0.4854 \approx 5.48 \times 10^4\text{ V/m} ]

The electric field strength is approximately $5.48 \times 10^4\text{ V/m}$, directed away from the disk.

5. Common Pitfalls to Avoid

When performing these types of calculations, keep the following in mind:

Neglecting Symmetry: Attempting to integrate the full 3D vector field without accounting for the cancellation of radial components is a common source of error and unnecessary complexity.
Incorrect Integration Limits: Ensure your integration variable (in this case, $r$) starts from the center ($0$) and ends at the physical boundary ($R$).
Unit Inconsistency: Always convert measurements to SI units (meters, Coulombs, Farads) before plugging them into the formula. Mixing centimeters and meters will lead to errors of several orders of magnitude.

6. Summary

Through this derivation, we have moved from the fundamental Coulomb's Law to a complete analytical model for the electric field of a charged disk. We have demonstrated how:

Symmetry simplifies vector problems into scalar ones.
Calculus allows us to sum the contributions of continuous charge distributions.
Limiting cases provide a sanity check for our mathematical models.

Understanding this model provides a foundation for more advanced topics, such as non-uniform charge distributions, the electric field of charged cylinders, and the relationship between electric fields and electrostatic potential.