Effective stability estimates close to resonances with applications to rotational dynamics
Pith reviewed 2026-05-13 18:30 UTC · model grok-4.3
The pith
Nekhoroshev-like estimates with optimized parameters and reduced perturbations provide effective stability for orbits near resonances in spin-orbit models.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In nearly-integrable Hamiltonian systems defined over non-resonant domains, Nekhoroshev-like estimates are applied in the neighborhood of resonances to bound the variation of the action variables over long times. An optimization algorithm selects parameters to maximize the stability time, and perturbation theory reduces the norm of the perturbing function. When implemented for the one-dimensional spin-orbit problem and the two-dimensional spin-spin-orbit model, this yields stability results for sequences of Diophantine frequencies converging to the main resonances of these rotational dynamics models.
What carries the argument
Nekhoroshev-like estimates applied near resonances, with parameter optimization and perturbation theory to minimize the perturbing function norm.
If this is right
- Explicit long-term bounds on action drift are obtained for the spin-orbit problem near its main resonances.
- The two-dimensional spin-spin-orbit model admits similar stability estimates close to its resonances.
- Stability holds for Diophantine frequencies arbitrarily close to resonant ones after optimization.
- The freedom in parameter choice allows tailoring the estimates to specific models.
Where Pith is reading between the lines
- The technique could be extended to other nearly-integrable systems in celestial mechanics, such as those with more degrees of freedom.
- Further reductions in perturbation norms might be achieved by higher-order perturbation theory, leading to even longer stability times.
- These bounds suggest that numerical explorations of rotational dynamics should focus on verifying stability within the predicted intervals.
- Applications might include assessing the long-term behavior of planetary spins in multi-body systems.
Load-bearing premise
The systems are nearly-integrable and the Nekhoroshev-like estimates apply in the resonance neighborhoods with parameters that can be freely chosen and optimized.
What would settle it
Direct numerical simulation of the spin-orbit Hamiltonian showing that the action variables change by more than the estimated bound within the computed stability time for a frequency close to resonance would falsify the result.
Figures
read the original abstract
We consider nearly-integrable Hamiltonian systems defined over a non-resonant domain. In the neighborhood of resonances, we use Nekhoroshev-like estimates to provide effective stability bounds for the action variables over long time. The applicability conditions of these estimates allow some freedom in the choice of parameters. Hence, we develop an optimization algorithm for choosing parameters that maximize the stability time. To further improve the stability estimates, we use perturbation theory to reduce the norm of the perturbing function. We implement this procedure (effective stability estimates and perturbation theory) to analyze the stability of sequences of irrational (Diophantine) frequencies converging to frequencies corresponding to resonances. We consider two applications to models describing problems of rotational dynamics in Celestial Mechanics: the spin-orbit problem, described by a 1D time-dependent Hamiltonian, and the spin-spin-orbit model, described by a 2D time-dependent Hamiltonian. We show stability results for orbits close to the main resonances associated with such models.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops effective stability estimates for action variables in nearly-integrable Hamiltonian systems near resonances by applying Nekhoroshev-like theorems, optimizing free parameters (such as ε and τ) to maximize stability times, and using perturbation theory to reduce the perturbing norm. The procedure is implemented for sequences of Diophantine frequencies approaching resonant ones and applied to the 1D spin-orbit and 2D spin-spin-orbit models, yielding stability results for orbits close to the main resonances in these celestial mechanics problems.
Significance. If the post-optimization validity conditions hold, the work supplies concrete, computable stability times for rotational dynamics near resonances, which is of practical value in celestial mechanics. The combination of parameter optimization with a preliminary perturbation reduction is a methodological strength that can improve upon direct Nekhoroshev application; the explicit algorithmic treatment and concrete model applications add to the utility.
major comments (3)
- [§4] §4 (Optimization procedure): The optimization algorithm for Nekhoroshev parameters is described at a high level and stated to maximize stability time, but no a-posteriori verification is supplied that the selected values keep the reduced perturbation norm below the smallness threshold required by the cited Nekhoroshev-like estimates (relative to resonance distance and Diophantine constants). This is load-bearing for the central claim.
- [§5] §5 (Spin-orbit application): The stability bounds for Diophantine frequency sequences converging to the main resonances are presented after optimization, yet the manuscript does not report the concrete optimized parameter values nor confirm they remain inside the domain of validity of the underlying theorem; without this, the reported stability times cannot be directly justified.
- [§3.1–3.2] §3.1–3.2 (Nekhoroshev estimates near resonances): The claim that applicability conditions 'allow some freedom' in parameter choice is used to justify optimization, but the text does not quantify how the optimization step alters the effective perturbation size or provide an explicit check against the theorem hypotheses after reduction; this leaves a moderate derivation gap between the abstract statement and the final bounds.
minor comments (3)
- [Notation] The notation for Diophantine constants, effective perturbation norms, and stability times is introduced piecemeal; a consolidated table of symbols with references to their first appearance would improve readability.
- [Figures] Figure captions for the frequency sequences and stability plots lack explicit indication of the optimized parameter values used to generate each curve.
- [§2] A few references to standard Nekhoroshev statements (e.g., the precise form of the stability time exponent) could be added for completeness in §2.
Simulated Author's Rebuttal
We thank the referee for the thorough reading and valuable comments, which help clarify the presentation of our optimization procedure and its verification. We address each major comment below and will revise the manuscript to incorporate the requested details and checks.
read point-by-point responses
-
Referee: [§4] §4 (Optimization procedure): The optimization algorithm for Nekhoroshev parameters is described at a high level and stated to maximize stability time, but no a-posteriori verification is supplied that the selected values keep the reduced perturbation norm below the smallness threshold required by the cited Nekhoroshev-like estimates (relative to resonance distance and Diophantine constants). This is load-bearing for the central claim.
Authors: We agree that an explicit a-posteriori verification is essential for the central claims. In the revised manuscript we will add a dedicated subsection (or appendix) that computes the reduced perturbation norm after optimization for each frequency sequence and directly compares it to the smallness threshold of the Nekhoroshev-like theorem, taking into account the resonance distance and the specific Diophantine constants. This verification will be performed numerically for the sequences considered in the applications. revision: yes
-
Referee: [§5] §5 (Spin-orbit application): The stability bounds for Diophantine frequency sequences converging to the main resonances are presented after optimization, yet the manuscript does not report the concrete optimized parameter values nor confirm they remain inside the domain of validity of the underlying theorem; without this, the reported stability times cannot be directly justified.
Authors: We will include a table (or set of tables) in the revised §5 that lists the concrete optimized values of the free parameters (ε, τ and any auxiliary constants) for each Diophantine frequency in the sequences approaching the main resonances. For each entry we will also report the corresponding reduced perturbation norm and explicitly confirm that it lies below the theorem threshold, thereby justifying the reported stability times. revision: yes
-
Referee: [§3.1–3.2] §3.1–3.2 (Nekhoroshev estimates near resonances): The claim that applicability conditions 'allow some freedom' in parameter choice is used to justify optimization, but the text does not quantify how the optimization step alters the effective perturbation size or provide an explicit check against the theorem hypotheses after reduction; this leaves a moderate derivation gap between the abstract statement and the final bounds.
Authors: We will expand §§3.1–3.2 to quantify the effect of optimization on the effective perturbation size (e.g., by tabulating or plotting the norm before and after the optimization step) and to supply explicit post-reduction checks against all hypotheses of the Nekhoroshev-like estimates. These additions will close the derivation gap and make the transition from the abstract theorem to the concrete bounds fully transparent. revision: yes
Circularity Check
Stability estimates derived from external Nekhoroshev theorems via algorithmic parameter optimization
full rationale
The paper applies Nekhoroshev-like estimates to resonant neighborhoods after an algorithmic optimization of free parameters (chosen to maximize the resulting stability time) and a preliminary perturbation-theory reduction of the perturbing norm. The stability bounds are direct consequences of the cited external theorems evaluated at the optimized values, without any redefinition of the target stability time in terms of itself, without fitting parameters to the final result, and without load-bearing self-citations that reduce the central claim to an unverified premise. The derivation chain remains self-contained against the stated applicability conditions of the theorems.
Axiom & Free-Parameter Ledger
free parameters (1)
- Nekhoroshev parameters (epsilon, tau, etc.)
axioms (1)
- domain assumption The Hamiltonian is nearly integrable and the frequency vector satisfies the non-resonance conditions outside a small neighborhood of the target resonances.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
we develop an optimization algorithm for choosing parameters that maximize the stability time... sequences of irrational (Diophantine) frequencies... Γ(k1,k2)z,s = −k2/k1 − s/(z+γ) with γ golden ratio
-
IndisputableMonolith/Foundation/BranchSelection.leanbranch_selection unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Theorem 1... E ≤ 1/(27ℓ) α r / K ... T ≡ s0 r / (5 E) exp(K s0 / 6)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
V. I. Arnol’d. Proof of a theorem of A. N. Kolmogorov on the invariance of quasi-periodic motions under small perturbations. Russian Math. Surveys, 18(5):9–36, 1963
work page 1963
-
[2]
G. Bou´ e. The two rigid body interaction using angular momentum theory formulae. Celestial Mechanics and Dynamical Astronomy, 128(2-3):261–273, June 2017
work page 2017
-
[3]
A. Bounemoura and J.-P. Marco. Improved exponential stability for near-integrable quasi-convex hamiltonians. Nonlinearity, 24(1):97, nov 2010
work page 2010
-
[4]
A. P. Bustamante, A. Celletti, and C. Lhotka. The dynamics of the spin–spin problem in Celestial Mechanics. Communications in Nonlinear Science and Numerical Simulations, 142:108548, Mar. 2025
work page 2025
-
[5]
J. Cassels. An introduction to Diophantine approximation. Cambridge University Press, 1957
work page 1957
- [6]
- [7]
-
[8]
A. Celletti, I. De Blasi, and C. Efthymiopoulos. Nekhoroshev estimates for the orbital stability of Earth’s satellites. Celestial Mechanics and Dynamical Astronomy, 135(2):10, Apr. 2023
work page 2023
-
[9]
A. Celletti, C. Falcolini, and U. Locatelli. On the break-down threshold of invariant tori in four dimensional maps. Regul. Chaotic Dyn., 9(3):227–253, 2004
work page 2004
-
[10]
A. Celletti and L. Ferrara. An Application of the Nekhoroshev Theorem to the Restricted Three- Body Problem. Celestial Mechanics and Dynamical Astronomy, 64(3):261–272, Sept. 1996
work page 1996
-
[11]
A. Celletti, J. Gimeno, and M. Misquero. The Spin-Spin Problem in Celestial Mechanics. Journal of NonLinear Science, 32(6):88, Dec. 2022
work page 2022
-
[12]
A. Celletti and A. Giorgilli. On the stability of the Lagrangian points in the spatial restricted problem of three bodies. Celestial Mechanics and Dynamical Astronomy, 50(1):31–58, Mar. 1990
work page 1990
- [13]
-
[14]
C. Efthymiopoulos. Canonical perturbation theory; stability and diffusion in Hamiltonian systems: applications in dynamical astronomy. Workshop Series of the Asociacion Argentina de Astronomia, 3:3–146, Jan. 2011. 38 A. CELLETTI, A. DOGKAS, AND A. FRANCESCA GUIDO
work page 2011
-
[15]
C. Efthymiopoulos, A. Giorgilli, and G. Contopoulos. Nonconvergence of formal integrals: II. Im- proved estimates for the optimal order of truncation. Journal of Physics A Mathematical General, 37(45):10831–10858, Nov. 2004
work page 2004
-
[16]
S. Ferraz-Mello. Canonical Perturbation Theories - Degenerate Systems and Resonance, volume
-
[17]
A. Giorgilli and C. Skokos. On the stability of the Trojan asteroids. Astron. & Astroph., 317:254– 261, Jan. 1997
work page 1997
-
[18]
P. Goldreich and S. Peale. Spin-orbit coupling in the solar system. Astronomical Journal, Vol. 71, p. 425 (1966), 71:425, 1966
work page 1966
-
[19]
J. M. Greene. A method for determining a stochastic transition. Jour. Math. Phys., 20:1183–1201, 1979
work page 1979
- [20]
- [21]
- [22]
-
[23]
A. N. Kolmogorov. On conservation of conditionally periodic motions for a small change in Hamil- ton’s function. Dokl. Akad. Nauk SSSR (N.S.), 98:527–530, 1954. English translation inStochastic Behavior in Classical and Quantum Hamiltonian Systems (Volta Memorial Conf., Como, 1977), Lecture Notes in Phys., 93, pages 51–56. Springer, Berlin, 1979
work page 1954
-
[24]
H. Lei. Spin–orbit coupling of the primary body in a binary asteroid system. Celestial Mechanics and Dynamical Astronomy, 136(5):37, Oct. 2024
work page 2024
-
[25]
P. Lochak. Canonical perturbation theory via simultaneous approximation. Russian Mathematical Surveys, 47(6):57, dec 1992
work page 1992
-
[26]
P. Lochak and A. I. Neishtadt. Estimates of stability time for nearly integrable systems with a quasiconvex Hamiltonian. Chaos, 2(4):495–499, Oct. 1992
work page 1992
-
[27]
A. J. Maciejewski. Reduction, Relative Equilibria and Potential in the Two Rigid Bodies Problem. Celestial Mechanics and Dynamical Astronomy, 63(1):1–28, Mar. 1995
work page 1995
- [28]
-
[29]
J. Moser. On invariant curves of area-preserving mappings of an annulus. Nachr. Akad. Wiss. G¨ ottingenMath.-Phys. Kl. II, 1962:1–20, 1962
work page 1962
-
[30]
N. N. Nekhoroshev. An exponential estimate of the time of stability of nearly integrable Hamil- tonian systems. Uspehi Mat. Nauk, 32(6(198)):5–66, 287, 1977. English translation: Russian Math. Surveys, 32(6):1–65, 1977
work page 1977
-
[31]
G. Pinzari, B. Scoppola, and M. Veglianti. Spin orbit resonance cascade via core shell model: application to Mercury and Ganymede. Celestial Mechanics and Dynamical Astronomy, 136(5):39, Oct. 2024
work page 2024
- [32]
-
[33]
E. T. Whittaker and G. N. Watson. A Course of Modern Analysis. Cambridge University Press, 5 edition, 2021
work page 2021
- [34]
-
[35]
J. Zhang and K. Zhang. Improved stability for analytic quasi-convex nearly integrable systems and optimal speed of arnold diffusion. Nonlinearity, 30(7):2918, jun 2017
work page 2017
-
[36]
K. Zhang. Speed of arnold diffusion for analytic hamiltonian systems. Inventiones Mathematicae, 186:255–290, 11 2011. EFFECTIVE STABILITY ESTIMATES CLOSE TO RESONANCES 39 Department of Mathematics, University of Roma Tor Vergata, Via della Ricerca Scientifica 1, 00133 Roma (Italy) Email address:celletti@mat.uniroma2.it Department of Mathematics, Universit...
work page 2011
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.