Abstract
We discuss a method for obtaining vector constants of motion for time-dependent central fields.
Vector Constants of Motion for Time-Dependent Kepler and Isotropic Harmonic Oscillator Potentials
O.M. Ritter* * e-mail: fsc1omr@fisica.ufsc.br ;
Departamento de Física,
Universidade Federal de Santa Catarina, Trindade
88040-900 Florianópolis SC, Brasil
F.C. Santos * e-mail: fsc1omr@fisica.ufsc.br , and A. C. Tort * e-mail: fsc1omr@fisica.ufsc.br
Instituto de Física, Universidade Federal do Rio de Janeiro
Cidade Universitária, Caixa Postal 68528
21945-970 Rio de Janeiro RJ, Brasil
Received on 8 June, 2000. Revised version received on 21 March, 2001.
We discuss a method for obtaining vector constants of motion for time-dependent central fields.
I Introduction
The search for first integrals has assumed an increasing importance in the determination of the integrability of a dynamical system. The notion of integrability is related to the existence of first integrals of motion. Several methods of finding first integrals are available in the literature. Some of them are based on the symmetries of the system [1], [2]. Another possibility relies on a direct search for invariants with a specified form (usually, a polynomial form on the velocities) [3], [4]. Still other methods to determine the integrability of a system are available as, for instance, the Painlevé Analysis [5]. One of the most aesthetically appealing and important problem in physics is the central field problem [6]. Energy and angular momentum associated with this type of field are well known conserved quantities. However, other vector and tensor conserved quantities have been associated with some particular central fields. The Laplace-Runge-Lenz vector [7] is a vector first integral of motion for the Kepler problem; the Fradkin tensor [8] is conserved for the case of the harmonic oscillator and for any central field it is possible to find a vector first integral of motion as was shown in [9]. In the general case these additional integrals of motion turn out to be complicated functions of the position r and linear momentum p of the particle probing the central field. When orbits are closed and degenerated with respect to the mechanical energy, however, we should expect these additional constants of motion to be simple functions of r and p. In this article we wish to exploit further this line of reasoning by determining the existence of such additional vector first integrals of motion for some special class of time-dependent Kepler and isotropic harmonic oscillator problems. In particular, we will show that for the time-dependent Kepler problem the existence of a vector constant of motion coupled to a simple transformation of variables turns the problem easily integrable.
II Constructing vector constants of motion
The force f(r, t) acting on a particle moving in a central but otherwise arbitrary and possible time-dependent field of force g(r, t) can be written as
where r = r(t) is the position vector with respect to the center of force, r is its magnitude, and t is the time. Assuming a constant vector j of the form
where p = p(t) = m (t) is the linear momentum, m is the reduced mass and A, B are arbitrary scalar functions of p, r and t, it follows that the functions A and B must satisfy
Eliminating B between (3) and (4) we obtain
It follows from (4) that j can be written in the form
If the field g(r, t) is known any solution of (5) will yield a vector constant of motion of the form given by (6). Equation (5), however, is a differential equation whose solution may turn out to be a hard task to accomplish. Nevertheless, we can make progress if instead of trying to tackle it directly we make plausible guesses concerning A thereby linking j to specific forms of the field g(r, t).
III Time-dependent fields
Let us consider time-dependent central force fields for which we can build more general vector first integrals of motion. As with the time-independent case there are of course several possibilities when it comes to the choice of a function A for a time-dependent central field. Here is one
Using (7) in (5) and assuming the additional condition
we arrive at
The vector first integral of motion associated with (9) is
where we have made use of the fact that the angular momentum is constant for any arbitrary central field whether it is time-independent or not. If in (10) we set f = 1 and C ¹ 0, we can recover the time-independent Kepler problem. If we set (m/f) d2f/dt2 = -k(t), that is, if f is an arbitrary function of the time, and also C = 0, we have the time-dependent isotropic harmonic oscillator field, F(r) = -k(t) r. In this case j is equal to the first term on the R.H.S. of (10). If (m/f)d2f/dt2 = -k and C = 0, we have the time-independent isotropic harmonic oscillator field. Finally, if (m/f)d2f/dt2 = 0, f ¹ 0 for any instant of time t and C ¹ 0, we have a time-dependent Kepler field. In this case
and the associated constant vector is a time-dependent generalization of the Laplace-Runge-Lenz vector
As another example let us consider the time-dependent isotropic harmonic oscillator and show how it is possible to generalize the Fradkin tensor for this case. Let the function A(r, p, t) be written as
Using (5) and assuming the additional condition
we obtain
We can solve (15) thoroughly if g(r, t) is a function of the time t only. In this case, as before, we end up with the time-dependent isotropic harmonic oscillator. From (13) and (14) we can derive the vector j associated with (13) which, in terms of components reads
where Fij is a generalized Fradkin tensor and can be written as
It is not hard to see that the trace of this generalized Fradkin tensor F(r, p, t) becomes the energy of the particle when g(r, t) is a constant field.
IV Obtaining explicit solutions: An alternative way
Now we will show how to take advantage of the vector constant j to obtain the solution for the Kepler and the isotropic harmonic oscillator potentials. But firstly we must establish some very general relationships. First notice that (6) can be recasted into the form
Integrating (18) we readily obtain
Let us now show how we can obtain the orbits in the case of the time-dependent Kepler problem. In what follows, we will be dealing with a special class of time-dependent Kepler problems, namely, those satisfying (9). In polar coordinates the angular momentum conservation law is written in the form
and this allows to rewrite (18) as
where we have introduced the polar unitary vectors and . The component of the vector j in the direction of is given by
In section 3 we determined a generalized Laplace-Runge-Lenz vector for the time-dependent Kepler problem. The scalar function A(r, p, t) associated with this vector was found to be
Considering A as a function of q we can write
Equations (22) and (24) lead to
where a is the angle between j and the axis. Upon integrating (25) we find
This is the orbit equation we were looking for.
It is possible to transform the time-dependent Kepler problem into the time-independent one. In order to see this we define a new position vector r¢ according to
and redefine our time parameter according to
We can recast the equation of motion for the time-dependent Kepler problem, namely
into a simpler form. Making use of the above transformations the equation of motion (29) can be written as
Equation (30) corresponds to the usual time-independent Kepler problem. Equations (27) and (28) show that the open solutions of (30) are transformed into the open solutions of (29) with the same angular size and that closed solutions of (30) are associated with spiraling solutions of (29). To conclude consider the total mechanical energy associated with (30)
Since p¢ and p are related by
and r¢ and r by (27) we easily obtain
which is a conserved quantity and can be interpreted as a generalization of the energy of the particle under the action of a time-dependent Kepler field.
V Conclusions
In this paper we have outlined a simple and effective method for treating problems related with time-dependent and time independent central force fields. In particular we have dealt with the Kepler problem and the isotropic harmonic oscillator fields. The difficulty in finding vector constants of motion for central fields stems from the fact that, in general, orbits for this type of problem are not closed, therefore any new ways to attack those problems are welcome. In our method this difficulty is transferred, so to speak, to the obtention for each possible central field, which can be time-dependent or not, of a certain scalar function of the position, linear momentum and time. For a given central field this scalar function is a solution of (5). In the general case, the obtention of the scalar function is a difficult task. Judicious guesses, however, facilitate the search for solution of (5) and this is what we have done here. Notice also that the method apllies to closed orbits f0 < 1 as well as open orbits f0 > 1. The case f0=1, the parabolic orbit is also contemplated.
References
[1] P.J. Olver, 1986 Applications of Lie groups to differential equations (New York: Springer Verlag).
[2] F. Cantrijn and W. Sarlet, Siam Review 23 467 (1981).
[3] E.T. Whittaker, A Treatise on the Analytical Dynamics of Particles and Rigid Bodies 4th ed. (Cambridge: Cambridge University Press), 1937.
[4] J. Hietarinta, Physics Reports 147, 87 (1987).
[5] M.J. Ablowitz, A. Ramani, and H. Segur Lett. Nuovo Cimento 23, 333 (1978).
[6] H. Goldstein Classical Mechanics 2nd ed (Reading MA: Addison-Wesley), 1980.
[7] Laplace P S 1799 Traité de Mécanique Celeste; Runge C Vektoranalysis V1 pg. 70 (Leipzig: S. Hirzel, Leipzig), 1919; W Pauli Z. Phyisik 36, 336-63 (1926); Lenz W Z. Phyisik 24, 197 (1926); Heintz W H Am. J. Phys. 42, 1078 (1974). See also Goldstein H Am. J. Phys. 43, 737 (1975); ibid. 44, 1123 (1976).
[8] Fradkin D M Am. J. Phys. 33, 207 (1965); Prog. Theor. Phys. 37, 798 (1967).
[9] C.C. Yan, J. Phys. A 24, 4731 (1991).
e-mail: filadelf@if.ufrj.br
e-mail:tort@if.ufrj.br
- [1] P.J. Olver, 1986 Applications of Lie groups to differential equations (New York: Springer Verlag).
- [2] F. Cantrijn and W. Sarlet, Siam Review 23 467 (1981).
- [3] E.T. Whittaker, A Treatise on the Analytical Dynamics of Particles and Rigid Bodies 4th ed. (Cambridge: Cambridge University Press), 1937.
- [4] J. Hietarinta, Physics Reports 147, 87 (1987).
- [5] M.J. Ablowitz, A. Ramani, and H. Segur Lett. Nuovo Cimento 23, 333 (1978).
- [6] H. Goldstein Classical Mechanics 2nd ed (Reading MA: Addison-Wesley), 1980.
- [7] Laplace P S 1799 Traité de Mécanique Celeste; Runge C Vektoranalysis V1 pg. 70 (Leipzig: S. Hirzel, Leipzig), 1919;
- W Pauli Z. Phyisik 36, 336-63 (1926);
- Lenz W Z. Phyisik 24, 197 (1926);
- Heintz W H Am. J. Phys. 42, 1078 (1974).
- See also Goldstein H Am. J. Phys. 43, 737 (1975); ibid 44, 1123 (1976).
- [8] Fradkin D M Am. J. Phys. 33, 207 (1965);
- Prog. Theor. Phys. 37, 798 (1967).
- [9] C.C. Yan, J. Phys. A 24, 4731 (1991).
Publication Dates
-
Publication in this collection
05 Mar 2002 -
Date of issue
June 2001
History
-
Reviewed
21 Mar 2001 -
Received
08 June 2000