trans-Neptunian Region Research Articles

A Near-infrared Survey of Candidate Haumea Family Members

The Haumea family is the only known dynamical family in the trans-Neptunian region. To date, 10 family members have been unambiguously identified using near-infrared (NIR) spectral or photometric data in combination with their strong dynamical proximity and the rest of the family. In this work, we build off previous empirically constructed models of the family to identify 39 candidate family members and follow up on eight of them using the Hubble Space Telescope (HST) to measure their visible and NIR colors. Six of the candidates have strong water-ice absorption features—consistent with family membership. Based on these initial findings, our sample of 39 candidate family members should contain about 20 more water-rich objects. Combining the HST visible and NIR photometry with past results, we find no evidence for significant color heterogeneity within the family. Of the six new family members, two have Δv ∼ 300 m s−1, well outside of the traditionally defined velocity dispersion limit of ∼150 m s−1. As evidence suggests they are not affected by any of Neptune’s resonances, we propose that these family members are the result of dynamical sculpting by Neptune during its outward migration. Further searches for far-flung family members will be able to further explore this hypothesis.

Testing MOND on Small Bodies in the Remote Solar System

Modified Newtonian dynamics (MOND), which postulates a breakdown of Newton's laws of gravity/dynamics below some critical acceleration threshold, can explain many otherwise puzzling observational phenomena on galactic scales. MOND competes with the hypothesis of dark matter, which successfully explains the cosmic microwave background and large-scale structure. Here we provide the first solar system test of MOND that probes the subcritical acceleration regime. Using the Bekenstein–Milgrom “aquadratic Lagrangian” (or AQUAL) formulation, we simulate the evolution of myriads of test particles (planetesimals or comets) born in the trans-Neptunian region and scattered by the giant planets over the lifetime of the Sun to heliocentric distances of 102–105 au. We include the effects of the Galactic tidal field and passing stars. While Newtonian simulations reproduce the distribution of binding energies of long-period and Oort-cloud comets detectable from Earth, MOND-based simulations do not. This conclusion is robust to plausible changes in the migration history of the planets, the migration history of the Sun, the MOND transition function, effects of the Sun's birth cluster, and the fading properties of long-period comets. For the most popular version of AQUAL, characterized by a gradual transition between the Newtonian and MOND regimes, our MOND-based simulations also fail to reproduce the orbital distribution of trans-Neptunian objects in the detached disk (perihelion q > 38 au). Our results do not rule out some MOND theories more elaborate than AQUAL, in which non-Newtonian effects are screened on small spatial scales, at small masses, or in external gravitational fields comparable in strength to the critical acceleration.

Wind erosion and transport on planetesimals

We consider the possibility that aeolian (wind blown) processes occur on small, 1 to 100 km diameter, planetesimals when they were embedded in the protosolar nebula. Drag from a headwind within a protostellar disk is sufficiently large to loft cm and smaller sized particles off the surface of a 10 km diameter asteroid in the inner solar system (at a few AU), and micron sized particles off the surface of a 10 km diameter object in the Transneptunian region. The headwind is sufficiently strong to overcome surface cohesion in the inner solar system, but not in the outer solar system. However, in the outer solar system, surface particles can be redistributed or escape due to impacts from particles that are in the protosolar disk’s wind. Based on scaling crater ejecta, we estimate that impacts from particles in the headwind will lead to erosion of mass rather than accretion for planetesimals below about 6 km in diameter. The erosion limit is independent of material strength but proportional to the wind velocity. We explore the sensitivity of splash particle trajectories to particle size, headwind velocity and Reynolds number. Winds from a protostellar disk could account for Kuiper Belt Object (486958) Arrokoth’s smooth undulating terrain but only during an epoch of high particle flux and low wind velocity. These conditions could have been present during and just after coalescence of Arrokoth’s building blocks.

Spectroscopy of the binary TNO Mors–Somnus with the JWST and its relationship to the cold classical and plutino subpopulations observed in the DiSCo-TNO project

Context. Trans-Neptunian objects (TNOs) are remnants of small icy bodies from planetary formation that orbit in the region beyond Neptune. Within the population of TNOs, Trans-Neptunian binaries (TNBs) provide a valuable opportunity to test the models of the formation and evolution of planetesimals in the trans-Neptunian region. Various theories have been proposed to describe the observed separations between binary components, their relative sizes, and other orbital parameters. The colors of TNOs have been used to trace the dynamical history of the outer Solar System and the colors of TNB components provide tests for formation theories. However, spectral information for the components of small TNBs, crucial information that could validate formation mechanisms, has until now remained elusive. Aims. The main goal of this work is to characterize the near-infrared spectral properties of the TNB plutino (341520) Mors–Somnus, the only TNB with resolved components in the James Webb Space Telescope (JWST) Large Cycle 1 General Observer program “DiSCo-TNOs” (PID 2418; PI: Pinilla-Alonso). The secondary goal is to use the surface compositions of the individual components of the Mors–Somnus system to probe formation and dynamical evolution in the outer Solar System through comparison to the surface properties of the cold classical and plutino (3:2 resonant) dynamical groups. Methods. To achieve these goals, we measured the spectral slope of the continuum and identified absorption bands in the individual spectra of Mors and Somnus, as well as in those of the cold classicals and plutinos obtained with the NIRSpec Integral Field Unit (IFU) and the PRISM/CLEAR disperser (0.6–5.3 μm), and compared these results to shed light on the dynamical evolution of the Mors–Somnus binary. Results. The spectra of Mors and Somnus are similar and indicate the presence of complex organic materials, CO2, CO, OH-compounds, and tentative nitrogen-rich materials. We find a high degree of compositional diversity in the plutino population, a group of TNOs that likely formed elsewhere and moved to their current orbits during the migration of Neptune, while the cold classical TNOs, which likely formed in situ, appear more homogeneous. Conclusions. The very wide separation between the components, their nearly equal sizes, and the high orbital inclination of the system suggest this plutino binary is a survivor of the primordial population of objects beyond 30 au. The similarities found between the spectral features of the plutinos Mors and Somnus and those of all of the cold classical TNOs in the DiSCo-TNOs sample as well as the high degree of compositional heterogeneity found in the plutino population provide compositional evidence for evaluation of Neptune’s migration in the trans-Neptunian region early on in the history of the Solar System.

The origin of an overpopulation of contact binary plutinos

ABSTRACT In the trans-Neptunian region, 10–25 per cent of binary objects are contact binaries. However, the amplitude of the light curves of a sample of plutinos indicates that this fraction could reach 50 per cent or even more. This remarkable difference cannot be explained by the single action of the Kozai cycles and tidal friction (KCTF), nor it can be assumed that it is due to a formation mechanism different from the one that gave rise to the binaries of the other hot populations. In this work, we present numerical simulations showing that a superabundance of contact binaries occurs as a natural consequence of the coupling between the dynamical evolution of the mutual orbit of the binaries subject to the action of KCTF and the variation of the heliocentric orbit of the binary centre of mass due to the perturbations of the giant planets. This same mechanism affects the other populations, but to a lesser extent, because the variations of the heliocentric orbits are much smaller for them than for objects that are within the 3:2 mean motion resonance with Neptune. This mechanism also increases the fraction of tight binaries in the plutino population, where it could reach $\sim 30~{{\ \rm per\ cent}}$.

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.

Steeper Scattered Disks Buckle Faster

Disks of low-mass bodies scattered by giant planets to large semimajor axis and constant periapsis orbits are vulnerable to a buckling instability. This instability exponentially grows orbital inclinations, raises periapsis distances, and coherently tilts orbits resulting in clustering of arguments of periapsis. The dynamically hot system is then susceptible to the formation of a lopsided mode. Here we show that the timescale of the buckling instability decreases as the radial surface density of the population becomes more centrally dense, i.e., steeper scattered disks buckle faster. Accounting for differential apsidal precession driven by giant planets, we find that ∼10 M ⊕ is sufficient for a primordial scattered disk in the trans-Neptunian region to have been unstable if .

A Dynamical Survey of the Trans-Neptunian Region. I. Mean-motion Resonances with Neptune

In this paper, we present a large-scale dynamical survey of the trans-Neptunian region, with particular attention to mean-motion resonances (MMRs). We study a set of 4121 trans-Neptunian objects (TNOs), a sample far larger than in previous works. We perform direct long-term numerical integrations that enable us to examine the overall dynamics of the individual TNOs as well as to identify all MMRs. For the latter purpose, we apply an own-developed FAIR method that allows the semi-automatic identification of even very high-order MMRs. Apart from searching for the more frequent eccentricity-type resonances that previous studies concentrated on, we set our method to allow the identification of inclination-type MMRs, too. Furthermore, we distinguish between TNOs that are locked in a given MMR throughout the whole integration time span (108 yr) and those that are only temporarily captured in resonances. For a more detailed dynamical analysis of the trans-Neptunian space, we also construct dynamical maps using test particles. Observing the fine structure of the 34–80 au region underlines the stabilizing role of the MMRs, with the regular regions coinciding with the positions of real TNOs.

Comet P/2021 HS (PANSTARRS) and the Challenge of Detecting Low-activity Comets

Jupiter-family comet (JFC) P/2021 HS (PANSTARRS) only exhibits a coma within a few weeks of its perihelion passage at 0.8 au, which is atypical for a comet. Here we present an investigation into the underlying cause using serendipitous survey detections and targeted observations. We find that the detection of the activity is caused by an extremely faint coma being enhanced by the forward scattering effect owing to the comet reaching a phase angle of ∼140°. The coma morphology is consistent with sustained, sublimation-driven activity produced by a small active area, ∼700 m2, one of the smallest values ever measured on a comet. The phase function of the nucleus shows a phase coefficient of 0.035 ± 0.002 mag deg−1, implying an absolute magnitude of H = 18.31 ± 0.04 and a phase slope of G = − 0.13, with color consistent with typical JFC nuclei. Thermal observations suggest a nucleus diameter of 0.6–1.1 km, implying an optical albedo of 0.04–0.23, which is higher than typical cometary nuclei. An unsuccessful search for dust trail and meteor activity confirms minimal dust deposit along the orbit, totaling ≲108 kg. As P/2021 HS is dynamically unstable, similar to typical JFCs, we speculate that it has an origin in the trans-Neptunian region and that its extreme depletion of volatiles is caused by a large number of previous passages to the inner solar system. The dramatic discovery of the cometary nature of P/2021 HS highlights the challenges of detecting comets with extremely low activity levels. Observations at high phase angle, where forward scattering is pronounced, will help identify such comets.

Tidally locked rotation of the dwarf planet (136199) Eris discovered via long-term ground-based and space photometry

The rotational states of the members in the dwarf planet-satellite systems in the trans-Neptunian region are determined by formation conditions and the tidal interaction between the components. These rotational characteristics serve as prime tracers of their evolution. A number of authors have claimed a very broad range of values for the rotation period for the dwarf planet Eris, ranging from a few hours to a rotation that is (nearly) synchronous with the orbital period (15.8 d) of its satellite, Dysnomia. In this Letter, we present new light curve data for Eris, taken with ∼1–2 m-class ground based telescopes and with the TESS and Gaia space telescopes. The TESS data did not provide a well-defined light curve period, but it could be used to constrain light curve variations to a maximum possible light curve amplitude of Δm ≤ 0.03 mag (1-σ) for P ≤ 24 h periods. Both the combined ground-based data and Gaia measurements unambiguously point to a light curve period equal to the orbital period of Dysnomia, P = 15.8 d, with a light curve amplitude of Δm ≈ 0.03 mag, indicating that the rotation of Eris is tidally locked. Assuming that Dysnomia has a collisional origin, calculations with a simple tidal evolution model show that Dysnomia must be relatively massive (mass ratio of q = 0.01–0.03) and large (radius of Rs ≥ 300 km) to have the potential to slow Eris down to a synchronised rotation. These simulations also indicate that (assuming tidal parameters usually considered for trans-Neptunian objects) the density of Dysnomia should be 1.8–2.4 g cm−3. This is an exceptionally high value among similarly sized trans-Neptunian objects, setting important constraints on their formation conditions.

OSSOS. XXVI. On the Lack of Catastrophic Collisions in the Present Kuiper Belt

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.

The Von Zeipel–Lidov–Kozai Effect inside Mean Motion Resonances with Applications to Trans-Neptunian Objects

Secular dynamics inside mean motion resonances (MMRs) plays an essential role in governing the dynamical structure of the trans-Neptunian region and sculpting the orbital distribution of trans-Neptunian objects (TNOs). In this study, semianalytical developments are made to explore the von Zeipel–Lidov–Kozai resonance inside MMRs. To this end, a semi-secular model is formulated from averaging theory and then a single-degree-of-freedom integrable model is achieved based on the adiabatic invariance approximation. In particular, we introduce a modified adiabatic invariant, which is continuous around the separatrices of MMRs. During long-term evolution, both the resonant Hamiltonian and the adiabatic invariant remain unchanged, thus phase portraits can be produced by plotting level curves of the adiabatic invariant with a given Hamiltonian. The phase portraits provide global pictures to predict long-term behaviors of the eccentricity, inclination, and argument of pericenter. Applications to some representative TNOs inside MMRs (2018 VO137, 2005 SD278, 2015 PD312, Pluto, 2004 HA79, 1996 TR66, and 2014 SR373) show good agreements between the numerically propagated trajectories under the full N-body model and the level curves arising in phase portraits. Interestingly, 2018 VO137 and 2005 SD278 exhibit switching behaviors during their long-term evolution and currently they are inside 2:5 MMR with Neptune.

OSSOS Finds an Exponential Cutoff in the Size Distribution of the Cold Classical Kuiper Belt

Abstract The cold main classical Kuiper Belt consists of the nonresonant small solar system bodies with low orbital inclinations and orbital semimajor axes between 42.4 and 47.7 au. These objects likely formed in situ, and the population has experienced minimal collisional modification since formation. Using the Outer Solar System Origins Survey ensemble sample and characterization, combined with constraints from deeper surveys and supported by evidence from the Minor Planet Center catalog and the Deep Ecliptic Survey, we determine the absolute magnitude H r distribution of the cold classical belt from H r ≃ 5 to 12 (roughly diameters of 400–20 km). We conclude that the cold population’s H r distribution exhibits an exponential cutoff at large sizes. Exponential cutoffs at large sizes are not a natural outcome of pairwise particle accretion, but exponentially tapered power-law size distributions are a feature of numerical simulations of planetesimal formation via a streaming instability. Our observation of an exponential cutoff agrees with previous observational inferences that no large objects (D ≳ 400 km) exist in the cold population. We note that the asymptotic slope of the H r distribution is consistent with α ∼ 0.4. This asymptotic slope is also found in streaming instability modeling of planetesimal formation and is thus not necessarily associated with achieving collisional equilibrium. Studies of the trans-Neptunian region are providing the parameters that will enable future streaming-instability studies to determine the initial conditions of planetesimal formation in the ≈45 au region of the Sun’s protoplanetary disk.

Long-term Dynamical Stability in the Outer Solar System. I. The Regular and Chaotic Evolution of the 34 Largest Trans-Neptunian Objects

Abstract We carried out an extensive analysis of the stability of the outer solar system, making use of the frequency analysis technique over short-term integrations of nearly 100,000 test particles, as well as a statistical analysis of 200 1 Gyr long numerical simulations, which consider the mutual perturbations of the giant planets and the 34 largest trans-Neptunian objects (we have called all 34 objects “dwarf planets,” DPs, even if probably only the largest of them are true DPs). From the frequency analysis, we produced statistical diffusion maps for a wide region of the a–e phase-space plane; we also present the average diffusion time for orbits as a function of perihelion. We later turned our attention to the 34 DPs, making an individualized analysis for each of them and producing a first approximation of their future stability. From the 200 distinct realizations of the orbital evolution of the 34 DPs, we classified the sample into three categories, including 17 stable, 11 unstable, and 6 resonant objects; we also found that, statistically, two objects from the sample will leave the trans-Neptunian region within the next gigayear, most likely being ejected from the solar system, but with a nonnegligible probability of going inside the orbit of Neptune, either to collide with a giant planet or even falling to the inner solar system, where our simulations are no longer able to resolve their continuous evolution.

Transneptunian Space

We provide a nonspecialist overview of the current state of understanding of the structure and origin of our Solar System's transneptunian region (often called the Kuiper Belt), highlighting perspectives on planetesimal formation, planet migration, and the contextual relationship with protoplanetary disks. We review the dynamical features of the transneptunian populations and their associated differences in physical properties. We describe aspects of our knowledge that have advanced in the past two decades and then move on to current issues of research interest (which thus still have unclear resolution). ▪ The current transneptunian population consists of both implanted and primordial objects. ▪ The primordial (aka cold) population is a largely unaltered remnant of the population that formed in situ. ▪ The reason for the primordial cold population's current outer edge is unexplained. ▪ The large semimajor-axis population now dynamically detached from Neptune is critical for understanding the Solar System's history. ▪ Observational constraints on the number and orbits of distant objects remain poor.

The Main Asteroid Belt: The Primary Source of Debris on Comet-like Orbits

Abstract Jupiter-family comets (JFCs) contribute a significant amount of debris to near-Earth space. However, telescopic observations of these objects seem to suggest that they have short physical lifetimes. If this is true, the material generated will also be short-lived, but fireball observation networks still detect material on cometary orbits. This study examines centimeter-to-meter-scale sporadic meteoroids detected by the Desert Fireball Network from 2014 to 2020 originating from JFC-like orbits. Analyzing each event’s dynamic history and physical characteristics, we confidently determined whether they originated from the main asteroid belt or the trans-Neptunian region. Our results indicate that <4% of sporadic meteoroids on JFC-like orbits are genetically cometary. This observation is statistically significant and shows that cometary material is too friable to survive in near-Earth space. Even when considering shower contributions, meteoroids on JFC-like orbits are primarily from the main belt. Thus, the presence of genuine cometary meteorites in terrestrial collections is highly unlikely.

OSSOS. XXI. Collision Probabilities in the Edgeworth–Kuiper Belt

Abstract Here, we present results on the intrinsic collision probabilities, P I , and range of collision speeds, V I , as a function of the heliocentric distance, r, in the trans-Neptunian region. The collision speed is one of the parameters that serves as a proxy for a collisional outcome (e.g., disruption and scattering of fragments, or formation of a crater, as both processes are related to the impact energy). We utilize an improved and debiased model of the trans-Neptunian object (TNO) region from the “Outer Solar System Origins Survey” (OSSOS). It provides a well-defined model of TNO orbital distribution, based on multiple opposition observations of more than 1000 bodies. We compute collisional probabilities for the OSSOS models of the main classical, resonant, detached+outer, and scattering TNO populations. The intrinsic collision probabilities and collision speeds are computed using Öpik’s approach, as revised and modified by Wetherill for noncircular and inclined orbits. The calculations are carried out for each of the dynamical TNO groups, allowing for inter-population collisions as well as collisions within each TNO population, resulting in 28 combinations in total. Our results indicate that collisions in the trans-Neptunian region are possible over a wide range in (r, V I ) phase space. Although collisions are calculated to happen within r ∼ 20–200 au and V I ∼ 0.1 km s−1 to as high as V I ∼ 9 km s−1, most of the collisions are likely to happen at low relative velocities V I < 1 km s−1 and are dominated by the main classical belt.

Prospects for Future Exploration of the Trans-Neptunian Region

Prospects for Future Exploration of the Trans-Neptunian Region

The origin of (47171) Lempo-like Kuiper belt triple systems during binary–binary interactions

ABSTRACT The trans-Neptunian triple Lempo–Paha–Hiisi is composed of a tight inner binary with components of similar size and an outer companion about half their size orbiting 10 times further away. Large trans-Neptunian objects like Pluto also have multiple small moons, but Lempo’s structure is unique in the Solar system, and the place and timing of its origin is still a subject of debate. We propose a new formation mechanism able to form a large number of systems like Lempo–Paha–Hiisi, which involves binary–binary close encounters in the primordial planetesimal disc at 30–40 au. Some of these systems were then implanted in different populations of the trans-Neptunian region during Neptune’s outward migration. Our results strongly support that the 4:7 resonant multiple object Manwë–Thorondor was once a triple system similar to Lempo–Paha–Hiisi, but the orbit of the inner binary evolved by tides, becoming a contact binary. As with Lempo–Paha–Hiisi, it should have originated in the planetesimal disc below 30–40 au. Triple systems like Lempo–Paha–Hiisi or Manwë–Thorondor could not have formed in situ and the existence of this kind of system is not expected in the cold classical Kuiper belt.

“TNOs are Cool”: A survey of the trans-Neptunian region

The goal of this work is to determine the physical characteristics of resonant, detached and scattered disk objects in the trans-Neptunian region, observed mainly in the framework of the “TNOs are Cool” Herschel open time key programme. Based on thermal emission measurements with the Herschel/PACS and Spitzer/MIPS instruments, we determine size, albedo, and surface thermal properties for 23 objects using radiometric modeling techniques. This is the first analysis in which the physical properties of objects in the outer resonances are determined for a notable sample. In addition to the results for individual objects, we compared these characteristics with the bulk properties of other populations of the trans-Neptunian region. The newly analyzed objects show a large variety of beaming factors, indicating a diversity of thermal properties, and in general they follow the albedo-color clustering identified earlier for Kuiper belt objects and Centaurs, further strengthening the evidence for a compositional discontinuity in the young Solar System.