Physics Library
 An open source physics library
Encyclopedia | Forums | Docs | Random  

Effective Potential and Radial Motion in Central Gravity

(Topic)

Effective Potential and Radial Motion in Central Gravity

The gravitational two-body problem is a central force problem.

Once angular momentum is conserved, the planar motion can be separated into:

  • radial motion, which changes the separation r,
  • angular motion, which changes the polar angle 𝜃.

The effective-potential method absorbs the angular part of the kinetic energy into an effective radial potential.

This reduces the orbital problem to a one-dimensional energy equation for

r (t).
(1)

The method is useful for:

  • locating periapsis and apoapsis,
  • identifying allowed and forbidden radial regions,
  • understanding the angular-momentum barrier,
  • finding circular orbits,
  • testing circular-orbit stability,
  • distinguishing radial from nonradial motion.

1 Velocity in plane polar coordinates

For planar motion,

|----------------|
|v = ˙r^er + r𝜃˙^e𝜃.|
------------------
(2)

Therefore

|--------------|
v2-=-r˙2-+-r2𝜃˙2.-
(3)

The kinetic energy is

|----------------------|
|     1-  2   1-   2 ˙2|
|T =  2m r˙ + 2 mr 𝜃  .|
-----------------------
(4)

PIC

Figure 1. The velocity separates into radial and transverse components. The effective-potential method uses angular-momentum conservation to replace the transverse kinetic-energy term by a function of radius.

2 Angular momentum in a central force

For a central force,

τ =  r × F = 0.
(5)

Therefore angular momentum is conserved.

For planar motion,

|-------2-˙------------|
-L-=-mr--𝜃-=-constant.-|
(6)

Define the specific angular momentum

|--------------|
|h = -L =  r2 ˙𝜃.
-----m---------|
(7)

Thus

|--------|
|˙   -h  |
|𝜃 = r2 .|
---------
(8)

3 Energy with angular momentum

The gravitational mechanical energy is

     1-   2   1-  2 ˙2  GM--m--
E  = 2 m ˙r +  2mr  𝜃 −    r   .
(9)

Use

       2 ˙
L = mr  𝜃.
(10)

Then

         L2
r2𝜃˙2 = --2-2.
       m  r
(11)

Therefore

|----------------------------|
|    1    2    L2     GM  m  |
E  = --m ˙r +  ----2 − ------.|
-----2--------2mr-------r-----
(12)

The entire radial problem is now contained in one equation.

4 The effective potential

Define

|------------2------------|
U   (r) = -L---−  GM--m-. |
--eff-------2mr2-------r-----
(13)

Then

|--------------------|
|     1-  2          |
|E =  2m ˙r +  Ueff(r).|
----------------------
(14)

The effective potential contains two terms:

|----------2--|
U    =  -L----|
--ang---2mr2--|
(15)

and

|--------------|
|       GM--m--|
Ug =  −   r   .|
----------------
(16)

The first term is positive and diverges as

r → 0
(17)

when

L  ⁄= 0.
(18)

The second term is the ordinary gravitational potential energy.

PIC

Figure 2. The effective potential is the sum of a positive angular-momentum term proportional to 1∕r2 and the negative gravitational term proportional to −1∕r. Their competition produces a minimum.

5 Specific effective potential

Divide the energy equation by

m.
(19)

Using

μ =  GM
(20)

and

h = L ∕m,
(21)

we obtain

|------------------|
|    1- 2          |
|𝜀 = 2 ˙r + Φe ff(r),|
-------------------
(22)

where

|-----------2------|
|Φeff(r) = -h--− μ-.|
----------2r2----r--
(23)

The units of

Φeff
(24)

are

m2 ∕s2,
(25)

the same as specific orbital energy.

6 The angular-momentum barrier

When

h ⁄= 0,
(26)

 h2
--2-→  + ∞
2r
(27)

as

r →  0.
(28)

Thus

|----------------------------------|
-Φeff(r)-→--+∞------(r-→--0,-h ⁄=-0).|
(29)

This is often called the centrifugal barrier or angular-momentum barrier.

It is not an additional physical gravitational interaction.

It arises because conservation of angular momentum forces the transverse speed

          h
v𝜃 = r˙𝜃 = --
          r
(30)

to grow as the radius decreases.

Consequently the transverse kinetic energy grows as

1∕r2.
(31)

7 Radial kinetic energy

From

𝜀 =  1˙r2 + Φeff,
     2
(32)

the radial kinetic energy per unit mass is

|------------------|
|1                 |
|--˙r2 = 𝜀 − Φe ff(r).
-2-----------------
(33)

Therefore

|-------------------|
r˙2 = 2 [𝜀 − Φeff (r )]. |
---------------------
(34)

Since

r˙2 ≥ 0,
(35)

the allowed radial region must satisfy

|------------|
|𝜀 ≥ Φe ff(r).|
-------------
(36)

Regions where

𝜀 < Φeff
(37)

are forbidden by energy conservation.

PIC

Figure 3. Radial motion is allowed only where the horizontal specific-energy level lies on or above the effective-potential curve. Intersections of the two are radial turning points.

8 Radial turning points

A radial turning point occurs when

˙r = 0.
(38)

Therefore

-------------
|            |
-𝜀 =-Φeff(rt).|
(39)

Substituting,

     h2    μ
𝜀 = ---2−  --.
    2rt    rt
(40)

Multiply by

   2
2r t :
(41)

|--2-----------2------|
2𝜀rt-+-2μrt-−-h--=-0.--
(42)

For

𝜀 ⁄= 0,
(43)

the roots are

|-----------∘------------|
|     −-μ ±---μ2-+-2𝜀h2- |
|rt =         2𝜀        .|
-------------------------
(44)

Only positive real roots are physical radii.

9 Bound nonradial motion

For

h ⁄= 0
(45)

and negative energy above the effective-potential minimum, there are two positive turning radii:

|r-<--r-.|
--p----a-|
(46)

The smaller radius

rp
(47)

is periapsis.

The larger radius

ra
(48)

is apoapsis.

The radial coordinate oscillates between these values:

|--------------|
|r ≤ r (t) ≤ r .|
--p----------a--
(49)

At each radial turning point,

˙r = 0,
(50)

but the transverse speed is generally nonzero.

PIC

Figure 4. A negative-energy level above the effective-potential minimum intersects the curve twice. These intersections give periapsis and apoapsis for bound nonradial motion.

10 Zero-energy turning point

For

𝜀 = 0,
(51)

the turning-point equation becomes

        2
2μr − h  =  0.
(52)

Thus the finite radial turning point is

|--------|
|     h2-|
rp =  2μ.|
----------
(53)

The other radial limit is

r → ∞.
(54)

This is the parabolic escape threshold.

11 Positive-energy motion

For

𝜀 > 0,
(55)

there is one positive finite turning point.

It is the periapsis of an unbound trajectory.

As

r → ∞,
(56)

Φe ff → 0,
(57)

so

1
-r˙2
2
(58)

approaches a nonzero value after the transverse speed has fallen toward zero.

The complete orbit is hyperbolic in the ideal inverse-square problem.

12 Circular orbit as an effective-potential equilibrium

A circular orbit has constant radius:

r˙= 0,     ¨r = 0.
(59)

Therefore it occurs at an extremum of the effective potential:

|----------|
|dΦeff      |
|-----=  0.|
--dr--------
(60)

Differentiate:

dΦeff-     h2-  μ-
 dr  =  − r3 + r2.
(61)

Set this equal to zero:

  h2    μ
− --3+  -2 = 0.
   rc   rc
(62)

Hence

|------2-|
|rc = h-.|
------μ---
(63)

In terms of total angular momentum,

|--------------|
|       L2     |
|rc = ------2. |
------GM--m----
(64)

13 Circular-orbit speed from the effective potential

Since

h =  rcv𝜃,
(65)

the circular condition

h2 =  μrc
(66)

gives

 2     h2-   μ-
vcirc = r2c =  rc.
(67)

Therefore

|------∘-----|
|vcirc =   μ-.|
|         rc |
--------------
(68)

The effective-potential method therefore reproduces the circular-orbit condition derived directly from centripetal acceleration.

14 Energy of the circular orbit

Evaluate the effective potential at

rc = h2 ∕μ.
(69)

Using

h2 = μrc,
(70)

Φeff(rc) = μrc-
2r2
  c −μ-
r
 c (71)
= − μ
---
2rc. (72)

Therefore

|--------------------|
|                 2  |
|𝜀c = − -μ- = − -μ--.|
--------2rc-----2h2--
(73)

This is the minimum possible energy for a fixed nonzero

h.
(74)

PIC

Figure 5. For fixed nonzero angular momentum, the circular orbit occurs at the minimum of the effective potential. Energies below this minimum are impossible for that value of angular momentum.

15 Stability of the circular orbit

An extremum is stable if the effective potential has a local minimum.

Differentiate again:

d2Φe-ff   3h2-   2μ-
 dr2   =  r4 −  r3 .
(75)

At

r = rc
(76)

with

 2
h  = μrc,
(77)

 2    |
d-Φeff-||
 dr2  |r c = 3μ-
r3c −2μ-
r3c (78)
= μ
-3
rc > 0. (79)

Thus the Newtonian circular orbit is stable against sufficiently small radial perturbations at fixed angular momentum.

16 Small radial oscillations about a circular orbit

Let

r = rc + ξ,    |ξ| ≪ rc.
(80)

Near the minimum,

                            |
                    1-d2Φeff-||  2
Φe ff(r) ≈ Φeff(rc) + 2  dr2  | ξ .
                            rc
(81)

The radial motion is approximately harmonic:

|-----2------|
-¨ξ +-ω-rξ-=-0,
(82)

with

|--------|
| 2   μ  |
ω r = -3.|
------rc--
(83)

The circular orbital angular frequency also satisfies

|--------|
ω2 =  μ-.|
--c---r3c--
(84)

Therefore

|--------|
-ωr-=-ωc-|
(85)

for Newtonian inverse-square gravity.

This equality is closely related to the closed nature of Keplerian ellipses.

17 Radial equation of motion

The effective potential also gives the radial force equation.

Since

E =  1m ˙r2 + Ueff(r),
     2
(86)

the effective radial force is

        dUeff
m ¨r = − -----.
         dr
(87)

Differentiate:

            2
dUeff-=  − L---+  GM--m-.
 dr       mr3      r2
(88)

Therefore

|--------------------|
|      L2     GM  m  |
m ¨r =  ---3 − ---2--.|
-------mr-------r-----
(89)

Divide by

m  :
(90)

|------------|
|    h2   μ  |
¨r =  -3-− -2.|
-----r----r---
(91)

18 Relation to the polar-coordinate acceleration equation

The radial acceleration in polar coordinates is

           ˙2
ar = ¨r − r𝜃 .
(92)

Newton’s radial equation is

  (       )
m  ¨r − r𝜃˙2  =  − GM--m-.
                   r2
(93)

Using

 ˙2   h2-
r𝜃  = r3 ,
(94)

we recover

|------------|
|    h2   μ  |
¨r =  -3-− -2.|
-----r----r---
(95)

Thus the effective-potential equation is not introducing a new force law.

It is a reformulation of the original Newtonian equations using the conserved angular momentum.

19 Changing angular momentum changes the barrier

For fixed central mass, the effective potential is

        2
Φ   =  h---− μ-.
  eff    2r2   r
(96)

Increasing

h
(97)

raises the small-radius barrier.

It also moves the circular-orbit minimum outward:

|------2-|
|rc = h-.|
------μ---
(98)

Thus larger angular momentum makes close approach more difficult for a trajectory with the same energy.

PIC

Figure 6. Increasing specific angular momentum raises the small-radius barrier and shifts the effective-potential minimum to a larger circular-orbit radius.

20 Special case: zero angular momentum

If

h = 0,
(99)

then

Φe ff = − μ.
         r
(100)

There is no angular-momentum barrier.

The energy equation becomes

|-------------|
|   1      μ  |
𝜀 = --˙r2 − -. |
----2------r---
(101)

Pure radial infall can therefore reach

r = 0
(102)

in the ideal point-mass model.

The singularity reflects the idealized Newtonian point source.

For a physical extended body, the interior gravitational field must replace the point-mass field.

21 Radial versus nonradial turning points

For

h = 0,
(103)

a turning point with

˙r = 0
(104)

also has total speed

v = 0.
(105)

For

h ⁄= 0,
(106)

a radial turning point has

˙r = 0
(107)

but

v =  h-⁄=  0.
 𝜃   r
(108)

Therefore

|------------|
|          h |
|vturning = r- |
-----------t-
(109)

for nonradial motion.

This is why periapsis and apoapsis are not points where the orbiting body stops.

PIC

Figure 8. A radial turning point means ṙ = 0. With zero angular momentum the total speed is zero, but with nonzero angular momentum the transverse speed remains finite.

22 Example 1: Earth orbit from a turning point

Take Earth’s gravitational parameter as

μE = 3.9860 × 1014 m3 ∕s2.
(110)

Suppose at

r0 = 7.00 × 106m
(111)

the velocity is purely tangential with magnitude

v0 = 8.00 × 103 m ∕s.
(112)

Because

˙r = 0,
(113)

this radius is a radial turning point.

The specific angular momentum is

h = r0v0 (114)
= 5.60 × 1010 m2∕s. (115)

The specific energy is

𝜀 = 1
--
2v02 −μE
---
 r0 (116)
≈−2.494 × 107 J∕kg. (117)

Since the energy is negative, the orbit is bound.

23 Example 1 continued: the second turning point

The turning-point equation is

2𝜀r2 + 2μEr − h2 =  0.
(118)

One root is the known turning radius

rp = 7.00 × 106 m.
(119)

The second positive root is

|------------------|
|ra ≈ 8.98 × 106 m. |
--------------------
(120)

Thus the radial coordinate oscillates between approximately

7000 km
(121)

and

8980 km.
(122)

Because the tangential speed at

7000 km
(123)

is larger than the circular speed there, that turning point is the periapsis.

PIC

Figure 7. In the numerical Earth example, the specific-energy level intersects the effective potential at two radii. The smaller is periapsis and the larger is apoapsis.

24 Example 2: circular orbit for a specified angular momentum

Suppose

              10   2
h =  5.00 × 10   m  ∕s
(124)

about Earth.

Then

rc = h2-
μ
 E (125)
≈ 6.27 × 106 m. (126)

The circular specific energy is

𝜀c = −μE
---
2rc (127)
≈−3.18 × 107 J∕kg. (128)

For that fixed

h,
(129)

no motion is possible with

𝜀 < 𝜀c.
(130)

25 Dimensionless form

It is useful to normalize radius by the circular-orbit radius

rc = h2 ∕μ.
(131)

Define

     r
x =  --.
     rc
(132)

Since

h2 = μrc,
(133)

the effective potential becomes

Φeff =  μ
---
2rc1
-2-
x −μ
--
rc 1
--
x. (134)

Thus

|----------------|
|Φeff     1     1 |
|-----= ---2 − -.|
-μ∕rc---2x-----x--
(135)

The minimum occurs at

x = 1
(136)

with value

− 1∕2.
(137)

This universal curve describes all Newtonian gravitational effective potentials with nonzero angular momentum after suitable scaling.

26 Discriminant and existence of turning points

The quadratic turning-point equation has discriminant

-----------------
|     2       2 |
D--=-μ--+-2𝜀h-.--
(138)

Real turning points require

D  ≥ 0.
(139)

Therefore

|-------2--|
𝜀 ≥  − μ--.|
-------2h2--
(140)

But

  μ2
− ----
  2h2
(141)

is exactly the minimum circular-orbit energy.

Thus the algebraic condition for real turning points is the same as the effective-potential requirement that the energy lie at or above its minimum.

27 Three energy regimes at fixed nonzero angular momentum

For fixed

h ⁄=  0 :
(142)

  • At the minimum energy, the two turning points coincide and the orbit is circular.
  • For negative energy above the minimum, there are two finite turning radii and the motion is bound.
  • For zero energy, there is one finite periapsis and the outer limit is infinity.
  • For positive energy, there is one finite periapsis and unbound outward motion is allowed.

The effective-potential graph displays all four cases without solving the full orbit equation.

28 What the effective potential does not determine

The effective potential determines the radial dynamics for a given

E
(143)

and

L.
(144)

It does not by itself give

r(𝜃)
(145)

or the geometric shape of the orbit.

To obtain the trajectory in the orbital plane, one must combine:

  • angular-momentum conservation,
  • the radial dynamics,
  • the relation between r and 𝜃.

The next articles develop that connection.

29 Common mistakes

  1. Treating the angular-momentum term as a new physical gravitational force.
  2. Forgetting that Ueff is used for radial motion at fixed angular momentum.
  3. Confusing total angular momentum L with specific angular momentum h = L∕m.
  4. Using h2∕(2r) instead of h2∕(2r2).
  5. Forgetting the negative sign of the gravitational term.
  6. Allowing motion in a region where 𝜀 < Φeff.
  7. Assuming a radial turning point means total speed is zero.
  8. Forgetting that periapsis and apoapsis have nonzero tangential speed when h≠0.
  9. Using the quadratic turning-point formula without checking for positive real roots.
  10. Forgetting that a circular orbit occurs at an extremum of the effective potential.
  11. Finding the extremum but failing to test stability with the second derivative.
  12. Assuming the effective potential has a small-radius barrier when h = 0.
  13. Using the radial-collision result for a physical extended body all the way to r = 0.
  14. Confusing the circular-orbit radius h2∕μ with the turning radius for arbitrary energy.
  15. Assuming the effective potential alone gives the full orbit shape r(𝜃).

30 Practice exercises

  1. Derive the polar-coordinate kinetic energy for planar motion.
  2. Starting from angular-momentum conservation, derive the effective potential.
  3. Derive the specific effective potential in terms of h and μ.
  4. Show that radial motion is allowed only where 𝜀 ≥ Φeff.
  5. Derive the quadratic equation for radial turning points.
  6. Analyze the number of positive turning points for negative, zero, and positive energy.
  7. Derive the circular-orbit radius rc = h2∕μ.
  8. Derive the circular-orbit energy 𝜀c = −μ2∕(2h2).
  9. Show that the circular orbit is stable by evaluating the second derivative.
  10. Derive the small-radial-oscillation frequency about a circular orbit.
  11. Recover the radial equation r = h2∕r3 − μ∕r2 from U eff.
  12. Explain why no angular-momentum barrier exists for h = 0.
  13. Reproduce the Earth turning-point example with r0 = 7000 km and v0 = 8.00 km/s.
  14. Reproduce the circular-orbit example for h = 5.00 × 1010 m2∕s.
  15. Derive the dimensionless effective-potential curve 1∕(2x2) − 1∕x.

31 Summary

For planar motion in Newtonian central gravity,

|----------------------|
|       2 ˙            |
-L-=-mr--𝜃-=-constant.-
(146)

The conserved energy can be written

----------------------
|     1              |
|E =  -m ˙r2 + Ueff(r),|
------2---------------
(147)

with

|-------------------------|
|         -L2---  GM--m-- |
Ueff (r) = 2mr2 −     r  . |
---------------------------
(148)

Per unit mass,

|--------------------|
|            2       |
|𝜀 = 1-˙r2 + h--−  μ-.|
-----2------2r2---r--
(149)

Allowed radial motion requires

|------------|
|𝜀 ≥ Φe ff(r).|
-------------
(150)

Radial turning points satisfy

|--------------------|
|2𝜀r2 + 2 μr − h2 = 0.|
----------------------
(151)

The circular orbit occurs at the effective-potential minimum:

|--------|
|     h2 |
|rc = μ-,|
----------
(152)

with

|----------|
|      -μ- |
𝜀c = − 2rc.|
------------
(153)

For Newtonian inverse-square gravity,

---------------------
| 2   ||             |
|d-Φeff|  =  μ- > 0, |
| dr2 |rc   r3c      |
---------------------
(154)

so the circular orbit is radially stable.

The effective-potential method converts a two-dimensional central force problem into a one-dimensional radial-energy problem while preserving the influence of angular momentum.

References

References

[1]   PhysicsLibrary, M07-05, Orbital Energy, Escape Speed, and Bound versus Unbound Motion.

[2]   PhysicsLibrary, M04-08, Conservation of Angular Momentum and Central Force Motion.

[3]   PhysicsLibrary, M01-10, Polar-Coordinate Particle Kinematics.

[4]   J. R. Taylor, Classical Mechanics, University Science Books, 2005.

[5]   D. Kleppner and R. Kolenkow, An Introduction to Mechanics, 2nd ed., Cambridge University Press, 2014.


"Effective Potential and Radial Motion in Central Gravity" is owned by bloftin.
(view preamble)
View style:
Other names:  M07-06
Also defines:  effective potential, specific effective potential, centrifugal barrier, radial kinetic energy, radial turning point, circular-orbit equilibrium, stable circular orbit
Keywords:  effective potential, central gravity, radial motion, angular momentum, turning point, centrifugal barrier, circular orbit, orbital stability, periapsis, apoapsis

Attachments:
GRE Physics Companion: Effective Potential and Radial Motion in Central Gravity (Example) by bloftin

Cross-references: two-dimensional, frequency, quadratic equation, formula, relation, orbit equation, graph, algebraic, circular-orbit energy, magnitude, parameter, gravitational field, radial equation, polar coordinates, acceleration, angular frequency, harmonic, centripetal acceleration, total angular momentum, turning point, energy conservation, mass, speed, forces, conservation of angular momentum, specific orbital energy, units, gravitational potential energy, mechanical energy, specific angular momentum, function, velocity, testing, energy, kinetic energy, angular momentum, central force
There are 8 references to this object.

This is version 1 of Effective Potential and Radial Motion in Central Gravity, born on 2026-10-31.
Object id is 1559, canonical name is EffectivePotentialAndRadialMotionInCentralGravity.
Accessed 7 times total.

Classification:
Physics Classification: 45.20.Dd (Newtonian mechanics)
 45.50.-j (Dynamics and kinematics of a particle and a system of particles)

Pending Errata and Addenda

None.

Discussion

No messages.

Interact