arxiv:physics/ v1 [physics.ed-ph] 19 Oct 2004

Transcription

1 Orbits in a central force field: Bounded orbits Subhankar Ray Dept of Physics, Jadavpur University, Calcutta 7 3, India J. Shaanna Physics Departent, Visva Bharati University, Santiniketan 7335, India (Dated: August, 3 arxiv:physics/449v [physics.ed-ph] 9 Oct 4 The nature of boundedness of orbits of a particle oving in a central force field is investigated. General conditions for circular orbits and their stability are discussed. In a bounded central field orbit, a particle oves clockwise or anticlockwise, depending on its angular oentu, and at the sae tie oscillates between a iniu and a axiu radial distance, defining an inner and an outer annulus. There are generic orbits suggested in popular texts displaying the general features of a central orbit. In this work it is deonstrated that soe of these orbits, seeingly possible at the first glance, are not copatible with a central force field. For power law forces, the general nature of boundedness and geoetric shape of orbits are investigated. I. INTRODUCTION The central force otion is one of the oldest and widely studied probles in classical echanics. Several failiar force-laws in nature, e.g., Newton s law of gravitation, Coulob s law, van-der Waals force, Yukawa interaction, and Hooke s law are all exaples of central forces. The central force proble gives an opportunity to test one s understanding of the Lagrange s equation, Hailton s equation, Hailton Jacobi ethod, and classical perturbation. It also serves as an introduction to the concept of integrals of otion and conservation laws. We need only to appeal to the principles of conservation of energy and angular oentu to describe the nature and geoetry of the possible trajectories in central force otion. Most books in classical echanics,,3,4,5,6, treatise 7,8, and advanced texts 9, discuss the central force proble. In this article we present soe interesting features of bounded orbits in a central field. A. Kepler s Laws One of the ost rearkable discoveries in the history of physics is that of Keplerian orbits. A treendous wealth of data on planetary positions was collected by Tycho Brahe and Johannes Kepler after detailed observation spread over several decades. After a thorough analysis of this data Johannes Kepler forulated three epirical laws that described and correlated the otion of the five planets then known:. Each planet oves in an elliptical orbit, with the sun at one of its foci.. The radius vector fro the sun to each planet sweeps out equal areas in equal ties. 3. The square of the periods (T of the planets are proportional to the cube of the lengths of the corresponding seiajor axes (a 3. B. Newtonian Synthesis Alost years later Newton realized that the planets go about in their nearly circular orbits around the sun under the influence of the sae force that causes an apple to fall to the ground, i.e., gravitation. Newton s law of gravitation gave a theoretical basis to Kepler s laws. Kepler s laws can be derived fro Newton s law of gravitation; this is often referred to as the Newtonian synthesis. Kepler s first and third laws are valid only in the specific case of inverse square force. There are, however, certain general features which are observed in all central field probles. They include (i certain conserved quantities (energy, and angular oentu, (ii planer nature of orbits, and (iii constancy of areal velocity (Kepler s second law. A large class of central forces allows circular orbits (stable or unstable, bounded orbits, and even closed and periodic orbits. Certain coon characteristics about the generic shapes of bounded orbits can also be ascertained. II. EQUATIONS OF MOTION AND THEIR FIRST INTEGRALS A. Central field orbits: confineent in a plane The central force otion between two bodies about their center of ass can be reduced to an equivalent one body proble in ters of their reduced ass and their relative radial distance r. Hence in this reduced syste, a body having the reduced ass oves about a fixed center of force. Consider the otion of a body under a central force, F = F(r = f(rˆr with the origin as its force center. The potential V (r fro which this force is derived is also a function of r alone, F = V, V V (r. On account of the central nature of the force, the echanical properties of the body do not vary under rotation in any anner around the center of force. Let

2 the body be rotated through an infinitesial angle δθ, where the agnitude δθ is the angle of rotation while the direction is that of the axis of rotation ˆn. The change in the radius vector fro the origin to the body is δr = rsin(θδθ, with δr being perpendicular to r and δθ. Hence δr = δθ r. The change in velocity is siilarly given by δv = δθ v. The Lagrangian of the syste is a function of r and ṙ. The otion is governed by the Lagrange s equation, ( d L L dt ṙ r = If we now require that the Lagrangian of the syste reain invariant under this rotation we obtain δl = L L δr + r ṙ δṙ =. One can define generalized oentu as, p = L ṙ Fro the Lagrange s equation we get, ṗ =. d ( L = L dt ṙ r Replacing L/ ṙ by p and L/ r by ṗ in the equation for δl we get, ṗ δθ r + p δθ ṙ =. δθ d (r p =. dt As δθ is arbitrary we conclude l = r p is a conserved quantity. l is called the angular oentu of the syste. Since l is a constant and is perpendicular to r it follows that the radius vector of the particle lies in a plane perpendicular to l. This iplies that the otion of the particle in a central field is confined to a plane. B. Lagrangian and equations of otion As the otion in a central force field is confined to a plane, it suffices to use plane polar coordinates. One ay write the Lagrangian of the particle as, L = ṙ V (r = (ṙ + r θ V (r. ( We assue the center of force to be at the origin. The coordinates of the body of ass undergoing the central field otion are given by (r, θ. The Lagrange s equation for the θ and r coordinates are given respectively by, d dt d dt ( L ( θ L ṙ L θ = ( L r = (3 d θ r Area = / (r ( r d θ FIG. : Area swept by radius vector C. First integrals and conservation laws r r d θ The canonical oentu corresponding to θ is called the angular oentu (or rather the agnitude of the angular oentu that we discussed before, p θ = L θ = r θ = l As θ is a cyclic coordinate, i.e., the Lagrangian is independent of θ, this angular oentu is conserved. This can be shown fro the Lagrange s equation for θ. d dt ( L θ L θ = p θ = d dt (r θ =. The corresponding integral of otion is, r θ = l (4 and this l can easily be shown to be the agnitude of the angular oentu vector l = r p. This conservation law is essentially equivalent to Kepler s nd Law: The eleentary triangular area swept by the radius vector in an infinitesial tie interval dt is, da = r(rdθ hence it follows fro Eq. (4 that the rate of areal sweep is a constant. da dt = r θ = l (5 There is another first integral of otion associated with the Lagrange s equation for the r coordinate. d dt (ṙ r θ + V r =

3 velocity 4 kinetic energy Angular oentu - total energy potential energy radial position FIG. : Constancy of angular oentu in an inverse square potential. Total energy is 6% of the iniu of effective potential Ṽ (r FIG. 3: Constancy of total energy in an inverse square potential. Total energy is 6% of the iniu of effective potential Ṽ (r. The force in ters of the potential (conservative force is given by, f(r = V/ r, whence the above equation becoes r r θ = f(r Using the first integral of otion, one can convert this second equation into an equation for r alone. r = d ( V + l dr r Integrating we get ṙ + + V (r = E (6 r l where E is a constant of integration, called the energy. This is the law of conservation of total echanical energy. It is interesting to note that for otion in a general force field, l = r θ reains invariant even though r and θ vary with tie (see Fig.. Siilarly E = ṙ / + l /(r +V (r reains a constant even though ṙ and r and hence ṙ / and V (r + l /(r each varies with tie (see Fig. 3. Lagrange s equations are two second order Ordinary Differential Equations (ODE in r and θ. However they decouple, i.e., each equation is expressible in ters of either r or θ. On integrating each equation once we get the first integrals of otion naely, the total echanical energy and the angular oentu. A further integration will yield the coplete solution to the proble. This second integration introduces two ore constants of integration naely, the initial radial (r and angular (θ positions. D. Equation for the orbit Fro the equation giving energy as a first integral of otion, we get the expression for radial velocity, ( ṙ = E V r (7 On integration we get, t = r r dr /(E V /(r. (8 This relation can be inverted to give r as a function of t, r = r(t. The other first integral gives, θ = l r. Substituting the expression for r as a function of t, and integrating, θ θ = l t dt r(t. (9 These two expressions for r(t and θ(t express the equation of the orbit for the particle in a central field in ters of tie t as a paraeter. Instead of expressing the orbit paraetrically in ters of t, one often wants to express the orbit directly as an equation connecting r and θ. Such an equation ay be obtained by eliinating t fro the above expressions for ṙ and θ. dθ dr = l/(r /(E V (r /(r

4 4 On integration this yields, θ θ = r l/(r r /(E V (r /(r dr For understanding the qualitative nature of otion in a central field one looks at the equivalent one diensional proble. ṙ = = (E V (r (E Ṽ (r. r We call Ṽ the effective potential, introduced to ake the proble siilar to that of a particle oving in a one diensional potential field. The effective radial force ( f(r is connected to the effective radial potential (Ṽ (r by the expected relation, Ṽ (r f(r = r ( At a point where the effective potential Ṽ (r equals the energy E, the radial velocity vanishes (ṙ =. In one diensional otion this corresponds to a particle coing oentarily to rest, and having zero kinetic energy. However in the case of central field, the otion is not really one diensional, and even for ṙ =, the particle is not at rest (v = r θˆθ, and it has a non-zero kinetic energy ((/r θ. For the inverse square force, as in the case of gravitation, we have f(r = k/r and V (r = k/r. Effective potential (see Fig. 4 is given by, Ṽ (r = k r + r ( The following properties of the effective potential are easily noted:. Ṽ (r = /r (l / kr, as r the ter within bracket is essentially l /( and hence li r Ṽ (r +. li r Ṽ (r = 3. Ṽ (r = /r(k /r and for large values of r, l /r is negligible copared to k, hence Ṽ (r has a negative value. 4. At r = l /k the function Ṽ (r intersects the r axis, i.e., Ṽ (r =. 5. Ṽ (r reaches a iniu at r = r = l /k. With Ṽ / r r = and Ṽ / r r >. For total energy E = Ṽ (r the particle has zero radial velocity at r = r, and no other radial position is physically accessible, since ṙ becoes iaginary at r r. This corresponds to circular otion ~ V(r l^/(r^ V(r ~ - /r (a circular orbit (b elliptic orbit (c elliptic orbit of higher eccentricity r FIG. 4: Effective potential for an inverse square force For > E > Ṽ (r there is a range of radial positions (r ax r r in for which E Ṽ (r. The particle can ove in this range of r with varying ṙ. The radial velocity ṙ vanishes at the end points r ax and r in where the energy E equals the effective potential Ṽ and ṙ reaches a axiu at r = r The points r ax and r in are called the turning points. The central field particle cannot ove beyond these points, as the energy becoes less than the effective potential, and the expression for radial velocity (ṙ turns iaginary. E < Ṽ (r is a physically ipossible situation, since no (radial position is physically allowed for the particle. For E we have an unbounded otion, where the particle can fly off to infinity. Thus in the case of inverse square force field we can have bounded (E < or unbounded (E otion depending on the energy of the particle. In particular for E = it is parabolic and for E > it is hyperbolic. III. EXISTENCE AND STABILITY OF CIRCULAR ORBITS FOR CENTRAL FORCES At all positions other than where the effective potential is a iniu or a axiu, we have a net effective force f(r = Ṽ / r. When the total energy E is not equal to the iniu of effective potential Ṽ (r, at the points of instantaneous zero radial velocity the particle is pushed away. If the effective potential has a iniu and energy is greater than that iniu then the radial distance has a lower and an upper bound and the particle oves between these radial liits. On the other hand if the effective potential has a axiu, and the energy is less than that axiu, the effective force pushes the (a (b (c

5 (a (b (c (a spiral in (b circle (c spiral out FIG. 5: Spiralling orbits, f(r /r 4 particle away fro the positions of zero radial velocity in an inward or an outward spiral. If the total energy is equal to the axiu or iniu of the effective potential, the syste can stay with zero radial velocity, and hence ove in a circular orbit. The condition for circular orbit is, whence we get, Ṽ r = V r r r f(r = V r r 3 = r = ( r 3 (3 The negative sign on the right hand side clearly shows that the force ust be attractive. The particle oves in a circular orbit, since the force of attraction due to the central field provides the necessary centripetal force. For a given attractive central force ( f(r and V (r given it is possible to have a circular orbit of radius r provided the angular oentu and energy of the particle are given by, l = r 3 f(r (4 E = V (r + r (5 However when the effective potential has a axiu the syste is in an unstable circular orbit. A sall deviation fro this radial position causes the orbit to be unbounded. The effective force that coes into play akes the particle ove away fro the position of zero radial velocity in an inward or outward spiral (Fig. 5. When the effective potential has a iniu, the effective forces cause the particle to reain in a bound orbit confined in an annular space. For sall deviations the annular radii are nearly the sae, and hence the orbit is close to a circle. One can understand the stability question by studying the forces (restoring or unsettling that coe into play when the syste is oved infinitesially fro the position of circular orbit. The circular orbit is stable if, That is equivalent to, Using Eq. (3, IV. Ṽ r = f r r < 3 r 4 Ṽ >. (6 r ( f r + 3 r 4 r > f r r < 3f(r r (7 NATURE OF BOUNDED ORBITS General bounded otion has both lower and upper bounds. It eans that the particle cannot approach nearer than soe iniu or ove farther than soe axiu distance. One has to reeber that the angular velocity has a constant sign, sae as that of the constant angular oentu, throughout the otion. However its agnitude decreases with increase in the radial distance ( /r. Together with the angular otion, the radial distance changes fro r ax to a r in, then back to r ax and so on. Fro this general nature of otion two failies of generic orbits are suggested in popular texts (Fig. 6, Fig. 7. However it can be shown that the generic types shown in the later figure (Fig. 7 are not feasible for any attractive potential. We consider below a general bounded orbit confined in an annular region. We would like to investigate the nature of the orbit close to the point where it touches the inner or the outer annulus. Let us choose the reference line (polar line of the coordinate syste such that the orbit touches the annular ring at θ =, or ore specifically at r = R, θ =. Since at this point P(r = R, θ = the orbit is at its closest (or farthest approach fro the pole, ( r/ t P =. As r is a function of θ, using Eq. (4 we get for otion along the trajectory, hence dr dt = l dr r dθ r ( = dr dθ = (8 θ=

6 FIG. 6: Generic orbits in central field (allowed FIG. 7: Generic orbits in central field (not allowed The equation of otion for r is given by, r r θ = f(r Fro which we get r (θ along the trajectory, l ( d l dr r dθ r r dθ r 3 = f(r d r dθ = r4 f(r l = r + r4 l f(r One can expand r(θ in a Taylor series in θ, and reeber that r ( =, r(θ traj = r( + r (θ + r ( θ! + h.o. = r( + r(4 θ f(r( l! + h.o. = r( + (r( + r(4 θ l f(r(! + h.o. Consider a tangent to the annulus at the point P, Fig.

7 7 Inner annulus r( θ tan θ r( θ traj trajectory (not allowed Outer annulus r( θ tan r( θ θ traj trajectory (not allowed rannul r annul trajectory (allowed trajectory (allowed tangent trajectory (not allowed tangent FIG. 8: Curvature of the trajectory in coparison to the tangent to the inner annulus FIG. 9: Curvature of the trajectory in coparison to the outer annulus 8, r( r(θ = cos(θ which can be expanded to give r(θ tan along the tangent, r(θ tan = r( + r( θ + h.o. Notice that the constant or the θ independent ters in r(θ traj and r(θ tan are equal, and the first leading power of θ is θ in both the cases. The coefficient of θ (curvature will deterine the nature of the trajectory in reference to the tangent. We first study the curvature of the orbit in reference to the tangent. The force being attractive f(r( <, for sall angular distance ( (θ away fro the point P, r( θ traj < r( θ tan (9 which eans that the orbit bends ore sharply than the tangent and stays closer to the inner annulus for nearby points (Fig. 8. Hence the second type of orbits shown in soe texts (Fig. 7 are not possible for central forces. For the outer annulus the analysis with respect to tangent is not very eaningful as one can set up an even stronger bound for its orbit, directly considering the potential Ṽ (r. In this case the effective potential has a positive slope with respect to the radius, and hence a negative effective force ( f(r <. The orbit should not only reain nearer than the tangent, but even nearer than the outer annulus. This condition is confired if we study the orbit near the point P where it touches the outer annulus. hence r( θ traj = r( + r(4 ( θ f(r l + h.o.! r( θ traj < r( = r ax ( This confirs that the outer annulus is indeed the outer bound of the trajectory. The nature of the orbit near the point of contact P at the outer annulus is shown in figure 9. V. BOUNDED ORBITS FOR THE POWER LAW CENTRAL FORCE A. Existence of stable circular orbit For the case of a power law central potential V (r = k r n, Fro the stability condition Eq. (7, we find, (n + nk r n+ f r r < 3f(r r ( < 3 nk f(r = r n+ ( k r n+ r n < (

8 8 4 For the point of extreu 3 l^/( r^ r n = an l. Finding the second derivative of Ṽ with respect to r at the point r, ~ V(r (a unstable circular orbit (b unbounded orbit d Ṽ an(n + dr = r n+ + 3 r 4 = ( an(n + r 4 r n + 3 = 3 ( (n r V(r ~ - /r^3 = r 4 l ( n < FIG. : Effective potential for f(r /r 4 force Hence the circular orbit is stable for an attractive power law potential that varies slower than inverse square (or a force that varies slower than inverse cube. Hence it is a point of axiu. For E > Ṽ (r we always have an unbounded orbit. For E < Ṽ (r the orbit is seibounded, bounded above or bounded below, according to its initial state. The particle either spirals in or spirals out. For E = Ṽ (r we get an unstable circular orbit.. V (r = a/r n where n = B. Study of boundedness of orbits. V (r = a/r n where n > Consider the case when n >. We have the effective potential, Ṽ (r = a r n + r (3. Ṽ (r = /r n (a l r n /(, and as r we ay neglect l r n /( in coparison to a, thus li r Ṽ (r =.. li r Ṽ (r =. 3. For large but finite r, Ṽ (r = /r (a/r n l /(, and a/r n is negligible copared to l /(, and hence Ṽ (r is positive. 4. Ṽ (r intersects the r axis at r where (r n = a/l. 5. Ṽ (r has a axiu at r where r n = an/l. For a axiu point of Ṽ we should have Ṽ / r r = and Ṽ / r r <. Ṽ r = an r n+ r 3 =. The effective radial potential is Ṽ (r = a/r + l /(r. This is essentially an attractive or repulsive inverse square ter. Ṽ (r = ã r, ã > if a > ã < if a < This potential cannot give circular orbit ever. If the effective potential is attractive and E < it has an upper bound of radial distance. A typical orbit would therefore be an inward spiral. For a repulsive effective potential we will have outward spiral oving to infinite radial distance. 3. V (r = a/r n where > n > In this case the effective potential is, Ṽ (r = a r n + r (4 We have the following properties of Ṽ (r. Ṽ (r = /r (l / a/r n and as r we ay neglect a/r n in coparison to l /, and thus li r Ṽ (r +.. li r Ṽ (r =.

9 9 3 3 l^/( r^ l^/( r^ ~ V(r - ~ V(r V(r ~ ln(r -3 V(r ~ -/r FIG. : Effective potential for f(r /r force FIG. : Effective potential for f(r /r force 3. Ṽ (r = /rn (a l /r n and for large but finite values of r, l /r n is negligible copared to a, hence Ṽ (r is negative. 4. Ṽ (r intersects the r axis at r where (r n = l /(a. 5. Ṽ (r reaches a iniu at r where r n = l /an. For a iniu point of Ṽ we should have Ṽ / r r = and Ṽ / r r >. Ṽ r = an r n+ r 3 =. For the point of extreu r n = an. Finding the second derivative of Ṽ with respect to r at the point r, d Ṽ an(n + dr = r n+ + 3 r 4 = ( an(n + r 4 r n + 3 = 3 ( (n r = l r 4 ( n > Hence it is a point of iniu. For this case the orbit can be bounded or unbounded depending on the total energy E of the particle. 4. V (r = a ln r This corresponds to one over r force f(r a/r. The effective potential is given by, Ṽ (r = a lnr + We have the following properties of Ṽ (r. li r Ṽ (r +.. li r Ṽ (r +. r (5 3. There is no point of intersection with the r axis, and Ṽ (r is always positive. 4. Ṽ (r reaches a iniu at r = l/ a. For a iniu point of Ṽ we should have Ṽ / r r = and Ṽ / r r >. Ṽ r = a r for the point of extreu r = a r 3 The second derivative of Ṽ at r = r, d Ṽ dr = a r + 3 r 4 = ( r a + 3 r = r a >

10 3 The second derivative of Ṽ with respect to r, 5 d Ṽ dr = an(n r n + 3 r 4 = ( r 4 an(n r n = l r 4 ( + n > 5 ~ V(r V(r ~ r^ l^/( r^ Hence it is a point of iniu. The orbit is always bounded. The findings for the general power law potential can be suarized as follows, V (r = sign(nar n n (7 f(r = abs(nar n (8 and V (r = b lnr (9 FIG. 3: Effective potential for f(r r force f(r = b r (3 Hence it is a point of iniu. The orbit is always bounded. TABLE I: Dependence of boundedness on power law. n of V (r ar n nature of boundedness, ( always unbounded 5. V (r = ar n where n > For this potential the effective potential is given by, spiralling orbit,( bounded or unbounded depending on E (lnr always bounded Ṽ (r = ar n + We have the following properties of Ṽ (r r n > (6,+ ( always bounded. li r Ṽ (r +.. li r Ṽ (r There is no point of intersection with the r axis, and Ṽ (r is always positive. 4. Ṽ (r has a iniu at r, r n+ = l /an. For a iniu point of Ṽ we should have Ṽ / r r = and Ṽ / r r >. Ṽ r = anrn for the point of extreu r n+ = an r 3 C. Stable bounded orbit, geoetric shape The stable bounded orbits in a power law central field can have the following fors. V (r = a/r n, < n. V (r = b log r 3. V (r = a r n, n > In all the above cases the derivative of V (r with respect to r is always positive, and hence the force is necessarily attractive (f(r = V (r/ r <. In the first case the effective one diensional potential Ṽ (r goes to infinity as r. Ṽ (r has a iniu at soe r = r, and it has a negative slope for all r < r, and positive slope for all r > r. For E = Ṽ (r we get stable circular orbit, and for Ṽ (r < E <. the orbits are still bounded. For cases ( and (3 the circular orbit is stable, and orbits are always bounded for any energy.

11 TABLE II: Dependence of stability of circular orbits on power law. n of V (r ar n stability of circular orbits, unstable, ( stable VI. CONCLUSION Existence of bounded orbit for a large class of attractive central field has been discussed. The generic nature of central field bounded orbits is analytically derived. Certain class of these orbits (Fig. 7 presented in popular texts, are shown to be non-feasible. Sall deviation fro circularity in the case of central field is often expressed in ters of inverse of radial distance (u = /r. u = u + a cos(βθ r r = + a r cos(βθ where r = /u. It is interesting to note that these orbits are soetie istakenly identified with diagras of the for shown in Fig. 7. The above expression in fact corresponds to a class of orbits that look generically like those shown in Fig. 6. The generic features of central force orbits discussed here have been verified by coputer siulation for a large class of central force fields. The figures shown here (Fig., 3, 5 and 6 were generated by these siulations. Students interested in studying and generating such orbits will find Ref. helpful. Acknowledgeent The authors wish to express their indebtedness to the well known texts by Goldstein, Landau and Arnold 9. They also acknowledge their teachers in related graduate courses at Stony Brook, Prof. Max Dresden, Prof. A. S. Goldhaber and Prof. Leon A. Takhtajan. Authors gratefully acknowledge the encourageent received fro Prof. Shyaal SenGupta of Presidency College, Calcutta. The aterial presented here was used in a graduate level classical echanics course at Jadavpur University during SR wishes to thank his students for stiulating discussions. Electronic address: subho@juphys.ernet.in H. Goldstein, Classical Mechanics, Addison-Wesley, 98. L.D. Landau, E.M. Lifshitz, Mechanics, Pergaon Press, K.R. Syon, Mechanics, Addison-Wesley, D.T. Greenwood, Classical Dynaics, Prentice Hall, J.L. Synge, and B.A. Griffith, Principles of Mechanics, Mc- GrawHill, A. Soerfeld, Mechanics, Acadeic Press, E.T. Whittaker, A Treatise on the Analytical Dynaics of Particles and Rigid Bodies, Dover, L.A. Pars, Introduction to Dynaics, Cabridge, V.I. Arnold, Matheatical Methods of Classical Mechanics, Springer Verlag, 989. R. Abraha, and J.E. Marsden, Foundations of Mechanics, Benjain-Cuings Publ.Co.Inc., 978. H. Goldstein, Fig. 3-7, Fig. 3-3 and Eq in Classical Mechanics, Addison-Wesley, 98. H. Gould, J. Tobochnik, An Introduction to Coputer Siulation Methods, Addison-Wesley, 988.