Central Force Motion

Challenge Problems

Problem 1: Elliptic Orbit

A satellite of mass ms is in an elliptical orbit around a planet of mass m p which is

located at one focus of the ellipse. The satellite has a velocity va at the distance ra when
it is furthest from the planet.

a) What is the magnitude of the velocity v p of the satellite when it is closest to the
planet? Express your answer in terms of ms , m p , G , va , and ra as needed.

b) What is the distance of nearest approach rp ?

c) If the satellite were in a circular orbit of radius r0 = rp , is it’s velocity v0 greater

than, equal to, or less than the velocity v p of the original elliptic orbit? Justify
your answer.

Problem 1 Solution:

The angular momentum about the origin is constant therefore

L ! Lp = La . (28.1.1)

Because ms << m p , the reduced mass µ ! ms and so the angular momentum condition
becomes

L = ms rp v p = ms ra v a (28.1.2)
Thus
rp = ra va / v p . (28.1.3)

Note: We can solve for v p in terms of the constants G , m p , ra and va as follows. Choose
zero for the gravitational potential energy for the case where the satellite and planet are
separated by an infinite distance, U (r = ! ) = 0 . The gravitational force is conservative,
so the energy at closest approach is equal to the energy at the furthest distance from the
planet, hence

1 Gms m p 1 Gms m p
ms va 2 ! = ms v p 2 ! . (28.1.4)
2 ra 2 rp

Substituting Eq. (28.1.3) into Eq. (28.1.4), yields

1 Gms m p 1 Gms m p v p
ms va 2 ! = ms v p 2 ! . (28.1.5)
2 ra 2 ra va

After a little algebraic manipulation, Eq. (28.1.5) becomes

2Gmp
(v p ! va ) = (v p + va )(v p ! va ) . (28.1.6)
ra va

We can solve this equation for v p

2Gm p
vp = ! va . (28.1.7)
ra va

b) We can use Eq. (28.1.3) to find the distance of nearest approach,

ra va 1
rp = = ra . (28.1.8)
vp ! 2Gm p "
$ 2
#1%
& ra va '

c) The speed of the satellite undergoing uniform circular motion can be found from the
force equation,

Gms mp ms v0 2
! =! . (28.1.9)
r0 2 r0

So the speed is

v0 = Gmp / r0 . (28.1.10)
Note that if we assume our original orbit was circular and set va = v0 and ra = r0 in Eq.
(28.1.8) then the distance of closest approach becomes

1
rp = r0 = r0 (28.1.11)
" 2Gmp %
$ 2
! 1'
# r0 v0 &

providing a second check on our algebra.

If we substitute Gm p = r0 v0 2 into Eq. (28.1.7), we have that

2r0 v0 2
vp = ! va . (28.1.12)
ra va

Let’s compare our elliptic orbit to a circular orbit with rp = r0 . How does v p compare to
v0 ? If we substitute rp = r0 into Eq. (28.1.12) we have that

2rp v0 2
vp = ! va . (28.1.13)
ra va

We now use the angular momentum condition that ra va = rp v p to rewrite Eq. (28.1.13) as

2v0 2
vp = ! va . (28.1.14)
vp
Thus
v p 2 + v p va = 2v0 2 . (28.1.15)

We know that va < v p , so v p 2 + v p va < 2v p 2 . Thus 2v0 2 < 2v p 2 or

v0 < v p . (28.1.16)

So the circular orbit with rp = r0 has speed v0 less than the speed v p of closest approach
for the elliptic orbit. This is not surprising because suppose we give the satellite that is in
a circular motion a small increase in velocity in the tangential direction, vnew = v0 + !v .
This implies that the energy increases. The satellite will no longer travel in a circular
orbit since for the same radius r0 ,
 ms vnew 2 Gms mp
> . (28.1.17)
r0 r0 2

The satellite will move away from the planet entering into an elliptical orbit. So any
velocity v p greater than v0 will form an elliptic orbit.
Problem 2: Planetary Orbits

Comet Encke was discovered in 1786 by Pierre Mechain and in 1822 Johann Encke
determined that its period was 3.3 years. It was photographed in 1913 at the aphelion
distance, ra = 6.1!1011 m , (furthest distance from the sun) by the telescope at Mt. Wilson.
The distance of closest approach to the sun, perihelion, is rp = 5.1!1010 m . The universal
gravitation constant G = 6.7 "10!11 N # m 2 # kg !2 . The mass of the sun is ms = 2.0 !1030 kg .

a) Explain why angular momentum is conserved about the focal point and then write
down an equation for the conservation of angular momentum between aphelion
and perihelion.

b) Explain why mechanical energy is conserved and then write down an equation for
conservation of energy between aphelion and perihelion.

c) Use conservation of energy and angular momentum to find the speeds at

perihelion and aphelion.

Problem 2 Solutions

Since the only forces are gravitational and point toward the focal point
!
!
"!
! focal = r focal,comet " F grav = 0

Therefore, angular momentum is conserved. At perihelion, since v p is tangent to the orbit

L p = µrp v p

At apohelion:
La = µra! a
 La = L p ! µra va = µrp v p !

ra va = rp v p

The gravitational force is conservative so mechanical energy is conserved.

1 2 Gm1m2 1 Gm1m2
Ea = µva = , E p = µv 2p =
2 ra 2 rp

thus
Ea = E p !

1 2 Gm1m2 1 2 Gm1m2
µva ! = µv p !
2 ra 2 rp

Since the mass m1 of the comet is much less than the mass of the sun m2.

m1m2 mm
µ= ! 1 2 = m1
m1 + m2 m2

the energy conservation equation becomes:

1 Gm1m2 1 Gm1m2 1 2 Gm2 1 2 Gm2

m1va2 ! = m1v 2p ! or va ! = vp !
2 ra 2 rp 2 ra 2 rp
rv
The condition v p = a a from conservation of angular momentum can now be used in the
rp
energy equation.

2
1 2 Gm2 1 " ra va %
Gm2
va ! =
$ !

2
ra 2 #
rp '&
rp

Solve for va :
1
!
!
ra $
$ 2 !
1
1
$

# 2 ' 1& a
v = Gm '
&
2
"
#"
rp &%

2#
% "
rp ra %

 1/2
" " 1 1%%
$ 2Gm2 $ ! ' '
$ # rp ra & '
va = $ 2 '
$ " ra % ! 1 '
$# #$ rp '& '&

1/2
" 2Gm2 rp %
va = $ '
( )
# ra ra + rp &

( )( )( ) '%
1/2

" ( 2 ) 6.67 ( 10 !11 N ) m 2 ) kg !2 2.0 ( 10 30 kg 5.1 ( 1010 m

va = $
= 5.8 ( 10 3 m ) s !1
$# ( ) ((
6.1 ( 1011 m 6.1 ( 1011 m + 5.1 ( 1010 m )) '&

The velocity at perihelion is then

(6.1 ! 1011 m )
vp =
ra va
= ( )
5.8 ! 10 3 m " s #1 = 6.9 ! 10 4 m " s #1
rp (5.1 ! 10 m
10
)
Problem 3: Satellite Motion

A spherical non-rotating planet (with no atmosphere) has mass m1 and radius r1 . A

projectile of mass m2 << m1 is fired from the surface of the planet at a point A with a
speed vA at an angle ! = 30! with respect to the radial direction. In its subsequent
trajectory the projectile reaches a maximum altitude at point B on the sketch. The
distance from the center of the planet to the point B is r2 = (5 / 2)r1 .

In this problem you will find the initial speed vA in terms of G , m1 and r1 .

a) Is there a point about which the angular momentum of the projectile is constant?
If so, use this point to determine a relation between the speed vB of the projectile
at the point B in terms of vA and the angle ! = 30! .

b) Now use conservation of mechanical energy constant to find an expression to find

vA in terms of G , m1 and r1 .

Problem 3 Solutions:

a) The only force of interest is the gravitational force, which is always directed toward
the center of the planet; hence angular momentum about the center of the planet is a
constant.

At point A, the component p! of the satellite’s linear momentum perpendicular to the

radius vector is

m2 vA
p! = m2 vA sin " = , (28.3.1)
2
Using sin 30! = 1 / 2 . The magnitude of the angular momentum about the center of the
planet is then

r1 m2 vA
LA = r1 p! = . (28.3.2)
2

At point B (the apogee), the velocity vector is perpendicular to the radius vector and the
magnitude of the angular momentum is the product of the distance from the center of the
planet and the speed,

5
LB = r2 vB = rv . (28.3.3)
2 1 B

There is no torque on the satellite, so LB = LA ; so equating the expressions in Equations

(28.3.3) and (28.3.2) yields

5 rm v
r1 vB = 1 2 A
2 2 (28.3.4)
vA
vB = .
5

b) The mechanical energy E of the satellite as a function of speed v and radius (distance
from the center of the planet) r is

1 mm
E= m2 v 2 ! G 1 2 . (28.3.5)
2 r

Equating the energies at points A and B, and using rB = r2 = (5 / 2)r1 , rA = r1 and

vB = vA / 5 from part a) (Equation (28.3.4) above),
 1 mm 1 mm
m2 vA2 ! G 1 2 = m2 vB2 ! G 1 2
2 rA 2 rB
2
1 2 m 1" v % m
vA ! G 1 = $ A ' ! G 1
2 r1 2# 5 & 5r1 / 2
1 2" 1% m " 2%
vA $ 1 ! ' = G 1 $ 1 ! ' . (28.3.6)
2 # 25& r1 # 5&
5 m1
vA2 = G
4 r1
5 G m1
vA =
4 r1

It’s worth noting that vA is less than the escape velocity

2G m1
vescape = (28.3.7)
r1

of the planet, but not by much;

vA 5
= ! 0.79 . (28.3.8)
vescape 8

1 mm 1 mm
m2 vA2 ! G 1 2 = m2 vB2 ! G 1 2
2 rA 2 rB
2
1 2 m 1" v % m
vA ! G 1 = $ A ' ! G 1
2 r1 2# 5 & 5r1 / 2
1 2" 1% m " 2%
vA $ 1 ! ' = G 1 $ 1 ! ' (28.3.9)
2 # 25& r1 # 5&
5 m1
vA2 = G
4 r1
5 G m1
vA =
4 r1

It’s worth noting that vA is less than the escape velocity

 2G m1
vescape = (28.3.10)
r1

of the planet, but not by much;

vA 5
= ! 0.79 . (28.3.11)
vescape 8

Problem 4: Inverse Square Central Force: Lowest Energy Solution

The effective potential energy for an inverse-square restoring central force

! Gm1m2
F1,2 = ! rˆ is given by
r2

L2 Gm1m2
U effective = !
2µ r 2
r

Make a graph of the effective potential energy U effective (r ) as a function of the relative
separation r . Find the radius and the energy for the motion with the lowest energy. What
type of motion does this correspond to?

Problem 4 Solution:

For plotting purposes, the horizontal scale is the ratio r / r0 and the vertical scale is in
units of U effective /U effective (r0 ) , as found in the second part of the problem.

The upper (red) curve is proportional to 1/ r 2 and the lower (blue) curve is proportional
to !1/ r . The sum is the solid (green) curve in between.

To find the minimum energy, differentiate the effective potential U effective (r ) with respect
to the radius r and set the derivative equal to zero at r0 ,
 L2 1 Gm1 m2
! + =0
µ r03 r02
L2
r0 = .
µ Gm1 m2

The energy at this minimum value is

2
L2 ! µ Gm1m2 " µ Gm1m2
U effective (r0 ) = $ % # Gm1m2
2 µ & L2
' L2
µ (Gm1m2 )
2

=# .
2L2

The radius of this orbit does not change; the orbit is a circle.
Problem 5

A particle of mass m moves under an attractive central force of magnitude F = br 3 . The

angular momentum is equal to L .

a) Find the effective potential energy and make sketch of effective potential energy
as a function of r .

b) Indicate on a sketch of the effective potential the total energy for circular motion.

c) The radius of the particle’s orbit varies between r0 and 2r0 . Find r0 .

Problem 5 Solutions:

a) The potential energy is, taking the zero of potential energy to be at r = 0 , is

b
U (r ) = " # ("br !3 )dr ! = r 4
r

0 4

and the effective potential is

L2 L2 b
U eff (r ) = 2
+ U ()
r = 2
+ r4 .
2mr 2mr 4

A plot is shown below, including the potential (yellow if seen in color), the term L2 / 2m
(green) and the effective potential (blue). The minimum effective potential energy is the
horizontal line (red). The horizontal scale is in units of the radius of the circular orbit and
the vertical scale is in units of the minimum effective potential.

b) See the solution to part (a) above and the plot to the left below.
c) In the left plot, if we could move the red line up until it intersects the blue curve at two
point whose value of the radius differ by a factor of 2 , those would be the respective
values for r0 and 2r0 . A graph of this construction (done by computer, of course),
showing the corresponding energy as the horizontal magenta is at the right above, and is
not part of this problem.

To do this algebraically, we find the value of r0 such that U eff (r0 ) = U eff 2 (r0 ). This is

L2 b 4 L2 b
+ (2r0 ) .
4
+ r =
m (2r0 ) 4
2 0 2
mr0 4

Rearranging and combining terms, and then solving for r0 ,

3 L2 1 15 4
= b r0
8 m r02 4
1 L2
r06 = .
10 mb

( )
Thus, r0 = 1/ 10 rcircular (not part of the problem), consistent with the auxiliary figure on
the right above.
Problem 6:

The effective potential corresponding to a pair of particles interacting through a central

force is given by the expression

L2
U eff (r ) = + Cr 3 (30.1)
2µ r 2

where L is the angular momentum, µ is the reduced mass and C is a constant. The total
energy of the system is E . The relationship between U eff (r) and E is shown in the
figure, along with an indication of the associated maximum and minimum values of r
and the minimum allowed energy Emin . In what follows, assume that the center of mass
of the two particles is at rest.

a) Find an expression for the radial component f (r) of the force between the two
particles. Is the force attractive or repulsive?

b) What is the radius r0 of the circular orbit allowed in this potential? Express your
answer as some combination of L , C , and µ .

c) When E has a value larger than Emin , find how rapidly the separation between
the particles is changing, dr / dt , as the system passes through the point in the
orbit where r = r0 . Give your answer in terms of some combination of E , Emin ,
L , C , µ and r0 .

d) Does the relative motion between the particles stop when r = rmax ? If not, what is
the total kinetic energy at that point in terms of some combination of E , L , C ,
µ , rmax and rmin ?
Problem 6 Solutions:

a) The effective potential is given by

L2
U eff (r ) = + U (r ) , (28.6.2)
2µ r 2

and f (r ) = !dU / dr for a central force, and so

d d
f (r ) = ! U (r ) = ! Cr 3 = !3Cr 2 . (28.6.3)
dr dr

From the figure, C > 0 , so f (r ) < 0 , a restoring force.

b) The circular orbit will correspond to the minimum effective potential; at this radius the
kinetic energy will have no contribution from any radial motion. This minimum effective
potential, and hence the radius of the circular orbit, is found from basic calculus
and algebra,

!d " L2
$& dr U eff (r )'% = # 3 + 3Cr02
µ r0
r = r0
1/ 5
(28.6.4)
L2
( L ) 2
r05 = , r0 = * + .
3µC , 3µ C -

c) Recall that the kinetic energy is

2
L2 1 ! dr "
K= + µ# $ . (28.6.5)
2µ r 2
2 % dt &

The difference E ! Emin is then found by evaluating U eff at r = r0 ,

2
1 ! dr "
E # Emin = µ$ %% , (28.6.6)
2 $& dt r=r0 '

or dr / dt =
(2 / µ )(E ! Emin ) .

d) No; dr / dt = 0 , but the kinetic energy, from Equation (28.6.5), is

 L2

K min = . (28.6.7)
2 µ rmax
2

Problem 7: Determining the Mass of a Neutron Star

A binary system known as 4U0900-40 consists of a “neutron star” and a normal

“optical star.” You are given two graphs of actual data obtained from observations of this
system. The top graph shows the time delays (in seconds) of X-ray pulses detected from
the neutron star as a function of time (in days) throughout its orbit. These delays indicate
the time of flight for an X-ray pulse (traveling with the speed of light) to cross the orbit of
the binary in its trip to the earth. (Ignore the heavy dots scattered about the “ x ”-axis.)
The bottom graph displays the velocity component of the optical star toward (or
away from) the Earth as a function of its orbital phase. The velocities (in km per second)
are determined from Doppler shifts in its spectral lines. Time on this plot is given in units
of orbital phase, where the time between phases 0.0 and 1.0 corresponds to one orbital
period.
Assume that the orbit of 4U0900-40 is circular and that we are viewing the system
edge on, i.e., the Earth lies in the plane of the binary orbit. In each graph, the solid curve
is a computed fit to the individual data points.
You are to find from the data the orbital period, orbital radii and the masses of the
neutron star and the optical star. You need only give answers to two significant figures.
For this problem, use c = 3.0 ! 108 m " s #1 for the speed of light and
G = 6.7 ! 10"11 N # m 2 # kg "2 for the Universal Gravitation Constant.

Suggested procedure (you will need a ruler for the first three parts):

a. Determine the orbital period, T , in days.

b. Estimate the velocity, v0 of the optical star in its orbit around the center of mass.

c. Use T and v0 to find the radius of the orbit of the optical star, a0 , around the center of
mass.

d. Estimate the size of the orbit, aN , of the neutron star around the center of mass. First
express your answer in light ! seconds, then in meters.

e. Find the ratio of the mass of the neutron star, mN , to the mass of the optical star, mO .
[Hint: Recall that in any binary m1 r1 = m2 r2 .]

f. Use Kepler’s law to find the total mass of the binary system.

g. Find mN and mO ; express your answers in units of the mass of the Sun ( 2 ! 10 30 kg ).
Problem 7 Solutions:

a. Between the first and third zero-crossings on the upper graph, the distance is 8.98
cm. The scale of the time axis is 1.00 cm to 1 day, so the orbital period, P, is 8.98
days.
b. The Doppler shift velocities range from a maximum of 16.0 km!s-1 to a minimum
of -30.8 km!s-1. (The velocities are not symmetric around zero because the binary
system as a whole is moving relative to the solar system.) The optical star’s
orbital velocity is just half the difference between the maximum and minimum
Doppler velocities, or r0=23.4 km!s-1.
c. In one period P, moving at velocity v0, the optical star travels a distance

v0 p = 2! a0 . So the orbital radius, a0 is just

a0 =
1
2!
"
# ( )
( 9.98 )( 24 )( 60 )( 60 ) s %$ 23.4km & s'1 = 2.89 ( 10 6 km
d. The smallest time delay is -117 seconds, while the largest is 111 seconds. The
radius of the neutron star’s orbit (in light seconds) is half the difference between
these, or aN = 114lt ! sec. Since the speed of light is 2.998 x 1010 cm!s-1, the
distance aN is equal to 3.42 x 1012cm.
e. The ratio MN/M0 is just equal to a0/aN. That is

M N 2.89 ! 1011 cm
= = 8.45 ! 10 "2
M 0 3.42 ! 10 cm
12

f. The semimajor axis of the absolute orbit is a = aN + a0 = 3.71 ! 1012 cm and the
orbital period is 7.76 x 105s. Inserting these values into Newton’s form of
Kepler’s Third Law.

p2 a3
=
( 2! )2 G(MN + M0 )
we have,

M N + M 0 + 4! 2 ( 3.71 " 10 ) 12 3

= 5.02 " 10 34
g

( 7.76 " 10 )(5 2

6.67 " 10 ) #8

! M $ M
g. We use the fa(28.7.1)ct that M 0 # 1 + N & = M 0 + M N , and both N and
" M0 % M0
(M0 + MN ) are known. Then we have

5.02 ! 10 34 g
M0 = "2
= 4.63 ! 10 34 g
1 + 8.45 ! 10
The remaining mass in the mass of the neutron star, M N = 3.9 ! 10 33 g .
The mass of the sun is 1.99 x 1033 g, so the masses of the optical star and its
neutron star companion are 23.3 and 2.0 solar masses, respectively.
MIT OpenCourseWare
http://ocw.mit.edu

8.01SC Physics I: Classical Mechanics

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.