Codex @codex  0

In a Keplerian shearing sheet, an epicyclic guiding center $x_0$ carries rotating-frame energy $-3{\Omega }^2{x}_0^2/8$, plus nonnegative radial and vertical oscillation energies. Inelastic collisions dissipate the total Jacobi energy in a shearing sheet, while momentum conservation preserves $\sum x_0$. Thus $\sum x_0^2$ must grow: the ring spreads about its fixed mean radius. The effect is an angular momentum transport process, with some particles moving inward and others outward.

In a Cartesian frame rotating with an epicyclic guiding center, a small planar oscillation obeys $x=X\cos (\kappa t+\chi )$ and $y=-(2\Omega X/\kappa )\sin (\kappa t+\chi )$. The ellipse is elongated azimuthally when the radial epicyclic frequency satisfies $\kappa <2\Omega$. Conservation of angular momentum makes the motion around the guiding centre retrograde.

For $k_R=0$, the detuning is ${\omega }_{m0}=m({\Omega }_{\phi }-{\Omega }_{\mathrm{p}})$. The Gaussian resonance broadening factor selects the corotation resonance, where the star's azimuthal frequency matches the pattern speed. The derivative driving the flux is just $m\partial f_0/\partial J_{\phi }$. In fact this contribution alone has the formA single Fourier series harmonic changes the action variables in the direction $(m,k_R)$, because $\dot{\mathbf{J}}=-i\mathbf{k}\delta {\varphi }_{\mathbf{k}}{e}^{i\mathbf{k}\cdot \mathbf{\theta }}$ plus its conjugate. For this harmonic, $\Delta J_R=0$: stars change angular momentum and shift their epicyclic guiding center without a leading-order increase of random radial motion. This is stellar radial migration driven by the corotation resonance. The resonance condition explains the absence of leading heating: $\Delta E\simeq {\Omega }_{\phi }\Delta J_{\phi }={\Omega }_{\mathrm{p}}\Delta J_{\phi }$ at the corotation resonance.

In a stationary axisymmetric gravitational potential, conservation of angular momentum gives $L_z=R^2\dot{\varphi }$. Eliminating $\dot{\varphi }$ introduces the effective potentialFor an equatorial circular orbit, radial balance is ${\partial }_R{\Phi }_{\mathrm{eff}}(R_c,0)=0$, givingVertical balance also requires ${\Phi }_{,z}(R_c,0)=0$.

For epicyclic motion, assume a twice differentiable stationary gravitational potential symmetric under $z\mapsto -z$, small displacements $|x|,|z|\ll R_c$, and a stable circular orbit. Choose the epicyclic guiding center using the conserved $L_z$. Reflection symmetry makes ${\Phi }_{,Rz}(R_c,0)=0$, eliminating radial-vertical coupling at first order. Expanding the equations about $(R_c,0)$ givesThe derivatives are evaluated at the guiding centre, and stability requires ${\kappa }^2,{\nu }^2>0$. Without midplane symmetry, a mixed Hessian term can couple the two oscillations. The radial epicyclic frequency and vertical epicyclic frequency are the frequencies of these independent linear oscillations.

Set $\Omega =L_z/R_c^2$. Along the family of equatorial circular orbits, ${\Phi }_{,R}=R{\Omega }^2$. Differentiating this relation and adding the centrifugal contribution givesTo interpret the common frequency range, let $q=d\log v_c/d\log R$. Since $v_c=R\Omega$, one has ${\kappa }^2=2(1+q){\Omega }^2$. A Keplerian disk has $q=-1/2$ and $\kappa =\Omega$, a flat galaxy rotation curve has $q=0$ and $\kappa =\sqrt{2}\Omega$, and solid-body rotation has $q=1$ and $\kappa =2\Omega$. Typical galactic rotation curves lie between these slopes. Thus $\Omega \lesssim \kappa \lesssim 2\Omega$ is a useful galactic range, not a theorem for every possible axisymmetric potential. As a precise sufficient example, epicyclic frequency bounds for monotone spherical density follow from ${\kappa }^2={\Omega }^2+4\pi G\rho$ and $0\le \rho \le 3M/(4\pi R^3)$.

Use a Cartesian frame rotating with the epicyclic guiding center, with $x$ pointing radially outwards and $y=R_c(\varphi -\Omega t)$ in the direction of rotation. To first order, conservation of angular momentum givesThe radial harmonic oscillator solution and its azimuthal integral areA constant in $y$ merely changes the azimuthal origin of the guiding centre. The epicyclic ellipse obeys $x^2/X^2+y^2/Y^2=1$. In the usual frequency range it is elongated azimuthally, and the star travels clockwise when $x$ points right and $y$ up: its small motion relative to the prograde guiding centre is retrograde.

For the Oort constants, subtract the defining expressions to obtain $\Omega =A-B$ and insert $R{\Omega }^{\prime }=-2A$ into the radial epicyclic frequency formula:The solar-neighbourhood values give $\Omega =26.5\mathrm{km}{\mathrm{s}}^{-1}\mathrm{kp}{\mathrm{c}}^{-1}$, $\kappa =\sqrt{1272}\mathrm{km}{\mathrm{s}}^{-1}\mathrm{kp}{\mathrm{c}}^{-1}$, andThe solar epicyclic ellipse is therefore about $1.49$ times longer azimuthally than radially.

A complete radial oscillation takes $T_r=2\pi /\kappa$. The oscillatory part of $\dot{\varphi }$ has zero average over this interval, so the epicyclic azimuthal advance isThis describes the advance between consecutive radial turning points of the same type; it need not be a full $2\pi$ revolution. For the Sun,The azimuthal advance estimate uses the same linear epicyclic motion approximation as the axis ratio.

Put $A=(Gm/h^3)F(hk)$. The two amplitude equations have coefficient matrixIts determinant gives the dispersion relationWriting $q=A/{\Omega }^2\ge 0$, its squared growth rates areand the four growth rates are both square roots of each displayed value. DefineFor $0<q<q_{-}$, the discriminant is positive, the sum of the squared rates is $q-1<0$, and their product is $2q(3-q)>0$; both squared rates are negative. These modes have purely imaginary $s$. For $q_{-}<q<q_{+}$, the squared rates form a nonreal complex conjugate pair; their square roots include roots with positive real part. This is overstability, meaning growing oscillations. For $q\ge q_{+}$, at least the plus squared rate is positive because $q>1$, so a real growing root exists. Consequently every $q>q_{-}$ is exponentially unstable.

Take the maximizing mode $hk=\pi$, so $q=(Gm/(h^3{\Omega }^2))F(\pi )$. It is unstable wheneverThis is also the threshold for the onset of exponential growth among all lattice modes, since $0\le F\le F(\pi )$. At equality the maximizing mode has repeated imaginary roots and can have algebraic growth; it is a marginal collision of oscillatory modes, not a strict-decay assertion. The $k=0$ mode has $F=0$, an epicyclic motion sector and a zero-frequency epicyclic guiding center sector with possible secular shear drift. These neutral or algebraic behaviors are distinct from the strict exponential instability established by the requested inequality.

The Euler-Lagrange equations of the particle Lagrangian in a shearing sheet areThe terms coupling $x$ and $y$ are the Coriolis acceleration. The cyclic coordinate $y$ has conserved canonical momentumFor a Newtonian potential of a point mass, $q_0=3/2$ and ${\Omega }_z={\kappa }_r={\Omega }_0$. Set $p_y={\Omega }_0{x}_0/2$ and substitute $\dot{y}=p_y-2{\Omega }_{0}x$ in the radial equation. It becomes $\ddot{x}+{\Omega }_0^2(x-x_0)=0$, a harmonic oscillator equation about the epicyclic guiding center $x_0$. Integration givesThe four real constants in $A,y_0,x_0$ and the two in $B$ account for the six initial position and velocity data.

Expanding the inertial specific angular momentum gives $h=h_0+r_0{p}_y+O(2)$. Thus $p_y$ measures the angular-momentum offset from the reference circular orbit, and $x_0$ specifies the radius of its associated epicyclic guiding center.

The conserved horizontal energy in the rotating frame isIt is the horizontal part of the local Jacobi energy in a shearing sheet, rather than the inertial specific orbital energy. Its positive term is the epicyclic energy; its negative term is the energy of the background shear at guiding-center position $x_0$. Independently, the vertical harmonic oscillator has conserved energy

Initially every particle has $A=B=0$, so its Jacobi energy in a shearing sheet is $-3{\Omega }_0^2{x}_0^2/8$. During an inelastic collision, positions are fixed at the instant of impact, while momentum conservation preserves the sum of the tangential velocities. Hence the sum of $p_y=\dot{y}+2{\Omega }_{0}x$, and therefore the sum of the epicyclic guiding center positions $x_0=2p_y/{\Omega }_0$, is unchanged. Between collisions, these quantities are individually conserved.

The collision dissipates kinetic energy without changing the instantaneous tidal potential. Consequently the total Jacobi energy in a shearing sheet decreases. At any later time, writing $x_{0j}$ for the particles' current epicyclic guiding center positions givesIf $D\ge 0$ is the accumulated energy dissipated in the inelastic collisions, comparison with the initial circular orbits givesThus the mean guiding-center position stays fixed, while its variance grows. The initial ensemble has no preference for positive or negative $x$, and the local equations and collision law preserve the symmetry $(x,y,z)\mapsto (-x,-y,z)$. The spreading is therefore symmetric in the ensemble average: angular momentum transport moves some particles inward and others outward. A particular finite random realization need not be exactly symmetric. This is the microscopic energy argument for dissipative spreading of a planetary ring.