A new expansion of planetary disturbing function and applications to interior, co-orbital and exterior resonances with planets

In this study, a new expansion of planetary disturbing function is developed for describing the resonant dynamics of minor bodies with arbitrary inclinations and semimajor axis ratios. In practice, the disturbing function is expanded around circular orbits in the first step and then, in the second step, the resulting mutual interaction between circular orbits is expanded around a reference point. As usual, the resulting expansion is presented in the Fourier series form, where the force amplitudes are dependent on the semimajor axis, eccentricity and inclination, and the harmonic arguments are linear combinations of the mean longitude, longitude of pericenter and longitude of ascending node of each mass. The resulting new expansion is valid for arbitrary inclinations and semimajor axis ratios. In the case of mean motion resonant configuration, the disturbing function can be easily averaged to produce the analytical expansion of resonant disturbing function. Based on the analytical expansion, the Hamiltonian model of mean motion resonances is formulated, and the resulting analytical developments are applied to Jupiter’s inner and co-orbital resonances and Neptune’s exterior resonances. Analytical expansion is validated by comparing the analytical results with the associated numerical outcomes.

Period Ratio Sculpting near Second-order Mean-motion Resonances

Second-order mean-motion resonances lead to an interesting phenomenon in the sculpting of the period-ratio distribution, due to their shape and width in period-ratio/eccentricity space. As the osculating periods librate in resonance, the time-averaged period ratio approaches the exact commensurability. The width of second-order resonances increases with increasing eccentricity, and thus more eccentric systems have a stronger peak at commensurability when averaged over sufficient time. The libration period is short enough that this time-averaging behavior is expected to appear on the timescale of the Kepler mission. Using N-body integrations of simulated planet pairs near the 5:3 and 3:1 mean-motion resonances, we investigate the eccentricity distribution consistent with the planet pairs observed by Kepler. This analysis, an approach independent from previous studies, shows no statistically significant peak at the 3:1 resonance and a small peak at the 5:3 resonance, placing an upper limit on the Rayleigh scale parameter, σ, of the eccentricity of the observed Kepler planets at σ = 0.245 (3:1) and σ = 0.095 (5:3) at 95% confidence, consistent with previous results from other methods.

A Pair of Warm Giant Planets near the 2:1 Mean Motion Resonance around the K-dwarf Star TOI-2202*

TOI-2202 b is a transiting warm Jovian-mass planet with an orbital period of P = 11.91 days identified from the Full Frame Images data of five different sectors of the TESS mission. Ten TESS transits of TOI-2202 b combined with three follow-up light curves obtained with the CHAT robotic telescope show strong transit timing variations (TTVs) with an amplitude of about 1.2 hr. Radial velocity follow-up with FEROS, HARPS, and PFS confirms the planetary nature of the transiting candidate (a b = 0.096 ± 0.001 au, m b = 0.98 ± 0.06 M Jup), and a dynamical analysis of RVs, transit data, and TTVs points to an outer Saturn-mass companion (a c = 0.155 ± 0.002 au, m c = 0.37 ± 0.10 M Jup) near the 2:1 mean motion resonance. Our stellar modeling indicates that TOI-2202 is an early K-type star with a mass of 0.82 M ⊙, a radius of 0.79 R ⊙, and solar-like metallicity. The TOI-2202 system is very interesting because of the two warm Jovian-mass planets near the 2:1 mean motion resonance, which is a rare configuration, and their formation and dynamical evolution are still not well understood.

The functional relation between three-body mean motion resonances and Yarkovsky drift speeds

ABSTRACT We examined the motion of asteroids across the three-body mean motion resonances (MMRs) with Jupiter and Saturn and with the Yarkovsky drift speed in the semimajor axis of the asteroids. The research was conducted using numerical integrations performed using the Orbit9 integrator with 84 000 test asteroids. We calculated time delays, dtr, caused by the seven three-body MMRs on the mobility of test asteroids with 10 positive and 10 negative Yarkovsky drift speeds, which are reliable for Main Belt asteroids. Our final results considered only test asteroids that successfully crossed over the MMRs without close approaches to the planets. We have devised two equations that approximately describe the functional relation between the average time 〈dtr〉 spent in the resonance, the strength of the resonance SR, and the semimajor axis drift speed da/dt (positive and negative) with the orbital eccentricities of asteroids in the range (0, 0.1). Comparing the values of 〈dtr〉 obtained from the numerical integrations and from the derived functional relations, we analysed average values of 〈dtr〉 in all three-body MMRs for every da/dt. The main conclusion is that the analytical and numerical estimates of the average time 〈dtr〉 are in very good agreement, for both positive and negative da/dt. Finally, this study shows that the functional relation we obtain for three-body MMRs is analogous to that previously obtained for two-body MMRs.

An integrable model for first-order three-planet mean motion resonances

Recent works on three-planet mean motion resonances (MMRs) have highlighted their importance for understanding the details of the dynamics of planet formation and evolution. While the dynamics of two-planet MMRs are well understood and approximately described by a one-degree-of-freedom Hamiltonian, little is known of the exact dynamics of three-body resonances besides the cases of zeroth-order MMRs or when one of the bodies is a test particle. In this work, I propose the first general integrable model for first-order three-planet mean motion resonances. I show that one can generalize the strategy proposed in the two-planet case to obtain a one-degree-of-freedom Hamiltonian. The dynamics of these resonances are governed by the second fundamental model of resonance. The model is valid for any mass ratio between the planets and for every first-order resonance. I show the agreement of the analytical model with numerical simulations. As examples of application, I show how this model could improve our understanding of the capture into MMRs as well as their role in the stability of planetary systems.

Relationship between Three-Dimensional Radiation Stress and Vortex-Force Representations

In numerical ocean models, the effect of waves on currents is usually expressed by either vortex-force or radiation stress representations. In this paper, the differences and similarities between those two representations are investigated in detail in conditions of both conservative and nonconservative waves. In addition, comparisons between different sets of equations of mean motion that apply different representations of wave-induced forcing terms are included. The comparisons are useful for selecting a suitable numerical ocean model to simulate the mean current in conditions of waves combined with currents.

Revisiting The Averaged Problem in The Case of Mean-Motion Resonances of The Restricted Three-Body Problem. Global Rigorous Treatment and Application To The Co-Orbital Motion.

A classical approach to the restricted three-body problem is to analyze the dynamics of the massless body in the synodic reference frame. A different approach is represented by the perturbative treatment: in particular the averaged problem of a mean-motion resonance allows to investigate the long-term behavior of the solutions through a suitable approximation that focuses on a particular region of the phase space. In this paper, we intend to bridge a gap between the two approaches in the specific case of mean-motion resonant dynamics, establish the limit of validity of the averaged problem, and take advantage of its results in order to compute trajectories in the synodic reference frame. After the description of each approach, we develop a rigorous treatment of the averaging process, estimate the size of the transformation and prove that the averaged problem is a suitable approximation of the restricted three-body problem as long as the solutions are located outside the Hill's sphere of the secondary. In such a case, a rigorous theorem of stability over finite but large timescales can be proven. We establish that a solution of the averaged problem provides an accurate approximation of the trajectories on the synodic reference frame within a finite time that depend on the minimal distance to the Hill's sphere of the secondary. The last part of this work is devoted to the co-orbital motion (i.e., the dynamics in 1:1 mean-motion resonance) in the circular-planar case. In this case, an interpretation of the solutions of the averaged problem in the synodic reference frame is detailed and a method that allows to compute co-orbital trajectories is displayed.