My final year project involves searching for stable orbits for the gravitational n-body problem where we use Newton’s point mass approximation to describe the force field around masses assuming that these are approximately spherical. Now given that this model uses point masses we can’t model rotations. In fact, it’s generally assumed that they are negligible.

However, I wondered whether this assumption can break down at some point. Consider the following:

- On the scale of millenia, the surface of our planet actually behaves like a fluid.
- If the distance between two planetary bodies is small and we assume zero rotation, this would imply a growing deviation from sphericity over time. It follows that the point-mass approximation would no longer make sense.
- In fact, our bodies would look something like blobs of gooey paint slowly dripping towards each other.

While it’s true that the probability of zero rotation is close to zero the challenge of dealing with non-spherical masses is something that Newton wouldn’t have handled very well with his point-mass approximation. However, this is not quite the end of Newton’s model.

Given that the gravitational potential in his model holds for point masses a natural solution would be to integrate over the volume enclosed by a compact body. This is exactly what Gauss did when he discovered the divergence theorem from which the Gauss Law of Gravity emerges as an immediate corollary. As this handles the case where our masses aren’t spherical the rest of my blog post will focus on the statement and proof of this theorem.

**Gauss Law of Gravitation:**

Given as the gravitational potential and as the gravitational field, Gauss’ Law of Gravitation states:

In order to arrive at this result, I’ll first need to introduce two useful definitions:

**flux:**The flux of a vector field through a surface is the integral of the scalar product of the field in a given point times the infinitesimal oriented surface at that point:

where denotes the gravitational field, the boundary of our mass and is our infinitesimal oriented surface.

**divergence:**The divergence represents volume density of the outward flux of a vector field from an infinitesimal volume around a given point:

The flux and divergence of a field are actually related in a local manner. Consider a box with dimensions . Using Cartesian coordinates, let’s suppose that denotes a corner on this box. Then the flux of projected onto is given by:

For very small ,

where we get equality in the limit as tends to zero.

As a result we have :

Now, we may deduce that the flux and divergence are related as follows:

If we take an arbitrary solid and slice it into compact volumes with associated surfaces we have:

In addition, using the relationship between flux and divergence we have:

By equating the last two expressions we conclude the proof of this marvellous result by Gauss which is succinctly summarised on Wolfram as follows:

in the absence of the creation or destruction of matter, the density within a region of space can change only by having it flow into or away from the region through its boundary

At this point it’s natural to ask how can we describe the gravitational field which is defined in terms of the gravitational potential. In general, for non-spherical masses we can no longer use the usual inverse radius gravitational potential due to Newton.

In some sense we’ll need an infinite sum of perturbations in the gravitational potential which can be done using spherical harmonics, another wonderful mathematical invention that has many applications. This shall be the subject of my next blog post.