Black-hole accretion

Black-hole accretion is the inflow of matter onto a black hole, converting some of its gravitational binding energy into heat, electromagnetic radiation and outflows. Matter with substantial angular momentum generally forms an accretion disk; internal stresses transfer angular momentum outward so that matter can move inward. The released energy can be radiated before the material crosses the event horizon. In a thin-disk model the innermost stable circular orbit marks the transition toward rapid plunging motion, although the detailed inner flow depends on the accretion state.^[ref]

The relation $L=\eta \dot{M}{c}^2$ defines the radiative efficiency of black-hole accretion, where $L$ is emitted luminosity, $\dot{M}$ the mass accretion rate and $c$ the speed of light. A radiatively inefficient accretion flow carries much of its available energy inward or into outflows instead of emitting it, so luminosity alone need not measure the mass supply.

An oscillatory feature with finite width in a source's variability power spectrum, rather than a perfectly periodic clock. Orbital, precession and disc-mode interpretations connect different frequencies to the same compact object. Identifying an observed peak with an innermost stable circular orbit orbital frequency is therefore a model assumption, not a direct black-hole mass or spin measurement.

Collisionless matter approaching a black hole generally cannot dissipate enough orbital energy to remain bound, so direct capture requires sufficiently small angular momentum. For nonrelativistic incident speed $v$, a Schwarzschild black hole has direct test-particle capture cross-section approximately $16\pi G^2{M}^2/(c^2{v}^2)$. A Newtonian pericentre estimate using the Schwarzschild radius instead gives $4\pi G^2{M}^2/(c^2{v}^2)$. Both show that capture at a given mass density is suppressed by order $(v/c{)}^2$ relative to supersonic Bondi--Hoyle--Lyttleton accretion. Tidal disruption events provide a distinct way to accrete some stellar debris.

For a singular isothermal sphere of dispersion $\sigma$ and gas fraction $f_g$, placing all gas inside radius $R$ onto a hole in time $R/\sigma$ gives ${\dot{M}}_{\mathrm{dyn}}\sim 2f_g{\sigma }^3/G$. This is an optimistic supply bound that omits angular momentum transport, star formation, and active-galactic-nucleus feedback. A numerical comparison with the Eddington accretion rate requires a gas fraction, reference radiative efficiency, opacity, and crossing-time convention.

The inflow time is the characteristic time for accreting matter to move radially by its current radius, $t_{\mathrm{in}}\sim R/|u_R|$.

The Paczyński-Wiita potential is a Newtonian pseudo-potential that reproduces the Schwarzschild innermost stable circular orbit at $r=3r_S$.

For $\beta =16c_s^2/c^2$, the exterior sonic point of the barotropic spherical accretion equation lies outside $r_G=2GM/c^2$. Its mass density is ${\rho }_c={\rho }_{\infty }\exp [4/(1+\sqrt{1+\beta })-1/2]$. The selected smooth transonic accretion rate is $4\pi r_c^2{\rho }_c{c}_s$. Subsonic formal solutions are not assigned this rate without transonic selection.

Let $r_S=2GM/c^2$. In the Paczyński-Wiita potential, circular orbit balance and the radial epicyclic frequency formula giveThus the innermost stable circular orbit lies at $r=3r_S$. Its specific orbital energy is $-c^2/16$, giving a zero-torque radiative efficiency of black-hole accretion of $1/16$. The model reproduces the Schwarzschild inner orbit exactly but approximates its relativistic energy and frequency. Abramowicz's derivation explains the potential's relation to the Schwarzschild effective potential.

The limiting inward departure from the marginal Paczyński-Wiita circular orbit at $r=3a$ conserves the displayed specific angular momentum and specific orbital energy. In the Paczyński-Wiita potential, its inward speed satisfies $v_r^2=GM(3a-r{)}^3/[4a{r}^2(r-a)]$. The exactly marginal circular orbit does not depart without a perturbation; the formula describes the separatrix limit of inward releases.

For a steady accretion disk with arbitrary rotation law and a zero-torque inner boundary condition at $R=3R_S$, the outward viscous torque in an accretion disk is $\mathcal{G}=\dot{M}(j-j_{\mathrm{in}})$. The dissipation per unit planar area summed over both disk faces is $D=-\mathcal{G}{\Omega }^{\prime }/(2\pi R)=\nu \Sigma (R{\Omega }^{\prime }{)}^2$. Integrating $2\pi RD$ gives $L=\dot{M}GM/(8R_S)=\dot{M}{c}^2/16$ in the Paczyński-Wiita potential. The radiative flux from one face is $D/2$, preventing a factor-of-two error in the luminosity.

The radiative efficiency is $\eta =L/(\dot{M}{c}^2)$. For a zero-torque thin disk it is $1-E_{\mathrm{ISCO}}$, the binding energy per unit rest mass released before matter crosses the innermost stable circular orbit.

The maximal radiative efficiency is the available binding energy per unit accreted rest-mass energy when all energy released outside the capture boundary escapes as radiation. A relativistic zero-torque disk has ${\eta }_{\max }=1-E_{\mathrm{ISCO}}$.

The innermost stable circular orbit is the smallest stable circular timelike orbit around a black hole. Its radius and orbital energy depend on black-hole spin and on whether the orbit is prograde or retrograde.

For equatorial test-particle circular orbits in the extremal-spin limit, the prograde Kerr black hole innermost stable circular orbit has Boyer-Lindquist radius tending to $GM/c^2$ and specific energy tending to $1/\sqrt{3}$. The retrograde limit has radius $9GM/c^2$ and specific energy $5/(3\sqrt{3})$. The prograde radius approaching the horizon coordinate does not make the massive limiting orbit identical to the null horizon generator.^[ref]

The Schwarzschild coordinate rotation frequency at the marginally stable circular orbit, measured per unit time of a stationary clock at infinity. It scales inversely with black-hole mass. A finite-radius static clock adds a gravitational time-dilation factor; a quasi-periodic brightness signal need not directly equal this orbital frequency.

For ideal zero-torque equatorial thin discs, efficiency is $1-E_{\mathrm{ISCO}}$. The formal extremal retrograde, Schwarzschild and extremal prograde limiting radii are respectively $9,6,1$ in units $GM/c^2$, giving efficiencies $1-\sqrt{25/27}$, $1-\sqrt{8/9}$ and $1-1/\sqrt{3}$. Photon capture and non-ideal accretion can change the realized efficiency.

Bondi accretion is steady spherical accretion from a uniform gas reservoir onto a compact gravitating object. In its adiabatic form, compression raises internal energy without providing the irreversible shear heating of a radiatively efficient disk.

The thermal sphere-of-influence scale $r_B=GM/c_s^2$ compares gravitational potential per unit mass with the ambient gas sound-speed squared. A convention based on escape speed instead introduces a factor of two. The scale is not a capture cross-section or an accretion rate without a flow and dissipation model.

With planet-to-star mass ratio $q$ and disk aspect ratio $h$, $r_B/a=q/h^2$ and $r_H/a=(q/3{)}^{1/3}$. Equating the two gives $q=h^3/\sqrt{3}$. The thermally limited radius grows faster with planet mass than the tidally limited radius; an envelope outer scale is often the smaller of the two.

For isentropic spherical accretion with this adiabatic index, define $\lambda =\dot{M}{c}_{\infty }^3/(4\pi G^2{M}^2{\rho }_{\infty })$ and $x=r{c}_{\infty }^2/(GM)$. The steady-flow equation is $\sqrt{\lambda }[{\mathcal{M}}^{3/2}/2+3{\mathcal{M}}^{-1/2}/2]=1+3x/2$. The bracket has minimum two at Mach one. Reaching arbitrarily small positive radii therefore requires $\lambda \le 1/4$. For smaller rates the outer-rest branch remains subsonic; at the maximal rate it approaches Mach one only as $r\to 0$, not at a finite sonic point. Supersonic outer branches fail the specified resting finite-density condition, and larger rates terminate at a positive minimum radius. An additional inner boundary prescription is needed to select a unique rate among the globally defined subsonic solutions.

Relativistic spherical accretion is steady radial gas inflow in Schwarzschild spacetime, with conserved rest-mass flux and relativistic Bernoulli function. A regular solution passes a finite sonic point outside the event horizon. For an isentropic gas with $\gamma =5/3$ and cold asymptotic sound speed $c_{\infty }\ll c$, the critical radius is approximately $3GM/(4c{c}_{\infty })$ and the leading mass accretion rate is $\pi G^2{M}^2{\rho }_{\infty }/c_{\infty }^3$. Relativistic corrections resolve the zero-radius critical point of the Newtonian $\gamma =5/3$ model without eliminating accretion.

Bondi--Hoyle--Lyttleton accretion is gravitational capture by a compact object moving through a gas. In the pressureless supersonic limit, ballistic streamlines focus into a downstream wake.

In a fixed central Newtonian gravitational potential, pressureless free fall has speed $u=a{R}^{-1/2}$, and randomly moving stars have typical gas-relative speed $V=b{R}^{-1/2}$. A steady, undepleted spherical mass flux $A$ gives $\rho =A{R}^{-3/2}/(4\pi a)$. The Bondi--Hoyle--Lyttleton accretion rate therefore becomes $\dot{M}=\lambda M^2$ with $\lambda =G^{2}A/(a{b}^3)$, independent of radius. Integrating $d(1/M)/dt=-\lambda$ gives the displayed solution before $\lambda M_{0}t=1$. Its finite-time divergence signals failure of the fixed potential, undepleted gas supply and negligible feedback assumptions; it is not a physically unlimited mass supply.

If stars are neither created nor lost and all follow quadratic accretion in a central Keplerian potential, conservation of number gives $F(M,t)dM=F_0(M_0)d{M}_0$. Inverting the trajectory gives $M_0=M/(1+\lambda tM)$ and its Jacobian determinant is $(1+\lambda tM{)}^{-2}$. This proves the displayed mass-density-in-mass-space transformation. In a broad high-mass range with $\lambda tM\gg 1$, whose inverse images remain in a narrow interval where $F_0$ varies little, the power law is $F\propto M^{-2}$. A strictly monochromatic initial population remains monochromatic. Finite initial support imposes a finite upper cutoff before the earliest formal blow-up, so the power law is an intermediate asymptotic subject to the accretion model's limits.

In a cold pressureless stream of asymptotic speed $v$, conservation of specific angular momentum and the eccentricity vector makes a streamline at impact parameter $d$ cross the downstream axis at $b=d^2{v}^2/(2GM)$. Its axial velocity there remains $v$, while opposing transverse velocities collide in a shock.

If transverse kinetic energy is dissipated at the downstream axis, the remaining specific orbital energy is $v^2/2-GM/b$. Binding requires $b<2GM/v^2$, corresponding to incident impact parameter $d<2GM/v^2$. Integrating the captured mass flux gives the Bondi--Hoyle--Lyttleton accretion rate.

In the ballistic supersonic limit, streamline collisions bind gas with impact parameter $b<2GM/v_{\infty }^2$, giving

Isothermal Bondi accretion has $p=c_s^2\rho$ with constant sound speed. Its regular transonic solution crosses the Bondi sonic point and approaches free fall near the accretor.

For steady Isothermal Bondi accretion from gas at rest at infinity, define $\lambda =\dot{M}{c}_s^3/(4\pi G^2{M}^2{\rho }_{\infty })$, $x=r/r_s$ and $y=u/c_s>0$, where $r_s=GM/(2c_s^2)$. The Bernoulli function and mass conservation give $y^2/2-\ln y=2\ln x+2/x-\ln (4\lambda )$. The left side has minimum $1/2$ at $y=1$, and the radius-dependent part on the right has minimum two at $x=1$. Existence over all positive radii therefore requires the displayed bound. Equality permits the smooth accelerating-inward transonic branch and gives the usual isothermal Bondi accretion rate. Smaller positive rates have an everywhere-subsonic reservoir-connected branch, and zero rate gives a hydrostatic atmosphere. An absorbing inner boundary inside the sonic radius selects the critical inflow; steady isothermal equations and the outer density alone do not specify that inner condition.

The regular spherical isothermal inflow from mass density ${\rho }_{\infty }$ at rest at infinity has sonic mass density $e^{3/2}{\rho }_{\infty }$. Combining it with the Bondi sonic point and mass conservation gives $\dot{M}=\pi e^{3/2}{G}^2{M}^2{\rho }_{\infty }/c_s^3$.

For inward speed $u>0$, steady spherical mass conservation and radial momentum balance give

With $x=r/r_s$ and $\mathcal{M}=|u_r|/c_s$, the radial isothermal equation integrates to ${\mathcal{M}}^2/2-\ln \mathcal{M}=2/x+2\ln x+C$. A regular transonic branch has $C=-3/2$, and the reservoir-fed accretion branch crosses $(x,\mathcal{M})=(1,1)$ with slope $-1$.

The isothermal Bondi sonic point is $r_s=GM/(2c_s^2)$, where the regular transonic solution has $u=c_s$.

At the critical point, both the coefficient of $du/dr$ and the remaining right-hand side of the Isothermal Bondi equation vanish, allowing a finite derivative through the sonic transition.

The mass accretion rate is the mass entering an accretor per unit time. In a steady spherical flow, $\dot{M}=4\pi r^2\rho u$ is independent of radius.

A steady full-sphere radial inflow has mass flux $\dot{M}=4\pi r^2\rho |v_r|$ by conservation of mass. With Newtonian free-fall speed from rest at infinity, $|v_r|=\sqrt{2GM/r}$, this becomes the displayed expression. The assumed velocity then requires $\rho \propto r^{-3/2}$ along the steady flow. A local density and speed do not determine the total mass accretion rate of a general nonspherical inflow, and this estimate does not impose the sonic boundary condition of Bondi accretion.

Consider steady spherical accretion of a gas obeying a polytropic equation of state in isentropic flow, with fixed background potential $\Phi =-A{r}^{-1/2}$, where $A=16\pi G{\Psi }_0$, with reservoir mass density ${\rho }_0$ and sound speed $c_0$ at infinity. At a regular sonic point, radial momentum balance requires $c_c^2=A/(4\sqrt{r_c})$. The Bernoulli function isSetting $u=c_s=c_c$ at the sonic point givesA finite positive critical sound speed therefore requires $1<\gamma <9/7$. Mass conservation gives $\dot{M}=4\pi r_c^2{\rho }_c{c}_c$, while the constant-entropy gas law gives ${\rho }_c/{\rho }_0=(c_c^2/c_0^2{)}^{1/(\gamma -1)}$. Since $r_c\propto A^2{c}_0^{-4}$ at fixed $\gamma$, these equations yield $\dot{M}\propto {\rho }_0{G}^4{\Psi }_0^4{c}_0^{-7}$. The steep sound-speed dependence comes from both the selected capture scale and the critical flow speed. The background potential is prescribed; gas self-gravity and finite outer boundaries require a different model.

For a background gravitating density ${\rho }_{\mathrm{dm}}={\Psi }_0{r}^{-5/2}$, the potential is $\Phi =-16\pi G{\Psi }_0{r}^{-1/2}$. A steady polytropic inflow has a regular sonic point only for $1<\gamma <9/7$, and its accretion rate scales as $\dot{M}\propto {\rho }_0{G}^4{\Psi }_0^4{c}_0^{-7}$.

The Newtonian free-fall speed from rest at infinity toward a point mass is $v_{\mathrm{ff}}=\sqrt{2GM/r}$.

For an adiabatic gas with $\gamma =5/3$, asymptotic density ${\rho }_0$, and sound speed $v_{s0}$ accreting onto a Newtonian point mass $M$, the largest globally subsonic steady rate is

A radiatively inefficient accretion flow carries a substantial fraction of locally generated energy with the gas instead of radiating it before accretion. It can occur at low density because emission is weak or at high optical depth because photons are trapped.

An advection-dominated accretion flow is a hot, geometrically thick, radiatively inefficient flow in which most dissipated energy remains in the gas and is carried inward rather than radiated locally.

A two-temperature accretion flow has ions and electrons at substantially different temperatures because their Coulomb equilibration time exceeds the inflow time. Ions can retain most viscous heat while electrons, which radiate more effectively, remain cooler.

The electron--proton thermal equilibration time is the timescale on which Coulomb collisions exchange energy between electron and proton populations. Its inverse dependence on number density helps make dilute accretion flows two-temperature and radiatively inefficient.

A slim accretion disk is an optically thick, moderately geometrically thick flow at a high accretion rate. Its inflow time can be shorter than its photon-diffusion time, so radial advection carries trapped radiation inward.

A slim accretion disk is optically thick and retains an important inward advective heat flux. It is distinct from the optically thin version of a radiatively inefficient accretion flow.

Photon trapping occurs when the radiative diffusion time across an optically thick flow exceeds its inward advection time. The trapped radiation is carried inward, reducing the emergent radiative efficiency.

The photon-trapping radius is where the outward radiative diffusion time equals the inward inflow time. Radiation generated inside it is predominantly advected toward the accretor.

A luminous hot accretion flow has radiative cooling greater than local viscous heating, so compressive inward advection supplies the difference. In the convention $q_{\mathrm{adv}}=q^{+}-q^{-}$, it has $q_{\mathrm{adv}}<0$.

The accretion-flow advection balance compares local volumetric heating $q^{+}$, radiative cooling $q^{-}$, and energy transport with the gas. Positive $q_{\mathrm{adv}}$ means advection removes locally generated heat; zero gives local radiative balance; negative means advection supplies heat that is radiated locally.

The Eddington ratio is the luminosity divided by the Eddington luminosity. At fixed radiative efficiency it also measures the accretion rate relative to the Eddington rate.

An Eddington accretion rate converts the Eddington luminosity into a mass supply rate using an adopted reference efficiency ${\eta }_0$: ${\dot{M}}_{\mathrm{Edd}}=L_{\mathrm{Edd}}/({\eta }_0{c}^2)$. The common choice ${\eta }_0=0.1$ gives ${\dot{M}}_{\mathrm{Edd}}=10L_{\mathrm{Edd}}/c^2$.

For radiatively efficient release of Newtonian surface binding energy, $L\sim GM\dot{M}/R_{\ast }$ and $L_E=4\pi GMc/\kappa$ imply the displayed rate. Its efficiency is $\eta \sim GM/(R_{\ast }{c}^2)$, so at fixed opacity the maximum rate scales with stellar radius. This differs from defining an Eddington supply rate with an arbitrarily fixed reference efficiency. Geometry, advection and mass loss can modify the simple limit.

The Salpeter time is the e-folding time for Eddington-limited black-hole growth. If $t_E={\sigma }_{T}c/(4\pi G{m}_p)$, then $t_{\mathrm{Sal}}=\eta t_E/[(1-\eta )f_{\mathrm{Edd}}]$.

For fixed efficiency, growth from $M_i$ to $M_f$ within $\Delta t$ at no more than the Eddington luminosity requires $[\eta /(1-\eta )]t_E\log (M_f/M_i)\le \Delta t$. Solving gives an upper bound, because greater radiative efficiency reduces retained mass growth at a fixed luminosity. A duty cycle below unity tightens the requirement; mergers and super-Eddington episodes invalidate its stated premises.

At constant radiative efficiency of black-hole accretion $\eta$ and the Eddington luminosity, a black hole grows exponentially with the Salpeter time $t_S=\eta t_E/(1-\eta )$, where $t_E={\sigma }_{T}c/(4\pi G{m}_p)$. Its doubling time is $t_S\log 2$. The factor $1-\eta$ accounts for the supplied rest energy that escapes as radiation instead of increasing the hole's mass.