# Black hole > [[PORTAL_Physics|Physics]] spine. <!-- MICROSIMGEN:BEGIN v1.7 — hand-placed to match siblings; regenerate with g08_place_microsims.py (§15) --> ## Microsims — three.js ### Black hole (three.js) <div class="microsim-player"> <!-- MICROSIM:PENDING_DEPLOY:BEGIN v1.7 g08 — embed target is not on the CDN; restore with g08 --undeploy-clear --> <p class="wt-pending"><strong>Microsim staged, not yet on the CDN.</strong> <code>Black_hole.html</code> is built and deploy-ready in <code>_3d_deploy_stage/</code>, but the Netlify project has not been redeployed yet. The player is disabled until the deploy lands; the explanatory text below is unchanged.</p> <!-- <iframe src="https://wikitube-3d-microsims.netlify.app/Black_hole.html" width="100%" height="620" frameborder="0" loading="lazy" sandbox="allow-scripts allow-same-origin" title="Black hole — three.js microsim"></iframe> --> <!-- MICROSIM:PENDING_DEPLOY:END --> </div> **Open it full-screen:** [Black_hole.html](https://wikitube-3d-microsims.netlify.app/Black_hole.html) · library `threejs` · route `microsim/threejs/` ### Related microsims Live sims on neighbouring articles: - [[Gravitational_lens]] - [[Spacetime]] - [[Stellar_structure]] - [[General_relativity]] - [[Gravity]] - [[Photon]] *Sim hosted off-article; the article owns the reference, not the runtime (WIKI_RULES §10.4).* <!-- MICROSIMGEN:END --> ## Overview Almost every popular image of a black hole gets the geometry backwards. The hole is drawn as an aperture — a puncture that things fall through — and its defining property is presented as an appetite. Neither is right. A non-rotating black hole is a piece of empty, static, spherically symmetric geometry specified entirely by one number, and the interesting things about it are not appetites but **radii**. Three of them matter: | radius | value | what happens there | |---|---|---| | event horizon | *r*<sub>s</sub> = 2*GM*/*c*² | the one-way causal boundary | | photon sphere | *r*<sub>ph</sub> = 3*GM*/*c*² | the only circular orbit light can have | | ISCO | *r*<sub>ISCO</sub> = 6*GM*/*c*² | innermost *stable* circular orbit for matter | Everything the microsim shows follows from firing test objects past those radii. The central fact is that the capture threshold is **sharp and calculable**. A photon arriving from far away is fully labelled by its impact parameter *b* — the perpendicular distance between the hole and the straight line the photon would have followed had gravity been absent — and there is a critical value > *b*<sub>crit</sub> = 3√3 *GM*/*c*² = 5.196 *GM*/*c*² = 2.598 *r*<sub>s</sub> Photons with *b* < *b*<sub>crit</sub> are swallowed; photons with *b* > *b*<sub>crit</sub> get out. No fuzzy transition, no probability: the boundary is a closed-form number. Because the beam is three-dimensional the threshold is not a line but a **cylinder** of radius *b*<sub>crit</sub> aimed at the hole, and the sim draws it explicitly. What makes this worth watching rather than merely reading is the behaviour *just outside* the threshold. A ray at *b* = 1.0001 *b*<sub>crit</sub> does not simply bend a lot — it spirals in to just above the photon sphere, winds around it two and a half times, and spirals back out, and the winding count grows without bound as *b* approaches the threshold. That looping is the direct cause of the bright ring in an Event Horizon Telescope image: light that would have missed the telescope is wrapped around the hole and delivered to it, piling up in an annulus just outside the dark patch. The dark patch is the shadow of the capture cylinder, angular radius *b*<sub>crit</sub>/*D*. Switch the test object to a massive particle and the ISCO appears for a reason you can watch: the effective potential grows a barrier whose two extrema — an unstable circular orbit and a stable one — merge and vanish at exactly 6*GM*/*c*². Below that there are no stable circular orbits at all, which is why an accretion disc has an inner edge. Two further things this article is careful about, because they are the errors that survive almost every popular account: 1. **A black hole is not a cosmic vacuum cleaner.** Outside the horizon, at large radius, its field is *exactly* the Newtonian field of the same mass. Replace the Sun with a solar-mass black hole and the Earth's orbit is unchanged to a part in 10⁸. 2. **The horizon is not a surface anybody notices.** An infalling observer crosses it in finite proper time and feels nothing local at the crossing. It is the *distant* observer for whom the crossing never happens, and that is a statement about causal structure and about a coordinate chart, not about a wall. Both are attacked directly in the sim: the first by a preset and a Newtonian overlay, the second by a falling clock drawn twice — once on its own proper time and once as a distant telescope would actually see it. ## The physics ### The metric, and what "radius" means Karl Schwarzschild found the exterior solution for a static, spherically symmetric mass within weeks of Einstein's field equations being published (Schwarzschild 1916). In the coordinates that bear his name, with *G* = *c* = 1 and *r*<sub>s</sub> = 2*M*, > d*s*² = −(1 − *r*<sub>s</sub>/*r*) d*t*² + (1 − *r*<sub>s</sub>/*r*)⁻¹ d*r*² + *r*²(d*θ*² + sin²*θ* d*φ*²) Three points govern everything below. First, *r* is **not** a distance from the centre: it is defined so that a sphere of coordinate *r* has proper area 4π*r*². Second, *t* is the proper time of a static observer infinitely far away, not anybody's local clock. Third, by **Birkhoff's theorem** (Birkhoff 1923) this is the *unique* spherically symmetric vacuum solution, so the exterior field of a solar-mass black hole and the exterior field of the Sun beyond its surface are the same geometry with the same *M*. That theorem is misconception 1, stated as mathematics. ### The orbit equation Spherical symmetry means every geodesic stays in a plane through the centre, so the problem is two-dimensional in (*r*, *φ*) and the sim places that plane in 3-D afterwards. Staticity and axisymmetry give two constants of motion per unit mass: the energy *E* and the angular momentum *h*. Substituting *u* = 1/*r* and eliminating the affine parameter turns the geodesic equation into one of the tidiest expressions in general relativity (Misner, Thorne & Wheeler 1973, ch. 25): > **photon** d²*u*/d*φ*² + *u* = 3*Mu*² > **matter** d²*u*/d*φ*² + *u* = *M*/*h*² + 3*Mu*² Delete the 3*Mu*² term and the matter equation becomes *u*″ + *u* = *M*/*h*², whose solutions are exactly the Newtonian conic sections — closed, non-precessing ellipses, parabolas and hyperbolas. So **the entire difference between Newton and Einstein for this problem is one cubic term**, and the sim's "Newtonian prediction" checkbox is literally that term switched off in the same integrator. The first integrals fix the initial slope, > (d*u*/d*φ*)² = 1/*b*² − *u*²(1 − 2*Mu*) (photon) > (d*u*/d*φ*)² = (*E*² − 1)/*h*² + 2*Mu*/*h*² − *u*² + 2*Mu*³ (matter) and the coordinate clock comes along for free as d*t*/d*φ* = *E r*²/[*h*(1 − 2*M*/*r*)]. The sim integrates the second-order form with RK4 at d*φ* = 0.003 rad from *r* = 400 *GM*/*c*², where the trajectory is straight to a quarter of a per cent, so *b* means the asymptotic *L*/*E* and not some local aiming offset. ### Where 3√3 comes from For a photon the radial equation is (d*r*/d*λ*)² = *E*² − *L*²*f*(*r*)/*r*² with *f* = 1 − 2*M*/*r*. The effective potential *L*²*f*/*r*² has a single maximum, > d/d*r* [(1 − 2*M*/*r*)/*r*²] = 0 ⟹ *r* = 3*M* which is the **photon sphere**: the one radius at which light can travel in a circle. It is a maximum, not a minimum, so the orbit is unstable — a knife edge, which is exactly why near-critical rays linger there and then leave. Setting the peak equal to *E*², with *b* = *L*/*E*: > *b*<sub>crit</sub>² = *r*³/(*r* − 2*M*) at *r* = 3*M* ⟹ *b*<sub>crit</sub>² = 27*M*² ⟹ *b*<sub>crit</sub> = 3√3 *M* The capture cross-section is therefore σ = π*b*<sub>crit</sub>² = 27π*G*²*M*²/*c*⁴, an exact result with no free parameters, and the sim prints it live in square metres for whatever mass you dial in. ### The strong-deflection limit and the photon ring The exact deflection is a quadrature over the cubic above, which the sim evaluates for the highlighted ray after removing the square-root endpoint singularity (substituting *u* = *u*<sub>p</sub>(1 − *s*²)). It has two clean limits. In the weak field it reduces to *α* → 4*GM*/(*c*²*b*), the Einstein value — 1.751 arcsec at the solar limb, confirmed at the 1919 eclipse (Dyson, Eddington & Davidson 1920). Approaching the threshold it diverges logarithmically (Darwin 1959; recast systematically by Bozza 2002): > *α*(*b*) ≈ −ln(*b*/*b*<sub>crit</sub> − 1) + ln[216(7 − 4√3)] − π The sim's integrator reproduces this to about 10⁻³ rad out to *b*/*b*<sub>crit</sub> − 1 = 10⁻⁶, which is a useful check that the geodesics on screen are the real ones. Each factor of *e* closer to the threshold buys roughly another 1/(2π) of a turn, so the "Photon ring" preset — which sets *b*/*b*<sub>crit</sub> − 1 = 1.5 × 10⁻⁴ — gets about 1.75 loops and a total deflection near 487°. Repeated windings produce an infinite nested sequence of ever-fainter photon subrings just outside the shadow edge; only the *n* = 1 ring contributes meaningfully to current EHT images. ### Matter: the effective potential and the ISCO For a massive test particle the radial equation is > (d*r*/d*τ*)² = *E*² − *V*(*r*), *V*(*r*) = (1 − 2*M*/*r*)(1 + *h*²/*r*²) Circular orbits sit at extrema of *V*, which gives > *r*<sub>±</sub> = (*h*²/2*M*)[1 ± √(1 − 12*M*²/*h*²)] with *r*<sub>−</sub> the unstable orbit (the barrier peak) and *r*<sub>+</sub> the stable one; the sim draws both as rings, orange and green, updated live. The whole ISCO story is in that square root: the two orbits exist only for *h* > 2√3 *M*, and at *h* = 2√3 *M* they merge at *r* = 6*M* and vanish. Below 6*GM*/*c*² there are no stable circular orbits at all — a particle nudged inward does not settle into a smaller orbit, it plunges. That is why a thin accretion disc has an inner edge. Two numbers fall out of the same algebra. The binding energy at the ISCO is 1 − √(8/9) = 5.72 %, the classic radiative efficiency of a thin Schwarzschild disc (hydrogen fusion releases 0.7 %). And the critical impact parameter for matter arriving with asymptotic speed *v*, which the sim gets by bisection, interpolates between > *v* → *c*: *b*<sub>crit</sub> → 3√3 *GM*/*c*² (the photon value) > *v* ≪ *c*: *b*<sub>crit</sub> → 4*GM*/(*c*²·*v*/*c*) The slow limit is worth a pause. A *Newtonian* body of radius *r*<sub>s</sub> would capture slow particles only inside *b* = 2*GM*/(*c*²·*v*/*c*), so general relativity swallows four times the area. For light the same comparison gives exactly 27/8 = 3.375, since *b*<sub>crit</sub>² = 27*M*² against the Michell–Laplace "dark star" value of 8*M*². The sim draws both circles at the entry plane, so the ratio is something you can see rather than a claim you have to accept. ### Misconception 1: not a vacuum cleaner The strongest form of the argument is Birkhoff's theorem: outside a spherical mass the geometry depends on *M* alone, so swapping the Sun for a one-solar-mass black hole changes *nothing* for the planets. The "Swap the Sun" preset makes it concrete and reports: - *r*<sub>s</sub> = 2.95 km, photon sphere 4.43 km, ISCO 8.86 km, *b*<sub>crit</sub> = 7.67 km. - Earth orbits at 1 AU = 5.07 × 10⁷ *r*<sub>s</sub>, where the post-Newtonian parameter *GM*/*rc*² is 9.9 × 10⁻⁹. - The entire relativistic effect on Earth's orbit is a perihelion drift of 6π*GM*/[*c*²*a*(1 − *e*²)] = 3.84 arcsec **per century** — the same formula that gives Mercury's famous 43 arcsec per century. - The hole is a *smaller* target than the Sun by a factor of 8.2 × 10⁹ in area: its light-capture disc has radius 7.67 km against the Sun's 696 000 km. That last number is the cleanest refutation available: replacing the Sun with a black hole makes the solar system's central object roughly ten billion times *harder* to hit. Black holes accrete efficiently only when matter arrives with almost no angular momentum, or when viscosity in a disc strips angular momentum away over many orbits. The geometry alone does not do it. Switch on the Newtonian overlay and drag the impact parameter inward to watch Newtonian intuition fail. The two curves are indistinguishable at *b* ≳ 10 *r*<sub>s</sub>, separate noticeably by 5 *r*<sub>s</sub>, and once the periapsis nears the photon sphere the Newtonian prediction is not merely inaccurate but qualitatively wrong — it has the ray escaping when in fact it is captured. For light that Newtonian curve is the 1784 Michell corpuscle, a historical curiosity rather than a theory (Newtonian gravity has nothing to say about massless particles); it gives half the correct deflection, 2*GM*/(*c*²*b*), which is also the value Einstein published in 1911. ### Misconception 2: the horizon is not a place The sim drops a clock from rest at the ISCO and draws it **twice**. The first marker follows the clock's own worldline. Radial free fall from rest at *R* has the exact cycloid solution *r* = (*R*/2)(1 + cos *η*), *τ* = √(*R*³/8*M*)(*η* + sin *η*), so the clock reaches *r* = 0 in **finite proper time**. Nothing marks the crossing: no discontinuity in *τ*, no jolt, no local measurement that could detect it. The only thing that is locally measurable is the tidal field, ~2*GM*/*r*<sub>s</sub>³ = *c*⁶/4*G*²*M*², which *decreases* with mass — at a 10 M<sub>☉</sub> hole the stretch across a 1.4 m body at the horizon is about 10⁷ *g* and lethal, at M87* about 3 × 10⁻¹¹ *g* and undetectable. The longest possible proper time between crossing the horizon and reaching *r* = 0 is π*GM*/*c*³: 15.5 μs for the Sun, 28 hours for M87*. The second marker shows where a distant static observer *sees* the clock: the coordinate time along the worldline plus the exact radial light-travel delay built from the tortoise coordinate, > *t*<sub>arr</sub> = *t*<sub>emit</sub> + (*R*<sub>obs</sub> − *r*) + *r*<sub>s</sub> ln[(*R*<sub>obs</sub> − *r*<sub>s</sub>)/(*r* − *r*<sub>s</sub>)] Both terms diverge logarithmically as *r* → *r*<sub>s</sub>, so on the distant clock the image asymptotes to the horizon and never crosses. It also fades. The received-to-emitted frequency ratio for a radially infalling emitter is (1 − *r*<sub>s</sub>/*r*)/(*E* + *V*), which equals the familiar static factor √(1 − *r*<sub>s</sub>/*r*) at the moment of release and goes to zero at the horizon. Since the observed flux of a point source scales as the *fourth* power of that ratio, and near the horizon *r* − *r*<sub>s</sub> ∝ e<sup>−*t*/*r*<sub>s</sub></sup>, the image does not merely redden — it goes dark exponentially, e-folding every *r*<sub>s</sub>/*c*: 9.8 μs for a solar-mass hole, 17.8 hours for M87*. "Frozen at the horizon forever" is technically true and observationally meaningless. ### What this looks like for real black holes | object | *M* | *r*<sub>s</sub> | shadow diameter 2*b*<sub>c</sub>/*D* | measured | |---|---|---|---|---| | Sun (hypothetical) | 1 M<sub>☉</sub> | 2.95 km | — | — | | Cygnus X-1 | 21.2 M<sub>☉</sub> | 62.6 km | — | — | | Sgr A* | 4.297 × 10⁶ M<sub>☉</sub> | 1.27 × 10⁷ km | 53.3 μas | 51.8 ± 2.3 μas | | M87* | 6.5 × 10⁹ M<sub>☉</sub> | 128 AU | 39.7 μas | 42 ± 3 μas | The last two rows are the point: the predicted shadow diameter follows from the mass and distance with no adjustable parameters, and the EHT's measured ring diameters agree. The sim computes exactly this number, live, whenever a preset supplies a distance. ## Controls -> what each maps to | control | symbol | range / units | what it does | |---|---|---|---| | mass **M** | *M* | 1 to 10¹⁰ M<sub>☉</sub>, log scale | Sets the physical scale only. Schwarzschild geometry is scale-free in *M*, so the picture is identical for every mass; the HUD's kilometres, microarcseconds and microseconds are what change. | | preset buttons | *M*, *D* | Sun / 10 M<sub>☉</sub> / Sgr A* / M87* | Loads a real mass and, where known, a real distance so the shadow's angular diameter can be shown. | | test object | — | photon \| massive | Switches the source term in the orbit equation between 3*Mu*² and *M*/*h*² + 3*Mu*². | | speed at infinity **v** | *v*/*c* | 0.40 to 0.999 (massive only) | Sets *E* = *γ* and *h* = *b γv*, and therefore moves *b*<sub>crit</sub> between 5.196 *GM*/*c*² and 4*GM*/(*c*²*v*/*c*). | | impact parameter **b** | *b* = *L*/*E* | 0 to ≈2.8 *b*<sub>crit</sub>; scale stretched at the threshold | The asymptotic perpendicular offset of the incoming trajectory. Sets the fate of the highlighted ray. | | rays in fan | *N* | 1 to 25 | Number of geodesics fired. Spread is 1.0 *b*<sub>crit</sub> about the primary, with azimuths on the golden angle so the bundle fills three dimensions. | | photon sphere | *r*<sub>ph</sub> | on/off | The 3*GM*/*c*² shell plus the closed photon orbit in the primary ray's plane. | | ISCO + accretion disc | *r*<sub>ISCO</sub> | on/off | The 6*GM*/*c*² ring plus 1300 instanced particles on exact circular geodesics, Ω = √(*GM*/*r*³). | | capture cross-section | *b*<sub>crit</sub> | on/off | The capture cylinder of radius *b*<sub>crit</sub>, with the Michell–Laplace √2 *r*<sub>s</sub> circle inside it for comparison. | | infalling clock + image | *τ*, *t* | on/off | Free fall from rest at the ISCO, drawn twice: proper time, and what a static observer at 30 *GM*/*c*² actually sees. | | Newtonian prediction | — | on/off | Re-integrates the primary trajectory with the 3*Mu*² term deleted and draws it dashed. | | camera | — | external \| down the beam \| following the faller | "Down the beam" looks back along the incoming rays, which is the geometry an EHT image records. | | Photon ring | *b* | preset | *b*/*b*<sub>crit</sub> = 1.00015 — about 1.75 windings of the photon sphere before escape. | | Swap the Sun | *M*, *b* | preset | Misconception 1, with the live arithmetic in the caption. | | Pause / Reset | — | — | Freezes the animation or restores every default. | | Kerr spin *a*/*M* | *a* | **not implemented** | This is Schwarzschild, *a* = 0. No frame dragging, no ergosphere. See below. | ## Learning objective **Replace the intuition of "a hole that sucks" with the intuition of "a geometry with three radii and one threshold".** After using this sim a reader should be able to say, without looking anything up, that a photon's fate is decided by a single number *b* compared against 3√3 *GM*/*c*²; that the sharpness of that threshold is what makes a black-hole image have an edge; that the bright ring in that image is light that went round and came back; and that the accretion disc stops at 6*GM*/*c*² because that is where stable circular orbits run out. Specific things to try: - Park the impact-parameter slider at 1.00000 and step outward one notch at a time: the winding count falls as ln(*b*/*b*<sub>crit</sub> − 1), and the deflection with it. - Set 25 rays and sweep *b* through the threshold. The fan splits into a red half and a green half along a boundary that does not blur. - Switch to a massive particle and drag *b* down. The orange and green circular-orbit rings converge and vanish at 6*GM*/*c*² as *h* falls to 2√3 *M*. - Press "Swap the Sun", read the caption, then turn on the Newtonian overlay and drag *b* inward until the two curves visibly disagree — at around 10 *r*<sub>s</sub>. - Watch the two infall markers separate. They are the same object; only the clock differs. ## Limits and connections This is a **test-particle treatment in a fixed background**, and the honest list of what that excludes is long. **No backreaction, no self-force.** Every trajectory is an exact geodesic of an unchanging metric. Real infalling matter carries energy and momentum, perturbs the geometry and radiates gravitationally. The approximation is excellent for photons and for anything much lighter than the hole; it fails completely for the comparable-mass mergers LIGO observes (Abbott et al. 2016). **No Kerr, stated plainly.** The spin parameter *a*/*M* is not implemented. Astrophysical black holes are expected to rotate, often rapidly, and rotation changes the numbers that matter: for a maximally spinning hole the prograde ISCO moves from 6*GM*/*c*² to *GM*/*c*², the retrograde ISCO to 9*GM*/*c*², and the thin-disc radiative efficiency rises from 5.7 % to 42 % (Bardeen, Press & Teukolsky 1972). Rotation also drags inertial frames, producing an **ergosphere** outside the horizon — the region *r* < *M* + √(*M*² − *a*²cos²*θ*), where no observer can remain static — from which the Penrose process can extract rotational energy. None of that appears here, and nothing in this sim should be read as evidence about spin. One thing does survive: a Kerr shadow is not a circle, but its mean angular size differs from the Schwarzschild value by only a few per cent, which is why EHT mass measurements are relatively insensitive to spin. **No accretion physics.** The disc particles ride exact circular geodesics with Ω = √(*GM*/*r*³) — Kepler's third law, unchanged, a genuine result — but there is no viscosity, no magnetorotational instability, no inflow, no radiative transfer, no Doppler beaming and no lensing of the far side of the disc over the top of the hole. A physically rendered thin disc looks very different, as Luminet (1979) showed in the first such calculation. The disc's *colours* are a decorative stand-in for the Shakura–Sunyaev *T* ∝ *r*<sup>−3/4</sup> run (Shakura & Sunyaev 1973). Only two things about the disc here carry physics: its inner edge sits on the ISCO, and it shears differentially. **Coordinates, not images.** The sim plots Schwarzschild coordinate positions; it is not a ray-traced view, so background stars are not lensed and there is no Einstein ring from sources behind the hole. Relatedly, the "distant observer sees it freeze" behaviour is a fact about the Schwarzschild chart and about causal structure — the coordinate singularity at *r* = *r*<sub>s</sub> is an artefact removed by better coordinates (Eddington–Finkelstein, Kruskal–Szekeres; see Finkelstein 1958). Inside the horizon *r* becomes a timelike coordinate, so drawing the faller's interior track as a shrinking radius is a schematic, not a picture of a place. **No quantum effects.** Hawking radiation is absent (Hawking 1975). For a stellar-mass hole the Hawking temperature is about 6 × 10⁻⁹ K and the evaporation time exceeds 10⁶⁸ years, so ignoring it is overwhelmingly correct for any astrophysical object. Whether information survives evaporation remains genuinely unsettled and is not touched here. **Numerics.** Geodesics use RK4 at d*φ* = 0.003 rad from *r* = 400 *GM*/*c*², capped at 13 windings and 34 000 steps; rays inside about *b*/*b*<sub>crit</sub> − 1 = 10⁻⁶ are reported as "orbiting" rather than resolved. The HUD's deflection angle comes from an exact quadrature, not from the drawn arc, so it is the asymptotic value. **Connections.** The weak-field end of the same equation is [[Gravitational_lens]], where 4*GM*/*c*²*b* governs everything from Einstein rings to cosmological mass mapping. The metric itself belongs to [[Spacetime]] and [[General_relativity]]; the Newtonian limit that the overlay makes visible is [[Gravity]]. The test object with *m* = 0 is a [[Photon]], and the reason stellar-mass black holes exist at all — the failure of degeneracy pressure at the end of a massive star's life — is [[Stellar_structure]]. ## References - Schwarzschild, K. (1916). "Über das Gravitationsfeld eines Massenpunktes nach der Einsteinschen Theorie." *Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften zu Berlin*, 189–196. English translation: S. Antoci and A. Loinger, arXiv:physics/9905030. - Einstein, A. (1916). "Die Grundlage der allgemeinen Relativitätstheorie." *Annalen der Physik* **354**(7), 769–822. DOI: 10.1002/andp.19163540702. - Misner, C. W., Thorne, K. S., and Wheeler, J. A. (1973). *Gravitation*. W. H. Freeman. ISBN 978-0-7167-0344-0. (Princeton University Press reprint, 2017, ISBN 978-0-691-17779-3.) Chapters 25 and 33 carry the orbit equation, the effective potentials and the radial-infall solutions used here. - Chandrasekhar, S. (1983). *The Mathematical Theory of Black Holes*. Oxford University Press. ISBN 978-0-19-850370-5. - Wald, R. M. (1984). *General Relativity*. University of Chicago Press. ISBN 978-0-226-87033-5. - Hartle, J. B. (2003). *Gravity: An Introduction to Einstein's General Relativity*. Addison-Wesley. ISBN 978-0-8053-8662-2. - Taylor, E. F., and Wheeler, J. A. (2000). *Exploring Black Holes: Introduction to General Relativity*. Addison Wesley Longman. ISBN 978-0-201-38423-9. - Birkhoff, G. D. (1923). *Relativity and Modern Physics*. Harvard University Press. (The uniqueness theorem was obtained independently by J. T. Jebsen in 1921.) - Michell, J. (1784). "On the Means of discovering the Distance, Magnitude, &c. of the Fixed Stars…" *Philosophical Transactions of the Royal Society of London* **74**, 35–57. The original "dark star" argument; Laplace gave an equivalent one in *Exposition du système du monde* (1796). - Dyson, F. W., Eddington, A. S., and Davidson, C. (1920). "A Determination of the Deflection of Light by the Sun's Gravitational Field, from Observations Made at the Total Eclipse of May 29, 1919." *Philosophical Transactions of the Royal Society A* **220**, 291–333. - Darwin, C. (1959). "The gravity field of a particle." *Proceedings of the Royal Society A* **249**, 180–194. The original strong-deflection logarithm. - Bozza, V. (2002). "Gravitational lensing in the strong field limit." *Physical Review D* **66**, 103001. DOI: 10.1103/PhysRevD.66.103001. - Shapiro, I. I. (1964). "Fourth Test of General Relativity." *Physical Review Letters* **13**, 789. DOI: 10.1103/PhysRevLett.13.789. The signal-delay effect that makes rays crawl near the hole in this sim. - Kerr, R. P. (1963). "Gravitational Field of a Spinning Mass as an Example of Algebraically Special Metrics." *Physical Review Letters* **11**, 237. DOI: 10.1103/PhysRevLett.11.237. - Bardeen, J. M., Press, W. H., and Teukolsky, S. A. (1972). "Rotating Black Holes: Locally Nonrotating Frames, Energy Extraction, and Scalar Synchrotron Radiation." *Astrophysical Journal* **178**, 347. Bibcode: 1972ApJ...178..347B. Source for the Kerr ISCO and efficiency figures quoted above. - Shakura, N. I., and Sunyaev, R. A. (1973). "Black holes in binary systems. Observational appearance." *Astronomy & Astrophysics* **24**, 337–355. Bibcode: 1973A&A....24..337S. - Luminet, J.-P. (1979). "Image of a spherical black hole with thin accretion disk." *Astronomy & Astrophysics* **75**, 228–235. Bibcode: 1979A&A....75..228L. - Oppenheimer, J. R., and Snyder, H. (1939). "On Continued Gravitational Contraction." *Physical Review* **56**, 455. DOI: 10.1103/PhysRev.56.455. - Finkelstein, D. (1958). "Past-Future Asymmetry of the Gravitational Field of a Point Particle." *Physical Review* **110**, 965. DOI: 10.1103/PhysRev.110.965. - Penrose, R. (1965). "Gravitational Collapse and Space-Time Singularities." *Physical Review Letters* **14**, 57. DOI: 10.1103/PhysRevLett.14.57. - Hawking, S. W. (1975). "Particle creation by black holes." *Communications in Mathematical Physics* **43**, 199–220. DOI: 10.1007/BF02345020. - Event Horizon Telescope Collaboration (2019). "First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole." *Astrophysical Journal Letters* **875**, L1. DOI: 10.3847/2041-8213/ab0ec7. - Event Horizon Telescope Collaboration (2019). "First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole." *Astrophysical Journal Letters* **875**, L6. DOI: 10.3847/2041-8213/ab1141. Source of the 42 ± 3 μas ring diameter and the (6.5 ± 0.7) × 10⁹ M<sub>☉</sub> mass. - Event Horizon Telescope Collaboration (2022). "First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way." *Astrophysical Journal Letters* **930**, L12. DOI: 10.3847/2041-8213/ac6674. Source of the 51.8 ± 2.3 μas ring diameter. - GRAVITY Collaboration (2021). "Mass distribution in the Galactic Center based on interferometric astrometry of multiple stellar orbits." *Astronomy & Astrophysics* **657**, L12. Bibcode: 2021A&A...657L..12G. Source of *M* = 4.297 × 10⁶ M<sub>☉</sub> and *R*₀ = 8277 pc for Sgr A*. - GRAVITY Collaboration (2020). "Detection of the Schwarzschild precession in the orbit of the star S2 near the Galactic centre massive black hole." *Astronomy & Astrophysics* **636**, L5. Bibcode: 2020A&A...636L...5G. - Miller-Jones, J. C. A., et al. (2021). "Cygnus X-1 contains a 21-solar mass black hole — Implications for massive star winds." *Science* **371**, 1046–1049. DOI: 10.1126/science.abb3363. - Abbott, B. P., et al. (LIGO Scientific Collaboration and Virgo Collaboration) (2016). "Observation of Gravitational Waves from a Binary Black Hole Merger." *Physical Review Letters* **116**, 061102. DOI: 10.1103/PhysRevLett.116.061102. - Prša, A., et al. (2016). "Nominal Values for Selected Solar and Planetary Quantities: IAU 2015 Resolution B3." *Astronomical Journal* **152**, 41. Source of the nominal solar mass parameter *GM*<sub>☉</sub> = 1.3271244 × 10²⁰ m³ s⁻² and radius used throughout the HUD. **On the spine:** [[Gravitational_lens]] · [[Spacetime]] · [[Stellar_structure]] · [[General_relativity]] · [[Gravity]] · [[Photon]] · [[Physics]]. ## Wikipedia : Wikitube **Strict pair:** [Wikipedia](https://en.wikipedia.org/wiki/Black_hole) : [Wikitube](https://en.wikitube.io/wiki/Black_hole) --- *PORTAL_Physics three.js batch · 2026-08-05 · sim staged in `_3d_deploy_stage/`.*