No Arabic abstract
We present N-body simulations of resonant planets with inclined orbits that show chaotically evolving eccentricities and inclinations that can persist for at least 10 Gyr. A wide range of behavior is possible, from fast, low amplitude variations to systems in which eccentricities reach 0.9999 and inclinations 179.9 degrees. While the orbital elements evolve chaotically, at least one resonant argument always librates. We show that the HD 73526, HD 45364 and HD 60532 systems may be in chaotically-evolving resonances. Chaotic evolution is apparent in the 2:1, 3:1 and 3:2 resonances, and for planetary masses from lunar- to Jupiter-mass. In some cases, orbital disruption occurs after several Gyr, implying the mechanism is not rigorously stable, just long-lived relative to the main sequence lifetimes of solar-type stars. Planet-planet scattering appears to yield planets in inclined resonances that evolve chaotically in about 0.5% of cases. These results suggest that 1) approximate methods for identifying unstable orbital architectures may have limited applicability, 2) the observed close-in exoplanets may be produced during the high eccentricity phases induced by inclined resonances, 3) those exoplanets orbital planes may be misaligned with the host stars spin axis, 4) systems with resonances may be systematically younger than those without, 5) the distribution of period ratios of adjacent planets detected via transit may be skewed due to inclined resonances, and 6) potentially habitable planets in resonances may have dramatically different climatic evolution than the Earth. The GAIA spacecraft is capable of discovering giant planets in these types of orbits.
We present preliminary though statistically significant evidence that shows that multiplanetary systems that exhibit a 2/1 period commensurability are in general younger than multiplanetary systems without commensurabilities, or even systems with other commensurabilities. An immediate possible conclusion is that the 2/1 mean-motion resonance in planetary systems, tends to be disrupted after typically a few Gyrs.
In circumstellar discs, collisional grinding of planetesimals produces second-generation dust. While it remains unclear whether this ever becomes a major component of the total dust content, the presence of such dust, and potentially the substructure within, it can be used to explore a discs physical conditions. A perturbing planet produces nonaxisymmetric structures and gaps in the dust, regardless of its origin. The dynamics of planetesimals, however, will be very different than that of small dust grains due to weaker gas interactions. Therefore, planetesimal collisions could create dusty disc structures that would not exist otherwise. In this work, we use N-body simulations to investigate the collision rate profile of planetesimals near mean-motion resonances. We find that a distinct bump or dip feature is produced in the collision profile, the presence of which depends on the libration width of the resonance and the separation between the peri- and apocenter distances of the edges of the resonance. The presence of one of these two features depends on the mass and eccentricity of the planet. Assuming that the radial dust emission traces the planetesimal collision profile, the presence of a bump or dip feature in the dust emission at the 2:1 mean-motion resonance can constrain the orbital properties of the perturbing planet. This assumption is valid, so long as radial drift does not play a significant role during the collisional cascade process. Under this assumption, these features in the dust emission should be marginally observable in nearby protoplanetary disks with ALMA.
The identification of mean motion resonances in exoplanetary systems or in the Solar System might be cumbersome when several planets and large number of smaller bodies are to be considered. Based on the geometrical meaning of the resonance variable, an efficient method is introduced and described here, by which mean motion resonances can be easily find without any a priori knowledge on them. The efficiency of this method is clearly demonstrated by using known exoplanets engaged in mean motion resonances, and also some members of different families of asteroids and Kuiper-belt objects being in mean motion resonances with Jupiter and Neptune respectively.
As the discoveries of more minor bodies in retrograde resonances with giant planets, such as 2015 BZ509 and 2006 RJ2, our curiosity about the Kozai-Lidov dynamics inside the retrograde resonance has been sparked. In this study, we focus on the 3D retrograde resonance problem and investigate how the resonant dynamics of a minor body impacts on its own Kozai-Lidov cycle. Firstly we deduce the action-angle variables and canonical transformations that deal with the retrograde orbit specifically. After obtaining the dominant Hamiltonian of this problem, we then carry out the numerical averaging process in closed form to generate phase-space portraits on a $e-omega$ space. The retrograde 1:1 resonance is particularly scrutinized in detail, and numerical results from a CRTBP model shows a great agreement with the our semi-analytical portraits. On this basis, we inspect two real minor bodies currently trapped in retrograde 1:1 mean motion resonance. It is shown that they have different Kozai-Lidov states, which can be used to analyze the stability of their unique resonances. In the end, we further inspect the Kozai-Lidov dynamics inside the 2:1 and 2:5 retrograde resonance, and find distinct dynamical bifurcations of equilibrium points on phase-space portraits.
Exoplanet systems with multiple planets in mean motion resonances have often been hailed as a signpost of disk driven migration. Resonant chains like Kepler-223 and Kepler-80 consist of a trio of planets with the three-body resonant angle librating and/or with a two-body resonant angle librating for each pair. Here we investigate whether close-in super-Earths and mini-Neptunes forming in situ can lock into resonant chains due to dissipation from a depleted gas disk. We simulate the giant impact phase of planet formation, including eccentricity damping from a gaseous disk, followed by subsequent dynamical evolution over tens of millions of years. In a fraction of simulated systems, we find that planets naturally lock into resonant chains. These planets achieve a chain of near-integer period ratios during the gas disk stage, experience eccentricity damping that captures them into resonance, stay in resonance as the gas disk dissipates, and avoid subsequent giant impacts, eccentricity excitation, and chaotic diffusion that would dislodge the planets from resonance. Disk conditions that enable planets to complete their formation during the gas disk stage enable those planets to achieve tight period ratios <= 2 and, if they happen to be near integer period ratios, lock into resonance. Using the weighting of different disk conditions deduced by MacDonald et al. (2020) and forward modeling Kepler selection effects, we find that our simulations of in situ formation via oligarchic growth lead to a rate of observable trios with integer period ratios and librating resonant angles comparable to observed Kepler systems.