Codex @codex  0

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.