Gauss's law for gravity

The two forms of Gauss's law for gravity are mathematically equivalent. The divergence theorem states:

${\displaystyle \oint _{\partial V}\mathbf {g} \cdot d\mathbf {A} =\int _{V}\nabla \cdot \mathbf {g} \ dV}$

where V is a closed region bounded by a simple closed oriented surface ∂V and dV is an infinitesimal piece of the volume V (see volume integral for more details). The gravitational field g must be a continuously differentiable vector field defined on a neighborhood of V.

Given also that

${\displaystyle M=\int _{V}\rho \ dV}$

we can apply the divergence theorem to the integral form of Gauss's law for gravity, which becomes:

${\displaystyle \int _{V}\nabla \cdot \mathbf {g} \ dV=-4\pi G\int _{V}\rho \ dV}$

which can be rewritten:

${\displaystyle \int _{V}(\nabla \cdot \mathbf {g} )\ dV=\int _{V}(-4\pi G\rho )\ dV.}$

This has to hold simultaneously for every possible volume V; the only way this can happen is if the integrands are equal. Hence we arrive at

${\displaystyle \nabla \cdot \mathbf {g} =-4\pi G\rho ,}$

which is the differential form of Gauss's law for gravity.

It is possible to derive the integral form from the differential form using the reverse of this method.

Although the two forms are equivalent, one or the other might be more convenient to use in a particular computation.

Gauss's law for gravity can be derived from Newton's law of universal gravitation, which states that the gravitational field due to a point mass is:

${\displaystyle \mathbf {g} (\mathbf {r} )=-GM{\frac {\mathbf {e_{r}} }{r^{2}}}}$

where

e_r is the radial unit vector,

r is the radius, |r|.

M is the mass of the particle, which is assumed to be a point mass located at the origin.

A proof using vector calculus is shown in the box below. It is mathematically identical to the proof of Gauss's law (in electrostatics) starting from Coulomb's law.[1]

Outline of proof: (Click [show] button to the right.)

g(r), the gravitational field at r, can be calculated by adding up the contribution to g(r) due to every bit of mass in the universe (see superposition principle). To do this, we integrate over every point s in space, adding up the contribution to g(r) associated with the mass (if any) at s, where this contribution is calculated by Newton's law. The result is: ${\displaystyle \mathbf {g} (\mathbf {r} )=-G\int \rho (\mathbf {s} ){\frac {(\mathbf {r} -\mathbf {s} )}{|\mathbf {r} -\mathbf {s} |^{3}}}d^{3}\mathbf {s} .}$ (d3s stands for dsxdsydsz, each of which is integrated from -∞ to +∞.) If we take the divergence of both sides of this equation with respect to r, and use the known theorem[1] ${\displaystyle \nabla \cdot \left({\frac {\mathbf {r} }{|\mathbf {r} |^{3}}}\right)=4\pi \delta (\mathbf {r} )}$ where δ(r) is the Dirac delta function, the result is ${\displaystyle \nabla \cdot \mathbf {g} (\mathbf {r} )=-4\pi G\int \rho (\mathbf {s} )\ \delta (\mathbf {r} -\mathbf {s} )\ d^{3}\mathbf {s} .}$ Using the "sifting property" of the Dirac delta function, we arrive at ${\displaystyle \nabla \cdot \mathbf {g} (\mathbf {r} )=-4\pi G\rho (\mathbf {r} )}$ which is the differential form of Gauss's law for gravity, as desired.

It is impossible to mathematically prove Newton's law from Gauss's law alone, because Gauss's law specifies the divergence of g but does not contain any information regarding the curl of g (see Helmholtz decomposition). In addition to Gauss's law, the assumption is used that g is irrotational (has zero curl), as gravity is a conservative force:

${\displaystyle \nabla \times \mathbf {g} =0}$

Even these are not enough: Boundary conditions on g are also necessary to prove Newton's law, such as the assumption that the field is zero infinitely far from a mass.

The proof of Newton's law from these assumptions is as follows:

Outline of proof

Start with the integral form of Gauss's law: ${\displaystyle \oint _{\partial V}\mathbf {g} \cdot d\mathbf {A} =-4\pi GM.}$ Apply this law to the situation where the volume V is a sphere of radius r centered on a point-mass M. It's reasonable to expect the gravitational field from a point mass to be spherically symmetric. (We omit the proof for simplicity.) By making this assumption, g takes the following form: ${\displaystyle \mathbf {g} (\mathbf {r} )=g(r)\mathbf {e_{r}} }$ (i.e., the direction of g is parallel to the direction of r, and the magnitude of g depends only on the magnitude, not direction, of r). Plugging this in, and using the fact that ∂V is a spherical surface with constant r and area ${\displaystyle 4\pi r^{2}}$, ${\displaystyle g(r)\oint _{\partial V}\mathbf {e_{r}} \cdot d\mathbf {A} =-4\pi GM}$ ${\displaystyle g(r)(4\pi r^{2})=-4\pi GM}$ ${\displaystyle g(r)=-GM/r^{2}}$ ${\displaystyle \mathbf {g} (\mathbf {r} )=-GM{\frac {\mathbf {e_{r}} }{r^{2}}}}$ which is Newton's law.

We can conclude (by using a "Gaussian pillbox") that for an infinite, flat plate (Bouguer plate) of any finite thickness, the gravitational field outside the plate is perpendicular to the plate, towards it, with magnitude 2πG times the mass per unit area, independent of the distance to the plate[2] (see also gravity anomalies).

More generally, for a mass distribution with the density depending on one Cartesian coordinate z only, gravity for any z is 2πG times the difference in mass per unit area on either side of this z value.

In particular, a parallel combination of two parallel infinite plates of equal mass per unit area produces no gravitational field between them.

In the case of an infinite uniform (in z) cylindrically symmetric mass distribution we can conclude (by using a cylindrical Gaussian surface) that the field strength at a distance r from the center is inward with a magnitude of 2G/r times the total mass per unit length at a smaller distance (from the axis), regardless of any masses at a larger distance.

For example, inside an infinite uniform hollow cylinder, the field is zero.

In the case of a spherically symmetric mass distribution we can conclude (by using a spherical Gaussian surface) that the field strength at a distance r from the center is inward with a magnitude of G/r² times only the total mass within a smaller distance than r. All the mass at a greater distance than r from the center has no resultant effect.

For example, a hollow sphere does not produce any net gravity inside. The gravitational field inside is the same as if the hollow sphere were not there (i.e. the resultant field is that of all masses not including the sphere, which can be inside and outside the sphere).

Although this follows in one or two lines of algebra from Gauss's law for gravity, it took Isaac Newton several pages of cumbersome calculus to derive it directly using his law of gravity; see the article shell theorem for this direct derivation.