Draft version October 6, 2017 Preprint typeset using LATEX style emulateapj v. 12/16/11
DYNAMICAL EVOLUTION INDUCED BY PLANET NINE Konstantin Batygin1, Alessandro Morbidelli2 1Division of Geological and Planetary Sciences, California Institute of Technology, 1200 E. California Blvd., Pasadena, CA 91125 and 2Laboratoire Lagrange, Universit Cˆoted’Azur, Observatoire de la Cˆoted’Azur, CNRS, CS 34229, 06304 Nice, France Draft version October 6, 2017
ABSTRACT The observational census of trans-Neptunian objects with semi-major axes greater than 250 AU exhibits unexpected orbital structure that is most readily attributed to gravitational perturbations∼ induced by a yet-undetected, massive planet. Although the capacity of this planet to (i) reproduce the observed clustering of distant orbits in physical space, (ii) facilitate dynamical detachment of their perihelia from Neptune, and (iii) excite a population of long-period centaurs to extreme inclinations is well established through numerical experiments, a coherent theoretical description of the dynamical mechanisms responsible for these effects remains elusive. In this work, we characterize the dynamical processes at play, from semi-analytic grounds. We begin by considering a purely secular model of orbital evolution induced by Planet Nine, and show that it is at odds with the ensuing stability of distant objects. Instead, the long-term survival of the clustered population of long-period KBOs is enabled by a web of mean-motion resonances driven by Planet Nine. Then, by taking a compact-form approach to perturbation theory, we show that it is the secular dynamics embedded within these resonances that regulates the orbital confinement and perihelion detachment of distant Kuiper belt objects. Finally, we demonstrate that the onset of large-amplitude oscillations of orbital inclinations is accomplished through capture of low-inclination objects into a high-order secular resonance and identify the specific harmonic that drives the evolution. In light of the developed qualitative under- standing of the governing dynamics, we offer an updated interpretation of the current observational dataset within the broader theoretical framework of the Planet Nine hypothesis.
1. INTRODUCTION a & 250 AU Kuiper belt objects (Figure1) is sculpted by Over the course of one hundred and seventy years that a m 10 m planet residing on an appreciably eccentric (e ∼0.3 0.⊕7), large semi-major axis (a 300 700 AU) followed the successful prediction of Neptune (LeVerrier ∼ − ∼ − 1846; Adams 1846), the presence of additional planets orbit, whose plane roughly coincides with the plane of within the solar system has been contemplated by an the distant bodies, and is characterized by perihelion extensive list of astronomers and celestial mechanicians direction that is anti-aligned with respect to the aver- (Hoyt 1980). Historically, the lines of evidence for the age apsidal orientation of the KBOs. In addition to existence of distant, massive bodies that orbit the sun (i) facilitating the orbital confinement of the aforemen- have ranged from the (apparently) anomalous motion of tioned population of long-period KBOs and (ii) providing Uranus (Pickering & Pickering 1909; Lowell 1915) and a physical mechanism for the perihelion detachment of unexpected orbital characteristics of long-period comets Sedna-type orbits from Neptune, the presence of Planet (Forbes 1880; Matese & Whitmire 1986), to the peculiar Nine entails a series of additional consequences for the structure of the solar system’s small body populations observed structure of the solar system (Brown & Baty- (Brunini & Melita 2002; Gladman & Chan 2006; Gomes gin 2016). In particular, it has been shown that the dy- et al. 2006; Lykawka & Mukai 2008; Trujillo & Sheppard namical effects of Planet Nine naturally explain (iii) the 2014; Volk & Malhotra 2017). The predicted physical existence of highly inclined, large semi-major axis Cen- and orbital properties of the putative planets have been taurs (Gomes et al. 2015; Batygin & Brown 2016a), (iv) equally as varied, with inferred masses and semi-major the six-degree obliquity of the sun (Bailey et al. 2016; arXiv:1710.01804v1 [astro-ph.EP] 4 Oct 2017 axes spanning the Mars-Jupiter range and tens to thou- Lai 2016; Gomes et al. 2017), as well as (v) the origins sands of astronomical units, respectively. of proximate (a < 100 AU) retrograde KBOs (Batygin & A recent addition to the aggregate of planetary pre- Brown 2016b). dictions within the solar system is the Planet Nine hy- A key characteristic that differentiates the various pothesis1 (Batygin & Brown 2016a). Within the frame- planetary proposals is the dynamical mechanism through work of this model, the observed orbital clustering of which the envisaged planet generates its observational signatures. In this regard, the Planet Nine (P9) hypoth- esis remains incomplete. While numerical simulations 1 The Planet Nine hypothesis was inspired by the work of Tru- reveal that synthetic models of the solar system that in- jillo & Sheppard(2014), who noted that the arguments of perihe- lion (the angle between the apsidal and nodal lines on an orbit) of clude Planet Nine can provide a good match to the ob- distant Kuiper belt objects are grouped together. In contrast with servational data (see however Nesvorny et al. 2017), the this finding, the primary aim of Planet Nine’s inferred influence is dynamical process through which physical confinement to explain the simultaneous clustering of the longitudes of perihe- of the distant orbits occurs remains poorly understood. lion (a proxy for the direction of the pericenter in physical space) and the longitudes of ascending node (orientation of the orbital The original study of Batygin & Brown(2016a) sug- plane). gested that mean-motion resonances (including those 2 of high order) are responsible for orbital clustering. Expanding on this idea, Malhotra et al.(2016) have semi-major axis > 250 AU suggested that each of the distant KBOs are cur- inclination < 40 deg rently trapped in N:1 and N:2 mean motion resonances perihelion > 30 AU 2015 KG163 (MMRs) with Planet Nine. More recently, Millholland & Laughlin(2017); Becker et al.(2017) have carried out a large-scale numerical exploration of the resonant hy- scaled Runge-Lenz pothesis in relation to the existing data, further demon- 0 AU vectors strating its viability. Meanwhile, Beust(2016) has shown 25 that a purely secular treatment of the dynamics pro- 2015 GT50 vides a good match to the simulation results of Baty- 2012 VP113 gin & Brown(2016a), suggesting that resonant dynam- ics may be irrelevant to the problem at hand. In light of Sedna these conflicting results, the dominant P9-KBO interac- tion mechanism remains elusive, and its identification is the primary purpose of this paper. Analytical characterization of P9-KBO coupling is im- 2014 FE72 portant for three reasons, the first being falsifiability. 2013 SY99 The specific perturbation mechanism (along with the data it aims to explain) draws the distinction between Fig. 1.— Current observational census of the distant Kuiper belt. various planetary hypotheses, and can be used to refute Thirteen known objects with semi-major axes greater than 250 AU a specific model upon confrontation with the data. As an and perihelion distance greater than 30 AU are shown in physical example, suppose that an observational survey discovers space, and are color-coded in accord with their dynamical class, as dictated by the Planet Nine hypothesis. Orbits belonging to the a planet at a radial separation of a few hundred AU, but primary longitude of perihelion cluster (inferred to be apsidally the gravitational effects of this planet do not facilitate anti-aligned with the orbit of Planet Nine) are shown in purple. a physical confinement among the distant KBO orbits Orbits that are diametrically opposed to the primary cluster (pre- through the envisioned dynamical process2. Such a dis- sumed to be apsidally aligned with the orbit of Planet Nine) are shown in green. The outlying object that does not correspond to covery would imply that the Planet Nine hypothesis, as either population is shown in gray. Each orbit is further labeled formulated, is incorrect. by its scaled Runge-Lenz vector, which is color-coded in the same A second, more practical motivation for studying P9- way. KBO interactions is the link between existing observa- tions and the predicted orbit of Planet Nine. That is, if markable dynamical states of highly eccentric, resonant all distant KBOs are presently trapped in MMRs with exoplanets (Lee 2004; Tan et al. 2013). P9, their mean longitudes may contain information about Before delving into calculations, we delineate a list of the location of Planet Nine on its orbit (Malhotra et al. specific questions we wish to answer in this work. 2016; Millholland & Laughlin 2017). If, on the other hand, the interactions are purely secular, such a connec- What role (if any) do resonant interactions play tion cannot be established3. Therefore, understanding • within the dynamical evolution induced by Planet the physics of P9-KBO coupling can offer a more com- Nine? If resonances are prevalent, what or- plete interpretation of the extant observational dataset, der/multiplet harmonics dominate the dynamics, and may yield an avenue towards further constraining and what are their characteristic widths? the astronomical search for Planet Nine. The third, and final reason for this study is purely What role (if any) do secular interactions play academic. While the exploration of the circular re- • within the dynamical evolution induced by Planet stricted three-body problem (as a paradigm for inter- Nine? If dominant, how are close encounters actions between planets and small bodies in the solar avoided on nearly co-planar, anti-aligned orbits? system) dates back multiple centuries (Laplace 1799; Moreover, if resonant interactions are relevant to Poincar´e 1902), its highly elliptical counterpart remains the Planet Nine hypothesis, why does the purely much more scarcely understood (see e.g. Beaug´eet al. secular phase-space portrait provide a good match 2006; Michtchenko et al. 2006; Pichierri et al. 2017). In- to the results of numerical simulations? deed, a well-formulated theory for gravitational coupling between a planet and a test particle in the severe orbit- What parameters determine the critical semi- crossing regime does not currently exist. As a result, a • major axis corresponding to the transition perturbative study of P9 dynamics can yield interesting between randomized and clustered longitudes of insights into the general mathematical structure of inter- perihelion? What physical effect controls this actions among highly excited orbits. Such a framework transition? would have considerable applications beyond the prob- lem at hand, including the characterization of the re- What is the qualitative behavior of inclination dy- • namics within the framework of P9-driven evolu- 2 We note that the existence of a trans-Neptunian planet at an orbital separation of a few hundred AU was first proposed by tion? What dynamical process allows some of the Forbes(1880). objects to acquire exceptionally high inclinations 3 The process of orbital averaging inherent to secular perturba- in the distant Kuiper belt? tion theory removes all information related to the mean anomalies of the interacting bodies (Morbidelli 2002). We take these questions as an approximate guide to 3
anti-aligned orbits a = 50 AU a = 150 AU a = 250 AU 1 q = a8 0.75 Neptune crossing orbits e 0.5 q = a8 collisional domain 0.25
0 a = 350 AU a = 450 AU a = 550 AU 1
0.75 e 0.5 collisional collisional collisional domain domain domain 0.25
0 0 90 180 270 360 0 90 180 270 360 0 90 180 270 360 $ (deg) $ (deg) $ (deg) aligned orbits
Fig. 2.— Purely secular dynamics induced upon Kuiper belt objects by Planet Nine. Each panel is labeled by the corresponding value of the particle’s semi-major axis, and depicts contours of the doubly averaged Hamiltonian (1) in eccentricity - longitude of perihelion (measured relative to the apsidal line of Planet Nine) space. In panels characterized by a > 250 AU, secular equilibria corresponding to both apsidally aligned (∆$ = 0) as well as anti-aligned (∆$ = 180 deg) configurations emerge, and are segregated by tangential collision curves, shown as thick curves. Note that in spite of the apparent smoothness of the level curves of the averaged Hamiltonian surrounding the anti-aligned equilibrium points, this region of phase-space describes orbital configurations that are not protected from close encounters, and thus entails a long-term unstable orbital evolution. Moreover, notice that only apsidally aligned configurations that are close to the ∆$ = 0 equilibrium points are protected from close encounters by the geometrical collinearity of the orbits. the logic of the paper, which is structured as follows. ther, we choose to work in a slowly rotating coordinate In section2, we revisit a purely secular description of frame, that co-precesses with Planet Nine’s perihelion P9-induced orbital evoluition, carrying out the averag- (the corresponding contact transformation is spelled out ing procedure in closed form. In section3, we consider in Batygin & Brown 2016a). The governing (doubly av- a pair of idealized N-body simulations and outline the eraged) Hamiltonian of a test particle under planetary key differences between numerical experiments and pure perturbations then has the form: secular theory. In section4, we present a semi-analytical 8 2 theory of resonant P9-KBO interactions, and elucidate 1 M 1 X mj aj
s = G / secular dynamics embedded within mean-motion reso- H −4 a 23 2 M a2 1 e j nances as the primary driver of apsidal confinement of − =5 p p 2 distant KBOs. Subsequently, we discuss the onset of +$ ˙ 9 M a 1 1 e large-scale inclination oscillations of long-period bodies IIG − − in section5. We then re-examine the existing observa- 1 m9 2 G dλ dλ9, (1) tional data in light of the aforementioned theoretical de- − 4π r r9 | − | velopments in section6. We summarize and discuss the where is the gravitational constant, M is the mass implications of our results in section7. A brief analy- of theG sun, r is the position vector, λ is the mean lon- sis of long-term stable, but apsidally unconfined orbits is gitude, a is semi-major axis, e is eccentricity, and $ is presented in the appendix. the longitude of perihelion. All quantities pertaining to the four canonical giant planets are labeled with indexes 2. PURELY SECULAR DYNAMICS 5 8, while the values corresponding to Planet Nine are We begin our study of Planet Nine-induced dynam- denoted− with the subscript 9. The unlabeled variables ics within the framework of purely secular perturbation correspond to the Kuiper belt object. Finally, Planet theory. In Batygin & Brown(2016a), our preliminary Nine’s orbit-averaged precession rate is given by the ex- exploration of secular dynamics relied on an octupole or- pression4: der expansion of the gravitational potential (Kaula 1962; s 8 2 Mardling 2013), that implicitly assumed that the or- 3 M 1 X mj aj $˙ 9 = G . (2) bits under consideration do not cross. Following Beust 3 22 2 4 a9 1 e M a9 (2016), here we abandon the series-expansion approach − 9 j=5 to modeling P9-KBO coupling, and carry out the phase- The three terms present in equation (1) have sim- averaging procedure in closed form. ple physical interpretations. The first term governs the As in Batygin & Brown(2016a); Beust(2016), we as- sume coplanarity, and model the mean-field effects of the 4 Note that there is a typo in the corresponding expression in known giant planets neglecting their eccentricities. Fur- Batygin & Brown(2016a). 4 secular advance of the KBO’s perihelion due to the the contours of s do not show any kinks within the apsi- phase-averaged gravitational potential of Jupiter, Sat- dally anti-alignedH domain is simply a consequence of the urn, Uranus, and Neptune (in direct parallel with equa- integrability of the non-tangential singularity (Thomas tion2). The second term accounts for the fact that the & Morbidelli 1996; Gronchi & Milani 1998), and does reference frame is slowly rotating. The third term gov- not mean that the system can elude collisions. Instead, erns secular P9-KBO interactions. Note that the indirect recalling that the physical setup of the orbits is planar, part of the disturbing potential is by default entirely av- it is trivial to demonstrate that under the assumption eraged out within the framework of secular theory, and of uncorrelated Keplerian motion, all objects entrained need not be accounted for (Morbidelli 2002). in the apsidally anti-aligned configuration with Planet The only pair of dynamical variables present in Hamil- Nine would suffer close-encounters on timescales much tonian (1) is (e, ∆$), meaning that the system is inte- shorter than the age of the solar system. Therefore, grable. In other words, contours of the numerically aver- despite giving the illusion of agreement with N-body aged function (1) fully encapsulate the accompanying or- simulations, pure secular theory predicts that the entire bital evolution. Figure (2) shows the secular phase-space distant Kuiper belt is dynamically unstable, and should portraits of the system, projected into e ∆$ space, have been cleared out on a timescale comparable to the for a = 50, 150, ..., 550 AU, adopting P9 parameters− from orbital precession time. Batygin & Brown(2016a) (specifically, a9 = 700 AU, Contrary to this view, published numerical experi- e9 = 0.6 and m9 = 10 m ). This figure can be readily ments reveal that particles residing deep within the cores compared with Figure (4)⊕ of Batygin & Brown(2016a), of the anti-aligned libration regions of Figure (2) remain and confirms that the purely secular portrait provides a stable over the multi-Gyr lifetime of the solar system good match to the numerically computed portraits in the (Batygin & Brown 2016a; Brown & Batygin 2016). This same dynamical regime (Beust 2016). disparity suggests that secular theory alone is unlikely On all panels denoted by a > 250 AU, e ∆$ space to represent a full dynamical description of P9-driven is characterized by two stable equilibrium points:− one at evolution, and an additional stabilizing mechanism is at ∆$ = 0 and another at ∆$ = 180 deg. On each dia- play in the simulations. Let us now examine this point gram, the two libration regions surrounding these fixed further. points are separated by a solid curve, which corresponds to a tangential configuration of the KBO and P9 orbits. 3. NUMERICAL SIMULATIONS Note further that on these tangential contact curves, the To quantify the discrepancy between published simula- tions and pure secular theory, we carry out a pair of sim- derivatives of the Hamiltonian are discontinuous (i.e. s is locally class 0), signaling a breakdown of the secularH plified numerical experiments that mirror those reported framework (GronchiC 2002). in (Batygin & Brown 2016a; Brown & Batygin 2016). In particular, we evolve an initially axisymmetric disk A simple examination of the contours of depicted Hs of 6,000 eccentric test-particles with a (150, 750) AU on Figure (2) reveals that KBOs initialized on nearly and q (30, 36) AU, under the influence∈ of the phase- Neptune-crossing orbits (immediately below horizontal averaged∈ potential of the four inner giants as well as lines labeled by q = a8) will suffer drastically different Planet Nine with (a9, e9) = (700 AU, 0.6) and (a9, e9) = secular evolutions depending on their starting value of (600 AU, 0.5). The test particles are initialized with a ∆$. Bodies initially close to ∆$ 0 will be driven ∼ null vertical velocity dispersion and thus remain confined onto the tangential orbit-crossing curves through apsi- to Planet Nine’s orbital plane throughout the simulation. dal precession, and will eventually be removed from the Unlike the simulation suite of Batygin & Brown system by recurrent close encounters with Planet Nine. (2016a); Brown & Batygin(2016); Millholland & Laugh- On the other hand, objects initialized at ∆$ π settle ∼ lin(2017); Becker et al.(2017) where Neptune was mod- onto secular trajectories that encircle that anti-aligned eled directly, or that of Batygin & Brown(2016b) where equilibrium, and never encounter the tangential collision the Keplerian motion of all four inner giants was resolved, curves. Thus, the dynamical lifetimes of apsidally anti- here we emulate the effects of Jupiter, Saturn, Uranus aligned objects can be na¨ıvely envisioned to be longer and Neptune with an effective quadrupolar gravitational than their aligned counterparts, and in time, a cluster of moment of the Sun: exclusively anti-aligned orbits should be carved out by 2 8 2 Planet Nine. 3 cos (i9) 1 X mj aj J = − , (3) We note further that there is an island of stable ap- 2 4 M 2 sidal libration around ∆$ = 0 that avoids crossing the j=5 R tangential configuration curve and is thus protected from setting the inner absorbing radius to = 20 AU. This close-encounters. However, the eccentricities along these idealization is employed specifically to avoidR contaminat- protected orbits are moderate and the corresponding per- ing P9-induced KBO evolution with chaotic dynamics ihelion distances are large. This means that even if KBOs that arise from scattering off of Neptune, and to yield somehow came to occupy these islands of stability, they the closest point of comparison between numerical and would be difficult to detect observationally. semi-analytical results. Moreover, here we will assume With this picture in mind, it is tempting to affirm that that the inclination of Planet Nine relative to the Laplace the agreement between theory and simulation is satis- plane of the canonical giant planets is sufficiently small factory, and proceed forward within the purely secular to approximate cos(i9) 1, while keeping in mind that framework. This is however a misconception, facilitated zˆ axis of our coordinate≈ system coincides with the or- by the apparent smoothness of the secular phase-space bital− plane of Planet Nine. portraits shown in Figure (2). In particular, the fact that To carry out the simulations, we utilized the mercury6 5
15
m9 = 10 m 7:1 28 10 a9 = 700 AU 8:1 e =0.6 N 9
5 apsidal circulation
0 6:1 5:1 4:1 7:2 3:1 5:2 2:1 7:4 5:3 3:2 4:3 5:4 6:5 1:1 15
m9 = 10 m 25 10 9:2 a9 = 600 AU e =0.5 N 11:2 9 5 apsidal libration anti-aligned 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
orbital period ratio ( P / P9 ) Fig. 3.— A histogram summarizing the orbital distribution of simulated particles with dynamical lifetimes that exceed 4 Gyr. Beyond a critical orbital period ratio (P/P9 & 0.1 for a9 = 700 AU, e9 = 0.6 and P/P9 & 0.15 for a9 = 600 AU, e9 = 0.5), all surviving members of the distant scattered disk reside in mean-motion resonances with Planet Nine, and derive their long-term stability from the associated phase-protection mechanism. The final orbital configurations of distant bodies sculpted by Planet Nine are similar in both numerical experiments, and show the onset of clustering in longitude of perihelion beyond the 3:1 mean motion resonance. Correspondingly, bins containing particles locked in a stable pattern of anti-aligned apsidal libration are shown in black, while bins containing objects with circulating longitudes of perihelion are shown in gray. The commensurabilities possessing the largest number of apsidally confined particles are indicated with blue triangles. gravitational dynamics software package (Chambers it can be clearly seen in Figure (3) that all long-lived 1999). The integrations were performed using the hybrid objects which exhibit perihelion confinement have semi- Wisdom-Holman/Bulirsch-Stoer algorithm (Wisdom & major axes that correspond to mean-motion commensu- Holman 1992; Press et al. 1992), with a time-step of rabilities with Planet Nine. Particularly, the 1:1, 3:2, 2:1, ∆t = 3, 100 days (i.e. 1/10 th of Uranus’ orbital period), and 5:2 resonances contain the largest populations of ap- and spanned 4 Gyr. Any particle that attained a radial sidally clustered KBOs in both simulations5. This point distance of r < or r > 10, 000 AU was removed from undercuts a key difference between the results of N-body the simulations. R simulations and pure secular theory, and demonstrates Within the context of these idealized numerical exper- that (at least within the framework of a planar physical iments, all long-term stable particles retain nearly con- setup) distant KBOs derive their long-term orbital sta- stant semi-major axes throughout the integration, and bility from the phase-protection mechanism inherent to Figure (3) depicts an orbital histogram of the surviv- mean-motion resonances. ing bodies. Bins shown in black correspond to apsidally The typical dynamical evolutions of bodies trapped in confined objects (defined by ∆$ 180 deg 6 90 deg the four aforementioned resonances are shown in Fig- throughout the simulation), while| − those shown| in gray ure (4). Specifically, the left and middle panels depict denote particles that are long-term stable, but experience the evolutions of orbital eccentricities and apsidal lines. apsidal circulation. We note that here, the gray bins are On timescales of order 0.1 1 Gyr, orbital eccentric- stacked on top of the black bins, such that the relative ities experience considerable∼ oscillations− in concert with size of gray and black components of each column is a librations of ∆$. This facilitates a periodic regression of measure of the contamination of apsidally confined pop- the KBO perihelion distance, generating dynamically de- ulation of objects (at a given semi-major axis) by those tached (Sedna-type) orbits (Gladman et al. 2002; Brown undergoing perihelion circulation. et al. 2004). A strong propensity towards apsidal clustering for Of particular importance is the behavior of the reso- P/P9 & 1/3 is clearly evident in both simulations. How- nant angles: ever, it is also important to note that the same orbital ψ = k λ ` λ (k `) $, periods (e.g. those corresponding to the 3:1, 2:1, and 1:1 res 9 − − − commensurabilities in Figure3) can be simultaneously φres = k λ9 ` λ (k `) $9, (4) occupied by apsidally circulating and librating orbits. − − − where k and ` are integers (note that for k ` > 1, addi- This means that even within the framework of highly | − | idealized treatment of P9-induced dynamics, the cluster- tional resonant harmonics related to these angles through ing of the longitude of perihelion beyond a critical semi- ∆$ exist). On the right panels of Figure (4), ψres is plot- major axis is not perfectly strict, and the existence of ted in green and φres is plotted in black. The actual li- objects that do not conform to the general anti-aligned brations of φres and ψres are short-periodic (as portrayed orbital pattern is an expected consequence of the model 5 We note that these simulations are not intended to represent (Shankman et al. 2017). a full exploration of the resonant capture probabilities within the In addition to the already established tendency to- context of the Planet Nine hypothesis. A more complete estimation wards orbital clustering with increasing orbital period, these probabilities is presented in a companion paper (Bailey et al. 2017). 6
eccentricity e apsidal angle $ (deg) resonant angles (deg) 1 360 180 res res 0.75 270 90
0.5 180 0
0.25 90 -90 1:1 MMR a = 700 AU resonant libration around 0 deg 0 0 -180 1 360 360
0.75 270 270
0.5 180 180
0.25 90 90 3:2 MMR a = 534 AU resonant libration around 180 deg 0 0 0 1 360 180
0.75 270 90
0.5 180 0
0.25 90 -90 2:1 MMR a = 441 AU resonant libration around 0 deg 0 0 -180 1 360 360
0.75 270 270
0.5 180 180
0.25 90 90 5:2 MMR a = 380 AU resonant libration around 180 deg 0 0 0 0 1 2 3 4 0 1 2 3 4 0 1 2 3 4 time (Gyr) time (Gyr) time (Gyr)
Fig. 4.— Orbital time series of four resonant objects drawn from the N-body simulation with m9 = 10 m⊕, a9 = 700 AU and e9 = 0.6. The left and middle panels show oscillations of the eccentricities and the relative longitudes of perihelion respectively, while the right panels depict the behavior of the resonant angles φres (black) and ψres (green). Generally, particles that exhibit stable anti-aligned apsidal libration (such as those shown in this Figure) are characterized by low-amplitude oscillation of φres with a period equal to that of libration of ∆$ (not to be confused with the actual libration amplitude of φres, which is short periodic and corresponds to the thickness of the black “band”, as shown in the inset of the top-right panel). On the contrary, orbital evolution driven primarily by low-amplitude libration of a resonant harmonic that contains $ (rather than only $9 as in φres) corresponds to circulation of the longitude of perihelion (examples shown in the appendix). Note that φres librates around 0 for the 1:1 & 2:1 MMRs and around 180 deg for the 3:2 & 5:2 MMRs. by the inset within the top-right panel of Figure4) and oscillations of the these angles are inversely correlated to their amplitudes correspond to the thickness of the green the relative strengths of the corresponding terms. Ac- and black “bands.” Meanwhile, the long-periodic oscilla- cordingly, in our analytic approach in the next section tions of these angles are mere reflections of the librations we will consider φres as a reference angle, in contrast to of ∆$, which modulate the locations of the resonant Malhotra et al.(2016) and Beust(2016), who considered equilibria associated with φres and ψres. ψres instead. The fact that the amplitude of long-periodic oscilla- Only containing the longitude of perihelion of Planet tions of φres is much smaller than that of ψres for orbits Nine and not that of the Kuiper belt objects themselves, that exhibit libration in ∆$, demonstrates that the reso- the (short-periodic) libration of φres drives the oscilla- nant multiplets containing φres represent a better approx- tion of the particle’s semi-major axis but does not af- imation of the real dynamics than those containing ψres. fect the evolution of its eccentricity. Thus, the libra- In particular, if the Hamiltonian depended only on φres, tion in φres, sometimes referred to as “corotation” res- this angle would have no long-periodic oscillations while onance, already implies a certain disconnect among the ψres (being equal to φres + (` k) ∆$) would oscillate degrees of freedom related to the particle’s semi-major with the amplitude and period− of (` k) ∆$. The op- axis and its eccentricity. Moreover, the striking sepa- | − | posite would be true if the Hamiltonian dependent only ration of timescales associated with resonant (φres a) on ψ . Thus the relative amplitudes of long-periodic dynamics and secular (e ∆$) dynamics motivates− the res − 7 construction of the semi-analytical model that will fol- 1:1 MMR low, based on the adiabatic approach. 710 The dominant dependence of the resonant dynamics on φres, as opposed to another harmonic that contains $ in its critical argument, is central to maintaining apsidally 705 anti-aligned libration of the orbits. In fact, long-term stable particles whose resonant dynamics is driven by any 700 angle other than φ exhibit circulation of the longitude (AU)
res a of perihelion, and correspond to objects denoted as gray bins in Figure (3). Although characterizing the dynamics 695 in this transitionary semi-major axis range is not the primary purpose of this study, we present a brief analysis 690 of this mode of orbital evolution in the appendix. 3:2 MMR 538 4. SEMI-ANALYTICAL THEORY As shown in the previous section, all surviving mem- bers of the synthetic scattered disk that exhibit persis- 536 tent apsidal anti-alignment with Planet Nine’s orbit are locked into mean-motion resonances with P9. The pur- 534 (AU)
pose of the following analysis is thus to explain this be- a havior from semi-analytic grounds. To achieve this goal, we consider the isolated resonant behavior first. As in 532 section (2), we will abandon traditional methods of ce- lestial mechanics based on series-expansions, and cast our 530 analysis of the governing Hamiltonian in closed form. 2:1 MMR 445 4.1. Mean Motion Resonances
In order to carry out the averaging process numerically 443 and elucidate the φres a resonant dynamics, we set up a rigorous Hamiltonian− approach to the problem. As 441 is typical for the planar restricted three body problem, (AU) the canonical Poincar´eaction-angle variables for the test a particle are (Murray & Dermott 1999): 439 p Λ = M a λ = + $
G M 437 p p 2 Γ = M a 1 1 e γ = $ (5) 5:2 MMR
G − − − 382 where is the mean anomaly. In terms of these vari- ables,M the Hamiltonian of the problem is (Morbidelli 2002) 381 2 1 M 1 r r9 = G m · . (6) 380 r 9 3 (AU) H −2 Λ − G r r9 − r9 | − | | | a
The vector r can be written as a function of the canon- 379 collision curves ical variables (5). Similarly, r9 can be expressed as a function of the parameters a9, e9, which are assumed to 378 be fixed in time, as well as the angles λ9 and $9, which 0 90 180 270 360 are assumed to advance with fixed frequencies, n9 and $˙ 9 respectively. Thus, Hamiltonian (6) depends on time res (deg) through these two planetary angles. In order to remove explicit time dependence from equa- Fig. 5.— Resonant phase-space portraits of the four most pop- ulated mean-motion commensurabilities, projected into a − φres tion (6), we extend the phase space by two degrees of space. The color scale and the black contours correspond to level freedom. That is, we consider λ and γ = $ as 9 9 − 9 curves of the singly averaged Hamiltonian (10), while the verti- independent variables, with conjugated actions Λ9 and cal black lines denote collision curves. Resonant contours of the Γ . Correspondingly, we add the term (n Λ $˙ Γ ) Hamiltonian that come into contact with the collision curves are 9 9 9 − 9 9 shown as thick purple lines, and thereby inform the widths of the to expression (6), making the Hamiltonian autonomous. corresponding resonances. These phase-space diagrams adopt our Of course, in doing so we do not alter the dynam- fiducial P9 parameters and assume that the Kuiper belt objects are ics in any way, since dλ9/dt = ∂ r/∂Λ9 = n9 and characterized by a fixed value of Γ that corresponds to q = 35 AU H at nominal resonance and ∆γ = π. The red curves that encircle d$9/dt = dγ9/dt = ∂ r/∂Γ9 =$ ˙ 9. − − H the elliptic equilibria in each panel represent the a − φres evolu- We now denote the mean motion of the particle by tion of the particles shown in Figure (3) near ∆$ = π, and signal 2 2 3 n = M /Λ and assume that there is a resonance excellent agreement with theory. of theG kind k n ` n = 0 for some integers k and `. 9 − 8
resonant domain occupied by stable apsidally anti-aligned orbits
3:2 MMR $ = $col $ = 100 deg $ = 140 deg $ = 180 deg 538
536
534 (AU) a
532
530 0 90 180 270 360 0 90 180 270 360 0 90 180 270 360 0 90 180 270 360