Next Article in Journal
A License Plate Recognition System with Robustness against Adverse Environmental Conditions Using Hopfield’s Neural Network
Next Article in Special Issue
A Swarming Meyer Wavelet Computing Approach to Solve the Transport System of Goods
Previous Article in Journal
Norm of Hilbert Operator’s Commutants
Previous Article in Special Issue
Various Solitons and Other Wave Solutions to the (2+1)-Dimensional Heisenberg Ferromagnetic Spin Chain Dynamical Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monoparametric Families of Orbits Produced by Planar Potentials

Department of Physics, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
Axioms 2023, 12(5), 423; https://doi.org/10.3390/axioms12050423
Submission received: 30 March 2023 / Revised: 23 April 2023 / Accepted: 25 April 2023 / Published: 26 April 2023
(This article belongs to the Special Issue Geometry and Nonlinear Computations in Physics)

Abstract

:
We study the motion of a test particle on the x y plane. The particle trajectories are given by a one-parameter family of orbits f ( x , y ) = c, where c = const. By using the tools of the 2D inverse problem of Newtonian dynamics, we find two-dimensional potentials that produce a pre-assigned monoparametric family of regular orbits f ( x , y ) = c that can be represented by the “slope function γ = f y f x uniquely. We apply a new methodology in order to find potentials depending on specific arguments, i.e., potentials of the form V ( x , y ) = P ( u ) where u = x 2 + y 2 , x y , x 3 y 3 , x y ( x , y 0). Then, we establish one differential condition for the family of orbits f ( x , y ) = c. If it is satisfied, it guarantees the existence of such a potential, generating the above family of planar orbits. Then, the potential function V = V ( x , y ) is found by quadratures. For known families of curves, e.g., ellipse, the logarithmic spiral, the lemniscate of Bernoulli, and circles, we find homogeneous and polynomial potentials that are compatible with this family of orbits. We offer pertinent examples that cover all of the cases, and we examine which of these potentials are integrable. We also study one-dimensional potentials. The families of straight lines in 2D space are also examined.

1. Introduction

The inverse problem of dynamics is an old problem of physics. In [1], Szebehely set the problem of finding all of the potentials V ( x , y ) that can produce a given monoparametric family of curves f ( x , y ) = c traced in the ( x , y )-Cartesian plane by a material point of unit mass. This equation was then studied by several authors (see [2,3,4]). Ref. [5] produced a linear second-order PDE in the unknown function V, which includes only the family of orbits and the potential—not the energy dependence. Two-dimensional homogeneous potentials that generate families of planar orbits were studied by [6], and inhomogeneous potentials were examined by [7] respectively. Family boundary curves were studied by [8], and the allowed region for the motion of the test particle was determined. A review on basic facts of the inverse problem in dynamics was made by [9]. Solvable cases of the planar inverse problem were given in [10]. In a series of papers ([11,12,13]), the authors presented methodologies; with the aid of these methodologies, one can obtain monoparametric families of planar orbits in the direct problem of Newtonian dynamics. Families of straight lines for planar potentials were also studied by [14]. Ref. [15] studied the monoparametric isoenergetic families of planar orbits f ( x , y ) = c created by homogeneous potentials V = V ( x , y ) . Ref. [16] studied the two basic equations of the inverse problem in dynamics by using another methodology. Moreover, Ref. [17] examined a solvable version of the planar inverse problem and found potentials of special type V = V ( γ ) , where γ = f y f x .
In the present work, we address the following question: Given a monoparametric family of regular curves f ( x , y ) = c, is there any potential V = P ( u ( x , y ) ) that generates this family of curves as trajectories? Many results were found in the past for the planar inverse problem; see, e.g., [9]. In order to give an appropriate answer to the above question, we select
(i)
u = x 2 + y 2 in order to find central potentials that have applications in Celestial Mechanics, e.g., the Newtonian one, and others,
(ii)
u = x y in order to determine 2D potentials that produce specific families of curves as orbits, e.g., the lemniscate of Bernoulli,
(iii)
u = x 3 + k y 3 for finding cubic potentials,
(iv)
u = x y in order to find homogeneous potentials of zero-degree and other results.
The paper is organized as follows. In Section 2, we present the basic facts of the 2D inverse problem of dynamics. In Section 3, we develop a new methodology for finding potentials of the special form V = P ( u ) , where P is an arbitrary function of C 2 class (u is defined above). Proceeding further, we find one compatibility condition for the orbital function γ = γ ( x , y ) that is related to the given family of orbits. If this condition is fulfilled, then such a potential exists, and it is found by quadratures. By using this methodology, someone can find two-dimensional or one-dimensional potentials that produce the given family of orbits. Indeed, in Section 4 we find central potentials of the form V = P ( x 2 + y 2 ) compatible with families of orbits, and we focus our interest to the Newtonian, cored, and logarithmic ones. Especially, the cored and logarithmic potentials have an absolute minimum and reflection symmetry with respect to both axes. In addition to that, they are very interesting in terms of problems of galactic dynamics and in terms of the models that are used in the study of elliptical galaxies. In Section 5, we study the potentials of the form V = P ( x y ) , and we find a family of curves, e.g., lemniscate of Bernoulli, produced by potentials as orbits. In Section 6, we study cubic potentials, and other results are given in Section 7. Furthermore, we find families of orbits compatible with integrable potentials (Section 8). We give figures of families of orbits in many cases. One-dimensional potentials are examined in Section 9. We refer to the case of straight lines (Section 10), and we make some concluding comments in Section 11.

2. The Basic Equation of the 2D Inverse Problem

We consider the monoparametric family of planar orbits
f ( x , y ) = c , c = c o n s t .
that is traced by a material point of unit mass under the action of the potential V = V ( x , y ) . The total energy is constant. In a two-dimensional frame, we consider monoparametric families of orbits given in the form (1). As indicated by ([5,6]), the family of orbits (1) can be represented by the “slope function γ = γ ( x , y )
γ ( x , y ) = f y f x ,
and
Γ = γ γ x γ y .
The following notation is useful:
κ = 1 γ γ , λ = 1 Γ Γ x + 1 γ Γ y , μ = λ γ + 3 Γ γ .
We note here the subindices “x, y” denote partial derivatives. There is a “one-to-one” correspondence between the slope function (2) and the family of orbits (1). This means that if the slope function γ is given in advance, then we can find the monoparametric family of orbits in the form (1) by analytically solving the ODE
d y d x = 1 γ ( x , y ) .
The potential V = V ( x , y ) has to satisfy a linear PDE of second order ([5,6])). Taking into account that Γ 0, this PDE reads
V x x + κ V x y + V y y = λ V x + μ V y .
Equation (6) can be written in polar coordinates. In doing so, [18] studied geometrically similar orbits and applied Yoshida’s criterion of non-integrability for some homogeneous potentials. The energy of family of orbits (1) is found to be ([11], p. 248):
E = V ( 1 + γ 2 ) ( V x + γ V y ) 2 Γ .
From (7), we ascertain that E V 0. As it was shown by [8], this requirement leads to the fact
A = V x + γ V y Γ 0 .
The inequality (8) defines the allowed region of motion of the test particle in the orbits of the family (1) in the x y plane, which are traced by a particle of unit mass in the presence of a 2D potential V = V ( x , y ) .
An interesting case of the family of curves on the x y plane is the straight lines. As it was shown by ([14], p. 4), the curvature of the orbits in (1) is given by
K = Γ 1 + γ 3 / 2 .
Consequently, the family of orbits (1) consists of straight lines if and only if Γ = 0. In view of (3), this condition can be written as
γ γ x γ y = 0 ,
or, equivalently, γ = γ x γ y . This case will be studied in Section 9.

3. The Methodology for the General Case

In this section, we find solutions of the form
V = P ( u ) ,
for the second order PDE (6). Here, u is an mathematical expression of the variables x , y , and P is an arbitrary function of C 2 −class. The case u = γ was studied extensively by [17] and will not be considered here.

One Condition on the Slope Function γ

We suppose that V = P ( u ) , and we find the derivatives of first and second order of the potential function V with respect to x , y , respectively. It is:
V x = P u x , V y = P u y ,
where P ( u ) = d P d u and
V x x = P u x 2 + P u x x , V x y = P u x u y + P u x y , V y y = P u y 2 + P u y y ,
and we insert them into Equation (6). Thus, we obtain the next relation
q 2 P ( u ) = q 1 P ( u ) ,
where
q 2 = u x 2 + κ u x u y + u y 2 , q 1 = λ u x + μ u y + u x x κ u x y u y y ,
If q 2 0, then we obtain from (13)
P P = q 1 q 2 = Q .
Now, we observe that P P depends only on the argument u. Consequently, the function Q must depend on the same argument. Then, we obtain
Q = Q ( u )
or, equivalently,
Q x u x = Q y u y .
This relation (17) is verified if and only if
Q x u y Q y u x = 0 .
The differential condition (18) is the differential condition for the slope function γ , which, if it is satisfied, ensures the existence of a potential (10). On the other hand, we consider that the condition (18) is satisfied by the slope function γ . Besides that, we have P ( u ) 0 . Then, we obtain P = P ( u ) , and we can determine the function P from (15). The result is
P ( u ) = d 1 e Q ( u ) d u d u + d 2 , d 1 , d 2 = c o n s t . 0 .
As a conclusion, by using the result (19), we can find the potential function V by quadratures. Now, we can formulate the next result.
Proposition 1.
If we consider that q 2 0 and q 1 0 and the differential condition (18) is satisfied for the given family of orbits (1), then a potential of the form V = P ( u ) 0 always exists, and it is determined uniquely from (19) up to two arbitrary constants.
Synthesis of the problem
(1)
We start with the family of orbits (1), and we determine the functions γ and Γ from (2) and (3), respectively.
(2)
If Γ = 0, then the family of orbits consists of straight lines, and the potential is found from the relation V x + γ V y = 0 ([14], p. 4).
(3)
If Γ 0, then we estimate the quantities q 1 , q 2 from (14), and we check the condition (18).
(4)
If the condition (18) is not satisfied, then we have no solution to the problem. On the other hand, if the condition (18) is satisfied, then we determine the potential of the form V = P ( u ) from (19).
We remark here that
  • Case 1. If q 2 0 and q 1 = 0, then we obtain P ( u ) = 0 from (13). Then, we have P ( u ) = d 1 u + d 2 , d 1 , d 2 = c o n s t .
  • Case 2. If q 2 = 0 and q 1 = 0, then the Equation (13) is identically zero. Thus, any function V = P ( u ) is the solution to our problem.
  • Case 3. If q 2 = 0 and q 1 0, we obtain P ( u ) = 0 from (13). Thus, P ( u ) = c o n s t . and this case is excluded from our study.

4. Central Potentials P = P ( x 2 + y 2 )

In this Section, we shall offer examples that cover the above cases. We select u = x 2 + y 2 , and we give the first
Example 1.
We consider the monoparametric family of orbits (ellipses, see Figure 1a)
f ( x , y ) = x + 2 x 2 + y 2 = c ,
which leads to the slope function
γ = 2 y 2 x + x 2 + y 2 , Γ 0 ,
Here, we have q 2 , q 1 0, and we check the condition (18). It is satisfied, and we determine the function Q ( u ) from (15). It is:
Q ( u ) = 3 2 u , u = x 2 + y 2 ,
and, from (19), we find
P ( u ) = 2 d 1 u + d 2 , d 1 , d 2 = c o n s t .
If we set d 1 = 1 2 and d 2 = 0, we conclude that the family of orbits (20) is created by the Newtonian potential
V ( x , y ) = 1 x 2 + y 2 .
The potential in (24) is non-integrable ([19]). The polynomial integrability of two-dimensional Hamiltonian systems with homogeneous potentials of degree k was also studied by [20]. The energy of the family of orbits (20) is found from (7) to be
E = 3 4 c = c o n s t .
and the allowed region is found from (8) to be A = x + 2 x 2 + y 2 0.
Example 2.
We regard the curves (logarithmic spirals, see Figure 1b) in the 2D Euclidean plane
x = a e b t cos t , y = a e b t sin t , a , b 0 .
Eliminating the parameter t from the Equation (26), we obtain the monoparametric family of orbits in x , y -coordinates
f ( x , y ) = ( x 2 + y 2 ) exp 2 b arctan ( y x ) = a 2 = c ,
which is represented by the slope function
γ = b x + y x + b y , Γ 0 ,
Here, it is q 2 , q 1 0, and we check the condition (18). It is satisfied, and we determine the function Q ( u ) from (15). It is:
Q ( u ) = 2 u , u = x 2 + y 2 ,
and, from (19), we find
P ( u ) = d 1 u + d 2 , d 1 , d 2 = c o n s t .
If we set d 1 = 1 and d 2 = 0, we conclude that the family of orbits (27) is produced by the potential (homogeneous of degree m = 2)
V ( x , y ) = 1 x 2 + y 2 .
The potential in (31) is completely integrable as it was shown by ([20], p. 3878). The energy is found to be (7)
E = 0 .
and the allowed region is the entire plane. We note here that the family of orbits (27) is trcaced isoenergetically ( E = 0) by a material point of unit mass under the action of potential (31). Such families of planar orbits were studied in detail by [15].

Special Cases

In this paragraph, we shall present one example that corresponds to the Case 2 of the general theory.
Example 3.
We take into account the monoparametric family of orbits (circles)
f ( x , y ) = x 2 + y 2 = c ,
which leads to the slope function
γ = y x , Γ = x 2 + y 2 x 3 0 .
Here, we have q 2 = q 1 = 0 (Case 2). Thus, any function V = P ( u ) is a solution to our problem. Thus, we have
P ( u ) = R ( u ) ,
where R is an arbitrary function of C 2 −class. The corresponding potential is
V ( x , y ) = R ( x 2 + y 2 ) .
It is remarkable to say that
a
If we select R ( u ) = 1 + u , then we obtain the cored potential V = 1 + x 2 + y 2 , and the energy of the family of orbits is
E = 2 + 3 c 2 1 + c
and the allowed region is the entire plane ( A = x 2 1 + x 2 + y 2 0).
b
If we select R ( u ) = log ( 1 + u ) , then we obtain the loagarithmic potential V = log ( 1 + x 2 + y 2 ) and the energy of the family of orbits is
E = c 1 + c + log ( 1 + c )
and the allowed region is the entire plane ( A = 2 x 2 1 + x 2 + y 2 0).
Both of them were studied in [21].

5. Potentials of the Form P = P ( x y )

At the beginning, we shall prove the following
Theorem 1.
All of the monoparametric families of orbits f ( x , y ) = M ( x 2 + y 2 ) = c where M is an arbitrary function of C 2 class are produced by the potential V ( x , y ) = 1 x y and the allowed region for a test particle is x y > 0 ( A = 2 x y ( x 2 + y 2 ) 0).
Proof. 
We consider the family of orbits
f ( x , y ) = M ( x 2 + y 2 ) = c ,
where M is an arbitrary function of C 2 class and the slope function is
γ = y x , Γ 0 .
Here, it is q 2 , q 1 0, and we check the condition (18). It is satisfied, and we determine the function Q ( u ) from (15). It is:
Q ( u ) = 2 u , u = x y ,
and, from (19), we find
P ( u ) = d 1 u + d 2 , d 1 , d 2 = c o n s t .
For d 1 = 1 , d 2 = 0, the potential function is
V ( x , y ) = 1 x y .
which is non-integrable in Liouville sense ([22], p. 80). Proceeding more, we present the following. □
Example 4.
We consider the curves (lemniscates of Bernoulli) on the 2D plane (see Figure 2a)
x 2 + y 2 = 2 a 2 ( x 2 y 2 ) , a 0 .
This family of curves can be written as follows:
f ( x , y ) = x 2 + y 2 x 2 y 2 = 2 a 2 = c ,
which is related to the slope function
γ = 3 x 2 y y 3 x 3 3 x y 2 , Γ 0 .
Here, it is q 2 , q 1 0, and we check the condition (18). It is satisfied, and we determine the function Q ( u ) from (15). It is:
Q ( u ) = 5 2 u , u = x y ,
and, from (19), we find
P ( u ) = 2 d 1 3 u 3 / 2 + d 2 , d 1 , d 2 = c o n s t .
For d 1 = 1 and d 2 = 0, the potential function is
V ( x , y ) = 2 3 ( x y ) 3 / 2 ,
The potential in (47) is integrable in the Liouville sense ([22], p. 80), and the energy of the family of orbits (43) is found to be E = 0. The allowed region for a test particle is x y > 0 ( A = 4 ( x 3 3 x y 2 ) 2 3 ( x y ) 3 / 2 ( x 2 + y 2 ) 3 0).
Other results are presented at Table 1.

Special Cases

In this paragraph, we shall present two examples that correspond to the Cases 1 and 3 of the general theory.
Example 5.
As we explained in Section 2, a monoparametric family of orbits f ( x , y ) = c can be represented by the slope function γ. Thus, here we shall consider a family of orbits studied by ([6], p. 470) that corresponds to
γ = r + s ,
where
r = 2 x y 1 x 2 + y 2 , s = r 2 1
Here, we have q 2 0 , q 1 = 0 (Case 1), and we check the condition (18). It is satisfied, and we find the function Q ( u ) from (15). It is: Q ( u ) = 0, and the function P = P ( u ) is found to be
P ( u ) = d 1 u + d 2 , d 1 , d 2 = c o n s t .
If we set d 1 = 1 and d 2 = 0, we find the two-dimensional potential
V ( x , y ) = x y .
The potential in (51) is non-integrable ([20,22]), and the energy of this family of orbits is E = 1 2 .
Example 6.(Counterexample)
We take into account the monoparametric family of orbits
f ( x , y ) = x y ( 1 + x y ) = c ,
which leads to the slope function
γ = x y , Γ = 2 x y 2 0 ,
Here, we have q 2 = 0 , q 1 0 (Case 3). Thus, P ( u ) = c o n s t . , which is excluded from our study. So, there is no solution to our problem.

6. Cubic Potentials

In this Section, we shall deal with polynomial potentials of degree m = 3. We select u = x 3 + k y 3 ( k Z , k 0 ), and we start with
Example 7.
We shall examine the family of orbits (see Figure 3a)
f ( x , y ) = x 3 + x 2 y = c , c = c o n s t .
which leads to the slope function
γ = x 3 x + 2 y , Γ 0 ,
and we shall search the potentials of the form
V ( x , y ) = P ( u ) , u = x 3 + k y 3 , k = c o n s t .
We check the condition (18). It is satisfied only for k = 1, and this leads to the case q 2 0 , q 1 = 0 (Case 1). Then, we determine the function Q ( u ) . It is as follows: Q ( u ) = 0 and the potential function is given by
P ( u ) = d 1 u + d 2 , d 1 , d 2 = c o n s t .
For d 1 = 1 and d 2 = 0, we obtain
V ( x , y ) = x 3 y 3 .
This result agrees with the result found by ([6], p. 471) by using another methodology. We note here that the potential in (58) is separable in Cartesian coordinates and, consequently, is integrable ([19], pp. 8614, 8619). The energy of the family of orbits is calculated from (7)
E = 13 2 c = c o n s t .
and the allowed region for a test particle is 3 x y 0 ( A = 1 2 ( 3 x y ) ( 3 x + 2 y ) 2 0).
Some interesting results are shown at Table 2 ( d 1 , d 2 > 0 ).

7. Other Results

Potentials of the Form V = P ( x y )

We shall prove the following
Theorem 2.
All the monoparametric families of orbits f ( x , y ) = f ( x , y ) = x y m = c ( m 0 , 1 ) are produced by the potential V ( x , y ) = y 2 2 m ( m x 2 y 2 ) .
Proof. 
We regard the family of orbits
f ( x , y ) = x y m = c , m 0 , 1 ,
which leads to the slope function
γ = m x y , Γ = m ( m + 1 ) x y ,
and we shall look for potentials of the form
V ( x , y ) = P ( u ) , u = x y n , n m .
We check the condition (18). It is satisfied only for n = 1. Thus, we take u = x y , and we are interested in potentials
V ( x , y ) = P ( x y ) .
The fact n = 1 leads to the case q 2 , q 1 0 . Then, we determine the function Q ( u ) . It is:
Q ( u ) = 1 + 3 m u 2 u m u 3 ,
and the potential function is given by
P ( u ) = d 1 2 m ( 1 m u 2 ) + d 2 , d 1 , d 2 = c o n s t .
For d 1 = 1 and d 2 = 0, we obtain
V ( x , y ) = y 2 2 m ( m x 2 y 2 ) ,
which is homogeneous of zero degree. The energy of the family of orbits is
E = 1 2 ( 1 + m ) = c o n s t .
and the allowed region for a test particle is determined from (8). It is:
A = y 2 m ( m + 1 ) ( m x 2 y 2 ) 0 .
From (68), we conclude that the area for the motion of test particle is A = m x 2 y 2 > 0 for m > 0, or, for m < 0. In the next step, we offer the following. □
Example 8.
We consider the monoparametric family of orbits (astroids curves, see Figure 4a)
f ( x , y ) = x 2 / 3 + y 2 / 3 = c ,
which leads to the slope function
γ = x 1 / 3 y 1 / 3 , Γ 0 , u = γ .
From (15) we obtain q 2 , q 1 0, and we check the condition (18). It is satisfied, and we determine the function Q ( u ) from (15). It is:
Q ( u ) = 1 + 2 u 2 + 3 u 4 u ( u 4 1 ) ,
and, from (19), we find
P ( u ) = d 1 1 u 4 2 ( 1 u 2 ) + d 2 , d 1 , d 2 = c o n s t .
and the allowed region is A = y 2 / 3 x 2 / 3 > 0.

8. Integrable Potentials

A Hamiltonian system with n degrees of freedom is integrable if the system admits n independent first integrals in involution (Liouville integrability). Let us consider a two-degree of freedom Hamiltonian
H = 1 2 ( p x 2 + p y 2 ) + V ( x , y ) .
We assume that Φ is a first integral of (73). Then, the Poisson bracket of Φ and H vanishes, i.e.,
[ Φ , H ] = 0 .
The question of integrability of 2D potentials was addressed by many authors in the past. More precisely, Ref. [23] studied a family of dynamical systems related to the motion of a test particle in two-dimensional space and found integrals of motion quadratic in velocities. The same authors extended their work to third- and fourth degree-polynomial potentials ([24]). In both cases, they used the weak-Pianlevé property as a criterion for integrability. Polynomial integrals of motion of a degree greater than 2 for planar systems have been found in the review paper of [25]. Integrable velocity-dependent potentials that have linear or quadratric invariants in the velocity have been studied by [26]; integrable velocity-dependent potentials with logarithmic integrals of motion were found by [27]. In [28], the author derived a new equation equivalent to the equation that follows from the condition of vanishing Poisson bracket for autonomous systems of two degrees of freedom. In addition to that, Ref. [29] formulated four necessary and sufficient conditions for a potential function U ( x , y ) in order to be integrable with the second integral of motion quartic in velocities. In the meanwhile, Ref. [30] combined the results of Darboux’s theory with the inverse problem of dynamics and found monoparametric families of orbits sufficient for the integrability of planar potentials having linear or quadratic invariants. Proceeding further, Ref. [31] found planar potentials with a second integral quadratic in the momenta by using Darboux’s integrability criterion. Ref. [19] gave a complete list of all of the integrable two-dimensional homogeneous polynomial potentials with a polynomial integral of order of at most four in the momenta. For algebraic homogeneous potentials of a non-zero rational homogeneity degree, necessary integrability conditions were constructed by [22]. The polynomial integrability of Hamiltonian systems with homogeneous potentials of degree k was examined by [20].
As an example, we shall study an integrable potential (see [19], p. 8614), a polynomial of the third degree, and we shall find a family of planar orbits compatible with it. This is
V ( x , y ) = x 3 + 3 16 x y 2 .
Example 9.
We shall consider the family of orbits
f ( x , y ) = f ( x , y ) = f ( x , y ) = x y m = c , m 0 , 1 ,
which leads to the slope function
γ = m x y , Γ 0 ,
and we shall search potentials of the form
V ( x , y ) = P ( u ) , u = x 3 + 3 16 x y 2 .
We check the condition (18). It is satisfied only for m = 4, and this leads to the case q 2 0 , q 1 = 0 (Case 1). Then, we determine the function Q ( u ) . It is: Q ( u ) = 0, and the potential function is given by
P ( u ) = d 1 u + d 2 , d 1 , d 2 = c o n s t .
For d 1 = 1 and d 2 = 0, we find the potential
V ( x , y ) = x 3 + 3 16 x y 2 .
which produces the family of orbits
f ( x , y ) = x y 4 = c .
and the allowed region is x < 0 ( A = 8 x 2 y 2 + y 4 64 x 0).

9. One-Dimensional Potentials

In this Section, we shall examine potentials of the form V = V ( y ) .
Example 10.
We consider the monoparametric family of orbits (parabolas)
f ( x , y ) = x + y 2 = c ,
which is represented by the slope function
γ = 2 y , Γ = 2 0 ,
and we shall find potentials of the form
V ( x , y ) = P ( u ) , u = γ .
From (14), we find the quantities q 1 , q 2 . These are: q 1 = 6 y , q 2 = 4 , and Q ( y ) = 3 2 y . The condition (18) is satisfied, and thus we have Q ( u ) = 3 u . From (19), we determine the potential function P = P ( u ) .
P = d 1 2 u 2 + d 2 , d 1 , d 2 0 .
If we set d 1 = 1 and d 2 = 0, then we obtain
V ( x , y ) = 1 8 y 2 .
The potential (85) is integrable, and it was studied in detail by ([31], p. 1635). The energy of the family of orbits (81) is E = 1 2 and the allowed region is the entire plane ( A = 1 8 y 2 0).

10. Families of Straight Lines

If Γ = 0, then we have to study a one-parameter family of straight lines (FSL) in 2-D space. As it was shown by [14], potentials that generate one-parametric families of straight lines on the x y plane have to satisfy the following necessary and sufficient differential condition:
V x V y ( V x x V y y ) = V x y ( V x 2 V y 2 ) .
We examined the following potentials:
I.
V = P ( x 2 + y 2 ) ,
II.
V = P ( c 1 x + c 2 y ) , c 1 , c 2 = c o n s t .
III.
V = P ( x y n ) , n 0
IV.
V = P ( x 3 + k y 3 ) , k 0
We replaced them into (86), and we found that only the first two cases, i.e., I and II, are applicable for the study of straight lines. Then, we find the corresponding family of staight lines
γ = V x V y .
For the first case, we obtain
γ = x y
and the family of straight lines is
f ( x , y ) = y x = c .
For the second case, we obtain
γ = c 1 c 2 = c 3 ,
or, equivalently,
f ( x , y ) = y = c 3 x + c 4 , c 3 , c 4 = c o n s t .

11. Conclusions

In the present paper, we studied four solvable versions of the 2D inverse problem of Newtonian dynamics, and we gave the reader a new idea on how to find new solutions to it by using the basic equations. This problem relates two-dimensional potentials V = V ( x , y ) with a preassigned monoparametric family of regular orbits f ( x , y ) = c = const.
We used the basic PDE (6) of the inverse problem of dynamics, taking into account that the quantity Γ is different from zero (Section 2). We imposed one differential condition on the slope function γ , i.e., Equation (18), in order to have solutions for our problem. We studied potentials of the form V ( x , y ) = P ( u ) where u = x 2 + y 2 , x y , x 3 y 3 , x y ( x , y 0). We focused our interest on central potentials that have many applications in physical problems. In order to obtain interesting results, we studied known curves in the x y plane, e.g., circles, ellipses, logarithmic spirals, lemniscates of Bernoulli, etc., which can be traced by a test particle of unit mass as orbits. Our results were not restricted only to central potentials, but we extended them to homogeneous potentials of degree m and to other known potentials from the bibliography, i.e., cubic potentials. We did not obtain only mathematical results; we also found potentials with applications in many areas of physics, e.g., galactic dynamics. Such potentials are the following ones: the Newtonian, the cored, and the logarithmic ones. Our aim was to find a suitable pair of orbits that are produced by these potentials. Furthermore, we found monoparametric families of orbits produced by integrable potentials. New figures with families of orbits were presented in all cases. All of the results are really new and original. We also studied the case of straight lines, which is a special category of orbits in 2D space.

Funding

The author, Thomas Kotoulas, received no funds for this work.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

I would like to thank three anonymous reviewers for their helpful comments, which improved the first version of this manuscript.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Szebehely, V. On the determination of the potential by satellite observations. In Proceedings of the International Meeting on Earth’s Rotation by Satellite Observation; Proverbio, G., Ed.; The University of Cagliari: Bologna, Italy, 1974; pp. 31–35. [Google Scholar]
  2. Bozis, G. Generalization of Szebehely’s Equation. Cel. Mech. 1983, 29, 329–334. [Google Scholar] [CrossRef]
  3. Bozis, G.; Tsarouhas, G. Conservative fields derived from two monoparametric families of planar orbits. Astron. Astrophys. 1985, 145, 215–220. [Google Scholar]
  4. Puel, F. Intrinsic formulation of the equation of Szebehely. Cel. Mech. 1984, 32, 209–216. [Google Scholar] [CrossRef]
  5. Bozis, G. Szebehely’s inverse problem for finite symmetrical material concentrations. Astron. Astrophys. 1984, 134, 360–364. [Google Scholar]
  6. Bozis, G.; Grigoriadou, S. Families of planar orbits generated by homogeneous potentials. Cel. Mech. Dyn. Astr. 1993, 57, 461–472. [Google Scholar] [CrossRef]
  7. Bozis, G.; Anisiu, M.-C.; Blaga, C. Inhomogeneous potentials producing homogeneous orbits. Astron. Nach. 1997, 318, 313–318. [Google Scholar] [CrossRef]
  8. Bozis, G.; Ichtiaroglou, S. Boundary Curves for Families of Planar Orbits. Cel. Mech. Dyn. Astr. 1993, 58, 371–385. [Google Scholar] [CrossRef]
  9. Bozis, G. The inverse problem of dynamics: Basic facts. Inverse Probl. 1995, 11, 687–708. [Google Scholar] [CrossRef]
  10. Grigoriadou, S.; Bozis, G.; Elmabsout, B. Solvable cases of Szebehely’s equation. Cel. Mech. Dyn. Astr. 1999, 74, 211–221. [Google Scholar] [CrossRef]
  11. Anisiu, M.-C.; Bozis, G.; Blaga, C. Special families of orbits in the direct problem of dynamics. Cel. Mech. Dyn. Astr. 2004, 88, 245–257. [Google Scholar] [CrossRef]
  12. Blaga, C.; Anisiu, M.-C.; Bozis, G. New solutions in the direct problem of dynamics. PADEU 2006, 17, 13. [Google Scholar]
  13. Bozis, G.; Anisiu, M.-C.; Blaga, C. A solavable version of the direct problem of dynamics. Rom. Astron. J. 2000, 10, 59–70. [Google Scholar]
  14. Bozis, G.; Anisiu, M.-C. Families of straight lines in planar potentials. Rom. Astron. J. 2001, 11, 27–43. [Google Scholar]
  15. Borghero, F.; Bozis, G. Isoenergetic families of planar orbits generated by homogeneous potentials. Meccanica 2002, 37, 545–554. [Google Scholar] [CrossRef]
  16. Anisiu, M.-C. An alternative point of view on the equations of the inverse problem of dynamics. Inverse Probl. 2004, 20, 1865–1872. [Google Scholar] [CrossRef]
  17. Bozis, G.; Anisiu, M.-C. A solvable version of the inverse problem of dynamics. Inverse Probl. 2005, 21, 487–497. [Google Scholar] [CrossRef]
  18. Bozis, G.; Meletlidou, E. Nonintegrability Detected from Geometrically Similar Orbits. Cel. Mech. Dyn. Astr. 1997, 68, 335–346. [Google Scholar] [CrossRef]
  19. Nakagawa, K.; Yoshida, H. A list of all integrable two-dimensional homogeneous polynomial potentials with a polynomial integral of order at most four in the momenta. J. Phys. A Math. Gen 2001, 34, 8611–8630. [Google Scholar] [CrossRef] [Green Version]
  20. Oliveira, R.; Valls, C. Polynomial integrability of Hamiltonian systems with homogeneous potentials of degree -k. Phys. Let. A 2016, 380, 3876–3880. [Google Scholar] [CrossRef]
  21. Jimenez-Lara, L.; Llibre, J. The cored and logarithmic potentials: Periodic orbits and integrability. J. Math. Phys. 2012, 53, 1–12. [Google Scholar] [CrossRef] [Green Version]
  22. Maciejewski, A.; Przybylska, M. Integrability of Hamiltonian systems with algebraic potentials. Phys. Let. A 2016, 380, 76–82. [Google Scholar] [CrossRef]
  23. Dorizzi, B.; Grammaticos, B.; Ramani, A. A new class of integrable systems. J. Math. Phys. 1983, 24, 2282–2288. [Google Scholar] [CrossRef]
  24. Dorizzi, B.; Grammaticos, B.; Ramani, A. Integrability of Hamiltonians with third- and fourth-degree polynomial potentials. J. Math. Phys. 1983, 24, 2288–2295. [Google Scholar]
  25. Hietarinta, J. Direct methods for the search of second invariants. Phys. Rep. 1987, 147, 87–154. [Google Scholar] [CrossRef]
  26. Dorizzi, B.; Grammaticos, B.; Ramani, A. Integrable Hamiltonian systems with velocity-dependent potentials. J. Math. Phys. 1985, 26, 3070–3079. [Google Scholar] [CrossRef]
  27. Icthiaroglou, S.; Voyatzis, G. Integrable potentials with logarithmic integrals of motion. J. Phys. A Math. Gen. 1988, 21, 3537–3546. [Google Scholar] [CrossRef]
  28. Bozis, G. A transformed equation for a vanishing Poisson bracket. J. Phys. A Math. Gen. 1989, 22, 1759–1764. [Google Scholar] [CrossRef]
  29. Bozis, G. Two-dimensional integrable potentials with quartic invariants. J. Phys. A Math. Gen. 1992, 25, 3329–3351. [Google Scholar] [CrossRef]
  30. Icthiaroglou, S.; Meletlidou, E. On monoparametric families of orbits sufficient for integrability of planar potentials with linear or quadratic invariants. J. Phys. A Math. Gen. 1990, 23, 3673–3679. [Google Scholar] [CrossRef]
  31. Grigoriadou, S. The inverse problem of dynamics and Darboux’s integrability criterion. Inverse Probl. 1999, 15, 1621–1637. [Google Scholar] [CrossRef]
Figure 1. (a) Some members of the family of orbits (20) (ellipses) for c = 3, 5, 7, 9. (b) Members of the family of orbits (26) for a = 1 and b = 0.25, 0.4, and 0.55 successively.
Figure 1. (a) Some members of the family of orbits (20) (ellipses) for c = 3, 5, 7, 9. (b) Members of the family of orbits (26) for a = 1 and b = 0.25, 0.4, and 0.55 successively.
Axioms 12 00423 g001
Figure 2. (a) Some members of the family of orbits (43) for c = 1, 2, 3, and 4. (b) A 3D plot of the potential (47), which produces the family of orbits (42).
Figure 2. (a) Some members of the family of orbits (43) for c = 1, 2, 3, and 4. (b) A 3D plot of the potential (47), which produces the family of orbits (42).
Axioms 12 00423 g002
Figure 3. (a) Some members of the family of orbits (54) for c = 1, 5, and 9. (b) Some members of the family of orbits (60) for c = 1, 5, 9, and m = 2.
Figure 3. (a) Some members of the family of orbits (54) for c = 1, 5, and 9. (b) Some members of the family of orbits (60) for c = 1, 5, 9, and m = 2.
Axioms 12 00423 g003
Figure 4. (a) Some members of the family of orbits (69) for c = 0.04, 0.36, and 1. (b) Some members of the family of orbits (80) for c = 0.25, 0.5, and 1.
Figure 4. (a) Some members of the family of orbits (69) for c = 0.04, 0.36, and 1. (b) Some members of the family of orbits (80) for c = 0.25, 0.5, and 1.
Axioms 12 00423 g004
Table 1. Families of orbits and potentials.
Table 1. Families of orbits and potentials.
Family of OrbitsPotential V(x,y)EnergyAllowed Region
f ( x , y ) = x 4 / 3 + y 4 / 3 = c 3 2 ( x y ) 2 / 3 3 c y > 0
f ( x , y ) = x y ( x + y ) = c d 1 6 x 6 y 6 ( d 1 < 0) 1 12 d 1 c 4 the entire plane
The first potential is integrable in Liouville sense ([22], p. 81), but the second one is non-integrable.
Table 2. Families of orbits and potentials.
Table 2. Families of orbits and potentials.
Family of OrbitsPotential V(x,y)EnergyAllowed Region
f ( x , y ) = x ( x + y ) 2 = c d 2 y ( 3 x 2 + y 2 ) 13 d 2 2 c x + 2 y < 0
f ( x , y ) = y ( x + y ) 2 = c d 1 x ( x 2 + 3 y 2 ) 13 d 1 2 c 2 x + y < 0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kotoulas, T. Monoparametric Families of Orbits Produced by Planar Potentials. Axioms 2023, 12, 423. https://doi.org/10.3390/axioms12050423

AMA Style

Kotoulas T. Monoparametric Families of Orbits Produced by Planar Potentials. Axioms. 2023; 12(5):423. https://doi.org/10.3390/axioms12050423

Chicago/Turabian Style

Kotoulas, Thomas. 2023. "Monoparametric Families of Orbits Produced by Planar Potentials" Axioms 12, no. 5: 423. https://doi.org/10.3390/axioms12050423

APA Style

Kotoulas, T. (2023). Monoparametric Families of Orbits Produced by Planar Potentials. Axioms, 12(5), 423. https://doi.org/10.3390/axioms12050423

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop