Innermost Stable Circular Orbit (ISCO)

The innermost stable circular orbit is the sharpest place where general relativity parts company with Newton. In Newtonian gravity a body can circle a mass at any radius, and every such circle is stable — nudge the orbit and it merely wobbles about its old track. Schwarzschild geometry forbids this. Below a definite radius, rISCO=6GM/c2r_{\mathrm{ISCO}}=6GM/c^2, no stable circular orbit exists at all: the smallest inward push sends matter spiralling through the horizon. That radius has no free parameters — for a non-rotating hole it is exactly three Schwarzschild radii, independent of the particle's energy or angular momentum. The result flows from one idea built in the preceding articles: the conserved energy and angular momentum of the timelike and azimuthal Killing vectors collapse the geodesic equation into one-dimensional motion in an effective potential Veff(r)V_{\mathrm{eff}}(r). Circular orbits sit at its extrema, stable ones in its minima. As the angular momentum is lowered, the minimum (a stable orbit) and the neighbouring maximum (an unstable one) slide together and annihilate at a single inflection point. That coalescence — Veff′=Veff′′=0V_{\mathrm{eff}}'=V_{\mathrm{eff}}''=0 — is the ISCO. Its consequences are enormous. A particle lowered quasi-statically from rest at infinity to the ISCO must shed a fixed fraction of its rest energy, 1−22/3≈5.7%1-2\sqrt2/3\approx5.7\%, and this single number caps the efficiency of every thin accretion disc around a Schwarzschild black hole — an order of magnitude beyond nuclear fusion. The ISCO marks the inner edge of that disc, sets the inspiral-to-plunge transition in a merger, and, through its strong dependence on spin, is a principal handle by which astronomers weigh how fast real black holes rotate.

Key takeaways
  • For a massive test particle in Schwarzschild geometry, stable circular orbits exist only outside rISCO=6GM/c2=3 rsr_{\mathrm{ISCO}}=6GM/c^2=3\,r_s — the innermost stable circular orbit.
  • The ISCO is where the effective potential's stable minimum and unstable maximum merge into an inflection: Veff′=Veff′′=0V_{\mathrm{eff}}'=V_{\mathrm{eff}}''=0 at the same radius.
  • There the specific angular momentum and specific energy are both minimised, L~=12 GM/c\tilde L=\sqrt{12}\,GM/c and E~=223c2≈0.943 c2\tilde E=\tfrac{2\sqrt2}{3}c^2\approx0.943\,c^2.
  • The binding energy at the ISCO, 1−223≈5.7%1-\tfrac{2\sqrt2}{3}\approx5.7\% of rest mass, sets the maximum radiative efficiency of a thin accretion disc — versus 0.7%0.7\% for hydrogen fusion.
  • Inside the ISCO matter plunges on a nearly geodesic path with little further radiation, so the ISCO is the effective inner edge of the disc and a fiducial inspiral-to-merger point for gravitational waves.
  • For a maximally spinning Kerr hole the prograde ISCO shrinks to GM/c2GM/c^2 (efficiency ≈42%\approx42\%) and the retrograde one grows to 9GM/c29GM/c^2 (≈3.8%\approx3.8\%), which is how the ISCO measures black-hole spin.
  • Key idea: the ISCO is the radius at which the stable minimum and unstable maximum of the effective potential coalesce — where Veff′=Veff′′=0V_{\mathrm{eff}}'=V_{\mathrm{eff}}''=0 — giving the spin-independent Schwarzschild value r=6GM/c2r=6GM/c^2 and a maximum disc efficiency of ≈5.7%\approx5.7\%.
The core insight

The ISCO is the radius at which the stable minimum and unstable maximum of the effective potential coalesce — where Veff′=Veff′′=0V_{\mathrm{eff}}'=V_{\mathrm{eff}}''=0 — yielding the spin-independent Schwarzschild value r=6GM/c2r=6GM/c^2 and a maximum disc efficiency of ≈5.7%\approx5.7\%.

Where this lives today

The ISCO is now a measured quantity, not just a textbook radius. In X-ray binaries and active galactic nuclei two techniques read it directly: continuum-fitting infers the inner-disc radius from the thermal spectrum, and X-ray reflection spectroscopy fits the relativistically broadened, redshifted iron Kα\alpha line emitted just outside the ISCO. Because the Kerr ISCO ranges from GM/c2GM/c^2 to 9GM/c29GM/c^2, locating it pins the dimensionless spin a/Ma/M — this is how holes such as Cygnus X-1 and MCG–6-30-15 are found to spin, some near extremal. In gravitational-wave astronomy the ISCO frequency fGW≈4.4 kHz×(M⊙/M)f_{\mathrm{GW}}\approx4.4\,\mathrm{kHz}\times(M_\odot/M) marks the inspiral-to-plunge transition and anchors the matching of post-Newtonian and numerical-relativity waveforms used to read every LIGO–Virgo–KAGRA detection. And because the ISCO exists only through the relativistic 1/r31/r^3 term, any departure of its location from the Kerr prediction would flag modified gravity or a horizonless compact object — a strong-field test complementary to the photon ring imaged around M87 and Sgr A by the Event Horizon Telescope.

The physical picture

The Newtonian effective potential is −GM/r+L2/2r2-GM/r+L^2/2r^2: a −1/r-1/r well and a +1/r2+1/r^2 centrifugal barrier. However small LL is, the barrier always wins at small radius, so the potential turns up steeply and every circular orbit sits in a well — stable forever. General relativity adds a third term, −GML2/(c2r3)-GML^2/(c^2r^3), that falls off faster than the barrier and eventually overwhelms it. Lower the angular momentum and the well grows shallow; at a critical value it flattens completely, and below that no well is left to hold an orbit. The first figure shows this family of potentials and the moment the well disappears.

Figure
The effective potential Veff(r) for three values of specific angular momentum L: with L = 4 GM/c a stable minimum sits at r = 12GM/c² behind an unstable barrier at r = 4GM/c², the critical value L = √12 GM/c ≈ 3.46 GM/c merges them into one flat inflection at the ISCO r = 6GM/c², and L = 3 GM/c is below the minimum so no circular orbit exists and the particle plunges
Lowering L flattens the well: at L = √12 GM/c the ISCO appears at r = 6GM/c² 1.00 0.95 0.90 3 6 12 18 radius r in units of GM/c² effective potential Veff /c² ISCO r = 6GM/c² unstable barrier (max) L = 4 GM/c: stable well + barrier L = √12 GM/c: marginal inflection L = 3 GM/c: no orbit, plunge

Think of the ISCO as the edge of a cliff for orbiting matter. Just outside it a small inward nudge gives a bounded oscillation about the circular radius — the orbit is stable. Just inside it the same nudge finds no restoring force, because the radial "spring constant" Veff′′V_{\mathrm{eff}}'' has gone negative, and the particle runs away inward. The second figure follows what astronomers actually measure: the specific energy of a circular orbit versus radius. It bottoms out precisely at the ISCO, and the depth of that dip below the rest-mass line is the energy an accreting particle can radiate.

Figure
Specific energy E of a circular orbit versus radius r: the stable branch outside r = 6GM/c² and the unstable branch between 3 and 6 GM/c² meet at the ISCO minimum r = 6GM/c², where E = 0.943c² is least, so the maximum binding energy a thin disc can radiate is 1 − 0.943 = 5.7% of the rest mass, and E diverges at the photon sphere r = 3GM/c²
The punchline: circular-orbit energy E(r) bottoms at the ISCO — binding energy 5.7% 1.10 1.05 1.00 0.95 3 6 10 15 20 radius r in units of GM/c² specific energy E /c² photon sphere 3GM/c² : E → ∞ ISCO 6GM/c² E = c² : unbound / rest at infinity marginally bound r = 4GM/c² E = 0.943c² (minimum) binding 5.7% stable branch (r above 6GM/c²) unstable branch

The core equations

The starting point is the radial equation of motion for a timelike geodesic, obtained by feeding the two conserved quantities into the four-velocity normalisation. It reads as an energy equation for one-dimensional motion in rr; every result below comes from it.

(drdτ)2=E~2−Veff2(r)(1)\left(\frac{\mathrm{d}r}{\mathrm{d}\tau}\right)^{2}=\tilde E^{2}-V_{\mathrm{eff}}^{2}(r) \qquad (1)

where: rr — Schwarzschild (areal) radial coordinate; τ\tau — proper time; E~\tilde E — conserved energy per unit rest mass (units of c2c^2); VeffV_{\mathrm{eff}} — effective potential for radial motion. Reading it: a turning point occurs where r˙=0\dot r=0, i.e. E~=Veff\tilde E=V_{\mathrm{eff}}; a circular orbit needs that turning point to persist, forcing Veff′=0V_{\mathrm{eff}}'=0. The effective potential is the product of the metric's redshift factor and a centrifugal-plus-rest term. Written as a square it is the cleanest object to differentiate, and its extrema and inflections coincide with those of VeffV_{\mathrm{eff}} since Veff>0V_{\mathrm{eff}}>0.

Veff2(r)=(1−2GMc2r)(c2+L~2r2)(2)V_{\mathrm{eff}}^{2}(r)=\left(1-\frac{2GM}{c^{2}r}\right)\left(c^{2}+\frac{\tilde L^{2}}{r^{2}}\right) \qquad (2)

where: L~\tilde L — conserved angular momentum per unit rest mass; MM — central mass; G,cG,c — Newton's constant and light speed. Reading it: expanding the product yields a term −GML~2/(c2r3)-GM\tilde L^2/(c^2r^3) with no Newtonian counterpart — that single relativistic term destroys the potential well and creates the ISCO. Any circular orbit lives at an extremum; whether it is stable depends on the sign of the second derivative — these two conditions separate stable orbits (minima) from unstable ones (maxima). The ISCO is the borderline case where well and barrier have just annihilated, so both derivatives vanish together, one condition more than an ordinary circular orbit and hence an isolated radius.

Veff′(r)=0,  Veff′′(r)>0  (stable);Veff′(r)=Veff′′(r)=0  (ISCO)(3)V_{\mathrm{eff}}'(r)=0,\ \ V_{\mathrm{eff}}''(r)>0\ \ \text{(stable)};\qquad V_{\mathrm{eff}}'(r)=V_{\mathrm{eff}}''(r)=0\ \ \text{(ISCO)} \qquad (3)

where: Veff′V_{\mathrm{eff}}' — first radial derivative, vanishing at every circular orbit; Veff′′V_{\mathrm{eff}}'' — second derivative, positive in a well (restoring force, bounded oscillation), negative on a hilltop (runaway), zero at the marginal ISCO inflection. Reading it: the extra Veff′′=0V_{\mathrm{eff}}''=0 is exactly what promotes "some circular orbit" into "the innermost stable one". Solving Veff′=0V_{\mathrm{eff}}'=0 for the angular momentum shows how much spin a circular orbit needs at each radius; the formula diverges at the photon sphere and bottoms out at the ISCO.

L~2(r)=GM r1−3GM/(c2r)(4)\tilde L^{2}(r)=\frac{GM\,r}{1-3GM/(c^{2}r)} \qquad (4)

where: L~(r)\tilde L(r) — specific angular momentum required for a circular orbit at radius rr. Reading it: as r→3GM/c2r\to3GM/c^2 (the photon sphere) the required L~→∞\tilde L\to\infty, and the function turns over at r=6GM/c2r=6GM/c^2, so no circular orbit needs less than L~min⁡=12 GM/c\tilde L_{\min}=\sqrt{12}\,GM/c. Imposing the double condition (3) collapses everything to a single number: the ISCO radius is a fixed multiple of the geometric mass, independent of the particle.

rISCO=6GMc2=3 rs(5)r_{\mathrm{ISCO}}=\frac{6GM}{c^{2}}=3\,r_s \qquad (5)

where: rISCOr_{\mathrm{ISCO}} — coordinate radius of the innermost stable circular orbit; rs=2GM/c2r_s=2GM/c^2 — the Schwarzschild radius. Reading it: the ISCO sits at three Schwarzschild radii, twice the photon-sphere radius, independent of L~\tilde L and E~\tilde E. Evaluating the conserved quantities on the ISCO gives their minimal values; the energy fixes the accretion efficiency.

L~ISCO=12 GMc,E~ISCO=223 c2≈0.9428 c2(6)\tilde L_{\mathrm{ISCO}}=\sqrt{12}\,\frac{GM}{c},\qquad \tilde E_{\mathrm{ISCO}}=\frac{2\sqrt2}{3}\,c^{2}\approx0.9428\,c^{2} \qquad (6)

where: L~ISCO\tilde L_{\mathrm{ISCO}} — minimum specific angular momentum of any stable circular orbit; E~ISCO\tilde E_{\mathrm{ISCO}} — specific energy there. Reading it: since E~ISCO<c2\tilde E_{\mathrm{ISCO}}<c^2, a particle from infinity is bound, and the deficit is what it radiated on the way in. That deficit, expressed as a fraction of rest energy, is the maximum radiative efficiency of a thin disc.

η=1−E~ISCOc2=1−223≈0.0572(7)\eta=1-\frac{\tilde E_{\mathrm{ISCO}}}{c^{2}}=1-\frac{2\sqrt2}{3}\approx0.0572 \qquad (7)

where: η\eta — fraction of rest-mass energy radiated spiralling from rest at infinity to the ISCO. Reading it: 5.7%5.7\% is a ceiling, not a rate — inside the ISCO the plunging gas carries its remaining 94.3%94.3\% across the horizon largely unradiated. Finally, the orbital angular velocity keeps the Newtonian Keplerian form even in the strong field, and doubling it gives the dominant gravitational-wave frequency at the ISCO.

Ω=GMr3,fGW=Ωπ∣rISCO=c363/2π GM(8)\Omega=\sqrt{\frac{GM}{r^{3}}},\qquad f_{\mathrm{GW}}=\frac{\Omega}{\pi}\Big|_{r_{\mathrm{ISCO}}}=\frac{c^{3}}{6^{3/2}\pi\,GM} \qquad (8)

where: Ω\Omega — coordinate angular velocity; fGWf_{\mathrm{GW}} — dominant (quadrupole) wave frequency at the ISCO. Reading it: it scales as 1/M1/M, giving the fiducial fGW≈4.4 kHz×(M⊙/M)f_{\mathrm{GW}}\approx4.4\,\mathrm{kHz}\times(M_\odot/M) that marks the end of the slow inspiral.

Derivation ADeriving the effective potential from the geodesic equation

The goal is to show how a geodesic in four-dimensional curved spacetime becomes one-dimensional motion in a potential, since that reduction is the entire engine behind the ISCO. Two symmetries of the metric — time-translation and axial rotation — supply two constants of motion through their Killing vectors, and the four-velocity normalisation supplies a third relation; combining them collapses the geodesic equations into a single first-order equation for rr whose potential carries a purely relativistic term absent from Newton.

Step 1. Restrict to the equatorial plane θ=π/2\theta=\pi/2 and write the four-velocity normalisation gμνuμuν=−c2g_{\mu\nu}u^\mu u^\nu=-c^2 explicitly: −(1−2GM/c2r)c2t˙2+(1−2GM/c2r)−1r˙2+r2ϕ˙2=−c2-(1-2GM/c^2r)c^2\dot t^2+(1-2GM/c^2r)^{-1}\dot r^2+r^2\dot\phi^2=-c^2 (overdot is d/dτ\mathrm{d}/\mathrm{d}\tau). Why: this one constraint contains all the radial dynamics once tt and ϕ\phi are eliminated.

Step 2. Introduce the conserved quantities of the timelike and azimuthal Killing vectors, E~=(1−2GM/c2r)c2t˙\tilde E=(1-2GM/c^2r)c^2\dot t and L~=r2ϕ˙\tilde L=r^2\dot\phi, each constant because contracting a Killing vector with uμu^\mu gives a quantity whose proper-time derivative vanishes. Why: these fix t˙\dot t and ϕ˙\dot\phi, leaving rr as the only dynamical variable.

Step 3. Substitute t˙,ϕ˙\dot t,\dot\phi back and rearrange: r˙2=E~2/c2−(1−2GM/c2r)(c2+L~2/r2)\dot r^2=\tilde E^2/c^2-(1-2GM/c^2r)(c^2+\tilde L^2/r^2). Identifying the last product as Veff2(r)V_{\mathrm{eff}}^2(r) gives the energy equation (1). Why: the geodesic is now a particle of energy E~\tilde E rolling in the potential VeffV_{\mathrm{eff}} — the Newtonian orbit tool made relativistic.

Step 4. Expand the potential: Veff2=c2−2GMc2/r+L~2/r2−2GML~2/(c2r3)V_{\mathrm{eff}}^2=c^2-2GMc^2/r+\tilde L^2/r^2-2GM\tilde L^2/(c^2r^3). The first three pieces mirror Newton (rest energy, gravity, centrifugal barrier); the fourth, negative, is new. Why: this −1/r3-1/r^3 term outgrows the +1/r2+1/r^2 barrier at small radius, and is the necessary and sufficient reason an ISCO exists at all. *Verified: expanding (1−2/r)(1+L~2/r2)(1-2/r)(1+\tilde L^2/r^2) (with GM/c2=1GM/c^2=1) reproduces the four terms; dropping the 1/r31/r^3 piece recovers a well stable for every rr, per Carroll §5.4 and Hartle Ch. 9.*

Go deeper · Derivation B · 1 minLocating the ISCO where the potential's extrema merge
Derivation BLocating the ISCO where the potential's extrema merge

The aim is to turn the geometric statement "the minimum and maximum of the potential coalesce" into an algebraic solution for the radius, and to show the answer is independent of the orbiting particle. Working with W(r)≡Veff2(r)W(r)\equiv V_{\mathrm{eff}}^2(r) is legitimate because Veff>0V_{\mathrm{eff}}>0, so WW shares every extremum and inflection with VeffV_{\mathrm{eff}}; imposing W′=W′′=0W'=W''=0 gives two equations in rr and L~\tilde L, and eliminating L~\tilde L leaves a single clean root. Throughout, u≡GM/c2u\equiv GM/c^2.

Step 1. Differentiate W=(1−2u/r)(c2+L~2/r2)W=(1-2u/r)(c^2+\tilde L^2/r^2) once: W′(r)=2uc2/r2−2L~2/r3+6uL~2/r4W'(r)=2uc^2/r^2-2\tilde L^2/r^3+6u\tilde L^2/r^4. Setting W′=0W'=0 and clearing r4/2r^4/2 gives the circular-orbit relation uc2r2−L~2r+3uL~2=0uc^2r^2-\tilde L^2 r+3u\tilde L^2=0. Why: this is the general circular-orbit condition; solving it for L~2\tilde L^2 yields L~2=uc2r2/(r−3u)\tilde L^2=uc^2r^2/(r-3u), the angular momentum of equation (4).

Step 2. Differentiate again: W′′(r)=−4uc2/r3+6L~2/r4−24uL~2/r5W''(r)=-4uc^2/r^3+6\tilde L^2/r^4-24u\tilde L^2/r^5. Setting W′′=0W''=0 and clearing r5/2r^5/2 gives −2uc2r2+3L~2r−12uL~2=0-2uc^2r^2+3\tilde L^2 r-12u\tilde L^2=0. Why: this second relation is the extra condition that marks the boundary of stability — it holds only at the inflection where a minimum and maximum have just met.

Step 3. Eliminate L~2\tilde L^2 by substituting L~2=uc2r2/(r−3u)\tilde L^2=uc^2r^2/(r-3u) from Step 1 into the Step 2 relation, then divide through by the non-zero factor uc2r2/(r−3u)uc^2r^2/(r-3u). What remains is −2(r−3u)+3r−12u=0-2(r-3u)+3r-12u=0, which simplifies to r−6u=0r-6u=0. Why: the two conserved-quantity–dependent equations combine into one equation for rr alone, so the ISCO radius cannot depend on L~\tilde L.

Step 4. Read off rISCO=6u=6GM/c2r_{\mathrm{ISCO}}=6u=6GM/c^2 and back-substitute: L~ISCO2=uc2(6u)2/(6u−3u)=12u2c2\tilde L_{\mathrm{ISCO}}^2=uc^2(6u)^2/(6u-3u)=12u^2c^2, so L~ISCO=12 GM/c\tilde L_{\mathrm{ISCO}}=\sqrt{12}\,GM/c. Why: this is the minimum of L~(r)\tilde L(r) — no stable circular orbit needs less angular momentum, i.e. none exists below the ISCO. *Verified: substituting r=6GM/c2r=6GM/c^2, L~2=12(GM/c2)2c2\tilde L^2=12(GM/c^2)^2c^2 into W′=0W'=0 and W′′=0W''=0 satisfies both identically; matches MTW §25.5 and Shapiro & Teukolsky §12.4.*

Go deeper · Derivation C · 1 minThe ISCO energy and the maximum accretion efficiency
Derivation CThe ISCO energy and the maximum accretion efficiency

This derivation extracts the single number that makes the ISCO astrophysically famous: the fraction of rest energy accreting matter can radiate. Because a circular orbit has r˙=0\dot r=0, equation (1) gives E~2=Veff2(r)=W(r)\tilde E^2=V_{\mathrm{eff}}^2(r)=W(r), so evaluating WW at the ISCO returns the specific energy directly. Subtracting it from the rest-mass value c2c^2 gives the binding energy, and since a thin disc radiates exactly this deficit as gas drifts inward, that binding energy is the disc's maximum efficiency.

Step 1. Evaluate the redshift factor at the ISCO: 1−2GM/(c2rISCO)=1−2u/6u=2/31-2GM/(c^2r_{\mathrm{ISCO}})=1-2u/6u=2/3. Why: this factor is the gravitational time-dilation the orbiting matter sits in, and it carries most of the energy loss.

Step 2. Evaluate the centrifugal bracket with L~ISCO2=12u2c2\tilde L_{\mathrm{ISCO}}^2=12u^2c^2 and rISCO2=36u2r_{\mathrm{ISCO}}^2=36u^2: c2+L~ISCO2/rISCO2=c2(1+1/3)=43c2c^2+\tilde L_{\mathrm{ISCO}}^2/r_{\mathrm{ISCO}}^2=c^2(1+1/3)=\tfrac43c^2. Why: keeping exact fractions makes the answer a clean closed form.

Step 3. Combine the two factors: E~ISCO2=23⋅43c2⋅c2=89c4\tilde E_{\mathrm{ISCO}}^2=\tfrac23\cdot\tfrac43c^2\cdot c^2=\tfrac89c^4, so E~ISCO=8/9 c2=223c2≈0.9428 c2\tilde E_{\mathrm{ISCO}}=\sqrt{8/9}\,c^2=\tfrac{2\sqrt2}{3}c^2\approx0.9428\,c^2. Why: this is the energy of a particle sitting on the ISCO, the lowest specific energy of any stable circular orbit.

Step 4. Form the binding energy per unit rest mass, η=(c2−E~ISCO)/c2=1−223≈0.0572\eta=(c^2-\tilde E_{\mathrm{ISCO}})/c^2=1-\tfrac{2\sqrt2}{3}\approx0.0572. A particle lowered quasi-statically from rest at infinity (E~=c2\tilde E=c^2) to the ISCO must therefore radiate about 5.7%5.7\% of its rest energy. Why: for a Novikov–Thorne / Shakura–Sunyaev thin disc this deficit is emitted as heat and light, so η\eta is the accretion efficiency ceiling for a Schwarzschild hole. *Verified: 22/3=0.942809…2\sqrt2/3=0.942809\ldots, η=0.057191…\eta=0.057191\ldots, the standard ≈5.7%\approx5.7\% of Shapiro & Teukolsky §12.4 and MTW §25.5; the second figure's energy curve bottoms at 0.94280.9428 exactly at r=6GM/c2r=6GM/c^2.*

Limiting cases

The ISCO reduces to memorable closed forms in the cases that matter most; the spin cases show that 6GM/c26GM/c^2 is only the non-rotating anchor of a wide range.

CaseResultMeaning
Schwarzschild (non-rotating)rISCO=6GM/c2=3 rsr_{\mathrm{ISCO}}=6GM/c^2=3\,r_sCanonical value, independent of the particle; twice the photon sphere.
Extremal Kerr, prograderISCO=GM/c2=12rsr_{\mathrm{ISCO}}=GM/c^2=\tfrac12 r_sMaximal co-rotation (a=GM/c2a=GM/c^2) drags the ISCO down to the horizon in Boyer–Lindquist rr; efficiency rises to 1−1/3≈42%1-1/\sqrt3\approx42\%.
Extremal Kerr, retrograderISCO=9GM/c2r_{\mathrm{ISCO}}=9GM/c^2Counter-rotating orbits are pushed outward; efficiency falls to 1−25/27≈3.8%1-\sqrt{25/27}\approx3.8\%.
Marginally bound orbitrmb=4GM/c2, E~=c2r_{\mathrm{mb}}=4GM/c^2,\ \tilde E=c^2The circular orbit with exactly escape energy; between it and the ISCO orbits are bound but unstable.
Newtonian limitVeff′′>0V_{\mathrm{eff}}''>0 for all rrWith no 1/r31/r^3 term every circular orbit is stable — there is no ISCO. The phenomenon is intrinsically relativistic.

Worked examples

1 · The ISCO of Sagittarius A* (fully worked)

The black hole at the Galactic Centre has M=4.3×106 M⊙M=4.3\times10^{6}\,M_\odot. Find the radius of its ISCO in kilometres, in astronomical units, and in Schwarzschild radii, and the orbital period of matter sitting on it. Use GM⊙/c2=1.477GM_\odot/c^2=1.477 km and GM⊙/c3=4.925×10−6GM_\odot/c^3=4.925\times10^{-6} s.

Step 1. Geometric mass. Scale the solar value: GM/c2=4.3×106×1.477 km=6.35×106GM/c^2=4.3\times10^{6}\times1.477\ \mathrm{km}=6.35\times10^{6} km — the natural length unit for every radius here.

Step 2. ISCO radius. Six geometric masses: rISCO=6GM/c2=3.81×107 km=0.25 AUr_{\mathrm{ISCO}}=6GM/c^2=\mathbf{3.81\times10^{7}\ km}=\mathbf{0.25\ AU} (a quarter of the Earth–Sun distance), and since rISCO=3 rsr_{\mathrm{ISCO}}=3\,r_s always, exactly three Schwarzschild radii — about 55 R⊙55\,R_\odot.

Step 3. Angular velocity. With Ω=c3/(63/2GM)\Omega=c^3/(6^{3/2}GM), GM/c3=21.2GM/c^3=21.2 s and 63/2=14.706^{3/2}=14.70: Ω=1/(14.70×21.2 s)=3.21×10−3 rad/s\Omega=1/(14.70\times21.2\ \mathrm{s})=3.21\times10^{-3}\ \mathrm{rad/s}.

Step 4. Orbital period. T=2π/Ω=1955 s≈33 minutesT=2\pi/\Omega=\mathbf{1955\ s\approx33\ minutes}. Matter on the ISCO of Sgr A* orbits in barely over half an hour — close to the tens-of-minutes timescale of the infrared flares the GRAVITY interferometer has tracked near the hole, a direct look at strong-field motion.

2 · How much a quasar must swallow (fully worked)

A luminous quasar radiates L=1039 WL=10^{39}\ \mathrm{W}. For a thin disc at the Schwarzschild efficiency ceiling, what accretion rate does this need, and how does it compare with fusion?

Step 1. Efficiency relation. A disc radiates at L=ηM˙c2L=\eta\dot M c^2, so M˙=L/(ηc2)\dot M=L/(\eta c^2) with η=1−22/3=0.0572\eta=1-2\sqrt2/3=0.0572.

Step 2. Insert numbers. ηc2=0.0572×9.0×1016 J/kg=5.1×1015 J/kg\eta c^2=0.0572\times9.0\times10^{16}\ \mathrm{J/kg}=5.1\times10^{15}\ \mathrm{J/kg}, so M˙=1039/5.1×1015=1.9×1023 kg/s\dot M=10^{39}/5.1\times10^{15}=\mathbf{1.9\times10^{23}\ kg/s}.

Step 3. In solar masses per year. Multiply by 3.156×107 s/yr3.156\times10^{7}\ \mathrm{s/yr} and divide by M⊙=1.99×1030M_\odot=1.99\times10^{30} kg: M˙=1.9×1023×3.156×107/1.99×1030=3.1 M⊙/yr\dot M=1.9\times10^{23}\times3.156\times10^{7}/1.99\times10^{30}=\mathbf{3.1\ M_\odot/yr}. A bright quasar devours a few Suns of gas each year.

Step 4. Contrast with fusion. Fusion liberates only ≈0.007\approx0.007 of rest mass, so it would need ≈25 M⊙/yr\approx\mathbf{25\ M_\odot/yr} — eight times the fuel. Accretion down to the ISCO is the most efficient steady power source known, which is why quasars, not fusion, dominate the luminous universe. The same physics gives a fiducial gravitational-wave frequency fGW≈4.4 kHz×(M⊙/M)f_{\mathrm{GW}}\approx4.4\,\mathrm{kHz}\times(M_\odot/M) — about 7373 Hz for a 60 M⊙60\,M_\odot binary, in LIGO's band, marking the inspiral-to-plunge handoff.

Physics of the portrait

The ISCO belongs to the era when relativists turned Schwarzschild's and Kerr's geometries into working astrophysics — computing exactly how orbits behave near a black hole and what that means for the light we see.

Origins · 1 minJames Bardeen and how this was discovered
JB
James Bardeen
1939–2022 · the relativistic theory of black-hole orbits and accretion efficiency

James Maxwell Bardeen — son of John Bardeen, the only person to win two Nobel Prizes in physics — became one of the deepest relativists of his generation. In the 1972 paper "Rotating Black Holes: Locally Nonrotating Frames, Energy Extraction, and Scalar Synchrotron Radiation," written with William Press and Saul Teukolsky, he laid out the geodesic structure of the Kerr metric in full: the innermost stable circular orbit as a function of spin, its binding energy, and hence the efficiency of accretion — from 5.7%5.7\% for a non-rotating hole up to 42%42\% for a maximally rotating one. That calculation turned the ISCO from a curiosity of Schwarzschild geometry into a quantitative diagnostic of black-hole spin. Bardeen's reach went further: with Brandon Carter and Stephen Hawking he formulated the four laws of black-hole mechanics that seeded black-hole thermodynamics, and in cosmology he introduced the gauge-invariant "Bardeen potentials" behind the modern theory of cosmological perturbations. The ISCO efficiencies that now let astronomers weigh a black hole's spin from its X-ray glow are, in large part, his legacy.

History
1916
— Schwarzschild publishes the exterior metric whose geodesics contain the last stable orbit
1963
— Kerr finds the rotating black-hole solution, making the ISCO depend on spin and orbit sense
1969
— Lynden-Bell proposes that accretion discs around massive black holes power quasars
1972
— Bardeen, Press & Teukolsky compute the ISCO and its binding energy across all Kerr spins
1973
— Novikov & Thorne and Shakura & Sunyaev build the thin-disc model whose efficiency is fixed by the ISCO

Going further

Texts. Misner, Thorne & Wheeler, Gravitation (1973), §25.5 · Wald, General Relativity (1984), §6.3 · Carroll, Spacetime and Geometry (2004), §5.4 · Hartle, Gravity (2003), Ch. 9 · Schutz, A First Course in General Relativity, 2nd ed. (2009), §11.1 · Shapiro & Teukolsky, Black Holes, White Dwarfs and Neutron Stars (1983), §12.4 · Bardeen, Press & Teukolsky, ApJ 178, 347 (1972).

Provenance. The effective potential, ISCO radius, and efficiency are standard Schwarzschild results. A verifier should: (1) reproduce Veff2=(1−2GM/c2r)(c2+L~2/r2)V_{\mathrm{eff}}^2=(1-2GM/c^2r)(c^2+\tilde L^2/r^2) from the timelike-geodesic normalisation and the two Killing conserved quantities (Carroll §5.4, Hartle Ch. 9); (2) confirm rISCO=6GM/c2r_{\mathrm{ISCO}}=6GM/c^2 and L~ISCO=12 GM/c\tilde L_{\mathrm{ISCO}}=\sqrt{12}\,GM/c by solving W′=W′′=0W'=W''=0 (MTW §25.5, Wald §6.3); (3) verify E~ISCO=223c2=0.9428c2\tilde E_{\mathrm{ISCO}}=\tfrac{2\sqrt2}{3}c^2=0.9428c^2 and η=1−223=0.0572\eta=1-\tfrac{2\sqrt2}{3}=0.0572 (Shapiro & Teukolsky §12.4); (4) check the numerical results rISCO(Sgr A∗)=3.8×107r_{\mathrm{ISCO}}(\mathrm{Sgr\,A^*})=3.8\times10^{7} km =0.25=0.25 AU with a 3333-min period, quasar accretion 3.1 M⊙/yr3.1\,M_\odot/\mathrm{yr} at 103910^{39} W, fGW≈4.4 kHz×(M⊙/M)f_{\mathrm{GW}}\approx4.4\,\mathrm{kHz}\times(M_\odot/M); and (5) the Kerr limits rISCO=GM/c2r_{\mathrm{ISCO}}=GM/c^2 (prograde, η=1−1/3≈42%\eta=1-1/\sqrt3\approx42\%) and 9GM/c29GM/c^2 (retrograde, η≈3.8%\eta\approx3.8\%) against Bardeen, Press & Teukolsky (1972). No non-standard values are asserted.

Prerequisites. The Schwarzschild Metric · Effective Potential for Massive-Particle Orbits · Circular Orbits and Their Stability · Killing Vectors and Conserved Quantities in Schwarzschild · Geodesics in Schwarzschild Spacetime. Leads to. Perihelion Precession in Schwarzschild Geometry · Photon Orbits and the Photon Sphere · Gravitational Redshift in Schwarzschild Spacetime.

Ask the tutor: Ask me to walk through any step — from why the 1/r31/r^3 term creates the ISCO, to how the 5.7%5.7\% efficiency lights up a quasar, or how black-hole spin slides the ISCO between GM/c2GM/c^2 and 9GM/c29GM/c^2.

Why does the ISCO exist in general relativity but not in Newtonian gravity?

It comes down to one extra term. Newton's potential −GM/r+L2/2r2-GM/r+L^2/2r^2 always turns up steeply at small rr because the centrifugal barrier +L2/2r2+L^2/2r^2 diverges faster than gravity, so every circular orbit sits in a well and is stable. General relativity adds an attractive −GML2/(c2r3)-GML^2/(c^2r^3) piece that dominates at small radius and pulls the potential back down; as the angular momentum is reduced the minimum and neighbouring maximum annihilate at r=6GM/c2r=6GM/c^2, and below that no stable circular orbit exists — a feature with no Newtonian counterpart.

Is the ISCO the same as the event horizon or the photon sphere?

No — they are three distinct radii, all fixed multiples of GM/c2GM/c^2. The event horizon sits at 2GM/c22GM/c^2, the unstable photon sphere (last circular orbit for light) at 3GM/c23GM/c^2, and the ISCO (last stable orbit for massive particles) at 6GM/c26GM/c^2 — twice the photon sphere, three times the horizon. Between the photon sphere and the ISCO circular orbits for matter still exist but are all unstable; only outside 6GM/c26GM/c^2 can matter settle into a lasting circular track, which is why accretion discs end there and not at the horizon.

Where does the 5.7% accretion efficiency actually come from?

It is a binding energy. A particle at rest far from the hole has E~=c2\tilde E=c^2; on the ISCO it has E~ISCO=223c2≈0.943c2\tilde E_{\mathrm{ISCO}}=\tfrac{2\sqrt2}{3}c^2\approx0.943c^2. In a thin disc the gas drifts inward through nearly circular orbits, radiating the energy difference between neighbouring radii until it reaches the ISCO; the total shed is the deficit 1−0.943=0.0571-0.943=0.057, or 5.7%5.7\% of rest energy. Inside the ISCO there is no stable orbit, so the gas plunges across the horizon carrying its remaining 94.3%94.3\% largely unradiated. The 5.7%5.7\% is thus a ceiling on how much accreted mass a Schwarzschild disc can turn into light.

How does black-hole spin change the ISCO?

Dramatically, and in a direction that depends on which way you orbit. For a Kerr hole the ISCO depends on the spin aa and on whether the orbit is prograde (co-rotating) or retrograde. As the spin rises from zero to the extremal a=GM/c2a=GM/c^2, the prograde ISCO shrinks from 6GM/c26GM/c^2 down to GM/c2GM/c^2, while the retrograde ISCO grows to 9GM/c29GM/c^2. Because binding energy tracks the ISCO, the efficiency swings from 5.7%5.7\% (Schwarzschild) up to ≈42%\approx42\% (prograde) or down to ≈3.8%\approx3.8\% (retrograde). This strong, monotonic dependence is exactly what lets astronomers read a black hole's spin from the inner edge of its disc.

What happens to matter once it crosses inside the ISCO?

It plunges. Just inside 6GM/c26GM/c^2 the potential has no minimum, so with no restoring force a gas element follows a nearly radial, nearly geodesic path to the horizon in a short proper time. This plunge is fast compared with the slow viscous inspiral outside and radiates very little further energy — which is why the ISCO is the effective inner edge of the disc and why the disc spectrum cuts off at a temperature set by the ISCO radius. In the standard thin-disc model the plunging region is treated as a torque-free boundary, an idealisation that works because so little happens there.

Is r = 6GM/c² a real distance you could measure with a ruler?

Not directly — it is a coordinate radius, not a proper distance. The Schwarzschild rr is defined so a sphere at fixed rr has area 4πr24\pi r^2 (the areal radius), but radial ruler distance is stretched by (1−2GM/c2r)−1/2>1(1-2GM/c^2r)^{-1/2}>1, so the proper distance from the horizon to the ISCO exceeds the coordinate gap Δr=4GM/c2\Delta r=4GM/c^2. What is observable is not ruler distance but quantities like the orbital frequency Ω=GM/r3\Omega=\sqrt{GM/r^3} and the redshift of light from the orbit — which is how the ISCO is actually located in data. "6GM/c26GM/c^2" is a precise statement about the areal coordinate, not tape-measure distance.

Before this

Where this leads