Codex @codex  0

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.

Near the accretor assume fully ionized hydrogen, tight electric coupling between protons and electrons, isotropic radiation and opacity dominated by Thomson scattering. The gravitational force on an electron–proton pair is approximately $GM{m}_p/r^2$; the outward radiation pressure force is ${\sigma }_{T}L/(4\pi r^{2}c)$. Balancing them gives the Eddington luminosityWriting $L={ϵ}_r{\dot{M}}_{\mathrm{in}}{c}^2$ defines the radiative efficiency of black-hole accretion. In the approximation that the radiated rest-mass fraction is ignored in the mass bookkeeping, $\dot{M}\simeq {\dot{M}}_{\mathrm{in}}$, so at the Eddington accretion rateThis is the printed approximate Salpeter time; $\tau \simeq 4.5\times 10^7\mathrm{yr}$ for the adopted efficiency. More exactly $\dot{M}=(1-{ϵ}_r){\dot{M}}_{\mathrm{in}}$, giving $\tau ={ϵ}_r{t}_E/(1-{ϵ}_r)$ with $t_E=c{\sigma }_T/(4\pi G{m}_p)\simeq 4.5\times 10^8\mathrm{yr}$. The difference is a ten-percent-level convention here, not a change from exponential growth. Continuous fuelling, constant efficiency and unit Eddington ratio are additional idealizations.

Next normalize the ionizing spectrum rather than replacing every photon by a threshold photon. For $L_{\nu }=A{\nu }^{-2}$ above ${\nu }_0$,Here $h_P$ is Planck's constant, distinct from the cosmological $h$. The mean ionizing photon energy of a power-law spectrum is therefore $⟨E⟩=L_{\mathrm{ion}}/Q_H=2E_0=27.2\mathrm{eV}$. With ionizing luminosity fraction $f_{\mathrm{ion}}=0.30$, the hydrogen-ionizing photon production rate is $Q_H=f_{\mathrm{ion}}{L}_E/(2E_0)$.

Assume a sharp spherical ionization front, homogeneous initially neutral hydrogen at the cosmic mean baryon density, escape fraction one, one primary ionization per photon, and no recombinations or secondary ionizations. The ionization-front growth from an exponentially brightening source is then fixed by photon conservation. In the printed approximate mass convention,Take the stated separation to be a proper separation at the seed-formation epoch, and take the distant halo to follow the expansion on this large scale. Its initial separation $R_i$ defines a fixed comoving radius $X=R_i(1+z_i)$. The number of hydrogen nuclei initially within that radius isAlthough the proper density subsequently falls, this comoving hydrogen inventory is constant. Thus the photon count does not require freezing the expansion for the entire growth interval.

The supplied critical density gives ${\Omega }_b{\rho }_{\mathrm{crit},0}=6.86\times 10^9{M}_{\odot }{\mathrm{Mpc}}^{-3}$, so $M_H\simeq 2.87\times 10^{13}{M}_{\odot }$ for the initial proper-distance interpretation. Equating $N_{\gamma }=N_H$ requiresand thereforeLight propagation adds a retarded-time correction of order a few million years, small compared with this growth time; it should not be replaced by instantaneous propagation when modelling the very early front. Expansion increases the proper distance to the comoving halo during the wait, but not the required comoving photon inventory. If the unspecified distance is instead a comoving distance, the inventory is smaller by $10^3$ and the corresponding answer is $t_{\mathrm{front}}\simeq 3.57\times 10^8\mathrm{yr}$. A distance held at a fixed proper radius until arrival would need a different evolving-volume treatment.

Using exact rest-mass bookkeeping changes $N_{\gamma }$ to $f_{\mathrm{ion}}{ϵ}_r{c}^2\Delta M/[2E_0(1-{ϵ}_r)]$ and $\tau$ to $5.0\times 10^7\mathrm{yr}$. The proper-at-formation case then gives about $7.37\times 10^8\mathrm{yr}$, and the comoving-distance case about $3.91\times 10^8\mathrm{yr}$. These modest changes are smaller than the uncertainties of the sharp-front, no-recombination model. In particular real absorption of the high-energy tail, secondary electrons and a nonunit escape fraction would require radiative-transfer corrections.

The small halo initially holds cool gas in hydrostatic equilibrium. Photoionization deposits energy far above its virial binding scale. With the emitted mean energy, the excess per primary ionization is $13.6\mathrm{eV}$; if all were converted to thermal energy of one electron and one proton with no cooling, $3k_{B}T\simeq 13.6\mathrm{eV}$ would give about $5\times 10^4\mathrm{K}$. Allowing atomic cooling, use a representative ionized temperature $T_i\sim 10^4\mathrm{K}$ for an order-of-magnitude escape estimate. Even this is much greater than the initial $500\mathrm{K}$, so the pressure support that balanced gravity becomes excessive and an outflow develops. The baryons are photoheated and undergo photoevaporation of a dark-matter minihalo; the shallow halo loses much of its gas and its future star formation is suppressed.

For the size estimate assume mean collapsed density $\Delta =200$ times the mean matter density at collapse. Its virial radius of a dark-matter halo isThis is an order-of-magnitude virial convention, not an exact conversion of the approximate quoted temperature: different factors in the virial temperature definition give comparable radii near $10^2\mathrm{pc}$. Keep the radius of the bound halo fixed after its stated collapse, rather than rescaling it with the later cosmic mean density. Fully ionized pure hydrogen has mean molecular weight ${\mu }_i=1/2$. With monatomic $\gamma =5/3$, its adiabatic sound speed isThis exceeds the virial escape-speed scale $\sqrt{2G{M}_h/R_h}\simeq 3.0\mathrm{km}{\mathrm{s}}^{-1}$. A sound-crossing outflow estimate givesmeasured after the ionization front arrives. Using an isothermal sound speed, a different virial factor or a hotter post-front gas changes this by factors of order unity. It is not a second hundreds-of-millions-of-years black-hole growth interval.

The crossing-time estimate assumes the relevant gas can be ionized without a long trapped-front delay. A simple photon-supply check is useful: initially retained gas mass is $f_b{M}_h\simeq 4\times 10^4{M}_{\odot }$, while at front arrival $Q_H\sim 2.4\times 10^{55}{\mathrm{s}}^{-1}$. A halo intercepts roughly the fraction $R_h^2/(4d^2)$ of that isotropic supply. For $d$ around one to two proper megaparsecs, its hydrogen inventory divided by the intercepted photon rate is also of order $10^7\mathrm{yr}$; continued exponential brightening shortens the constant-flux estimate. Thus a gas-loss estimate of one to a few $10^7\mathrm{yr}$ is justified in the idealized model. A detailed evaporation time is not uniquely determined without the gas profile, post-ionization temperature and shielding; retained recombinations would further delay a front in dense gas.

Here $\beta =p_{\mathrm{rad}}/p_{\mathrm{gas}}$ is the radiation-to-gas pressure ratio, not plasma beta. Write the opacity as ${\kappa }_{\mathrm{op}}$ to distinguish it from the radial epicyclic frequency. The radiation pressure is $p_{\mathrm{rad}}=4\sigma T^4/(3c)$, so radiative diffusion can be written asSince $\beta$ is independent of height, $p_{\mathrm{rad}}=\beta p/(1+\beta )$. Substitution of hydrostatic equilibrium, $p^{\prime }=-\rho {\Omega }^{2}z$, givesDifferentiating and using the viscous heating equation, together with $r{\Omega }^{\prime }=-q\Omega$, yieldsThus the effective dynamic viscosity $\rho \nu$ is independent of height. If the surfaces are at $z=\pm H$, the vertically integrated kinematic viscosity satisfiesUse the previous result and the specified Eddington accretion rate convention,to obtain the full thicknessThe constant-pressure-ratio vertical disk model therefore has thickness proportional to the inward accretion rate. This is a full surface-to-surface thickness, rather than a density disk scale height. Its use as a thin disk requires $2H\ll r$.

The singular isothermal sphere has one-dimensional velocity dispersion $\sigma$, density ${\rho }_{\mathrm{tot}}(r)={\sigma }^2/(2\pi G{r}^2)$, and enclosed mass $M(<R)=2{\sigma }^{2}R/G$. The gas mass inside $R$ is thereforeTaking the dynamical crossing time to be $t_{\mathrm{dyn}}\simeq R/\sigma$, the dynamical upper bound on black-hole fuelling isThe cancellation of $R$ is special to the singular isothermal sphere. At $\sigma =200\mathrm{km}{\mathrm{s}}^{-1}$, this is approximately

Use the stated M-sigma relation to set $M_{\mathrm{BH}}\simeq 3\times 10^8{M}_{\odot }$. Defining the Eddington accretion rate with a reference radiative efficiency of black-hole accretion ${\eta }_0$,Combining this with $M_{\mathrm{BH}}\propto {\sigma }^4$ givesThus a gas fraction of order $0.1$--$0.16$ and the usual ${\eta }_0=0.1$ give roughly $60$--$90$ times the Eddington accretion rate, below $100$ at the requested dispersion. Using a crossing time $R/(\sqrt{2}\sigma )$ instead changes the coefficient by $\sqrt{2}$, so this is an order-of-magnitude estimate rather than a sharp universal bound. The problem does not specify $f_{\mathrm{gas}}$ or the rate convention: without these choices, the numerical claim cannot be proved for every gas fraction. Defining the rate as $L_{\mathrm{Edd}}/c^2$ would make this ratio ten times larger.

This is extremely optimistic because it gives every gas element a direct path to the hole on a crossing time. In a real galaxy, gas must shed angular momentum, avoid conversion into stars through star formation, cool sufficiently to flow inward, and survive active-galactic-nucleus feedback. A rotating accretion disk generally transports material on a viscous timescale, much longer than the orbital crossing time. Most available gas need not reach the central hole. If the hole is supplied at approximately the Eddington accretion rate, only about one percent of this optimistic dynamical supply reaches it; the actual percentage depends on its luminosity and duty cycle, so the bound alone is not a measurement of the fuelling efficiency.

For the temperature jump, use Bondi accretion as the local hot-gas estimate:Assuming the initial gas had $c_s\sim \sigma$, its density has not yet changed, and its polytropic coefficient remains comparable, heating to $c_{s,\mathrm{feed}}\sim 10\sigma$ givesThe capture radius $G{M}_{\mathrm{BH}}/c_s^2$ also shrinks by about $100$. Subsequent expansion can lower the density and suppress the mass accretion rate further. The result illustrates thermal suppression of black-hole accretion and negative active-galactic-nucleus feedback: vigorous activity can shut off its own gas supply. The factor $10^{-3}$ is conditional on the fixed-density comparison; the previous dynamical upper bound is not itself an actual pre-heating Bondi rate.

The usual Pressure-opacity zones of a Shakura--Sunyaev thin disk, in order of increasing cylindrical radius, are an inner region dominated by radiation pressure and electron-scattering opacity, a middle region dominated by gas pressure and electron-scattering opacity, and an outer region dominated by gas pressure and the Kramers' opacity law, often approximated by free-free opacity. The boundaries depend on mass, mass accretion rate, and the alpha disk viscosity parameter. At sufficiently large radius, accretion-disk self-gravity or changes in ionization can invalidate this three-zone model.

For the inner region writeHere $F$ is the standard thin-disk dissipation flux through one face, and $R_{\ast }$ is the innermost stable circular orbit treated as a zero-torque inner boundary. The given factor $1/2$ in the dissipation rate is appropriate to one face. For Keplerian rotation, it is $(9/8)\nu \Sigma {\Omega }_K^2$.

Let $H$ denote the disk scale height or vertical half-thickness, not the full thickness. In vertical hydrostatic equilibrium, the leading vertical acceleration is ${\Omega }_K^{2}z$, so a one-zone estimate givesIn the inner zone, $P_c\simeq P_{\mathrm{rad},\mathrm{c}}=a_{\mathrm{rad}}{T}_c^4/3$. For an optically thick medium with constant electron-scattering opacity, radiative diffusion gives the emergent one-face fluxEquivalently, combine $F=-c({\kappa }_{\mathrm{es}}\rho {)}^{-1}d{P}_{\mathrm{rad}}/dz$ with vertical pressure balance. Equating the radiative flux to the dissipated flux yields the radiation-supported inner-disk height:Order-unity vertical-profile factors depend on the precise height convention. The central mass cancels at fixed physical mass accretion rate, because both vertical gravity and local dissipation are proportional to $G{M}_{\mathrm{BH}}/R^3$.

If ${\dot{m}}_{\mathrm{Edd}}=L_{\mathrm{Edd},\mathrm{es}}/({\eta }_0{c}^2)$ and $r_g=G{M}_{\mathrm{BH}}/c^2$, this can also be writtenThe radiation-supported height rises from the formal zero at $R_{\ast }$ toward the approximately constant value $H_0=3{\kappa }_{\mathrm{es}}\dot{m}/(8\pi c)$. That zero is not a reliable description of the physical plunging region: the radiation-dominated approximation and the Newtonian zero-torque formula fail very close to the inner edge.

Beyond the inner zone, the height grows slowly with radius rather than staying on that plateau forever. Far from the inner boundary, the standard gas-pressure scalings give $H\propto R^{21/20}$ in the electron-scattering zone and $H\propto R^{9/8}$ in the Kramers zone. These follow from $H\propto T_c^{1/2}{R}^{3/2}$ with $T_c\propto R^{-9/10}$ and $R^{-3/4}$, respectively. A schematic profile is

vertical half-thickness H
  ^
  |                                               / outer gas/Kramers
  |                                          ____/
  |                                     ____/ middle gas/scattering
  |              ______________________/
  |          ___/ inner radiation/scattering plateau
  |       __/
  |     _/
  +----|----------------------------------------------> radius R
      R_*              zone boundaries are schematic

The height formula also tests the thin disk assumption. For a nonrotating hole with $R_{\ast }=6r_g$, the largest $H/R$ in this inner approximation occurs at $R=(9/4)R_{\ast }$. At $\dot{m}={\dot{m}}_{\mathrm{Edd}}$ and ${\eta }_0=0.1$, it is about $0.37$, so a disc this close to the conventional Eddington accretion rate is only marginally thin. Slim accretion disk effects and advected heat can then matter.

Well inside the Bondi accretion radius, pressure is negligible and the inflow approaches the free-fall speedSteady spherical mass conservation givesEstimating the outward optical depth across a radial scale $r$ with electron-scattering opacity $\kappa ={\sigma }_T/m_p$,The integral optical depth to infinity differs only by a factor of two in this free-fall approximation.

At the Schwarzschild radius, $\tau =1$ whenIf the conventional Eddington accretion rate is ${\dot{M}}_{\mathrm{Edd}}=L_{\mathrm{Edd}}/({\eta }_0{c}^2)$, then ${\dot{M}}_{\mathrm{crit}}=2{\eta }_0{\dot{M}}_{\mathrm{Edd}}$. This comparison depends on the reference efficiency ${\eta }_0$, while ${\dot{M}}_{\mathrm{crit}}$ itself follows directly from optical depth.

Writing $\dot{m}=\dot{M}/{\dot{M}}_{\mathrm{crit}}$ givesFor $\dot{m}>1$, the inner flow is optically thick and its photosphere lies at

The diffusion speed is $v_{\mathrm{diff}}\sim c/\tau$. Henceand equality defines the photon-trapping radiusInside $r_t$, inward advection outruns outward diffusion, so photon trapping in an accretion flow carries most released radiation into the black hole.

Only the binding energy released outside $r_t$ escapes efficiently. Its order of magnitude isThe radiative efficiency of black-hole accretion therefore decreases inversely with supply rate:Order-one coefficients depend on the optical-depth and inner-boundary conventions, but the luminosity saturation and $eta\propto {\dot{M}}^{-1}$ are robust.

The inner Shakura--Sunyaev thin disk is supported mainly by radiation pressure and has nearly constant electron-scattering opacity. Equating radiative diffusion with local dissipation and imposing vertical hydrostatic equilibrium gives the characteristic half-thicknessThe radiation-supported height approaches a plateau away from the inner boundary. The formal zero at that boundary is outside the reliable radiation-dominated approximation. At an Eddington accretion rate defined using efficiency $0.1$, this height can make the inner flow only marginally thin.