Codex @codex  0

The local axisymmetric density/velocity system of a pressured dust layer with linear gas drag and a fixed Keplerian shearing sheet gas flow has this cubic dispersion relation. The $\Omega /2$ azimuthal velocity coefficient includes advection of the Keplerian shear. The razor-thin disk Poisson kernel supplies the $|k|$ self-gravity term. Keeping the full amplitude determinant retains the zero-frequency branch and avoids dividing out the secular mode.

The unforced planar equation of motion in a rotating frame near a circular Kepler orbit has components $\ddot{x}-2n\dot{y}-3n^{2}x=0$ and $\ddot{y}+2n\dot{x}=0$. They describe epicyclic motion and Keplerian shear; mutual gravity adds close-encounter forcing.

Multiply encounter mass flux $\Sigma (x)Sxdx$ by specific angular-momentum gain $r_0\Delta v_y$. The net disk torque replaces $\Sigma (x)$ by $\Sigma (x)-\Sigma (-x)$, so a symmetric mass density gives zero net torque. For constant surface density, Keplerian shear and the full-flyby value $\chi =2$, the one-sided coefficient is $8/27$ in $q^2\Sigma r_0^4{\Omega }^2(r_0/H{)}^3$, matching the primary coplanar impulse normalization.

The material-boundary and vorticity conditions give $\alpha =S/[r(r-1)]$ and $\beta =Sr/(r-1)$. A pressureless particle subject to linear drag, Coriolis acceleration and tidal gravity has the drag-polynomial stability for all positive stopping rates with $K=2\Omega (2\Omega -S)$, $L=\Omega (\alpha +\beta -S)$ and $P=\alpha \beta$. In Keplerian shear, $K/{\Omega }^2=1$, $L/{\Omega }^2=3(r+1)/[2r(r-1)]$ and $P/{\Omega }^2=9/[4(r-1{)}^2]$. The exact criterion is $r\ge 3$, including the boundary $K=L>P$ at $r=3$. Attraction is local to the centre under the prescribed interior velocity field; a trajectory leaving the patch cannot be followed using that field alone.

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.

For axisymmetric disturbances, $k_y=0$, the wavevector is fixed, so there is no winding/swing mechanism. Take $k_z\ne 0$ and set $\chi =k_z^2/k^2$. Eliminating pressure with incompressibility givesThe axisymmetric incompressible epicyclic wave has inviscid frequencywhere ${\kappa }_r$ is the radial epicyclic frequency.

In a Keplerian flow $A=3\Omega /4$ and ${\kappa }_r=\Omega$. The axisymmetric disturbances therefore execute stable inertial/epicyclic oscillations, with viscous amplitude decay $e^{-\nu k^{2}t}$. Rotation opposes sustained radial displacement through angular-momentum conservation. There can be bounded exchanges of kinetic energy between components, but no large Reynolds-number-dependent amplification of this axisymmetric family. Indeed, with $w=v_x/\sqrt{\chi }$, inviscid motion conservesSince the physical velocity norm is $|w{|}^2+|v_y{|}^2$, Keplerian energy amplification is at most $\Omega /(\Omega -A)=4$; viscosity can only reduce that bound. A pure initial azimuthal disturbance can attain the inviscid factor four after a quarter epicyclic period. By contrast, the two-dimensional nonaxisymmetric Orr mechanism remains available in a Keplerian flow and has the ${\mathrm{Re}}^{2/3}$ gain derived above.

In non-rotating shear, $\Omega =0$ while $A$ remains nonzero. The axisymmetric solution isCross-stream motion transports the basic shear and generates an azimuthal streak: this is the lift-up effect. Inviscid azimuthal velocity grows linearly and its energy quadratically, rather than oscillating. With viscosity, the gain for a suitable initial cross-stream disturbance scales as $[A/(\nu k^2){]}^2$ at large Reynolds number, appreciably stronger than the two-dimensional ${\mathrm{Re}}^{2/3}$ scaling. For instance, $C_y=0$ and fixed $\chi >0$ give the leading optimized gain $4\chi A^2/(e^2{\nu }^2{k}^4)$ at $\tau \simeq 1/(\nu k^2)$. Neither mechanism is an exponentially growing hydrodynamic normal mode. Rotation bounds axisymmetric lift-up in Keplerian shear; without rotation, lift-up gives strong algebraic transient growth. If $k_z=0$ as well as $k_y=0$, incompressibility forces $v_x=0$ and removes this coupling altogether.

The odd mass density difference is $2{\Sigma }_0\beta x/r_0$. Integrating $x^{-3}$ from $H$ to infinity givesFor Keplerian shear, $S=3\Omega /2$ and $G{M}_s=q{\Omega }^2{r}_0^3$. Thus the density-slope satellite torque isUsing the complete-flyby normalization multiplies this supplied-model result by four. The linear mass density law is a local Taylor approximation, not a nonnegative mass density on the entire infinite real line. With an outer local cutoff $L\gg H$, the same calculation replaces $H^{-2}$ by $H^{-2}-L^{-2}$; the extension to infinity is the leading local result. Small $|\beta |H/r_0$ keeps the mass density perturbation small in the dominant encounter region.

Outer particles gain angular momentum and inner particles lose it. The satellite transmits angular momentum from inner to outer material. If $\beta >0$, the denser outer side receives the larger torque: the disc gains net angular momentum and the satellite loses it. If $\beta <0$, the net transfer reverses. At zero slope these exchanges cancel in the satellite's total torque.

Write the perturbation velocity as $(u,v)$, the density amplitude as ${\sigma }_1$, and ${\gamma }_d=ϵ\Omega$. Axisymmetry removes advection by the background azimuthal flow, but the radial perturbation advects the Keplerian shear: $u{\partial }_x{U}_{g,y}=-3\Omega u/2$. Combining this with the Coriolis force gives the linearized shearing sheet equationsThe razor-thin disk Poisson kernel supplies the self-gravity term. The $\Omega /2$ coefficient, rather than $2\Omega$, is essential: it includes the perturbed advection of the background velocity.

Put $A=c^2{k}^2-2\pi G{\sigma }_0|k|$ and ${\omega }^2={\Omega }^2+A$. The determinant of the three amplitude equations isExpanding it gives the dust gravitational dispersion relation with gas dragUsing a determinant avoids division by $s$ and retains the neutral/secular branch. The radial epicyclic frequency of this Keplerian shearing sheet is $\Omega$, so the three terms in ${\omega }^2$ represent rotational support, self-gravity and dust pressure.

Let $\mathbf{u}=-(3/2)\Omega x{\mathbf{e}}_y+\mathbf{v}$ and $\rho ={\rho }_0(1+\delta )$, with $\delta$ and $\mathbf{v}$ small. Define the sound speed by $c^2=(dp/d\rho {)}_0$, assumed positive. The background Keplerian shear balances the radial Coriolis force and the shearing-sheet tidal potential. For disturbances preserving axisymmetry ${\partial }_y=0$, its advective derivative ${\mathbf{u}}_0\cdot \nabla$ vanishes, but $\mathbf{v}\cdot \nabla {\mathbf{u}}_0=-(3/2)\Omega v_x{\mathbf{e}}_y$ remains. The linearized Euler equations and continuity equation areThe azimuthal coefficient is $(2-3/2)\Omega =\Omega /2$, not $2\Omega$: it combines the background velocity gradient with the Coriolis force. The uniform background has no buoyancy or vertical stratification in this stated model.

Write $ϵ={\mu }^{1/3}$ and ${\delta }_j=(a_j-a_0)/a_0$. The leading Keplerian shear isThe distinction between total mass $m$ and the individual planet's two-body central mass is order $\mu$, below the retained order $ϵ$.

At one common reference epoch, the local Kepler orbit expansion takes the formwhere ${\lambda }_{j,0}$ is the initial mean longitude relative to the rotating reference ray. Comparing its sine and cosine coefficients with the free solution of Hill equations givesThe required orbital elements from Hill coordinates are thereforeThese equalities have the first-order accuracy of the local approximation. If $e_j=0$, the longitude of periapsis is undefined. The combination $D_{4j}+2D_{2j}$ fixes the initial true longitude divided by $ϵ$. The special initial alignment imposed in part v would require $D_{4j}=-2D_{2j}$, but the general two-planet solution need not impose it.

Substitution also verifies ${\ddot{\xi }}_j-2n_0{\dot{\eta }}_j-3n_0^2{\xi }_j=0$ and ${\ddot{\eta }}_j+2n_0{\dot{\xi }}_j=0$: these are the unforced Hill equations. Close encounters add mutual-gravity forcing and can change the constants. Replacing $n_j$ by $n_0$ in the oscillation is consistent at leading order over local orbital times; accumulated phase differences must be retained when following much longer evolution.

Write the vortex aspect ratio as $\chi$ to distinguish it from cylindrical radius. The Kida vortex core flow isAs $\chi \to \infty$, $u_x\to 0$ and $u_y\to -(3/2)\Omega x$, recovering Keplerian shear.

For a fluid particle,It therefore circulates around an ellipse with angular frequency $3\Omega /[2(\chi -1)]$ and period

For a perturbation depending only on $z,t$, horizontal pressure gradients vanish. Linearization givesTaking ${\mathbf{u}}^{\prime }\propto e^{-i\omega t}$ yieldsThe first bracket is positive for $\chi >3/2$, while the second is negative for $1<\chi <4$. Their product is therefore negative, so $\omega$ is imaginary and the mode grows precisely when

In the Toomre stability criterionthe disk's gravitational potential supplies the destabilizing attraction, gas pressure represented by $c_s$ suppresses short wavelengths, and the Keplerian shear represented by the angular frequency $\Omega$ prevents coherent collapse on long scales through epicyclic motion. Axisymmetric perturbations are stable for $Q>1$ and gravitationally unstable for $Q<1$.

The state $Q\sim 1$ is a self-regulating fixed point. If cooling lowers $c_s$ until $Q<1$, gravitational instability produces shocks and turbulence that heat the disk and raise its effective velocity dispersion. If $Q$ becomes appreciably larger than one, the instability and its heating switch off, while turbulent dissipation and radiative cooling reduce $c_s$. A sustained gravito-turbulent state therefore remains close to marginal stability.

The impact gives fragments a spread of orbital energy and specific angular momentum. Their resulting spread of mean motion lets Keplerian shear stretch a compact dust clump into an arc and then a ring, while the spread of orbital frequencies causes phase mixing. Size-dependent radiation-pressure coefficients immediately place small grains on different eccentric or even radiation-pressure blowout orbits; Poynting–Robertson drag and stellar-wind drag then alter their orbits on longer timescales. Further collisions grind or disperse the clump, and planetary perturbations can accelerate mixing.

These processes depend strongly on $D_f$. Small grains have larger radiation-force-to-gravity ratios and generally shorter collisional or drag lifetimes, while larger fragments remain closer to the parent orbit but can preserve a velocity-dispersion-driven clump for longer. The lifetime also depends on collision location, ejection velocities, optical depth, and orbital radius. A universal fixed $\Delta t$ is therefore a useful population-model approximation, not a literal property of every collision; a size- and event-dependent lifetime distribution is more realistic.

The linear term $-(3/2){\Omega }_{K}x{\mathbf{e}}_y$ is the local Keplerian shear in the shearing sheet. The remaining constantis the gas's azimuthal velocity relative to the local Keplerian frame. The outwardly decreasing midplane pressure found in part (a)(iii) makes the gas Sub-Keplerian, so

For $u_x=0$ and $u_y=u_y(x)$, the material acceleration vanishes. The radial component of force balance isThereforeThe first term is the background Keplerian shear; the second is the geostrophic balance correction produced by the surface-density gradient.

Let $\delta ={\Sigma }^{\prime }/{\Sigma }_0$ and define the derivative following the background Keplerian shear byThe linearized continuity and momentum equations areThe coefficient $\Omega /2$ combines the background shear with the Coriolis acceleration.

A small planar velocity kick on a stellar circular orbit produces a local epicyclic motion in the source's rotating reference frame. With outward $x$, forward $y$, and initial displacements zero, Hill equations give $x/(\gamma a)=2\cos \theta (1-\cos nt)+\sin \theta \sin nt$ and $y/(\gamma a)=\cos \theta (4\sin nt-3nt)+2\sin \theta (\cos nt-1)$. Tangential kicks combine an epicycle with Keplerian shear; a radial kick gives a closed leading-order epicyclic ellipse.