A Novel Potential Featuring Off-Center Circular Orbits
In Book 1, Proposition 7, Problem 2 of his 1687 Philosophiae Naturalis Principia Mathematica, Isaac Newton poses and answers the following question: Let the orbit of a particle moving in a central force field be an off-center circle. How does the magnitude of the force depend on the position of the...
Saved in:
| Published in: | Symmetry, Integrability and Geometry: Methods and Applications |
|---|---|
| Date: | 2023 |
| Main Author: | |
| Format: | Article |
| Language: | English |
| Published: |
Інститут математики НАН України
2023
|
| Online Access: | https://nasplib.isofts.kiev.ua/handle/123456789/211883 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Journal Title: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Cite this: | A Novel Potential Featuring Off-Center Circular Orbits, Maxim Olshanii, SIGMA 19 (2023), 001, 8 pages |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1859825407615827968 |
|---|---|
| author | Olshanii, Maxim |
| author_facet | Olshanii, Maxim |
| citation_txt | A Novel Potential Featuring Off-Center Circular Orbits, Maxim Olshanii, SIGMA 19 (2023), 001, 8 pages |
| collection | DSpace DC |
| container_title | Symmetry, Integrability and Geometry: Methods and Applications |
| description | In Book 1, Proposition 7, Problem 2 of his 1687 Philosophiae Naturalis Principia Mathematica, Isaac Newton poses and answers the following question: Let the orbit of a particle moving in a central force field be an off-center circle. How does the magnitude of the force depend on the position of the particle onthat circle? In this article, we identify a potential that can produce such a force, only at zero energy. We further map the zero-energy orbits in this potential to finite-energy free motion orbits on a sphere; such a duality is a particular instance of a general result by Goursat, from 1887. The map itself is an inverse stereographic projection, and this fact explains the circularity of the zero-energy orbits in the system of interest. Finally, we identify an additional integral of motion—an analogue of the Runge-Lenz vector in the Coulomb problem—that is responsible for the closeness of the zero-energy orbits in our problem.
|
| first_indexed | 2026-03-16T13:28:44Z |
| format | Article |
| fulltext |
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 19 (2023), 001, 8 pages
A Novel Potential Featuring Off-Center Circular
Orbits
Maxim OLSHANII
Department of Physics, University of Massachusetts Boston, Boston Massachusetts 02125, USA
E-mail: maxim.olchanyi@umb.edu
URL: https://sites.google.com/view/integrability/home
Received October 04, 2022, in final form January 03, 2023; Published online January 07, 2023
https://doi.org/10.3842/SIGMA.2023.001
Abstract. In Book 1, Proposition 7, Problem 2 of his 1687 Philosophiæ Naturalis Principia
Mathematica, Isaac Newton poses and answers the following question: Let the orbit of a par-
ticle moving in a central force field be an off-center circle. How does the magnitude of the
force depend on the position of the particle on that circle? In this article, we identify a po-
tential that can produce such a force, only at zero energy. We further map the zero-energy
orbits in this potential to finite-energy free motion orbits on a sphere; such a duality is
a particular instance of a general result by Goursat, from 1887. The map itself is an inverse
stereographic projection, and this fact explains the circularity of the zero-energy orbits in
the system of interest. Finally, we identify an additional integral of motion—an analogue
of the Runge–Lenz vector in the Coulomb problem—that is responsible for the closeness of
the zero-energy orbits in our problem.
Key words: off-center circular orbits; integrals of motion
2020 Mathematics Subject Classification: 37J35; 68-03
1 Introduction
In Book 1, Proposition 7, Problem 2 of his 1687 Philosophiæ Naturalis Principia Mathematica [9],
Isaac Newton poses the following question:
Let the orbit of a particle moving in a central force field be an off-center circle. We assume
that the orbit encloses the potential center. How does the magnitude of the force depend
on the position on the circle?
He produces the answer that the force magnitude must depend on the position on the circle as
follows:
F (θ) = − 4α
r2(θ) chord3(θ)
, (1.1)
where
r(θ) =
√
R2 + l2 + 2Rl cos θ (1.2)
is the distance between the particle and the force center, and
chord(θ) = 2R(R+ l cos θ)/r(θ)
is the length of the chord of the orbital circle with the particle as one of its ends and the force
center as one of its points; both are given as functions of the particle position on the orbit,
mailto:maxim.olchanyi@umb.edu
https://sites.google.com/view/integrability/home
https://doi.org/10.3842/SIGMA.2023.001
2 M. Olshanii
x
y
r
r′
R
v⃗
θϕ
l
cho
rd
R
R
F⃗
Figure 1. The geometry of the problem. An unknown central force F , acting on a point particle of
mass m, generates a circular orbit with the center at (+l, 0) and a radius R. In Book 1, Proposition 7,
Problem 2 of his 1687 Philosophiæ Naturalis Principia Mathematica [9] Newton shows that on the orbit,
the force must have the form F (θ) = − 4α
r2(θ) chord3(θ)
. Here, (r, ϕ) are the polar coordinates of the particle,
and the function chord(θ) is the length of the chord that has the particle as one of its ends and that
crosses the force center. The distance r′ is the complement of the distance r to the full chord length.
2R is the length of the shortest chord that passes through the origin. The angle θ provides a convenient
parametrization for the points on the orbit.
angle θ (see Figure 1). Here and below, R is the radius of the orbit, l is the distance between
the force center and the center of the orbit, and α is an arbitrary constant.
Further, in Corollary I to [9, Book 1, Proposition 7, Problem 2], Newton looks at the par-
ticular case of orbit passing through the force center (i.e., the case of l = R). In this case,
chord(θ) = r(θ). It leads to the force
F0(θ) = − 4α
r(θ)5
. (1.3)
One may further pose the question of which central potential field can produce the force (1.3)
on the orbit. One obvious surmise is
V0(r) = − α
r4
(1.4)
everywhere.
However, to our knowledge, the problem of finding a potential that induces the force (1.1)
has never been posed. In this article, we find such a potential.
Likewise, the mathematical reason for the appearance of off-center circular orbits has been
identified so far only in the limiting case (1.3)–(1.4). It turns out that the potential (1.4), at
a given energy, is self-dual under the so-called Bohlin–Arnold–Vassiliev (BAV) duality [2, 4, 5,
7, 8], with the dual being at a different energy and featuring a different coupling constant. In
the particular case (1.4), the map between the two is a circular inversion. At zero energy, the
A Novel Potential Featuring Off-Center Circular Orbits 3
potential (1.4) maps to a free problem, at a particular non-zero energy. Circle inversion is known
to map circles and straight lines to circles or straight lines; in the present case, straight-line orbits
of the dual system get mapped to circular orbits of the original problem.
The BAV duality has a long history that starts from a 1742 book by MacLaurin [8, Book II,
Chapter V, Section 875]. His result is described is the review [1, Section 3]. Historically, a map
between a harmonic oscillator and a Coulomb problem has been the principal object of interest.
The timeline of the development is presented in [1, Section 10].
Note that the conventional BAV duality requires a conformal transformation from a Euclidean
space onto itself, with the potential energy of the problem of interest playing the role of a metric
tensor. However, the explanation for the appearance of zero-energy off-center circular orbits
in our case requires a generalization of the BAV scheme to a conformal map from a curved
to an Euclidean space, first introduced by Goursat in 1887 [6]. The map turns out to be
a stereographic projection [11] from a sphere to its equatorial plane. Similarly to the inverse
quartic case, zero-energy orbits of our potential map to a free motion orbits of the dual problem:
those are represented by the grand circles on a sphere. Finally, according to a known property
of the stereographic projection, these circles map back to the circular orbits of the potential we
found.
The circularity of the orbits aside, the zero-energy orbits of (2.1) are also closed curves. We
identify a maximal superintegrability set of conserved quantities responsible for the phenomenon.
In the final section, we summarize our results, put them in a context, and suggest future
projects. In particular, we propose to look at the case where the potential center lies outside of
a circular orbit; such orbits are conjectured to be generated by a free motion on a hyperbolic
surface.
2 Main result: the potential that generates off-center
circular orbits
Theorem 2.1 (a potential description of the force (1.1)). The expression (1.1) for the radial
component of a central force that supports a circular orbit (1.2) admits a potential description
with
V (r) = − α
(r2 +R2)2
, (2.1)
as the potential, where R =
√
R2 − l2 is a half of the shortest chord of the orbit circle that
crosses the force center (see Figure 1).
Proof. According to the intersecting chords theorem, r× r′ = R×R (Figure 1). On the other
hand, the radial component of the force produced by the potential (2.1) is
∂
∂r
V (r) = − 4αr(
r2 +R2
)3 .
Thus, Newton’s expression for the force on the orbit, equation (1.1), can be transformed as
F = − 4α
r2 chord3
= − 4α
r2(r + r′)3
= − 4α
r2
(
r + R2
r
)3 = − 4αr(
r2 +R2
)3 =
∂
∂r
V (r). ■
Remark 2.2 (orbits (1.2) are zero-energy orbits). The circular orbits (1.2) are only a small
subset of the orbits generated by the potential (2.1). In fact, they are the zero-energy orbits for
4 M. Olshanii
this potential. Indeed, consider the aphelion (the rightmost point of the orbit of Figure 1) and
perihelion (the leftmost point there). Energy at these two points will read
E =
L2
z
2m(R+ l)2
− α(
(R+ l)2 +R2
)2 ,
E =
L2
z
2m(R− l)2
− α(
(R− l)2 +R2
)2 ,
respectively, where Lz = mr2ϕ̇ = const is the angular momentum. Solving for E and Lz we get
E = 0,
and Lz = ± 1√
2
√
mα
R2 .
Conversely, every zero-energy orbit of the potential (2.1) has a form (1.2), up to a rotation
about the potential center. This property follows from Remark 3.4.
3 Why zero-energy orbits in a 1/
(
r2 +R2
)2
potential
well are circular
The potential (2.1) has been specifically constructed to support circular orbits at zero energy.
Nonetheless, one would want to know the deeper mathematical reason for this property. In this
section, we will present an example of such a reason.
Let us embed the xy-plane we are working with in a three-dimensional space with coordinates
x, y, and z, and introduce a radius R sphere there with its center at the origin. Consider
a stereographic projection from the sphere to the xy-plane of interest (Figure 2):x
y
0
=
s′ − (ez · s′)ez
1− (ez · s′)
R. (3.1)
Its inverse reads
s′x =
2R2
r2 +R2
x
R
, s′y =
2R2
r2 +R2
y
R
, s′z =
r2 −R2
r2 +R2
. (3.2)
Here, s′ is the unit vector directed to a point on the sphere and ez is the unit vector in the
z-direction.
Theorem 3.1 (correspondence between the motion in a 1/
(
r2 + R2
)2
potential well and
a free motion on a sphere). The stereographic projection (3.1) maps free motion orbits on the
sphere (3.2) to the orbits in the potential (1.4).
Proof. An energy E orbit Ω in the potential (2.1) extremises the Maupertuis–Jacobi functional
S[P] =
∫
P
√
2m(E − V (r)) dl,
where P is any path in the coordinate space that shares the beginning and the end points with
the orbit [3, p. 247]. At zero energy, the Maupertuis–Jacobi action becomes
S[P] =
√
2mE′
∫
P
dl′, (3.3)
A Novel Potential Featuring Off-Center Circular Orbits 5
R
Figure 2. A stereographic projection from a sphere of radius R to its equatorial plane. The spherical
geodesics—the great circles—become, in the plane, circles that cross the sphere’s equator at two dia-
metrically opposite points. In the main text, we show that the free-motion orbits on the sphere become
orbits in the potential (2.1).
with
dl′ =
2R2
r2 +R2
dl and E′ =
α
4R4
. (3.4)
The length element dl′ is the distance traveled on the sphere (3.2) by the stereographic prototype
of the point that traveled a distance dl in the xy-plane. The whole functional (3.3) is then
proportional to the length of the prototype of the path P on the sphere. Thus, a path in xy
that extremises P (i.e., an actual orbit) will be a stereographic image of a free motion path on
the prototype sphere. ■
Remark 3.2. Theorem 3.1 is a particular case of a much more general result that covers any
conformal map between two Riemannian manifolds, at finite as well as at zero energy [6].
Remark 3.3 (why orbits are circles). A stereographic projection from a sphere to its equatorial
plane maps circles (that do not cross the north pole) on the sphere to circles in the equatorial
plane [11]. This explains why the zero energy orbits of (2.1) are circles.
Remark 3.4 (geometry of the orbits, for a given potential). For a given potential (2.1), the
family of zero-energy orbits consists of circles that cross an origin-centered circle of a radius R
at two points diametrically opposite to each other (see Figure 3). Note that this radius R circle
is the intersection between the xy-plane and the stereographic sphere.
Remark 3.5 (relationship between trajectories). Consider a trajectory of a free particle of an
energy E′ given by (3.4), on the sphere (3.2). One can show that r(s′(t′(t))) will be a valid
zero-energy trajectory in the potential (2.1), where the inverse stereographic map r(s′) is given
by (3.2) and t′(t) obeys
dt′
dt
=
4R4(
r(t)2 +R2
)2 . (3.5)
Remark 3.6 (the known limit of a 1/r4 potential well). The well-studied r ≫ R limit leads to
the inverse quartic potential (1.4). There, the inverse stereographic projection (3.2) degenerates
onto a circle inversion, with
R0 =
√
2R
being the radius of the inversion circle, itself centered at the origin.
6 M. Olshanii
n⃗
R
Figure 3. Three representative zero-energy orbits in the potential (2.1), for a given displacement ori-
entation n (solid lines). The origin-centered circle of radius R (dashed-line) is the intersection of the
stereographic sphere with the coordinate plane. The thin dotted line segment is circle’s diameter.
4 Why zero-energy orbits in a 1/
(
r2 +R2
)2
are closed:
an additional integral of motion
Circular or not, the zero-energy orbits of our model are also closed. Such a phenomenon is usually
associated with a presence of 2N − 1 functionally independent integrals of motion where N is
the number of degrees of freedom (the so-called maximal super-integrability set). In our two-
dimensional system, we need to find one additional constant of motion besides the Hamiltonian
H(r,p) =
p2
2m
− α(
r2 +R2
)2
and the angular momentum
Lz(r,p) = xpy − ypx.
Here, p = m d
dtr is the momentum vector.
Theorem 4.1 (an additional integral of motion). At zero energy, a unit vector
n ≡ − sign(Lz)
Ixy
|Ixy|
, (4.1)
with
Ixy ≡
(
Ix(r,p)
Iy(r,p)
)
= − 1
2R
(
rLz(r,p) + (ez × r)Q(r,p) + (ez × p)R2
)
, (4.2)
is conserved, and it is functionally independent from both the Hamiltonian and the angular
momentum. Here,
Q(r,p) = r · p.
Proof. A straightforward calculation shows that
{Ixy, H} = − 2
R
(
−y
+x
)
H;
A Novel Potential Featuring Off-Center Circular Orbits 7
here and below, {A,B} ≡ ∂A
∂r
∂B
∂p − ∂A
∂p
∂B
∂r are the Poisson brackets. Indeed, this commutator
vanishes at zero energy.
For an off-center zero-energy orbit, |Ixy| is strictly greater than zero: indeed, one can show
that the difference between the aphelion (rmax) and perihelion (rmin) distances is
rmax − rmin = 2
√
|Ixy|
|Lz|
R.
Consider such an orbit. If one rotates it about the origin, both the energy H and the angular
momentum Lz will remain unchanged, but the orientation of the vector Ixy changes together
with the orbit. This proves the functional independence of n of H and Lz, only for the off-center
orbits.
In the case of a circular orbit centered at the origin, the functional independence of n is not
needed for a proof of the orbit closeness. ■
Remark 4.2 (magnitude of Ixy). The magnitude of the vector Ixy,
|Ixy| =
√
mα
2R2
− L2
z,
is a function of Lz.
Remark 4.3 (the physical meaning of the integrals of motion (4.2)). Note that each of the
three components of the angular momentum on the sphere,
I ′i = mR2
(
s′ × d
dt′
s′
)
i
,
is conserved. Here, i = 1, 2, 3 correspond to x, y, z, with a cyclic convention 2+2
·
= 1, 3+1
·
= 1,
3 + 2
·
= 2. One can now express these conserved quantities through the variables of the original
problem, using the connections (3.2) and (3.5). It can be shown that the 2D Ix,y(t) = I ′x,y(t
′(t)),
while the 2D angular momentum Lz = I ′z.
Remark 4.4 (the analogy with the Runge–Lenz vector). The vector Ixy, together with the
angular momentum Lz controls the displacement of the orbit center from the force center:
ln = − R
Lz
Ixy.
The vector n will determines the orientation of the position of the orbit center, ln, with respect
to the origin; as such, it is analogous to the Laplace–Runge–Lenz vector of the Coulomb problem.
5 Summary, interpretation of the results, and outlook
This article is inspired by Proposition 7, Problem 2 of Book 1 of Newton’s Principia. There,
Newton gives a formula for how a force should depend on the point on the orbit in order for the
orbit to be an off-center circle. We identify a potentialthat can induce such a force (2.1). All
zero-energy orbits in this potential are, indeed, off-center circles.
Next, we show that the motion in the potential (2.1) is a stereographic image of free motion
on a sphere. This connection can be regarded as a manifestation of the Bohlin–Arnold–Vassiliev
(BAV) duality [2, 4, 5, 7, 8] introduced by MacLaurin, in 1742, and its non-Euclidean gener-
alization by Goursat [6], from 1887. At large distances from the potential center, where our
potential becomes an inverse quartic well, our map degenerates to a circle inversion; the latter
8 M. Olshanii
relates the motion in the inverse quartic potential to a free motion in a flat space, a well-known
manifestation of the BAV duality [2] (a self-duality, in this particular case).
Since the zero-energy orbits in the potential (2.1) are closed curves, we were also compelled
to identify an additional integral of motion, complementary to the Hamiltonian and the angular
momentum. This integral is an analogue of the Runge–Lenz vector of the Coulomb problem.
For future projects, it appears promising to search for other problems that map to a free
motion on curved surfaces with solvable geodesics.
In particular, Newton’s result can be generalized to the case where the potential center lies
outside of the orbit. This motion leads to a potential
Vhyperbolic(r) = − α(
r2 − R̃2
)2 , with R̃ =
√
l2 −R2,
where R and l are the radius of the orbit and the distance between its center and the force
center respectively. We expect that this potential would induce a metric of the Poincaré disk
model of hyperbolic geometry [10, p. 45]. There, the geodesics are circle segments intersecting
the border of the Poincaré disk at a right angle, inside the disk. We would like to point out that
in this case, it can be predicted that the analogue of the integral of motion (4.1) will only be
conserved in between the subsequent collisions with the singularity at r = R̃; angular momentum
conservation is likely to require sudden jumps in n, at the moment of such a collision.
Acknowledgements
We thank Francisco J. Jauffred and Bala Sundaram for inspiring this project and Vanja Dunjko
and Steven Jackson for the helpful discussions. We are indebted to the anonymous referee
who inspired us to drastically reorganize the original draft and provided the absolutely crucial
references [8] and [6]. This work was supported by NSF grant PHY-1912542.
References
[1] Albouy A., Zhao L., Darboux inversions of the Kepler problem, Regul. Chaotic Dyn. 27 (2022), 253–280,
arXiv:2201.00808.
[2] Arnold V., Huygens and Barrow, Newton and Hooke, Birkhäuser, Basel, 1990.
[3] Arnold V.I., Mathematical methods of classical mechanics, 2nd ed., Springer, New York, 1989.
[4] Arnold V.I., Vasilev V.A., Newton’s Principia read 300 years later, Notices Amer. Math. Soc. 36 (1989),
1148–1154, Addendum, Notices Amer. Math. Soc. 37 (1990), 144–144.
[5] Bohlin M., Note sur le problème des deux corps et sur une intégration nouvelle dans le problème des trois
corps, Bull. Astron. 28 (1911), 113–119.
[6] Goursat E., Les transformation isogonales en méchanique, C. R. Math. Acad. Sci. Paris 108 (1887), 446–450.
[7] Horvathy P., Zhang P., The laws of planetary motion, derived from those of a harmonic oscillator (following
Arnold), arXiv:1404.2265.
[8] MacLaurin C., A treatise of fluxions, Ruddimans, Edinburgh, 1742.
[9] Newton I., Mathematical principles of natural philosophy, Daniel Adee, New York, 1846.
[10] Penrose R., The road to reality: a complete guide to the laws of the universe, Alfred A. Knopf, Inc., New
York, 2005.
[11] Rosenfeld B., Sergeeva N., Stereographic projection, Mir, Moscow, 1977.
https://doi.org/10.1134/s1560354722030017
https://arxiv.org/abs/2201.00808
https://doi.org/10.1007/978-3-0348-9129-5
https://doi.org/10.1007/978-1-4757-2063-1
https://arxiv.org/abs/1404.2265
1 Introduction
2 Main result: the potential that generates off-center circular orbits
3 Why zero-energy orbits in a 1/(r^2+R^2)^2 potential well are circular
4 Why zero-energy orbits in a 1/(r^2+R^2)^2 are closed: an additional integral of motion
5 Summary, interpretation of the results, and outlook
References
|
| id | nasplib_isofts_kiev_ua-123456789-211883 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 1815-0659 |
| language | English |
| last_indexed | 2026-03-16T13:28:44Z |
| publishDate | 2023 |
| publisher | Інститут математики НАН України |
| record_format | dspace |
| spelling | Olshanii, Maxim 2026-01-14T09:57:13Z 2023 A Novel Potential Featuring Off-Center Circular Orbits, Maxim Olshanii, SIGMA 19 (2023), 001, 8 pages 1815-0659 2020 Mathematics Subject Classification: 37J35; 68-03 arXiv:2207.09606 https://nasplib.isofts.kiev.ua/handle/123456789/211883 https://doi.org/10.3842/SIGMA.2023.001 In Book 1, Proposition 7, Problem 2 of his 1687 Philosophiae Naturalis Principia Mathematica, Isaac Newton poses and answers the following question: Let the orbit of a particle moving in a central force field be an off-center circle. How does the magnitude of the force depend on the position of the particle onthat circle? In this article, we identify a potential that can produce such a force, only at zero energy. We further map the zero-energy orbits in this potential to finite-energy free motion orbits on a sphere; such a duality is a particular instance of a general result by Goursat, from 1887. The map itself is an inverse stereographic projection, and this fact explains the circularity of the zero-energy orbits in the system of interest. Finally, we identify an additional integral of motion—an analogue of the Runge-Lenz vector in the Coulomb problem—that is responsible for the closeness of the zero-energy orbits in our problem. We thank Francisco J. Jauffred and Bala Sundaram for inspiring this project, and Vanja Dunjko and Steven Jackson for the helpful discussions. We are indebted to the anonymous referee who inspired us to reorganize the original draft drastically and provided the absolutely crucial references [8] and [6]. This work was supported by NSF grant PHY-1912542. en Інститут математики НАН України Symmetry, Integrability and Geometry: Methods and Applications A Novel Potential Featuring Off-Center Circular Orbits Article published earlier |
| spellingShingle | A Novel Potential Featuring Off-Center Circular Orbits Olshanii, Maxim |
| title | A Novel Potential Featuring Off-Center Circular Orbits |
| title_full | A Novel Potential Featuring Off-Center Circular Orbits |
| title_fullStr | A Novel Potential Featuring Off-Center Circular Orbits |
| title_full_unstemmed | A Novel Potential Featuring Off-Center Circular Orbits |
| title_short | A Novel Potential Featuring Off-Center Circular Orbits |
| title_sort | novel potential featuring off-center circular orbits |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/211883 |
| work_keys_str_mv | AT olshaniimaxim anovelpotentialfeaturingoffcentercircularorbits AT olshaniimaxim novelpotentialfeaturingoffcentercircularorbits |