The Yukawa and Exponential Potentials for a Galaxy’s Spherical Central Bulge* ()
1. Introduction
In Book III of Principia, Newton [1] published proofs of two theorems that afford an easy calculation of the gravitational potential of any spherically symmetric distribution of matter and to calculate excellent approximations of the potential for approximately spherical distributions of matter such as the central bulge of a spiral galaxy. Thus, since 1687 we have known that the Newtonian acceleration of an orbiting object depends only on the mass inside the orbit. Any uniform spherical shell of matter outside the orbit has no effect on the orbital dynamics. These useful Newtonian theorems are not valid for other gravitational potentials proposed over the past decade to explain anomalous galactic dynamics. For example, with a linear potential [2] , the centripetal acceleration of a disk star is sensitive to the halo mass outside its orbit and, even more inconveniently, to the mass of the entire Universe.
The foci of our exposition are two other non-Newtonian gravitational models that are more tractable than the linear potential. The first is a Yukawa potential [3] which, for a point mass M at
, has the form
where G is the Newtonian gravitational constant,
is a dimensionless constant, and
is the wave number (
is the Compton wavelength). The second is an exponential potential [4] [5] which is actually the superposition of two Yukawa potentials with coupling constants
that are identical except for their signs: an attractive potential,
, and a repulsive potential,
. Their wave numbers are different, but almost identical:
and
. Then
(1)
In this paper, we present the mathematical development that has allowed us to resolve for the exponential potential the same problem that Newton solved for the Newtonian potential. That is, given the mass distribution of a spherically symmetric central bulge and halo, then what are their contributions to the specific forces on disk matter that are due to the exponential potential? With this distribution, how much does the exponential potential contribute to the orbital velocity of disk matter? In fact, we provide a route to the potential, and then to the acceleration, that is often more convenient than the iteration by Eckhardt [5] . This model assumes that the central bulge of a galaxy is spherically symmetric, that the mass of the galaxy disk is negligible or is treated separately [6] , and that there is an exponential potential as well as a Newtonian potential. To make it clear, the same approach is used to model the shell for the Newtonian potential, then the Yukawa potential, and finally the exponential potential. Although the solution for the Newtonian potential is trivial, it is presented as a guide for the approach for solving the other potentials.
2. Newtonian Potential
The governing equation for the [1] potential
, which is mediated by the massless graviton, is Poisson’s equation:
The potential of a point source with mass
at a distance R from the source is
. (That is, the appropriate Green’s function that vanishes at
is
.) Thus the potential at
of a shell at
with mass
is
Because of the spherical symmetry, Laplace’s Equation (
) is
The two solutions to this equation are of the form
and
. The solution that remains finite at
is
, so for
, we have
(2)
The solution that remains finite at
is
, so for
, we have
(3)
3. Differential Equation for the Exponential Potential
We define the operators
and
With spherical symmetry, the differential equations for the Yukawa potentials
and
are then
and
In addition to
, we define
. We then have
and
It is obvious that
, so we neglect the
term and then define
and finally derive
(4)
and
(5)
The exponential potential satisfies the fourth order differential equation
4. The Yukawa and Exponential Potentials of a Thin Spherical Shell
For a thin spherical shell of mass M at
, the Yukawa potential (coupling constant
) at the origin is
In free space, the Yukawa potential satisfies the homogeneous equation
, whose solutions are linear combinations of
. The only combination that is finite and equal to unity at the origin is
, so the Yukawa potential inside the shell is
(6)
The solution interval is closed at the shell because the potential (but not its derivative) is continuous there. Thus
The homogeneous solution that is equal to unity at the shell boundary and remains finite outside the shell is
, so the Yukawa potential outside the shell is
(7)
Note that Equation (6) becomes Equation (7) on interchanging
and
.
The derivative of
is discontinuous at
. Let
be the shell thickness, and let
represent
just inside the shell, while
represents r just outside the shell. Then
For
, this is the Gauss divergence theorem.
Setting
in Equation (1), the exponential potential is
(8)
and, on interchanging
and
and then rearranging terms,
(9)
The operators
and
commute, so
Thus, unlike
, the derivative of
is continuous at
. This is the reason why
cannot be modeled with a second order differential equation. Note that in the far field, where
is negligible, Equation (9) becomes
and, for
,
All the variables in Equations (8) and (9) are linear combinations of
and
. Thus, because
and
the equation
is satisfied wherever
.
5. Specific Forces and Galaxy Rotation Velocities
The centripetal specific force (acceleration) due to the Newtonian and exponential potentials of the spherically symmetric bulge is
where
From Equations (3) and (9),
and from Equations (2) and (8),
Thus
The sum of the two last terms (dominated by the last one), which is negative, represents a very small outward force due to the exponential potential of matter further from the origin than r; they vanish at the origin. In modeling spiral galaxies with a combination of Newtonian and exponential potentials [7] , we have estimated that
and
, roughly one order of magnitude larger than the radius of the central bulge, so a good approximation is
where
If a star is in a circular orbit a distance
from the origin with velocity
, then
. Setting
, we then have
6. Conclusions
Whereas the Poisson equation for a spherically symmetric Newtonian potential is of second order in
, the corresponding equation for the exponential potential is of fourth order. We have derived rigorous and precise mathematical solutions to calculate the exponential potential of a spherically symmetric distribution of mass such as a galactic bulge or halo. Moreover, we have derived a simple and efficient algorithm to closely approximate the associated radial gradients in galaxies. Using the Abel transform, all integrations of spherically symmetric models are relatively simple.
Combined with our technique for modeling the gravitational field of a galactic disk [6] , we have found it easy to use in evaluating the effects of Newtonian and exponential potentials for numerous spiral galaxies [7] . Furthermore, it is evident that our approach will be convenient to apply at other distance scales, smaller and larger, such as globular clusters, elliptical galaxies and galaxy clusters.
NOTES
*This article is dedicated to the fond memory of Don Eckhardt, who passed away in the morning of 21st October 2023. He was my friend and collaborator, and a true visionary.