# Does the gravitational attraction near the surface of dense celestial objects diverge from inverse square?

Does the gravitational attraction near the surface of dense celestial objects (neutron star, white dwarf itc) diverge (to infinity) from inverse square?

This question is inspired by the similarity between EM and gravity (inverse square force.). The paper by John Lekner here (doi:10.1098/rspa.2012.0133) shows that there is an electrostatic attraction between charged spheres no matter the polarities of the charges and that it diverges at close separation until electrical short for almost all charge ratios. I am wondering if there is a similar kind of gravitational inverse square divergence for anything other than a black hole.

Actually make this for a black hole as well, although I know that a black hole is not thought to have a normal surface.

• Lekner's paper deals with spheres where charges redistribute themselves on the surface. The corresponding gravity thing would redistribute mass and deform the objects. That definitely complicates things compared to 1/r^2 – Anders Sandberg Mar 15 '20 at 0:12
• I'm not expert enough to answer properly, but I think it does, due to general relativity, but to describe how it does, you need to first define your radius and time coordinates, and there is no single particularly privileged choice. – Steve Linton Mar 15 '20 at 1:22
• Classical departure from inverse square happens when the "charges" (electrostatic or otherwise) are no longer spherically symmetric and Newton's Shell theorem no longer applies. Any time one or both gravitational bodies depart from spherical symmetric mass distribution the force between them can no longer be assumed to be inverse square. The bodies may naturally depart from spherical symmetry due to rotation, but they can also induce it. There will be a very tiny effect on the Moon's orbit due to the the motion of Earth's ocean tides it induces. – uhoh Mar 15 '20 at 1:25
• @uhoh I clarified the question above with the two words (to infinity). – DMac Mar 16 '20 at 17:24

Black holes produce another complication: since spacetime nearby is curved and expanded the meaning of the distance in the inverse square law becomes problematic. The Paczyński–Wiita potential is an approximation of the potential, and it deviates from the $$U=-GM/r$$ as $$U_{PW}=-GM/(r-R_S)$$ (where $$R_S$$ is the Schwarzschild radius). It makes the force increase faster than the classical potential as we approach $$r=R_S$$.
The result is indeed that the force increases faster than $$1/r^2$$ as the bodies approach each other, since they elongate and eventually merge (a bit before this they will deviate from my ellipsoidal assumption). If one multiplies the force by the squared distance the product should be constant for pure $$1/r^2$$ forces, but it starts to increase as they approach enough. Note that this is a non-rotating model: with rotation the numbers will change and the ellipsoids will become tri-axial, but I suspect the qualitative behaviour stays the same.