No Arabic abstract
An episode of dynamical instability is thought to have sculpted the orbital structure of the outer solar system. When modeling this instability, a key constraint comes from Jupiters fifth eccentric mode (quantified by its amplitude M55), which is an important driver of the solar systems secular evolution. Starting from commonly-assumed near-circular orbits, the present-day giant planets architecture lies at the limit of numerically generated systems, and M55 is rarely excited to its true value. Here we perform a dynamical analysis of a large batch of artificially triggered instabilities, and test a variety of configurations for the giant planets primordial orbits. In addition to more standard setups, and motivated by the results of modern hydrodynamical simulations of the giant planets evolution within the primordial gaseous disk, we consider the possibility that Jupiter and Saturn emerged from the nebular gas locked in 2:1 resonance with non-zero eccentricities. We show that, in such a scenario, the modern Jupiter-Saturn system represents a typical simulation outcome, and M55 is commonly matched. Furthermore, we show that Uranus and Neptunes final orbits are determined by a combination of the mass in the primordial Kuiper belt and that of an ejected ice giant.
We find an interesting fact that fictitious retrograde co-orbitals of Saturn, or small bodies inside the retrograde 1:1 resonance with Saturn, are highly unstable in our numerical simulations. It is shown that in the presence of Jupiter, the retrograde co-orbitals will get ejected from Saturns co-orbital space within a timescale of 10 Myr. This scenario reminds us of the instability of Saturn Trojans caused by both the Great Inequality and the secular resonances. Therefore, we carry out in-depth inspections on both mechanisms and prove that the retrograde resonance overlap, raised by Great Inequality, cannot serve as an explanation for the instability of retrograde co-orbitals, due to the weakness of the retrograde 2:5 resonance with Jupiter at a low eccentricity. However, we discover that both $ u_5$ and $ u_6$ secular resonances contribute to the slow growth of the eccentricity, therefore, are possibly the primary causes of the instability inside Saturns retrograde co-orbital space.
Using astrometric observations spanning more than a century and including a large set of Cassini data, we determine Saturns tidal parameters through their current effects on the orbits of the eight main and four coorbital moons. We have used the latter to make the first determination of Saturns Love number, $k_2=0.390 pm 0.024$, a value larger than the commonly used theoretical value of 0.341 (Gavrilov & Zharkov, 1977), but compatible with more recent models (Helled & Guillot, 2013) for which $k_2$ ranges from 0.355 to 0.382. Depending on the assumed spin for Saturns interior, the new constraint can lead to a reduction of up to 80% in the number of potential models, offering great opportunities to probe the planets interior. In addition, significant tidal dissipation within Saturn is confirmed (Lainey et al., 2012) corresponding to a high present-day tidal ratio $k_2/Q=(1.59 pm 0.74) times 10^{-4}$ and implying fast orbital expansions of the moons. This high dissipation, with no obvious variations for tidal frequencies corresponding to those of Enceladus and Dione, may be explained by viscous friction in a solid core, implying a core viscosity typically ranging between $10^{14}$ and $10^{16}$ Pa.s (Remus et al., 2012). However, a dissipation increase by one order of magnitude at Rheas frequency could suggest the existence of an additional, frequency-dependent, dissipation process, possibly from turbulent friction acting on tidal waves in the fluid envelope of Saturn (Ogilvie & Li, 2004). Alternatively, a few of Saturns moons might themselves experience large tidal dissipation.
We investigate the resonant rotation of co-orbital bodies in eccentric and planar orbits. We develop a simple analytical model to study the impact of the eccentricity and orbital perturbations on the spin dynamics. This model is relevant in the entire domain of horseshoe and tadpole orbit, for moderate eccentricities. We show that there are three different families of spin-orbit resonances, one depending on the eccentricity, one depending on the orbital libration frequency, and another depending on the pericenters dynamics. We can estimate the width and the location of the different resonant islands in the phase space, predicting which are the more likely to capture the spin of the rotating body. In some regions of the phase space the resonant islands may overlap, giving rise to chaotic rotation.
We investigate the properties of the hydrodynamic flow around eccentric protoplanets and compare them with the often assumed case of a circular orbit. To this end, we perform a set of 3D hydrodynamic simulations of protoplanets with small eccentricities ($eleq 0.1$). We adopt an isothermal equation of state and concentrate resolution on the protoplanet to investigate flows down to the scale of the protoplanets circumplanetary disk (CPD). We find enhanced prograde rotation exterior to the CPD for low planet masses undergoing subsonic eccentric motion. If the eccentricity is made large enough to develop a bow shock, this trend reverses and rotation becomes increasingly retrograde. The instantaneous eccentric flow field is dramatically altered compared to circular orbits. Whereas the latter exhibit a generic pattern of polar inflow and midplane outflow, the flow geometry depends on orbital phase in the eccentric case. For even the modest eccentricities tested here, the dominant source of inflow can come from the midplane instead of the poles. We find that the amount of inflow and outflow increases for higher $e$ and lower protoplanet masses, thereby recycling more gas through the planets Bondi radius. These increased fluxes may increase the pebble accretion rate for eccentric planets up to several times that of the circular orbit rate. In response to eccentric motion, the structure and rotation of the planets bound CPD remains unchanged. Because the CPD regulates the eventual accretion of gas onto the planet, we predict little change to the gas accretion rates between eccentric and circular planets.
In a recent paper we proposed that the giant planets primordial orbits may have been eccentric (~0.05), and used a suite of dynamical simulations to show outcomes of the giant planet instability that are consistent with their present-day orbits. In this follow-up investigation, we present more comprehensive simulations incorporating superior particle resolution, longer integration times, and eliminating our prior means of artificially forcing instabilities to occur at specified times by shifting a planets position in its orbit. While we find that the residual phase of planetary migration only minimally alters the the planets ultimate eccentricities, our work uncovers several intriguing outcomes in realizations where Jupiter and Saturn are born with extremely large eccentricities (~0.10 and ~0.25, respectively). In successful simulations, the planets orbits damp through interactions with the planetesimal disk prior to the instability, thus loosely replicating the initial conditions considered in our previous work. Our results therefore suggest an even wider range of plausible evolutionary pathways are capable of replicating Jupiter and Saturns modern orbital architecture.