arxiv: 2605.05470 · v1 · submitted 2026-05-06 · 🌌 astro-ph.EP

Recognition: unknown

Resonant Hamiltonian Dynamics in the CR3BP: Bistability and Stochastic Resonance in Binary Planetary Systems

R. Capuzzo-Dolcetta

Authors on Pith no claims yet

Pith reviewed 2026-05-08 15:23 UTC · model grok-4.3

classification 🌌 astro-ph.EP

keywords mean-motion resonancecircular restricted three-body problembistable effective potentialstochastic resonancebinary star systemsexoplanet dynamicsHamiltonian averagingfinite mass ratio

0 comments

The pith

The effective potential for mean-motion resonances in binary star systems becomes bistable when the ratio of the first two harmonic amplitudes exceeds one quarter.

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

The paper constructs a single Hamiltonian description that covers both planets orbiting one star and planets orbiting both stars in a binary, with the stellar companion carrying finite mass. After transforming to resonant variables and averaging away the fast orbital motion, the reduced effective potential develops two stable wells once the second resonant harmonic grows stronger than one quarter the strength of the first. This bistable landscape supplies the necessary condition for stochastic resonance, in which weak noise can drive the planet to switch between the two wells. A reader would care because the same mechanism could govern long-term orbital stability or migration in real binary systems, and the derived scaling laws point to specific extreme binaries where the effect might be observable.

Core claim

Starting from the full CR3BP Hamiltonian in both S-type and P-type configurations, canonical transformation to resonant action-angle variables followed by averaging over fast angles produces a one-degree-of-freedom Hamiltonian whose effective potential exhibits bistability precisely when |ε₂/ε₁| > 1/4. Leading-order Fourier scaling laws for ε₁ and ε₂ are obtained as functions of binary mass ratio and semimajor-axis ratio; current observed systems fall below the threshold, while sufficiently close, near-equal-mass binaries are predicted to enter the bistable regime.

What carries the argument

The reduced one-degree-of-freedom resonant Hamiltonian obtained by averaging the perturbing potential over fast orbital motion, whose shape is controlled by the relative amplitudes of the first two Fourier harmonics.

If this is right

Bistable resonances become possible only in extreme P-type configurations with a/a_b ≲ 1.5 and near-equal stellar masses.
All currently catalogued binary-planet systems lie below the |ε₂/ε₁| = 1/4 threshold and therefore lack the bistable setting.
The same averaged Hamiltonian framework applies uniformly to both circumstellar and circumbinary resonant orbits once the finite-mass perturber is retained.
Stochastic resonance can in principle operate once the bistability condition is met, allowing noise to induce transitions between the two resonant wells.

Where Pith is reading between the lines

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

Detection of orbital-element switching on timescales set by the noise strength in future extreme-binary observations would constitute direct evidence of the predicted stochastic resonance.
The same reduction technique could be applied to higher-order mean-motion resonances or to systems with small eccentricities to test whether bistability persists.
If the bistability threshold is crossed in nature, the resulting two-state dynamics would impose a new constraint on the long-term survival of planets in tight binaries.

Load-bearing premise

Systematic averaging over fast orbital motion together with the leading-order Fourier scaling laws fully capture the resonant dynamics for finite-mass-ratio binaries without significant contamination from higher-order terms or non-averaged effects.

What would settle it

Numerical integration of the unaveraged CR3BP equations for a binary-planet configuration where the calculated |ε₂/ε₁| exceeds 1/4, followed by inspection of the long-term orbital elements to check whether two distinct stable resonant states appear and can be switched by added noise.

Figures

Figures reproduced from arXiv: 2605.05470 by R. Capuzzo-Dolcetta.

**Figure 1.** Figure 1: Schematic geometry of S-type (circumprimary) and P view at source ↗

**Figure 2.** Figure 2: Effective potential V(θ) for 3 choices of |ǫ2/ǫ1| as labelled. The case |ǫ2/ǫ1| = 0.4 > 1/4 shows bistability with a barrier height ∆V = 0.11. The two minima are indicated by the red dots and the maximum by the brown square. and the barrier height is ∆V = V(π) − V(θmin) view at source ↗

**Figure 3.** Figure 3: Resonant phase portraits showing the bifurcation of view at source ↗

read the original abstract

Context: The Circular Restricted Three-Body Problem provides a fundamental framework for understanding resonant dynamics in binary star systems. Aims: We develop a unified Hamiltonian formulation for mean-motion resonances that encompasses both circumstellar and circumbinary planetary orbits within the CR3BP. Unlike the Solar System case where the perturbing body is a planet of negligible mass, here the perturber (a stellar companion) has a non-negligible, finite mass, a crucial difference that we fully incorporate. Methods: Starting from the full Hamiltonian in each configuration, we perform canonical transformations to resonant action angle variables and derive reduced one-degree-of-freedom Hamiltonians through systematic averaging over the fast orbital motion. Leading-order scaling laws for the Fourier coefficients of the resonant perturbation are obtained, revealing their dependence on the binary mass ratio and the planet's orbital distance. Results: The resulting effective potential is shown to exhibit bistability under the well-defined condition |epsilon2/epsilon1| > 1/4, where epsilon1 and epsilon2 are the amplitudes of the first two resonant harmonics. This bistability creates the essential dynamical setting for stochastic resonance. Scaling laws for the Fourier coefficients are derived for both S-type and P-type configurations. Estimates for known binary-planet systems show that while currently observed systems lie below the bistability threshold, the theory predicts that extreme configurations (a/a_b <~ 1.5 for P-type, almost equal mass binary) could host bistable resonances accessible to future observations. Conclusions: This work provides a natural Hamiltonian framework for studying stochastic resonance in binary planetary systems, bridging analytical celestial mechanics and the nonlinear dynamics of exoplanetary systems subject to realistic perturbations.

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.

Desk Editor's Note private letter to a colleague

The paper derives scaling laws and a |ε2/ε1| > 1/4 bistability threshold in averaged resonant Hamiltonians for finite-mass binaries, but the claim rests on unverified truncation that may not survive in strong-perturbation P-type cases.

read the letter

The core advance here is extending the standard CR3BP resonant averaging to binaries where the perturber has finite mass, then extracting explicit mass-ratio and separation scalings for the leading Fourier terms in both S-type and P-type geometries. That produces a clean condition for bistability in the reduced effective potential and ties it to stochastic resonance, which is a reasonable next step beyond the usual negligible-mass treatments used for Solar-System resonances. The derivations appear internally consistent on the abstract level and avoid obvious circularity by starting from the full Hamiltonian and performing the canonical transformations and averaging in the usual way. Credit for that systematic approach. The soft spot is exactly the one flagged in the stress test. For P-type orbits with mass ratio near 0.5 and a/ab ≲ 1.5 the perturbation is order-one, so the neglected higher harmonics and any breakdown of the fast-angle averaging can easily flatten or erase the double-well structure that the leading-order ε1, ε2 terms are supposed to create. The abstract states the bistability result but supplies no error estimates, no comparison to the unaveraged equations, and no numerical integrations that would show whether the threshold survives in the full system. Without those checks the central claim stays provisional. This is the kind of analytical framework that people modeling circumbinary planets or tight S-type systems might want to build on, provided the truncation holds. It is coherent enough and grounded enough in the existing CR3BP literature to deserve a serious referee, but any review should insist on at least one set of direct numerical tests near the quoted threshold before the bistability prediction is taken as reliable. I would not cite it yet and would not bring it to a reading group until the validation appears.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a unified Hamiltonian formulation for mean-motion resonances in the CR3BP applied to binary star systems with planets in S-type and P-type orbits. It starts from the full CR3BP Hamiltonian, performs canonical transformations to resonant action-angle variables, and derives reduced 1DOF Hamiltonians via systematic averaging over fast orbital motion. Leading-order scaling laws for the resonant Fourier coefficients ε1 and ε2 (depending on mass ratio μ and a/ab) are obtained. The resulting effective potential is shown to exhibit bistability when |ε2/ε1| > 1/4, which sets the stage for stochastic resonance. Scaling laws are derived separately for S-type and P-type cases, and estimates for known binary-planet systems are provided, indicating that observed systems lie below threshold while extreme configurations (e.g., a/ab ≲ 1.5 for P-type, near-equal masses) could reach bistability.

Significance. If the averaging and truncation approximations hold, the work supplies a useful analytical framework for resonant dynamics in binaries that incorporates finite perturber mass, unlike standard Solar-System treatments. The explicit bistability threshold |ε2/ε1| > 1/4 and the derived scaling laws constitute concrete, falsifiable predictions that could guide future observations or simulations. The unified treatment of S-type and P-type configurations is a clear strength.

major comments (2)

[Results / derivation of effective potential] The bistability result (Results section) is obtained from the averaged 1DOF effective Hamiltonian after canonical transformation. For P-type configurations with μ ≈ 0.5 and a/ab ≲ 1.5 the disturbing function is order-1, so the neglected higher harmonics and non-resonant secular terms could modify or destroy the double-well structure near the |ε2/ε1| = 1/4 threshold. The manuscript provides no error estimates on the averaging or direct numerical comparisons between the reduced Hamiltonian and the full CR3BP equations of motion for these regimes.
[Methods / scaling laws for Fourier coefficients] The leading-order Fourier scaling laws for ε1 and ε2 (Methods) are used to evaluate the bistability condition. The paper does not quantify the magnitude of the next-order corrections to these coefficients or demonstrate that they remain small enough to preserve the sign of the ratio near threshold for finite-μ binaries.

minor comments (2)

[Estimates for known systems] A table summarizing the computed |ε2/ε1| values and mass-ratio/a/ab parameters for the estimated known systems would improve clarity and allow readers to reproduce the threshold comparisons.
[Canonical transformations] Notation for the resonant angles and actions after the canonical transformation could be introduced with an explicit table or diagram to aid readers unfamiliar with resonant Hamiltonian reductions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments and for recognizing the value of the unified Hamiltonian framework. We respond to each major comment below.

read point-by-point responses

Referee: [Results / derivation of effective potential] The bistability result (Results section) is obtained from the averaged 1DOF effective Hamiltonian after canonical transformation. For P-type configurations with μ ≈ 0.5 and a/ab ≲ 1.5 the disturbing function is order-1, so the neglected higher harmonics and non-resonant secular terms could modify or destroy the double-well structure near the |ε2/ε1| = 1/4 threshold. The manuscript provides no error estimates on the averaging or direct numerical comparisons between the reduced Hamiltonian and the full CR3BP equations of motion for these regimes.

Authors: We agree that for P-type configurations with μ ≈ 0.5 and a/ab ≲ 1.5 the disturbing function reaches order-1, so higher harmonics and secular terms omitted by the averaging could alter the effective potential near the |ε₂/ε₁| = 1/4 threshold. The bistability condition is derived strictly within the leading-order averaged 1DOF Hamiltonian and therefore constitutes a necessary condition inside that approximation. The submitted manuscript contains neither quantitative error estimates for the averaging nor direct numerical integrations against the full CR3BP equations. In the revised version we will add a subsection that supplies order-of-magnitude bounds on the neglected terms and presents a limited set of numerical comparisons for representative extreme cases to test whether the double-well structure survives. revision: yes
Referee: [Methods / scaling laws for Fourier coefficients] The leading-order Fourier scaling laws for ε1 and ε2 (Methods) are used to evaluate the bistability condition. The paper does not quantify the magnitude of the next-order corrections to these coefficients or demonstrate that they remain small enough to preserve the sign of the ratio near threshold for finite-μ binaries.

Authors: The scaling laws for ε₁ and ε₂ are obtained at leading order from the Fourier expansion of the disturbing function. The manuscript does not evaluate the size of the next-order corrections. For the observed systems we examine, |ε₂/ε₁| lies sufficiently far below 1/4 that moderate corrections are unlikely to change the sign. For the extreme configurations that approach the threshold, such corrections could in principle affect the ratio. We will revise the Methods section to include rough estimates of the relative magnitude of the next-order terms and to state explicitly the regime in which the leading-order sign of the ratio is expected to remain reliable. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the Hamiltonian reduction and bistability derivation

full rationale

The paper starts from the standard CR3BP Hamiltonian, applies canonical transformations to resonant action-angle variables, and performs systematic averaging over fast angles to obtain a reduced 1DOF effective Hamiltonian. The bistability condition |ε2/ε1| > 1/4 follows directly as a mathematical property of the resulting effective potential containing the first two resonant harmonics; the scaling laws for the Fourier coefficients ε1 and ε2 are derived from leading-order perturbation theory in the mass ratio and orbital distance. No self-definitional loops, fitted parameters renamed as predictions, or load-bearing self-citations appear in the derivation chain. The central claims remain independent of the paper's own inputs and are self-contained within the averaged model.

Axiom & Free-Parameter Ledger

0 free parameters · 3 axioms · 0 invented entities

The model rests on standard CR3BP assumptions and averaging approximations without introducing new free parameters or postulated entities.

axioms (3)

domain assumption The system obeys the Circular Restricted Three-Body Problem with a planet of negligible mass orbiting two finite-mass stars.
Explicitly stated as the foundational framework in the context and methods.
standard math Canonical transformations to resonant action-angle variables preserve the Hamiltonian structure.
Standard technique invoked to reach the reduced one-degree-of-freedom model.
domain assumption Averaging over fast orbital motion isolates the slow resonant dynamics without loss of essential behavior.
Core step used to obtain the effective potential and bistability condition.

pith-pipeline@v0.9.0 · 5610 in / 1571 out tokens · 111810 ms · 2026-05-08T15:23:04.788857+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

16 extracted references

[1]

1982, Tellu s, 34, 10

Benzi, R., Parisi, G., Sutera, A., & Vulpiani, A. 1982, Tellu s, 34, 10

1982
[2]

1981, Journal of Physic s A: Mathematical and General, 14, L453

Benzi, R., Sutera, A., & Vulpiani, A. 1981, Journal of Physic s A: Mathematical and General, 14, L453

1981
[3]

R., Carter, J

Doyle, L. R., Carter, J. A., Fabrycky, D. C., et al. 2011, Scie nce, 333, 1602

2011
[4]

Freitas, U. S. & Ozorio de Almeida, A. M. 2005, Physica D: Nonl inear Phenom- ena, 202, 129

2005
[5]

1998,Reviews of Modern Physics, 70, 223

Gammaitoni, L., Hänggi, P ., Jung, P ., & Marchesoni, F. 1998,Reviews of Modern Physics, 70, 223

1998
[6]

& Lemaitre, A

Henrard, J. & Lemaitre, A. 1983, Celestial Mechanics, 30, 19 7

1983
[7]

B., Orosz, J

Kostov, V . B., Orosz, J. A., Feinstein, A. D., et al. 2020, The Astronomical Jour- nal, 159, 253

2020
[8]

& Boué, G

Laskar, J. & Boué, G. 2010, Astronomy and Astrophysics, 522, A60

2010
[9]

Martin, D. V . 2019, Living Reviews in Computational As- trophysics, 5, 1, regularly updated online version at http://www.astro.twam.info/cbp-catalogue/

2019
[10]

2002, Nature, 419, 47

Morbidelli, A. 2002, Nature, 419, 47

2002
[11]

Murray, C. D. & Dermott, S. F. 1999, Solar System Dynamics (Cambridge: Cam- bridge University Press), ﬁrst edition, reprinted with cor rections 2000, 2008 NASA Exoplanet Archive. 2024, NASA Exoplanet Archive, acce ssed: 2024-12- 06

1999
[12]

2018, in Handbo ok of Exoplan- ets, ed

Schwarz, R., Bazsó, Á., Funk, B., & Dvorak, R. 2018, in Handbo ok of Exoplan- ets, ed. H. J. Deeg & J. A. Belmonte (Cham: Springer Internati onal Publish- ing), 1–25

2018
[13]

2021, Space Sci ence Reviews, 217, 43

Schwarz, R., Bazsó, Á., Funk, B., & Dvorak, R. 2021, Space Sci ence Reviews, 217, 43

2021
[14]

2016, A&A, 594 , A39

Schwarz, R., Bazsó, Á., Funk, B., & Zechner, R. 2016, A&A, 594 , A39

2016
[15]

Shevchenko, I. I. 2011, Solar System Research, 45, 90 The TESS Science Team & The Kepler Science Team. 2023, Astron omy & As- trophysics Review, 31, 5, compilation of results from Keple r, TESS, and RV surveys through 2023

2011
[16]

Winter, O. C. & Murray, C. D. 1997, Celestial Mechanics and Dy namical As- tronomy, 66, 1 Article number, page 8

1997