Gauss's Law: Differential Form and Physical Meaning

Hello! Welcome back to our study of electrostatics.

In the last lesson, you mastered Gauss's Law in its integral form, $\oint_{S} E \cdot d S = Q_{e n c} / ϵ_{0}$ . We saw how this "global" law, which relates the electric field over a closed surface to the total charge inside that surface, provides a powerful method for calculating electric fields in situations with high symmetry.

Today, we delve deeper to uncover the "local" version of this law. Our learning outcome is to derive the differential form of Gauss's Law ( $\nabla \cdot E = ρ / ϵ_{0}$ ) and understand its physical meaning. This equation is arguably even more fundamental than the integral form. It tells us how charge at a single point in space generates the electric field at that very same point. This local perspective is essential for tackling more complex problems where symmetry is absent and for building the complete theory of electromagnetism.

We will explore two complementary paths to this result. The first is a swift, elegant derivation using the Divergence Theorem you've already mastered. The second is a more physically intuitive derivation from first principles, which will build a concrete understanding of what the "divergence" of a field truly represents.

1. The Mathematical Path: Applying the Divergence Theorem

Let's begin with the most direct route from our last lesson. This approach showcases the power of vector calculus as a language for physics.

We start with Gauss's Law in integral form:
$\oint_{S} E \cdot d S = \frac{Q _{e n c}}{ϵ _{0}}$

The total charge enclosed, $Q_{e n c}$ , can be expressed as the volume integral of the charge density, $ρ (r)$ , over the volume $V$ enclosed by the surface $S$ :
$Q_{e n c} = ∭_{V} ρ (r) d V$

Substituting this into Gauss's Law gives:
$\oint_{S} E \cdot d S = \frac{1}{ϵ _{0}} ∭_{V} ρ (r) d V$

Now, we invoke the Divergence Theorem for the left-hand side of the equation, which connects the surface integral of a vector field to the volume integral of its divergence:
$\oint_{S} E \cdot d S = ∭_{V} (\nabla \cdot E) d V$

By equating these two expressions for the flux, we arrive at a single equation involving volume integrals:
$∭_{V} (\nabla \cdot E) d V = ∭_{V} \frac{ρ ( r )}{ϵ _{0}} d V$
$⟹ ∭_{V} (\nabla \cdot E - \frac{ρ}{ϵ _{0}}) d V = 0$

This equation must hold true for any volume $V$ we choose, no matter how large or small. The only way for an integral to be zero for every possible integration domain is if the integrand itself is zero everywhere. This powerful argument leads us directly to our destination.

Gauss' Law - Differential Form

The following resource provides a concise summary of this elegant derivation using the Divergence Theorem.

Please read the short section titled "Derivation via the Divergence Theorem". It presents the exact logical steps we just discussed.

This leaves us with the first of Maxwell's equations in differential form:
$\nabla \cdot E = \frac{ρ}{ϵ _{0}}$

2. The Physical Path: Flux from an Infinitesimal Volume

The derivation above is efficient, but it might feel a bit abstract. To build a solid physical intuition for what divergence means, let's derive the same result from scratch by considering the flux out of an infinitesimally small box.

Imagine a tiny rectangular box at a point $(x, y, z)$ with side lengths $Δ x$ , $Δ y$ , and $Δ z$ . We will calculate the total electric flux out of its six faces.

The key insight, which your calculus background makes you well-equipped to handle, is that for a very small $Δ x$ , the change in the electric field can be approximated using a partial derivative:
$E_{x} (x + Δ x) \approx E_{x} (x) + \frac{\partial E _{x}}{\partial x} Δ x$

By applying this idea to pairs of opposite faces on the box, we can calculate the net flux and relate it to the charge inside.

Lecture 26 Maxwell Equations - The Full Story

This video provides an excellent, step-by-step walkthrough of this derivation. The professor carefully calculates the flux through each pair of faces and assembles the result, providing a clear visual and conceptual foundation for the divergence operator.

Please watch from 19:04 to 30:01. Focus on how the net flux through opposite faces is related to the partial derivative of the field component perpendicular to those faces. This process physically constructs the divergence term.

As the video demonstrates, the total flux $Φ_{E}$ out of the box is:
$Φ_{E} \approx (\frac{\partial E _{x}}{\partial x} + \frac{\partial E _{y}}{\partial y} + \frac{\partial E _{z}}{\partial z}) Δ x Δ y Δ z$

The term in the parentheses is the definition of the divergence of E, written as $\nabla \cdot E$ , in Cartesian coordinates. The product of the side lengths is the volume of the box, $Δ V$ . So, $Φ_{E} \approx (\nabla \cdot E) Δ V$ .

From Gauss's integral law, we also know $Φ_{E} = Q_{e n c} / ϵ_{0}$ . For our tiny box, the enclosed charge is simply $Q_{e n c} = ρ Δ V$ . Equating our two expressions for flux gives:
$(\nabla \cdot E) Δ V = \frac{ρ Δ V}{ϵ _{0}}$
$\nabla \cdot E = \frac{ρ}{ϵ _{0}}$
We have arrived at the same result, but this time with a clear physical picture: the divergence of the electric field at a point is the net flux flowing out of an infinitesimal volume at that point, per unit volume.

3. The Physical Meaning of Divergence

So, what does $\nabla \cdot E = ρ / ϵ_{0}$ really tell us?

If $ρ > 0$ at a point, then $\nabla \cdot E > 0$ . This means there is a net outward flux from that point. Electric field lines "originate" or "diverge" from positive charges. The point is a source of the field.
If $ρ < 0$ at a point, then $\nabla \cdot E < 0$ . There is a net inward flux. Electric field lines "terminate" or "converge" on negative charges. The point is a sink of the field.
If $ρ = 0$ at a point (in empty space), then $\nabla \cdot E = 0$ . There is no net flux originating or terminating. Field lines just pass through that point without starting or stopping.

This concept of "sourciness" is the core physical meaning of divergence.

ESF04: Differential Form of Gauss's Law for Electric Field

For a memorable and intuitive explanation of this idea, watch the following segment.

Watch from 14:37 to 21:48. The instructor introduces the "sourciness" analogy, which is an extremely effective way to visualize and remember what the divergence of a vector field represents.

4. Divergence in Other Coordinate Systems

The expression $\nabla \cdot E = \frac{\partial E _{x}}{\partial x} + \frac{\partial E _{y}}{\partial y} + \frac{\partial E _{z}}{\partial z}$ is only valid in Cartesian coordinates. Since many problems in physics have spherical or cylindrical symmetry, you must use the correct form of the divergence operator for those systems.

Electromagnetism I: Electrostatics

This Olympiad-level handout gives the formula for divergence in spherical coordinates and then immediately applies it to an important problem, checking for consistency between the integral and differential forms of Gauss's Law.

Read the section for "Example 4" (page 4). Note the formula provided for divergence in spherical coordinates. Then, follow the example to see how the field of a spherically symmetric charge distribution, which you found using the integral law in the last lesson, is shown to satisfy the differential law.

For reference, here are the divergence formulas:

Cylindrical $(ρ, ϕ, z)$ :
$\nabla \cdot F = \frac{1}{ρ} \frac{\partial}{\partial ρ} (ρ F_{ρ}) + \frac{1}{ρ} \frac{\partial F _{ϕ}}{\partial ϕ} + \frac{\partial F _{z}}{\partial z}$
Spherical $(r, θ, ϕ)$ :
$\nabla \cdot F = \frac{1}{r ^{2}} \frac{\partial}{\partial r} (r^{2} F_{r}) + \frac{1}{r s i n θ} \frac{\partial}{\partial θ} (F_{θ} sin θ) + \frac{1}{r s i n θ} \frac{\partial F _{ϕ}}{\partial ϕ}$

A common mistake is to naively apply the Cartesian formula to the components in other coordinate systems. This is incorrect because the unit vectors themselves ( $r^$ , $θ^$ , etc.) change direction from point to point, a complexity that these formulas correctly account for.

Test your understanding!

Suppose the electric field in a region of space is given in cylindrical coordinates by $E = k ρ^{2} ρ^$ , where $k$ is a constant. What is the volume charge density $ρ_{c ha r g e} (ρ, ϕ, z)$ in this region?

Show answer

We use the differential form of Gauss's Law, $ρ_{c ha r g e} = ϵ_{0} (\nabla \cdot E)$ , with the divergence operator in cylindrical coordinates.
Here, $E_{ρ} = k ρ^{2}$ , and $E_{ϕ} = E_{z} = 0$ .

$\nabla \cdot E = \frac{1}{ρ} \frac{\partial}{\partial ρ} (ρ E_{ρ}) + 0 + 0$
$\nabla \cdot E = \frac{1}{ρ} \frac{\partial}{\partial ρ} (ρ \cdot k ρ^{2}) = \frac{1}{ρ} \frac{\partial}{\partial ρ} (k ρ^{3})$
$\nabla \cdot E = \frac{1}{ρ} (3 k ρ^{2}) = 3 k ρ$

Therefore, the charge density is $ρ_{c ha r g e} = 3 k ϵ_{0} ρ$ . The charge density increases linearly with the distance from the z-axis.

5. A Deeper Look: The Role of Boundary Conditions

Your goal is to delve into theoretical physics, which often involves scrutinizing the assumptions behind our equations. The differential form of Gauss's Law is a first-order partial differential equation. As with any differential equation, its solution (the electric field $E$ ) is not unique without specifying boundary conditions.

Consider a universe filled uniformly with charge density $ρ_{0}$ . What is the electric field?

A symmetry argument might suggest $E = 0$ , since there's no preferred direction. But this gives $\nabla \cdot E = 0$ , which violates Gauss's Law if $ρ_{0} \neq = 0$ .
The field $E = (ρ_{0} /3 ϵ_{0}) r$ works mathematically ( $\nabla \cdot E = ρ_{0} / ϵ_{0}$ ).
The field $E = (ρ_{0} / ϵ_{0}) x x^$ also works!

Which is correct? The question is ill-posed because for an infinite charge distribution, we cannot apply the usual boundary condition that $E \to 0$ at infinity. The answer depends on how you imagine "building" this infinite distribution, which is a way of imposing a symmetry that isn't required by the local physics.

Electromagnetism I: Electrostatics

This is a subtle but profound point. For a challenging look into this issue, which puzzled Newton himself in the gravitational analogue, read the following problem.

Read 'Problem 6' and its solution on pages 5-6. This will challenge you to think about the limitations of physical laws and the crucial role of boundary conditions and symmetry assumptions, a key theme in advanced theoretical physics.

Conclusion

In this lesson, we have transformed the global, integral form of Gauss's Law into its local, differential form, which is the first of the four fundamental Maxwell's Equations.

Key Takeaways:

Gauss's Law in Differential Form: $\nabla \cdot E = ρ / ϵ_{0}$ . This is a local law relating the structure of the electric field at a point to the charge density at that same point.
Two Paths to Derivation: We showed it can be derived elegantly using the Divergence Theorem or more physically by considering the flux out of an infinitesimal volume.
Physical Meaning of Divergence: $\nabla \cdot E$ represents the "sourciness" of the electric field. It measures the net outflow of flux per unit volume. Positive charges are sources ( $\nabla \cdot E > 0$ ), and negative charges are sinks ( $\nabla \cdot E < 0$ ).
Coordinate Systems Matter: The mathematical expression for the divergence operator is different in Cartesian, cylindrical, and spherical coordinates. Using the correct form is crucial for problem-solving.

Preview of the Next Lesson:
We have now established the "divergence" of the electrostatic field. But what about its "curl" ( $\nabla \times E$ )? In the next lesson, we will investigate this and discover another of Maxwell's fundamental equations: $\nabla \times E = 0$ . We will see that this property means the electrostatic field is conservative, which allows us to define a scalar quantity that is often much easier to work with than the electric field itself: the electric potential.

Recommended resources