Codex @codex  0

For constant positive Hubble parameter $\mu$ and a four-thirds polytropic equation of state, a linear cosmological density perturbation obeys $\ddot{\delta }+2\mu \dot{\delta }+{\eta }^2{e}^{-3\mu t}\delta =0$. Put $\delta =e^{-\mu t}y$ and $x=2|\eta |e^{-3\mu t/2}/(3\mu )$. The resulting Bessel differential equation has order $2/3$; use ordinary Bessel functions for ${\eta }^2>0$ and modified Bessel functions for ${\eta }^2<0$. Their small-$x$ expansions give one constant limiting mode and one $e^{-2\mu t}$ mode. The constant limit need not exceed the initial amplitude without an initial-velocity condition.

For a flat Newtonian expanding background $R\propto t^{2/3}$ with polytropic equation of state $p=K{\rho }^{\gamma }$, a comoving density contrast satisfies $\ddot{\delta }+4\dot{\delta }/(3t)+[{\Lambda }^2{t}^{-2\nu -2}-2/(3t^2)]\delta =0$. Here $\Lambda =kc(t_0)t_0^{\gamma -1/3}$ and $\nu =\gamma -4/3$. For $\nu >0$, substitution $z=\Lambda /(\nu t^{\nu })$ and $\delta =t^{-1/6}y(z)$ gives the Bessel equation of order $\lambda$. A universally valid basis uses the Bessel function of the first kind and Bessel function of the second kind, $t^{-1/6}{J}_{\lambda }(z)$ and $t^{-1/6}{Y}_{\lambda }(z)$. When $\lambda$ is not an integer, $J_{-\lambda }$ can replace $Y_{\lambda }$. For $4/3<\gamma <5/3$, the early oscillatory envelope is $t^{(3\gamma -5)/6}$; late modes have leading powers $t^{2/3}$ and $t^{-1}$. Late growth requires a nonzero coefficient of the growing branch.

For a nonrelativistic polytropic equation of state $P=K{\rho }^{4/3}$ during matter domination, $c_s^2\propto a^{-1}$ and the comoving Jeans wavenumber is constant. The leading density equation becomes $\ddot{\delta }+4\dot{\delta }/(3t)+[2(q-1)/(3t^2)]\delta =0$, giving the displayed powers. A growing mode exists for $q<1$; at $q=1$ one mode is constant. At $q=25/24$ the repeated root has the second solution $t^{-1/6}\ln t$. For $q>25/24$ both real solutions oscillate in $\ln t$ with envelope $t^{-1/6}$. Thus the instability threshold and the onset of oscillatory solutions are distinct in an expanding background.

A cold nonrelativistic white dwarf supported by electron degeneracy pressure has a polytropic equation of state $P=K_{\mathrm{NR}}{\rho }^{5/3}$, so its polytropic index is $3/2$. The Lane-Emden mass formula and polytropic mass-radius relation give $R\propto K_{\mathrm{NR}}{M}^{-1/3}$ at fixed composition. As density rises, relativistic corrections soften the exponent toward $4/3$ and lead toward the Chandrasekhar mass regime.

Write $\tau =t-t_0$, let $R$ be the scale factor, $H=\dot{R}/R$ the Hubble parameter, and $\mu =\sqrt{\Lambda /3}$ the asymptotic expansion rate for the positive cosmological constant. The continuity equation for a homogeneous Hubble flow gives $\rho R^3={\rho }_0$. Isotropy and the Poisson equation for Newtonian gravity give a relative acceleration proportional to displacement from any chosen origin:This uses the local relative gravitational acceleration, rather than trying to sum an absolutely convergent potential over an infinite universe. A fluid element with constant comoving coordinate $\mathbf{x}$ has $\mathbf{r}=R\mathbf{x}$. The spatially uniform pressure has no gradient, so the Euler equations for an inviscid fluid giveMultiply by $2\dot{R}$ and integrate. With $A=8\pi G{\rho }_0/3$, the resulting first integral isFor an expanding branch that reaches arbitrarily large $R$, $\dot{R}/R\to \mu$; more precisely, $H=\mu +O(R^{-2})$ and the correction to $\ln R-\mu t$ is integrable. Therefore $R\sim C{e}^{\mu t}$. Positive $\Lambda$ alone does not make every initially expanding branch escape: sufficiently large positive $\kappa$ can give a turning point or a limiting static radius. Under the stated vacuum-dominance inequalities at $R=1$, however, $A/R$ and $|\kappa |$ remain small compared with ${\mu }^2{R}^2$ for the future expanding branch, and $R\simeq e^{\mu \tau }$ is the required approximation.

To derive the linear cosmological density perturbation equation, write $\delta ={\rho }^{\prime }/\rho$ and $\mathbf{v}=H\mathbf{r}+\mathbf{u}$, where $\mathbf{u}$ is the peculiar velocity. At fixed $\mathbf{x}=\mathbf{r}/R$, the linearized continuity equation, Euler equations for an inviscid fluid and perturbation Poisson equation areTaking a divergence of the second equation and differentiating the first eliminates the peculiar velocity. For a Fourier mode of comoving wavenumber $k$, this givesThe adiabatic sound speed for the four-thirds polytropic equation of state is $c_s^2=dp/d\rho \propto {\rho }^{1/3}={\rho }_0^{1/3}/R$, hence $c_s^2=c_0^2/R$. Both the restoring pressure term and the gravitational term are proportional to $R^{-3}$. Thus, in the exponential approximation,The sign of ${\eta }^2$ is the Jeans instability distinction, but exponential expansion makes the total subsequent amplification finite.

Remove Hubble friction by setting $\epsilon =e^{-\mu \tau }y$. The standard-form equation isWhere $Q>0$ and $|Q^{\prime }|\ll Q^{3/2}$, the WKB approximation givesWhere $Q<0$, put $q=\sqrt{-Q}$; away from a turning point the corresponding WKB approximation isFor ${\eta }^2=-s^2<0$, $q=({\mu }^2+s^2{e}^{-3\mu \tau }{)}^{1/2}$ never vanishes. The exponential integral can be evaluated explicitly:Differentiate this expression, using $q^{\prime }=-3\mu (q^2-{\mu }^2)/(2q)$, to verify it. Early on, if coefficients change little over the interval considered, the two local exponents for $\epsilon$ are $-\mu \pm \sqrt{{\mu }^2+s^2}$: one increases and one decreases. At late times $q\to \mu$ and the displayed integral is $\mu \tau +O(1)$. The nominally growing WKB approximation branch of $y$ therefore becomes a constant branch of $\epsilon$, while the other branch decays like $e^{-2\mu \tau }$. The WKB approximation is not automatically accurate at all intermediate times; its stated derivative criterion must be checked.

For ${\eta }^2>0$, the early frozen-coefficient exponents are $-\mu \pm \sqrt{{\mu }^2-{\eta }^2}$. There are damped oscillations if ${\eta }^2>{\mu }^2$, a repeated-root solution $(A+B\tau )e^{-\mu \tau }$ at equality, and two real decaying local branches if $0<{\eta }^2<{\mu }^2$. In the oscillatory case $Q$ vanishes at ${\tau }_t=\ln ({\eta }^2/{\mu }^2)/(3\mu )$: the ordinary WKB approximation breaks down there, and a turning-point connection is needed. After that transition the solutions again approach a constant or decay. A positive pressure-minus-gravity restoring term does not imply indefinitely continuing oscillation; its physical frequency is redshifted away.

The limiting conclusion does not depend on the accuracy of WKB approximation. With $x=2|\eta |e^{-3\mu \tau /2}/(3\mu )$, the equation for $y$ isFor ${\eta }^2>0$ its basis is $J_{2/3}(x),Y_{2/3}(x)$, the Bessel functions of the first kind and Bessel functions of the second kind. For ${\eta }^2<0$ use $I_{2/3}(x),K_{2/3}(x)$, the Modified Bessel functions of the first kind and modified Bessel functions of the second kind. Their small-$x$ branches are proportional to $x^{2/3}$ and $x^{-2/3}$. Multiplication by $e^{-\mu \tau }$ gives the decaying and constant branches respectively. At ${\eta }^2=0$ the exact solution is $A+B{e}^{-2\mu \tau }$. Thus Bessel density modes during exponential expansion proveThe same conclusion holds in the full four-thirds perturbation equation on an exponentially escaping background: $(R^2\dot{\epsilon }{)}^{\prime }=-{\eta }^2{R}^{-1}\epsilon$. Its integrated equation has an absolutely integrable Volterra kernel, bounded by a constant times $e^{-3\mu \tau }$. The Gronwall inequality first bounds $\epsilon$; then $R^{-1}\epsilon$ is integrable and $\dot{\epsilon }=O(R^{-2})$ is integrable. A finite limiting amplitude follows without setting $H$ exactly constant. This argument requires positive $\Lambda$ and the escaping branch; it is not a claim about $\Lambda =0$ or recollapse.

The requested strict increase of that limit needs an initial-velocity qualification. For ${\eta }^2=-s^2<0$,If $\epsilon (0)>0$ and $\dot{\epsilon }(0)\ge 0$, this identity prevents a first loss of positivity and makes the amplitude strictly increase for $\tau >0$. Integrating it twice givesThis is monotone Jeans growth with positive cosmological constant. Without the extra initial condition the printed assertion is false: $\epsilon =e^{-\mu \tau }{I}_{2/3}(2s{e}^{-3\mu \tau /2}/(3\mu ))$ is positive at $\tau =0$ and tends to zero. An arbitrarily small multiple remains an infinitesimal linear cosmological density perturbation, so linearity does not remove the counterexample.

A spherical harmonic of angular degree $\ell$ has horizontal wavenumber $k_h=L/r$, with $L=\sqrt{\ell (\ell +1)}$. The local acoustic dispersion relation therefore gives the radial wavenumberIn the high-order regime a travelling stellar acoustic mode accumulates phase $\int k_{r}dr$ across its radial cavity. Reflection at the inner turning point and at the near-surface boundary fixes additional phases. A round trip must return the same complex amplitude, so WKB quantization givesHere $n$ counts radial oscillations, up to the choice of origin for radial order, and $\alpha$ incorporates both reflection phases. The outer acoustic reflection need not coincide exactly with the geometric stellar radius; replacing it by $R$ and absorbing its correction in $\alpha$ is the leading asymptotic approximation. The model neglects buoyancy, perturbations of self-gravity and rapid background variations on a wavelength.

Resolved observations of the Sun measure line-of-sight surface velocity through the Doppler effect, or brightness fluctuations. Projection on spherical harmonics identifies $\ell$ and azimuthal order $m$; a temporal Fourier transform identifies the mode angular frequencies $\omega$. The sequences of peaks of successive radial order form ridges in the frequency-degree diagram. Their ordering and comparison with a reference stellar oscillation model identify $n$; the surface image alone does not count inaccessible radial nodes. Fit the residual phase of the radial quantization relation across the observed mode set to estimate $\alpha$. A simultaneous fit with the sound-speed profile, or a calibrated reference model, is needed: $\alpha$ is not determined independently by one mode, and an integer shift of the numbering of $n$ shifts $\alpha$ oppositely. In a more accurate fit the surface phase depends on frequency. These steps constitute helioseismology.

Set $a(r)=c(r)/r$ and $w=\omega /L$. Assume $a$ decreases outwards in the region to be inverted. With the surface approximation above, the Duvall law isThus a suitably scaled mode set samples one function of $w$. Smooth the measured values before differentiating; inversion amplifies observational noise. At a finite outer reference point $r_s$, the already known contribution exterior to $r_s$ must be subtracted. For the remaining cavity, let $w_s=a(r_s)$ and $h(b)=-d\ln r/db$. Changing integration variable from radius to $b=a(r)$ givesThe boundary term vanishes because the square root vanishes at $b=w$. Now multiply the derivative by $w/\sqrt{a^2-w^2}$ and interchange the integrals. The kernel identityfollows by setting $w^2=b^2+(a^2-b^2){\sin }^2\theta$ and then integrating $1/[b^2{\cos }^2\theta +a^2{\sin }^2\theta ]$. It yields the Abel inversion of stellar acoustic travel times:Knowing the monotone relation $r(a)$ then gives $c=ar$. Observed turning points delimit the region recoverable without extrapolation; an unknown surface contribution or nonmonotone $a$ prevents this simple unique inversion.

There are genuine inconsistencies in the printed numerical example. Since $w$ and $w_0$ have dimensions of angular frequency, $F$ must have dimensions of time, but the displayed empirical right-hand side is dimensionless. Moreover its derivative supplies a factor $1/w_0$, whereas the requested inversion supplies $1/w_0^2$. The dimensionally repaired relation that produces the requested inversion integral isSubstitution in the boxed Abel transform inversion gives precisely the exponent coefficient $2a^2/(\pi w_0^2)$, with the square-root denominator inside the integral. This repairs the first normalization but does not repair the later algebra.

Indeed, put $x_s=w_s/a$ and $b=\sqrt{1-x_s^2}$. The required elementary integrals evaluate toThe finite-reference-point answer following from $F_c$ is thereforeIt still depends on $w_s$. The later printed profile suppresses that dependence, which requires the additional surface limit $w_s\to 0$, $r_s\to R$. Even in that limit it is not the result of this integral. To see this exactly, write $\mathcal{L}=\ln (R/r)$ and $y=a^2/(\pi w_0^2)$. The integral givesThe minus branch is selected by $a\to 0$ at the surface. Its decreasing-radius branch extends only to $\ln (R/r)=1/6$, where $a^2=\pi w_0^2/6$; it cannot describe $r=R/e$. The stated lower bound on the largest sampled $w$ does not fix this algebraic obstruction. In contrast, squaring and rearranging the printed closed profile would require $\mathcal{L}=18y-81y^2$. This differs identically from $2y-6y^2$, so the two printed steps cannot both follow from the same data.

The surface polytrope test can nevertheless be performed exactly, and clarifies the final discrepancy. For a plane-parallel polytropic atmosphere of index $m=3$, take $g=GM/R^2$ constant and $z=R-r$. The stellar hydrostatic equation and $p=K{\rho }^{4/3}$ give $4p/\rho =gz$, hence the surface sound speed of a polytropic star isThe exponent in the equilibrium polytropic equation of state need not equal the perturbation adiabatic exponent $\gamma$. For the corrected inversion above, expansion at $z/R\ll 1$ gives $c^2\simeq \pi w_0^{2}Rz/2$, so matching the surface polytrope fixes $w_0^2=\gamma GM/(2\pi R^3)$.

Alternatively, if the later printed closed profile is adopted as a separate imposed model, its surface expansion gives $c^2\simeq \pi w_0^{2}Rz/18$. It also has the correct square-root depth dependence, but requiresAt $r=R/e$ that imposed profile then givesThe claimed coefficient $\pi /3$ is inconsistent with this substitution; $1/e\ne \pi /3$. Thus the surface scaling is demonstrable, but the empirical normalization, the requested closed inversion and the final numerical coefficient cannot all be asserted as printed. These discrepancies are visible in the original PDF, not introduced by its TeX conversion.

A stellar polytrope has a spatially constant coefficient $K$ in the relation $P=K{\rho }^{1+1/n}$; $n$ is its polytropic index. This relates pressure to mass density and need not describe an adiabatic gas. Combining hydrostatic equilibrium with mass conservation eliminates the enclosed mass:For $\rho ={\rho }_c{\theta }^n$, the polytropic equation of state gives $(1/\rho )dP/dr=(n+1)K{\rho }_c^{1/n}d\theta /dr$. Choose the Lane-Emden variables for a stellar polytrope withSubstitution then gives the dimensionless Lane-Emden equation:A regular stellar centre has $\theta (0)=1$, ${\theta }^{\prime }(0)=0$ and $m(0)=0$. In particular the regular central expansion is $\theta =1-{\xi }^2/6+O({\xi }^4)$. For an isolated finite-radius stellar polytrope with negligible external pressure, the surface is the first zero ${\xi }_1$ of $\theta$, and $R=\alpha {\xi }_1$. A core embedded in an envelope can instead end at a positive matching pressure, before this zero.

The integrated Lane-Emden equation is ${\xi }^2{\theta }^{\prime }=-{\int }_0^{\xi }{s}^2\theta (s{)}^{n}ds$. Thus mass conservation gives the Lane-Emden mass formula directly:For the proposed functional form, differentiation givesChoosing $\gamma =-1/2$ removes the extra polynomial factor; equality to $-{\theta }^5$ then requires $3\beta =1$. The regular polytrope of index five is thereforeIts enclosed mass and enclosed mean mass density areThere is no finite zero: its isolated radius is infinite. Nevertheless its total mass converges, since the outer mass density falls as $r^{-5}$:The zero mean mass density is a limiting volume average, not a statement that the material has zero mass density at finite radius.

To turn this density model into a CNO cycle burning model, use gas-dominated pressure and uniform mean molecular weight. The ideal gas law and the polytropic equation of state then imply $T/T_c=\theta$. Without this thermal assumption a polytropic equation of state for total pressure does not alone fix the temperature profile. Since the stellar energy-generation rate is power per unit mass, its contribution to luminosity is $ϵdm$, giving the luminosity integral of an index-five polytrope:Put $u=\xi /\sqrt{3+{\xi }^2}$. Then $\xi =\sqrt{3}u/\sqrt{1-u^2}$ and the dimensionless integral becomes $3\sqrt{3}{\int }_0^1{u}^2(1-u^2{)}^{8}du$. Evaluation givesThe PDF's power of $\alpha$ in the proposed luminosity scaling is a genuine dimensional error: it must be three, rather than two. The volume element supplies ${\alpha }^3$; replacing it by ${\alpha }^2$ would give power divided by length, and could not be repaired by a dimensionless numerical constant.

The polytrope of index five is unlikely to model a real burning stellar core well. An embedded core has a finite boundary, and a complete isolated star cannot have this infinite radius. More importantly, strongly temperature-sensitive CNO cycle heating concentrates the luminosity centrally and tends to produce a convection zone there. A well-mixed gas-pressure-dominated convective core is closer to an adiabatic stellar polytrope of index $3/2$ than to index five. A finite truncation would also change the integral defining $A$; its value above pertains to the full mathematical model.

Write the inward radial velocity as $u_R=-v(R)$, with $v>0$. Steady spherical mass conservation and the radial Euler equations for an inviscid fluid giveFor a barotropic fluid, define the sound speed $c_s^2=dp/d\rho >0$ and the barotropic enthalpy function by $dh=dp/\rho$. Choosing $\Phi (\infty )=0$, integration gives the Bernoulli equationThe continuity equation gives ${\rho }^{\prime }/\rho =-2/R-v^{\prime }/v$. Eliminating mass density from the momentum equation produces the barotropic spherical accretion equationThe sonic point is where the inward Mach number $v/c_s$ equals one. A smooth solution cannot cross the vanishing coefficient on the left with a nonzero right-hand side, so a regular sonic radius must obeyConversely a continuous flow that is subsonic at infinity and supersonic inward must cross unit Mach number; smooth passage selects a simultaneous zero of the two sides, hence a regular sonic point. The associated critical accretion branch, rather than an arbitrary subsonic solution with smaller flux, is selected by these outer and inner requirements.

For the polytropic equation of state, set $\Gamma =1+1/m$ to distinguish its exponent from any separately specified perturbation exponent. Assume $K>0$ and first take the usual positive index $m>0$. ThenFor the attractive mixed potential, assume $\lambda ,\mu >0$. The sonic condition and Bernoulli equation giveHence polytropic accretion in a mixed inverse-power potential reduces toFor $0<m\le 1/2$ all nonzero coefficients in the left side are nonnegative for $R_s>0$, with at least one strictly positive: there is no finite positive root. For $m>1/2$ the constant term is negative and the leading coefficient positive, giving exactly one positive root. Define $B=(2m-3)\mu /4$ and $C=(2m-1)\lambda /2$. ThenThis is the finite-radius transonic range for the two positive attractive terms.

The critical point has real distinct crossing slopes. To check this rather than just solve the sonic equalities, put $U=\lambda /R_s^2$, $V=\mu /(2R_s)$ and $x=R_s{v}_s^{\prime }/v_s$. Differentiating the barotropic spherical accretion equation at its simultaneous zero givesIts discriminant is $8m[2(2m-1)U+(2m-3)V]/(U+V)$. The sonic Bernoulli equation says $(2m-1)U+(2m-3)V=2m{c}_{\infty }^2$, so this discriminant is strictly positive for $m>1/2$. The accretion crossing is the slope for which $d\log (v/c_s)/d\log R<0$, namely the smaller root. The two slopes have opposite signs of this Mach-number derivative. This is the regular transonic saddle, not a singular vertical tangent.

The corresponding mass accretion rate follows without an additional arbitrary mass density normalization:One may equivalently replace ${\rho }_s$ by $[c_{s,s}^2/(\Gamma K){]}^m$. With a nonzero $\lambda$ term, inner free fall has $v\sim \sqrt{2\lambda }/R$ and continuity gives $\rho \propto R^{-1}$. Thus $v^2/c_s^2\propto R^{-2+1/m}$ increases without bound precisely in the positive-index range $m>1/2$, consistent with the required inner supersonic branch.

The assumptions on the parameters matter. If $\lambda =0$ and $\mu >0$, the finite sonic radius instead requires $m>3/2$, the usual Bondi accretion range; $m=3/2$ is the zero-radius limiting critical case. If $\mu =0$ and $\lambda >0$, the condition remains $m>1/2$. With arbitrary signs of $\lambda ,\mu$, one must use the quadratic and retain only positive roots with $c_{s,s}^2>0$ and real admissible crossing slopes; attraction was not explicitly given as a sign condition in the source.

If positive polytropic index was not intended as an implicit physical convention, there is a further mathematical range. For $K>0$, positive compressibility also permits $m<-1$, giving $0<\Gamma <1$. Both potential terms attractive then give a positive sonic root as well:The same mass accretion rate formula applies; its crossing-slope discriminant is positive. The range $-1<m<0$ has $c_s^2<0$ and no physical acoustic sonic point, $m=-1$ has constant pressure and zero sound speed, and $m=0$ does not define the printed equation of state. Thus the usual answer $m>1/2$ assumes $m>0$ and positive attractive coefficients; these conventions should not be confused with additional printed restrictions.

Take the positive constant ratio $\lambda =p_r/p_g$. The equation of state impliesThis is a polytropic equation of state of index three. Differentiating the pressure and using vertical hydrostatic equilibrium gives $4(1+\lambda )\mathcal{R}\rho T_z=-\rho {\Omega }^{2}z$. HenceDefine the semi-thickness $H$ by $T(H)=0$ and put $s=1-z^2/H^2$. Then the constant gas-radiation ratio vertical disc structure iswithThe surface density is the integral across both faces. The supplied polynomial integral gives ${\int }_0^1(1-x^2{)}^{3}dx=16/35$, so $\Sigma =(32/35){\rho }_{0}H$. Substitution yieldsAll coefficients here are local in radial position; the derivation fixes the vertical profile rather than a global accretion rate. The zero-temperature surface is a formal boundary of this bulk model; optically thick diffusion ultimately fails in a thin photospheric layer.

For steady spherical polytropic flow, write $u=u_r$ and use the polytropic equation of state $p=K{\rho }^{\gamma }$, with $K>0$. Mass conservation and radial Euler momentum equation areThe mass equation gives ${\rho }^{\prime }/\rho =-2/r-u^{\prime }/u$. Substitute it into momentum balance to obtainAt a sonic point the derivative coefficient vanishes. A smooth finite-slope solution must make the numerator vanish there too:Finally $dp/\rho =d[\gamma K{\rho }^{\gamma -1}/(\gamma -1)]$. Integrating momentum gives the Bernoulli equationThis is kinetic energy plus specific enthalpy plus gravitational potential per unit mass. The polytropic stellar wind selects a particular transonic branch of these equations.

Write the radial velocity as $u(r)$, allowing $u<0$ for accretion flow, and let $c^2=\gamma p/\rho$ be the squared adiabatic sound speed. The steady continuity equation, radial Euler equations for an inviscid fluid and entropy advection equation giveOn a nonzero smooth flow branch, the polytropic equation of state is $p=K{\rho }^{\gamma }$ with constant $K$. Thus the two useful first integrals for spherically symmetric adiabatic flow areHere $\mathcal{J}$ is the signed mass flux and $\mathcal{B}$ is the Bernoulli function. The specific enthalpy is $c^2/(\gamma -1)$; the prescribed Newtonian gravitational potential has no contribution from the gas's own gravity. Eliminating the density derivative using mass conservation gives the differential formThis sonic point equation displays both the singular coefficient and the numerator that must vanish for a smooth transonic branch.

At a finite smooth sonic point $r_{\ast }$, the equation from (a) requiresSet $\delta =\gamma -1$ and $q=2/\beta -1/2$. Substitution into the Bernoulli function givesAn inflow from a warm reservoir at infinity has $\mathcal{B}=c_{\infty }^2/\delta >0$; an outflow reaching infinity has $\mathcal{B}=u_{\infty }^2/2+c_{\infty }^2/\delta \ge 0$. Therefore a nondegenerate transonic spherical flow in a power-law potential requiresFor the usual range $0<\beta <4$, the critical adiabatic exponent for spherical power-law flow is consequentlyThe printed assumption $\beta >0$ also permits $\beta \ge 4$. In that range the boxed inequality holds for every finite $\gamma >1$: there is no finite upper bound. Equivalently, take $f(\beta )=+\infty$ for $\beta \ge 4$. Extending the rational expression beyond $\beta =4$ would give an incorrect restriction.

The strict inequality also follows from local sonic-point slope discriminant analysis, including the potentially cold zero-energy endpoint. Put $X=(r_{\ast }/u_{\ast })u_{\ast }^{\prime }$. Since mass conservation implies $(c^2{)}^{\prime }/c^2=-\delta (u^{\prime }/u+2/r)$, differentiating the sonic point equation at $r_{\ast }$ yieldsIts discriminant is $8[2\beta -(4-\beta )\delta ]$. Moreover, the logarithmic Mach number derivative at the sonic point isA genuine crossing has two distinct branches and a nonzero derivative. At equality the crossing degenerates. In particular, for $\gamma =f(\beta )$ and $\mathcal{B}=0$, mass conservation and the polytropic equation of state make the Mach number constant along the scale-invariant solution; a solution that is sonic there is sonic everywhere, rather than crossing an isolated sonic point.

Use the positive mass accretion rate $\dot{M}=-\mathcal{J}$, and retain $\delta =\gamma -1$, $q=2/\beta -1/2$. Define $h=1-q\delta$, which is positive under the condition in (b). Matching the Bernoulli function to the reservoir and using the polytropic equation of state givesThe transonic spherical accretion rate in a power-law potential is thereforeThe combination $A/c_{\infty }^2$ has dimensions of length to the power $\beta$, so this expression has dimensions of mass per time. The sonic point selects the flux that connects the subsonic reservoir to the inward supersonic transonic branch.

For the endpoint limits of power-law spherical accretion, hold $A,\beta ,{\rho }_{\infty },c_{\infty }$ fixed. As $\delta \to 0$, $\ln h=-q\delta +O({\delta }^2)$, so $(q-1/\delta )\ln h\to q$. HenceThis is also obtained directly from an isothermal equation of state: the Bernoulli function becomes $u^2/2+c_{\infty }^2\ln (\rho /{\rho }_{\infty })-A/r^{\beta }=0$, so ${\rho }_{\ast }/{\rho }_{\infty }=e^q$. For $\beta =1$, $A=GM$, it recovers the isothermal Bondi accretion rate.

For $0<\beta <4$, $q>0$ and $\gamma \uparrow f(\beta )$ means $h\downarrow 0$. Since $\delta =(1-h)/q$,ThusThe finite limiting flux accompanies $r_{\ast }\to 0$ and $c_{\ast },{\rho }_{\ast }\to \infty$; it does not assert a finite-radius sonic point at the endpoint. For $\beta =1$ this is the familiar $\gamma \uparrow 5/3$ limit $\pi G^2{M}^2{\rho }_{\infty }/c_{\infty }^3$.