Newtonian Shirokov Effect: Epicyclic Frequency Splitting from Mass Multipoles

Technical Summary: Newtonian Shirokov Effect: Epicyclic Frequency Splitting from Mass Multipoles

Problem Statement
The paper addresses the origin of epicyclic frequency splitting—the difference between radial ($\Omega_r$) and vertical ($\Omega_\theta$) oscillation frequencies of a test particle in a nearly circular orbit. In General Relativity (GR), this splitting (the Shirokov effect) arises from spacetime curvature even for a static, spherically symmetric source (Schwarzschild metric). However, for realistic astrophysical bodies, the gravitational potential includes Newtonian multipole corrections. The authors investigate the "classical counterpart" of the Shirokov effect: specifically, which Newtonian mass multipoles (monopole, dipole, quadrupole, octupole, etc.) induce a splitting of these frequencies in a static, axisymmetric potential, and how this splitting relates to the source's physical properties versus coordinate choices.

Methodology
The authors analyze small oscillations of a test particle in an axisymmetric Newtonian potential expanded in mass multipoles.

1. Potential Expansion: The external potential $\Phi(r)$ is expanded to include the monopole, dipole ($D_i$), trace-free quadrupole ($Q_{ij}$), and trace-free octupole ($O_{ijk}$). The dipole is explicitly linked to the displacement of the coordinate origin from the center of mass ($r_{CM}$).
2. Effective Potential & Linearization: The effective potential $V_{eff}$ is constructed by eliminating azimuthal motion via conservation of specific angular momentum. The system is linearized around a circular equilibrium $(r_0, \theta_0)$.
3. Hessian Analysis: The authors compute the full Hessian matrix of $V_{eff}$ at the true equilibrium. Crucially, for multipoles that break reflection symmetry (dipole, octupole), the equilibrium is not at the equatorial plane ($\theta = \pi/2$) but on a tilted cone ($\theta_0 \neq \pi/2$).
4. Eigenvalue Solution: The coupled two-mode oscillator equations are solved exactly to find the normal-mode frequencies $\omega_\pm$. The splitting is defined as $\omega_+^2 - \omega_-^2$ (or $\Omega_\theta^2 - \Omega_r^2$ when coupling vanishes).
5. Verification: Analytic results are confirmed via direct numerical integration of orbits in the exact fields of displaced point masses and specific multipole configurations.

Key Contributions and Results

• The Quadrupole Splitting (Genuine Newtonian Shirokov Effect):
  The paper establishes that a genuine mass quadrupole (associated with oblateness $J_2$) splits the radial and vertical epicyclic frequencies at first order. The splitting is given by:
  $$\Omega_\theta^2 - \Omega_r^2 = -\frac{3GQ}{r_0^5} = \frac{6GMJ_2R^2}{r_0^5}$$
  This result is identified as the true Newtonian analogue of the relativistic Shirokov splitting. It corresponds to the classical fact that an oblate body produces different rates for apsidal precession ($\Omega_\phi - \Omega_r$) and nodal regression ($\Omega_\phi - \Omega_\theta$).

• The Dipole Does Not Split Frequencies:
  Contrary to the intuition that a dipole (representing an off-center mass) might break symmetry, the authors prove it produces no frequency splitting.

  • Physical Reason: The dipole moment $D = Mr_{CM}$ encodes only the choice of coordinate origin relative to the center of mass. It is a "pure gauge" quantity.
  • Mechanism: While a naive calculation at the equator suggests a first-order coupling term ($C \neq 0$), evaluating the Hessian at the true tilted equilibrium (where the orbit naturally tilts by $\delta\theta_0 \approx -r_{CM}/r_0$) causes this coupling to vanish exactly. Furthermore, any second-order residual is cancelled by the induced quadrupole of the shifted source.
  • Conclusion: The frequencies remain degenerate ($\Omega_r = \Omega_\theta$) for a displaced point mass, regardless of the coordinate origin.

• The Octupole Splitting:
  A genuine octupole (measuring north-south asymmetry or "pear-shape") does split the frequencies. Unlike the dipole, the octupole is a genuine degree of freedom of the mass distribution. The splitting occurs at first order in the octupole moment $O$:
  $$\omega_+^2 - \omega_-^2 \approx \frac{6G|O|}{r_0^6}$$
  The splitting depends on the magnitude $|O|$ and is independent of the sign (direction of asymmetry).

• Selection Rule:
  The paper corrects the naive expectation that splitting depends on parity (even vs. odd multipoles). The actual rule is: Every mass multipole splits the epicyclic frequencies except the dipole. The dipole is unique because it is removable by a coordinate translation, whereas the quadrupole, octupole, and higher genuine moments are intrinsic properties of the source.

• Observable Probes:
  The results yield two complementary, non-overlapping probes for an axisymmetric source:

  1. Probe A (Geometry): The tilt of the orbital plane ($\delta\theta_0$) measures the center-of-mass offset ($r_{CM}$). This is a geometric observable, not a frequency observable.
  2. Probe B (Frequencies): The frequency splitting measures the intrinsic oblateness ($J_2$).

• Secular Drift:
  The authors translate the quadrupole splitting into Shirokov's original observable: the secular transverse drift ($\xi_\theta$) of a test particle after $n$ orbits. For the Newtonian quadrupole, the drift is:
  $$\xi_\theta(n) \approx \xi_{\theta,0} \pi n \left( \frac{6J_2R^2}{r_0^2} \right)$$
  Estimates for the Sun suggest this drift is $\sim 10^{-8}$ cm at 1 au, comparable to the relativistic Schwarzschild estimate only at very close distances ($\sim 0.1$ au). For Earth satellites, the Newtonian quadrupole drift dominates the relativistic effect by a factor of $\sim 10^6$.

Significance and Claims
The paper claims to isolate the purely Newtonian content of the Shirokov effect, demonstrating that frequency splitting is a probe of the source's intrinsic structure (shape) rather than coordinate artifacts.

• It clarifies that the "Shirokov effect" in a Newtonian context is driven by the quadrupole moment, not the dipole.
• It distinguishes between coordinate-dependent quantities (center-of-mass offset, detectable only via orbital tilt) and coordinate-independent dynamical quantities (frequency splitting, detectable only via oblateness).
• It provides a rigorous analytic and numerical confirmation that the dipole, despite breaking reflection symmetry, does not affect the eigenfrequencies of the system when the equilibrium is correctly identified.
• The work serves as a foundation for interpreting orbital data from solar-system bodies (e.g., determining solar $J_2$ from orbital data) and highlights that for bodies with significant oblateness (like Earth), Newtonian quadrupole effects can vastly exceed relativistic corrections in epicyclic frequency splitting.