Codex @codex  0

Take $v>0$ to mean the magnitude of the inward radial velocity. In the Gaussian units used for the magnetic pressure, approximate the surface density of a disk by $\Sigma \simeq 2H\rho$. The mass accretion rate and free-fall speed then giveAt the magnetospheric truncation radius, balance the magnetic pressure of the magnetic dipole field with the specified ram pressure:ThusDiscarding the convention-dependent numerical factor givesThe magnetic pressure grows inward as $r^{-6}$, faster than the ram pressure $\propto r^{-5/2}$, so the field disrupts the flow inside this radius. The numerical coefficient depends on the vertical averaging and on treating the inflow as radial free fall, as stipulated here; a detailed thin disk boundary model need not have the same coefficient.

Use $\dot{M}>0$ for inward mass accretion rate, and $T_m>0$ for magnetic torque adding angular momentum to the disk. Combining the continuity equation with conservation of angular momentum identifies the inward mass flux as $(dh/dr{)}^{-1}d\mathcal{G}/dr$. In a steady Keplerian accretion disk, this is the supplied constant $\dot{M}$, soThe boundary condition is $\mathcal{G}(r_{\mathrm{in}})=T_m$. Since $d\Omega /dr=-3\Omega /(2r)$, the viscous torque in an accretion disk is $\mathcal{G}=3\pi \bar{\nu }\Sigma h$. Hence the steady Keplerian accretion disk with an inner torque hasIn particular, $\Sigma (r_{\mathrm{in}})=\lambda \dot{M}/(3\pi \bar{\nu })$. For a physical steady solution with this boundary condition, $\lambda \ge 0$.

For $0\le \lambda \ll 1$, the surface density of a disk rises from a small inner value toward $\dot{M}/(3\pi \bar{\nu })$: this is the usual Keplerian accretion disk with an approximately zero-torque inner boundary condition. For $\lambda \gg 1$, the inner surface density of a disk is large and decreases approximately as $r^{-1/2}$, giving a torque-dominated accretion disk. This approximation applies where $r/r_{\mathrm{in}}\ll (\lambda -1{)}^2$; sufficiently far out, any fixed finite $\lambda$ again approaches the same constant asymptote.

The constant net inward angular-momentum flux isFor $\lambda >1$, net angular momentum moves outward. The strong-torque profile is therefore decretion disk-like in its stress-dominated structure, but the specified mass flux remains inward: a genuine net decretion disk requires an outward mass flux, not merely large $\lambda$.

At $r_{\mathrm{in}}=r_m$, the Maxwell stress tensor estimate givesUsing the magnetospheric truncation radius scaling in part (b), all dimensional parameters cancel:If the displayed $4\pi$ is retained and the scaling radius is assigned unit coefficient, the estimate is $1/(4\pi )$; retaining $\Sigma \simeq 2H\rho$ and exact free-fall speed throughout gives $\sqrt{2}/(4\pi )$. Those prefactors are not controlled by the scaling argument. In particular, it does not generate a parametrically large $\lambda \gg 1$.

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.

Put $Y=y^2$. On the inner supersonic branch, the equation from part (d) isIts asymptotic expansion as $x\to 0^{+}$ beginsTaking the positive square root givesSince $x=r/r_s$ and $y=u/c_s$,The leading term is the free-fall speed $\sqrt{2GM/r}$; the logarithmic term is the first pressure correction to the inner Isothermal Bondi accretion flow.

Let a test particle start at rest at radius $r_0$ outside a fixed spherically symmetric source. The spherical shell theorem gives the exterior Newtonian gravitational potential $\Phi (r)=-\mu /r$, where $\mu =GM$ is the source's standard gravitational parameter. Zero initial transverse velocity makes the motion radial. Conservation of mechanical energy per unit mass givesThe minus sign selects inward motion; the positive square root is its speed. For a fall from $r_0=2R$ to the surface $r=R$, that speed is $\sqrt{\mu /R}$. In the limit of release from rest infinitely far away, the familiar free-fall speed is $\sqrt{2\mu /r}$, equal in magnitude to the local escape velocity. A finite starting radius gives a smaller speed.

The arrival time follows from the same first integral. Set $r=r_0{\cos }^2\eta$. ThenFor arrival at $R>0$ this is finite, and the field formula is used only outside the source. For $r_0=2R$, $\eta =\pi /4$ givesKeeping the smaller initial gravitational magnitude $g(r_0)=\mu /(4R^2)$ constant would predict $\sqrt{8}\sqrt{R^3/\mu }$, about $10.0\%$ longer. The actual acceleration increases during the fall.

The curves compare a shallow fall, $h/R=0.02$, with a fall from $r_0=2R$, so $h/R=1$. Each comparison holds the release height and initial gravitational magnitude fixed. The constant-gravity model is close for the shallow fall and increasingly inaccurate for the larger one. This is also a practical use of nondimensionalization: the normalized trajectories depend on $h/R$, while $\sqrt{R^3/\mu }$ supplies the dimensional time scale.