28.08.01 · astronomy / high-energy-astrophysics

High-energy astrophysics — compact objects, accretion, and cosmic explosions

shipped3 tiersLean: none

Anchor (Master): Longair, High Energy Astrophysics (3e), full treatment; Shapiro & Teukolsky, Black Holes, White Dwarfs, and Neutron Stars (1983), Ch. 2-4, 8-10, 14; Frank, King & Raine, Accretion Power in Astrophysics

Intuition Beginner

Most of the light our eyes can see comes from stars: gentle, warm, and steady. But the universe also produces far more violent radiation, in X-rays and gamma rays that are invisible to us and blocked by Earth's atmosphere. The objects that make this high-energy radiation are the most extreme things in nature: the corpses of dead stars.

When a star runs out of fuel, gravity crushes what is left into a compact object. A star like the Sun becomes a white dwarf, a glowing ember the size of Earth but with the mass of a star. More massive stars collapse harder, crushing their cores into neutron stars: objects with the mass of the Sun squeezed into a sphere only about 20 kilometres across. The heaviest stars collapse all the way into black holes, where matter is crushed out of existence.

These objects are extreme in every way. A teaspoon of neutron-star material would weigh about six billion tonnes. Their magnetic fields can be a trillion times stronger than Earth's. Some spin hundreds of times per second. Around them, gas spiralling inward is heated to millions of degrees, glowing in X-rays. This process, called accretion, is the most efficient energy source known, far more efficient than nuclear fusion.

Gravity is the engine. When gas falls onto a compact object, gravity releases enormous energy, just as a falling rock picks up speed. Because a neutron star or black hole is so small and dense, gas can fall incredibly deep into its gravitational well, converting up to about forty percent of its own rest-mass energy into radiation and outflows. This is why accreting black holes outshine entire galaxies.

Rotation provides a second engine. A spinning magnetised neutron star acts like a colossal dynamo. Its rotating magnetic field whips charged particles around at nearly the speed of light, sweeping beams of radio waves across space like a lighthouse. When such a beam crosses Earth, we see a pulse: a pulsar. The fastest pulsars spin hundreds of times each second and keep time more accurately than atomic clocks.

The biggest black holes sit at the centres of galaxies and can weigh as much as billions of Suns. When they accrete gas, they blaze as quasars, the brightest steady objects in the universe, and fire narrow jets of plasma outward at nearly the speed of light. Even more violent are gamma-ray bursts, brief flashes that in a few seconds release more energy than the Sun will emit in its entire ten-billion-year life.

High-energy astrophysics is the study of these extremes. It tells us what happens at the edge of what physics allows: the densest matter, the strongest gravity, the largest magnetic fields, and the most powerful explosions. These objects are rare, but they shape galaxies, seed space with heavy elements, and let us test physics that no laboratory on Earth could ever reproduce.

Visual Beginner

The compact-object family and its power sources, plotted by density and accretion luminosity.

Object Mass (solar) Radius Density Power source
White dwarf up to 1.4 Earth-sized (~10,000 km) ~tonne per cm³ Residual heat; accretion (novae)
Neutron star 1.4 to ~2.3 ~10-20 km Nuclear (~10¹⁷ kg/m³) Rotation; accretion (X-ray binaries)
Stellar black hole 3 to ~50 event horizon ~3 km per solar mass (collapsed) Accretion; jets
Supermassive black hole 10⁶ to 10¹⁰ up to size of Solar System (collapsed) Accretion (quasars, AGN)

The key contrast is size. The Sun has roughly the same mass as a neutron star, but the neutron star is a hundred thousand times smaller. Falling that much deeper into gravity is what releases the enormous energies of high-energy astrophysics.

Worked example Beginner

The Eddington luminosity is the largest luminosity an object can have before its own radiation blows away the gas falling onto it. Beyond this limit, the outward push of light overwhelms gravity and accretion stops. It sets the maximum steady brightness of an accreting star.

The formula is

where is the object's mass and is the Sun's mass. So the Eddington limit scales directly with mass: heavier objects can shine brighter before blowing their food away.

Take a neutron star of mass . Its Eddington luminosity is

That is about times the Sun's luminosity (the Sun emits W). A single neutron star, accreting at its limit, can outshine fifty thousand Suns in X-rays.

A second useful number is the Schwarzschild radius, the distance from a black hole's centre to its event horizon (the point of no return):

For a black hole of , the event horizon sits at km — a sphere smaller than a city, into which ten Suns have vanished. Gas spiralling through a disc down to this horizon converts roughly ten percent of its mass into energy, which is why stellar-mass black holes in X-ray binaries are so bright.

Check your understanding Beginner

Formal definition Intermediate+

Compact objects

A compact object is the end product of stellar evolution in which matter has been compressed to a density comparable to or exceeding atomic (white dwarfs), nuclear (neutron stars), or infinite (black holes) density. Their defining property is compactness, the dimensionless ratio of Schwarzschild radius to actual radius:

For a white dwarf ; for a neutron star ; for a black hole at the horizon . As grows, General Relativity (GR) becomes progressively more important for the object's structure and the radiation it produces.

A white dwarf is supported against collapse by electron degeneracy pressure. Its maximum mass, derived by Chandrasekhar in 1931 [Chandrasekhar1931], is

where is the mean molecular weight per electron and . A neutron star is supported by neutron degeneracy pressure and nuclear forces. Its maximum mass, derived by Oppenheimer and Volkoff in 1939 using the relativistic TOV equation [OppenheimerVolkoff1939], lies between roughly and , depending on the unknown nuclear equation of state. Above this limit no known pressure can resist collapse, and a black hole forms.

A non-rotating black hole is described by the Schwarzschild metric with event horizon at the Schwarzschild radius ; a rotating one by the Kerr metric with horizon radius , where is the spin parameter.

Accretion and its efficiency

Accretion is the gravitationally driven capture of gas by a compact object. The luminosity released is a fraction of the mass-energy accreted:

The radiative efficiency equals the binding energy of gas at the inner edge of the disc, expressed as a fraction of rest-mass energy. For a Newtonian disc reaching the surface of a neutron star ( km, ), (about ten percent). For a Schwarzschild black hole, whose innermost stable circular orbit (ISCO) lies at , the standard value is . For a maximally spinning (extremal) Kerr black hole, rises to about . Compare these with nuclear fusion of hydrogen to helium, which releases only of rest-mass energy: accretion onto compact objects is up to sixty times more efficient.

The high-energy spectrum

The temperature of the inner accretion disc follows from equating its radiated flux to a blackbody. For a steady-state thin disc, the effective temperature at radius scales as

so is largest at small . Around a neutron star or stellar-mass black hole this temperature reaches K, peaking in X-rays (Wien's law: nm for K). Around a supermassive black hole the ISCO is much larger, so the inner disc is cooler (often ultraviolet), and most observed X-rays are produced by a hot corona Compton-upscattering disc photons.

Key derivation: accretion and the Eddington limit Intermediate+

The Eddington luminosity

Consider ionised gas of hydrogen plasma accreting onto a compact object of mass . Each electron-proton pair is held to the object by gravity (acting mainly on the massive proton) but pushed outward by Thomson scattering of radiation off the electron. The Eddington luminosity is the value of at which these two forces exactly balance for an optically thin plasma.

Outward radiation force on one electron: the Thomson cross-section is , the radiation flux at radius is , and the momentum flux absorbed is . Inward gravitational force on the attached proton is . Equating,

the radius cancels, giving

In solar units, . Above the outward force wins, radiation drives a wind, and steady spherical accretion is choked off. The Eddington limit is therefore the natural scale for the brightest steady sources in the universe. Sub-Eddington accretion produces a cool, geometrically thin disc; super-Eddington accretion (seen in some ultraluminous X-ray sources and gamma-ray bursts) drives massive outflows and thick radiation-dominated funnels.

Why accretion beats fusion

Combining with the Eddington ceiling fixes a maximum accretion rate that can be radiated steadily,

about /yr per solar mass for . Over a billion years, a supermassive black hole can swallow enough gas at this rate to reach billions of solar masses, which is the basic argument for how quasars grow.

Pulsars as rotating magnetic dipoles

A spinning magnetised neutron star of moment of inertia kg m² and angular velocity radiates away rotational energy through magnetic dipole radiation. The spin-down luminosity is

with the stellar radius, the polar surface field, and the angle between magnetic and rotation axes. The observed slowdown of a pulsar fixes its surface field, typically G for ordinary pulsars and up to G for magnetars. The total rotational energy reservoir for the Crab pulsar ( ms) is about J, of which roughly W is currently being radiated — comparable to the Eddington luminosity of a object.

Bridge. This Eddington argument builds toward the entire theory of active galactic nuclei and X-ray binaries, and it appears again in the cosmology of quasar growth and black-hole seeding 28.04.01; the central insight is that the brightest steady sources in the universe are bounded by a single mass-dependent luminosity, this is exactly why quasar luminosity functions map onto black-hole mass functions, and the foundational reason that accretion, not fusion, powers the high-energy sky. Putting these together, the same formula governs a neutron star and a supermassive black hole, and the bridge is that gravitational compactness sets both the energy budget and the brightness ceiling of every high-energy source.

Exercises Intermediate+

Lean formalization Intermediate+

This unit has no Lean formalization (lean_status: none). High-energy astrophysics rests on the relativistic structure equations of compact objects (Tolman-Oppenheimer-Volkoff for neutron stars, the Kerr metric for rotating black holes), on radiative transfer through optically thick and thin plasmas, and on magnetohydrodynamic disc and jet solutions. Mathlib currently provides neither a general-relativistic fluid equilibrium nor radiation magnetohydrodynamics, and the physical constants needed here (Thomson cross-section, specific nuclear equations of state) are only partially in the library. The Eddington luminosity itself is a one-line dimensional identity that could in principle be stated once the supporting GR infrastructure exists, but the structure equations, the ISCO, and the spin-down formula remain out of reach. Until that infrastructure is built, this unit stays a prose-and-calculation treatment validated against the primary literature.

Advanced results Master

Neutron-star structure and the equation of state

A neutron star in hydrostatic equilibrium in General Relativity obeys the Tolman-Oppenheimer-Volkoff (TOV) equations, the relativistic generalisation of Newtonian hydrostatic balance:

with the enclosed gravitational mass. The TOV equations stiffen the pressure gradient relative to the Newtonian case, lowering the maximum mass below the naive Chandrasekhar-style estimate. The maximum mass depends sensitively on the equation of state (EOS) of matter above nuclear saturation density ( kg/m³), which is unknown because it cannot be probed in terrestrial laboratories.

Two observational constraints bracket the EOS. The pulsars PSR J1614-2230 and PSR J0348+0432 demand an EOS stiff enough to support at least two solar masses, ruling out many exotic "soft" proposals (free quarks, kaon condensates). The gravitational-wave event GW170817, a neutron-star merger observed by LIGO/Virgo, constrained the tidal deformability , favouring radii of about km. The combination of a high maximum mass and a relatively small radius leaves a narrow window for the EOS, and ongoing NICER X-ray pulse-profile measurements of pulsars like J0030+0451 and J0740+6620 further shrink it.

Black-hole accretion discs

The standard thin-disc model (Shakura & Sunyaev 1973) solves the steady-state viscous relativistic fluid equations on circular Keplerian orbits, with "viscosity" parameterised by representing magnetorotational turbulence. It predicts a multi-temperature blackbody spectrum peaking in the band set by , a surface temperature profile in the outer disc, and most of the luminosity emerging from the innermost few gravitational radii. Deviations from the thin-disc model classify observationally distinct states: a hard state dominated by a hot ( K) Comptonising corona, a soft state dominated by the thermal disc, and a steep power-law state.

A black hole's spin is measured by the ISCO location. The relativistic iron K fluorescent line at keV, produced in the inner disc and redshifted and broadened by GR, has a skewed profile whose red wing directly encodes and hence the dimensionless spin . Measurements of the stellar-mass black hole in Cygnus X-1 give , near maximal, confirming that accretion can spin up black holes to extremality and raise the radiative efficiency toward .

Relativistic jets and the Blandford-Znajek mechanism

The narrow, light-year-long jets of active galactic nuclei and the smaller jets of some X-ray binaries move outward at bulk Lorentz factors . Their composition (electron-proton plasma versus electron-positron pairs) and launching mechanism remain open problems, but two leading extraction processes dominate theory. Blandford-Znajek extracts rotational energy from the ergosphere of a spinning black hole through magnetic field lines threading the horizon, functioning as a colossal circuit with the hole's rotation as the EMF. Blandford-Payne extracts binding energy from the accretion disc via magnetocentrifugal acceleration along field lines anchored in the disc. Both require large-scale ordered magnetic flux and reproduce the observed correlation between jet power and accretion luminosity.

The apparent superluminal motion of jet knots is a geometric projection effect: a blob moving at Lorentz factor toward the observer at small angle has an apparent transverse speed , which can exceed without violating relativity. This effect, combined with Doppler beaming (the boosting of emission along the direction of motion by a factor ), explains why gamma-ray bursts and blazars are visible across the observable universe.

Gamma-ray bursts

Gamma-ray bursts (GRBs) divide into long ( s) and short ( s) classes. Long GRBs are tied to the collapsar model: the core of a massive, rapidly rotating Wolf-Rayet star collapses to a black hole, accreting through a debris torus at rates of /s, far above Eddington. A bipolar relativistic jet punches through the stellar envelope, and its internal shocks produce the prompt gamma-ray emission while the resulting supernova is sometimes visible (e.g. SN 1998bw). Short GRBs are produced by neutron-star mergers (confirmed by the gravitational-wave event GW170817 and its kilonova counterpart AT2017gfo), the same events that forge the heaviest r-process elements.

The energy budget of a long GRB, correcting for jet opening angles of a few degrees, is of order J, a sizeable fraction of a supernova's total output concentrated into a few seconds and channelled into a narrow cone. The afterglow, seen across the electromagnetic spectrum from X-rays to radio, arises as the decelerating jet sweeps up external medium and produces synchrotron radiation, allowing the burst's energy, ambient density, and jet opening angle to be reconstructed.

Synthesis. High-energy astrophysics builds toward a unified picture in which the Eddington luminosity, accretion efficiency, and relativistic gravity together explain objects from city-sized neutron stars to billion-solar-mass quasars; the foundational reason is that gravitational compactness is the single dimensionless parameter setting every energy scale, this is exactly why the same formula bounds a pulsar and a AGN, the central insight is that rotation and magnetic field convert gravitational binding energy into observable radiation through accretion discs and pulsar magnetospheres, the framework generalises from steady X-ray binaries to transient collapsars and merging neutron stars, appears again in the cosmological growth of supermassive black holes 28.04.01 and in the structure of galaxies hosting AGN 28.03.03 pending, and putting these together the bridge is that every high-energy source is a gravitationally powered engine whose luminosity, spectrum, and variability are set by the mass, spin, and accretion rate of its central compact object.

Full proof set Master

Proposition (Eddington luminosity is the maximum for steady spherical accretion)

Statement. Consider spherically symmetric, steady, optically thin accretion of fully ionised hydrogen onto a point mass , with the radiation field produced by the accretion itself carrying luminosity isotropically. Then steady inflow is possible only if .

Proof. In steady state each electron-proton pair moves radially. By spherical symmetry and the assumption of optically thin coupling, the radiation flux at radius is and is directed radially outward. A single electron absorbs momentum from this flux at a rate

directed radially outward. The attached proton feels the Newtonian gravitational attraction of the central mass,

directed radially inward. The electron's negligible mass means the pair is dynamically coupled: any net force on the pair is dominated by these two terms. For steady inflow at constant (or non-accelerating) sub-relativistic speed the net radial force must be non-positive (inward), i.e. . Substituting,

Both denominators equal and cancel, leaving a radius-independent bound,

If the net force is outward, the gas is accelerated away from the source, and steady spherical inflow cannot be maintained; a radiation-driven wind forms instead. Hence is the maximum luminosity compatible with steady spherical accretion of pure ionised hydrogen.

Remark. Three assumptions are load-bearing: pure hydrogen composition (a mean molecular weight per electron ; heavier donors raise by ); geometrically isotropic emission (beaming lowers the true Eddington ceiling in the beam); and optically thin coupling via Thomson scattering (at very high optical depth, radiation is advected and the simple balance fails, allowing super-Eddington sources such as GRB central engines).

Proposition (Accretion radiative efficiency at the ISCO)

Statement. In Newtonian gravity, the radiative efficiency of a thin accretion disc extending to an inner radius around a compact object of mass is .

Proof sketch. A gas parcel of unit mass in a Keplerian orbit at radius has specific binding energy , the standard virial result: half the gravitational potential energy appears as orbital kinetic energy, and the other half, released by viscous torques as the parcel drifts inward, is radiated. Starting at large radius (zero binding energy) and drifting to , the parcel radiates an energy

per unit mass. The rest-mass energy per unit mass is , so the radiative efficiency is

For the Schwarzschild ISCO, , giving ; for a neutron-star surface km with , ; and for an extremal Kerr ISCO () at the Newtonian level (full GR gives after accounting for the last stable orbit's specific energy).

Connections Master

  • Stars and stellar endpoints 28.02.01. Every compact object is the remnant of a star, and the mass thresholds set in stellar evolution 28.02.01 fix which remnant forms: cores below the Chandrasekhar mass become white dwarfs, cores between roughly and become neutron stars, and heavier cores become black holes. The bifurcation of stellar fates is the input condition for all of high-energy astrophysics.

  • Black holes in General Relativity 13.06.01. The Schwarzschild and Kerr metrics studied in the GR unit 13.06.01 supply the spacetime geometry in which all black-hole accretion occurs. The event horizon, the ISCO, the ergosphere, and the no-hair theorem are the geometric inputs on which the Eddington argument, the iron-line spin diagnostics, and the Blandford-Znajek jet mechanism all depend.

  • Galaxies and active galactic nuclei 28.03.03 pending. Supermassive black holes and their accretion discs are the central engines of AGN and quasars 28.03.03 pending, and the observed correlations between black-hole mass and host-galaxy bulge properties imply that AGN feedback co-evolves with galaxies, regulating star formation and shaping the galaxy luminosity function.

  • Cosmology and black-hole growth 28.04.01. The growth of supermassive black holes across cosmic time, the luminosity function of quasars, and the role of AGN feedback in reheating the intergalactic medium are framed by cosmology 28.04.01; the Eddington limit sets the minimum time to grow a billion-solar-mass quasar by redshift , the classic argument for massive seed black holes in the early universe.

  • Gravitational-wave and multi-messenger astronomy. Mergers of neutron-star and black-hole binaries, detected as gravitational waves by LIGO/Virgo/KAGRA and in the electromagnetic band for GW170817, are simultaneously high-energy astrophysical sources and precision tests of the TOV equation of state, connecting this unit to relativistic multi-messenger astronomy and the r-process nucleosynthesis of heavy elements.

Historical & philosophical context Master

The history of high-energy astrophysics is a history of objects predicted long before they were seen, and of radiation bands that required leaving the atmosphere. The Chandrasekhar limit was derived on the boat trip from India to Cambridge in 1930 [Chandrasekhar1931]; Eddington's public rejection of the result at the 1935 Paris meeting — that there ought to be a law of nature "to prevent a star from behaving in this absurd way" — set the field back, but the conclusion stood. Oppenheimer and Volkoff's 1939 paper [OppenheimerVolkoff1939] showed that General Relativity places an even lower ceiling on neutron-star masses, and Oppenheimer with Snyder the same year showed that continued collapse produces what Wheeler would later christen the black hole.

The observational side opened in 1949 when radio astronomers identified the first cosmic radio sources, but the decisive step into the high-energy universe came on 12 June 1962, when Riccardo Giacconi's Aerobee rocket, launched to look for X-rays from the Moon, instead discovered a far brighter source in Scorpius [Giacconi1962]. Scorpius X-1 was the first extrasolar X-ray source, and its optical counterpart proved to be a faint blue star, implying an X-ray-to-optical luminosity ratio of about a thousand — direct evidence that accretion onto a compact object was the mechanism. Giacconi received the 2002 Nobel Prize for this work.

The serendipitous discovery of pulsars followed in July 1967, when Jocelyn Bell Burnell, a graduate student in Antony Hewish's group at Cambridge, noticed a regularly repeating signal in radio data at 81.5 MHz with a period of 1.337 seconds [HewishBell1968]. The precision of the signal earned it the half-serious designation LGM-1 ("Little Green Men"). Thomas Gold's 1968 identification of pulsars as rotating magnetised neutron stars vindicated the thirty-year-old theoretical prediction of Baade and Zwicky (1934) and opened pulsar timing as a discipline; Hewish received the 1974 Nobel Prize, while Bell Burnell's exclusion remains a long-standing controversy.

Quasars were identified as cosmologically distant, superluminous sources in the early 1960s through Maarten Schmidt's decoding of the redshifted Balmer lines of 3C 273, which placed it at a recession velocity of nearly . Their enormous luminosities were explained by Edwin Salpeter (1964) and Yakov Zel'dovich (1964) as accretion onto supermassive black holes, and Donald Lynden-Bell (1969) argued that such black holes should reside in the nuclei of essentially all massive galaxies — a prediction confirmed decades later by stellar-dynamical measurements of the Galactic Centre (the 2020 Nobel-recognised work of Genzel and Ghez on Sgr A*).

The last quarter of the twentieth century added relativistic jets, the gamma-ray-burst phenomenon (discovered by the Vela satellites in the late 1960s and shown to be cosmological by BATSE in the 1990s), and X-ray timing that revealed kilohertz quasi-periodic oscillations from accreting neutron stars. The twenty-first century closed several loops: NICER measured pulsar radii, LIGO directly detected merging neutron stars and black holes, and the Event Horizon Telescope imaged the shadow of the supermassive black hole in M87 and in the Galactic Centre.

Philosophically, high-energy astrophysics forces a confrontation with the limits of physical law. Compact objects are the only places in the universe where General Relativity, quantum mechanics at nuclear density, and electrodynamics at super-strong fields all matter at once. They are the cleanest tests of GR in the strong-field regime and the strongest constraints on matter beyond nuclear density. The existence of black holes — regions from which no information can return — also sharpens long-standing questions about determinism, information, and the role of the observer, and the resolution of Hawking's information paradox remains one of the deepest open problems in theoretical physics.

Bibliography Master

@article{Chandrasekhar1931,
  author  = {Chandrasekhar, S.},
  title   = {The Maximum Mass of Ideal White Dwarfs},
  journal = {Astrophysical Journal},
  volume  = {74},
  pages   = {81--82},
  year    = {1931},
  doi     = {10.1086/143316}
}

@article{OppenheimerVolkoff1939,
  author  = {Oppenheimer, J. R. and Volkoff, G. M.},
  title   = {On Massive Neutron Cores},
  journal = {Physical Review},
  volume  = {55},
  pages   = {374--381},
  year    = {1939},
  doi     = {10.1103/PhysRev.55.374}
}

@article{BaadeZwicky1934,
  author  = {Baade, W. and Zwicky, F.},
  title   = {On Super-novae},
  journal = {Proceedings of the National Academy of Sciences},
  volume  = {20},
  pages   = {254--259},
  year    = {1934}
}

@article{Giacconi1962,
  author  = {Giacconi, R. and Gursky, H. and Paolini, F. R. and Rossi, B. B.},
  title   = {Evidence for X-rays From Sources Outside the Solar System},
  journal = {Physical Review Letters},
  volume  = {9},
  pages   = {439--443},
  year    = {1962},
  doi     = {10.1103/PhysRevLett.9.439}
}

@article{HewishBell1968,
  author  = {Hewish, A. and Bell, S. J. and Pilkington, J. D. H. and Scott, P. F. and Collins, R. A.},
  title   = {Observation of a Rapidly Pulsating Radio Source},
  journal = {Nature},
  volume  = {217},
  pages   = {709--713},
  year    = {1968},
  doi     = {10.1038/217709a0}
}

@article{Gold1968,
  author  = {Gold, T.},
  title   = {Rotating Neutron Stars as the Origin of the Pulsating Radio Sources},
  journal = {Nature},
  volume  = {218},
  pages   = {731--732},
  year    = {1968}
}

@article{Salpeter1964,
  author  = {Salpeter, E. E.},
  title   = {Accretion of Interstellar Matter by Massive Objects},
  journal = {Astrophysical Journal},
  volume  = {140},
  pages   = {796--800},
  year    = {1964}
}

@article{LyndenBell1969,
  author  = {Lynden-Bell, D.},
  title   = {Galactic Nuclei as Collapsed Old Quasars},
  journal = {Nature},
  volume  = {223},
  pages   = {690--694},
  year    = {1969}
}

@article{ShakuraSunyaev1973,
  author  = {Shakura, N. I. and Sunyaev, R. A.},
  title   = {Black holes in binary systems. Observational appearance},
  journal = {Astronomy \& Astrophysics},
  volume  = {24},
  pages   = {337--355},
  year    = {1973}
}

@article{BlandfordZnajek1977,
  author  = {Blandford, R. D. and Znajek, R. L.},
  title   = {Electromagnetic extraction of energy from Kerr black holes},
  journal = {Monthly Notices of the Royal Astronomical Society},
  volume  = {179},
  pages   = {433--456},
  year    = {1977}
}

@book{ShapiroTeukolsky1983,
  author    = {Shapiro, S. L. and Teukolsky, S. A.},
  title     = {Black Holes, White Dwarfs, and Neutron Stars: The Physics of Compact Objects},
  publisher = {Wiley-Interscience},
  year      = {1983}
}

@book{LongairHEA2011,
  author    = {Longair, M. S.},
  title     = {High Energy Astrophysics},
  edition   = {3rd},
  publisher = {Cambridge University Press},
  year      = {2011}
}

@book{FrankKingRaine2002,
  author    = {Frank, J. and King, A. and Raine, D. J.},
  title     = {Accretion Power in Astrophysics},
  edition   = {3rd},
  publisher = {Cambridge University Press},
  year      = {2002}
}

@book{CarrollOstlie2017,
  author    = {Carroll, B. W. and Ostlie, D. A.},
  title     = {An Introduction to Modern Astrophysics},
  edition   = {2nd},
  publisher = {Cambridge University Press},
  year      = {2017}
}

@article{AbbottGW1708172017,
  author  = {{LIGO Scientific Collaboration and Virgo Collaboration} and {Fermi-GBM} and {INTEGRAL}},
  title   = {Multi-messenger Observations of a Binary Neutron Star Merger},
  journal = {Astrophysical Journal Letters},
  volume  = {848},
  pages   = {L12},
  year    = {2017}
}

@article{EventHorizonTelescope2019,
  author  = {{Event Horizon Telescope Collaboration}},
  title   = {First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole},
  journal = {Astrophysical Journal Letters},
  volume  = {875},
  pages   = {L1},
  year    = {2019}
}