Lessons · circuits

The field of a charged disk: between a point and a sheet

Regime 2 — calculus does what algebra cannot. The integrand isn't constant, so no single algebra product gives the answer; the accumulated quantity is the area under the curve — and SymPy proves that area is exactly the closed-form result.

Regime 2 · calculus does more Math machine-derived & checked by SymPy 1 modeling assumptions (author-asserted)

A disk of radius R=0.3R = 0.3 m carries a uniform surface charge density σ=100\sigma = 100 nC/m². We want the electric field a distance z=0.1z = 0.1 m out along its axis. Algebra hands you two extremes — a point charge kQ/z2kQ/z^2 and an infinite sheet σ/2ε0\sigma/2\varepsilon_0 — but a real disk is neither, and there is no algebra formula for one. Calculus is the only way in: slice the disk into concentric rings dq=σ2πrdrdq = \sigma\,2\pi r\,dr, each contributing kdqz/(z2+r2)3/2k\,dq\,z/(z^2+r^2)^{3/2} on the axis, and integrate. The shaded area under those ring contributions is the field E=σ2ε0(1z/z2+R2)E = \tfrac{\sigma}{2\varepsilon_0}(1 - z/\sqrt{z^2+R^2}) — which collapses to the point charge far away and to the infinite sheet for a huge disk.

Total field E=σ2ε0 ⁣(1zz2+R2)E = \dfrac{\sigma}{2\varepsilon_0}\!\left(1 - \dfrac{z}{\sqrt{z^2+R^2}}\right)
3862 N/C
Point-charge approx kQ/z2kQ/z^2 (exact only far away; overestimates up close)
25420 N/C
Infinite-sheet limit σ/2ε0\sigma/2\varepsilon_0 (overestimate)
5649 N/C
Slice the disk into rings and add their on-axis fields
dE=kdqz(z2+r2)3/2,dq=σ(2πrdr)dE = \frac{k\,dq\,z}{(z^2+r^2)^{3/2}}, \qquad dq = \sigma\,(2\pi r\,dr)

Cut the disk into concentric rings of radius rr and width drdr, each carrying charge dq=σ2πrdrdq = \sigma\,2\pi r\,dr. By symmetry only the axial part of each ring's field survives, kdqz/(z2+r2)3/2k\,dq\,z/(z^2+r^2)^{3/2}. Summing over the rings is an integral — the only way to handle a continuous distribution.

check dE/dr=2πkσzr/(z2+r2)3/2dE/dr = 2\pi k\sigma z\,r/(z^2+r^2)^{3/2}; check E=0rfdrE = \int_0^r f\,dr'; far limit kQ/z2kQ/z^2 (zz\to\infty); large-disk limit σ/2ε0\sigma/2\varepsilon_0 (RR\to\infty)

  • The accumulated field's slope is the ring contribution: dEdr=2πkσzr(z2+r2)3/2\dfrac{dE}{dr} = \dfrac{2\pi k\sigma z\,r}{(z^2+r^2)^{3/2}}. [structural]
  • The field is the area under the ring contributions: E(r)=0r2πkσzr(z2+r2)3/2drE(r) = \int_0^r \dfrac{2\pi k\sigma z\,r'}{(z^2+r'^2)^{3/2}}\,dr'. [structural]
  • Far away the disk looks like a point charge: EkQ/z2E \to kQ/z^2 with Q=σπR2Q = \sigma\pi R^2 as zz \to \infty. [structural]
  • A very large disk is an infinite sheet: Eσ/2ε0=2πkσE \to \sigma/2\varepsilon_0 = 2\pi k\sigma, independent of zz, as RR \to \infty. [structural]

Dimensional homogeneity: checked by SymPy (holds).

Common misconception: “A charged disk is basically a flat sheet, so its field is just σ/2ε0\sigma/2\varepsilon_0 — or, if you prefer, lump all its charge at the centre and use kQ/z2kQ/z^2.

Both are limits, not the answer. σ/2ε0\sigma/2\varepsilon_0 is the field of an infinite sheet (or the disk in the limit z0z \to 0); for a finite disk at finite distance it is an overestimate, because the missing charge beyond radius RR would only have added more field. kQ/z2kQ/z^2 is the far-field limit (zRz \gg R); up close it overestimates badly, because the inverse-square point model piles all the charge at one distance instead of spreading it from r=0r=0 out to RR. The true field is the integral E=σ2ε0(1z/z2+R2)E = \tfrac{\sigma}{2\varepsilon_0}(1 - z/\sqrt{z^2+R^2}), which lies below both extremes and only matches one of them in its own regime. Drag the cursor to sweep how much of the disk is included; the area under the ring-contribution curve is the field.

Modeling assumptions — author-asserted, disclosed not discharged
  • A flat disk with uniform surface charge density σ\sigma, and the field evaluated on its central axis — so by symmetry the off-axis components of each ring cancel and only the axial field survives.

The dE/dr–r graph, fully annotated

A static rendering (Matplotlib): the shaded area under dE/dr is the accumulated integral E, and the slope of E is dE/dr. The interactive version with a draggable cursor is in the Graph tab above.

2026-06-27T13:58:23.701412 image/svg+xml Matplotlib v3.11.0, https://matplotlib.org/