Canceling Effects of Conjunctions Render Higher Order Mean Motion Resonances Weak

Mean motion resonances (MMRs) are a key phenomenon in orbital dynamics. The traditional disturbing function expansion in celestial mechanics shows that, at low eccentricities, $p$:$p-q$ MMRs exhibit a clear hierarchy of strengths, scaling as $e^q$, where $q$ is the order of the resonance. This explains why first-order MMRs (e.g., 3:2 and 4:3) are important, while the infinite number of higher order integer ratios are not. However, this relationship derived from a technical perturbation series expansion provides little physical intuition. In this paper, we provide a simple physical explanation of this result for closely spaced orbits. In this limit, interplanetary interactions are negligible except during close encounters at conjunction, where the planets impart a gravitational “kick” to each other’s mean motion. We show that while first-order MMRs involve a single conjunction before the configuration repeats, higher order MMRs involve multiple conjunctions per cycle, whose effects cancel out more precisely the higher the order of the resonance. Starting from the effects of a single conjunction, we provide an alternate, physically motivated derivation of MMRs’ $e^q$ strength scaling.

💡 Research Summary

The paper tackles a long‑standing question in celestial mechanics: why does the strength of a mean‑motion resonance (MMR) scale as e^q at low eccentricities, where q is the resonance order? Traditional derivations rely on a perturbative expansion of the disturbing function, which mathematically yields the e‑dependence but offers little physical intuition. The authors present a simple, physically motivated picture that works in the limit of closely spaced (Hill‑limit) coplanar orbits. In this regime the mutual gravitational interaction between two planets is negligible except during the brief moments of conjunction, when the planets pass each other at a small separation and impart an impulsive “kick” to each other’s mean motion.

The analysis begins by mapping the general two‑planet problem onto the circular restricted three‑body problem (CR3BP) in the Hill limit. The inner planet is taken to be on a fixed circular orbit, while the outer test particle follows an eccentric orbit characterized by a normalized eccentricity ẽ = e/e_c, where e_c is the eccentricity at which the orbits would cross. The conjunction angle is defined as θ = λ_conj − ϖ, i.e., the angular offset of the conjunction from the line of closest approach.

A key result is the Fourier series for the fractional change in the test particle’s mean motion at a conjunction:

δn/n_p ≈ 2π μ e_c⁻² ∑{j=1}^∞ A{2j} ẽ^j sin(jθ),

where μ = m/M_* is the planet‑to‑star mass ratio, and the coefficients A_{2j} are order‑unity numbers tabulated by Tamayo & Hadden (2025). This expression shows that a single kick contains contributions from all harmonics sin(jθ), with higher‑j terms rapidly diminishing for small ẽ.

The authors then treat the evolution of θ as a continuous variable, deriving a second‑order differential equation by relating the time derivative of the resonant angle to the mean‑motion change. Smoothing the discrete kicks over many conjunctions yields

¨θ = (p q) δn/t_conj ≈ C_q sin(qθ),

which is the equation of a simple pendulum. For a first‑order resonance (q = 1) the coefficient

C₁ = A₁ μ ẽ/e_c

matches the coefficient of the leading cosine term in the classical disturbing function. The pendulum Hamiltonian H = ½ θ̇² + C₁ cosθ describes the familiar libration dynamics of a first‑order MMR.

For higher‑order resonances (q > 1) the situation is more subtle. In a p : (p − q) resonance the system experiences q conjunctions per resonant cycle. Each conjunction contributes a kick with a different phase θ, and the Fourier expansion shows that only terms with j = k q (k = 1, 2, …) survive after summing over the q kicks; all other harmonics cancel to leading order. Consequently the net effect after a full resonant cycle scales as ẽ^q, because the surviving term is proportional to ẽ^{q} sin(qθ). This cancellation explains why higher‑order resonances are intrinsically weaker at low eccentricities: the multiple kicks destructively interfere, leaving only a residual that is higher‑order in ẽ.

The theoretical predictions are validated with a suite of REBOUND N‑body integrations. The authors initialize many simulations with a single conjunction occurring at various θ values, measure the resulting fractional change in mean motion, and compare the results to the truncated Fourier series. The numerical data align closely with the analytic curves, confirming that the first harmonic dominates for small ẽ and that adding higher harmonics reproduces the exact kick profile.

In conclusion, the paper provides a clear physical mechanism for the e^q scaling of MMR strength: in the Hill limit, resonant dynamics are governed by impulsive kicks at conjunctions, and the number of conjunctions per resonant cycle determines the degree of cancellation among the Fourier components of those kicks. This framework reproduces the classic disturbing‑function results while offering intuitive insight, and it can be extended to more complex systems (multiple planets, disks, satellite systems) where close orbital spacing makes the impulsive‑kick picture applicable. The work bridges the gap between formal perturbation theory and a tangible dynamical picture, deepening our understanding of why only low‑order resonances dominate planetary system architecture.