Codex @codex  0

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.

Write $h=r^2\dot{\theta }$ for the conserved specific angular momentum. Zero specific orbital energy in the inverse-square Newtonian gravitational field givesDefine $r_0=h^2/(2GM)$, the periapsis radius for a nonradial parabolic trajectory. ThenOn the outgoing branch, integrate from the periapsis time $t_0$ and put $s=\sqrt{r-r_0}$:The incoming branch has the opposite sign. Squaring gives the radius-time relation for a parabolic Kepler orbit, a version of the Barker equation:Here $r_0$ and $t_0$ are fixed by angular momentum and the time origin at periapsis. The relation applies on both branches, with $r\ge r_0$; it is not necessary to solve for the polar angle.

For an attractive central force, the transverse acceleration vanishes. In polar coordinates,Thus the specific angular momentum $h=r^2\dot{\theta }$ is constant. Assume $h\ne 0$, so angle can parameterize the orbit, and put $u=1/r$. Primes denote angle derivatives. The chain rule givesThe outward radial component of the attractive force is $-f(u)$. Its Newton's second law equation is therefore $m(\ddot{r}-r{\dot{\theta }}^2)=-f(u)$, giving the Binet equationThe sign reflects $f$ being the positive inward force magnitude. Zero angular momentum produces radial motion and is outside this angular parametrization.

Take constant kinematic viscosity $\nu$ and the inner specific angular momentum found above. The Angular velocity of a circular Schwarzschild geodesic gives $-{\Omega }^{\prime }=3c\sqrt{m}/(2r^{5/2})$. Hence the inward radial component of four-velocity in the stated viscous model isTo convert proper time to Schwarzschild time, neglect ${\dot{r}}^2$ in the normalization as prescribed and insert the circular value of $h$. One obtainsThe change in observed Schwarzschild time equals this coordinate-time change when differential light travel time is neglected. Since $dr/d\tau =-u$,This time is nonnegative: $r^2-12m(r-3m)=(r-6m{)}^2\ge 0$ on the integration interval.

Put $x=r/(3m)$ and $X=r_1/(3m)$. The integral becomesPolynomial division gives $x^2/(x-1)=x+1+1/(x-1)$. Separately, writing $y=x-1$ gives $\int 2x/\sqrt{x-1}dx=(4/3)(x-1{)}^{3/2}+4\sqrt{x-1}$. Therefore the antiderivative inside the integral isIts value at $x=2$ is $-4/3$. Subtracting this lower endpoint yields the formal viscous inspiral time in Schwarzschild spacetimeThe logarithm has a positive coefficient. The minus sign in the printed final expression is inconsistent with its preceding integral. Differentiating the corrected bracket gives $6X[X-2\sqrt{X-1}]/(X-1)$, precisely the positive integrand after rescaling. The bracket vanishes at $X=2$; replacing $+6\log (X-1)$ by its negative would give derivative $-12$ there and a negative elapsed time for $X$ just above $2$. This is a sign error, not a choice of potential or time convention.

The result assumes constant kinematic viscosity; if $\nu$ varies with radius it remains inside the time integral. There is also a physical limitation to the stipulated approximation: $h-h_0$ vanishes quadratically at $6m$, so its formal $u$ diverges and eventually violates the neglected-radial-motion condition. The finite integral is the formal extrapolation of the nearly circular viscous model. A real accretion disk must instead match to its plunging region of a black-hole accretion disk; this approximation cannot describe that transition arbitrarily closely.

Let $\mathcal{G}$ be the outward viscous torque in an accretion disk. For a declining angular velocity, the vertically integrated viscous shear stress givesIn a steady accretion disk without radiated angular momentum, the net inward angular-momentum flux $Fh-\mathcal{G}$ is constant and equals the swallowed flux $F{h}_0$. Thus conservation of angular momentum gives $F(h-h_0)=\mathcal{G}$. Substituting $F=2\pi r\Sigma u$ cancels the surface density and provesThis uses the stipulated phenomenological viscous shear law; it does not replace it by a different fully relativistic stress prescription.

For the orbital calculation let a dot denote a derivative with respect to proper time. Substituting the geodesic energy and specific angular momentum first integrals into the timelike normalization in Schwarzschild spacetime givesDifferentiate with $h$ and $ϵ$ constant along the geodesic. Where $\dot{r}\ne 0$, division by $2\dot{r}$ givesThis equation also holds at turning points and on circular orbits. To justify that without dividing by zero, put $A=1-2m/r$ in the equatorial geodesic Lagrangian $L=(A{c}^2{\dot{t}}^2-A^{-1}{\dot{r}}^2-r^2{\dot{\varphi }}^2)/2$. Its radial Euler-Lagrange equation is $\ddot{r}=(A^{\prime }/2A){\dot{r}}^2-A{A}^{\prime }{c}^2{\dot{t}}^2/2+Ar{\dot{\varphi }}^2$. Substitution of the timelike normalization cancels the ${\dot{r}}^2$ terms and yields $\ddot{r}=-A^{\prime }{c}^2/2+(A/r-A^{\prime }/2)h^2/r^2$, exactly the displayed radial equation, including when $\dot{r}=0$.

A circular orbit requires $\ddot{r}=\dot{r}=0$. Thereforechoosing the positive sense of rotation. Differentiating the square is simpler than differentiating $h$:It is negative for $3m<r<6m$ and positive for $r>6m$, so the global minimum over timelike circular orbits isAs a stability check, the second radial derivative of $W=(1-2m/r)(c^2+h^2/r^2)$ at a circular orbit, holding $h$ fixed for the perturbation, is $W^{\prime \prime }=2m{c}^2(r-6m)/[r^3(r-3m)]$. Thus the exterior branch is stable, and $6m$ is the innermost stable circular orbit. It is the natural matching radius for the idealized zero-torque inner boundary condition.

Assume axisymmetry in a thin disk orbiting a fixed central mass $M$. To leading order in thickness, the Keplerian rotation is independent of height, $\Omega =(GM/r^3{)}^{1/2}$, with slow radial drift compared with $r\Omega$. Neglect radial pressure support in the orbital balance, changes in the central potential, mass loss through the disk surfaces and applied surface torques. The Cauchy stress tensor can represent viscous or effective turbulent/magnetic transport; only its vertically integrated azimuthal component enters the conservation of angular momentum equation.

Define the surface density and its mass-weighted radial velocity byand set $W_{r\varphi }=\int T_{r\varphi }dz$. Integrating mass conservation vertically givesThe azimuthal equation in cylindrical coordinates, multiplied by $r$ and combined with mass conservation, is the conservative equation for specific angular momentum $l=r^2\Omega$. The cylindrical stress tensor divergence contributes $r^{-1}{\partial }_r(r^2{T}_{r\varphi })+{\partial }_z(r{T}_{\varphi z})$. With vanishing surface torque, vertical integration therefore givesBecause $l(r)$ is fixed in time, subtracting $l$ times the mass equation yieldsDefine the effective viscosity defined by disk stress byThis sign convention treats $T_{r\varphi }$ as the shear stress appearing with $+\nabla \cdot \mathbf{T}$. For the Newtonian fluid stress tensor it reduces to the density-weighted viscosity of a disk, $\bar{\nu }={\Sigma }^{-1}\int \rho \nu dz$, rather than an unweighted height average.

Using $l=\sqrt{GMr}$ and $W_{r\varphi }=-(3/2)\bar{\nu }\Sigma \Omega$ in the conservation of angular momentum equation givesSubstitution into mass conservation proves the Keplerian viscous diffusion equation:The outward viscous torque in an accretion disk is $\mathcal{G}=-2\pi r^2{W}_{r\varphi }=3\pi \sqrt{GM}\sqrt{r}\bar{\nu }\Sigma$.

For $\bar{\nu }=Ar$ with $A>0$, make the heat-equation transform for linear-radius disk viscosityWriting $y=\sqrt{r}$, the diffusion equation becomes $g_t=(3A/4)g_{yy}$; the displayed rescaling of $y$ gives exactly $\boxed{g_t=g_{xx}}$, with the original time unchanged. The torque is $3\pi \sqrt{GM}g$, so the zero-torque inner boundary condition at the origin becomes $g(0,t)=0$ on the half-line $x>0$.

A positive similarity solution obeying this boundary condition isFor example, substitution of $g=t^{-1}F(x/\sqrt{t})$ reduces the heat equation to $F^{\prime \prime }+(\eta /2)F^{\prime }+F=0$, and $F=C\eta e^{-{\eta }^2/4}$ satisfies it and vanishes at zero. A replacement of $t$ by $t+t_0$ gives a finite-mass initial profile at time zero if desired.

To determine the global integrals, transform the mass and angular momentum:The Gaussian integrals giveThus the accreting similarity disk with viscosity proportional to radius loses mass as $t^{-1/2}$ while conserving angular momentum; its characteristic radius grows proportional to $At$. The accretion rate at the origin is $-{\dot{M}}_d=2\pi C\sqrt{3/A}{t}^{-3/2}$. Its radial velocity ${\bar{u}}_r=r/t-3A/2$ shows inward flow at small radius and outward spreading at large radius. The integrable inner density singularity is part of the idealized extension to $r=0$.

For the other profile, direct differentiation givesIt also satisfies the heat equation, but has $g(0,t)=t^{-1/2}\ne 0$ and $g_x(0,t)=0$. It therefore has a nonzero inner torque and zero mass flow at the origin, not the previous zero-torque inner boundary condition. HereAn arbitrary positive amplitude multiplies these quantities. This constant-mass disk with viscosity proportional to radius spreads outward under a central torque that supplies angular momentum, with no ongoing central mass supply. It has the reflecting Neumann boundary condition for the heat equation and can be regarded as the spreading of an initially central mass distribution driven by an applied torque. It is distinct from the accreting zero-torque solution; the ideal origin does not model a finite stellar surface or a relativistic inner disk.

Write the gravitational potential with spherical symmetry as $\Phi (r)$ and let $\mathbf{n}$ be the unit vector normal to a circular orbit, oriented in the direction of its angular velocity. Put $\mathbf{\Omega }=\Omega \mathbf{n}$ and denote the specific angular momentum by $\mathbf{h}=h\mathbf{n}$. Circular force balance and the specific orbital energy giveDifferentiating along the family of circular orbits,A change in orbital orientation contributes $hd\mathbf{n}$ to $d\mathbf{h}$, but $\mathbf{n}\cdot d\mathbf{n}=0$. ConsequentlyThis circular-orbit energy gradient accounts for radial and inclination changes together; spherical symmetry makes the orbital energy independent of the plane.

For either variable-mass particle, $\mathbf{H}=M\mathbf{h}$ and $E=M\epsilon$. Thus $d\mathbf{H}=Md\mathbf{h}+\mathbf{h}dM$, soConservation of mass and conservation of angular momentum impose $d{M}_2=-d{M}_1$, $d{\mathbf{H}}_2=-d{\mathbf{H}}_1$. Adding the two particle contributions provesTo interpret the signs, label the inner orbit 1 and the outer orbit 2. For $f(r)=\epsilon -\Omega h$ the first identity gives $f^{\prime }=-h{\Omega }^{\prime }$. The assumed outward decrease of angular velocity therefore makes $f_1<f_2$, so the mass-exchange part lowers the energy when $d{M}_1>0$: mass is transferred inward. In a coplanar, co-rotating pair, ${\Omega }_1>{\Omega }_2$ and an outward angular momentum transfer has $d{\mathbf{H}}_1=-dH\mathbf{n}$ with $dH>0$. Its contribution is $-({\Omega }_1-{\Omega }_2)dH<0$.

The vector difference also identifies the inclination-reducing exchange. For a fixed small magnitude of admissible angular momentum exchange, its energy contribution is most negative whenby the Cauchy-Schwarz inequality. This gives $d{E}_{\mathrm{angular}}=-s|{\mathbf{\Omega }}_1-{\mathbf{\Omega }}_2{|}^2$. Let $q={\mathbf{n}}_1\cdot {\mathbf{n}}_2=\cos \iota$ and $H_i=|{\mathbf{H}}_i|$, taking fixed particle masses for this exchange. Projecting $d{\mathbf{H}}_i$ perpendicular to its own direction givesThereforeThe downhill exchange reduces the relative orbital inclination. Its inner angular momentum component is ${\mathbf{n}}_1\cdot d{\mathbf{H}}_1=-s({\Omega }_1-{\Omega }_{2}q)<0$; in the aligned limit the outer orbit receives that angular momentum directly. For strongly inclined orbits, part of the exchange instead cancels differently directed angular momenta as the planes align. Together these gradients explain the tendencies toward inward mass transfer, outward angular momentum transport and coplanarity.

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.

The marginal circular orbit has conserved specific angular momentum and specific orbital energyAt fixed $j_0$, the derivative of the effective potential factors asAn inward displacement therefore produces inward radial acceleration; no inner circular equilibrium or centrifugal barrier halts the plunge. In the infinitesimal-displacement limit, energy conservation givesSelecting the inward branch and the original positive sense of rotation yields the Paczyński-Wiita ballistic plungeThe azimuthal speed remains positive while the radius decreases, so the trajectory spirals inward. From any fixed starting radius below $3a$, the time to reach $a$ is finite: $1/|v_r|$ behaves as $(r-a{)}^{1/2}$ near $a$, which is integrable.

The exactly marginal circular state with exactly zero displacement remains a circular solution. A finite inward release from rest changes the energy by a term cubic in the displacement; the displayed expression is its limiting separatrix with $E=E_0$. On that separatrix $|v_r|\propto (3a-r{)}^{3/2}$ near $3a$, so departure from exactly $3a$ takes an infinite limiting time. This distinguishes the infinitesimal-release formula from spontaneous motion of the unperturbed marginal orbit.

For a circular mid-plane orbit, radial force balance in the Newtonian gravitational potential is $r{\Omega }^2={\Phi }_{,r}$. HenceAxisymmetry conserves the specific angular momentum $j=r^2\dot{\varphi }$. At fixed $j$, radial motion obeys $\ddot{r}=-{\partial }_r{\Phi }_{\mathrm{eff}}$, where the effective potential is ${\Phi }_{\mathrm{eff}}=\Phi +j^2/(2r^2)$. Linearize around a circular radius $r_0$, writing $r=r_0+\xi$. Since $j^2=r_0^3{\Phi }_{,r}(r_0,0)$ there,This is the radial epicyclic frequency; equivalently ${\kappa }^2=r^{-3}d(r^4{\Omega }^2)/dr$.

Reflection symmetry gives ${\Phi }_{,z}(r,0)=0$ and ${\Phi }_{,rz}(r,0)=0$, so vertical and radial displacements decouple at linear order. The vertical equation is $\ddot{z}=-{\Phi }_{,zz}(r_0,0)z$, giving the vertical epicyclic frequencyPositive squared frequencies give oscillatory stability; a negative squared frequency gives exponential instability rather than a real oscillation frequency.

For Roche-lobe overflow, $R_2=R_L$. Cubing the Roche lobe approximation and using the donor's mass-radius relation givesKepler's third law therefore yields the Roche-lobe period-mass relation for a linear donor radius law:The component distances from the centre of mass are $a_1=a{M}_2/M$ and $a_2=a{M}_1/M$. Summing their angular momenta gives the circular-binary orbital angular momentumDuring conservative binary mass transfer, $\dot{M}=0$ and ${\dot{M}}_1=-{\dot{M}}_2$. Logarithmic differentiation at constant angular momentum gives $\dot{a}/a=-2(1-q){\dot{M}}_2/M_2$. Thus the Roche-lobe radius response exponent givesTheir difference is $d\log (R_L/R_2)/dt=(2q-8/3){\dot{M}}_2/M_2$. Since the donor loses mass, this is positive for $q<4/3$: the donor shrinks relative to its Roche lobe, so transfer shuts off without an external driver. The Roche lobe itself expands only for $q<5/6$; for $5/6<q<4/3$ it shrinks more slowly than the donor.

Gravitational-wave emission from a binary system provides such a driver: the orbit emits gravitational waves carrying energy and angular momentum, and the shrinking Roche lobe can maintain contact. Allowing $\dot{J}<0$, the binary mass-transfer contact equation becomesIt gives a negative transfer rate for $q<4/3$.

For the classical nova, assume isotropic re-emission from a binary star: ejecta carry the accretor's specific angular momentum $j_1=a_1^2\Omega$. The eruption lasts many orbital periods, so the orbit responds approximately adiabatically and remains nearly circular; an instantaneous asymmetric kick would be a different model. To first order,Differentiating the circular-binary orbital angular momentum formula during the ejection givesThe accompanying Roche lobe change is thereforeThe donor does not change mass in the eruption and moves inside its expanded Roche lobe: this is nova-induced binary detachment.

Put $\Gamma =-\dot{J}>0$ for the constant external loss. During detachment the component masses are fixed, so the Roche lobe shrinks at fractional rate $2\Gamma /J$. Closing the fractional gap $4\delta m/(3M)$ takes $t_d=2J\delta m/(3M\Gamma )$. During the subsequent semidetached binary phase, the binary mass-transfer contact equation gives $|{\dot{M}}_2|=3M_2\Gamma /[J(4-3q)]$, so accumulating the next eruption's mass takes $t_s=J(4-3q)\delta m/(3M_2\Gamma )$. ConsequentlyBoth times are evaluated to first order in $\delta m/M$, with the same equilibrium donor mass-radius relation; changes of $q,J$ and the transfer rate within one cycle affect higher orders.

More than three integrals occur when additional symmetry or degeneracy restricts motion beyond a generic three-frequency invariant torus. In a time-independent central gravitational potential, energy and the three components of specific angular momentum are four generically independent isolating integrals. They fix the orbital plane as well as the radial-motion parameters. A nonclosed planar rosette typically has two-dimensional orbit closure, rather than a three-dimensional invariant torus.

The Kepler orbit is a stronger example. In addition to $E$ and $\mathbf{L}$, the conserved Laplace-Runge-Lenz vectorfixes the apsidal direction. Its constraints $\mathbf{A}\cdot \mathbf{L}=0$ and $A^2={\mu }^2+2E{L}^2$ leave five independent integrals among $E,\mathbf{L},\mathbf{A}$. It is a maximally superintegrable Hamiltonian system, with closed nondegenerate bound orbits.

More than three integrals need not form a set of Poisson-commuting functions; the three-integral limit applies to an independent commuting set, rather than to every conserved orbit constraint. At a nonstationary point in six-dimensional phase space, integral gradients annihilate the nonzero flow vector, so at most five are independent. This explains the isolating integrals and orbital closure dimension count. A resonance on one isolated orbit does not by itself prove an additional smooth global integral on a neighborhood.

In cylindrical coordinates the particle's radial and vertical equations are $\ddot{r}-r{\dot{\varphi }}^2=-{\Phi }_r$ and $\ddot{z}=-{\Phi }_z$. Axisymmetry gives conservation of angular momentum, with specific angular momentum $h=r^2\dot{\varphi }$. Eliminating the cyclic azimuthal motion therefore givesThus $\boxed{{\Phi }_1=\Phi +h^2/(2r^2)}$ governs the meridional motion. Its conserved energy is $\frac{1}{2}({\dot{r}}^2+{\dot{z}}^2)+{\Phi }_1$. The additional term is the centrifugal barrier appropriate to fixed angular momentum.

For a point mass, $\Phi =-GM/\sqrt{r^2+z^2}$. If $h\ne 0$, the vertical derivative vanishes only at $z=0$, and the radial stationary condition there gives $r_0=h^2/(GM)$. This is the sole stationary point with $r>0$, at energy $-GM/(2r_0)$. Its Hessian matrix isConsequently the local contours are closed, approximately circular curves about a strict minimum. Finite negative-energy contours above the minimum remain closed; the centrifugal barrier prevents approach to the axis and the negative energy prevents escape to infinity. The zero contour satisfies $z^2=4r^4/r_0^2-r^2$, while positive-energy contours can be unbounded. The left panel illustrates these features; portions of the outer contours extend outside the plotting window.

Each panel uses its own circular radius as length unit and $GM/r_0$ as potential unit. Linearizing near the minimum givesThe radial and vertical epicyclic frequencies both equal the circular orbital frequency. Circular Kepler orbits are stable to small meridional perturbations of a freely moving particle. With no dissipation the oscillations persist rather than decay. A change in h shifts the circular radius to a nearby orbit; azimuthal phase itself is neutral by rotational symmetry. If $h=0$, there is no stationary point or centrifugal barrier and this circular-orbit stability discussion does not apply.

Use outward radial mass flux $\mathcal{F}=2\pi r\Sigma u_r$ and outward stress-carried angular-momentum flux $\mathcal{G}$. The PDF includes a second conservation equation omitted from the converted TeX:Together with mass conservation, ${\Sigma }_t+(2\pi r{)}^{-1}{\mathcal{F}}_r=-S$, it fixes the signs of the wind terms. Subtract h times the mass equation, noting that the prescribed $h(r)$ is independent of time. One obtainsThe viscous torque in an accretion disk in this outward-flux convention is $\mathcal{G}=-2\pi \bar{\nu }\Sigma r^3{\Omega }^{\prime }$. Thus, where $h^{\prime }\ne 0$,Putting this result back into mass conservation givesThe derivation requires a fixed rotation law; it would acquire additional terms if h evolved independently in time.

For Keplerian rotation, $h=\sqrt{GMr}$, $h^{\prime }=h/(2r)$, and $r^3{\Omega }^{\prime }=-(3/2)h$. With $T=qSh$, the radial flux becomesConsequently the wind-loaded Keplerian disk evolution equation isThe final term removes mass locally. The extra flux term describes the disc's radial response when the outflow removes more or less than the local specific angular momentum. In particular, $q=1$ cancels this additional transport term without canceling the local mass sink.

Inside the uniform sphere, $V_c^2=3{\sigma }^2{r}^2/r_0^2$. Thus, with $X=V_c/(\sqrt{2}\sigma )$, part (b) givesThe Gaussian integral needed for the lowered Maxwellian velocity distribution isSubstitution $v=\sqrt{2}\sigma y$ therefore yieldsSince $V_c<v_e$ throughout the core, the integrand is positive and so is $I(X)$ for $X>0$. Using $m{n}_0={\rho }_0$, part (a)'s massive-test-body drag magnitude isThe slow drift of a circular orbit is determined by torque, not by equating its changing circular speed to the tangential drag. Here $V_c=\omega r$, with constant $\omega =\sqrt{4\pi G{\rho }_0/3}$. The specific angular momentum $j=r{V}_c=\omega r^2$ obeysHence the positive local decay timescale isThere is a genuine factor-of-two discrepancy with the requested denominator $32{\pi }^2$. That denominator would give $V_c/(2D)$, which is the timescale for a flat galaxy rotation curve, not this uniform-density core. It also follows from the incorrect replacement ${\dot{V}}_c=-D$, which drops the $V_c\dot{r}$ term in $d(r{V}_c)/dt$. The literal signed definition $r/(2\dot{r})$ is negative for an inward drift; it equals $-{\tau }_{+}$. A positive decay time requires the absolute value or a minus sign.