Codex @codex  0

Assume a thin disk has height-independent opacity ${\kappa }_{\mathrm{op}}$ and radiation-to-gas pressure ratio $\beta$, with vertical gravity $-{\Omega }_z^{2}z$. Hydrostatic equilibrium and radiative diffusion implyBalancing its derivative against viscous heating $\rho \nu q^2{\Omega }^2$ fixes the dynamic viscosityIt is independent of height under these assumptions. For a spherically symmetric potential, ${\Omega }_z=\Omega$. Integrating $\rho \nu$ through the full surface-to-surface thickness $2H$ gives $\bar{\nu }\Sigma =2H\rho \nu$, which determines the thickness from the radial accretion rate solution.

Let $M_{\ast }$ be the central mass, so ${\Omega }^2=G{M}_{\ast }/r^3$ and the specific orbital energy is $\epsilon =-G{M}_{\ast }/(2r)$. Define the inward accretion rate $\dot{M}>0$. The rate at which the gas acquires binding energy in moving from $b$ to $a$ isThe vertically integrated viscous heating per unit horizontal area, including the entire column, is $\bar{\nu }\Sigma (r{\Omega }^{\prime }{)}^2$. Since $r{\Omega }^{\prime }=-3\Omega /2$, the dissipated power in the annulus isHere the column integral already includes dissipation associated with both disk faces; adding another factor of two would double-count it. Performing the elementary integrals givesFor the whole disk, take $a=r_{\mathrm{in}}$ and $b\to \infty$. Both expressions giveIf $a\gg r_{\mathrm{in}}$, the bracket $1-\sqrt{r_{\mathrm{in}}/r}$ is uniformly close to one throughout the annulus. Comparing the positive integrals showswith the correction negative. Thus far from the inner edge the viscous dissipation is approximately three times the local binding-energy acquisition.

There is no violation of conservation of energy because the viscous torque in an accretion disk also transports mechanical energy. The outward viscous torque in an accretion disk isand its outward energy flux is $\Omega \mathcal{G}=G{M}_{\ast }\dot{M}[r^{-1}-\sqrt{r_{\mathrm{in}}}{r}^{-3/2}]$. The explicit expressions satisfyEnergy transported from the inner disk supplies the extra dissipation at larger radii. The torque power vanishes at the zero-torque inner boundary condition and tends to zero at infinity, leaving global balance. This is torque transport of energy in a steady accretion disk.

Both profiles have integrable central cusps. Directly integrating the surface density of a disk and its specific angular momentum $j=\sqrt{GMr}$ givesFor $a=1$, substitution of the similarity functions givesThe inward accretion rate is $\dot{M}(r)=-2\pi r\Sigma {\bar{u}}_r$. At the origin it tends to $M_d/(5t)=-d{M}_d/dt$. The outward viscous torque in an accretion disk isFor this profile $\mathcal{G}\to 0$ as $r\to 0$, and the advected angular-momentum flux $\dot{M}j$ also vanishes there. Both outer fluxes vanish at $R$. Thus disk mass decreases through central accretion while disk angular momentum is conserved. Including the accreted central mass restores total mass conservation for the disk-plus-sink system.

For $a=5/4$, instead,Here the mass flux tends to zero at the origin, despite positive radial velocity everywhere else. The inner torque has the nonzero limitThus the nonaccreting similarity disk supplied by an inner torque conserves disk mass but gains angular momentum from an external central torque. Its vanishing edge flux cannot supply that gain. The disk is not globally isolated in angular momentum; the combined disk and torque source may of course conserve it.

These profiles illustrate different boundary-controlled similarity classes of the Keplerian viscous diffusion equation. The $a=1$ solution is compatible with a zero-torque inner boundary condition and conserved finite angular momentum, and can describe a late-time spreading accretion disk after the detailed initial profile has been forgotten. Its mass diverges as $t\downarrow 0$, so the formal zero-radius origin is not a finite-mass initial disk. The $a=5/4$ solution has finite initial mass concentrated at the origin and initially zero disk angular momentum, but requires the specified inner torque for $t>0$; that torque's $t^{-3/4}$ singularity is time-integrable. It may describe a spreading nonaccreting disk for compatible torque-driven boundary conditions, but cannot be an isolated angular-momentum-conserving disk's generic asymptotic state. Similarity scaling suggests their possible asymptotic roles; it does not establish convergence for arbitrary initial and boundary data.

Take $\dot{M}>0$ to mean inward accretion rate, and let $S$ count loss through both faces per unit projected area of the accretion disk. An annulus of width $dR$ loses $2\pi RSdR$ to the wind. In steady state, its incoming mass flux at its outer edge equals the inward flux at its inner edge plus this loss. Hence mass conservation givesIf a loss rate is instead defined per individual face, it must be doubled before substituting for $S$ here.

Use Keplerian rotation ${\Omega }^2=GM/R^3$. The vertically integrated viscous dissipation, summed over both disk faces, isHere $\nu$ is the kinematic viscosity, and $\Sigma$ is the surface density of a disk; for height-dependent viscosity the product means $\int \rho \nu dz$. Material already in a circular orbit has specific orbital energy $v_K^2/2-GM/R=-GM/(2R)$. Giving it zero total binding energy therefore costs an additional $GM/(2R)$ per unit mass. The total launch speed is the escape velocity $\sqrt{2GM/R}$, but one must not charge the wind for the orbital kinetic energy it already possesses. Equating wind power to $f{Q}^{+}$ givesCombining this with mass conservation yields the wind-powered mass loss from a steady accretion disk relation

Define the outward viscous torque in an accretion disk asThe net inward angular momentum flux is $\dot{M}h-\mathcal{G}$. Because each unit mass lost to the wind carries the local specific angular momentum $h$, annular conservation of angular momentum readsIn the wind region, $\mathcal{G}=Rh{\dot{M}}^{\prime }/(3f)$. For constant $f$, differentiation and cancellation of $h$ givesThus the requested Cauchy-Euler differential equation is

For $f=1/3$, a power-law trial $\dot{M}\propto R^p$ gives $p^2+p/2-1/2=0$, with roots $p=1/2,-1$. Therefore $\dot{M}=A\sqrt{R}+B/R$ in $R_{\ast }\le R\le R_0$. A zero inner viscous torque in an accretion disk implies ${\dot{M}}^{\prime }(R_{\ast })=0$, so $B=A{R}_{\ast }^{3/2}/2$. There is no point mass sink at $R_0$, so the accretion rate is continuous there and equals ${\dot{M}}_0$. The complete result isIts derivative is nonnegative throughout the wind region, as required for a nonnegative loss rate $S$. The derivative can jump to zero outside $R_0$: the energy relation applies only inside the wind zone and does not impose ${\dot{M}}^{\prime }(R_0)=0$ on its inner limit. The viscous torque in an accretion disk remains continuous and extends outside via ${\mathcal{G}}^{\prime }={\dot{M}}_0{h}^{\prime }$.

At the inner boundary the exact accreted fraction isThe last expression is a leading-order asymptotic equality, not the exact finite-radius fraction. Most of the supplied mass is expelled before reaching the star when $R_{\ast }/R_0$ is small.

Take the accretion rate $\dot{M}>0$ to mean inward flow. In a steady state, conservation of mass makes the inward mass flux independent of radius, so $\mathcal{F}=\dot{M}$. Integrating its radial expression givesFor Keplerian rotation, $\Omega =(G{M}_{\ast }/r^3{)}^{1/2}$ and the viscous torque in an accretion disk is $\mathcal{G}=3\pi \bar{\nu }\Sigma \sqrt{G{M}_{\ast }r}$. The zero-torque inner boundary condition therefore sets $\bar{\nu }\Sigma =0$ at $r_{\ast }$. It fixes $C=-\dot{M}\sqrt{r_{\ast }}/(3\pi )$, giving the steady density profileThis is the Keplerian accretion disk solution on $r_{\ast }<r<r_{\mathrm{out}}$. It determines the kinematic viscosity–surface density product; a separate viscosity closure is needed to turn it into an explicit power law. If the kinematic viscosity is finite and nonzero at the inner edge, the surface density tends to zero there.

The standard thin-disk dissipation flux, summed over both faces, follows by substituting the Keplerian accretion disk profile into the viscous heating rate:The Stefan–Boltzmann law gives $\mathcal{C}=2\sigma T^4$ because there are two emitting faces. Local radiative equilibrium consequently gives the effective-temperature profile of a zero-torque diskHere $\sigma$ is the Stefan-Boltzmann constant. The scale $T_{\mathrm{in}}$ is not the temperature exactly at the inner boundary: the zero-torque inner boundary condition makes that formal temperature zero. Maximizing $T^4$ shows that the maximum effective temperature of a zero-torque disk occurs at $r=49r_{\ast }/36$, with $T_{\max }\simeq 0.488T_{\mathrm{in}}$.

At equal central mass and accretion rate, characteristic effective temperatures scale as $r_{\ast }^{-3/4}$. Thus, comparing corresponding values of $r/r_{\ast }$,The neutron-star disk is about 180 times hotter. The Planck law and Wien displacement law move its characteristic emission to about 180 times higher frequency, or 180 times shorter wavelength. A white dwarf disk commonly emits in optical and ultraviolet bands, while the hotter neutron star disk can emit in X-rays. Absolute bands require an actual accretion rate; the relative shift follows directly from the stated scaling.

A cataclysmic variable is a close semidetached binary in which a white dwarf accretes from a cool, usually low-mass donor star filling its Roche lobe. In an ordinary hydrogen-rich system the donor star is often near the lower main sequence. In a nonmagnetic cataclysmic variable (a disk-fed CV), the white dwarf's magnetic field is too weak to control the flow over the disk. This leaves a characteristic disk-fed geometry, summarized in the original schematic below. NASA's introduction to cataclysmic variables describes the basic components.

Roche-lobe overflow passes through the inner Lagrange point $L_1$. A nearly ballistic stream bends in the rotating frame and strikes the outer accretion disk, producing an accretion-disk stream-impact hot spot. Its retained angular momentum prevents direct radial infall. Viscous evolution of an accretion disk transports angular momentum outward while gas moves inward through a nearly Keplerian accretion disk. Close to the white dwarf, gas slows from orbital rotation toward the stellar rotation in an accretion-disk boundary layer. The disk, accretion-disk stream-impact hot spot, accretion-disk boundary layer, white dwarf and donor star all contribute to the spectrum and, where the orbital inclination permits, to the eclipses of an eclipsing binary. In particular, “nonmagnetic” describes the accretor's control of the flow; it does not imply that the donor star cannot sustain a magnetic field.

The gravitational power available from accretion is approximately $L_{\mathrm{acc}}=G{M}_{\mathrm{WD}}{\dot{M}}_{\mathrm{acc}}/R_{\mathrm{WD}}$. For a slowly rotating white dwarf and a thin steady Keplerian accretion disk, the specific energy changes from approximately zero far out to $-G{M}_{\mathrm{WD}}/(2R_{\mathrm{WD}})$ at the inner disk. Roughly half the available power is radiated by the accretion disk; the rest is released in the accretion-disk boundary layer as the orbital kinetic energy is dissipated. Stellar rotation and departures from a steady thin disk change this partition.

A classical nova has a different energy source. Transferred hydrogen accumulates on the white dwarf; compression heats the base of its envelope until temperature-sensitive hydrogen burning accelerates. Electron degeneracy pressure initially weakens the expansion response to heating, helping a thermonuclear runaway develop. The envelope subsequently expands and ejects material, producing a large optical outburst followed by a decline as the ejecta expand and residual burning ends. The white dwarf usually survives, so continued accretion can build another fuel layer. A rough recurrence scale is the ignition-envelope mass divided by the mean accretion rate; both this mass and the rate vary strongly among systems. The event is an envelope eruption, and the retained fraction is not automatically unity. Starrfield, Iliadis and Hix's nova calculations explains this nuclear mechanism. A classical nova is powered by unstable nuclear burning on the white dwarf.

A dwarf nova undergoes recurrent, shorter brightenings powered principally by enhanced gravitational accretion. The hydrogen-ionization disk instability creates cold, mostly neutral and hot, ionized branches of the disk's accretion-disk thermal S-curve, separated by unstable equilibria. In quiescence the cool disk stores matter because inward transport is slow. Once a critical surface density is reached, a heating transition puts the disk into a hotter, more state with higher effective viscosity: the inward accretion rate and luminosity rise and the disk drains. A cooling transition returns it to quiescence, completing the cycle. A persistent increase in the donor star's transfer rate is not required. Sufficiently high transfer rates can keep the disk on its hot stable branch, giving a nova-like variable rather than ordinary disk cycles. Lasota's disk-instability analysis and Hameury's disk-instability review develop this picture. A dwarf-nova outburst is a disk instability, not a white-dwarf thermonuclear explosion. The names classify mechanisms and need not identify permanently distinct binaries: a nova-producing binary can also possess an unstable disk between nuclear eruptions.

The cataclysmic-variable orbital-period distribution is not smooth. For ordinary hydrogen-rich cataclysmic variables, prominent features are a cataclysmic-variable period gap around two to three hours, a cataclysmic-variable period minimum near eighty minutes, and an accumulation near that minimum. These are population features rather than absolute exclusions. Selection effects matter: luminous high-accretion rate systems are easier to find than faint evolved systems. Helium-transferring binaries have a different period range and are not described by the hydrogen-rich minimum. Gänsicke and collaborators' period-minimum study documents the observed accumulation.

The Roche-lobe-filling period-density relation makes the orbital period a measure of donor structure. Combining $R_d\simeq 0.46a(M_d/M{)}^{1/3}$ with Kepler's third law givesAs the donor star loses mass, its stellar radius response exponent $\zeta =d\log R_d/d\log M_d$ impliesA donor with $\zeta >1/3$ evolves toward shorter orbital periods. When its effective response falls below $1/3$, continued mass loss instead lengthens the orbital period: this is the cataclysmic-variable period bounce. A very low-mass donor may be substellar and increasingly affected by electron degeneracy pressure; the ideal degenerate scaling $R_d\propto M_d^{-1/3}$ illustrates the reversal. The precise cataclysmic-variable period minimum depends on thermal disequilibrium and the strength of orbital angular momentum loss. Near a turning point $|\dot{P}|$ is small; for an approximately steady evolutionary flow of systems, the number per period interval scales as $N(P)\propto 1/|\dot{P}|$, explaining the accumulation. Knigge, Baraffe and Patterson's donor-based evolutionary study relates the donor star's response to these period features.

Long-term transfer is driven mainly by losses of orbital angular momentum, rather than by disk outbursts. Magnetic braking of a binary star removes the cool donor star's spin through a magnetized stellar wind. Tidal synchronization makes the orbit replenish that spin, so the wind extracts orbital angular momentum. This is usually the dominant standard driving mechanism above the cataclysmic-variable period gap. Gravitational-wave emission from a binary system supplies a baseline loss, especially important below the gap. For a weak-field, slowly moving circular orbit, the circular gravitational-wave inspiral givesBecause ordinary transfer from the lighter donor star tends to expand its Roche lobe if orbital angular momentum is conserved, an external loss is needed to sustain contact. In the conservative contact approximation, the binary mass-transfer contact equation readswhere the stellar radius response exponent must match the evolutionary timescale and additional donor expansion has been neglected. The positive denominator on the stable branch makes ${\dot{J}}_{\mathrm{ext}}<0$ drive ${\dot{M}}_d<0$. Nova ejecta can introduce additional nonconservative binary mass transfer. Knigge's evolutionary discussion describes the standard loss mechanisms and their limitations.

In the disrupted magnetic braking model, relatively rapid mass loss above the cataclysmic-variable period gap keeps the donor star inflated relative to stellar thermal equilibrium. When the donor approaches the fully convective star transition, the model postulates a substantial reduction in magnetic braking of a binary star. The donor star can contract within its Roche lobe, suppressing Roche-lobe overflow near the upper edge of the gap. Gravitational-wave emission from a binary system continues to shrink the detached binary; near the lower edge the Roche lobe again reaches the donor radius and transfer resumes. No mass transfer is needed during the detached crossing. The density relation predicts a radius ratio $(3/2{)}^{2/3}\simeq 1.31$ between contact at three and two hours if the masses remain nearly fixed, illustrating the required inflation before detachment. The torque reduction is a model ingredient, not a claim that all fully convective stars lose their magnetic fields. Zorotovic and collaborators' detached-binary study tests the predicted detached population in the gap.

The standard formation channel starts with an initially wider binary star containing two main sequence stars. The initially more massive component evolves first and becomes a giant. Unstable Roche-lobe overflow can engulf the companion in a common envelope. Drag causes inward orbital motion, releasing two-body orbital energy and transferring angular momentum to the envelope. If the envelope is expelled before merger, a close detached binary survives, containing the exposed core, which becomes a white dwarf, and the lower-mass companion. This explains how a binary becomes much tighter than the giant progenitor's radius would have allowed. The outcome depends on envelope binding and the efficiency of energy deposition, summarized approximately by the common-envelope energy formalism; ejection is not guaranteed. Ivanova and collaborators' common-envelope analysis discusses the relevant physics and uncertainties.

Subsequent magnetic braking of a binary star and gravitational-wave emission from a binary system shrink the detached binary until its companion fills its Roche lobe. Stable Roche-lobe overflow, for a suitable binary mass ratio and stellar radius response exponent, then creates a cataclysmic variable. Its later nuclear eruptions, disk cycles and secular orbital period evolution occur on different timescales. The formation sequence is a wide binary, envelope ejection, a close detached white-dwarf binary, and angular-momentum-driven contact.

Interpret the outer Dirac delta function as a boundary source carrying its full mass, or as the limit of a narrow source lying just inside $r_2$. Then ${\int }_u^{r_2}vS(v)dv=r_2{S}_0$ for $r_1<u<r_2$, andThus the interior is the usual steady Keplerian accretion disk with a constant inward accretion rate. The product $\bar{\nu }\Sigma$ vanishes at the stress-free inner edge, and the surface density of a disk does too if the viscosity has a finite nonzero boundary value; far from that edge $\bar{\nu }\Sigma$ approaches $\dot{M}/(3\pi )$. Its density need not itself be constant because the viscosity may vary.

The outer-edge source of a viscous disk requires care about the boundary trace. The flux just inside the injecting edge is $\mathcal{F}(r_2^{-})=\dot{M}$, while the no-flux value outside the source is $\mathcal{F}(r_2^{+})=0$; the source produces the intervening jump. A simultaneously continuous interior condition $\mathcal{F}(r_2^{-})=0$ would contradict nonzero steady injection. A symmetrically regularized delta integrated only up to its center can instead give half its full weight; that would mean a different supplied mass, not a different interior solution for the same $\dot{M}$. The physically normalized one-sided limiting source removes this convention ambiguity.

Take $0<r_1<r_2$ and positive mean kinematic viscosity $\bar{\nu }$. The steady state of the supplied viscous evolution of an accretion disk equation satisfies ${\mathcal{F}}^{\prime }=-2\pi rS(r)$. The outer no-flux condition therefore givesPositive $S$ produces positive inward accretion rate in this sign convention. Set $Y=r^{1/2}\bar{\nu }\Sigma$. In a Keplerian disk, the viscous torque in an accretion disk is proportional to $Y$, so the zero-torque inner boundary condition gives $Y(r_1)=0$. Since $Y^{\prime }=\mathcal{F}/(6\pi r^{1/2})$, a second integration yieldsThis is the steady viscous disk supplied by a radial source. It does not require $\bar{\nu }$ to be independent of radius: the equation directly determines $\bar{\nu }\Sigma$. If viscosity depends on the local thermodynamic state, the displayed formula is an implicit relation completed by a separate closure. For nonnegative integrable $S$ it gives nonnegative surface density of a disk, and the inner accretion rate equals the total supplied mass per unit time, $2\pi {\int }_{r_1}^{r_2}rS(r)dr$. Differentiating both nested integrals recovers the flux and the differential equation.

Use a gas-pressure-supported, optically thick alpha disk, with $\mathcal{R}=k/({\mu }_m{m}_p)$ and positive parameters. Here $\sigma$ is the Stefan–Boltzmann constant in the Stefan–Boltzmann law, and $F_{+}$ and $F_{-}$ below are heating and cooling flux estimates per disk face. Hydrostatic equilibrium and the ideal gas relation give $\mathcal{R}T\sim {\Omega }^2{H}^2$, while $\rho \sim \Sigma /H$. Integrating viscous heating over one semi-thickness and estimating radiative diffusion giveHere the stipulated negative hydrogen ion opacity exponent is $\kappa ={\kappa }_0{\rho }^{1/3}{T}^{10}$, which is a local approximation distinct from other empirical opacity fits. Substituting $T\sim {\Omega }^2{H}^2/\mathcal{R}$ and equating heating and cooling givesDropping numerical constants of order one is essential here. This is the negative-hydrogen-opacity alpha-disk thickness scaling. Its angular thickness estimate isFor a central mass $M$, Keplerian rotation gives $\Omega =(GM/r^3{)}^{1/2}$, so $H/r\propto r^{53/82}\Sigma (r{)}^{-7/41}$. The isothermal sound speed $c_s^2=\mathcal{R}T$ isAn adiabatic sound speed differs by the constant factor $\sqrt{\gamma }$ if the adiabatic index is fixed. For $\Sigma \propto r^{-p}$, constant $\alpha$, and fixed composition, the disk aspect ratio grows as $r^{(53+14p)/82}$, so there is disk flaring when $p>-53/14$; $c_s\propto r^{(6+7p)/41}$. The problem supplies no numerical parameters or specified $\Sigma (r)$, so no unique numerical $H/r$ or unconditional radial power can be inferred.

If the additional standard closure of a constant inward accretion rate well outside a zero-torque edge is intended, $\bar{\nu }\sim \alpha \Omega H^2$ and $\bar{\nu }\Sigma \simeq \dot{M}/(3\pi )$ imply $\Sigma \propto (\alpha \Omega H^2{)}^{-1}$. Eliminating it gives $H^9\propto {\Omega }^{-38/3}$ at fixed $\dot{M},\alpha$, hence $H/r\propto r^{10/9}$ and $c_s\propto r^{11/18}$. These stronger powers use that extra steady-flow assumption. All these local scalings cease to apply if the thin disk or the stated opacity regime fails.

In a steady state, mass conservation makes $\mathcal{F}$ independent of radius. Taking the accretion rate $\dot{M}>0$ for inward flow gives $\mathcal{F}=-\dot{M}$. The conservation of angular momentum equation then makes $-\dot{M}h+\mathcal{G}$ constant. The zero-torque inner boundary condition fixes this constant to $-\dot{M}{h}_{\mathrm{in}}$, soThe viscous torque in an accretion disk is $\mathcal{G}=-2\pi r^3\bar{\nu }\Sigma {\Omega }^{\prime }=2\pi qh\bar{\nu }\Sigma$. Combining these expressions gives the steady accretion disk with arbitrary rotation lawFor $x=r/R_g$, the Paczyński-Wiita circular orbit formulas giveConsequentlyIn particular, $f(6)=0$ and $f\to 1$ far from the black hole, recovering the outer Keplerian accretion disk.

Inside the innermost stable circular orbit, the gas enters the plunging region of a black-hole accretion disk. Its rapid inward motion leaves little time for stresses to exchange angular momentum, motivating the zero-torque inner boundary condition. This is a thin-disc approximation; a strong magnetic stress could change it.

Matter supplied from very large radius has negligible specific orbital energy, while matter crossing the inner edge carries ${\epsilon }_{\mathrm{in}}=-\eta c^2$. With zero inner torque, no energy is supplied by a stress at that edge. If the heat released outside it escapes by radiative transfer, rather than being lost through inward advection, the integrated conservation of energy balance gives

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 standard Shakura--Sunyaev thin disk assumes: a steady state; axial symmetry; a geometrically thin disk $H/R\ll 1$; Newtonian, nearly Keplerian circular motion outside the inner edge; subsonic radial drift $|u_R|\ll u_{\varphi }$; vertical hydrostatic equilibrium; negligible disk self-gravity; an optically thick, locally thermal spectrum; local radiative balance between viscous heating and cooling; and a local alpha-viscosity stress with a zero-torque inner boundary near the innermost stable circular orbit. These assumptions also exclude dynamically dominant winds and large-scale external torques from the standard solution.

Start from the continuity equationIntegrate vertically, use axial symmetry and steadiness, and define the surface density of a disk $\Sigma =\int \rho dz$. Thenso the inward-positive accretion rate is

With $\mathcal{G}=0$ and $\mathcal{T}=2\pi v_m{h}^{\prime }\Sigma$, the evolution equation reduces to the radial continuity equationor, for the mass per unit radius $M_r=2\pi r\Sigma$,The outward mass flux is $F_M=2\pi r{v}_m\Sigma$. If $\dot{M}>0$ denotes the inward accretion rate, steady mass conservation requires $F_M=-\dot{M}$. HenceThis is positive because the magnetic drift velocity satisfies $v_m<0$.

Write $y=\nu \Sigma$. A steady inward accretion rate has $F_M=-\dot{M}$, so the integrated angular-momentum equation isFor $\beta =0$, $\mathcal{T}=0$ and thereforeafter imposing the zero-torque inner boundary condition at $r_{\mathrm{in}}$. In a Keplerian accretion disk, $h\propto r^{1/2}$ and $\mathcal{G}=3\pi hy$, so