Past exam of the mathematics course of the University of Cambridge / 2007 / iii / Paper 70 (Accretion Discs)

For an equatorial circular orbit in the Paczyński-Wiita potential, radial force balance givesIts specific angular momentum is $j=R^2\Omega =\sqrt{GM}{R}^{3/2}/(R-R_S)$. Perturbing the circular radius at fixed $j$ gives the radial epicyclic frequency squared ${\kappa }^2=R^{-3}d(j^2)/dR$. ThereforeCircular motion is radially stable for $R>3R_S$, marginal at $3R_S$, and unstable between $R_S$ and $3R_S$. Thus the inner edge of the nearly circular thin accretion disk is $R_{\mathrm{in}}=3R_S$, joining a plunging region of a black-hole accretion disk inside the innermost stable circular orbit. This conclusion concerns the circular thin-disk model, not an assertion that all gas disappears inside this radius.

Take the inward mass accretion rate $\dot{M}$ as positive, and write the positive outward viscous torque in an accretion disk asSteady conservation of angular momentum makes the net inward flux $C=\dot{M}j-\mathcal{G}$ independent of radius. Impose a zero-torque inner boundary condition at $3R_S$, appropriate to an approximately stress-free transition to plunging motion. ThenThis is the steady accretion disk with arbitrary rotation law angular-momentum balance.

The local viscous dissipation integrated through the full disk thickness is $D=\nu \Sigma (R{\Omega }^{\prime }{)}^2=-\mathcal{G}{\Omega }^{\prime }/(2\pi R)$. It is the power per unit planar area summed over both radiating faces. Sincesubstitution gives the two-face dissipation of a Paczyński-Wiita disk:In particular $D$ vanishes at the stress-free inner edge and is nonnegative outside it because $j$ increases there. A one-face radiative flux is $D/2$; the total accretion luminosity is therefore $L={\int }_{3R_S}^{\infty }2\pi RDdR$, with no additional factor of two.

For a direct evaluation put $x=R/R_S$. ThenThe first integral is $7/4$, using $(3x-1)/(x-1{)}^3=3/(x-1{)}^2+2/(x-1{)}^3$. For the second,so it equals $2/[\sqrt{3}(3-1)]=1/\sqrt{3}$. HenceThis is the radiative efficiency of black-hole accretion within this pseudo-Newtonian model.

Independently, the specific orbital energy of a circular particle isThus $E(\infty )=0$ and $E(3R_S)=-GM/(8R_S)=-c^2/16$. The energy released per unit accreted mass is exactly $c^2/16$, confirming the integrated accretion luminosity. The Paczyński-Wiita circular orbit reproduces the Schwarzschild inner radius but gives an approximate binding energy; the efficiency derived here is for the stated potential.

Far from the inner edge of a steady Keplerian accretion disk, the inner angular-momentum flux is negligible compared with $\dot{M}j$. Since $\Omega \propto R^{-3/2}$, the viscous torque in an accretion disk becomes $\mathcal{G}=3\pi \nu \Sigma j$. The steady balance $\mathcal{G}\simeq \dot{M}j$ therefore gives $\nu \Sigma =\dot{M}/(3\pi )$. With the stated surface density of a disk,Steady conservation of mass gives $\dot{M}=-2\pi R\Sigma v_R$, and consequentlyBoth relations use the far-inner-edge approximation, rather than extrapolating the stress-free boundary profile to the inner edge.

For the trace passive scalar, put $\sigma =\Sigma C$. Substitute this into the contaminant conservation equation and expand its advective flux:The braces multiplying $C$ vanish by background conservation of mass. Hence the concentration satisfiesThe small contaminant fraction ensures it does not alter the ambient disk dynamics. The surface-density model already averages vertically; an axisymmetric background and ring-like injection give $C=C(R,t)$. More generally, rapid differential orbital motion winds azimuthal structure and facilitates mixing, but radial-only dependence remains a symmetry or mixing assumption, not a property of every initial injection.

For the axisymmetric advection-diffusion equation, cylindrical geometry givesHere $k=\zeta \nu$ and $k\Sigma =\zeta {\nu }_0{\Sigma }_0$ are independent of radius. ThusLet $Y=\ln (R/R_{\ast })$, with any fixed reference length $R_{\ast }$. Since $C_Y=R{C}_R$ and $C_{YY}=R^2{C}_{RR}+R{C}_R$, this is the tracer diffusion in a power-law accretion disk equationAssume $\zeta >0$, ${\nu }_0>0$, and write $v=3{\nu }_0/2$, $D=\zeta {\nu }_0$. The change $X=Y+vt$ turns this into the heat equation $C_t=D{C}_{XX}$. The radial Dirac delta function must be transformed with its Jacobian:Convolution with the Gaussian heat kernel givesThis full-line logarithmic solution extends the far-edge power-law model to $0<R<\infty$. An actual finite inner boundary would require its own boundary condition and modified kernel.

The fraction must be computed with the contaminant mass measure, rather than with $CdR$. IndeedThe total contaminant mass is the conserved quantity $M_c=2\pi {\Sigma }_0{C}_0/R_0$. Its normalized logarithmic radius is therefore a normal distribution with mean $Y_0-vt$ and variance $2Dt$. Integrating this tail gives, for $r>0$,The complementary error function accounts for the one-time mass fraction outside the specified radius.

For an outward target $r_1>1$, put $\ell =\ln r_1>0$. Since $\mathrm{erfc}$ is decreasing, maximize $F$ by minimizing $z(t)=(\ell +vt)/\sqrt{4Dt}$. Its derivative has the sign of $vt-\ell$, so the unique maximum occurs at $t_{\ast }=\ell /v=2\ln r_1/(3{\nu }_0)$. At this time $z(t_{\ast })=\sqrt{v\ell /D}$. Thus the printed expression is correctly derived as the maximum instantaneous exterior fraction:

It is not, however, the probability of ever reaching the radius. A tracer can pass through the target and later return; its history is not counted by a one-time exterior fraction. The underlying logarithmic diffusion process iswhere $B_t$ is standard Brownian motion. To compute its actual hitting probability, first put absorbing levels at $-A$ and $\ell$ relative to $Y_0$. Conditioning on the first small time step gives the backward equation $D{h}^{\prime \prime }-v{h}^{\prime }=0$ for the probability of reaching the upper level first, with $h(-A)=0$, $h(\ell )=1$. Solving givesLetting the lower level tend to $-\infty$ exhausts all finite-time upward-hitting paths. Therefore the correct eventual probability in this idealized model isFor example, if $\ln r_1=2\zeta /3$, the maximum exterior fraction is $\frac{1}{2}\mathrm{erfc}(1)\simeq 0.07865$, whereas the actual hitting probability is $e^{-1}\simeq 0.36788$. This is an explicit counterexample to the last inference in the question and an instance of maximum occupation is not a hitting probability.

For $r_1=1$, the target is the initial radius and the hitting probability is one, while ${\sup }_{t>0}F(1,t)=1/2$ with no attained maximum. For $0<r_1<1$, the inward drift ensures eventual hitting with probability one, and ${\sup }_{t>0}F(r_1,t)=1$ as $t\downarrow 0$. The square-root expression involving $\ln r_1$ is therefore only meaningful for the outward-target case, and even there it describes occupancy rather than first passage.

Reflection symmetry gives ${\Phi }_z(R,0)=0$. In a vertically stable thin disk, expand the vertical force as ${\Phi }_z(R,z)={\Omega }_z^{2}z+O(z^3)$, where ${\Omega }_z^2={\Phi }_{zz}(R,0)$ is the vertical epicyclic frequency squared. For local hydrostatic equilibrium and height-independent isothermal sound speed, $p=c_s^2\rho$ andThus $\rho ={\rho }_0\exp [-z^2/(2H^2)]$ with $H=c_s/{\Omega }_z$. Integrating the pressure and mass density through the disk gives the vertically integrated pressure of an isothermal disk:The force expansion is approximate in $H/R$, while $P=c_s^2\Sigma$ is exact for a height-independent sound speed.

In a nearly Keplerian disk, the vertical oscillation and radial epicyclic motion frequencies are close to the orbital frequency. The exact leading orbital differences are $\Omega -{\Omega }_z$ and $\Omega -\kappa$, soThus ${\omega }_1$ is the nodal precession rate of a slightly tilted orbit, and ${\omega }_2$ is the apsidal precession rate of a slightly eccentric orbit. Their signs distinguish prograde from retrograde precession. This also agrees with the sign convention in the complex tilt: a freely precessing ring obeys $W_t=i{\omega }_{1}W$.

For exactly Keplerian motion, ${\Omega }_z=\kappa =\Omega$, so ${\omega }_1={\omega }_2=0$ and $P/\Sigma =c_s^2$. In the local approximation, freeze the slowly varying coefficients and take $W,G\propto e^{i(kR-\omega t)}$. The inviscid warp equations of a nearly Keplerian disk becomeEliminating $G$ gives ${\omega }^2=k^{2}P/(4\Sigma )=k^2{c}_s^2/4$. HenceThese are bending waves of a Keplerian disk, carrying changes in orbital tilt. The local approximation needs wavelength small compared with $R$; the secular thin-disk warp equations also describe the long-vertical-wavelength regime $H\ll \lambda \ll R$.

For a nonspinning hole, $a=0$ gives ${\omega }_1=0$. A steady disk warp then satisfies $G_R=0$. The zero inner torque forces $G=0$ throughout, and the second warp equation, with positive $P$, forces $W_R=0$. The outer orientation fixes $W(R)=W_{\infty }$ everywhere. An apsidal precession term by itself therefore does not require a tilt warp.

For nonzero spin, set the time derivatives to zero. The two equations giveEliminating the internal torque yieldsUse ${\omega }_2=3\Omega /r$, ${\omega }_1=2a\Omega r^{-3/2}$, ${\Omega }^2{R}^3=GM$, and $P\simeq {\Omega }^2\Sigma H^2$. The first coefficient becomes $GM\Sigma H^{2}r/12$ and the second becomes $2aGM\Sigma r^{-3/2}$. ThereforeThis is the stationary disk warp equation under the stated weak-precession approximations. The inner condition is $W_R(R_{\mathrm{in}})=0$ when its coefficient is finite and nonzero, because $G=iGM\Sigma H^{2}r{W}_R/12$.

Let $R_g=GM/c^2$, so $R=R_{g}r$, and impose $H=ϵR$, $\Sigma \propto R^{-3/4}$. Converting both radial derivatives carefully to $r$ and dividing out the common positive factor givesor $W_{rr}+9W_r/(4r)+(24a/{ϵ}^2)r^{-9/2}W=0$. For $a>0$, introduceSince $x_r$ is proportional to $r^{-9/4}$, the first-derivative term cancels exactly, leaving $W_{xx}+W=0$. Thus $W=A\cos x+B\sin x$.

The outer boundary $r\to \infty$ is $x\to 0$, giving $A=W_{\infty }$. The specified inner-radius relation is exactly $x_{\mathrm{in}}=\pi$. The inner zero-torque condition $W_x(\pi )=0$ then gives $B=0$. Hence the stationary bending warp of a Kerr disk isThe tilt magnitude is $\beta (r)=|W_{\infty }||\cos x|$. Far out it tends to the imposed tilt. Moving inward, it decreases to zero atthen grows again with opposite signed tilt, reaching $W(r_{\mathrm{in}})=-W_{\infty }$. All tilt vectors lie on the same line in the complex $W$ plane. Their azimuth is $\arg W_{\infty }$ outside the node and differs by $\pi$ inside; azimuth is undefined at the aligned node itself. The tilt vector remains smooth there. Thus this is an oscillatory standing bending warp, with equal outer and inner tilt magnitudes and reversed inner orientation, rather than monotone inner alignment.

The prescribed real inner radius assumes $a>0$: for $a<0$, its square root is not real, so that particular boundary prescription does not apply. With a separately chosen positive inner radius and negative spin, put $x=4\sqrt{24|a|}{r}^{-5/4}/(5ϵ)$. The equation becomes $W_{xx}-W=0$, and the same outer tilt and zero inner torque giveThis negative-spin solution decreases smoothly toward the inner disk without a sign reversal. The specified positive-spin solution is the cosine profile above; the $a=0$ case is the flat solution already derived.