Gravitational potential

In classical mechanics, the gravitational potential at a location is equal to the work (energy transferred) per unit mass that would be needed to move an object to that location from a fixed reference location. It is analogous to the electric potential with mass playing the role of charge. The reference location, where the potential is zero, is by convention infinitely far away from any mass, resulting in a negative potential at any finite distance.

In mathematics, the gravitational potential is also known as the Newtonian potential and is fundamental in the study of potential theory. It may also be used for solving the electrostatic and magnetostatic fields generated by uniformly charged or polarized ellipsoidal bodies.^[1]

Potential energy

The gravitational potential (V) at a location is the gravitational potential energy (U) at that location per unit mass:

V={\frac {U}{m}},

where m is the mass of the object. Potential energy is equal (in magnitude, but negative) to the work done by the gravitational field moving a body to its given position in space from infinity. If the body has a mass of 1 kilogram, then the potential energy to be assigned to that body is equal to the gravitational potential. So the potential can be interpreted as the negative of the work done by the gravitational field moving a unit mass in from infinity.

In some situations, the equations can be simplified by assuming a field that is nearly independent of position. For instance, in a region close to the surface of the Earth, the gravitational acceleration, g, can be considered constant. In that case, the difference in potential energy from one height to another is, to a good approximation, linearly related to the difference in height:

\Delta U\approx mg\Delta h.

Mathematical form

The gravitational potential V at a distance x from a point mass of mass M can be defined as the work W that needs to be done by an external agent to bring a unit mass in from infinity to that point:^[2]^[3]^[4]^[5]

V(\mathbf {x} )={\frac {W}{m}}={\frac {1}{m}}\int \limits _{\infty }^{x}\mathbf {F} \cdot d\mathbf {x} ={\frac {1}{m}}\int \limits _{\infty }^{x}{\frac {GmM}{x^{2}}}dx=-{\frac {GM}{x}},

where G is the gravitational constant, and F is the gravitational force. The potential has units of energy per mass, e.g., J/kg in the MKS system. By convention, it is always negative where it is defined, and as x tends to infinity, it approaches zero.

The gravitational field, and thus the acceleration of a small body in the space around the massive object, is the negative gradient of the gravitational potential. Thus the negative of a negative gradient yields positive acceleration toward a massive object. Because the potential has no angular components, its gradient is

\mathbf {a} =-{\frac {GM}{x^{3}}}\mathbf {x} =-{\frac {GM}{x^{2}}}{\hat {\mathbf {x} }},

where x is a vector of length x pointing from the point mass toward the small body and ${\hat {\mathbf {x} }}$ is a unit vector pointing from the point mass toward the small body. The magnitude of the acceleration therefore follows an inverse square law:

|\mathbf {a} |={\frac {GM}{x^{2}}}.

The potential associated with a mass distribution is the superposition of the potentials of point masses. If the mass distribution is a finite collection of point masses, and if the point masses are located at the points x₁, ..., x_n and have masses m₁, ..., m_n, then the potential of the distribution at the point x is

V(\mathbf {x} )=\sum _{i=1}^{n}-{\frac {Gm_{i}}{|\mathbf {x} -\mathbf {x_{i}} |}}.

If the mass distribution is given as a mass measure dm on three-dimensional Euclidean space R³, then the potential is the convolution of −G/|r| with dm.^[6] In good cases this equals the integral

V(\mathbf {x} )=-\int _{\mathbf {R} ^{3}}{\frac {G}{|\mathbf {x} -\mathbf {r} |}}\,dm(\mathbf {r} ),

where |x − r| is the distance between the points x and r. If there is a function ρ(r) representing the density of the distribution at r, so that dm(r)= ρ(r)dv(r), where dv(r) is the Euclidean volume element, then the gravitational potential is the volume integral

V(\mathbf {x} )=-\int _{\mathbf {R} ^{3}}{\frac {G}{|\mathbf {x} -\mathbf {r} |}}\,\rho (\mathbf {r} )dv(\mathbf {r} ).

If V is a potential function coming from a continuous mass distribution ρ(r), then ρ can be recovered using the Laplace operator, Δ:

\rho (\mathbf {x} )={\frac {1}{4\pi G}}\Delta V(\mathbf {x} ).

This holds pointwise whenever ρ is continuous and is zero outside of a bounded set. In general, the mass measure dm can be recovered in the same way if the Laplace operator is taken in the sense of distributions. As a consequence, the gravitational potential satisfies Poisson's equation. See also Green's function for the three-variable Laplace equation and Newtonian potential.

The integral may be expressed in terms of known transcendental functions for all ellipsoidal shapes, including the symmetrical and degenerate ones.^[7] These include the sphere, where the three semiaxes are equal; the oblate (see reference ellipsoid) and prolate spheroids, where two semiaxes are equal; the degenerate ones where one semiaxis is infinite (the elliptical and circular cylinder) and the unbounded sheet where two semiaxes are infinite. All these shapes are widely used in the applications of the gravitational potential integral (apart from the constant G, with 𝜌 being a constant charge density) to electromagnetism.

Spherical symmetry

A spherically symmetric mass distribution behaves to an observer completely outside the distribution as though all of the mass was concentrated at the center, and thus effectively as a point mass, by the shell theorem. On the surface of the earth, the acceleration is given by so-called standard gravity g, approximately 9.8 m/s², although this value varies slightly with latitude and altitude: The magnitude of the acceleration is a little larger at the poles than at the equator because Earth is an oblate spheroid.

Within a spherically symmetric mass distribution, it is possible to solve Poisson's equation in spherical coordinates. Within a uniform spherical body of radius R and density ρ, the gravitational force g inside the sphere varies linearly with distance r from the center, giving the gravitational potential inside the sphere, which is^[8]

V(r)={\frac {2}{3}}\pi G\rho (r^{2}-3R^{2}),\qquad r\leq R,

which differentiably connects to the potential function for the outside of the sphere (see the figure at the top).

General relativity

In general relativity, the gravitational potential is replaced by the metric tensor. When the gravitational field is weak and the sources are moving very slowly compared to light-speed, general relativity reduces to Newtonian gravity, and the metric tensor can be expanded in terms of the gravitational potential.^[9]

Multipole expansion

The potential at a point x is given by

V(\mathbf {x} )=-\int _{\mathbb {R} ^{3}}{\frac {G}{|\mathbf {x} -\mathbf {r} |}}\ dm(\mathbf {r} ).

The potential can be expanded in a series of Legendre polynomials. Represent the points x and r as position vectors relative to the center of mass. The denominator in the integral is expressed as the square root of the square to give

{\begin{aligned}V(\mathbf {x} )&=-\int _{\mathbb {R} ^{3}}{\frac {G}{\sqrt {|\mathbf {x} |^{2}-2\mathbf {x} \cdot \mathbf {r} +|\mathbf {r} |^{2}}}}\,dm(\mathbf {r} )\\{}&=-{\frac {1}{|\mathbf {x} |}}\int _{\mathbb {R} ^{3}}G\,\left/\,{\sqrt {1-2{\frac {r}{|\mathbf {x} |}}\cos \theta +\left({\frac {r}{|\mathbf {x} |}}\right)^{2}}}\right.\,dm(\mathbf {r} )\end{aligned}}

where, in the last integral, r = |r| and θ is the angle between x and r.

The integrand can be expanded as a Taylor series in Z = r/|x|, by explicit calculation of the coefficients. A less laborious way of achieving the same result is by using the generalized binomial theorem.^[10] The resulting series is the generating function for the Legendre polynomials:

\left(1-2XZ+Z^{2}\right)^{-{\frac {1}{2}}}\ =\sum _{n=0}^{\infty }Z^{n}P_{n}(X)

valid for |X| ≤ 1 and |Z| < 1. The coefficients P_n are the Legendre polynomials of degree n. Therefore, the Taylor coefficients of the integrand are given by the Legendre polynomials in X = cos θ. So the potential can be expanded in a series that is convergent for positions x such that r < |x| for all mass elements of the system (i.e., outside a sphere, centered at the center of mass, that encloses the system):

{\begin{aligned}V(\mathbf {x} )&=-{\frac {G}{|\mathbf {x} |}}\int \sum _{n=0}^{\infty }\left({\frac {r}{|\mathbf {x} |}}\right)^{n}P_{n}(\cos \theta )\,dm(\mathbf {r} )\\{}&=-{\frac {G}{|\mathbf {x} |}}\int \left(1+\left({\frac {r}{|\mathbf {x} |}}\right)\cos \theta +\left({\frac {r}{|\mathbf {x} |}}\right)^{2}{\frac {3\cos ^{2}\theta -1}{2}}+\cdots \right)\,dm(\mathbf {r} )\end{aligned}}

The integral $\int r\cos \theta dm$ is the component of the center of mass in the x direction; this vanishes because the vector x emanates from the center of mass. So, bringing the integral under the sign of the summation gives

V(\mathbf {x} )=-{\frac {GM}{|\mathbf {x} |}}-{\frac {G}{|\mathbf {x} |}}\int \left({\frac {r}{|\mathbf {x} |}}\right)^{2}{\frac {3\cos ^{2}\theta -1}{2}}dm(\mathbf {r} )+\cdots

This shows that elongation of the body causes a lower potential in the direction of elongation, and a higher potential in perpendicular directions, compared to the potential due to a spherical mass, if we compare cases with the same distance to the center of mass. (If we compare cases with the same distance to the surface, the opposite is true.)

Numerical values

The absolute value of gravitational potential at a number of locations with regards to the gravitation from ^{[clarification needed]} the Earth, the Sun, and the Milky Way is given in the following table; i.e. an object at Earth's surface would need 60 MJ/kg to "leave" Earth's gravity field, another 900 MJ/kg to also leave the Sun's gravity field and more than 130 GJ/kg to leave the gravity field of the Milky Way. The potential is half the square of the escape velocity.

Location	W.r.t. Earth	W.r.t. Sun	W.r.t. Milky Way
Earth's surface	60 MJ/kg	900 MJ/kg	≥ 130 GJ/kg
LEO	57 MJ/kg	900 MJ/kg	≥ 130 GJ/kg
Voyager 1 (17,000 million km from Earth)	23 J/kg	8 MJ/kg	≥ 130 GJ/kg
0.1 light-year from Earth	0.4 J/kg	140 kJ/kg	≥ 130 GJ/kg

Compare the gravity at these locations.

Notes

^ Solivérez, C.E. (2016). Electrostatics and magnetostatics of polarized ellipsoidal bodies: the depolarization tensor method (1st English ed.). Free Scientific Information. ISBN 978-987-28304-0-3.
^ Marion, J.B.; Thornton, S.T. (1995). Classical Dynamics of particles and systems (4th ed.). Harcourt Brace & Company. p. 192. ISBN 0-03-097302-3.
^ Arfken, George B.; Weber, Hans J. (2005). Mathematical Methods For Physicists International Student Edition (6th ed.). Academic Press. p. 72. ISBN 978-0-08-047069-6.
^ Sang, David; Jones, Graham; Chadha, Gurinder; Woodside, Richard; Stark, Will; Gill, Aidan (2014). Cambridge International AS and A Level Physics Coursebook (illustrated ed.). Cambridge University Press. p. 276. ISBN 978-1-107-69769-0.
^ Muncaster, Roger (1993). A-level Physics (illustrated ed.). Nelson Thornes. p. 106. ISBN 978-0-7487-1584-8.
^ Vladimirov 1984, §7.8
^ MacMillan, W.D. (1958). The Theory of the Potential. Dover Press.
^ Marion & Thornton 2003, §5.2
^ Grøn, Øyvind; Hervik, Sigbjorn (2007), Einstein's General Theory of Relativity: With Modern Applications in Cosmology, Springer Science & Business Media, p. 201, ISBN 978-0-387-69200-5
^ Wylie, C. R., Jr. (1960). Advanced Engineering Mathematics (2nd ed.). New York: McGraw-Hill. p. 454 [Theorem 2, Section 10.8].{{cite book}}: CS1 maint: multiple names: authors list (link)

References

Vladimirov, V. S. (1971), Equations of mathematical physics, Translated from the Russian by Audrey Littlewood. Edited by Alan Jeffrey. Pure and Applied Mathematics, vol. 3, New York: Marcel Dekker Inc., MR 0268497.
Wang, W. X. (1988). "The potential for a homogeneous spheroid in a spheroidal coordinate system. I. At an exterior point". J. Phys. A: Math. Gen. 21: 4245-4250. Bibcode:1988JPhA...21.4245W.
Milon, T. (1990). "A note on the potential of a homogenous ellipsoid in ellipsoidal coordinates". J. Phys. A: Math. Gen. 23: 581–584. doi:10.1088/0305-4470/23/4/027.
Rastall, Peter (1991). Postprincipia: Gravitation for Physicists and Astronomers. World Scientific. pp. 7ff. ISBN 981-02-0778-6.
Template:Cite article
Template:Cite article
Thornton, Stephen T.; Marion, Jerry B. (2003), Classical Dynamics of Particles and Systems (5th ed.), Brooks Cole, ISBN 978-0-534-40896-1.
Zhu, Lupeia (1988). "Gravity and Earth's Density Structure". Department of Earth and Atmospheric Sciences. EAS-437 Earth Dynamics. Saint Louis University. California Institute of Technology. Retrieved 2009-03-25.
Charles D. Ghilani (2006-11-28). "The Gravity Field of the Earth". The Physics Fact Book. Penn State Surveying Engineering Program. Archived from the original on 2011-07-18. Retrieved 2009-03-25.
Fukushima, Toshio (2014). "Prolate spheroidal harmonic expansion of gravitational field". Astrophys. J. 147 (6): 152. Bibcode:2014AJ....147..152F. doi:10.1088/0004-6256/147/6/152.

[1] Solivérez, C.E. (2016). Electrostatics and magnetostatics of polarized ellipsoidal bodies: the depolarization tensor method (1st English ed.). Free Scientific Information. ISBN 978-987-28304-0-3.

[2] Marion, J.B.; Thornton, S.T. (1995). Classical Dynamics of particles and systems (4th ed.). Harcourt Brace & Company. p. 192. ISBN 0-03-097302-3.

[3] Arfken, George B.; Weber, Hans J. (2005). Mathematical Methods For Physicists International Student Edition (6th ed.). Academic Press. p. 72. ISBN 978-0-08-047069-6.

[4] Sang, David; Jones, Graham; Chadha, Gurinder; Woodside, Richard; Stark, Will; Gill, Aidan (2014). Cambridge International AS and A Level Physics Coursebook (illustrated ed.). Cambridge University Press. p. 276. ISBN 978-1-107-69769-0.

[5] Muncaster, Roger (1993). A-level Physics (illustrated ed.). Nelson Thornes. p. 106. ISBN 978-0-7487-1584-8.

[6] Vladimirov 1984, §7.8

[7] MacMillan, W.D. (1958). The Theory of the Potential. Dover Press.

[8] Marion & Thornton 2003, §5.2

[Newtonian_or_gravitoelectric_potential-9] Grøn, Øyvind; Hervik, Sigbjorn (2007), Einstein's General Theory of Relativity: With Modern Applications in Cosmology, Springer Science & Business Media, p. 201, ISBN 978-0-387-69200-5

[AEM-10] Wylie, C. R., Jr. (1960). Advanced Engineering Mathematics (2nd ed.). New York: McGraw-Hill. p. 454 [Theorem 2, Section 10.8].{{cite book}}: CS1 maint: multiple names: authors list (link)

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]