Gauss's Law (Integral Form) from Coulomb
Statement
For an electrostatic field built by superposing Coulomb fields, the outward flux of E through any closed surface S equals the net charge it encloses divided by ε0: ∮S E·dA = Qenc/ε0. The surface may have any shape and size; only the charge inside it survives, and each interior charge contributes the same amount independent of where inside it sits.
Why it matters
Coulomb's law is a statement about the field of one point charge; Gauss's law is a statement about the field of everything at once, encoded as a flux balance. Recasting Coulomb in this form is what makes electrostatics tractable: whenever a charge distribution has enough symmetry, a wisely chosen Gaussian surface reduces a hard vector integral to a one-line algebraic solve for E.
The derivation also exposes why the law is true. Nothing in it depends on the details of the surface — the whole result is geometry, driven by one fact: the Coulomb field falls off exactly as 1/r2, so its flux through an area element is precisely the solid angle that element subtends. Gauss's law is inverse-square made global.
Assumptions
Derivation
Result
Reading. The net number of field lines threading out of any closed surface is fixed by the charge trapped inside it and by nothing else — not the shape of the surface, not where inside the charges sit, not any charge outside. Positive enclosed charge is a net source (outward flux); negative charge is a net sink.
Units check. [E·dA] = (N·C−1)(m2) = N·m2·C−1. On the right, [ε0] = C2·N−1·m−2, so [q/ε0] = C / (C2·N−1·m−2) = N·m2·C−1. Both sides carry units of N·m2·C−1 (equivalently V·m). ✓
Limiting cases
- Concentric sphere, single charge. With S a sphere of radius R about q, symmetry gives Φ = E·4πR2; setting this to q/ε0 recovers E = q/4πε0R2 — the derivation runs backwards to Coulomb.
- Charge fully outside. Qenc = 0, so Φ = 0 for any surface, however strong the external field piercing it.
- Charge on the axis of symmetry but off-centre. Flux is unchanged (still q/ε0); only the distribution of E over S shifts.
- Neutral enclosure. Equal positive and negative charge inside ⇒ Qenc = 0 and zero net flux, even though E is nonzero everywhere on S.
- Small surface shrinking onto a point of density ρ. Φ/V → ρ/ε0, the seed of the differential form ∇·E = ρ/ε0.
Breaks when
- A charge lies on S. The subtended solid angle is neither 4π nor 0 but the fractional angle the surface presents there (2π for a locally flat patch → flux q/2ε0, and still less at an edge or vertex). Qenc becomes ambiguous and the equality fails as stated.
- The force is not exactly inverse-square. A massive-photon (Proca/Yukawa) field E ∝ (1+μr)e−μr r̂/r2 no longer makes (r̂·n̂)dA/r2 a pure solid angle; step 4 collapses and the flux acquires a volume term depending on the field throughout the enclosed region, not just the charge.
- Time-varying / radiating sources. Once E is not a static Coulomb superposition, this Coulomb-based route is invalid (retardation matters). Gauss's law survives — but as an independent Maxwell equation, not as a consequence of this argument.
- Unbounded charge or a singular field on S. If total enclosed charge diverges, or a point charge sits arbitrarily close to S, the flux integral fails to converge absolutely and the tidy result breaks.
Failure modes
- Position-dependence error. Believing the flux depends on where inside the surface the charge sits. It does not — only on whether it is enclosed.
- Counting external charges. Adding nearby outside charges into Qenc. They contribute exactly zero net flux (step 6, outside case), no matter how close or how strong.
- Dropping the sign of cos θ. Writing |cos θ| or |dΩ| in step 6, which destroys the entry/exit cancellation and wrongly gives nonzero flux for an exterior charge.
- Confusing E with Φ. Concluding that if Qenc = 0 then E = 0 on S. The field on the surface depends on all charges; only its net flux is set by the enclosed ones.
- Open surfaces. Applying Φ = q/ε0 to a surface that is not closed (a hemisphere alone, a disk). The 4π result needs a closed surface; an open one subtends a partial solid angle.
- Using symmetry that isn't there. Pulling E out of ∮E dA = E·A when the field is not constant and normal over S. The flux law always holds; the algebraic shortcut needs a symmetry-matched surface.
Discussion
The physical content lives entirely in step 4. The Coulomb field spreads its 4π steradians of flux outward and thins as 1/r2, exactly matching the way a fixed solid angle’s cross-section grows as r2. The two powers cancel, so every closed surface around a charge catches the same total flux — the field lines are conserved between charges, terminating only on charge. Gauss's law is the accountant's statement of that conservation.
Logically, Gauss's law is stronger than Coulomb's law even though we derived it from Coulomb. Coulomb's law is the electrostatic special case; Gauss's law is one of the four Maxwell equations and holds even for wildly time-dependent fields, where no simple Coulomb expression exists. What this derivation establishes is the equivalence in electrostatics: given superposition, "inverse-square point field" and "flux = enclosed charge / ε0" say the same thing.
Applying the divergence theorem, ∮S E·dA = ∫V (∇·E) dV, and matching integrands with Qenc = ∫V ρ dV gives the local form ∇·E = ρ/ε0 — the integral law holding for every surface is exactly the content of the differential law at every point.
The exact 1/r2 is not a convenience; it is a probe of deep physics. In distributional language ∇·(r̂/r2) = 4π δ3(r), so the point charge is a clean delta-function source and the 4π is the surface area of the unit sphere in three dimensions. Any deviation of the exponent from 2 would signal a nonzero photon mass; null-flux experiments (Cavendish's concentric spheres, refined by Williams, Faller and Hill) look for a would-be charge on the inner sphere and, finding none, bound the exponent to 2 to within parts in 1016 and the photon mass below ~10−14 eV/c2.
Common misconceptions. Gauss's law does not say the field is easy to compute — it says the flux is. The field E at a point on S generally depends on all charges everywhere; it is only the surface-integrated normal component that forgets the outside world. Solving for E pointwise from Gauss's law is legitimate only when symmetry guarantees E is uniform and normal (or tangential, contributing nothing) across a matched surface.
Worked examples
Example 1 — Flux is independent of the surface radius. A point charge q = 5.0 nC sits at the centre of a sphere of radius R = 0.10 m. Find the outward flux, and confirm it does not change if R is doubled.
Reading. Doubling R to 0.20 m quarters E but quadruples the area, leaving Φ = 565 N·m2·C−1 unchanged — the flux belongs to the charge, not the sphere.
Example 2 — Symmetry and solid angle: flux through one face of a cube. A charge q = 12 nC sits at the centre of a cube. Find the flux through a single face.
Reading. No integration over the awkward square face was needed: the solid-angle bookkeeping of the derivation does the work, giving q/6ε0 directly. Total over six faces returns 12×10−9/8.854×10−12 = 1355 N·m2·C−1 = q/ε0.
Problems
- A point charge q = 2.0 nC sits at the centre of a sphere of radius r. What is the flux, and how does it change when r is doubled?
Solution
Flux depends only on enclosed charge: Φ = q/ε0 = (2.0×10−9)/(8.854×10−12) = 226 N·m2·C−1. Doubling r leaves it unchanged: 226 N·m2·C−1.
- A closed surface encloses two charges, +3.0 nC and −8.0 nC; a third charge +5.0 nC lies just outside it. Find the net flux.
Solution
Only enclosed charge counts: Qenc = +3.0 − 8.0 = −5.0 nC; the external +5.0 nC contributes zero net flux. Φ = Qenc/ε0 = (−5.0×10−9)/(8.854×10−12) = −565 N·m2·C−1. The negative sign means net inward flux.
- A charge q = 6.0 nC sits at the centre of the flat circular base of a hemispherical surface. Find the flux through the curved (hemispherical) part alone.
Solution
From a point on the base plane, the curved surface caps the entire upper half-space and subtends solid angle 2π sr. Using step 5, Φ = (q/4πε0)(2π) = q/2ε0. Numerically Φ = (6.0×10−9)/[2(8.854×10−12)] = 339 N·m2·C−1. (This is a charge on the boundary case — note the flat base carries no flux since E is parallel to it.)
- A charge q = 8.0 nC sits at one corner of a cube. Find the flux through the whole cube.
Solution
At a corner the cube occupies one octant of the surrounding space, subtending 4π/8 = π/2 sr. Hence Φ = (q/4πε0)(π/2) = q/8ε0. Numerically Φ = (8.0×10−9)/[8(8.854×10−12)] = 113 N·m2·C−1. (Seven other cubes would be needed to surround the charge fully, together giving q/ε0.)
- A uniform charge density ρ = 1.0 μC·m−3 fills a ball of radius a = 0.20 m. Find the flux through a concentric spherical Gaussian surface of radius r = 0.10 m (inside the ball).
Solution
The enclosed charge is only that within r: Qenc = ρ·(4/3)πr3 = (1.0×10−6)·(4/3)π(0.10)3 = 4.19×10−9 C. Then Φ = Qenc/ε0 = (4.19×10−9)/(8.854×10−12) = 473 N·m2·C−1. (For any r > a the flux would saturate at the total charge ρ(4/3)πa3/ε0.)