Arxiv:Physics/0410149V1 [Physics.Ed-Ph] 19 Oct 2004 Ftepsil Rjcoisi Eta Oc Oin Most Motion
![Arxiv:Physics/0410149V1 [Physics.Ed-Ph] 19 Oct 2004 Ftepsil Rjcoisi Eta Oc Oin Most Motion](http://data.docslib.org/img/5092806/arxiv-physics-0410149v1-physics-ed-ph-19-oct-2004-ftepsil-rjcoisi-eta-oc-oin-most-motion.jpg)
arXiv:physics/0410149v1 [physics.ed-ph] 19 Oct 2004 ftepsil rjcoisi eta oc oin Most motion. mechanics force classical central in in books trajectories geometry and possible and nature the energy the of of describe conservation to of momentum principles only angular the need We to of appeal laws. concept conservation to the and to motion introduction an of as integrals serves pertur- also classical It Hamilton’s and bation. equation, method, Jacobi one’s Lagrange’s Hamilton test equation, The the to opportunity forces. of an central understanding gives of problem examples force all central are gravitation, law of Hooke’s interaction, Yukawa and law force, Waals Newton’s van-der law, e.g., Coulomb’s nature, familiar Several in mechanics. force-laws classical in problems studied ril epeetsm neetn etrso bounded field. of central features a interesting in some orbits present we article texts vanced asta ecie n orltdtemto ftefive the of motion known: the then correlated empirical planets and three analysis described formulated thorough that Kepler a laws Johannes After observation data this detailed Tycho decades. of by after several collected over Kepler was spread Johannes positions and planetary wealth Brahe tremendous on A data orbits. of Keplerian of that is physics h eta oc oini n fteods n widely and oldest the of one is motion force central The n ftems eakbedsoeisi h itr of history the in discoveries remarkable most the of One .Tesur ftepros( periods planet the of each square to The 3. sun the from vector radius The the 2. with orbit, elliptical an in moves planet Each 1. wesoteulaesi qa times. equal in areas equal out sweeps foci. its of one at sun rprinlt h ueo h egh ftecor- the of ( lengths axes the semimajor of responding cube the to proportional 9,10 fbuddesadgoercsaeo risaeinvestigat Fo are orbits field. of t force of shape central some geometric a that and with demonstrated boundedness compatible is of not it are work glance, this popula first In in d suggested orbit. radial orbits central maximum generic a a are and There minimum annulus. a outer depend between anticlockwise, oscillates or stability time clockwise their same moves and particle orbits a circular orbit, for conditions General h aueo onens fobt fapril oigi c a in moving particle a of orbits of boundedness of nature The .INTRODUCTION I. ics h eta oc rbe.I this In problem. force central the discuss .Kpe’ Laws Kepler’s A. hsc eatet iv hrt nvriy Santiniket University, Bharati Visva Department, Physics risi eta oc ed one orbits Bounded field: force central a in Orbits 1,2,3,4,5,6 eto hsc,Jdvu nvriy acta7002 Ind 032, 700 Calcutta University, Jadavpur Physics, of Dept T a 2 3 treatise , ftepaesare planets the of ) ). 7,8 Dtd uut1 2003) 1, August (Dated: n ad- and , uhna Ray Subhankar .Shamanna J. hi etro ascnb eue oa qiaetone mass distance equivalent reduced radial an their relative of to terms reduced in be problem can body mass of center their eitc bu h eei hpso one riscan orbits bounded ascertained. of be shapes also generic charac- the common Certain about and teristics orbits. allows orbits, periodic forces bounded and central unstable), closed of even or class (stable large orbits A (Ke- circular velocity law). areal of second constancy pler’s (iii) planer conserved and (ii) certain orbits, momentum), of cen- (i) nature angular all include and They in (energy, observed quantities problems. are field however, which are, tral the features There in general force. syn- only square valid Newtonian certain inverse are the of laws case as third specific and to first referred grav- Kepler’s of often thesis. law is laws. Newton’s this Kepler’s from derived to itation; law be basis Newton’s can theoretical laws gravitation. Kepler’s a ap- i.e., gave an ground, causes gravitation the that of to force fall same to sun the ple the of around influence orbits circular the nearly under their in about go ets oyhvn h eue asmvsaotafixed a F about moves mass reduced force. the of center having body a potential ucinof function aini n anraon h etro oc.Let force. of ro- center under the vary around not manner do any body in the tation of properties chanical h eta oc oinbtentobde about bodies two between motion force central The lot10yasltrNwo elzdta h plan- the that realized Newton later years 100 Almost osdrtemto fabd ne eta force, central a under body a of = motion the Consider nacuto h eta aueo h oc,teme- the force, the of nature central the of account On .Cnrlfil ris ofieeti plane a in confinement orbits: field Central A. I QAIN FMTO N THEIR AND MOTION OF EQUATIONS II. F r icse.I one eta field central bounded a In discussed. are ∗ n nisaglrmmnu,ada the at and momentum, angular its on ing ( et ipaigtegnrlfaue of features general the displaying texts r r ed. = ) oe a ocs h eea nature general the forces, law power r eeobt,seigypsil tthe at possible seemingly orbits, hese V ( r f r sac,dfiiga ne n an and inner an defining istance, alone, rmwihti oc sdrvdi loa also is derived is force this which from ) ( nrlfrefil sinvestigated. is field force entral r .NwoinSynthesis Newtonian B. ) ˆ r IS INTEGRALS FIRST n713,India 731235, an ihteoii sisfrecne.The center. force its as origin the with F = r ec nti eue system, reduced this in Hence . − ia ∇ ,V V, ≡ V ( r ). m n their and 2 the body be rotated through an infinitesimal angle δθ, where the magnitude δθ is the angle of rotation while the direction is that of the axis of rotation ˆn. The change in the radius vector from the origin to the body is δr = rsin(θ)δθ, with δr being perpendicular to r r θ and| δθ.| Hence δr = δθ r. The change in velocity is r d similarly given by δv = δ×θ v. The Lagrangian of the system is a function of r and×r˙. The motion is governed θ by the Lagrange’s equation, Area = 1/2 (r) ( r d ) d ∂L ∂L =0 d θ dt ∂r˙ − ∂r If we now require that the Lagrangian of the system re- main invariant under this rotation we obtain r ∂L ∂L δL = δr + δr˙ =0 . FIG. 1: Area swept by radius vector ∂r · ∂r˙ · One can define generalized momentum as, C. First integrals and conservation laws ∂L p = ∂r˙ The canonical momentum corresponding to θ is called From the Lagrange’s equation we get, the angular momentum (or rather the magnitude of the angular momentum that we discussed before), . d ∂L ∂L p˙ = = dt ∂r˙ ∂r ∂L p = = mr2θ˙ = l θ ˙ Replacing ∂L/∂r˙ by p and ∂L/∂r by p˙ in the equation ∂θ for δL we get, As θ is a cyclic coordinate, i.e., the Lagrangian is inde- p˙ δθ r + p δθ r˙ =0 . pendent of θ, this angular momentum is conserved. This · × · × can be shown from the Lagrange’s equation for θ. d δθ (r p)=0 . d ∂L ∂L · dt × = 0 dt ∂θ˙ − ∂θ As δθ is arbitrary we conclude l = r p is a conserved × d 2 quantity. l is called the angular momentum of the sys- p˙θ = (mr θ˙) = 0. tem. Since l is a constant and is perpendicular to r it dt follows that the radius vector of the particle lies in a The corresponding integral of motion is, plane perpendicular to l. This implies that the motion of the particle in a central field is confined to a plane. mr2θ˙ = l (4)
and this l can easily be shown to be the magnitude of the B. Lagrangian and equations of motion angular momentum vector l = r p. This conservation law is essentially× equivalent to Ke- As the motion in a central force field is confined to a pler’s 2nd Law: plane, it suffices to use plane polar coordinates. One may The elementary triangular area swept by the radius write the Lagrangian of the particle as, vector in an infinitesimal time interval dt is,
1 2 1 2 2 2 1 L = mr˙ V (r)= m(r ˙ + r θ˙ ) V (r) . (1) dA = r(rdθ) 2 − 2 − 2 We assume the center of force to be at the origin. The hence it follows from Eq. (4) that the rate of areal sweep coordinates of the body of mass m undergoing the central is a constant. field motion are given by (r, θ). The Lagrange’s equation for the θ and r coordinates dA 1 l = r2θ˙ = (5) are given respectively by, dt 2 2m d ∂L ∂L = 0 (2) There is another first integral of motion associated dt ˙ − ∂θ with the Lagrange’s equation for the r coordinate. ∂θ d ∂L ∂L = 0 (3) d ∂V dt ∂r˙ − ∂r (mr˙) mrθ˙2 + =0 dt − ∂r 3
11 60
10 40 9 kinetic energy velocity 8 20
7 Angular momentum 0 6 total energy 5 -20
4
3 -40 potential energy
2 -60 radial position 1
0 -80 0 1 2 3 4 5 6 0 1 2 3 4 5 6
FIG. 2: Constancy of angular momentum in an inverse square FIG. 3: Constancy of total energy in an inverse square poten- potential. Total energy is 60% of the minimum of effective tial. Total energy is 60% of the minimum of effective potential potential V˜ (r). V˜ (r).
The force in terms of the potential (conservative force) D. Equation for the orbit is given by, f(r)= ∂V/∂r, whence the above equation − becomes From the equation giving energy as a first integral of motion, we get the expression for radial velocity, mr¨ mrθ˙2 = f(r) − 2 l2 Using the first integral of motion, one can convert this r˙ = E V (7) m − − 2mr2 second equation into an equation for r alone. s d 1 l2 On integration we get, mr¨ = V + −dr 2 mr2 r dr t = . (8) 2 2 Integrating we get r0 2/m(E V l /(2mr )) Z − − 2 This relation can bep inverted to give r as a function of t, 1 2 1 l mr˙ + + V (r)= E (6) r = r(t). The other first integral gives, 2 2 mr2 l where E is a constant of integration, called the energy. θ˙ = . This is the law of conservation of total mechanical energy. mr2 It is interesting to note that for motion in a general Substituting the expression for r as a function of t, and force field, l = mr2θ˙ remains invariant even though r and integrating, ˙ 2 θ vary with time (see Fig. 2). Similarly E = mr˙ /2+ t 2 2 l dt l /(2mr )+V (r) remains a constant even thoughr ˙ and r θ θ = . (9) 2 2 2 0 2 and hence mr˙ /2 and V (r)+ l /(2mr ) each varies with − m 0 r(t) Z time (see Fig. 3). These two expressions for r(t) and θ(t) express the equa- Lagrange’s equations are two second order Ordinary tion of the orbit for the particle in a central field in terms Differential Equations (ODE) in r and θ. However they of time t as a parameter. decouple, i.e., each equation is expressible in terms of ei- Instead of expressing the orbit parametrically in terms ther r or θ. On integrating each equation once we get of t, one often wants to express the orbit directly as an the first integrals of motion namely, the total mechanical equation connecting r and θ. Such an equation may be energy and the angular momentum. A further integra- obtained by eliminating t from the above expressions for tion will yield the complete solution to the problem. This r˙ and θ˙. second integration introduces two more constants of in- dθ l/(mr2) tegration namely, the initial radial (r0) and angular (θ0) = positions. dr 2/m(E V (r) l2/(2mr2)) − − p 4
On integration this yields, 30
r 2 l^2/(2mr^2) (a) circular orbit l/(mr ) (b) elliptic orbit θ θ0 = dr 20 2 2 (c) elliptic orbit of higher eccentricity − r0 2/m(E V (r) l /(2mr )) Z − −
For understandingp the qualitative nature of motion in 10 a central field one looks at the equivalent one dimensional ~ problem. V(r) 0 2 l2 r˙ = (E V (r) ) m − − 2mr2 r -10 (c) 2 = (E V˜ (r)). (b) m − r (a) -20 We call V˜ the effective potential, introduced to make the problem similar to that of a particle moving in a -30 one dimensional potential field. The effective radial force V(r) ~ - 1/r (f˜(r)) is connected to the effective radial potential (V˜ (r)) by the expected relation, -40 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 ∂V˜ (r) r f˜(r)= (10) − ∂r FIG. 4: Effective potential for an inverse square force At a point where the effective potential V˜ (r) equals the energy E, the radial velocity vanishes (r ˙ = 0). In one For 0 >E> V˜ (r ) there is a range of radial positions dimensional motion this corresponds to a particle com- 0 (r r r ) for which E V˜ (r). The particle can ing momentarily to rest, and having zero kinetic energy. max min move in≥ this≥ range of r with varying≥ r ˙. The radial velocity However in the case of central field, the motion is not r˙ vanishes at the end points r and r where the really one dimensional, and even forr ˙ = 0, the particle max min energy E equals the effective potential V˜ andr ˙ reaches a is not at rest (v = rθ˙θˆ), and it has a non-zero kinetic maximum at r = r The points r and r are called 2 ˙2 0 max min energy ((1/2m)r θ ). the turning points. The central field particle cannot move For the inverse square force, as in the case of gravita- beyond these points, as the energy becomes less than the tion, we have f(r)= k/r2 and V (r)= k/r. Effective − − effective potential, and the expression for radial velocity potential (see Fig. 4) is given by, (r ˙) turns imaginary. ˜ k l2 E < V (r0) is a physically impossible situation, since V˜ (r)= + (11) no (radial) position is physically allowed for the particle. − r 2mr2 For E 0 we have an unbounded motion, where the The following properties of the effective potential are particle can≥ fly off to infinity. Thus in the case of inverse easily noted: square force field we can have bounded (E < 0) or un- 2 2 bounded (E 0) motion depending on the energy of the 1. V˜ (r) = 1/r (l /2m kr), as r 0 the term ≥ within bracket is essentially− l2/(2→m) and hence particle. In particular for E = 0 it is parabolic and for E > 0 it is hyperbolic. limr→0 V˜ (r) + → ∞ 2. limr→∞ V˜ (r)=0 III. EXISTENCE AND STABILITY OF 3. V˜ (r) = 1/r(k l2/2mr2) and for large values of CIRCULAR ORBITS FOR CENTRAL FORCES − − r, l2/2mr2 is negligible compared to k, hence V˜ (r) has a negative value. At all positions other than where the effective potential ∗ is a minimum or a maximum, we have a net effective force 4. At r = l2/2mk the function V˜ (r) intersects the r ∗ f˜(r) = ∂V˜ /∂r = 0. When the total energy E is not axis, i.e., V˜ (r )=0. − 6 equal to the minimum of effective potential V˜ (r), at the ∗ 2 5. V˜ (r) reaches a minimum at r0 = 2 r = l /mk. points of instantaneous zero radial velocity the particle ˜ 2 ˜ 2 · is pushed away. If the effective potential has a minimum With ∂V /∂r r0 = 0 and ∂ V /∂ r r0 > 0. | | and energy is greater than that minimum then the radial For total energy E = V˜ (r0) the particle has zero radial distance has a lower and an upper bound and the particle velocity at r = r0, and no other radial position is phys- moves between these radial limits. On the other hand if ically accessible, sincer ˙ becomes imaginary at r = r0. the effective potential has a maximum, and the energy is This corresponds to circular motion. 6 less than that maximum, the effective force pushes the 5
2 the effective potential has a minimum, the effective forces (a) spiral in (b) circle cause the particle to remain in a bound orbit confined (c) spiral out 1.5 in an annular space. For small deviations the annular radii are nearly the same, and hence the orbit is close to
1 a circle. One can understand the stability question by studying the forces (restoring or unsettling) that come into play when the system is moved infinitesimally from 0.5 the position of circular orbit. The circular orbit is stable if, 0 ∂2V˜ > 0. (16) ∂r2 -0.5 (a) That is equivalent to, (b) -1 ∂2V˜ ∂f 3l2 = + > 0 (c) 2 4 ∂r − ∂r mr r0 -1.5 2 ∂f 3l r0 < 4 -2 ∂r | mr0 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 Using Eq. (13), FIG. 5: Spiralling orbits, f(r) ∼−1/r4 ∂f 3f(r0) r0 < (17) ∂r | − r0 particle away from the positions of zero radial velocity in an inward or an outward spiral. If the total energy IV. NATURE OF BOUNDED ORBITS is equal to the maximum or minimum of the effective potential, the system can stay with zero radial velocity, and hence move in a circular orbit. General bounded motion has both lower and upper The condition for circular orbit is, bounds. It means that the particle cannot approach nearer than some minimum or move farther than some ∂V˜ ∂V l2 maximum distance. One has to remember that the angu- = 3 = 0 (12) lar velocity has a constant sign, same as that of the con- ∂r ∂r r0 − mr0 r0 stant angular momentum, throughout the motion. How-
ever its magnitude decreases with increase in the radial whence we get, distance ( 1/r2). Together with the angular motion, the radial distance∼ changes from r to a r , then back ∂V l2 max min f(r0)= = (13) to rmax and so on. From this general nature of motion − ∂r −mr3 r0 0 two families of generic orbits are suggested in popular 11 texts (Fig. 6, Fig. 7). However it can be shown that The negative sign on the right hand side clearly shows the generic types shown in the later figure (Fig. 7) are that the force must be attractive. The particle moves in not feasible for any attractive potential. a circular orbit, since the force of attraction due to the We consider below a general bounded orbit confined in central field provides the necessary centripetal force. an annular region. We would like to investigate the na- For a given attractive central force ( f(r) and V (r) ture of the orbit close to the point where it touches the given ) it is possible to have a circular orbit of radius inner or the outer annulus. Let us choose the reference r0 provided the angular momentum and energy of the line (polar line) of the coordinate system such that the particle are given by, orbit touches the annular ring at θ = 0, or more specifi- 2 3 cally at r = R,θ = 0. Since at this point P (r = R,θ = 0) l = mr0f(r0) (14) the orbit is at its closest (or farthest) approach from the − 2 l pole, (∂r/∂t)P =0. As r is a function of θ, using Eq. (4) E = V (r0)+ 2 (15) 2mr0 we get for motion along the trajectory, dr l dr However when the effective potential has a maximum = the system is in an unstable circular orbit. A small devi- dt mr2 dθ ation from this radial position causes the orbit to be un- hence bounded. The effective force that comes into play makes the particle move away from the position of zero radial ′ dr r (0) = = 0 (18) velocity in an inward or outward spiral (Fig. 5). When dθ θ=0
6
2 1.5
1.5 1
1
0.5 0.5
0 0
-0.5 -0.5
-1
-1 -1.5
-2 -1.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 -1.5 -1 -0.5 0 0.5 1 1.5 1.5 1.5
1 1
0.5 0.5
0 0
-0.5 -0.5
-1 -1
-1.5 -1.5 -1.5 -1 -0.5 0 0.5 1 1.5 -1.5 -1 -0.5 0 0.5 1 1.5
FIG. 6: Generic orbits in central field (allowed) FIG. 7: Generic orbits in central field (not allowed)
The equation of motion for r is given by, One can expand r(θ) in a Taylor series in θ, and remem- ber that r′(0) = 0, mr¨ mrθ˙2 = f(r) − ′′ ′ ′′ θ From which we get r (θ) along the trajectory, r(θ) = r(0) + r (0)θ + r (0) + h.o. traj 2! l d l dr l2 mr(0)4 θ2 m mr = f(r) = r(0) + f˜(r(0)) + h.o. mr2 dθ mr2 dθ − mr3 2 l 2! d2r mr4 mr(0)4 θ2 = f˜(r) = r(0) + r(0) + f(r(0)) + h.o. 2 2 l2 2! dθ l mr4 = r + f(r) l2 Consider a tangent to the annulus at the point P , Fig. 7
r( θ ) r( θ ) tan tan r( θ ) Inner traj Outer r( θ ) annulus trajectory annulus traj trajectory θ (not allowed) θ (not allowed)
trajectory r (allowed) rannul annul
tangent tangent trajectory trajectory (allowed) (not allowed)
FIG. 8: Curvature of the trajectory in comparison to the FIG. 9: Curvature of the trajectory in comparison to the outer tangent to the inner annulus annulus
8, we study the orbit near the point P where it touches the outer annulus. r(0) 4 2 = cos(θ) mr(0) ˜ (∆θ) r(θ) r(∆θ)traj = r(0) + f(r) + h.o. l2 2! which can be expanded to give r(θ)tan along the tangent, hence
θ2 r(∆θ)traj < r(0) = rmax (20) r(θ)tan = r(0) + r(0) + h.o. 2 This confirms that the outer annulus is indeed the outer Notice that the constant or the θ independent terms in bound of the trajectory. The nature of the orbit near the r(θ)traj and r(θ)tan are equal, and the first leading power point of contact P at the outer annulus is shown in figure of θ is θ2 in both the cases. The coefficient of θ2 (cur- 9. vature) will determine the nature of the trajectory in reference to the tangent. We first study the curvature of the orbit in reference V. BOUNDED ORBITS FOR THE POWER LAW to the tangent. The force being attractive f(r(0)) < 0, CENTRAL FORCE for small angular distance (∆(θ)) away from the point P , A. Existence of stable circular orbit
r(∆θ)traj < r(∆θ)tan (19) For the case of a power law central potential k nk which means that the orbit bends more sharply than the V (r)= , f(r)= (21) tangent and stays closer to the inner annulus for nearby rn −rn+1 points (Fig. 8). Hence the second type of orbits shown From the stability condition Eq. (17), in some texts (Fig. 7) are not possible for central forces. For the outer annulus the analysis with respect to tan- ∂f 3f(r0) gent is not very meaningful as one can set up an even r0 < ∂r | − r0 stronger bound for its orbit, directly considering the po- tential V˜ (r). In this case the effective potential has a we find, positive slope with respect to the radius, and hence a (n + 1)nk k 1 negative effective force (f˜(r) < 0). The orbit should not n+2 < 3 n+1 r − −r r0 only remain nearer than the tangent, but even nearer 0 0 than the outer annulus. This condition is confirmed if n < 2 (22) 8
40 For the point of extremum
n−2 amn r0 = 2 . 30 l l^2/(2m r^2) Finding the second derivative of V˜ with respect to r at the point r0, 20 d2V˜ an(n + 1) 3l2 (a) unstable circular orbit = + dr2 − rn+2 mr4 10 2 (b) unbounded orbit 1 an(n + 1) 3l = 4 n−2 + ~ r − r m V(r) 1 l2 3l2 0 = (n + 1) + r4 − m m 0 1 l2 V(r) ~ - 1/r^3 -10 = 4 (2 n) < 0 r0 m −
Hence it is a point of maximum. For E > V˜ (r0) we -20 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 always have an unbounded orbit. For E < V˜ (r0) the orbit is semibounded, bounded above or bounded below, 4 FIG. 10: Effective potential for f(r) ∼−1/r force according to its initial state. The particle either spirals in or spirals out. For E = V˜ (r0) we get an unstable circular orbit. 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). 2. V (r)= −a/rn where n = 2
The effective radial potential is V˜ (r) = a/r2 + B. Study of boundedness of orbits l2/(2mr2). This is essentially an attractive or− repulsive inverse square term. 1. V (r)= −a/rn where n> 2 2 ˜ a˜ l V (r)= 2 , a>˜ 0 if a> Consider the case when n > 2. We have the effective −r 2m l2 potential, a<˜ 0 if a< 2m a l2 V˜ (r)= + (23) This potential cannot give circular orbit ever. If the ef- rn 2mr2 − fective potential is attractive and E < 0 it has an upper 1. V˜ (r) = 1/rn(a l2 rn−2/(2m)), and as r 0 bound of radial distance. A typical orbit would therefore we may− neglect l2− rn−· 2/(2m) in comparison to→a, be an inward spiral. For a repulsive effective potential · we will have outward spiral moving to infinite radial dis- thus limr→0 V˜ (r)= . −∞ tance.
2. limr→∞ V˜ (r)=0. n − 3. − where 3. For large but finite r, V˜ (r) = 1/r2(a/rn 2 V (r)= a/r 2 >n> 0 l2/(2m)), and a/rn−2 is negligible− compared to− l2/(2m), and hence V˜ (r) is positive. In this case the effective potential is, 2 4. V˜ (r) intersects the r axis at r∗ where (r∗)n−2 = ˜ a l V (r)= n + 2 (24) 2am/l2. −r 2mr ˜ ˜ n−2 2 We have the following properties of V (r) 5. V (r) has a maximum at r0 where r0 = anm/l . For a maximum point of V˜ we should have 1. V˜ (r)=1/r2(l2/2m a/rn−2) and as r 0 we may ˜ 2 ˜ 2 n−2 − 2 → ∂V /∂r r0 = 0 and ∂ V /∂r r0 < 0. neglect a/r in comparison to l /2m, and thus | | limr→0 V˜ (r) + . ∂V˜ an l2 → ∞ = =0 . ˜ ∂r rn+1 − mr3 2. limr→∞ V (r)=0. 9
30 30
20 20
l^2/(2m r^2) l^2/(2m r^2)
10 10 ~ V(r)
0 0
-10 ~ -10 V(r)
-20 -20 V(r) ~ ln(r)
V(r) ~ -1/r -30 -30
-40 -40 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
FIG. 11: Effective potential for f(r) ∼−1/r2 force FIG. 12: Effective potential for f(r) ∼−1/r force
− 3. V˜ (r) = 1/rn(a l2/2mr2 n) and for large but 4. V (r)= a ln r finite values− of r, l−2/2mr2−n is negligible compared to a, hence V˜ (r) is negative. This corresponds to one over r force f(r) a/r. The effective potential is given by, ∼− 4. V˜ (r) intersects the r axis at r∗ where (r∗)2−n = l2/(2am). 2 ˜ l V (r)= a ln r + 2 (25) 2−n 2mr 5. V˜ (r) reaches a minimum at r0 where r0 = l2/amn. We have the following properties of V˜ (r) For a minimum point of V˜ we should have ˜ 2 ˜ 2 1. limr→0 V˜ (r) + . ∂V /∂r r0 = 0 and ∂ V /∂r r0 > 0. | | → ∞ ˜ ∂V˜ an l2 2. limr→∞ V (r) + . = =0 . → ∞ ∂r rn+1 − mr3 3. There is no point of intersection with the r axis, ˜ For the point of extremum and V (r) is always positive.
2 4. V˜ (r) reaches a minimum at r0 = l/√am. 2−n l r0 = . ˜ amn For a minimum point of V we should have ˜ 2 ˜ 2 ∂V /∂r r0 = 0 and ∂ V /∂r r0 > 0. Finding the second derivative of V˜ with respect to | | ∂V˜ a l2 r at the point r0, = ∂r r − mr3 d2V˜ an(n + 1) 3l2 = + for the point of extremum dr2 − rn+2 mr4 1 an(n + 1) 3l2 l2 = + r2 = r4 − rn−2 m 0 am 2 2 1 l 3l ˜ = (n + 1) + The second derivative of V at r = r0, r4 − m m 0 2 d2V˜ a 3l2 1 l = + = 4 (2 n) > 0 dr2 −r2 mr4 r0 m − 1 3l2 = a + Hence it is a point of minimum. For this case the orbit r2 − mr2 can be bounded or unbounded depending on the total 1 = 2a> 0 energy E of the particle. r2 · 10
30 The second derivative of V˜ with respect to r, 2 ˜ 2 d V n−2 3l 25 = an(n 1)r + dr2 − mr4 1 3l2 = an(n 1)rn+2 + 20 r4 − m 1 l2 = 4 (2 + n) > 0 15 r m Hence it is a point of minimum. The orbit is always
10 V(r) ~ r^2 bounded. ~ V(r) The findings for the general power law potential can be summarized as follows, 5 V (r)= sign(n)arn n = 0 (27) l^2/(2m r^2) − 6 f(r)= abs(n)arn 1 (28) 0 − and
-5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 V (r)= b ln r (29) b f(r)= (30) FIG. 13: Effective potential for f(r) ∼−r force −r
Hence it is a point of minimum. The orbit is always TABLE I: Dependence of boundedness on power law. bounded. n of V (r) ∼ arn nature of boundedness
) −∞, −2( always unbounded
−2 spiralling orbit 5. V (r)= arn where n> 0 ) − 2, 0( bounded or unbounded depending on E For this potential the effective potential is given by, 0 (ln r) always bounded l2 V˜ (r)= arn + n> 0 (26) )0, +∞( always bounded 2mr2
We have the following properties of V˜ (r) C. Stable bounded orbit, geometric shape 1. limr→0 V˜ (r) + . → ∞ The stable bounded orbits in a power law central field 2. limr→∞ V˜ (r) + . can have the following forms → ∞ 1. V (r)= a/rn, 0 orbits are sometime mistakenly identified with diagrams TABLE II: Dependence of stability of circular orbits on power of the form shown in Fig. 711. The above expression in law. ∼ n fact corresponds to a class of orbits that look generically n of V (r) ar stability of circular orbits like those shown in Fig. 6. The generic features of central force orbits discussed ) −∞, −2) unstable here have been verified by computer simulation for a large ) − 2, ∞( stable class of central force fields. The figures shown here (Fig. 2, 3, 5 and 6) were generated by these simulations. Stu- dents interested in studying and generating such orbits VI. CONCLUSION will find Ref.12 helpful. Existence of bounded orbit for a large class of attrac- Acknowledgement tive 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 The authors wish to express their indebtedness to the 1 2 9 texts1, are shown to be non-feasible. well known texts by Goldstein , Landau and Arnold . Small deviation from circularity in the case of central They also acknowledge their teachers in related graduate field is often expressed in terms of inverse of radial dis- courses at Stony Brook, Prof. Max Dresden, Prof. A. S. tance (u =1/r). Goldhaber and Prof. Leon A. Takhtajan. Authors gratefully acknowledge the encouragement re- u = u0 + a cos(βθ) ceived from Prof. Shyamal SenGupta of Presidency Col- · r lege, Calcutta. The material presented here was used in r = 0 1+ a r0 cos(βθ) a graduate level classical mechanics course at Jadavpur · · University during 1998-2001. SR wishes to thank his stu- where r0 = 1/u0. It is interesting to note that these dents for stimulating discussions. ∗ Electronic address: [email protected] 8 L.A. Pars, Introduction to Dynamics, Cambridge, 1953. 1 H. Goldstein, Classical Mechanics, Addison-Wesley, 1980. 9 V.I. Arnold, Mathematical Methods of Classical Mechanics, 2 L.D. Landau, E.M. Lifshitz, Mechanics, Pergamon Press, Springer Verlag, 1989. 1976. 10 R. Abraham, and J.E. Marsden, Foundations of Mechan- 3 K.R. Symon, Mechanics, Addison-Wesley, 1971. ics, Benjamin-Cummings Publ.Co.Inc., 1978. 4 D.T. Greenwood, Classical Dynamics, Prentice Hall, 1977. 11 H. Goldstein, Fig. 3-7, Fig. 3-13 and Eq. 3-45 in Classical 5 J.L. Synge, and B.A. Griffith, Principles of Mechanics, Mc- Mechanics, Addison-Wesley, 1980. GrawHill, 1970. 12 H. Gould, J. Tobochnik, An Introduction to Computer 6 A. Sommerfeld, Mechanics, Academic Press, 1952. Simulation Methods, Addison-Wesley, 1988. 7 E.T. Whittaker, A Treatise on the Analytical Dynamics of Particles and Rigid Bodies, Dover, 1944.