* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download RHIG - Wayne State University
Atomic theory wikipedia , lookup
Path integral formulation wikipedia , lookup
Photon polarization wikipedia , lookup
Fictitious force wikipedia , lookup
Relativistic quantum mechanics wikipedia , lookup
Analytical mechanics wikipedia , lookup
Old quantum theory wikipedia , lookup
Centrifugal force wikipedia , lookup
Brownian motion wikipedia , lookup
Classical mechanics wikipedia , lookup
N-body problem wikipedia , lookup
Seismometer wikipedia , lookup
Lagrangian mechanics wikipedia , lookup
Relativistic angular momentum wikipedia , lookup
Routhian mechanics wikipedia , lookup
Hunting oscillation wikipedia , lookup
Laplace–Runge–Lenz vector wikipedia , lookup
Rigid body dynamics wikipedia , lookup
Relativistic mechanics wikipedia , lookup
Theoretical and experimental justification for the Schrödinger equation wikipedia , lookup
Centripetal force wikipedia , lookup
Newton's theorem of revolving orbits wikipedia , lookup
Newton's laws of motion wikipedia , lookup
Central-Force Motion Chapter 8 Prof. Claude A Pruneau Physics and Astronomy Department Wayne State University 8.1 Introduction • Consider motion of two particles affected by a force connecting the center of the two bodies. • One of few problems that can be solved completely. • Historically important e.g motion of planets, alpha-particle scattering on nuclei 8.2 Reduced Mass • Description of a two-particle system • Discussion restricted to frictionless (conservative) systems. 8.2 Reduced Mass (cont’d) • Assume a force is acting between the two particles along a line joining them. • Particle positions described in some arbitrary reference frame as positions r1 and r2, or in terms of the CM frame position, R, and the relative position vector r = r1 - r2 . ! r1 !m r1 m1 ! r CM ! R ! r2 Arbitrary Frame 1 ! R=0 ! r2 CM m2 m2 CM Frame Lagrangian for a two-body system • Assume the potential energy is only a function of the distance between the particles, U=U(r) • Lagrangian may be written: 2 2 L = m1 r!1 + m2 r!2 ! U(r) 1 2 1 2 • Translational motion of the system uninteresting; – Use R=0. Two-body CM Coordinates ! ! m1r1 + m2 r2 = 0 • We have Solving for r1 and r2 : ! r1 = m2 ! r m1 + m2 ! ! ! r = r1 ! r2 ! r2 = ! m1 ! r m1 + m2 Substitute in the Lagrangian r1 and r2 : 2 2 L = 12 m1 r!1 + 12 m2 r!2 ! U(r) 2 2 m2 !" m1 !" r + 12 m2 r ! U(r) m1 + m2 m1 + m2 m22 m12 !" 2 1 !" 2 1 L = 2 m1 r + m r ! U(r) 2 2 2 2 ( m1 + m2 ) ( m1 + m2 ) m2 m1 !" 2 1 !" 2 L = 12 m1m2 r + m m r ! U(r) 1 2 2 2 2 ( m1 + m2 ) ( m1 + m2 ) L = 12 m1 L= 1 m1m2 !" 2 r ! U(r) 2 m1 + m2 !" 2 L = µ r ! U(r) 1 2 m1m2 µ= m1 + m2 ! reduced mass 2 to 1 reduction 8.3 Conservation Theorems • Particle of mass µ in a central force field described by the potential function U(r). • Symmetry implies conservation of angular momentum. ! ! L = r ! p = constant Radius vector and momentum lie in a plane normal to the angular momentum vector L. The problem reduced to 2 dimensions: i.e. along “r” and “q”. ! L ! r ! p • Lagrangian ( ) L = 12 µ r! 2 + r 2!! 2 " U (r) Lagrangian cyclic in q implies: Angular momentum, pq, conjugate to q, is a conserved quantity. "L d "L p!! = =0= "! dt "!! #L 2! p! " = µr ! = constant #!! First integral of motion • The system’s symmetry permits the integration of one equation of motion. • pθ is the first integral of the motion. • Denote it l = µr !! = constant 2 Note • l can be negative or positive Interpretation of l as Areal velocity • The radius vector sweeps out an area dA in a time interval dt. dA = r d! 1 2 ! r (t1 ) d! rd! ! r (t 2 ) The areal velocity is thus dA 1 2 d! 1 2 ! = 2r = 2r ! dt dt l = = constant 2µ 2 ! r (t) Kepler’s 2nd law of planetary motion dA = constant dt • No particular assumption made about the form of U(r) implies: • This result is NOT limited to an inverse-square law force but is valid for all central forces. • Since the motion of the CM is not interesting, only one degree of freedom remains to be considered. • Linear momentum conservation adds nothing new here… • Energy conservation provides the only remaining equation of motion. T + U = E = constant E= 12 µv 2 + U (r) 2 2 !2 1 E= 2 µ r! + r ! + U (r) ( Total Energy ) 2 1 l E= 12 µr! 2 + + U (r) 2 2 µr 8.4 Equations of Motion • Assume U(r) is specified. Solve for dr/dt: 1 l2 2 1 E= 2 µr! + + U (r) 2 2 µr dr 2 l2 r! = =± E ! U (r) ! 2 2 dt µ µr ( ) Solving for dt, and integrate to get a solution t = t(r). Invert it to get r = r(t) Alternatively, obtain q= q(r), starting with d! dt !! d! = dr = dr dt dr r! dr 2 l2 r! = =± E ! U (r) ! 2 2 dt µ µr ( ) ! (r) = ) d! dt !! d! = dr = dr dt dr r! ( l = µr 2!! = constant ) ± l / r 2 dr # l2 & 2 µ % E " U (r) " 2 µ r 2 (' $ • Inversion of the result (if possible) yields the standard form (general) solution r = r(t). • Because l is constant, dθ/dt is a monotonic function of time. • The above integral is in practice possible only for a limited number of cases… Remarks • with F(r) = rn, solutions may be expressed in terms of elliptic integrals for certain integers and fractional values of “n”. • Solution may be expressed in terms of circular functions for n=1, -2, and –3. • Case n = 1 is the harmonic oscillator. • Case n = -2 is the inverse square law. Solution using Lagrange equations !L d !L " =0 !r dt !r! Lagrange equation for “r” ( ) #U 2 ! !! µ r ! r" = ! = F(r) #r 1 Use variable change u = r du 1 dr 1 dr dt 1 r! =" 2 =" 2 =" 2 ! d! r d! r dt d! r ! Remember l = µr !! 2 l ! != 2 µr Compute… l ! != 2 µr du 1 r! µ = " 2 ! = " r! d! l r ! d 2u d # µ & dt d # µ & µ r!! = " r!( = " r!( = " ! 2 % % d! $ l ' d! dt $ l ' l! d! d 2u µ r!! µ2 2 =" = " 2 r r!! 2 l # l & d! l %$ µr 2 (' • Solving l 2 2 d 2u r!! = ! 2 u µ d" 2 2 l r"! 2 = 2 u 3 µ ( ) #U 2 ! = F(r) • Substitute back into µ r!! ! r" = ! #r d 2 " 1% 1 µr 2 + =( 2 F r 2 $ ' d! # r & r l () • Which is useful if one wishes to find the force law that produces a particular orbit r=r(θ). Example 8.1 – Log-spiral Find the force law for a central-force field that allows a particle to !" move in a logarithmic spiral orbit given by, r = ke , where k and α are constants. d " 1% d " e()! % () e()! Solution = = d! $# r '& d! $# k '& k Calculate Now use To find: d 2 " 1 % ) 2 e()! ) 2 = = 2 $ ' k r d! # r & d 2 " 1% 1 µr 2 + =( 2 F r d! 2 $# r '& r l () !1 !1 " µr 2 % " d 2 " 1 % 1 % " µr 2 % " ) 2 1 % F r = $! 2 ' $ 2 $ ' + ' = $! 2 ' $ + ' r& # l & # d( # r & r & # l & # r () ( ) l2 = ! 3 )2 +1 µr Force is Attractive and Inverse cube! Example 8.2 – r(t), θ(t) Determine the functions r(t) and θ(t) for the problem in Ex 8.1. Solution: Start with: l l ! ! = 2 = 2 2"! µr µk e Rearrange, integrate: e 2!" l d" = dt 2 µk Answer: #e 2!" l d" = # 2 dt µk & 1 # 2" lt ! (t) = ln % + C( 2 2" $ µ k ' e2!" lt = + C' 2 2! µk Similarly for r(t), remember And write Answer (2): l l = 2 µr µ k 2 e2!" r2 2! lt 2!" =e = +C 2 2 k µk r(t) = 2! lt + k 2C µ Where l and C are determined by the initial conditions Example 8.3 – Total Energy What is the total energy of the orbit of the previous two examples? Solution: Need U… ( ) +l 2 2 U r = ! " Fdr = # +1 µ () ( ) l2 1 2 U (r) = ! # +1 2 2µ r l l ! ! = 2 = 2 2"! µr µk e r = ke!" d! d! dr l !! = = = 2 dt dr dt µr dr l l "l "! r! = = 2 " ke = 2 d! µr µr µr dr " r3 () lim U r = 0 r !" 2 1 " !l % l E = µ$ ' + ( 2 2 # µr & 2µr Given the reference 2 ( 2µr 2 () lim U r = 0 r !" )=0 l2 ! 2 + 1 8.5 Orbits in a central field • Radial velocity of a particle in central field dr 2 l2 r! = =± E ! U (r) ! 2 2 dt µ µr ( • ) Vanishes at the roots of the radical 2 l2 E ! U (r) ! 2 2 = 0 µ µr ( ) l2 E ! U (r) ! =0 2 2µr l2 E ! U (r) ! =0 2 2µr • Vanishing of dr/dt implies turning points • Two roots in general: rmin and rmax. • Motion confined to an annular region between rmin and rmax. • Certain combinations of E and l may lead to a single root: one then has a circular motion, and dr/dt = 0 at all times. • Periodic motion in U(r) implies the orbit is closed; I.e. loops on itself after a certain number of excursions about the center of force. • The change in θ while going from rmin to rmax is a function of the potential and need not be 180o. • It can be calculated! • Because the motion is symmetric in time: !" = * rmax rmin ( ) ± l / r 2 dr $ l2 ' 2 µ & E # U (r) # 2) 2 µ r % ( • Path closed only if Dq is a rational fraction of 2p. • Dq= 2p(a/b) where a and b are integers. • In this case, after b periods the particle will have completed a revolutions and returned to its original position. n+1 • For U (r) ! r a closed noncircular path exists only for n = -2 or +1. 8.6 Centrifugal Energy and Effective Potential • In dr/dt, dθ/dt, …, we have l2 E ! U (r) ! 2µr 2 • Where each term has the dimension of energy. • Remember that • Write 2 l = µr !! l 1 2 !2 = µr ! 2 2 2µr 2 l2 1 2 !2 = µr ! • Interpret 2 2 2µr l2 energy” U c = 2µr 2 as a “potential "U c l2 2 ! F = ! = = µ r # • The associated force is: c "r µr 3 • Traditionally called a centrifugal force. – Although it is, STRICTLY SPEAKING, NOT A FORCE – but rather a pseudo-force. – We continue to use the term nonetheless… • The term l 2 2 µr 2 can then be interpreted as the centrifugal potential energy, and included with U(r) to define an effective potential energy. 2 l V (r) ! U (r) + 2µr 2 • V(r) is a fictitious potential that combines the real or actual potential U(r) with the energy term associated with the angular motion about the center of force. • For an inverse-square law central-force motion, one gets: 2 "k l V (r) ! + 2 r 2µr l2 1 2! = µr ! 2 2 2µr V (r) k ! r V (!) = 0 Energy V (r) Turning point(s) (apsidal distances) unbound r3 r1 r2 bound r4 E1 1 2 µr! 2 r E2 E3 • Values of E less than Vmin = ! µ k 2l do not result in physically real motion; given velocity is imaginary. 2 2 • Techniques illustrated here are used in modern atomic, molecular and nuclear physics (but in the context of QM). 8.7 Planetary Motion – Kepler’s Problem • Consider the specific case of an inverse-square force law. ± l / r 2 dr ! (r) = ) + constant # k l2 & 2µ % E + " r 2 µr 2 (' $ ( ) • Integral soluble with variable substitution u=1/r. • Define the origin of θ so r is a minimum. cos! = l2 1 "1 µk r 2El 2 1+ µk 2 cos! = l2 1 "1 µk r 2El 2 1+ µk 2 • Define constants l2 != µk • Then one can re-write: 2El 2 ! = 1+ µk l2 1 2El 2 = 1+ 1+ cos! 2 µk r µk • To get the equation of a conic section with one focus at the origin ! = 1 + " cos# r ! = 1 + " cos# r • The quantity, ε, is called eccentricity, and • 2α is termed the latus rectum of the orbit. • Conic sections are formed by the intersection of a plane and a cone. • More specifically … by the loci of points (formed by a plane) where the ratio of the distance from a fixed point (the focus) to a fixed line (called the directrix) is a constant. Parabola, ε=1 Hyperbola, ε>1 Ellipse, 0<ε<1 Directrix For parabola Circle, ε=0 ! = 1 + " cos# r • q=0 corresponds to a pericenter, i.e. rmin, whereas rmax corresponds to the apocenter. • The general term for turning points is apsides. • Planetary Motion: • Major axis ! k a= • Minor axis b= 1" # 2 ! 1" # 2 = = 2E l 2µ E # rmin = a 1 ! " = 1+ " # rmax = a 1 + " = 1! " ( ) ( ) a b P aε P α • Period of elliptic motion: 2µ dt = dA l T A 2µ ! dt = l ! dA 0 0 2µ T= A Aab!= l 2µ 2µ k • The area of an ellipse is: T = ! ab = ! l l 2E • The period is then…. • Noting b = ! a • One also finds: l 2µ E µ "3/ 2 T = !k E 2 2 4 ! µ 3 2 T = a k Kepler’s Third Law • Given the gravitational force: F(r) = ! Gm1m2 r2 k =! 2 r 2 3 2 3 4 ! a µ 4 ! a • The square of the period: T 2 = " Gm2 G m1 + m2 ( ) • Where the last approx is realized for m1 << m2. • Kepler’s statement is correct only if the mass m1 of a planet can be neglected with respect to the mass m2 of the sun. • Correction needed for Jupiter given that it is 1/1000 of the mass of the Sun. Kepler’s Laws 1. Planets move in elliptical orbits about the sun with the sun at one focus. 2. The area per unit time swept out by a radius vector from the sun to a planet is constant. 3. The square of a planet’s period is proportional to the cube of the major axis of the planet’s orbit. Example 8.4 Halley’s comet, which passed around the sun early in 1986, moves in a highly elliptical orbit with an eccentricity of 0.967 and a period of 76 years. Calculate its minimum and maximum distances from 1 the sun. 2 2 $ ' Solution: ! = Gmsun" & ) 2 % 4# ( 2 ,$ ' $ 365day 24hr 3600s ' +11 Nm 30 . & 6.67 * 10 1.99 * 10 kg & 76 yr 2 ) yr day hr )( % kg ( .% =. 2 4 # . . ! = 2.68 * 1012 m ( We thus find: ) ( ) m (1 + 0.967 ) = 5.27 ! 10 rmin = 2.68 ! 1012 m 1 " 0.967 = 8.8 ! 1010 m rmax = 2.68 ! 1012 10 m 2 / 1 1 1 1 1 0 1 2 8.8 Orbital Dynamics