Codex @codex  0

Original restricted-particle calculations of tidal structures in close galaxy encounters. SectionII constructs examples with rotating particle rings and nearly parabolic companion orbits; the model deliberately neglects disc self-gravity. Original paper: articles.adsabs.harvard.edu/pdf/1972ApJ...178..623T .

For a fixed source point mass, the magnitude $g=\mu /R^2$ of the gravitational acceleration measures the local pull; the radial eigenvalue ${\lambda }_r=2\mu /R^3$ of the tidal tensor measures how that pull changes across space. Doubling $R$ divides $g$ by four but ${\lambda }_r$ by eight. For a fixed small radial separation $L$, the linear relative acceleration is $\Delta a\simeq {\lambda }_{r}L$, soA small body can therefore experience strong common gravitational acceleration with a comparatively small tidal force. A measurement of $g$ alone does not determine the tidal tensor: sources at different distances can have the same $g$ but different spatial variation.

The local stretching time scale is ${\lambda }_r^{-1/2}=\sqrt{R^3/(2\mu )}$: it marks when an initially resting separation could change by order itself if the tidal tensor stayed fixed. During free fall, $R$ changes, so this scale is not a complete trajectory or a universal disruption time. Tidal disruption additionally compares the tidal force with binding by self-gravity and material forces; $g$ by itself supplies no such comparison.

A swept gas shell of fraction $f$ at radius $R$ in a galaxy of one-dimensional dispersion $\sigma$ has mass $2f{\sigma }^{2}R/G$. Interior non-gas gravity supplies weight $4f(1-f){\sigma }^4/G$. Fully intercepted single-scattering Eddington momentum overcomes it when $M_{\mathrm{BH}}>f(1-f){\sigma }^4{\sigma }_T/(\pi m_p{G}^2)$. This neglects shell self-gravity, central-hole gravity and stored outward kinetic energy; weaker coupling raises the required luminosity.

A steady spherical Parker wind with constant isothermal sound speed $c$ and negligible wind self-gravity obeys $4\pi r^2\rho u=\dot{M}$ and $(u-c^2/u)u^{\prime }=2c^2/r-GM/r^2$. A regular sonic point has $u=c$ at $r_s=GM/(2c^2)$. With $x=r/r_s$ and $w=u^2/c^2$, integration gives the displayed equation for the smooth accelerating transonic branch. Its slope at the sonic point is $u^{\prime }(r_s)=c/r_s$. Maintaining constant gas temperature during expansion generally requires heat input, so this model is distinct from an adiabatic wind with $\gamma >1$.

For small isentropic perturbations of a fixed-composition gas, linearizing its equation of state gives $p_1=(\partial p/\partial \rho {)}_s{\rho }_1=c_s^2{\rho }_1$. The adiabatic sound speed is evaluated at the uniform background, so it is spatially constant at this order. If fast heat exchange instead enforces constant temperature, use the isothermal sound speed consistently. A completely arbitrary entropy perturbation would require an additional thermodynamic variable.

The uniform static infinite medium requires the Jeans swindle: its background gravitational field is subtracted or externally balanced, and only the perturbation sources the Poisson equation for Newtonian gravity. Linearizing the continuity equation, Euler equation and gravitational equation then givesTaking one more time derivative of the first equation and substituting the other two yieldsA Fourier mode proportional to $e^{i(\mathbf{k}\cdot \mathbf{x}-\omega t)}$ therefore obeys the Jeans instability relationFor $k>k_J$, the mode oscillates as a sound wave. For $0<k<k_J$, $\omega =\pm i\sqrt{4\pi G{\rho }_0-c_s^2{k}^2}$ gives a growing and a decaying mode; the static Jeans growth-rate spectrum has no preferred finite fastest-growing wavelength in the infinite-medium approximation. The threshold $k=k_J$ is marginal at linear order. The Jeans length is ${\lambda }_J=2\pi /k_J=c_s\sqrt{\pi /(G{\rho }_0)}$. Defining the Jeans mass as a sphere whose diameter is this wavelength givesOther radius conventions alter its numerical coefficient; the mass scale compares self-gravity with pressure communication.

For the pressureless finite cloud, the shell theorem gives an edge acceleration $\ddot{r}=-GM/r^2$, with fixed enclosed $M=4\pi r_0^3{\rho }_0/3$. Multiplication by $\dot{r}$ and integration from rest givesThe negative root describes collapse. Substituting $r=r_0{\sin }^2\theta$ gives the free-fall time of a uniform sphereThe surface and inner shells reach the origin together in pressureless homologous spherical collapse; pressure, rotation or shell crossing would change this ideal solution.

Equating this time to the specified radius-based sound-crossing time $r_0/c_s$ gives $r_0=c_s\sqrt{3\pi /(32G{\rho }_0)}$. The resulting crossing-time cloud support mass isThis is the same dimensional mass scale as $M_J$, with a different coefficient because the adopted length and time criteria differ. In particular, $t_{\mathrm{ff}}=r_0/c_s$ is not the exact marginal dispersion-relation condition.

With the requested isotropic support approximation, replace $c_s$ by the Alfvén speed $v_A=B/\sqrt{{\mu }_0{\rho }_0}$. In the crossing-time convention the Magnetic Jeans mass isUsing the wavelength-based Jeans mass convention instead replaces $C_{\mathrm{cross}}$ by ${\pi }^{5/2}/6$. The scaling, rather than a convention-independent coefficient, is the shared prediction.

For a fixed cloud mass undergoing isothermal homologous contraction, $\rho /{\rho }_0=(r_0/r{)}^3$ and $c_s$ is constant. Its number of thermal Jeans masses therefore grows asA cloud initially close to the support threshold has $N_{J,0}$ of order one, so smaller pieces can subsequently become independently unstable: this is isothermal fragmentation.

For approximately isotropic contraction with magnetic flux freezing, $B{r}^2$ stays constant, so $B\propto {\rho }^{2/3}$. Hence $M_B\propto B^3/{\rho }^2$ stays constant andThis ideal magnetic support therefore does not create progressively more unstable masses in the way thermal isothermal support does. The same conclusion follows from the conserved mass-to-flux ratio. A magnetically subcritical cloud cannot initiate unrestricted cross-field collapse under ideal flux freezing alone; flux loss or initially supercritical conditions are needed. A real magnetic field also gives anisotropic support: motion and fragmentation along its field lines are not excluded by this scalar-speed estimate.

Hydrostatic equilibrium with self-gravity. A resting fluid obeys $\nabla p=-\rho \nabla \Psi$. Spherical symmetry and the Poisson equation for Newtonian gravity giveEliminate the Newtonian gravitational potential derivative to obtainEquivalently, the shell theorem gives ${\Psi }^{\prime }=GM(r)/r^2$ and $M^{\prime }(r)=4\pi r^2\rho$, which leads to the same hydrostatic equilibrium equation.

A uniform-density liquid. With constant $\rho$, $M(r)=4\pi \rho r^3/3$. There is no central point mass, so the interior gravitational acceleration magnitude is $4\pi G\rho r/3$. ThusUse $p(R)=0$ at a vacuum boundary, neglecting surface tension. The uniform-density planetary pressure profile isThe central value can also be writtenIt has units pascals and decreases smoothly from its maximum at the centre to zero at the surface.

An isothermal gas. Here $c$ is the constant isothermal sound speed, with $p=c^2\rho$. Try $p=A/r^2$, so $\rho =A/(c^2{r}^2)$. The two sides of the equilibrium equation becomeThus $A=c^4/(2\pi G)$, giving the singular isothermal sphere,This does not describe a finite gaseous sphere in vacuum. Both pressure and mass density diverge at the centre, and $p(R)>0$ for every finite outer radius $R$, so it cannot meet $p=0$ there. Extending it to infinity removes the finite-radius pressure jump but makes the total mass infinite. It is useful as an idealized profile over an intermediate radial range, or as a truncated model with nonzero external pressure. The vacuum boundary obstruction for an isothermal sphere also applies to the regular, nonsingular isothermal solution; removing the central singularity does not produce a finite vacuum surface.

Steady outflow. In the wind calculation, take the central mass $M$ to be fixed and neglect the wind's own self-gravity. Constant mass flow rate and the isothermal equation of state giveThe steady radial Euler equation isSubstituting the density derivative givesThis is the isothermal Parker wind equation. A smooth sonic point requires both sides to vanish when $u=c$, soMaintaining constant temperature during expansion generally requires heat input; it is an isothermal model, not the adiabatic blast-wave model of Question 5.

The escaping branch at large radius. Integrate the wind equation, using dimensionless logarithm arguments:Set $x=r/r_s$ and $w=u^2/c^2$. ThenFor the accelerating transonic branch, $w(1)=1$, which fixes $C_0=-3$. At $x\gg 1$ on its supersonic branch, $w\to \infty$ and $\ln w$ is smaller than $w$. ConsequentlyA first correction makes the ordering explicit:This is the large-radius asymptotics of an isothermal Parker wind. The leading expression also holds for other supersonic branches with finite $C_0$.

The source implicitly selects an escaping wind. Positive steady outward flow alone does not determine that branch: a subsonic isothermal breeze also satisfies the differential equation. Its $w\ll 1$ limit gives $w\sim e^{-C_0}{x}^{-4}$, hence $u\propto r^{-2}$ and $\rho \to {\rho }_{\infty }>0$. Such a breeze needs nonzero pressure at infinity. Selecting a regular outflow from a subsonic base into vanishing ambient pressure instead gives the accelerating transonic wind used above.

The asymptotic speed has the same units as $c$; the logarithm is dimensionless. Its unlimited growth is a formal feature of indefinite isothermal heating and a prescribed point-mass potential. Real heating, surrounding pressure, and eventually the neglected wind mass limit that extrapolation.

The first argument mixes a material fluid element with a fixed spatial volume. For a material volume element following the fluid, its volume is generally time dependent:Therefore the material derivative cannot be taken outside the product with $\delta V$ as in the proposed cancellation. Instead, mass conservation statesThe momentum balance for a material fluid element, using Newton's second law, is consequentlyFor a fixed control volume, fluid crosses its boundary and carries momentum; a momentum flux term must be included. Its local conservative form is ${\partial }_t(\rho \mathbf{u})+\nabla \cdot (\rho \mathbf{u}\otimes \mathbf{u})=-\nabla p$, equivalent to the displayed Euler equations after using the continuity equation. The erroneous equation omits the term $\rho \mathbf{u}\nabla \cdot \mathbf{u}$ needed to convert a material derivative of momentum density into this conservative form.

For the second task, take a homogeneous fluid with constant mass density $\rho$, prescribed gravitational acceleration $\mathbf{g}$, incompressible flow and inviscid flow. Its equations areUnder the replacement $\rho \mapsto \tilde{\rho }=C\rho$, the same velocity solves these equations if $p\mapsto \tilde{p}=Cp$, up to a spatially constant pressure offset. This is density scaling invariance of incompressible Euler flow:The comparison assumes the same geometry and compatible velocity or gravitational free surface conditions; an independently specified nonuniform boundary pressure would supply additional information. The gravitational field is taken as prescribed. If the fluid's own self-gravity must be included, its field changes with mass density, so the invariance need not hold. A single freely moving particle's density-independent gravitational acceleration does not justify neglecting that distinction.

Using the supplied first critical value ${\xi }_c\simeq 6.5$, the critical radius of a pressure-confined isothermal sphere isNo radius in metres can be specified without $K$ and a central mass density. Here $c_s$ is the isothermal sound speed, consistent with the fixed-temperature variations in part (d). At a given central mass density, equilibria truncated at $\xi <{\xi }_c$ lie on the first stable branch, while increasing the dimensionless outer radius beyond this first threshold gives an unstable configuration. Numerical integration gives ${\xi }_c\simeq 6.451$, consistent with the paper's rounded $6.5$.

For a transparent comparison with the given Jeans length, set the uniform comparison density ${\rho }_0={\rho }_c$ and use the same speed of sound. ThenBoth scales express the competition between pressure support and self-gravity, with characteristic length $c_s/\sqrt{G\rho }$. Their numerical closeness under this density choice does not mean that the models impose the same threshold: $r_c$ is a boundary radius, whereas ${\lambda }_J$ is the wavelength of a perturbation. Comparing $r_c$ with ${\lambda }_J/2$ would instead give a factor about $2.07$. Using the sphere's mean or boundary density instead of its central density changes the coefficient again, so the reference density must always be stated.

The pressure-confined isothermal sphere gives an actual finite, centrally concentrated hydrostatic equilibrium and an explicit confining pressure, making its radial stability threshold more self-consistent for a roughly spherical isothermal cloud. The uniform Jeans instability model gives a simpler local criterion, includes perturbations at different wavelengths, and does not require constructing a spherical equilibrium. Its ideal uniform static infinite background requires the Jeans swindle or an explicitly specified treatment of background gravity. Neither model includes rotation, magnetic fields, turbulent support or a realistic temperature evolution, so neither is a universal collapse criterion.

differential rotation winds a material arm because $d\Omega /dr\ne 0$. A spiral density wave instead has a pattern speed different from the material speed: stars and gas can pass through its crests. Collective self-gravity, orbital resonances and shear can sustain or repeatedly amplify non-axisymmetric patterns; a visible spiral need not be a single permanent rigid material structure.

For a razor-thin gas disk, local tightly wound disturbances satisfyRotation stabilizes long wavelengths and pressure stabilizes short ones. Minimizing the axisymmetric right-hand side over $|k|$ gives $|k{|}_{\ast }=\pi G\Sigma /c_s^2$ andfor local axisymmetric stability. A collisionless stellar disk instead has $Q_{\ast }={\sigma }_R\kappa /(3.36G{\Sigma }_{\ast })$; the radial random velocity dispersion is not literally a gas sound speed.^[ref]

Local axisymmetric stability does not guarantee global disk stability. Gas and stars gravitate together, thickness modifies the short-scale force, and a disk with $Q>1$ may still have spiral or bar responses. Dissipation in gas can generate shock waves and concentrate star formation near a spiral crest; stellar random motions tend to reduce the collective response.

A dense cloud core first collapses when self-gravity overwhelms its effective support. While radiation escapes readily, compression can remain approximately isothermal and the Jeans mass decreases as density rises, encouraging gravitational fragmentation. Angular momentum usually leads to a disk rather than unrestricted spherical infall.

As the central optical depth increases, cooling becomes slower than compression and pressure produces a first hydrostatic core. molecular hydrogen dissociation near $2000\mathrm{K}$ absorbs compressional energy and lowers the effective adiabatic response, permitting a second collapse of a protostar. After dissociation the gas stiffens and a much more compact stellar core forms. An accretion shock and surrounding disk then convert infall energy into radiation. A small central core forms before most of the envelope arrives, as demonstrated in radiation-hydrodynamic collapse calculations.^[ref]

The isothermal inside-out collapse of a singular isothermal sphere model provides one idealized supply law, $\dot{M}\sim a_T^3/G$, but rotation, magnetic fields and an unsteady reservoir modify real accretion. Its luminosity scale is $L_{\mathrm{acc}}\sim GM\dot{M}/r_{\ast }$. Dust can absorb the emerging short-wavelength radiation and re-emit it in the infrared, so an embedded source need not have a directly visible stellar photosphere.

As envelope accretion declines, gravitational contraction becomes important on the Kelvin-Helmholtz cooling time $t_{\mathrm{KH}}\sim G{M}^2/(r_{\ast }L)$. deuterium-burning energy-generation rate can temporarily influence the radius; sustained core hydrogen burning eventually establishes the main sequence if the final mass is sufficient. A low-mass convective object approaches along a Hayashi-like track, while more massive objects need not share that track. accretion, contraction and nuclear burning are distinct energy sources, and a protostar is not generally a fixed-mass object.

Radial epicyclic motion. At fixed specific angular momentum, use ${\Phi }_{\mathrm{eff}}=\Phi +{\ell }^2/(2r^2)$ and $\ddot{r}=-{\Phi }_{\mathrm{eff}}^{\prime }$. A circular orbit at $r_c$ obeys ${\Phi }^{\prime }(r_c)={\ell }^2/r_c^3=r_c{\Omega }^2$. Expanding $r=r_c+\xi$ gives $\ddot{\xi }+{\kappa }^2\xi =0$, whereFor ${\kappa }^2>0$ the radial epicyclic frequency describes simple harmonic motion. The angular correction follows from conservation of angular momentum, $\dot{\varphi }\simeq \Omega (1-2\xi /r_c)$, and describes an epicyclic motion around the guiding circular orbit.

Stabilizing terms. In the given local fluid dispersion relation, ${\kappa }^2$ is rotational support, important at long wavelength; $a_0^2{k}^2$ is random-motion or pressure support, important at short wavelength. The term $-2\pi G\Sigma |k|$ is destabilizing self-gravity. Although the wording refers to stellar dispersion, this algebra uses the supplied fluid approximation; a collisionless stellar treatment has a reduction factor and a different numerical normalization.^[ref]

Axisymmetric threshold. For $n=0$, set $q=|k|\ge 0$. Completing the square givesThus $Q>1$ gives local axisymmetric oscillatory stability; $Q<1$ gives a growing band, withAt $Q=1$, the minimizing mode has zero frequency and is marginal, not a strictly restoring oscillation. Local axisymmetric stability does not establish stability to every global nonaxisymmetric mode.

long-wave and short-wave branches at marginal Toomre stability. For $Q=1$ and $n\ge 1$,The long-wavelength branch, using the minus sign, requiresThese endpoints are the inner and outer Lindblad resonance conditions. The stated dispersion relation does not restrict every wave solution to that interval. Its plus-sign short-wavelength branch has $q>0$ even outside it: for example, ${\Omega }_P=\Omega +2\kappa /n$ gives $q=3\kappa /a_0$ and satisfies the equation. Additional branch selection or resonant absorption physics would be needed for an absolute confinement claim. No such condition is supplied.

Toomre stability radius of a flat-rotation disk. With $\Omega =v_0/r$, $\kappa =\sqrt{2}{v}_0/r$. Using $G=4.3009\times 10^{-3}\mathrm{pc}(\mathrm{km}{\mathrm{s}}^{-1}{)}^2{M}_{\odot }^{-1}$,The supplied local model therefore predicts stable radii below about $4.19\mathrm{kpc}$, marginality at that radius, and axisymmetric instability farther out while its assumed constant surface density and thin-disk approximation remain applicable. An infinite constant-density, flat-rotation disk is an idealization, not a global self-consistent galaxy model.

Small primordial fluctuations seed later structure. For adiabatic cosmological perturbations on scales outside the horizon, the uniform-density curvature perturbation is conserved when the non-adiabatic pressure perturbation and the decaying mode can be neglected. A density contrast outside the horizon is gauge-dependent, so its apparent power-law growth should not be confused with directly measured local collapse.

For subhorizon, pressureless matter, the linearized continuity equation, Euler equation and Poisson equation combine intoHubble expansion damps peculiar motion; self-gravity amplifies overdensities. During radiation domination, radiation pressure drives acoustic oscillations and cold dark matter growth after horizon entry is only logarithmic. During matter domination the growing and decaying solutions are $\delta \propto a$ and $\delta \propto a^{-3/2}$. Dark-energy domination later suppresses growth relative to the matter-era law.

For a pressure-supported fluid the Fourier-mode equation also contains $+c_s^2{k}^2\delta /a^2$, giving a Jeans length separating acoustic support from gravitational growth. Before recombination, baryons are tightly coupled to photons and oscillate rather than freely following dark matter. Silk damping from photon diffusion damps small-scale photon-baryon fluctuations; after decoupling, baryons can fall into dark-matter potential wells. free streaming smooths species-dependent scales and distinguishes cold from warm or hot dark matter.

Schematically, the linear matter power spectrum can be written $P_{\delta }(k,a)=A{k}^{n_s}{T}^2(k)D^2(a)$ under the usual primordial-curvature normalization. The cosmological transfer function records horizon entry, matter-radiation equality, acoustic physics and damping; the linear growth factor records later amplification. density contrasts approaching unity require nonlinear collapse and halo formation beyond this approximation. cosmic microwave background anisotropies and galaxy clustering test the same initial perturbations through different transfer histories; galaxy bias, redshift-space distortions and observational selection enter the galaxy comparison.

With $C_{\mathbf{k}}$ and $a$ from e, inversion givesFor $t>0$, deform the contour downwards. Closing below is clockwise, so the residue theorem multiplies the sum of enclosed residues by $-i$. The forcing pole has residue $i{C}_{\mathbf{k}}{e}^{-iat}/{ϵ}_{\mathbf{k}}(a+i{0}^{+})$ and therefore contributesA simple collective root ${\omega }_j$ of the continued gravitational dielectric function contributesThese are the Landau mode transients excited by introducing the mass abruptly. For a strictly damped background they decay, together with any remaining phase-mixing contour contribution under the usual regularity and decay assumptions. ThusThe causal boundary value is what the paper abbreviates by ${ϵ}_{\mathbf{k}}(\mathbf{k}\cdot {\mathbf{v}}_p)$. Physically, the moving mass polarizes the stellar population and carries a gravitational wake. Each spatial harmonic is changed in amplitude and phase by $1/{ϵ}_{\mathbf{k}}$ at the forcing frequency. Self-gravity can amplify the wake. The induced part alone is $(1/{ϵ}_{\mathbf{k}}-1)\delta {\Phi }_{\mathbf{k}}^{\mathrm{ext}}$; the wake's back-reaction on the perturber produces dynamical friction.

The waiting time is controlled by the slowest excited collective damping, not just an orbital or crossing time. If ${\omega }_j={\omega }_{rj}-i{\Gamma }_j$ with ${\Gamma }_j>0$, then a modal amplitude is suppressed by $e^{-{\Gamma }_{j}t}$, and one needs $t{\Gamma }_j\gg 1$ for every appreciably excited mode, together with decay of the phase-mixing remainder. Large transient prefactors can require additional time. Near marginal stability, a relevant ${\Gamma }_j$ tends to zero and the required waiting time diverges; the stationary wake approximation is then not uniform in the stability parameter.^[ref]

Qualification of stability. If stable means only non-growing, a persistent neutral mode invalidates the asserted asymptote. A concrete example is a homogeneous waterbag distribution with parallel marginal $F(u)=\rho 1_{|u|<V}/(2V)$. Its gravitational dielectric function isFor $k^2{V}^2>4\pi G\rho$, the roots $\omega =\pm \sqrt{k^2{V}^2-4\pi G\rho }$ are real, so the switch-on transients persist although that harmonic has no growing mode. Taking the smallest box wavenumber to satisfy the same inequality makes the whole waterbag model free of growing harmonics. The claimed late-time formula therefore requires damping of all excited modes and ${ϵ}_k(a)\ne 0$; exact forcing at an undamped root instead produces a secular response.

The first equation is vertical hydrostatic equilibrium: the upward pressure gradient balances the central mass's vertical gravitational force. For $|z|\ll r$, the point-mass acceleration is $-GMz/(r^2+z^2{)}^{3/2}\simeq -{\Omega }_K^{2}z$, where Keplerian rotation gives ${\Omega }_K^2=GM/r^3$. The disc's self-gravity and vertical inertia are neglected.

The second equation is local thermal balance. The divergence of the vertical radiative flux equals viscous dissipation. In the convention of this alpha disk, $\nu =\alpha p/(\rho \Omega )$, and the dominant Keplerian shear has $rd\Omega /dr=-3\Omega /2$. ThereforeThe radiative flux increases from zero at the symmetric midplane to a positive outgoing value at the upper surface. Radial heat transport, thermal advection and rapid storage of thermal energy are excluded from this local balance.

The third equation is radiative diffusion, equivalentlyIt applies to optically thick material, with $\kappa$ the opacity per unit mass and $\sigma$ the Stefan-Boltzmann constant. Here constant electron-scattering opacity represents Thomson scattering; the formula uses a diffusive, nearly thermal radiation field. Positive upward flux requires temperature to decrease upwards.

The final equation is the ideal gas pressure law, with mean molecular weight ${\mu }_m$ in proton-mass units. Writing $\mathcal{R}=k/({\mu }_m{m}_p)$ gives $p=\mathcal{R}\rho T$. The pressure is gas pressure because radiation pressure is assumed negligible. The four equations express vertical support, local heating balance, diffusive cooling and the gas equation of state.

Estimate characteristic bulk terms. Put $ϵ=H/r\ll 1$ and assume smooth radial variation on scale $r$, secular evolution, and $\alpha \lesssim 1$. Vertical support gives $c_s^2\sim p/\rho \sim {\Omega }^2{H}^2$, while the alpha disk prescription and slow accretion implyThese velocity estimates follow from the viscous radial evolution scale and mass conservation. They are estimates for a slowly evolving thin disc, not a claim that geometric thinness excludes independently forced rapid vertical motion.

For vertical balance, the exact central force expands asAt $z\sim H$, the neglected correction has magnitude $\rho {\Omega }^{2}H{ϵ}^2$, compared with the retained force $\rho {\Omega }^{2}H$. Radial pressure support similarly changes the rotation from its Keplerian value byVertical meridional inertia scales as $u_r{u}_z/r$ or $u_z^2/H$, both $O({\alpha }^2{ϵ}^4{\Omega }^{2}H)$. Meridional viscous accelerations can scale as $\nu u_z/H^2$ or $\nu u_r/(rH)$, both $O({\alpha }^2{ϵ}^2{\Omega }^{2}H)$. Thus the gravity correction and possible viscous correction giveSome omitted terms are smaller than this leading order.

For the heating equation, the retained dissipation is $q_0^{+}\sim \rho \nu {\Omega }^2\sim \alpha p\Omega$. Radial pressure support and the finite-height correction change the dominant shear by a relative $O({ϵ}^2)$. Smooth vertical variation of the rotation gives ${\partial }_z(r\Omega )=O(ϵ\Omega )$, whose additional squared shear contributes $O({ϵ}^2{q}_0^{+})$. Meridional shear of order $u_r/H$ gives an additional relative $O({\alpha }^2{ϵ}^2)$ contribution.

Thermal advection and compressional work are of order $\rho u_r{c}_s^2/r$ and $p{u}_r/r$, respectively. Relative to $\alpha p\Omega$ they are $O({ϵ}^2)$. Thermal storage on the viscous evolution timescale $t_{\nu }\sim r^2/\nu$ has the same order, since $p/t_{\nu }\sim \alpha p\Omega {ϵ}^2$. For radiative diffusion with the same local diffusivity in different directions, smooth gradients give $F_r/F_z\sim H/r$, soConsequentlyThis is the thin-accretion-disk truncation error for the secular, radially smooth approximation. Neglect of disc self-gravity and radiation pressure, optical thickness, the opacity model and the idealized zero surface boundary have their own validity conditions; their errors do not become small solely because $H/r$ does. Near a sharp radial boundary layer or during thermal evolution on the heating timescale, the above estimates need not apply.

Consider the nontrivial gravitational mode with $k\ne 0$. At ${\omega }^2=0$, horizontal momentum requires $ikW=0$, hence $W=0$ andSubstituting into the Poisson equation for Newtonian gravity and using $4\pi G{\rho }_0/c_s^2=2/H^2$ givesSet $\tau =\tanh (z/H)$ and $\nu =kH$. Since ${\partial }_z=(1-{\tau }^2){\partial }_{\tau }/H$, this becomes the associated Legendre equation of degree one and order $\nu$:

To verify the two explicit solutions, put $\eta =z/H$ and $T=\tanh \eta$. They arebecause $e^{2\eta }=(1+\tau )/(1-\tau )$. Differentiating $F_{+}$ twice giveswhich cancels identically against $[2(1-T^2)-{\nu }^2]F_{+}$. Reflection $\eta \mapsto -\eta$ proves the same for $F_{-}$. Their Wronskian is $2\nu (1-{\nu }^2)$, evaluated at $\eta =0$; thus they are generally independent but coincide or become dependent at the exceptional orders $0,\pm 1$.

Use $n=|k|H>0$ to analyze decay, which is unaffected by the sign of $k$. The solution decaying at $+\infty$ is proportional to $e^{-n\eta }(n+\tanh \eta )$. At $-\infty$ its leading term is $(n-1)e^{n|\eta |}$, so it also decays there only if $n=1$. At this valueTherefore the marginal fragmentation mode of an isothermal slab hasAt $n=1$ the second independent solution obtained by reduction of order grows at infinity, so the coincidence of the two explicit formulas does not supply an extra admissible mode. At $k=0$ a uniform vertical translation gives another neutral disturbance, and mass-preserving relabeling displacements can have $\delta \rho =\varphi =0$. Neither is the nonzero-wavenumber fragmentation threshold just calculated.

The unstable fragmentation modes have wavelengths greater than $2\pi H$. Gas-pressure support becomes less effective relative to self-gravity at long wavelength. The side of the threshold can also be checked from the energy identity: at $|k|=1/H$ the marginal mass density disturbance has zero quadratic energy. Keep that mass density disturbance and decrease $|k|$; the pressure integral stays fixed, while the gravitational term becomes more negative because the positive inverse operator $(-{\partial }_z^2+k^2{)}^{-1}$ increases as $k^2$ decreases. Thus the quadratic energy becomes negative, giving an unstable mode on the longer-wavelength side. This argument concerns the gravitational sector and its stated decay conditions.

Assume a nonrotating, nonmagnetic equilibrium, inviscid adiabatic perturbations, and no dissipation. Let $g={\varphi }^{\prime }(R)>0$, so hydrostatic equilibrium is $p^{\prime }=-\rho g$. Use $\mathbf{\xi }{e}^{i\omega t}$ for the fluid displacement. To avoid a source ambiguity, write Eulerian fluid perturbations as $p_1,{\rho }_1,{\varphi }_1$, and use $\Delta$ for quantities following the displaced material element, as in Eulerian and Lagrangian fluid perturbations. They are related by $\Delta p=p_1+{\xi }_R{p}^{\prime }$ and similarly for mass density and potential.

The printed identity needs two corrections: its compression denominator must be $\gamma p$, not $\gamma \rho$, and the pressure and potential appearing in this form of the identity must be Eulerian fluid perturbations. The printed denominator has the wrong dimensions for an energy density. The claimed Lagrangian interpretation does not produce the displayed buoyancy-compression split; the following derivation specifies the consistent version.

Linearized mass conservation, the adiabatic relation, and the Poisson equation giveConsequently $p_1=-\gamma p\nabla \cdot \mathbf{\xi }-{\xi }_R{p}^{\prime }$. Linearizing the Euler equations for an inviscid fluid about the static star gives the linear adiabatic stellar oscillation equationsRequire regularity at the center, a vanishing Lagrangian surface pressure, and the usual matching to a decaying exterior gravitational potential. Because $p$ and $\rho$ vanish at the surface, the pressure boundary term vanishes for regular displacements; compactly supported trials also suffice below. The perturbed Newtonian gravitational potential solves Laplace's equation outside the star and must be included there when its field energy is computed.

Multiply the momentum equation by $-{\mathbf{\xi }}^{\ast }$ and integrate. The pressure term is $-{\int }_V{p}_1\nabla \cdot {\mathbf{\xi }}^{\ast }dV$. Let $d=\nabla \cdot \mathbf{\xi }$. The adiabatic relation and hydrostatic equilibrium implyThe pressure and equilibrium-gravity terms then combine pointwise aswhere $N^2=g[(1/\gamma )p^{\prime }/p-{\rho }^{\prime }/\rho ]$ is the stellar buoyancy frequency. The potentially troublesome cross terms cancel exactly.

For the perturbed self-gravity term, integration by parts givesThus the corrected stellar displacement energy identity isThe right side is a real quadratic form. With these conservative boundary conditions, the force operator is a self-adjoint operator. Its Rayleigh quotient characterizes squared mode frequencies, and a negative admissible trial quotient proves an unstable spectral component; it does not require the trial displacement itself to be an exact eigenfunction.

For a generic compressive disturbance of scale $\ell$ in a region with slowly varying equilibrium coefficients, estimate ${\rho }_1\sim \rho \xi /\ell$, $p_1\sim \gamma p\xi /\ell$ and ${\varphi }_1\sim 4\pi G{\rho }_1{\ell }^2$. The self-gravity and compression energy densities have orderswith adiabatic sound speed $c_s^2=\gamma p/\rho$. Under the plausible stellar estimate $c_s^2\sim G\bar{\rho }{a}^2$, global modes can have comparable gravity and pressure contributions, while their ratio for small-scale compressive modes is of order $(\rho /\bar{\rho })(\ell /a{)}^2$. Thus perturbations with $\ell \ll c_s/\sqrt{4\pi G\rho }$ admit the Cowling approximation. This estimate compares two nonzero compressive terms; it cannot alone justify neglecting gravity for a trial deliberately chosen to have $p_1=0$.

For an explicit convective trial, choose $c<d<e<b$ and a nonzero smooth bumpLet $Y$ be a normalized spherical harmonic of degree $\ell \ge 1$, let $L=\ell (\ell +1)$, and let ${\nabla }_{\Omega }$ be the unit-sphere gradient. The pressure-free localized stellar convective trial isIt is supported away from the center and surface, so it satisfies every needed mechanical boundary condition. Since ${\Delta }_{\Omega }Y=-LY$,Therefore $p_1=-\gamma p\nabla \cdot \mathbf{\xi }+\rho g{\xi }_R=0$ exactly. Its buoyancy energy is ${\int }_d^e\rho N^2{R}^2|f{|}^{2}dR<0$, whereas its kinetic norm ${\int }_d^e\rho R^2(|f{|}^2+L|h{|}^2)dR$ is positive. In the Cowling approximation this gives a negative Rayleigh quotient, and hence $N^2<0$ in an interior shell implies convective instability. Retaining perturbed self-gravity only adds a nonpositive term, so the same trial proves instability in the full conservative problem too.