Newtonian Shirokov Effect: Epicyclic Frequency Splitting from Mass Multipoles

We analyze small oscillations of nearly circular orbits in an axisymmetric Newtonian potential expanded in mass multipoles, as the classical counterpart of the relativistic Shirokov effect. Computing the full Hessian of the effective potential at the true (possibly tilted) equilibrium and solving the coupled two-mode oscillator exactly, we obtain a complete picture. (i) A quadrupole splits the radial and vertical epicyclic frequencies, theta 2- r 2=-3GQ/r 0 5=6GMJ 2R 2/r 0 5, at first order in J 2; the Newtonian analogue of the Shirokov splitting, equivalent to the classical statement that an oblate body's apsidal and nodal rates differ. (ii) A gravitational dipole produces no splitting: it equals M r CM, is removable by re-centering at the center of mass, and cannot appear in any coordinate independent frequency; the apparent first order coupling cancels at the true tilted equilibrium, any residual absorbed by the induced quadrupole of the shifted source, confirmed by direct orbit integration. (iii) A genuine octupole does split the frequencies, omega+ 2-omega-2 about6GO/r 0 6. The selection rule is thus not even/odd parity: every multipole splits the frequencies except the dipole. These yield two complementary probes of an axisymmetric source: the frequency splitting measures the oblateness J 2, while the orbital plane tilt, delta theta 0-r CM /r 0, measures the center of mass offset r CM, an orbital geometry observable rather than a frequency one. We give solar system estimates for both. Carried through to Shirokov's original observable -- the secular transverse drift after n orbits -- the quadrupole effect gives xi theta= xi 0 theta, pi n,(6J 2R 2/r 0 2), of order 10-8 cm at 1 au and about10-6 cm near 0.1 au, comparable to Shirokov's Schwarzschild estimate.