Processing math: 100%
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts
  • You do not have permission to view this page - please try signing in.

11.8: Inverse-square, two-body, central force

( \newcommand{\kernel}{\mathrm{null}\,}\)

The most important conservative, two-body, central interaction is the attractive inverse-square law force, which is encountered in both gravitational attraction and the Coulomb force. This force F(r) can be written in the form

F(r)=kr2ˆr

The force constant k is defined to be negative for an attractive force and positive for a repulsive force. In S.I. units the force constant k=Gm1m2 for the gravitational force and k=+q1q24πϵ0 for the Coulomb force. Note that this sign convention is the opposite of what is used in many books which use a negative sign in Equation ??? and assume k to be positive for an attractive force and negative for a repulsive force.

The conservative, inverse-square, two-body, central force is unique in that the underlying symmetries lead to four conservation laws, all of which are of pivotal importance in nature.

  1. Conservation of angular momentum: Like all conservative central forces, the inverse-square central two-body force conserves angular momentum as proven in chapter 11.3.
  2. Conservation of energy: This conservative central force can be represented in terms of a scalar potential energy U(r) as given by equation (11.1.2), where for this central force

    U(r)=kr

    Moreover, equation (11.6.3) showed that the center-of-mass Hamiltonian is conserved, that is, Hcm=Ecm

  3. Gauss’ Law: For a conservative, inverse-square, two-body, central force, the flux of the force field out of any closed surface is proportional to the algebraic sum of the sources and sinks of this field that are located inside the closed surface. The net flux is independent of the distribution of the sources and sinks inside the closed surface, as well as the size and shape of the closed surface. Chapter 2.14.5 proved this for the gravitational force field.
  4. Closed orbits: Two bodies interacting via the conservative, inverse-square, two-body, central force follow closed (degenerate) orbits as stated by Bertrand’s Theorem. The first consequence of this symmetry is that Kepler’s laws of planetary motion have stable, single-valued orbits. The second consequence of this symmetry is the conservation of the eccentricity vector defined in Equation ???.

Observables that depend on Gauss’s Law, or on closed planetary orbits, are extremely sensitive to addition of even a miniscule incremental exponent ξ to the radial dependence r(2±ξ) of the force. The statement that the inverse-square, two-body, central force leads to closed orbits can be proven by inserting Equation ??? into the orbit differential equation,

d2udψ2+u=μl21u2ku2=μkl2

Using the transformation

yu+μkl2

the orbit equation becomes

d2ydψ2+y=0

A solution of this equation is

y=Bcos(ψψ0)

Therefore

u=1r=μkl2[1+ϵcos(ψψ0)]

This is the equation of a conic section. For an attractive, inverse-square, central force, Equation ??? is the equation for an ellipse with the origin of r at one of the foci of the ellipse that has eccentricity ϵ, defined as ϵBl2μk

Equation ??? is the polar equation of a conic section. Equation ??? also can be derived with the origin at a focus by inserting the inverse square law potential into equation (11.6.10) which gives

ψ=±du2μEcml2+2μkl2uu2+ constant

The solution of this gives

u=1r=μkl2[1+1+2Ecml2μk2cos(ψψ0)]

Equations ??? and ??? are identical if the eccentricity ϵ equals

ϵ=1+2Ecml2μk2

The value of ψ0 merely determines the orientation of the major axis of the equivalent orbit. Without loss of generality, it is possible to assume that the angle ψ is measured with respect to the major axis of the orbit, that is ψ0=0. Then the equation can be written as

u=1r=μkl2[1+ϵcos(ψ)]=μkl2[1+1+2Ecml2μk2cos(ψ)]

This is the equation of a conic section where ϵ is the eccentricity of the conic section. The conic section is a hyperbola if ϵ>1, parabola if ϵ=1, ellipse if ϵ<1, and a circle if ϵ=0. All the equivalent one-body orbits for an attractive force have the origin of the force at a focus of the conic section. The orbits depend on whether the force is attractive or repulsive, on the conserved angular momentum l, and on the center-of-mass energy Ecm.

Bound orbits

Closed bound orbits occur only if the following requirements are satisfied.

  1. The force must be attractive, (k<0) then Equation ??? ensures that r is positive.
  2. For a closed elliptical orbit. the eccentricity ϵ<1 of the equivalent one-body representation of the orbit implies that the total center-of-mass energy Ecm<0, that is, the closed orbit is bound.

Bound elliptical orbits have the center-of-force at one interior focus F1 of the elliptical one-body representation of the orbit as shown in Figure 11.8.1.

9.8.1.PNG
Figure 11.8.1: Bound elliptical orbit.

The minimum value of the orbit r=rmin occurs when ψ=0, where

rmin=l2μk[1+ϵ]

This minimum distance is called the periapsis1.

The maximum distance, r=rmax, which is called the apoapsis, occurs when ψ=180o

rmax=l2μk[1ϵ]

Remember that since k<0 for bound orbits, the negative signs in equations ??? and ??? lead to r>0. The most bound orbit is a circle having ϵ=0 which implies that Ecm=μk2l2.

The shape of the elliptical orbit also can be described with respect to the center of the elliptical equivalent orbit by deriving the lengths of the semi-major axis a and the semi-minor axis b shown in Figure 11.8.1. a=12(rmin+rmax)=12(l2μk[1+ϵ]+l2μk[1ϵ])=l2μk[1ϵ2]b=a1ϵ2=l2μk[1ϵ2]

Remember that the predicted bound elliptical orbit corresponds to the equivalent one-body representation for the two-body motion as illustrated in Figure (11.2.2). This can be transformed to the individual spatial trajectories of the each of the two bodies in an inertial frame.

Kepler’s laws for bound planetary motion

Kepler’s three laws of motion apply to the motion of two bodies in a bound orbit due to the attractive gravitational force for which k=Gm1m2.

  1. Each planet moves in an elliptical orbit with the sun at one focus
  2. The radius vector, drawn from the sun to a planet, describes equal areas in equal times
  3. The square of the period of revolution about the sun is proportional to the cube of the major axis of the orbit.

Two bodies interacting via the gravitational force, which is a conservative, inverse-square, two-body central force, is best handled using the equivalent orbit representation. The first and second laws were proved in chapters 11.8 and 11.3. That is, the second law is equivalent to the statement that the angular momentum is conserved. The third law can be derived using the fact that the area of an ellipse is

A=πab=πa21ϵ2=πlμka32

Equations (11.3.7) and (11.3.8) give that the rate of change of area swept out by the radius vector is

dAdt=12r2˙ψ=l2μ

Therefore the period for one revolution τ is given by the time to sweep out one complete ellipse

τ=A(dAdt)=2π(μk)12a32

This leads to Kepler’s 3rd law τ2=4π2μka3

Bound orbits occur only for attractive forces for which the force constant k is negative, and thus cancel the negative sign in Equation ???. For example, for the gravitational force k=Gm1m2.

Note that the reduced mass μ=m1m2m1+m2 occurs in Kepler’s 3rd law. That is, Kepler’s third law can be written in terms of the actual masses of the bodies to be

τ2=4π2G(m1+m2)a3

In relating the relative periods of the different planets Kepler made the approximation that the mass of the planet m1 is negligible relative to the mass of the sun m2.

The eccentricity of the major planets ranges from ϵ=0.2056 for Mercury, to ϵ=0.0068 for Venus. The Earth has an eccentricity of ϵ=0.0167 with rmin=91106 miles and rmax=95106 miles. On the other hand, ϵ=0.967 for Halley’s comet, that is, the radius vector ranges from 0.6 to 18 times the radius of the orbit of the Earth.

The orbit energy can be derived by substituting the eccentricity, given by Equation ???, into the semi-major axis length a, given by Equation 11.8.15, which leads to the center-of-mass energy of

Ecm=k2a

However, the Hamiltonian, given by equation (11.6.3), implies that Ecm is

Ecm=12μv2+(kr)=k2a

For the simple case of a circular orbit, a=r then the velocity v equals v=kμr

For a circular orbit, the drag on a satellite lowers the total energy resulting in a decrease in the radius of the orbit and a concomitant increase in velocity. That is, when the orbit radius is decreased, part of the gain in potential energy accounts for the work done against the drag, and the remaining part goes towards increase of the kinetic energy. Also note that, as predicted by the Virial Theorem, the kinetic energy always is half the potential energy for the inverse square law force.

Unbound orbits

Attractive inverse-square central forces lead to hyperbolic orbits for ϵ>1 for which Ecm>0, that is, the orbit is unbound. In addition, the orbits always are unbound for a repulsive force since U=kr is positive as is the kinetic energy Tcm, thus Ecm=Tcm+Ucm>0. The radial orbit equation for either an attractive or a repulsive force is

r=l2μk[1+ϵcosψ]

For a repulsive force k is positive and l2 always is positive. Therefore to ensure that r remain positive the bracket term must be negative. That is

[1+ϵcosψ]<0k>0

For an attractive force k is negative and since l2 is positive then the bracket term must be positive to ensure that r is positive. That is, [1+ϵcosψ]>0k<0

9.8.2.PNG
Figure 11.8.2: Hyperbolic two-body orbits for a repulsive (left) and attractive (right) inverse-square, central two-body forces. Both orbits have the angular momentum vector pointing upwards out of the plane of the orbit

Figure 11.8.2 shows both branches of the hyperbola for a given angle ψ for the equivalent two-body orbits where the center of force is at the origin. For an attractive force, k<0, the center of force is at the interior focus of the hyperbola, whereas for a repulsive force the center of force is at the exterior focus. For a given value of |ψ| the asymptotes of the orbits both are displaced by the same impact parameter b from parallel lines passing through the center of force. The scattering angle, between the outgoing direction of the scattered body and the incident direction, is designated to be θ, which is related to the angle ψ by θ=1802ψ.

Eccentricity vector

Two-bodies interacting via a conservative two-body central force have two invariant first-order integrals, namely the conservation of energy and the conservation of angular momentum. For the special case of the inverse-square law, there is a third invariant of the motion, which Hamilton called the eccentricity vector2, that unambiguously defines the orientation and direction of the major axis of the elliptical orbit. It will be shown that the angular momentum plus the eccentricity vector completely define the plane and orientation of the orbit for a conservative inverse-square law central force.

Newton’s second law for a central force can be written in the form

˙p=f(r)ˆr

Note that the angular moment L=r×p is conserved for a central force, that is ˙L=0. Therefore the time derivative of the product p×L reduces to

ddt(p×L)=˙p×L=f(r)ˆr×(r×μ˙r)=f(r)μr[r(r˙r)r2˙r]

This can be simplified using the fact that

r˙r=12ddt(rr)=r˙r

thus

f(r)μr[r(r˙r)r2˙r]=μf(r)r2[˙rrr˙rr2]=μf(r)r2ddt(rr)

This allows Equation ??? to be reduced to

ddt(p×L)=μf(r)r2ddt(rr)

Assume the special case of the inverse-square law, Equation ???, then the central force Equation ??? reduces to

ddt(p×L)=ddt(μkˆr)

or

ddt[(p×L)+(μkˆr)]=0

Define the eccentricity vector A as

A(p×L)+(μkˆr)

then Equation ??? corresponds to

dAdt=0

This is a statement that the eccentricity vector A is a constant of motion for an inverse-square, central force.

The definition of the eccentricity vector A and angular momentum vector L implies a zero scalar product,

AL=0

Thus the eccentricity vector A and angular momentum L are mutually perpendicular, that is, A is in the plane of the orbit while L is perpendicular to the plane of the orbit. The eccentricity vector A, always points along the major axis of the ellipse from the focus to the periapsis as illustrated on the left side in Figure 11.8.3. As a consequence, the two orthogonal vectors A and L completely define the plane of the orbit, plus the orientation of the major axis of the Kepler orbit, in this plane. The three vectors A, p×L, and (μkˆr) obey the triangle rule as illustrated in the left side of Figure 11.8.3.

9.8.3.PNG
Figure 11.8.3: The elliptical trajectory and eccentricity vector A for two bodies interacting via the inversesquare, central force for eccentricity ϵ=0.75. The left plot shows the elliptical spatial trajectory where the semi-major axis is assumed to be on the x-axis and the angular momentum L=lˆz, is out of the page. The force centre is at one foci of the ellipse. The vector coupling relation A(p×L)+(μkˆr) is illustrated at four points on the spatial trajectory. The right plot is a hodograph of the linear momentum p for this trajectory. The periapsis is denoted by the number 1 and the apoapsis is marked as 3 on both plots. Note that the eccentricity vector A is a constant that points parallel to the major axis towards the perapsis.

Hamilton noted the direct connection between the eccentricity vector A and the eccentricity ϵ of the conic section orbit. This can be shown by considering the scalar product

Ar=Arcosψ=r(p×L)+μkr

Note that the triple scalar product can be permuted to give

r(p×L)=(r×p)L=LL=l2

Inserting Equation ??? into ??? gives 1r=μkl2(1Aμkcosψ)

Note that equations ??? and ??? are identical if ψ0=0. This implies that the eccentricity ϵ and A are related by

ϵ=Aμk

where k is defined to be negative for an attractive force. The relation between the eccentricity and total center-of-mass energy can be used to rewrite Equation ??? in the form A2=μ2k2+2μEcml2

The combination of the eccentricity vector A and the angular momentum vector L completely specifies the orbit for an inverse square-law central force. The trajectory is in the plane perpendicular to the angular momentum vector L, while the eccentricity, plus the orientation of the orbit, both are defined by the eccentricity vector A. The eccentricity vector and angular momentum vector each have three independent coordinates, that is, these two vector invariants provide six constraints, while the scalar invariant energy E, adds one additional constraint. The exact location of the particle moving along the trajectory is not defined and thus there are only five independent coordinates governed by the above seven constraints. Thus the eccentricity vector, angular momentum, and center-of-mass energy are related by the two equations ??? and ???.

Noether’s theorem states that each conservation law is a manifestation of an underlying symmetry. Identification of the underlying symmetry responsible for the conservation of the eccentricity vector A is elucidated using Equation ??? to give

(μkˆr)=A(p×L) Take the scalar product

(μkˆr)(μkˆr)=(μk)2=p2L2+A22L(p×L)

Choose the angular momentum to be along the z-axis, that is, L=lˆz, and, since p and A are perpendicular to L, then p and A are in the ˆxˆy plane. Assume that the semimajor axis of the elliptical orbit is along the x-axis, then the locus of the momentum vector on a momentum hodograph has the equation

p2x+(pyAL)2=(μkL)2

Equation ??? implies that the locus of the momentum vector is a circle of radius |μkL| with the center displaced from the origin at coordinates (0,AL) as shown by the momentum hodograph on the right side of an Figure 11.8.3. The angle β and eccentricity ϵ are related by,

cosβ=A/Lμk/L=Aμk=ϵ

The circular orbit is centered at the origin for ϵ=Aμk=0, and thus the magnitude |p| is a constant around the whole trajectory.

The inverse-square, central, two-body, force is unusual in that it leads to stable closed bound orbits because the radial and angular frequencies are degenerate, i.e. ωr=ωψ. In momentum space, the locus of the linear momentum vector p is a perfect circle which is the underlying symmetry responsible for both the fact that the orbits are closed, and the invariance of the eccentricity vector. Mathematically this symmetry for the Kepler problem corresponds to the body moving freely on the boundary of a four-dimensional sphere in space and momentum. The invariance of the eccentricity vector is a manifestation of the special property of the inverse-square, central force under certain rotations in this four-dimensional space; this O(4) symmetry is an example of a hidden symmetry.


1The greek term apsis refers to the points of greatest or least distance of approach for an orbiting body from one of the foci of the elliptical orbit. The term periapsis or pericenter both are used to designate the closest distance of approach, while apoapsis or apocenter are used to designate the farthest distance of approach. Attaching the terms "perí-" and "apo-" to the general term "-apsis" is preferred over having different names for each object in the solar system. For example, frequently used terms are "-helion" for orbits of the sun, "-gee" for orbits around the earth, and "-cynthion" for orbits around the moon.

2The symmetry underlying the eccentricity vector is less intuitive than the energy or angular momentum invariants leading to it being discovered independently several times during the past three centuries. Jakob Hermann was the first to indentify this invariant for the special case of the inverse-square central force. Bernoulli generalized his proof in 1710. Laplace derived the invariant at the end of the 18th century using analytical mechanics. Hamilton derived the connection between the invariant and the orbit eccentricity. Gibbs derived the invariant using vector analysis. Runge published the Gibb’s derivation in his textbook which was referenced by Lenz in a 1924 paper on the quantal model of the hydrogen atom. Goldstein named this invariant the "Laplace-Runge-Lenz vector", while others have named it the "Runge-Lenz vector" or the "Lenz vector". This book uses Hamilton’s more intuitive name of "eccentricity vector".


This page titled 11.8: Inverse-square, two-body, central force is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Douglas Cline via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?