Two‐dimension potentials which generate spatial families of orbits

Abstract

We consider the following case of the 3D inverse problem of dynamics: Given a spatial two-parametric family of curves \(f(x,y,z)=c_{1},g(x,y,z)=c_{2}\), find possibly existing two-dimension potentials under whose action the curves of the family are trajectories for a unit mass particle. First we establish the
conditions which must be fulfilled by the family so that potentials of the form \(w(y,z)\) give rise to the curves of the family, and we present some applications. Then we examine briefly the existence of potentials depending on \((x,z)\), respectively \((x,y)\), which are compatible with the given family.

Authors

Mira Cristiana Anisiu
T. Popoviciu Institute of Numerical Analysis,  Cluj-Napoca, Romania

George Bozis
Aristotle University of Thesaloniki, GR-54006, Greece

Keywords

celestial mechanics – stellar dynamics

PDF

Cite this paper as:

Anisiu, M-C., Bozis, G, Two-dimension potentials which generate spatial families of orbits, Astronomische Nachrichten, vol.330, issue 4,  2009, pag.411-415

About this paper

Journal

Online Library

Publisher Name
Print ISSN

Not available yet.

Online ISSN

Not available yet.

Google Scholar Profile

References

Anisiu, M.-C.: 2005, Inverse Probl. Sci. Eng. 13, 545
Anisiu, M.-C., Bozis, G.: 2008, Didactica Mathematica 26, 9
Borghero, F., Melis, A.: 1990, CMDA 49, 273
Bozis, G.: 1983, CeMec 31, 129
Bozis, G., Kotoulas, T.A.: 2004, Rendiconti del Seminario della
Facolta di Scienze dell’ Universita di Cagliari 74, 83
Bozis, G., Kotoulas, T.A.: 2005, Inverse Problems 21, 343
Bozis, G., Nakhla, A.: 1986, CeMec 38, 357
Erdi, B.: 1982, CeMec 28, 209 ´
Kotoulas, T.A., Bozis, G.: 2006, JPhA 39, 9223
Melis, A., Borghero, F.: 1986, Mecc 21, 71
Mertens, R.: 1981, ZAMM 61, 252
Puel, F.: 1992, CMDA 53, 207
Shorokhov, S.G.: 1988, CeMec 44, 193
Szebehely, V.: 1974, in: G. Proverbio (ed.), Proc. of the Internat.
Meeting on Earth’s Rotation by Satellite Observation, p. 31
Varadi, F., Erdi, B.: 1983, CeMec 32, 395

Paper (preprint) in HTML form

2009-Anisiu_Astron-Two-dimension-potentials

Two-dimension potentials which generate spatial families of orbits

M.-C. Anisiu 1 , 1 , ^(1,***){ }^{1, \star}1, and G. Bozis 2 2 ^(2){ }^{2}2 1 1 ^(1){ }^{1}1 T. Popoviciu Institute of Numerical Analysis, P. O. Box 68, 400110 Cluj-Napoca, Romania 2 2 ^(2){ }^{2}2 Aristotle University of Thesaloniki, GR-54006, Greece

Received 2008 Apr 18, accepted 2008 Aug 14
Published online 2009 Apr 17
Key words celestial mechanics - stellar dynamics
We consider the following case of the 3D inverse problem of dynamics: Given a spatial two-parametric family of curves f ( x , y , z ) = c 1 , g ( x , y , z ) = c 2 f ( x , y , z ) = c 1 , g ( x , y , z ) = c 2 f(x,y,z)=c_(1),g(x,y,z)=c_(2)f(x, y, z)=c_{1}, g(x, y, z)=c_{2}f(x,y,z)=c1,g(x,y,z)=c2, find possibly existing two-dimension potentials under whose action the curves of the family are trajectories for a unit mass particle. First we establish the conditions which must be fulfilled by the family so that potentials of the form w ( y , z ) w ( y , z ) w(y,z)w(y, z)w(y,z) give rise to the curves of the family, and we present some applications. Then we examine briefly the existence of potentials depending on ( x , z ) ( x , z ) (x,z)(x, z)(x,z), respectively ( x , y ) ( x , y ) (x,y)(x, y)(x,y), which are compatible with the given family.

1 Introduction

The three-dimensional version of the inverse problem of dynamics has already a history of more than 25 years. Following Szebehely's (1974) presentation of the problem in two dimensions, Érdi (1982) was the first to address the 3D problem by considering a monoparametric family of spatial orbits. Seen from this viewpoint, the treatment of the problem is essentially similar to the study of the motion of a particle describing a monoparametric family of (spatial) orbits on a fixed surface (Mertens 1981). There followed studies which put the problem in a more accurate perspective: A two-parametric family of spatial curves is given and the potential giving rise to these curves is required (Bozis 1983; Váradi & Érdi 1983; Bozis & Nakhla 1986; Shorokhov 1988; Puel 1992). At almost the same period, the problem was generalized to account for holonomic systems with n n nnn degrees of freedom (Melis & Borghero 1986; Borghero & Melis 1990).
Assuming that a two-parametric set of orbits
(1) f ( x , y , z ) = c 1 , g ( x , y , z ) = c 2 , (1) f ( x , y , z ) = c 1 , g ( x , y , z ) = c 2 , {:(1)f(x","y","z)=c_(1)","quad g(x","y","z)=c_(2)",":}\begin{equation*} f(x, y, z)=c_{1}, \quad g(x, y, z)=c_{2}, \tag{1} \end{equation*}(1)f(x,y,z)=c1,g(x,y,z)=c2,
in the Oxyz space can be traced by a material point in the presence of an unknown potential V = V ( x , y , z ) V = V ( x , y , z ) V=V(x,y,z)V=V(x, y, z)V=V(x,y,z) with energy dependence function E = E ( f , g ) E = E ( f , g ) E=E(f,g)E=E(f, g)E=E(f,g) which is considered to be known, one aims to finding the potential. From its mathematical viewpoint, the problem consists in solving a linear system of two PDEs in the unique unknown function V = V ( x , y , z ) V = V ( x , y , z ) V=V(x,y,z)V=V(x, y, z)V=V(x,y,z). It was soon made clear that this problem has no solution unless the 'given' energy function E = E ( c 1 , c 2 ) E = E c 1 , c 2 E=E(c_(1),c_(2))E=E\left(c_{1}, c_{2}\right)E=E(c1,c2) satisfies certain conditions depending on the given functions f ( x , y , z ) f ( x , y , z ) f(x,y,z)f(x, y, z)f(x,y,z) and g ( x , y , z ) g ( x , y , z ) g(x,y,z)g(x, y, z)g(x,y,z) (Bozis & Nakhla
1986; Shorokhov 1988). In fact, if some specific conditions are satisfied, the problem admits of a unique solution V = V ( x , y , z ) V = V ( x , y , z ) V=V(x,y,z)V=V(x, y, z)V=V(x,y,z), up to an additive and a multiplicative constant.
In a manner analogous to that used in the 2D case, Bozis & Kotoulas (2005) and also Anisiu (2005) eliminated the energy and produced a system of two linear in V ( x , y , z ) V ( x , y , z ) V(x,y,z)V(x, y, z)V(x,y,z) PDEs, one of the first and one of the second order (Eqs. (5) and (6) below). If this system is compatible and can be solved for a given family of orbits (1), then the corresponding energy with which each member of the family is being traced can be found (Eq. (7) below).
In this paper we shall deal with these two energy-free PDEs, taking into account that, in spite of their linearity, the above system cannot be treated in a straightforward manner. Further assumptions, regarding either the form of the given orbits or the form of the required potential or both may simplify the problem. This practice was followed by Bozis & Kotoulas (2005) and by Kotoulas & Bozis (2006) and will be followed in the present paper also: We look for two-dimension potentials which can possibly give rise to 3D families of orbits (1), given in advance. We study in detail the case of potentials of the form V = w ( y , z ) V = w ( y , z ) V=w(y,z)V=w(y, z)V=w(y,z).
In Sect. 2 we give the general energy-free PDEs of the inverse problem. In Sect. 3 we write down these equations for the case at hand and we pursue their solution by establishing the necessary and sufficient conditions on the 'given' family so that such solutions w = w ( y , z ) w = w ( y , z ) w=w(y,z)w=w(y, z)w=w(y,z) do exist. If all types of 2 D 2 D 2D2 D2D potentials creating a given family are required, the cases w = w ( x , z ) w = w ( x , z ) w=w(x,z)w=w(x, z)w=w(x,z) and w = w ( x , y ) w = w ( x , y ) w=w(x,y)w=w(x, y)w=w(x,y) must be studied separately. This can be done, of course, in a manner quite similar to that followed here. In Sect. 4 we examine some special cases exempted at a first step, so that the analysis could continue. In Sect. 5 we apply the theory to two rather broad sets of spatial curves and enrich the library of the ex-
isting examples. Finally, in Sect. 6 we discuss shortly the case w = w ( x , z ) w = w ( x , z ) w=w(x,z)w=w(x, z)w=w(x,z) of 2D potentials and we comment on the case w = w ( x , y ) w = w ( x , y ) w=w(x,y)w=w(x, y)w=w(x,y).

2 The energy-free equations of the spatial inverse problem of dynamics

In a three-dimensional frame we deal with two-parametric families of orbits written in the form (1), which are in an one-to-one correspondence with a pair ( α , β α , β alpha,beta\alpha, \betaα,β ) of 'slope functions' defined by
α = f z g x f x g z f y g z f z g y , β = f x g y f y g x f y g z f z g y α = f z g x f x g z f y g z f z g y , β = f x g y f y g x f y g z f z g y alpha=(f_(z)g_(x)-f_(x)g_(z))/(f_(y)g_(z)-f_(z)g_(y)),quad beta=(f_(x)g_(y)-f_(y)g_(x))/(f_(y)g_(z)-f_(z)g_(y))\alpha=\frac{f_{z} g_{x}-f_{x} g_{z}}{f_{y} g_{z}-f_{z} g_{y}}, \quad \beta=\frac{f_{x} g_{y}-f_{y} g_{x}}{f_{y} g_{z}-f_{z} g_{y}}α=fzgxfxgzfygzfzgy,β=fxgyfygxfygzfzgy.
The indices denote partial derivatives. Equations (1) then constitute the general solution of the ODE system
d y d x = α ( x , y , z ) , d z d x = β ( x , y , z ) d y d x = α ( x , y , z ) , d z d x = β ( x , y , z ) (dy)/((d)x)=alpha(x,y,z),quad(dz)/((d)x)=beta(x,y,z)\frac{\mathrm{d} y}{\mathrm{~d} x}=\alpha(x, y, z), \quad \frac{\mathrm{d} z}{\mathrm{~d} x}=\beta(x, y, z)dy dx=α(x,y,z),dz dx=β(x,y,z).
There exist in this case two energy-free PDEs (one of the first and one of the second order) relating families ( α , β α , β alpha,beta\alpha, \betaα,β ) and potentials V ( x , y , z ) V ( x , y , z ) V(x,y,z)V(x, y, z)V(x,y,z) (Bozis & Kotoulas 2005; Anisiu 2005).
Let us assume that α 0 0 α 0 0 alpha_(0)!=0\alpha_{0} \neq 0α00 and adopt the notation
ϵ ¯ = ( 1 , α , β ) , α 0 = ϵ ¯ grad α , β 0 = ϵ ¯ grad β ϵ ¯ = ( 1 , α , β ) , α 0 = ϵ ¯ grad α , β 0 = ϵ ¯ grad β bar(epsilon)=(1,alpha,beta),quadalpha_(0)= bar(epsilon)grad alpha,quadbeta_(0)= bar(epsilon)grad beta\bar{\epsilon}=(1, \alpha, \beta), \quad \alpha_{0}=\bar{\epsilon} \operatorname{grad} \alpha, \quad \beta_{0}=\bar{\epsilon} \operatorname{grad} \betaϵ¯=(1,α,β),α0=ϵ¯gradα,β0=ϵ¯gradβ,
Θ = 1 + α 2 + β 2 , n = Θ α 0 , n 0 = ϵ ¯ grad n Θ = 1 + α 2 + β 2 , n = Θ α 0 , n 0 = ϵ ¯ grad n Theta=1+alpha^(2)+beta^(2),quad n=(Theta)/(alpha_(0)),quadn_(0)= bar(epsilon)grad n\Theta=1+\alpha^{2}+\beta^{2}, \quad n=\frac{\Theta}{\alpha_{0}}, \quad n_{0}=\bar{\epsilon} \operatorname{grad} nΘ=1+α2+β2,n=Θα0,n0=ϵ¯gradn.
The two equations satisfied by the potential in the 3D inverse problem are
( α β 0 α 0 β ) V x β 0 V y + α 0 V z = 0 α β 0 α 0 β V x β 0 V y + α 0 V z = 0 (alphabeta_(0)-alpha_(0)beta)V_(x)-beta_(0)V_(y)+alpha_(0)V_(z)=0\left(\alpha \beta_{0}-\alpha_{0} \beta\right) V_{x}-\beta_{0} V_{y}+\alpha_{0} V_{z}=0(αβ0α0β)Vxβ0Vy+α0Vz=0,
(5) α V x x + ( α 2 1 ) V x y + α β V x z α V y y β V y z + (5) α V x x + α 2 1 V x y + α β V x z α V y y β V y z + {:(5)alphaV_(xx)+(alpha^(2)-1)V_(xy)+alpha betaV_(xz)-alphaV_(yy)-betaV_(yz)+:}\begin{equation*} \alpha V_{x x}+\left(\alpha^{2}-1\right) V_{x y}+\alpha \beta V_{x z}-\alpha V_{y y}-\beta V_{y z}+ \tag{5} \end{equation*}(5)αVxx+(α21)Vxy+αβVxzαVyyβVyz+
1 n ( ( 2 + α n 0 + α 0 n ) V x + ( 2 α n 0 ) V y + 2 β V z ) = 0 1 n 2 + α n 0 + α 0 n V x + 2 α n 0 V y + 2 β V z = 0 (1)/(n)((2+alphan_(0)+alpha_(0)n)V_(x)+(2alpha-n_(0))V_(y)+2betaV_(z))=0\frac{1}{n}\left(\left(2+\alpha n_{0}+\alpha_{0} n\right) V_{x}+\left(2 \alpha-n_{0}\right) V_{y}+2 \beta V_{z}\right)=01n((2+αn0+α0n)Vx+(2αn0)Vy+2βVz)=0.
For any compatible pair of potential V ( x , y , z ) V ( x , y , z ) V(x,y,z)V(x, y, z)V(x,y,z) and orbit ( α , β α , β alpha,beta\alpha, \betaα,β ) the energy is given by
E ( f , g ) = V + Θ ( α V x V y ) / 2 α 0 E ( f , g ) = V + Θ α V x V y / 2 α 0 E(f,g)=V+Theta(alphaV_(x)-V_(y))//2alpha_(0)E(f, g)=V+\Theta\left(\alpha V_{x}-V_{y}\right) / 2 \alpha_{0}E(f,g)=V+Θ(αVxVy)/2α0,
and real motion is allowed in the region
α V x V y α 0 0 α V x V y α 0 0 (alphaV_(x)-V_(y))/(alpha_(0)) >= 0\frac{\alpha V_{x}-V_{y}}{\alpha_{0}} \geq 0αVxVyα00.
The potentials producing families of straight lines, for which both α 0 α 0 alpha_(0)\alpha_{0}α0 and β 0 β 0 beta_(0)\beta_{0}β0 are identically zero, must satisfy (Bozis & Kotoulas 2004) the equations
α V x V y = 0 , β V x V z = 0 α V x V y = 0 , β V x V z = 0 alphaV_(x)-V_(y)=0,quad betaV_(x)-V_(z)=0\alpha V_{x}-V_{y}=0, \quad \beta V_{x}-V_{z}=0αVxVy=0,βVxVz=0.
If α 0 = 0 α 0 = 0 alpha_(0)=0\alpha_{0}=0α0=0 but β 0 0 β 0 0 beta_(0)!=0\beta_{0} \neq 0β00 equation (6) is replaced by a second-order differential equation presented in Anisiu (2005) and Kotoulas & Bozis (2006).

3 Potentials V = w ( y , z ) V = w ( y , z ) V=w(y,z)V=w(y, z)V=w(y,z) which generate spatial families of orbits

We shall consider exclusively potentials of the form
V = w ( y , z ) V = w ( y , z ) V=w(y,z)V=w(y, z)V=w(y,z).
The material point moving under the action of such a potential will have a uniform motion in x x xxx, i. e. x ˙ ( t ) = x ˙ ( t ) = x^(˙)(t)=\dot{x}(t)=x˙(t)= const. As a consequence, the orbits cannot be closed. For V V VVV given in (10) the system of the two PDEs (5) and (6) becomes
(11) w z = G w y (11) w z = G w y {:(11)w_(z)=Gw_(y):}\begin{equation*} w_{z}=G w_{y} \tag{11} \end{equation*}(11)wz=Gwy
and
Θ ( α + β G ) w y y + Ψ w y = 0 Θ ( α + β G ) w y y + Ψ w y = 0 Theta(alpha+beta G)w_(yy)+Psiw_(y)=0\Theta(\alpha+\beta G) w_{y y}+\Psi w_{y}=0Θ(α+βG)wyy+Ψwy=0,
where
(13) G = β 0 α 0 (13) G = β 0 α 0 {:(13)G=(beta_(0))/(alpha_(0)):}\begin{equation*} G=\frac{\beta_{0}}{\alpha_{0}} \tag{13} \end{equation*}(13)G=β0α0
and
Ψ = β Θ G y + α 0 ( n 0 2 ( α + β G ) ) Ψ = β Θ G y + α 0 n 0 2 ( α + β G ) Psi=beta ThetaG_(y)+alpha_(0)(n_(0)-2(alpha+beta G))\Psi=\beta \Theta G_{y}+\alpha_{0}\left(n_{0}-2(\alpha+\beta G)\right)Ψ=βΘGy+α0(n02(α+βG)).
Equation (11) comes directly from (5) but, to obtain (12) from (6), we used (11) and, with its help, we expressed w y z w y z w_(yz)w_{y z}wyz and w z z w z z w_(zz)w_{z z}wzz in terms of w y w y w_(y)w_{y}wy and w y y w y y w_(yy)w_{y y}wyy.
It is seen from Eq. (11) that the function G G GGG must be independent of x x xxx, i. e.
G x = 0 G x = 0 G_(x)=0G_{x}=0Gx=0.
We assume that α + β G 0 α + β G 0 alpha+beta G!=0\alpha+\beta G \neq 0α+βG0 and we put
H = Ψ Θ ( α + β G ) H = Ψ Θ ( α + β G ) H=(Psi)/(Theta(alpha+beta G))H=\frac{\Psi}{\Theta(\alpha+\beta G)}H=ΨΘ(α+βG).
For the PDE (12), now written as
w y y + H w y = 0 w y y + H w y = 0 w_(yy)+Hw_(y)=0w_{y y}+H w_{y}=0wyy+Hwy=0,
to have a solution of the form (10) it must be
H x = 0 H x = 0 H_(x)=0H_{x}=0Hx=0.
Solving Eq. (12) for w y w y w_(y)w_{y}wy we obtain
(19) w y = D ( z ) exp ( y H ( u , z ) d u ) (19) w y = D ( z ) exp y H ( u , z ) d u {:(19)w_(y)=D(z)exp(-int^(y)H(u,z)du):}\begin{equation*} w_{y}=D(z) \exp \left(-\int^{y} H(u, z) \mathrm{d} u\right) \tag{19} \end{equation*}(19)wy=D(z)exp(yH(u,z)du)
(with D ( z ) D ( z ) D(z)D(z)D(z) momentarily arbitrary function of z z zzz ), and from Eq. (11) we obtain w z w z w_(z)w_{z}wz. The compatibility condition ( w y z = w z y ) w y z = w z y (w_(yz)=:}{:w_(zy))\left(w_{y z}=\right. \left.w_{z y}\right)(wyz=wzy) for w y w y w_(y)w_{y}wy and w z w z w_(z)w_{z}wz, as these are given by Eqs. (19) and (11) respectively, is
D ( z ) D ( z ) J = 0 D ( z ) D ( z ) J = 0 D^(')(z)-D(z)J=0D^{\prime}(z)-D(z) J=0D(z)D(z)J=0,
where the function J J JJJ is given by
J = G y ( y , z ) G ( y , z ) H ( y , z ) + y H z ( u , z ) d u J = G y ( y , z ) G ( y , z ) H ( y , z ) + y H z ( u , z ) d u J=G_(y)(y,z)-G(y,z)H(y,z)+int^(y)H_(z)(u,z)duJ=G_{y}(y, z)-G(y, z) H(y, z)+\int^{y} H_{z}(u, z) \mathrm{d} uJ=Gy(y,z)G(y,z)H(y,z)+yHz(u,z)du.
The function J J JJJ must depend merely on z z zzz, i. e. it must be
J y = 0 J y = 0 J_(y)=0J_{y}=0Jy=0.
This last condition (22) can also be written as
G y y G y H G H y + H z = 0 G y y G y H G H y + H z = 0 G_(yy)-G_(y)H-GH_(y)+H_(z)=0G_{y y}-G_{y} H-G H_{y}+H_{z}=0GyyGyHGHy+Hz=0.
For a given family (1) the 'integration function' D ( z ) D ( z ) D(z)D(z)D(z) is specified from (20). It is
D ( z ) = d exp ( J ( z ) d z ) D ( z ) = d exp J ( z ) d z D(z)=d exp(int J(z)dz)D(z)=d \exp \left(\int J(z) \mathrm{d} z\right)D(z)=dexp(J(z)dz),
where J ( z ) J ( z ) J(z)J(z)J(z) is given by (21) and d d ddd is a constant. From Eqs. (7) and (8) we conclude that real motion is allowed in the region defined by the inequality
w y α 0 0 w y α 0 0 (w_(y))/(alpha_(0)) <= 0\frac{w_{y}}{\alpha_{0}} \leq 0wyα00,
and the particle moves with total energy
E = w ( y , z ) Θ w y 2 α 0 E = w ( y , z ) Θ w y 2 α 0 E=w(y,z)-(Thetaw_(y))/(2alpha_(0))E=w(y, z)-\frac{\Theta w_{y}}{2 \alpha_{0}}E=w(y,z)Θwy2α0.
We summarize by the following
Proposition 1. For α 0 0 , α + β G 0 α 0 0 , α + β G 0 alpha_(0)!=0,alpha+beta G!=0\alpha_{0} \neq 0, \alpha+\beta G \neq 0α00,α+βG0 and for any family ( α , β ) ( α , β ) (alpha,beta)(\alpha, \beta)(α,β) satisfying the conditions (15), (18) and (22) (equivalently 23), there exists a two-dimension compatible potential w = w ( y , z ) w = w ( y , z ) w=w(y,z)w=w(y, z)w=w(y,z). The potential is given by the (compatible) Eqs. (19) and (11) with D ( z ) D ( z ) D(z)D(z)D(z) given by (24). The family is lying in the region (25) and is traced with total energy given by Eq. (26).

4 Special cases

a. α 0 = 0 , β 0 = 0 α 0 = 0 , β 0 = 0 quadalpha_(0)=0,beta_(0)=0\quad \alpha_{0}=0, \beta_{0}=0α0=0,β0=0
The case is trivial. Only if the given orbits are families of straight lines we are led to α 0 = β 0 = 0 α 0 = β 0 = 0 alpha_(0)=beta_(0)=0\alpha_{0}=\beta_{0}=0α0=β0=0 and the corresponding potential obtained from (9) is w ( y , z ) = w ( y , z ) = w(y,z)=w(y, z)=w(y,z)= const.
b. α 0 = 0 , β 0 0 α 0 = 0 , β 0 0 quadalpha_(0)=0,beta_(0)!=0\quad \alpha_{0}=0, \beta_{0} \neq 0α0=0,β00
The Eq. (5) is still valid but Eq. (6) must be replaced by another second order PDE, as mentioned at the end of Sect. 2. Here we write down directly the pertinent system of the Eqs. (5) and (6) for 2D potentials of the form w = w ( y , z ) w = w ( y , z ) w=w(y,z)w=w(y, z)w=w(y,z). This system reads
w y = 0 , w z z = 2 β n ~ 0 n ~ β w z w y = 0 , w z z = 2 β n ~ 0 n ~ β w z w_(y)=0,w_(zz)=(2beta- widetilde(n)_(0))/(( widetilde(n))beta)w_(z)w_{y}=0, w_{z z}=\frac{2 \beta-\widetilde{n}_{0}}{\widetilde{n} \beta} w_{z}wy=0,wzz=2βn~0n~βwz,
where
n ~ = Θ β 0 , n ~ 0 = n ~ x + α n ~ y + β n ~ z n ~ = Θ β 0 , n ~ 0 = n ~ x + α n ~ y + β n ~ z widetilde(n)=(Theta)/(beta_(0)), widetilde(n)_(0)= widetilde(n)_(x)+alpha widetilde(n)_(y)+beta widetilde(n)_(z)\widetilde{n}=\frac{\Theta}{\beta_{0}}, \widetilde{n}_{0}=\widetilde{n}_{x}+\alpha \widetilde{n}_{y}+\beta \widetilde{n}_{z}n~=Θβ0,n~0=n~x+αn~y+βn~z.
The system (27) has as solution a 1D potential w = w ( z ) w = w ( z ) w=w(z)w= w(z)w=w(z), provided that the ratio ( 2 β n ~ 0 ) / ( n ~ β ) 2 β n ~ 0 / ( n ~ β ) (2beta- widetilde(n)_(0))//( widetilde(n)beta)\left(2 \beta-\widetilde{n}_{0}\right) /(\widetilde{n} \beta)(2βn~0)/(n~β) is independent of x x xxx and y y yyy.
c. α 0 0 , α + β G = 0 , Ψ 0 α 0 0 , α + β G = 0 , Ψ 0 quadalpha_(0)!=0,alpha+beta G=0,Psi!=0\quad \alpha_{0} \neq 0, \alpha+\beta G=0, \Psi \neq 0α00,α+βG=0,Ψ0
The case is trivial. From (12) we obtain w y = 0 w y = 0 w_(y)=0w_{y}=0wy=0, and from (11) w z = 0 w z = 0 w_(z)=0w_{z}=0wz=0 i. e. w = w = w=w=w= const.
d. α 0 0 , α + β G = 0 , Ψ = 0 α 0 0 , α + β G = 0 , Ψ = 0 quadalpha_(0)!=0,alpha+beta G=0,Psi=0\quad \alpha_{0} \neq 0, \alpha+\beta G=0, \Psi=0α00,α+βG=0,Ψ=0
Then Eq. (12) is satisfied identically. The potential is found from (11), provided that the condition (15) holds, and it will not be uniquely determined.
Although at first sight questionable, the equations
(29) G x = 0 , α + β G = 0 , n β G y + n 0 = 0 (29) G x = 0 , α + β G = 0 , n β G y + n 0 = 0 {:(29)G_(x)=0","quad alpha+beta G=0","quad n betaG_(y)+n_(0)=0:}\begin{equation*} G_{x}=0, \quad \alpha+\beta G=0, \quad n \beta G_{y}+n_{0}=0 \tag{29} \end{equation*}(29)Gx=0,α+βG=0,nβGy+n0=0
may hold simultaneously; for example, for the family f = y sin x + z cos x = c 1 , g = z sin x + y cos x = c 2 f = y sin x + z cos x = c 1 , g = z sin x + y cos x = c 2 f=y sin x+z cos x=c_(1),g=-z sin x+y cos x=c_(2)f= y \sin x+z \cos x=c_{1}, g=-z \sin x+y \cos x=c_{2}f=ysinx+zcosx=c1,g=zsinx+ycosx=c2 with the slope functions α = z α = z alpha=z\alpha=zα=z and β = y β = y beta=-y\beta=-yβ=y, it is α α 0 + β β 0 = 0 α α 0 + β β 0 = 0 alphaalpha_(0)+betabeta_(0)=0\alpha \alpha_{0}+\beta \beta_{0}=0αα0+ββ0=0 and Ψ = 0 Ψ = 0 Psi=0\Psi=0Ψ=0. Equation (11) can be easily integrated and gives V = F ( y 2 + z 2 ) V = F y 2 + z 2 V=F(y^(2)+z^(2))V=F\left(y^{2}+z^{2}\right)V=F(y2+z2), with F F FFF an arbitrary function. The energy is E = ( 1 + f 2 + g 2 ) F ( y 2 + z 2 ) + F ( y 2 + z 2 ) E = 1 + f 2 + g 2 F y 2 + z 2 + F y 2 + z 2 E=(1+f^(2)+g^(2))F^(')(y^(2)+z^(2))+F(y^(2)+z^(2))E=\left(1+f^{2}+g^{2}\right) F^{\prime}\left(y^{2}+z^{2}\right)+F\left(y^{2}+z^{2}\right)E=(1+f2+g2)F(y2+z2)+F(y2+z2), and the allowed region F ( y 2 + z 2 ) 0 F y 2 + z 2 0 F^(')(y^(2)+z^(2)) >= 0F^{\prime}\left(y^{2}+z^{2}\right) \geq 0F(y2+z2)0.
Proposition 2. If one of the conditions a a a\boldsymbol{a}a or c c c\boldsymbol{c}c from above holds, the problem has only the trivial solution w ( y , z ) = w ( y , z ) = w(y,z)=w(y, z)=w(y,z)= const.
If condition b b b\boldsymbol{b}b is fulfilled and the ratio ( 2 β n ~ 0 ) / ( n ~ β ) 2 β n ~ 0 / ( n ~ β ) (2beta- widetilde(n)_(0))//( widetilde(n)beta)\left(2 \beta-\widetilde{n}_{0}\right) /(\widetilde{n} \beta)(2βn~0)/(n~β) is independent of x x xxx and y y yyy, the potential will depend only on the z z zzz-variable, otherwise no potential can be found.
Finally, if the given family satisfies condition d d d\boldsymbol{d}d and G x = 0 G x = 0 G_(x)=0G_{x}=0Gx=0, we shall obtain a family of potentials w ( y , z ) w ( y , z ) w(y,z)w(y, z)w(y,z); if G x 0 G x 0 G_(x)!=0G_{x} \neq 0Gx0 there will be no such potential.

5 Applications

Application 1 Let f 1 , f 2 , g 1 , g 2 f 1 , f 2 , g 1 , g 2 f_(1),f_(2),g_(1),g_(2)f_{1}, f_{2}, g_{1}, g_{2}f1,f2,g1,g2 be constants such that
(30) f 1 ( 1 + g 2 ) = g 1 ( 1 + f 2 ) (30) f 1 1 + g 2 = g 1 1 + f 2 {:(30)f_(1)(1+g_(2))=g_(1)(1+f_(2)):}\begin{equation*} f_{1}\left(1+g_{2}\right)=g_{1}\left(1+f_{2}\right) \tag{30} \end{equation*}(30)f1(1+g2)=g1(1+f2)
and
g 1 ( f 2 g 2 1 ) ( 1 + g 1 + g 2 ) 0 g 1 f 2 g 2 1 1 + g 1 + g 2 0 g_(1)(f_(2)g_(2)-1)(1+g_(1)+g_(2))!=0g_{1}\left(f_{2} g_{2}-1\right)\left(1+g_{1}+g_{2}\right) \neq 0g1(f2g21)(1+g1+g2)0.
For any specific values of these constants (satisfying Eq. 30) we consider the family of curves in the inertial frame Oxyz
f = x f 1 y z f 2 = c 1 , g = x g 1 y g 2 z = c 2 f = x f 1 y z f 2 = c 1 , g = x g 1 y g 2 z = c 2 f=x^(f_(1))yz^(f_(2))=c_(1),quad g=x^(g_(1))y^(g_(2))z=c_(2)f=x^{f_{1}} y z^{f_{2}}=c_{1}, \quad g=x^{g_{1}} y^{g_{2}} z=c_{2}f=xf1yzf2=c1,g=xg1yg2z=c2,
where c 1 , c 2 c 1 , c 2 c_(1),c_(2)c_{1}, c_{2}c1,c2 are parameters. It can be checked that (32) constitutes a spatial two-parametric set of orbits produced by the 2D potential field
(33) w ( y , z ) = d 1 ( y 2 + z 2 ) 1 + g 1 + g 2 g 1 (33) w ( y , z ) = d 1 y 2 + z 2 1 + g 1 + g 2 g 1 {:(33)w(y","z)=d_(1)(y^(2)+z^(2))^((1+g_(1)+g_(2))/(g_(1))):}\begin{equation*} w(y, z)=d_{1}\left(y^{2}+z^{2}\right)^{\frac{1+g_{1}+g_{2}}{g_{1}}} \tag{33} \end{equation*}(33)w(y,z)=d1(y2+z2)1+g1+g2g1
traced with total energy
E = d 1 ( 1 + g 2 g 1 ) 2 ( f k 0 g 2 g k 0 + f k 0 g k 0 f 2 ) 1 + g 2 g 1 E = d 1 1 + g 2 g 1 2 f k 0 g 2 g k 0 + f k 0 g k 0 f 2 1 + g 2 g 1 E=-d_(1)((1+g_(2))/(g_(1)))^(2)(f^(k_(0)g_(2))g^(-k_(0))+f^(-k_(0))g^(k_(0)f_(2)))^((1+g_(2))/(g_(1)))E=-d_{1}\left(\frac{1+g_{2}}{g_{1}}\right)^{2}\left(f^{k_{0} g_{2}} g^{-k_{0}}+f^{-k_{0}} g^{k_{0} f_{2}}\right)^{\frac{1+g_{2}}{g_{1}}}E=d1(1+g2g1)2(fk0g2gk0+fk0gk0f2)1+g2g1,
where
k 0 = 2 f 2 g 2 1 k 0 = 2 f 2 g 2 1 k_(0)=(2)/(f_(2)g_(2)-1)k_{0}=\frac{2}{f_{2} g_{2}-1}k0=2f2g21,
and d 1 d 1 d_(1)d_{1}d1 is an arbitrary constant.
Remarks:
  1. In the formulae (32), the variable y y yyy in f f fff and the variable z z zzz in g g ggg were taken to be raised to the power 1. Obviously, if the nonzero exponents are not so, they can always be done equal to one.
  2. The condition g 1 0 g 1 0 g_(1)!=0g_{1} \neq 0g10 does not reduce the generality of the family (32), because the x x xxx-variable must appear at least in one of f f fff and g g ggg, hence we can suppose that g 1 g 1 g_(1)g_{1}g1 is nonzero (otherwise we may perform the change f g f , y z y f g f , y z y f rarr g rarr f,y rarr z rarr yf \rightarrow g \rightarrow f, y \rightarrow z \rightarrow yfgf,yzy ).
  3. The condition f 2 g 2 1 0 f 2 g 2 1 0 f_(2)g_(2)-1!=0f_{2} g_{2}-1 \neq 0f2g210 ensures the fact that α , β α , β alpha,beta\alpha, \betaα,β, α 0 , β 0 α 0 , β 0 alpha_(0),beta_(0)\alpha_{0}, \beta_{0}α0,β0 and k 0 k 0 k_(0)k_{0}k0 are meaningfull, that is the denominator of these expressions is nonzero. It also guarantees that the two surfaces in (32) do not coincide.
  4. For the family (32) we get (taking into account condition 31)
    α 0 = r ( r 1 ) y x 2 , β 0 = r ( r 1 ) z x 2 α 0 = r ( r 1 ) y x 2 , β 0 = r ( r 1 ) z x 2 alpha_(0)=r(r-1)(y)/(x^(2)),quadbeta_(0)=r(r-1)(z)/(x^(2))\alpha_{0}=r(r-1) \frac{y}{x^{2}}, \quad \beta_{0}=r(r-1) \frac{z}{x^{2}}α0=r(r1)yx2,β0=r(r1)zx2,
    where
    r = g 1 f 2 f 1 1 f 2 g 2 r = g 1 f 2 f 1 1 f 2 g 2 r=(g_(1)f_(2)-f_(1))/(1-f_(2)g_(2))r=\frac{g_{1} f_{2}-f_{1}}{1-f_{2} g_{2}}r=g1f2f11f2g2.
It is r 0 r 0 r!=0r \neq 0r0 because of g 1 ( f 2 g 2 1 ) 0 g 1 f 2 g 2 1 0 g_(1)(f_(2)g_(2)-1)!=0g_{1}\left(f_{2} g_{2}-1\right) \neq 0g1(f2g21)0 (if g 1 f 2 f 1 = 0 g 1 f 2 f 1 = 0 g_(1)f_(2)-f_(1)=0g_{1} f_{2}-f_{1}=0g1f2f1=0, then f 2 g 2 1 = ( f 1 g 2 g 1 ) / g 1 = ( g 1 f 2 f 1 ) / g 1 = 0 f 2 g 2 1 = f 1 g 2 g 1 / g 1 = g 1 f 2 f 1 / g 1 = 0 f_(2)g_(2)-1=(f_(1)g_(2)-g_(1))//g_(1)=(g_(1)f_(2)-f_(1))//g_(1)=0f_{2} g_{2}-1=\left(f_{1} g_{2}-g_{1}\right) / g_{1}=\left(g_{1} f_{2}-f_{1}\right) / g_{1}=0f2g21=(f1g2g1)/g1=(g1f2f1)/g1=0, contradiction). It is also r 1 r 1 r!=1r \neq 1r1 because of g 1 ( 1 + g 1 + g 2 ) 0 g 1 1 + g 1 + g 2 0 g_(1)(1+g_(1)+:}{:g_(2))!=0g_{1}\left(1+g_{1}+\right. \left.g_{2}\right) \neq 0g1(1+g1+g2)0, since in view of (31) f 2 ( g 1 + g 2 ) ( f 1 + 1 ) = ( g 1 f 2 f 1 ) ( 1 + g 1 + g 2 ) / g 1 f 2 g 1 + g 2 f 1 + 1 = g 1 f 2 f 1 1 + g 1 + g 2 / g 1 f_(2)(g_(1)+g_(2))-(f_(1)+1)=(g_(1)f_(2)-f_(1))(1+g_(1)+g_(2))//g_(1)f_{2}\left(g_{1}+g_{2}\right)-\left(f_{1}+1\right)= \left(g_{1} f_{2}-f_{1}\right)\left(1+g_{1}+g_{2}\right) / g_{1}f2(g1+g2)(f1+1)=(g1f2f1)(1+g1+g2)/g1.
5. The case f 1 = 0 f 1 = 0 f_(1)=0f_{1}=0f1=0 is not excluded, and in view of (30) it must be f 2 = 1 f 2 = 1 f_(2)=-1f_{2}=-1f2=1. The pertinent family f = y z 1 = c 1 , g = x g 1 1 y g 2 z = c 2 f = y z 1 = c 1 , g = x g 1 1 y g 2 z = c 2 f=yz^(-1)=c_(1),g=x^(g_(1_(1)))y^(g_(2))z=c_(2)f=y z^{-1}= c_{1}, g=x^{g_{1_{1}}} y^{g_{2}} z=c_{2}f=yz1=c1,g=xg11yg2z=c2 leads to potential and energy found directly from (33) and (34) for f 1 = 0 , f 2 = 1 f 1 = 0 , f 2 = 1 f_(1)=0,f_(2)=-1f_{1}=0, f_{2}=-1f1=0,f2=1.
6 . The exponents f 1 , f 2 f 1 , f 2 f_(1),f_(2)f_{1}, f_{2}f1,f2 do not appear in the potential (33), so, in view of (30) also, we understand that the family (32) is essentially three-parametric.
7. As seen from (25), real motion is allowed all over the space for negative d 1 d 1 d_(1)d_{1}d1 and nowhere for positive d 1 d 1 d_(1)d_{1}d1.
Application 2 It can be checked that, for any values of the constants f 1 , g 1 , g 2 f 1 , g 1 , g 2 f_(1),g_(1),g_(2)f_{1}, g_{1}, g_{2}f1,g1,g2 such that not both f 1 f 1 f_(1)f_{1}f1 and g 1 g 1 g_(1)g_{1}g1 are zero, the family
(38) f = f 1 x 2 + y 2 + z 2 = c 1 , g = g 1 x + g 2 y + z = c 2 (38) f = f 1 x 2 + y 2 + z 2 = c 1 , g = g 1 x + g 2 y + z = c 2 {:(38)f=f_(1)x^(2)+y^(2)+z^(2)=c_(1)","quad g=g_(1)x+g_(2)y+z=c_(2):}\begin{equation*} f=f_{1} x^{2}+y^{2}+z^{2}=c_{1}, \quad g=g_{1} x+g_{2} y+z=c_{2} \tag{38} \end{equation*}(38)f=f1x2+y2+z2=c1,g=g1x+g2y+z=c2
is produced by the potential
(39) w ( y , z ) = d 2 ( y g 2 z ) 2 (39) w ( y , z ) = d 2 y g 2 z 2 {:(39)w(y","z)=(d_(2))/((y-g_(2)z)^(2)):}\begin{equation*} w(y, z)=\frac{d_{2}}{\left(y-g_{2} z\right)^{2}} \tag{39} \end{equation*}(39)w(y,z)=d2(yg2z)2
and is traced with energy
E = d 2 ( f 1 1 ) c 1 g 1 2 + f 1 ( c 1 + c 1 g 2 2 c 2 2 ) E = d 2 f 1 1 c 1 g 1 2 + f 1 c 1 + c 1 g 2 2 c 2 2 E=(d_(2)(f_(1)-1))/(c_(1)g_(1)^(2)+f_(1)(c_(1)+c_(1)g_(2)^(2)-c_(2)^(2)))E=\frac{d_{2}\left(f_{1}-1\right)}{c_{1} g_{1}^{2}+f_{1}\left(c_{1}+c_{1} g_{2}^{2}-c_{2}^{2}\right)}E=d2(f11)c1g12+f1(c1+c1g22c22),
d 2 d 2 d_(2)d_{2}d2 being an arbitrary constant.
Remarks:
  1. The first of Eqs. (38) is a family of surfaces of revolution around the x x xxx-axis. (For f 1 > 0 , c 1 > 0 f 1 > 0 , c 1 > 0 f_(1) > 0,c_(1) > 0f_{1}>0, c_{1}>0f1>0,c1>0 we have an ellipsoid of revolution, for f 1 < 0 f 1 < 0 f_(1) < 0f_{1}<0f1<0 we have a hyperboloid of one or two sheets, depending on whether c 1 > 0 c 1 > 0 c_(1) > 0c_{1}>0c1>0 or c 1 < 0 c 1 < 0 c_(1) < 0c_{1}<0c1<0 respectively.)
  2. The (integrable) potential (39) is independent of f 1 f 1 f_(1)f_{1}f1 and g 1 g 1 g_(1)g_{1}g1, so, the family (38) is essentially four-parametric. Notice, however, that, besides c 1 c 1 c_(1)c_{1}c1 and c 2 c 2 c_(2)c_{2}c2, the total energy depends on the constants g 1 g 1 g_(1)g_{1}g1 and g 2 g 2 g_(2)g_{2}g2.
  3. Depending on the values of the constants and the parameters involved and also on the sign in front of the potential (39), real motion of the massive particle takes place either everywhere (nowhere) in the 3D space or inside (outside) a cylinder whose generatrice is parallel to the x x xxx-axis. This remark is valid for all potentials of the form w = w ( y , z ) w = w ( y , z ) w=w(y,z)w=w(y, z)w=w(y,z) which we study here.

6 Concluding remarks

For any family of spatial curves (1), we established the conditions which, if fulfilled, allow for the existence of a compatible 2D potential V = w ( y , z ) V = w ( y , z ) V=w(y,z)V=w(y, z)V=w(y,z) of the form (10). These are
(41) G x = 0 , H x = 0 , J y = 0 (41) G x = 0 , H x = 0 , J y = 0 {:(41)G_(x)=0","quadH_(x)=0","quadJ_(y)=0:}\begin{equation*} G_{x}=0, \quad H_{x}=0, \quad J_{y}=0 \tag{41} \end{equation*}(41)Gx=0,Hx=0,Jy=0
with the functions G , H , J G , H , J G,H,JG, H, JG,H,J given by (13), (16), (21) respectively, in terms of the given family.
For a specific family (1), the above conditions (41) are not expected to be satisfied, of course, meaning that no compatible potential of the form (10) exists. Yet, for that specific family, a compatible 2D potential may exist e. g. of the form
V = w ( x , z ) V = w ( x , z ) V=w^(**)(x,z)V=w^{*}(x, z)V=w(x,z).
Working with (42) and following the same steps as in Sect. 3, we come up with the following three conditions for the family (1):
G y = 0 , H y = 0 , J x = 0 G y = 0 , H y = 0 , J x = 0 G_(y)^(**)=0,quadH_(y)^(**)=0,quadJ_(x)^(**)=0G_{y}^{*}=0, \quad H_{y}^{*}=0, \quad J_{x}^{*}=0Gy=0,Hy=0,Jx=0,
where
G = β α β 0 α 0 G = β α β 0 α 0 G^(**)=beta-alpha(beta_(0))/(alpha_(0))G^{*}=\beta-\alpha \frac{\beta_{0}}{\alpha_{0}}G=βαβ0α0,
H = Θ α β G x + α 0 ( 2 + α n 0 + α 0 n + 2 β G ) α Θ ( 1 + β G ) H = Θ α β G x + α 0 2 + α n 0 + α 0 n + 2 β G α Θ 1 + β G H^(**)=(Theta alpha betaG_(x)^(**)+alpha_(0)(2+alphan_(0)+alpha_(0)n+2betaG^(**)))/(alpha Theta(1+betaG^(**)))H^{*}=\frac{\Theta \alpha \beta G_{x}^{*}+\alpha_{0}\left(2+\alpha n_{0}+\alpha_{0} n+2 \beta G^{*}\right)}{\alpha \Theta\left(1+\beta G^{*}\right)}H=ΘαβGx+α0(2+αn0+α0n+2βG)αΘ(1+βG),
J = G x G H + x H z ( u , z ) d u J = G x G H + x H z ( u , z ) d u J^(**)=G_(x)^(**)-G^(**)H^(**)+int^(x)H_(z)^(**)(u,z)duJ^{*}=G_{x}^{*}-G^{*} H^{*}+\int^{x} H_{z}^{*}(u, z) \mathrm{d} uJ=GxGH+xHz(u,z)du.
The two PDEs (5) and (6) become
w z = G w x w z = G w x w_(z)^(**)=G^(**)w_(x)^(**)w_{z}^{*}=G^{*} w_{x}^{*}wz=Gwx,
(45) α Θ ( 1 + β G ) w x x + ( Θ α β G x (45) α Θ 1 + β G w x x + Θ α β G x {:(45)alpha Theta(1+betaG^(**))w_(xx)^(**)+(Theta alpha betaG_(x)^(**):}:}\begin{equation*} \alpha \Theta\left(1+\beta G^{*}\right) w_{x x}^{*}+\left(\Theta \alpha \beta G_{x}^{*}\right. \tag{45} \end{equation*}(45)αΘ(1+βG)wxx+(ΘαβGx
+ α 0 ( 2 + α n 0 + α 0 n + 2 β G ) ) w x = 0 + α 0 2 + α n 0 + α 0 n + 2 β G w x = 0 {:+alpha_(0)(2+alphan_(0)+alpha_(0)n+2betaG^(**)))w_(x)^(**)=0\left.+\alpha_{0}\left(2+\alpha n_{0}+\alpha_{0} n+2 \beta G^{*}\right)\right) w_{x}^{*}=0+α0(2+αn0+α0n+2βG))wx=0
Let us apply the above reasoning to the same set of orbits (32) for which we found the potential (33), under the restriction (30) for the four exponents f 1 , f 2 , g 1 , g 2 f 1 , f 2 , g 1 , g 2 f_(1),f_(2),g_(1),g_(2)f_{1}, f_{2}, g_{1}, g_{2}f1,f2,g1,g2. In place of (30) we now put the restriction
( f 1 + f 2 ) g 2 = 1 + g 1 f 1 + f 2 g 2 = 1 + g 1 (f_(1)+f_(2))g_(2)=1+g_(1)\left(f_{1}+f_{2}\right) g_{2}=1+g_{1}(f1+f2)g2=1+g1.
The corresponding potential is
w ( x , z ) = d 3 ( x 2 + z 2 ) f 1 + f 2 + 1 w ( x , z ) = d 3 x 2 + z 2 f 1 + f 2 + 1 w^(**)(x,z)=d_(3)(x^(2)+z^(2))^(f_(1)+f_(2)+1)w^{*}(x, z)=d_{3}\left(x^{2}+z^{2}\right)^{f_{1}+f_{2}+1}w(x,z)=d3(x2+z2)f1+f2+1,
and each orbit is traced with total energy
E = d 3 ( f 1 + f 2 ) 2 ( c 1 0 c 2 0 f 2 + c 1 0 g 1 c 2 0 f 1 ) f 1 + f 2 E = d 3 f 1 + f 2 2 c 1 0 c 2 0 f 2 + c 1 0 g 1 c 2 0 f 1 f 1 + f 2 E=-d_(3)(f_(1)+f_(2))^(2)(c_(1)^(ℓ_(0))c_(2)^(-ℓ_(0)f_(2))+c_(1)^(-ℓ_(0)g_(1))c_(2)^(ℓ_(0)f_(1)))^(f_(1)+f_(2))E=-d_{3}\left(f_{1}+f_{2}\right)^{2}\left(c_{1}^{\ell_{0}} c_{2}^{-\ell_{0} f_{2}}+c_{1}^{-\ell_{0} g_{1}} c_{2}^{\ell_{0} f_{1}}\right)^{f_{1}+f_{2}}E=d3(f1+f2)2(c10c20f2+c10g1c20f1)f1+f2,
where
0 = 2 f 1 f 2 g 1 0 = 2 f 1 f 2 g 1 ℓ_(0)=(2)/(f_(1)-f_(2)g_(1))\ell_{0}=\frac{2}{f_{1}-f_{2} g_{1}}0=2f1f2g1,
and d 3 d 3 d_(3)d_{3}d3 is an arbitrary constant.
Let us check now if the family (32), this time under the restriction (30), which can be described in the presence of the potential w ( y , z ) w ( y , z ) w(y,z)w(y, z)w(y,z) given by (33), is compatible with a potential w ( x , z ) w ( x , z ) w^(**)(x,z)w^{*}(x, z)w(x,z). It can be easily seen that G = 0 G = 0 G^(**)=0G^{*}=0G=0 and the one-variable potential
(51) w ( x ) = d 4 x 2 ( 1 + g 1 + g 2 ) 1 + g 2 (51) w ( x ) = d 4 x 2 1 + g 1 + g 2 1 + g 2 {:(51)w^(**)(x)=d_(4)x^((2(1+g_(1)+g_(2)))/(1+g_(2))):}\begin{equation*} w^{*}(x)=d_{4} x^{\frac{2\left(1+g_{1}+g_{2}\right)}{1+g_{2}}} \tag{51} \end{equation*}(51)w(x)=d4x2(1+g1+g2)1+g2
will be compatible with the given family, traced with the energy
E = d 4 ( g 1 1 + g 2 ) 2 ( f k 0 g 2 g k 0 + f k 0 g k 0 f 2 ) E = d 4 g 1 1 + g 2 2 f k 0 g 2 g k 0 + f k 0 g k 0 f 2 E=-d_(4)((g_(1))/(1+g_(2)))^(2)(f^(k_(0)g_(2))g^(-k_(0))+f^(-k_(0))g^(k_(0)f_(2)))E=-d_{4}\left(\frac{g_{1}}{1+g_{2}}\right)^{2}\left(f^{k_{0} g_{2}} g^{-k_{0}}+f^{-k_{0}} g^{k_{0} f_{2}}\right)E=d4(g11+g2)2(fk0g2gk0+fk0gk0f2),
with k 0 k 0 k_(0)k_{0}k0 given by (35). Examples of one-dimension potentials producing planar families of curves are to be found in Anisiu & Bozis (2007).
With the aid of the criteria (41) and (43) we can decide whether a given family ( α , β α , β alpha,beta\alpha, \betaα,β ) is or is not derived from a 2D potential of the form (10) or (42), respectively. If we want to find all the two-dimension potentials, there remains to us the obligation to check whether the family is derived from a potential of the form
V = w ( x , y ) V = w ( x , y ) V=w^(****)(x,y)V=w^{* *}(x, y)V=w(x,y).
To face the question, we recall that the transformation x x , y z y x x , y z y x rarr x,y rarr z rarr yx \rightarrow x, y \rightarrow z \rightarrow yxx,yzy implies that α β α α β α alpha rarr beta rarr alpha\alpha \rightarrow \beta \rightarrow \alphaαβα and α 0 β 0 α 0 α 0 β 0 α 0 alpha_(0)rarrbeta_(0)rarralpha_(0)\alpha_{0} \rightarrow \beta_{0} \rightarrow \alpha_{0}α0β0α0. So, instead of checking ( α ( x , y , z ) , β ( x , y , z ) ) ( α ( x , y , z ) , β ( x , y , z ) ) (alpha(x,y,z),beta(x,y,z))(\alpha(x, y, z), \beta(x, y, z))(α(x,y,z),β(x,y,z)), we check the family
( β ( x , z , y ) , α ( x , z , y ) ) ( β ( x , z , y ) , α ( x , z , y ) ) (beta(x,z,y),alpha(x,z,y))(\beta(x, z, y), \alpha(x, z, y))(β(x,z,y),α(x,z,y)).
If the conditions (43) are fulfilled for (54), we understand that there exists a potential V = w ( x , y ) V = w ( x , y ) V=w^(****)(x,y)V=w^{* *}(x, y)V=w(x,y) creating the given family. The above trick may be applied also in case that the given family leads to α 0 = 0 , β 0 0 α 0 = 0 , β 0 0 alpha_(0)=0,beta_(0)!=0\alpha_{0}=0, \beta_{0} \neq 0α0=0,β00. If we look for a potential as (53) for the family (32) under the restriction (30), we find again (51).
As for application 2, it can be checked that for f 1 = 1 f 1 = 1 f_(1)=1f_{1}=1f1=1, beside the potential w ( z , y ) w ( z , y ) w(z,y)w(z, y)w(z,y) given by (39), the problem admits of solutions
w ( x , z ) = d 5 ( g 1 z x ) 2 (55) w ( x , y ) = d 6 ( g 2 x g 1 y ) 2 w ( x , z ) = d 5 g 1 z x 2 (55) w ( x , y ) = d 6 g 2 x g 1 y 2 {:[w^(**)(x","z)=(d_(5))/((g_(1)z-x)^(2))],[(55)w^(****)(x","y)=(d_(6))/((g_(2)x-g_(1)y)^(2))]:}\begin{align*} & w^{*}(x, z)=\frac{d_{5}}{\left(g_{1} z-x\right)^{2}} \\ & w^{* *}(x, y)=\frac{d_{6}}{\left(g_{2} x-g_{1} y\right)^{2}} \tag{55} \end{align*}w(x,z)=d5(g1zx)2(55)w(x,y)=d6(g2xg1y)2
the orbits being traced isoenergetically.
Acknowledgements. It is a pleasure for us to thank Prof. F. Puel for his very careful refereeing of the paper and for his valuable suggestions. The work of the first author was financially supported by the scientific program 2CEEX0611-96 of the Romanian Ministry of Education and Research.

References

Anisiu, M.-C.: 2005, Inverse Probl. Sci. Eng. 13, 545
Anisiu, M.-C., Bozis, G.: 2008, Didactica Mathematica 26, 9
Borghero, F., Melis, A.: 1990, CMDA 49, 273
Bozis, G.: 1983, CeMec 31, 129
Bozis, G., Kotoulas, T.A.: 2004, Rendiconti del Seminario della Facolta di Scienze dell' Universita di Cagliari 74, 83
Bozis, G., Kotoulas, T.A.: 2005, Inverse Problems 21, 343
Bozis, G., Nakhla, A.: 1986, CeMec 38, 357
Érdi, B.: 1982, CeMec 28, 209
Kotoulas, T.A., Bozis, G.: 2006, JPhA 39, 9223
Melis, A., Borghero, F.: 1986, Mecc 21, 71
Mertens, R.: 1981, ZAMM 61, 252
Puel, F.: 1992, CMDA 53, 207
Shorokhov, S.G.: 1988, CeMec 44, 193
Szebehely, V.: 1974, in: G. Proverbio (ed.), Proc. of the Internat. Meeting on Earth's Rotation by Satellite Observation, p. 31
Váradi, F., Érdi, B.: 1983, CeMec 32, 395

2009

Related Posts