skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Suppression of the inclination instability in the trans-Neptunian Solar system
ABSTRACT The trans-Neptunian scattered disc exhibits unexpected dynamical structure, ranging from an extended dispersion of perihelion distance to a clustered distribution in orbital angles. Self-gravitational modulation of the scattered disc has been suggested in the literature as an alternative mechanism to Planet nine for sculpting the orbital architecture of the trans-Neptunian region. The numerics of this hypothesis have hitherto been limited to N < O(103) superparticle simulations that omit direct gravitational perturbations from the giant planets and instead model them as an orbit-averaged (quadrupolar) potential, through an enhanced J2 moment of the central body. For sufficiently massive discs, such simulations reveal the onset of collective dynamical behaviour – termed the ‘inclination instability’ – wherein orbital circularisation occurs at the expense of coherent excitation of the inclination. Here, we report N = O(104) GPU-accelerated simulations of a self-gravitating scattered disc (across a range of disc masses spanning 5–40 M⊕) that self-consistently account for intraparticle interactions as well as Neptune’s perturbations. Our numerical experiments show that even under the most favourable conditions, the inclination instability never ensues. Instead, due to scattering, the disc depletes. While our calculations show that a transient lopsided structure can emerge within the first few hundreds of Myr, the terminal outcomes of these calculations systematically reveal a scattered disc that is free of any orbital clustering. We conclude thus that the inclination instability mechanism is an inadequate explanation of the observed architecture of the Solar system.  more » « less
Award ID(s):
2109276
PAR ID:
10528117
Author(s) / Creator(s):
;
Publisher / Repository:
OUP
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
523
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
6103 to 6113
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract There is a complex inclination structure present in the trans-Neptunian object (TNO) orbital distribution in the main classical-belt region (between orbital semimajor axes of 39 and 48 au). The long-term gravitational effects of the giant planets make TNO orbits precess, but nonresonant objects maintain a nearly constant “free” inclination (Ifree) with respect to a local forced precession pole. Because of the likely cosmogonic importance of the distribution of this quantity, we tabulate free inclinations for all main-belt TNOs, each individually computed using barycentric orbital elements with respect to each object’s local forcing pole. We show that the simplest method, based on the Laplace–Lagrange secular theory, is unable to give correct forcing poles for objects near theν18secular resonance, resulting in poorly conservedIfreevalues in much of the main belt. We thus instead implemented an averaged Hamiltonian to obtain the expected nodal precession for each TNO, yielding significantly more accurate free inclinations for nonresonant objects. For the vast majority (96%) of classical-belt TNOs, theseIfreevalues are conserved to < 1° over 4 Gyr numerical simulations, demonstrating the advantage of using this well-conserved quantity in studies of the TNO population and its primordial inclination profile; our computed distributions only reinforce the idea of a very coplanar surviving “cold” primordial population, overlain by a largeI-width implanted “hot” population. 
    more » « less
  2. Abstract A giant-impact origin for the Moon is generally accepted, but many aspects of lunar formation remain poorly understood and debated. Ćuk et al. proposed that an impact that left the Earth–Moon system with high obliquity and angular momentum could explain the Moon’s orbital inclination and isotopic similarity to Earth. In this scenario, instability during the Laplace Plane transition, when the Moon’s orbit transitions from the gravitational influence of Earth’s figure to that of the Sun, would both lower the system’s angular momentum to its present-day value and generate the Moon’s orbital inclination. Recently, Tian & Wisdom discovered new dynamical constraints on the Laplace Plane transition and concluded that the Earth–Moon system could not have evolved from an initial state with high obliquity. Here we demonstrate that the Earth–Moon system with an initially high obliquity can evolve into the present state, and we identify a spin–orbit secular resonance as a key dynamical mechanism in the later stages of the Laplace Plane transition. Some of the simulations by Tian & Wisdom did not encounter this late secular resonance, as their model suppressed obliquity tides and the resulting inclination damping. Our results demonstrate that a giant impact that left Earth with high angular momentum and high obliquity (θ> 61°) is a promising scenario for explaining many properties of the Earth–Moon system, including its angular momentum and obliquity, the geochemistry of Earth and the Moon, and the lunar inclination. 
    more » « less
  3. Abstract Axisymmetric disks of eccentric orbits in near-Keplerian potentials are unstable and undergo exponential growth in inclination. Recently, Zderic et al. showed that an idealized disk then saturates to a lopsided mode. Here we show, using N -body simulations, that this apsidal clustering also occurs in a primordial Scattered Disk in the outer solar system, which includes the orbit-averaged gravitational influence of the giant planets. We explain the dynamics using Lynden-Bell's mechanism for bar formation in galaxies. We also show surface density and line-of-sight velocity plots at different times during the instability, highlighting the formation of concentric circles and spiral arms in velocity space. 
    more » « less
  4. Abstract The most distant known trans-Neptunian objects (TNOs), those with perihelion distance above 38 au and semimajor axis above 150 au, are of interest for their potential to reveal past, external, or present but unseen perturbers. Realizing this potential requires understanding how the known planets influence their orbital dynamics. We use a recently developed Poincaré mapping approach for orbital phase space studies of the circular planar restricted three-body problem, which we have extended to the case of the 3D restricted problem withNplanetary perturbers. With this approach, we explore the dynamical landscape of the 23 most distant TNOs under the perturbations of the known giant planets. We find that, counter to common expectations, almost none of these TNOs are far removed from Neptune’s resonances. Nearly half (11) of these TNOs have orbits consistent with stable libration in Neptune’s resonances; in particular, the orbits of TNOs 148209 and 474640 overlap with Neptune’s 20:1 and 36:1 resonances, respectively. Five objects can be ruled currently nonresonant, despite their large orbital uncertainties, because our mapping approach determines the resonance boundaries in angular phase space in addition to semimajor axis. Only three objects are in orbital regions not appreciably affected by resonances: Sedna, 2012 VP113 and 2015 KG163. Our analysis also demonstrates that Neptune’s resonances impart a modest (few percent) nonuniformity in the longitude of perihelion distribution of the currently observable distant TNOs. While not large enough to explain the observed clustering, this small dynamical sculpting of the perihelion longitudes could become relevant for future, larger TNO data sets. 
    more » « less
  5. Abstract We investigate different conditions, including the orbital and size–frequency distribution (SFD) of the early Kuiper Belt, that can trigger catastrophic planetesimal destruction. The goal of this study is to test if there is evidence for collisional grinding in the Kuiper Belt that has occurred since its formation. This analysis has important implications for whether the present-day SFD of the cold classical trans-Neptunian objects (TNOs) is a result of collisional equilibrium or if it reflects the primordial stage of planetesimal accretion. As an input to our modeling, we use the most up-to-date debiased OSSOS++ ensemble sample of the TNO population and orbital model based on the present-day architecture of the Kuiper Belt. We calculate the specific impact energies between impactor–target pairs from different TNO groups and compare our computed energies to catastrophic disruption results from smoothed particle hydrodynamics simulations. We explore different scenarios by considering different total primordial Kuiper Belt masses and power slopes of the SFD and allowing collisions to take place over different timescales. The collisional evolution of the Kuiper Belt is a strong function of the unknown initial mass in the trans-Neptunian region, where collisional grinding of planetesimals requires a total primordial Kuiper Belt mass of M > 5 M ⊕ , collision speeds as high as 3 km s −1 , and collisions over at least 0.5 Gyr. We conclude that presently, most of the collisions in the trans-Neptunian region are in the cratering rather than disruption regime. Given the low collision rates among the cold classical Kuiper Belt objects, their SFD most likely represents the primordial planetesimal accretion. 
    more » « less