Since my blog claims to talk about physics sometimes and I just finished teaching multivariable calculus, I thought I’d do a post on one form of Gauss’ law. As a teacher of the course, I found this to be an astonishingly beautiful “application” of the divergence theorem. It turned out to be a touch too difficult for my students (and I vaguely recall being extremely confused about this when I took the class myself).

First, I’ll remind you what some of this stuff is if you haven’t thought about these concepts for awhile. Let’s work in for simplicity. Consider some subset . Let be a vector field. Mathematically this is just assigning a vector to each point of . For calculus we usually put some fairly restrictive conditions on , such as all partial derivatives exist and are continuous.

The above situation is ubiquitous in classical physics. The vector field could be the gravitational field or the electric field or it could describe velocity of a flowing fluid or … One key quantity you might want to know about your field is what is the flux of the field through a given surface ? This measures the net change of the field flowing through the surface. If is just a sphere, then it is easy to visualize the flux as the amount leaving the sphere minus the amount flowing in.

Let’s suppose is a smooth surface bounding a solid volume (e.g. the sphere bounding the solid ball). In this case we have a well-defined “outward normal” direction. Define to be the unit vector field in this direction at all points of . Just by definition the flux of through must be “adding up” the values of over , because this dot product just tells us how much is pointing in the outward direction.

Thus we define the flux (using Stewart’s notation) to be:

Note the second integral is integrating a scalar valued function with respect to surface area “dS.” Now recall that the divergence theorem says that in our situation (given that extends to a vector field on an open set containing ) we can calculate this rather tedious surface integral by converting it to a usual triple integral:

If you’re advanced, then of course you could just work this out as a special case of Stoke’s theorem using the musical isomorphisms and so on. Let’s now return to our original problem. Suppose I have a charge inside some surface and I want to compute the flux of the associated electric field through .

From my given information this would seem absolutely impossible. If can be anything, and can be located anywhere inside, then of course there are just way too many variables to come up with a reasonably succinct answer. Surprisingly, Gauss’ law tells us that no matter what is and where is located, the answer is always the same, and it is just a quick application of the divergence theorem to prove it.

First, let’s translate everything so that is located at the origin. Since flux is translation invariant, this will not change our answer. We first need to know what the electric field is, and this is essentially a direct consequence of Coloumb’s law:

If we care about higher dimensions, then we might want to note that the value only depends on the radial distance from the origin and write it in the more succinct way , where is just some constant that depends on the textbook/units you are working in. Let’s first compute the partial of the first coordinate with respect to (ignoring the constant factor for now):

You get similar things for taking the other derivatives involved in the divergence except the minus sign moves to and respectively. When you add all these together you get in the numerator . Thus the divergence is everywhere and hence by the divergence theorem the flux must be too, right? Wrong! And that’s where I lost most of my students.

Recall that pesky hypothesis that can be extended to a vector field on an open neighborhood of . Our can’t even be defined at all to extend continuously across the origin. Thus we must do something different. Here’s the idea, we just change our region . Since is open and contains the origin, we can find a small sphere of radius and centered at whose interior is properly contained in , say .

Let be the region between these two surfaces. Effectively this “cuts out” the bad point of and now we are allowed to apply the divergence theorem to where our new boundary is oriented outwards and oriented inward (negatively). We already calculated that , thus one side of the equation is . This gives us

This is odd, because it says that no matter how bizarre or gigantic was we can just compute the flux through a small sphere and get the same answer. At this point we’ve converted the problem to something we can do because the unit normal is just . Direct computation gives us

Plugging this all in we get that the flux through is

That’s Gauss’ Law. It says that no matter the shape of or the location of the charge inside , you can always compute the flux of the electric field produced by through as a constant multiple of the amount of charge! In fact, most books use where $\varepsilon_0$ is the “permittivity of free space” which kills off practically all extraneous symbols in the answer.