Arxiv:1111.3682V1 [Astro-Ph.EP] 15 Nov 2011 Ant Planets, As Well As Their Dynamical Architecture (I.E
Total Page:16
File Type:pdf, Size:1020Kb
Draft version November 17, 2011 Preprint typeset using LATEX style emulateapj v. 11/10/09 INSTABILITY-DRIVEN DYNAMICAL EVOLUTION MODEL OF A PRIMORDIALLY 5 PLANET OUTER SOLAR SYSTEM Konstantin Batygin1, Michael E. Brown1 & Hayden Betts2 1Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 91125 and 2Polytechnic School, Pasadena, CA 91106 Draft version November 17, 2011 ABSTRACT Over the last decade, evidence has mounted that the solar system's observed state can be favorably reproduced in the context of an instability-driven dynamical evolution model, such as the \Nice" model. To date, all successful realizations of instability models have concentrated on evolving the four giant planets onto their current orbits from a more compact configuration. Simultaneously, the possibility of forming and ejecting additional planets has been discussed, but never successfully implemented. Here we show that a large array of 5-planet (2 gas giants + 3 ice giants) multi-resonant initial states can lead to an adequate formation of the outer solar system, featuring an ejection of an ice giant during a phase of instability. Particularly, our simulations demonstrate that the eigenmodes which characterize the outer solar system's secular dynamics can be closely matched with a 5-planet model. Furthermore, provided that the ejection timescale of the extra planet is short, orbital excitation of a primordial cold classical Kuiper belt can also be avoided in this scenario. Thus the solar system is one of many possible outcomes of dynamical relaxation and can originate from a wide variety of initial states. This deems the construction of a unique model of solar system's early dynamical evolution impossible. Subject headings: Planets and satellites: dynamical evolution and stability; Kuiper belt: general 1. INTRODUCTION cal population of the Kuiper belt (Batygin et al. 2011). In the last two decades, discoveries of the Kuiper belt Finally, the presence of Jupiter's and Neptune's Trojan (Jewitt & Luu 1993), as well as planets orbiting stars asteroids has been attributed to chaotic capture of plan- other than the Sun (Mayor & Queloz 1995), have sup- etesimals during the instability (Morbidelli et al. 2005; plied the centuries-old quest to understand the formation Nesvorn´yet al. 2007). of the solar system with fresh constraints and insights All successful realizations of the Nice model to date into physical processes at play. Among a multitude of have been comprised exclusively of the four currently newly proposed formation scenarios, the \Nice" model present giant planets. However, there exists no strong (Tsiganis et al. 2005; Gomes et al. 2005; Morbidelli et al. evidence that suggests that additional planets were not 2005) is particularly notable, as it has attained a con- present in the solar system, at the epoch of the dispersion siderable amount of success in reproducing the various of the nebula. In fact, theoretical arguments, presented observed features of the solar system. Within the con- by Goldreich et al. (2004) point to a possibility of ini- text of the scenario envisioned by the Nice model, giant tially forming as many as five ice giants, three of which planets start their post-nebular evolution in a compact, get subsequently removed via ejections (see however Lev- multi-resonant configuration, and following a brief period ison & Morbidelli (2007)). The dynamical sensibility of of dynamical instability, scatter onto their current orbits such a scenario is further strengthened by the fact that a (Morbidelli et al. 2007; Batygin & Brown 2010; Walsh et considerable fraction of standard Nice model simulations al. 2011). result in an ejection of an ice-giant after an encounter The first success of the Nice model lies in its ability with at least one of the gas giants. to quantitatively reproduce the observed orbits of the gi- In this paper, we explore an instability-driven dynam- arXiv:1111.3682v1 [astro-ph.EP] 15 Nov 2011 ant planets, as well as their dynamical architecture (i.e. ical evolution of a 5-planet system (2 gas giants + 3 ice secular eigenmodes of the system) (Tsiganis et al. 2005; giants) with an eye towards identifying a pathway to- Morbidelli et al. 2009). Simultaneously, the brief insta- wards reproduction of solar system-like dynamical archi- bility, inherent to the model, provides a natural trigger tecture. In principle, the realm of possibility available to to the Late Heavy Bombardment (Gomes et al. 2005), as this study is enormous. Consequently, rather than per- well as a transport mechanism for emplacement of dy- forming a comprehensive parameter-search, here we limit namically \hot" Kuiper belt objects (KBOs) from inside ourselves to systems that contain an additional Uranus- of ∼ 35AU (Levison et al. 2008). Meanwhile, it has been like planet, with the aim of presenting a few proof-of- recently demonstrated that survival of a dynamically concept numerical experiments. The plan of the paper is \cold" primordial population between Neptune's current as follows. In section 2, we describe our numerical setup. 3:2 and 2:1 exterior mean-motion resonances (MMRs) is In section 3, we show that in our model, the planetary fully consistent with a Nice model-like evolution of the orbits, their secular eigenmodes, as well as various popu- planets, implying an in-situ formation of the cold classi- lations of the Kuiper belt are approximately reproduced. We conclude and discuss our results in section 4. [email protected] 2 Batygin Brown Betts 40 40 B 3:2 2:1 5:4 5:4 initial condition A 3:2 3:2 2:1 4:3 initial condition 30 30 (AU) (AU) 20 20 a, q, Q q, a, a, q, Q q, a, innermost ice giant ejection innermost ice giant ejection 10 10 0 0 0 1.107 2.107 3.107 4.107 5.107 0 1.107 2.107 3.107 4.107 5.107 t (years) t (years) Fig. 1.| Orbital evolution of planets. Each planet's semi-major axis, as well as perihelion and apohelion distances are shown as functions of time. The actual perihelion and apohelion distances of the planets are also shown for comparison as black error bars. In both cases, the innermost ice-giant ejects during the transient phase of instability, leaving behind 4 planets, whose orbits resemble that of the solar system. See main text for a description of the initial conditions. 2. NUMERICAL EXPERIMENTS described above). Each integration was performed using the mercury6 integration software package (Chambers The numerical setup of the simulations performed here 1 was qualitatively similar to those presented by Batygin 1999) and spanned 50Myr . The calculations were per- & Brown (2010) and Batygin et al. (2011). Particularly, formed on Caltech's PANGU super-computer. the five giant planets were initialized in a compact, multi- After their completion, simulations that were deemed resonant initial condition, surrounded by a massive plan- successful were reintegrated with the use of tracer simu- etesimal disk that extended between its immediate sta- lations (see Levison et al. (2008); Batygin et al. (2011)), bility boundary and 30AU. to address the dynamical evolution of a locally formed Two of the three ice-giants were taken to have the same population of KBOs. In particular, each run was sup- mass as Uranus and were initially placed next to Saturn plemented with an additional disk of mass-less particles and as the outermost planet respectively. The middle that resided in the cold classical region of the Kuiper ice-giant was taken to have Neptune's mass. In all sim- belt (i.e. between the final exterior 3:2 and 2:1 MMRs of ulations, Jupiter and Saturn started out in a 3:2 MMR, Neptune). in accord with the results of hydrodynamical simulations 3. RESULTS of convergent migration of the planets in the solar neb- ula (Masset & Snellgrove 2001; Morbidelli & Crida 2007; Out of the 810 integrations that were initially per- Pierens & Nelson 2008). The ice giants were also se- formed, 214 (∼ 25%) cases featured an ejection of a single quentially assembled into first order MMRs by apply- ice-giant, yielding a system composed of 4 planets. Per- ing dissipative forces, designed to mimic the presence of haps unsurprisingly, in most cases the ejected planet is the nebula (Lee & Peale 2002). Following resonant lock- the ice-giant that neighbors Saturn. Of the 750 simu- ing, each assembled multi-resonant initial condition was lations that were performed following the elimination of evolved in isolation for 10Myr, as an immediate test of initial conditions, 33 evolutions resulted in orbits remi- orbital stability. This procedure yielded a total of 81 niscent of the solar system. Specifically, we searched for stable multi-resonant initial conditions. solutions where the Saturn-Jupiter period ratio exceed The search for adequate dynamical evolutions was per- 2, while the final semi-major axes of the ice-giants were formed in two steps. First, we evolved 10 permutations of within 3AU of their observed counterparts. No strong each initial condition, with planetesimal disks composed requirements were placed on the planetary eccentricities of N = 1000 planetesimals. Disk masses were chosen and inclinations. min max It is noteworthy that evolving the giant planets onto so- randomly between Mdisk = 25M⊕ and Mdisk = 100M⊕. The density profiles followed a power-law distribution, lar system-like orbits is insufficient for a simulation to be Σ / rk where the power-law index, k, was chosen ran- deemed successful for indeed, there are additional con- domly between kmin = 1 and kmax = 2. The plan- straints that must be satisfied. The first orbital con- etesimals were initialized on near-circular orbits (e ∼ straint is the reproduction of the secular architecture of sin i ∼ 10−3). To reduce the already substantial com- the system.