pith. sign in

arxiv: 2606.23912 · v1 · pith:6X4Z2DSZnew · submitted 2026-06-22 · 🌀 gr-qc · astro-ph.EP· astro-ph.HE· astro-ph.SR

Newtonian Shirokov Effect: Epicyclic Frequency Splitting from Mass Multipoles

Anuar Idrissov , Kuantay Boshkayev , Abylaikhan Tlemissov , Hernando Quevedo This is my paper

Pith reviewed 2026-06-26 07:05 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.EPastro-ph.HEastro-ph.SR
keywords epicyclic frequenciesmass multipolesNewtonian gravityquadrupole momentorbital oscillationsShirokov effectaxisymmetric potentialcenter-of-mass offset
0
0 comments X

The pith

Quadrupole and higher multipoles split radial and vertical epicyclic frequencies in axisymmetric Newtonian potentials, except the removable dipole.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that in an axisymmetric Newtonian gravitational potential expanded in mass multipoles, nearly circular orbits experience splitting between the radial and vertical epicyclic frequencies due to all multipoles except the dipole. This follows from evaluating the full Hessian of the effective potential at the true, possibly tilted, equilibrium position and exactly solving the resulting two-mode coupled oscillator. A quadrupole produces the splitting Ω_θ^{2} − Ω_r^{2} = −3 G Q / r_{0}^{5} = 6 G M J_{2} R^{2} / r_{0}^{5} at first order in J_{2}, while a genuine octupole produces ω_{+}^{2} − ω_{-}^{2} ≈ 6 G |O| / r_{0}^{6}; the dipole yields no splitting because it equals a center-of-mass shift removable by re-centering. The result supplies a classical counterpart to the relativistic Shirokov effect and two independent observables: frequency splitting for oblateness and orbital tilt for center-of-mass offset, together with explicit solar-system estimates for the secular transverse drift.

Core claim

In an axisymmetric Newtonian potential expanded in multipoles, every multipole except the dipole splits the epicyclic frequencies of nearly circular orbits. A quadrupole gives Ω_θ^{2} − Ω_r^{2} = −3 G Q / r_{0}^{5} = 6 G M J_{2} R^{2} / r_{0}^{5} at leading order. A genuine octupole gives ω_{+}^{2} − ω_{-}^{2} ≈ 6 G |O| / r_{0}^{6}. The dipole produces no splitting because it is equivalent to a center-of-mass displacement that can be removed by coordinate shift; any apparent first-order coupling vanishes once the true tilted equilibrium is used. The secular transverse drift after n orbits is ξ^θ = ξ_{0}^θ π n (6 J_{2} R^{2} / r_{0}^{2}), of order 10^{-8} cm at 1 au.

What carries the argument

The Hessian of the effective potential evaluated at the true (possibly tilted) equilibrium, which supplies the exact frequencies of the coupled two-mode radial-vertical oscillator.

If this is right

  • Frequency splitting directly measures the quadrupole moment J_{2} of the central body.
  • Orbital-plane tilt δ heta_{0} ≈ −r_CM / r_{0} measures the center-of-mass offset as a geometric observable.
  • The secular transverse drift scales linearly with J_{2} and with the number of orbits n.
  • Solar-system estimates give drifts of order 10^{-8} cm at 1 au and 10^{-6} cm near 0.1 au.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Frequency splitting measurements in multi-planet systems could constrain higher multipoles without reference to relativistic corrections.
  • The removal of dipole effects by coordinate shift implies that any coordinate-independent observable must be insensitive to small center-of-mass displacements.
  • The same Hessian approach could be applied to galactic or stellar-cluster potentials to test whether observed epicyclic splittings require non-axisymmetric terms.

Load-bearing premise

The gravitational potential is axisymmetric and expandable in mass multipoles, with the orbit remaining nearly circular so the Hessian at the true equilibrium fully describes the small-oscillation dynamics.

What would settle it

Numerical integration of test-particle orbits in an axisymmetric potential containing only a known octupole term, after shifting to the true equilibrium, that shows identical radial and vertical frequencies would falsify the octupole splitting result.

read the original abstract

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, $\Omega_\theta^2-\Omega_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_{\rm 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\approx6G|O|/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\simeq-r_{\rm CM}/r_0$, measures the center-of-mass offset $r_{\rm 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 $\sim10^{-6}$ cm near $0.1$ au, comparable to Shirokov's Schwarzschild estimate.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The paper analyzes small oscillations of nearly circular orbits in an axisymmetric Newtonian potential expanded in mass multipoles, as the classical counterpart to the relativistic Shirokov effect. It computes the Hessian of the effective potential V_eff at the true (possibly tilted) equilibrium, solves the resulting coupled two-mode oscillator exactly, and derives frequency splittings: quadrupole gives Ω_θ² − Ω_r² = −3 G Q / r₀⁵ (= 6 G M J₂ R² / r₀⁵), dipole produces none after re-centering at the center of mass (with apparent coupling canceling and confirmed by direct numerical integration), octupole splits with ω₊² − ω₋² ≈ 6 G |O| / r₀⁶; every multipole except dipole splits frequencies. It also maps to secular transverse drift ξ^θ after n orbits and gives solar-system estimates.

Significance. If the derivations hold, this supplies a parameter-free Newtonian derivation of epicyclic splitting from multipoles, with explicit strengths in the exact oscillator solution and numerical validation of dipole cancellation. It yields two independent, falsifiable probes of axisymmetric sources (frequency splitting for oblateness J₂; orbital tilt δθ₀ for r_CM) and a concrete mapping to Shirokov's secular-drift observable, with order-of-magnitude estimates at 1 au and 0.1 au. The selection rule (all multipoles except dipole) is a clear outcome of the standard effective-potential construction.

minor comments (2)
  1. [Abstract] Abstract and § on quadrupole: the two expressions for the splitting are equated without an intermediate step showing how Q is normalized to J₂; adding one line relating the quadrupole moment to the standard J₂ definition would remove ambiguity.
  2. [Dipole section] Numerical confirmation paragraph: the direct orbit integration for the dipole is cited as confirmation, but no details are given on integrator, timestep, or convergence criterion; a brief methods sentence would aid reproducibility.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary, significance assessment, and recommendation of minor revision. No major comments were provided in the report, so we have no specific points requiring response or revision.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper derives the epicyclic frequency splitting directly from the Newtonian multipole-expanded potential by constructing the axisymmetric effective potential V_eff = Φ_multipole + L_z²/(2r²sin²θ), locating the (possibly tilted) equilibrium via first derivatives, and extracting the eigenvalues of the 2×2 Hessian at that point. All splitting expressions (quadrupole: Ω_θ²−Ω_r² = −3GQ/r₀⁵; octupole: ω₊²−ω₋² ≈ 6G|O|/r₀⁶; dipole cancellation after re-centering) follow from this standard small-oscillation analysis without any fitted parameters renamed as predictions, self-citation load-bearing steps, or ansatzes imported from prior author work. The selection rule and secular-drift mapping are obtained by explicit expansion and eigenvalue solution, remaining independent of the target results.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on the standard Newtonian gravitational potential and the assumption that it admits a multipole expansion for an axisymmetric source; no new entities are introduced and the multipole coefficients are treated as given inputs rather than fitted parameters.

axioms (2)
  • standard math Newtonian gravitational potential
    The entire analysis is performed in the Newtonian limit.
  • domain assumption Axisymmetric mass distribution admitting multipole expansion
    Invoked to justify the potential form used for the effective potential and its Hessian.

pith-pipeline@v0.9.1-grok · 5972 in / 1339 out tokens · 30993 ms · 2026-06-26T07:05:19.509420+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

38 extracted references · 3 linked inside Pith

  1. [1]

    The selection rule is thus not even/odd parity: every multipole splits the frequenciesexceptthe dipole. These results yield two complementary, non-overlapping probes of an axisymmetric source from orbital data: the frequency splittingmeasures the intrinsic oblatenessJ2, while theorbital-plane tilt,δθ 0 ≃ −r CM/r0 relative to the symmetry axis, measures th...

    Pith/arXiv arXiv 2026

  2. [2]

    This is the genuine Newtonian analogue of the Shirokov splitting

    The two frequencies arenot shifted equally: the quadrupole splits them at first order inQ(equivalently inJ 2). This is the genuine Newtonian analogue of the Shirokov splitting. It is equivalent to the classical statement that an oblate body produces apsidal precession at the rateΩϕ −Ωr and nodal regression at the rateΩ ϕ−Ωθ, withthetworatesunequalprecisel...

  3. [3]

    dipole on a monopole

    Their expansion runs over the even multipolesJ 2m alone, so it contains neither the dipole nor the octupole treated below. IV. THE DIPOLE DOESNOTSPLIT THE FREQUENCIES We now keeprCM ̸= 0(withO= 0) and ask whether the dipole adds any physical splitting. The answer is no, and the reason is best stated physically before any calculation. A. The physical argum...

  4. [4]

    dipoleonly

    The shift scales asr−2 0 and is independent of the particle mass. D. Numerical estimate for the Sun Taking the solar values of Table II (J⊙ 2 = 2×10 −7, R⊙ = 6.96×10 8m) with Shirokov’s choicesn= 10,ξ0 = TABLE I. Measured epicyclic ratioΩr/Ωθ from direct orbit integration (G=M=r 0 = 1). Unity means no splitting. Source (isolated)Ω r/Ωθ Pure monopole (base...

  5. [5]

    The two are genuinely disjoint: Probe A reads the plane orientation (set by the dipole, i.e

    Probe B is coordinate-independent and needs only the orbital and equatorial radii as inputs. The two are genuinely disjoint: Probe A reads the plane orientation (set by the dipole, i.e. byrCM), Probe B reads the frequency splitting (set by the quadrupole, i.e. byJ 2). The dipole contributes only to A, the quadrupole only to B. Their agreement, where both ...

    2025

  6. [6]

    M. F. Shirokov, On one new effect of the Einsteinian theory of gravitation, General Relativity and Gravitation 4, 131 (1973)

    1973

  7. [7]

    Idrissov, K

    A. Idrissov, K. Boshkayev, S. Momynov, H. Quevedo, D. Utepova, A. Urazalina, and B. Baitimbetova, Shi- rokov and Shapiro Effects in the Hartle-Thorne Space- time, arXiv e-prints (2026), arXiv:2605.01639 [gr-qc]

    Pith/arXiv arXiv 2026

  8. [8]

    E. Y. Melkumova and V. I. Khlebnikov, Calculation of the Shirokov effect for circular orbits in an axisymmetric static gravitational field, Soviet Physics Journal33, 349 (1990)

    1990

  9. [9]

    C. M. Will, The Confrontation between General Relativ- ity and Experiment, Living Reviews in Relativity17, 4 (2014)

    2014

  10. [10]

    C. W. F. Everittet al., Gravity Probe B: Final Results of a Space Experiment to Test General Relativity, Physical Review Letters106, 221101 (2011)

    2011

  11. [11]

    Ciufoliniet al., A test of general relativity using the LARES and LAGEOS satellites and a GRACE Earth gravity model, European Physical Journal C76, 120 (2016)

    I. Ciufoliniet al., A test of general relativity using the LARES and LAGEOS satellites and a GRACE Earth gravity model, European Physical Journal C76, 120 (2016)

    2016

  12. [12]

    J. P. Rozelot and C. Damiani, History of solar oblate- ness measurements and interpretation, European Physi- cal Journal H36, 407 (2011)

    2011

  13. [13]

    Müller, R

    D. Müller, R. G. Marsden, O. C. St. Cyr, H. R. Gilbert, C. Garcia-Marirrodriga, P. Geithner, B. Schmieder, M. Accorsi, L. Ali, S. Aljanabi,et al., The Solar Orbiter mission. Science overview, Astronomy & Astrophysics 642, A1 (2020)

    2020

  14. [14]

    N. M. Günter,Potential Theory and Its Applications to Basic Problems of Mathematical Physics(Frederick Un- gar Publishing Co., New York, 1967) translated from Russian

    1967

  15. [15]

    Quevedo, Multipole Moments in General Relativity - Static and Stationary Vacuum Solutions, Fortschritte der Physik38, 733 (1990)

    H. Quevedo, Multipole Moments in General Relativity - Static and Stationary Vacuum Solutions, Fortschritte der Physik38, 733 (1990)

    1990

  16. [16]

    Binney and S

    J. Binney and S. Tremaine,Galactic Dynamics, 2nd ed. (Princeton University Press, Princeton, 2008)

    2008

  17. [17]

    C. D. Murray and S. F. Dermott,Solar System Dynamics (Cambridge University Press, Cambridge, 1999)

    1999

  18. [18]

    R. M. Wald,General Relativity(The University of Chicago Press, Chicago, 1984)

    1984

  19. [19]

    Kluźniak and M

    W. Kluźniak and M. Rosińska, Orbital and epicyclic fre- quencies of a spheroidal mass distribution, Monthly No- tices of the Royal Astronomical Society434, 2825 (2013)

    2013

  20. [20]

    M. A. Abramowicz and W. Kluźniak, A precise determi- nation of black hole spin in GRO J1655-40, Astronomy and Astrophysics374, L19 (2001)

    2001

  21. [21]

    Boshkayev, A

    K. Boshkayev, A. Idrissov, O. Luongo, and M. Muccino, Quasiperiodic oscillations for spherically symmetric reg- ular black holes, Physical Review D108, 044063 (2023)

    2023

  22. [22]

    Utepova, K

    D. Utepova, K. Boshkayev, S. Momynov, A. Idrissov, G. Baimbetova, H. Quevedo, and A. Urazalina, Con- straining quadrupole deformations with relativistic ef- fects, European Physical Journal C85, 987 (2025)

    2025

  23. [23]

    C. W. Misner, K. S. Thorne, and J. A. Wheeler,Gravi- tation(W. H. Freeman, San Francisco, 1973)

    1973

  24. [24]

    Toshmatov, D

    B. Toshmatov, D. Malafarina, and N. Dadhich, Harmonic oscillations of neutral particles in theγmetric, Physical Review D100, 044001 (2019), arXiv:1905.01088 [gr-qc]

    arXiv 2019

  25. [25]

    Iorio, Post-Keplerian corrections to the orbital periods of a two-body system and their measurability, Monthly Notices of the Royal Astronomical Society460, 2445 (2016)

    L. Iorio, Post-Keplerian corrections to the orbital periods of a two-body system and their measurability, Monthly Notices of the Royal Astronomical Society460, 2445 (2016)

    2016

  26. [26]

    L. Iorio, Analytically calculated post-keplerian range and range-rate perturbations for planetary radar ranging: the solar lense–thirring effect and bepicolombo, Monthly No- tices of the Royal Astronomical Society476, 1811 (2018)

    2018

  27. [27]

    R. S. Park, W. M. Folkner, J. G. Williams, and D. H. Boggs, The JPL Planetary and Lunar Ephemerides DE440 and DE441, The Astronomical Journal161, 105 (2021)

    2021

  28. [28]

    W. M. Folkner, J. G. Williams, D. H. Boggs, R. S. Park, and P. Kuchynka, The planetary and lunar ephemerides de430 and de431, Interplanetary Network Progress Re- port42-196, 1 (2014)

    2014

  29. [29]

    R. N. Manchester, The international pulsar timing array, Classical and Quantum Gravity30, 224010 (2013)

    2013

  30. [30]

    J. P. W. Verbiestet al., The International Pulsar Timing Array: First data release, Monthly Notices of the Royal Astronomical Society458, 1267 (2016)

    2016

  31. [31]

    D. J. Championet al., Measuring the Mass of Solar System Planets Using Pulsar Timing, The Astrophysi- cal Journal Letters720, L201 (2010)

    2010

  32. [32]

    E. V. Pitjeva and N. P. Pitjev, Mass and the location of mass concentration in the Kuiper belt, Celestial Mechan- ics and Dynamical Astronomy130, 57 (2018). 12

    2018

  33. [33]

    Agazieet al.(NANOGrav), The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Back- ground, The Astrophysical Journal Letters951, L8 (2023)

    G. Agazieet al.(NANOGrav), The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Back- ground, The Astrophysical Journal Letters951, L8 (2023)

    2023

  34. [34]

    M. Vallisneriet al.(NANOGrav), Modeling the Un- certainties of Solar System Ephemerides for Robust Gravitational-wave Searches with Pulsar-timing Arrays, The Astrophysical Journal893, 112 (2020)

    2020

  35. [35]

    Benkhoffet al., BepiColombo - Mission Overview and Science Goals, Space Science Reviews217, 90 (2021)

    J. Benkhoffet al., BepiColombo - Mission Overview and Science Goals, Space Science Reviews217, 90 (2021)

    2021

  36. [36]

    di Stefano, P

    I. di Stefano, P. Cappuccio, and L. Iess, The Bepi- Colombo MORE gravimetry and rotation experiments with the ORBIT14 software, Classical and Quantum Gravity38, 055002 (2021)

    2021

  37. [37]

    Iorio, Constraining the angular momentum of the sun with planetary orbital motions and general relativity, So- lar Physics244, 125 (2007)

    L. Iorio, Constraining the angular momentum of the sun with planetary orbital motions and general relativity, So- lar Physics244, 125 (2007)

    2007

  38. [38]

    196 (Springer, Cham, 2019) p

    D.Philipp, D.Puetzfeld,andC.Lämmerzahl,OntheAp- plicability of the Geodesic Deviation Equation in General Relativity, inRelativistic Geodesy, Fundamental Theo- ries of Physics, Vol. 196 (Springer, Cham, 2019) p. 419, arXiv:1604.07173 [gr-qc]

    Pith/arXiv arXiv 2019