Stability of Circular Orbits

Every satellite you rely on — GPS, weather, the ISS — lives on a nearly circular orbit that has survived years of small kicks: thruster burns, drag, lunar tugs. That survival is not automatic. A circular orbit exists wherever gravity exactly supplies the needed centripetal force, but existing and persisting are different questions. Nudge the radius slightly: does the orbit breathe gently around the circle, or does the perturbation grow until the body spirals away or crashes inward? Classical mechanics answers with one clean idea. Conservation of angular momentum reduces the motion to a single radial coordinate rolling in the effective potential Veff(r)=V(r)+L2/2μr2V_{\mathrm{eff}}(r)=V(r)+L^2/2\mu r^2, and a circular orbit is simply a radius where that landscape is flat. Stability is the familiar valley-versus-hilltop test: if Veff′′(r0)>0V_{\mathrm{eff}}''(r_0)>0 the orbit sits in a valley and small perturbations oscillate harmlessly; negative curvature and they run away. In force language the criterion reads F′(r0)/F(r0)+3/r0>0F'(r_0)/F(r_0)+3/r_0>0 — for power laws F=−k rnF=-k\,r^n, simply n>−3n>-3. Gravity (n=−2n=-2) and the spring force (n=+1n=+1) pass with room to spare; any attraction steeper than inverse-cube confines nothing. The same curvature also sets the frequency ωr\omega_r of the radial breathing, and comparing ωr\omega_r to the orbital frequency opens the door to apsidal precession and Bertrand's theorem. Push the analysis into general relativity and the valley itself disappears below a critical radius — the innermost stable circular orbit, the inner edge of black-hole accretion disks.

Key takeaways
  • A circular orbit exists wherever Veff(r)=V(r)+L2/2μr2V_{\mathrm{eff}}(r)=V(r)+L^2/2\mu r^2 is stationary: F(r0)=−L2/μr03F(r_0)=-L^2/\mu r_0^3.
  • It is stable only if that stationary point is a minimum: Veff′′(r0)>0V_{\mathrm{eff}}''(r_0)>0, equivalently F′(r0)F(r0)+3r0>0\dfrac{F'(r_0)}{F(r_0)}+\dfrac{3}{r_0}>0.
  • For power-law forces F=−k rnF=-k\,r^n the criterion is n>−3n>-3; the frequency ratio is ωr2/ωφ2=3+n\omega_r^2/\omega_\varphi^2=3+n (Kepler n=−2n=-2: ratio 1; Hooke n=+1n=+1: ratio 4).
  • Stable perturbations oscillate at ωr=Veff′′(r0)/μ\omega_r=\sqrt{V_{\mathrm{eff}}''(r_0)/\mu}; the ratio ωr/ωφ\omega_r/\omega_\varphi controls precession and closure (Bertrand).
  • In general relativity the effective potential loses its minimum below r=6GM/c2r=6GM/c^2 — the ISCO, the inner edge of accretion disks.
  • Key idea: a circular orbit is the equilibrium of the radial coordinate in VeffV_{\mathrm{eff}}; stability means that equilibrium is a genuine minimum, which for the force reads F′/F+3/r0>0F'/F+3/r_0>0 — power laws need n>−3n>-3.
The core insight

A circular orbit is the equilibrium of the radial degree of freedom in the effective potential Veff(r)=V(r)+L2/(2μr2)V_{\mathrm{eff}}(r)=V(r)+L^2/(2\mu r^2). Stability ⇔ that equilibrium is a minimum ⇔ Veff′′(r0)>0V_{\mathrm{eff}}''(r_0)>0, which for the force translates into F′(r0)/F(r0)+3/r0>0F'(r_0)/F(r_0)+3/r_0>0.

Where this lives today

The ISCO is a working number in black-hole astronomy. The thermal and iron-line X-ray spectra of accreting black holes are fitted with disk models truncated at the innermost stable circular orbit; because a spinning (Kerr) hole drags its ISCO inward from 6GM/c26GM/c^2 toward GM/c2GM/c^2, the fitted inner radius is a standard way of measuring black-hole spin across X-ray binaries and active galactic nuclei. Gravitational-wave astronomy hears the same boundary: an inspiral signal sweeps up in frequency until the orbit reaches the ISCO region and transitions to plunge — the 440 Hz of the worked example is why LIGO's sweet spot sits at hundreds of hertz for stellar-mass pairs. Closer to home, the criterion is bread-and-butter in accelerator and trap design: radial confinement must stay shallower than 1/r31/r^3, and the frequency ratio ωr/ωφ\omega_r/\omega_\varphi — the "tune" — is engineered away from low-order rationals to avoid resonant kicks.

The physical picture

Fix the angular momentum and the problem collapses to one coordinate: the radius, sliding in the landscape Veff(r)V_{\mathrm{eff}}(r). The centrifugal term walls off small radii; the attraction pulls the landscape down. Whether their sum has a genuine valley — and whether your orbit sits in it — is the entire question:

Figure
Valley or hilltop — the same circular-orbit condition sits at a minimum for gentle forces and a maximum for steep ones
gentle force, n greater than −3 genuine valley: nudged orbit oscillates, radius breathes at frequency ω r — STABLE steep force, n less than −3 the extremum is a hilltop: any nudge grows — spiral to the centre or escape. UNSTABLE radius radius

Why is −3-3 the dividing exponent? The centrifugal barrier contributes an outward force scaling as 1/r31/r^3. An attraction shallower than that loses the tug-of-war at small radius — the barrier wins and the body bounces back out, so a valley forms. An attraction steeper than 1/r31/r^3 overwhelms the barrier at small rr while fading faster at large rr: the curves cross at a hilltop instead, and the balance point repels rather than restores. Once the valley exists, the physics near its floor is a harmonic oscillator: the radius breathes at ωr\omega_r set by the curvature, while the body revolves at ωφ\omega_\varphi. The two clocks generally disagree, and their ratio — 3+n\sqrt{3+n} for a power law — is the bridge from stability to the apsidal-precession and Bertrand's-theorem articles.

Circular orbits: existence

With LL conserved, the radial coordinate obeys Newton's law in the effective potential:

μr¨=F(r)+L2μr3=−dVeffdr,Veff(r)=V(r)+L22μr2(1)\mu\ddot r = F(r)+\frac{L^2}{\mu r^3} = -\frac{\mathrm{d}V_{\mathrm{eff}}}{\mathrm{d}r},\qquad V_{\mathrm{eff}}(r)=V(r)+\frac{L^2}{2\mu r^2} \qquad (1)

where: μ\mu — the reduced mass; F(r)=−dV/drF(r)=-\mathrm{d}V/\mathrm{d}r — the central force, negative when attractive; L2/μr3L^2/\mu r^3 — the outward centrifugal force. The angular motion has been traded for the centrifugal term, so bead-in-landscape intuition now applies literally. A circular orbit is a radius that stays constant forever. Setting the right side of (1) to zero gives the existence condition:

F(r0)=−L2μr03⟺Veff′(r0)=0(2)F(r_0)=-\frac{L^2}{\mu r_0^3}\qquad\Longleftrightarrow\qquad V_{\mathrm{eff}}'(r_0)=0 \qquad (2)

where: r0r_0 — the circular-orbit radius. The right side is strictly negative for L≠0L\neq0, so the force must be attractive at r0r_0 — a repulsive force has no circular orbits. For gravity, (2) is the familiar μv2/r0=GMμ/r02\mu v^2/r_0=GM\mu/r_0^2; existence holds at every radius for any attractive power law — existence is cheap, stability is not.

The stability criterion

Whether the stationary point is a valley or a hilltop is second-derivative business. Derivation A linearizes the radial equation and shows a perturbation η=r−r0\eta=r-r_0 obeys η¨+ωr2η=0\ddot\eta+\omega_r^2\eta=0 with ωr2=Veff′′(r0)/μ\omega_r^2=V_{\mathrm{eff}}''(r_0)/\mu, so everything hinges on the sign of the curvature:

Veff′′(r0)>0⟺F′(r0)F(r0)+3r0>0(3)V_{\mathrm{eff}}''(r_0)>0 \qquad\Longleftrightarrow\qquad \frac{F'(r_0)}{F(r_0)}+\frac{3}{r_0}>0 \qquad (3)

where: F′F' — the radial derivative of the force; the second form follows by differentiating (1) and eliminating LL with (2) — Derivation A, Step 4. Read it as a competition of logarithmic slopes: the force may fall with radius, but not faster than the centrifugal force's own 1/r31/r^3. If the curvature is negative the "orbit" is a pencil balanced on its tip — an exact solution no real system can hold. For the workhorse case of power-law forces the criterion collapses to a single number. Substituting F=−k rnF=-k\,r^{n} gives rF′/F=nrF'/F=n, and with it both the stability condition and the frequency of the radial breathing:

F(r)=−k rn:ωr2ωφ2=3+n  >0  ⟺  n>−3(4)F(r)=-k\,r^{n}:\qquad \frac{\omega_r^2}{\omega_\varphi^2}=3+n\;>0\;\Longleftrightarrow\; n>-3 \qquad (4)

where: ωφ=L/μr02\omega_\varphi=L/\mu r_0^2 — the orbital angular frequency, with ωφ2=−F(r0)/μr0\omega_\varphi^2=-F(r_0)/\mu r_0 from (2); ωr\omega_r — the radial oscillation frequency. Gravity (n=−2n=-2) gives ratio 11: the radius breathes exactly once per revolution, which is why a perturbed circular orbit is just a small Kepler ellipse. Hooke (n=+1n=+1) gives ratio 44, i.e. ωr=2ωφ\omega_r=2\omega_\varphi. Convention alert: many texts write F=−k/rpF=-k/r^{p} with p=−np=-n, turning the criterion into p<3p<3 — same physics, flipped sign. The same machinery runs beyond Newton. General relativity adds an attractive 1/r41/r^4-type force term to the Schwarzschild orbit equation, steepening the effective potential at small radius — and below a critical angular momentum the valley vanishes entirely:

Figure
The Schwarzschild effective potential — lower the angular momentum and the valley disappears at the ISCO
high L: valley (stable) plus inner hilltop (unstable) critical L: valley and hilltop merge — ISCO, r = 6GM/c² low L: no extremum at all — nothing to orbit on, plunge radius Schwarzschild effective potential, three values of L

Derivation B works this out: the inner hilltop and outer valley merge when L2=12 G2M2μ2/c2L^2=12\,G^2M^2\mu^2/c^2, at radius rISCO=6GM/c2r_{\mathrm{ISCO}}=6GM/c^2 — the innermost stable circular orbit. Inside it, no amount of angular momentum stabilizes a circle; matter spirals in. That single radius sets the inner edge of black-hole accretion disks.

Derivation AStability from linearizing the radial equation

The goal is the master criterion (3): a circular orbit is stable exactly when Veff′′(r0)>0V_{\mathrm{eff}}''(r_0)>0, with the perturbation oscillating at ωr2=Veff′′(r0)/μ\omega_r^2=V_{\mathrm{eff}}''(r_0)/\mu, and the equivalent force form F′/F+3/r0>0F'/F+3/r_0>0. The strategy is the universal recipe for stability of any equilibrium: expand the equation of motion to first order in the displacement and inspect the sign of the restoring coefficient. Watch the last step — converting to the force form divides by F(r0)F(r_0), which is negative, flipping the inequality.

Step 1. Perturb the radius. Write r(t)=r0+η(t)r(t)=r_0+\eta(t) with ∣η∣≪r0|\eta|\ll r_0, at fixed LL. The radial equation (1) becomes μη¨=−Veff′(r0+η)\mu\ddot\eta=-V_{\mathrm{eff}}'(r_0+\eta). Why: LL is conserved, so a radial nudge cannot change it — the perturbed motion lives in the same effective potential.

Step 2. Expand to first order. Veff′(r0+η)=Veff′(r0)+Veff′′(r0) η+O(η2)V_{\mathrm{eff}}'(r_0+\eta)=V_{\mathrm{eff}}'(r_0)+V_{\mathrm{eff}}''(r_0)\,\eta+\mathcal O(\eta^2), and the first term vanishes by the circular-orbit condition (2), leaving μη¨=−Veff′′(r0) η\mu\ddot\eta=-V_{\mathrm{eff}}''(r_0)\,\eta. Why: equilibrium kills the zeroth order; the curvature is the leading survivor and acts as a spring constant.

Step 3. Read off the verdict. If Veff′′(r0)>0V_{\mathrm{eff}}''(r_0)>0, η\eta oscillates harmonically at ωr=Veff′′(r0)/μ\omega_r=\sqrt{V_{\mathrm{eff}}''(r_0)/\mu} — bounded, stable. If negative, η∝e+∣ωr∣t\eta\propto e^{+|\omega_r|t} grows exponentially — unstable. Zero is marginal and needs higher orders. Why: this is the harmonic-oscillator trichotomy; the borderline case is decided by the first nonvanishing higher derivative.

Step 4. Convert to the force form. Differentiating Veff′=−F−L2/μr3V_{\mathrm{eff}}'=-F-L^2/\mu r^3 gives Veff′′(r0)=−F′(r0)+3L2μr04V_{\mathrm{eff}}''(r_0)=-F'(r_0)+\dfrac{3L^2}{\mu r_0^4}, and eliminating L2L^2 with (2): Veff′′(r0)=−F′(r0)−3F(r0)r0V_{\mathrm{eff}}''(r_0)=-F'(r_0)-\dfrac{3F(r_0)}{r_0}. Requiring this positive and dividing by −F(r0)>0-F(r_0)>0 yields F′(r0)F(r0)+3r0>0\dfrac{F'(r_0)}{F(r_0)}+\dfrac{3}{r_0}>0 — criterion (3). Why: the circular condition trades the orbit-specific LL for the force law itself, making the criterion a property of FF alone; the division is by a positive number, so the inequality survives intact.

Step 5. Specialize to power laws. For F=−k rnF=-k\,r^{n}, rF′/F=nrF'/F=n, so stability is 3+n>03+n>0, and with ωφ2=−F/μr0\omega_\varphi^2=-F/\mu r_0 the ratio is ωr2/ωφ2=3+n\omega_r^2/\omega_\varphi^2=3+n — equation (4). Why: one exponent now encodes stability (n>−3n>-3), the breathing rate, and — via the apsidal angle π/3+n\pi/\sqrt{3+n} — the precession story of the neighbouring articles. *Verified: matches Goldstein §3.3 and Taylor Ch. 8; n=−2n=-2 gives ratio 1 (perturbed orbit = small Kepler ellipse), n=+1n=+1 gives ωr=2ωφ\omega_r=2\omega_\varphi, both standard.*

Go deeper · Derivation B · 2 minThe innermost stable circular orbit of Schwarzschild geometry
Derivation BThe innermost stable circular orbit of Schwarzschild geometry

The goal is the most famous stability boundary in physics: rISCO=6GM/c2r_{\mathrm{ISCO}}=6GM/c^2, below which no stable circular orbit exists around a non-rotating black hole. The strategy reuses this article's toolkit unchanged — find where Veff′=0V_{\mathrm{eff}}'=0 has solutions, then ask when the valley survives — applied to the Schwarzschild effective potential, which differs from Newton's by one extra attractive 1/r31/r^3 term (we keep Newtonian-style notation). Watch the structure: the circular-orbit condition becomes a quadratic in rr, so orbits come in pairs — stable outer, unstable inner — and the ISCO is where the pair merges.

Step 1. Write the effective potential. For Schwarzschild geodesics the radial problem has Veff(r)=−GMμr+L22μr2−GML2μc2r3V_{\mathrm{eff}}(r)=-\dfrac{GM\mu}{r}+\dfrac{L^2}{2\mu r^2}-\dfrac{GML^2}{\mu c^2 r^3}. Why: the first two terms are exactly Newton's; the third is the relativistic correction — an attraction steepening like 1/r31/r^3, precisely the kind of term our criterion flags as dangerous.

Step 2. Find the circular orbits. Veff′(r)=0V_{\mathrm{eff}}'(r)=0 gives GMμr2−L2μr3+3GML2μc2r4=0\dfrac{GM\mu}{r^2}-\dfrac{L^2}{\mu r^3}+\dfrac{3GML^2}{\mu c^2 r^4}=0; multiplying by μr4/GMμ2\mu r^4/GM\mu^2 turns it into the quadratic r2−L2GMμ2 r+3L2μ2c2=0r^2-\dfrac{L^2}{GM\mu^2}\,r+\dfrac{3L^2}{\mu^2c^2}=0. Why: Newton's version lacks the last term and has a single root at every LL; the correction makes the equation quadratic, so circular orbits now exist only in pairs — or not at all.

Step 3. Solve and interpret the pair. r±=L22GMμ2[1±1−12 G2M2μ2L2c2]r_\pm=\dfrac{L^2}{2GM\mu^2}\left[1\pm\sqrt{1-\dfrac{12\,G^2M^2\mu^2}{L^2c^2}}\right]: r+r_+ is the valley floor (stable), r−r_- the inner hilltop (unstable). Why: the outer extremum is the minimum, the inner the maximum — exactly the high-LL curve of the figure.

Step 4. Merge the pair. The square root vanishes — valley and hilltop coalesce into an inflection, Veff′′=0V_{\mathrm{eff}}''=0 — when L2=12 G2M2μ2/c2L^2=12\,G^2M^2\mu^2/c^2. Substituting back: rISCO=L22GMμ2=6GMc2r_{\mathrm{ISCO}}=\dfrac{L^2}{2GM\mu^2}=\dfrac{6GM}{c^2}. Why: marginal stability is precisely the double-root condition; below this LL the discriminant is negative, no circular orbit of any kind exists, and matter plunges.

Step 5. Attach the numbers. 6GM/c26GM/c^2 is three Schwarzschild radii: 88.588.5 km for a 10 M⊙10\,M_\odot black hole, about 0.80.8 AU for the 4×106 M⊙4\times10^6\,M_\odot hole at the galactic centre. Why: the ISCO scales linearly with mass — stellar-mass holes pack their disk edges into tens of kilometres, making their orbital frequencies audible to gravitational-wave detectors. *Verified: rISCO=6GM/c2r_{\mathrm{ISCO}}=6GM/c^2 and LISCO=12 GMμ/cL_{\mathrm{ISCO}}=\sqrt{12}\,GM\mu/c are the standard Schwarzschild results (e.g. Hartle Ch. 9); the Newtonian limit c→∞c\to\infty removes the quadratic's second root, restoring stable circles at all radii.*

Limiting cases

CaseResultMeaning
Kepler, n=−2n=-2ωr=ωφ\omega_r=\omega_\varphi, ratio 3+n=13+n=1Stable; perturbed circle = small closed ellipse; no precession.
Hooke, n=+1n=+1ωr=2ωφ\omega_r=2\omega_\varphi, ratio 44Stable; centred ellipse, closes after half a revolution.
Inverse-cube force, n=−3n=-3ωr2=0\omega_r^2=0Marginal: no restoring force at linear order; any nudge destroys the orbit.
Steeper, n<−3n<-3ωr2<0\omega_r^2<0Extremum is a maximum: perturbations grow exponentially; spiral in or escape.
Schwarzschild, r<6GM/c2r<6GM/c^2no minimum of VeffV_{\mathrm{eff}}Relativity's 1/r31/r^3 potential term kills the valley: the ISCO, inner edge of accretion disks.

Worked examples

1 · Nudging the ISS (confidence builder)

The ISS orbits at r0=6771r_0=6771 km (400400 km altitude), where GM⊕=3.986×1014 m3/s2GM_\oplus=3.986\times10^{14}\ \mathrm{m^3/s^2}. Verify its circular frequency, check stability, and describe what a small radial kick does.

Step 1. Orbital frequency from the circular condition (2): ωφ=GM⊕/r03=3.986×1014/(6.771×106)3=1.13×10−3 rad/s\omega_\varphi=\sqrt{GM_\oplus/r_0^3}=\sqrt{3.986\times10^{14}/(6.771\times10^6)^3}=\boxed{1.13\times10^{-3}\ \mathrm{rad/s}}, i.e. a period 2π/ωφ=92.4 min2\pi/\omega_\varphi=\boxed{92.4\ \mathrm{min}} and speed ωφr0=7.67\omega_\varphi r_0=7.67 km/s.

Step 2. Stability: gravity has n=−2n=-2, and 3+n=1>03+n=1>0 — comfortably stable.

Step 3. Radial breathing frequency from (4): ωr=ωφ3+n=ωφ\omega_r=\omega_\varphi\sqrt{3+n}=\omega_\varphi. A small radial kick makes the radius oscillate with the same 92.4-minute period as the orbit itself.

Step 4. One radial cycle per revolution means the kicked orbit is a slightly eccentric Kepler ellipse that closes on itself — the station just trades a circle for a fixed ellipse. Stable, non-precessing.

2 · The last stable circle around a black hole (exam level)

A 10 M⊙10\,M_\odot black hole accretes gas that spirals inward through ever-smaller stable circles. Find the radius of the last one, its orbital frequency, and the gravitational-wave frequency a blob there would radiate. Use GM⊙=1.327×1020 m3/s2GM_\odot=1.327\times10^{20}\ \mathrm{m^3/s^2}, c=2.998×108c=2.998\times10^8 m/s.

Step 1. Gravitational length scale: GM/c2=10×1.327×1020/(2.998×108)2=1.475×104 m=14.75GM/c^2=10\times1.327\times10^{20}/(2.998\times10^8)^2=1.475\times10^4\ \mathrm{m}=14.75 km.

Step 2. ISCO radius from Derivation B: rISCO=6GM/c2=6×14.75=88.5 kmr_{\mathrm{ISCO}}=6GM/c^2=6\times14.75=\boxed{88.5\ \mathrm{km}}. Inside this the effective potential has no valley: gas cannot circle, only plunge.

Step 3. Orbital frequency (the Kepler form Ω=GM/r3\Omega=\sqrt{GM/r^3} holds exactly in Schwarzschild coordinates): Ω=1.327×1021/(8.85×104)3=1384 rad/s\Omega=\sqrt{1.327\times10^{21}/(8.85\times10^4)^3}=\boxed{1384\ \mathrm{rad/s}} — the last stable circle is traversed 220220 times per second.

Step 4. Orbital frequency f=Ω/2π=220f=\Omega/2\pi=220 Hz; a mass orbiting there radiates gravitational waves at twice the orbital frequency, fGW=440 Hzf_{\mathrm{GW}}=\boxed{440\ \mathrm{Hz}} — squarely in the LIGO band, which is why stellar-mass ISCOs are audible.

Origins · 1 minJoseph-Louis Lagrange and how this was discovered
JL
Joseph-Louis Lagrange
1736–1813 · Turin, Berlin and Paris

Lagrange rebuilt mechanics as analysis: the Mécanique analytique (1788) famously contains not a single diagram. Stability questions run through his greatest work — his 1772 essay on the three-body problem found the five equilibrium points that bear his name, and his analysis of when motion near equilibrium stays bounded seeded the method this article uses: expand about the equilibrium, keep the linear term, examine the restoring coefficient. The stable pair, L4 and L5, passed his test where the collinear three failed — and two centuries later the Trojan asteroids and the James Webb Space Telescope (parked near the unstable Sun–Earth L2, needing station-keeping exactly as the analysis predicts) live out the distinction daily.

History
1687
— Newton's Principia derives circular-orbit conditions and probes apsidal motion via the revolving-orbits theorem.
1772
— Lagrange finds the five three-body equilibrium points and analyzes their stability; L4 and L5 pass.
1892
— Poincaré's Méthodes nouvelles turns orbital stability into a general mathematical discipline.
1916
— Schwarzschild solves Einstein's equations; the solution's effective potential hides an innermost stable orbit.
1973
— Novikov and Thorne build the standard accretion-disk model with its inner edge pinned at the ISCO.

Going further

References. Morin, Introduction to Classical Mechanics, ch. 7 (Central Forces) · Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., §3.3, §3.6 (Bertrand's theorem) · Landau & Lifshitz, Mechanics, 3rd ed., §14, §15 · Taylor, Classical Mechanics, ch. 8 · Hartle, Gravity, Ch. 9 (Schwarzschild orbits and the ISCO)

Provenance. Verifier should recompute: (i) Veff′′(r)=−F′(r)+3L2/(μr4)V_{\mathrm{eff}}''(r)=-F'(r)+3L^2/(\mu r^4) and the resulting force-form criterion F′/F+3/r0>0F'/F+3/r_0>0; (ii) the power-law identity ωr2/ωθ2=3−n\omega_r^2/\omega_\theta^2=3-n giving the n<3n<3 stability bound; (iii) ωθ2=−F(r0)/(μr0)\omega_\theta^2=-F(r_0)/(\mu r_0) and the apsidal angle Φ=π/3−n\Phi=\pi/\sqrt{3-n} with Φ=π\Phi=\pi (Kepler) and Φ=π/2\Phi=\pi/2 (Hooke). Cross-check against Goldstein §3.3 & §3.6, Landau-Lifshitz Mechanics §14–15, Taylor ch. 8, and Morin ch. 7. (The provenance identities use the inverse-power convention F=−k/rnF=-k/r^{n}; the article text uses F=−k rnF=-k\,r^{n}, so n→−nn\to-n: ratio 3+n3+n, bound n>−3n>-3, Φ=π/3+n\Phi=\pi/\sqrt{3+n}.) Added for this upgrade and recomputed: ISS ω=1.13×10−3\omega=1.13\times10^{-3} rad/s, T=92.4T=92.4 min, v=7.67v=7.67 km/s from GM⊕=3.986×1014GM_\oplus=3.986\times10^{14}, r=6.771×106r=6.771\times10^6 m; ISCO r=6GM/c2=88.5r=6GM/c^2=88.5 km, Ω=1384\Omega=1384 rad/s, f=220f=220 Hz, fGW=440f_{\mathrm{GW}}=440 Hz for 10 M⊙10\,M_\odot; LISCO2=12G2M2μ2/c2L_{\mathrm{ISCO}}^2=12G^2M^2\mu^2/c^2 (Hartle Ch. 9).

Builds on: The Effective Potential · Reduction to the Radial Equation · Angular Momentum Conservation in Central Fields · Central Force Fundamentals Leads to: The Kepler Problem · Classification of Orbits by Energy · The Orbit Equation and Binet's Equation · Kepler's Third Law: Period-Semimajor Relation

Where does the 3 in the stability criterion actually come from?

From the centrifugal force's own steepness. The barrier term in the radial equation is L2/μr3L^2/\mu r^3 — an outward force falling as 1/r31/r^3, logarithmic slope −3-3. Stability is a tug-of-war of slopes: displaced inward, the body needs the outward barrier to grow faster than the attraction does, which is exactly rF′/F>−3rF'/F>-3. So the 3 is not numerology; it is the exponent of the centrifugal force, which traces to angular-momentum conservation via L=μr2φ˙L=\mu r^2\dot\varphi. Any physics that modified the centrifugal balance would move the 3 — and general relativity effectively does, which is how the ISCO appears.

My textbook says stability needs n less than 3 — this article says greater than −3. Who is wrong?

Nobody; the two statements are the same inequality in different conventions. Many texts write the force as F=−k/rpF=-k/r^{p} with the exponent in the denominator, so Kepler is p=2p=2 and stability reads p<3p<3. This article (matching its sibling orbit articles) writes F=−k rnF=-k\,r^{n}, so Kepler is n=−2n=-2 and the same criterion reads n>−3n>-3; evidently n=−pn=-p. The invariant statement, safe in any convention, is logarithmic: dln⁡∣F∣/dln⁡r>−3\mathrm{d}\ln|F|/\mathrm{d}\ln r>-3 — the attraction must fall off more slowly than the inverse-cube centrifugal force. When in doubt, compute rF′/FrF'/F and add 3.

Is an unstable circular orbit a real solution or just a mathematical artifact?

It is a perfectly exact solution — start with precisely the right radius and velocity and the body circles forever. What it is not is realizable: like a pencil on its tip, any perturbation grows exponentially at rate ∣ωr∣|\omega_r|, and perturbations are never absent. Unstable circular orbits still matter: they organize the dynamics around them. In Schwarzschild geometry the unstable branch r−r_- governs how close a flyby can approach before plunging, and the related unstable photon circular orbit at 3GM/c23GM/c^2 shapes the black-hole shadow imaged by the Event Horizon Telescope.

Why does a kicked satellite trace a closed ellipse instead of drifting or precessing?

Because for inverse-square gravity the two frequencies coincide: ωr2/ωφ2=3+n=1\omega_r^2/\omega_\varphi^2=3+n=1 at n=−2n=-2. The radius completes exactly one in-out breath per revolution, so the point of closest approach stays at a fixed direction — the perturbed orbit is a genuine small-eccentricity Kepler ellipse, not a drifting spiral. Any other exponent breaks the coincidence: the perturbed orbit is still bounded, but its apsides rotate each radial cycle, tracing a rosette. That ratio, and which force laws make it rational everywhere, is the subject of the apsidal-angles and Bertrand's-theorem articles.

Newtonian gravity gives stable circles at every radius — why does relativity suddenly forbid them inside the ISCO?

Because general relativity adds an extra attractive term to the effective potential that steepens like 1/r31/r^3 — beyond the inverse-cube danger threshold this article establishes. Far from the hole the term is negligible and the Newtonian valley survives; close in, it overwhelms the centrifugal barrier. The valley and the inner hilltop approach each other as LL decreases and merge at L2=12G2M2μ2/c2L^2=12G^2M^2\mu^2/c^2, radius 6GM/c26GM/c^2 — Derivation B. Inside that radius the landscape is all downhill: circular motion is not merely unstable but nonexistent, and gas at the disk's inner edge plunges within a few dynamical times.

Does the stability of a circular orbit depend on how much angular momentum the body has?

For power-law forces, no — and that is worth internalizing. The criterion 3+n>03+n>0 contains only the exponent: changing LL moves the location of the circular orbit but not the sign of the curvature there, so every circular orbit of a given power law is stable or every one is unstable. This changes the moment the force is not a pure power: then rF′/FrF'/F varies with radius, and different LL values park the orbit in stable or unstable zones of the same potential. Schwarzschild is the canonical example — stability holds for r0>6GM/c2r_0>6GM/c^2 and fails inside, which is precisely why the ISCO is defined by an angular-momentum threshold, L2≥12G2M2μ2/c2L^2\geq12G^2M^2\mu^2/c^2.

Before this

Where this leads